def test_UHF_SVP(self):
test_dir = 'UHF/SVP/'
with open(self.base_dir + test_dir + self.molden_file, "r") as f:
psi4 = molden2qmc.PSI4(f)
psi4.gwfn()
self.assertTrue(
filecmp.cmp(self.base_dir + test_dir + 'gwfn.data', 'gwfn.data'))
def test_PSI4_RHF_SVP(self):
test_dir = 'PSI4/RHF/SVP/'
with open(self.base_dir + test_dir + self.molden_file, "r") as f:
psi4 = molden2qmc.PSI4(f, qmcpack_normalization=True)
psi4.qmcpack()
res = os.system('h5diff -d 0.000001 Mol.orbs.h5 ' + self.base_dir +
test_dir + 'Mol.orbs.h5')
os.remove("Mol.orbs.h5")
self.assertTrue(res == 0)
def test_RHF_SVP(self):
test_dir = 'RHF/SVP/'
with open(self.base_dir + test_dir + self.molden_file, "r") as f:
psi4 = molden2qmc.PSI4(f)
psi4.gwfn()
self.assertTrue(
filecmp.cmp(self.base_dir + test_dir + 'gwfn.data', 'gwfn.data'))
with open('test/N4/ORCA4/RHF/SVP/N4.molden.input', "r") as f:
orca = molden2qmc.Orca(f)
self.assertTrue(
np.allclose(mo_matrix(psi4), mo_matrix(orca), atol=0.0001))