[gpaw-users] Re: dipole momentet från en gammal räkning?

[Jens Jørgen Mortensen]

On Mon, 2010-05-31 at 15:44 +0200, Anders Hellman wrote:
> Hej JJ,
> 
> Är det något speciallet som jag måste tänka på om jag vill ta ut dipole momentet från en gammal räkning?
> 
> Att skriva 
> 
> atoms,calc=restart('out.gpw') 
> calc.get_dipole_moment()
> 
> verkar inte fungera. 
> 
> >>> calc.get_dipole_moment()
> Traceback (most recent call last):
>   File "<stdin>", line 1, in ?
>   File "/home/a/andhe/Public/apps/lib64/python/gpaw/aseinterface.py", line 644, in get_dipole_moment
>     return self.density.finegd.calculate_dipole_moment(rhot_g) * Bohr
>   File "/home/a/andhe/Public/apps/lib64/python/gpaw/grid_descriptor.py", line 432, in calculate_dipole_moment
>     rho_01 = rho_g.sum(axis=2)
> AttributeError: 'NoneType' object has no attribute 'sum'
> 
> Det är bara ett problem om jag startar om, inte annars.

Try to do "calc.set_positions()" before "calc.get_dipole_moment()".

JJ

> mvh
> Anders