Using Scalapack does not work
This example works without Scalapack and gives strange numbers with (on 4 cores):
from ase.build import molecule
from gpaw import GPAW
system = molecule('H2O')
system.center(2.0)
system.calc = GPAW(mode='pw',
nbands=6,
parallel=dict(domain=(1, 1, 1),
band=4,
sl_default=(2, 1, 2)))
system.get_potential_energy()