[gpaw-users] AssertionError in transformers.py
Jens Jørgen Mortensen
jensj at fysik.dtu.dk
Thu Feb 2 08:56:08 CET 2012
On 01-02-2012 10:40, Falco Hüser wrote:
> Hi,
> can anyone tell me why I get this error in a regular groundstate
> calculation:
>
> GPAW CLEANUP (node 13): exceptions.AssertionError occurred. Calling
> MPI_Abort!
> File "/home/niflheim/falh/mygpaw/gpaw/transformers.py", line 48, in
> __init__
> assert ((gdin.n_c[0] + 2 * nn - 1) * (gdin.n_c[1] + 2 * nn - 1)<=
> AssertionError
This happens when there are too few grid points per cpu - I think. Try
with smaller h or fewer cpu's.
Jens Jørgen
> the script is:
>
> import numpy as np
> from ase import Atom, Atoms
> from ase.structure import bulk
> from ase.units import Hartree, Bohr
> from gpaw import GPAW, FermiDirac
>
> a = 3.567
>
> atoms = bulk('C', 'diamond', a=a)
> calc = GPAW(
> h=0.20,
> kpts=(9,9,9),
> xc='LDA',
> basis='dzp',
> txt='groundstate.txt',
> nbands=300,
> eigensolver='cg',
> occupations=FermiDirac(0.001),
> convergence={'bands':250}
> )
>
> atoms.set_calculator(calc)
> atoms.get_potential_energy()
> calc.write('C_kpt9.gpw','all')
>
> I used exactly the same script with 7x7x7 kpoints and had no problems.
> The same happens also for other structures.
>
> Thanks,
> Falco
> _______________________________________________
> gpaw-users mailing list
> gpaw-users at listserv.fysik.dtu.dk
> https://listserv.fysik.dtu.dk/mailman/listinfo/gpaw-users
More information about the gpaw-users
mailing list