calc = a.calc = GPAW(kpts=kpts,
eigensolver='rmm-diis',
symmetry={'point_group': sym},
mode='pw',
width=0.001,
txt=name + '.txt')
e = a.get_potential_energy()
#calc.diagonalize_full_hamiltonian(nbands=100)
calc.write(name, 'all')
calc = GPAW(name, txt=None, communicator=serial_comm)
chiold = CHI(calc,
w=omega,
q=q_c,
ecut=100,
hilbert_trans=False,
xc='RPA',
G_plus_q=True,
txt=name + '.logold')
chiold.initialize()
chiold.calculate()
chi0old_wGG = chiold.chi0_wGG
chi = Chi0(calc,
omega,
hilbert=False,
ecut=100,
txt=name + '.log')
pd, chi0_wGG, _, _ = chi.calculate(q_c)
assert abs(chi0_wGG - chi0old_wGG).max() < 1e-15
