diff --git a/tests/test_evgw.py b/tests/test_evgw.py
index ada5bd40..fdb98ca7 100644
--- a/tests/test_evgw.py
+++ b/tests/test_evgw.py
@@ -71,8 +71,8 @@ def _test_regression(self, xc, kwargs, nmom_max, ip, ea, name=""):
         gw.kernel(nmom_max)
         gf = gw.gf.physical(weight=0.1)
         self.assertTrue(gw.converged)
-        self.assertAlmostEqual(gf.occupied().energies[-1], ip, 7, msg=name)
-        self.assertAlmostEqual(gf.virtual().energies[0], ea, 7, msg=name)
+        self.assertAlmostEqual(gf.occupied().energies[-1], ip, 6, msg=name)
+        self.assertAlmostEqual(gf.virtual().energies[0], ea, 6, msg=name)
 
     def test_regression_simple(self):
         ip = -0.278612876943
@@ -90,7 +90,7 @@ def test_regression_g0w(self):
         self._test_regression("hf", dict(g0=True, damping=0.5), 1, ip, ea, "g0w")
 
     def test_regression_pbe_fock_loop(self):
-        ip = -0.281393616901
+        ip = -0.28139405242856597
         ea = 0.007257181880
         self._test_regression("pbe", dict(), 1, ip, ea, "pbe")
 
diff --git a/tests/test_gw.py b/tests/test_gw.py
index 89eabdf7..d9d9c2b7 100644
--- a/tests/test_gw.py
+++ b/tests/test_gw.py
@@ -201,8 +201,8 @@ def _test_regression(self, xc, kwargs, nmom_max, ip, ea, name=""):
         gw = GW(mf, **kwargs)
         gw.kernel(nmom_max)
         gf = gw.gf.physical(weight=0.1)
-        self.assertAlmostEqual(gf.occupied().energies[-1], ip, 7, msg=name)
-        self.assertAlmostEqual(gf.virtual().energies[0], ea, 7, msg=name)
+        self.assertAlmostEqual(gf.occupied().energies[-1], ip, 6, msg=name)
+        self.assertAlmostEqual(gf.virtual().energies[0], ea, 6, msg=name)
 
     def test_regression_simple(self):
         ip = -0.277578450082
@@ -210,7 +210,7 @@ def test_regression_simple(self):
         self._test_regression("hf", dict(), 3, ip, ea, "simple")
 
     def test_regression_pbe(self):
-        ip = -0.233369739990
+        ip = -0.2333686494079075
         ea = 0.002658170914
         self._test_regression("pbe", dict(), 3, ip, ea, "pbe")
 
@@ -220,8 +220,8 @@ def test_regression_fock_loop(self):
         self._test_regression("hf", dict(fock_loop=True), 1, ip, ea, "fock loop")
 
     def test_diagonal_pbe0(self):
-        ip = -0.261876372990
-        ea = 0.008159826670
+        ip = -0.2618765321330203
+        ea = 0.00815559003481674
         self._test_regression("pbe0", dict(diagonal_se=True), 5, ip, ea, "diagonal pbe0")
 
     def test_regression_tda(self):
diff --git a/tests/test_kthcgw.py b/tests/test_kthcgw.py
index 5a5f8222..a8ba21f3 100644
--- a/tests/test_kthcgw.py
+++ b/tests/test_kthcgw.py
@@ -173,12 +173,12 @@ def _test_regression(self, xc, kwargs, nmom_max, ip, ea, name=""):
         for k in range(len(kpts)):
             gf = kgw.gf[k].physical(weight=0.1)
             self.assertTrue(kgw.converged)
-            self.assertAlmostEqual(gf.occupied().energies[-1], ip[k], 7, msg=name)
-            self.assertAlmostEqual(gf.virtual().energies[0], ea[k], 7, msg=name)
+            self.assertAlmostEqual(gf.occupied().energies[-1], ip[k], 6, msg=name)
+            self.assertAlmostEqual(gf.virtual().energies[0], ea[k], 6, msg=name)
 
     def test_regression_pbe_fock_loop(self):
         ip = [
-            -0.48261234253482393,
+            -0.48261216353603703,
             -0.5020081305060984,
             -0.5020430194398028,
             -0.5098036236597017,
diff --git a/tests/test_qsgw.py b/tests/test_qsgw.py
index 4ce868de..ac495801 100644
--- a/tests/test_qsgw.py
+++ b/tests/test_qsgw.py
@@ -35,8 +35,8 @@ def _test_regression(self, xc, kwargs, nmom_max, ip, ea, ip_full, ea_full, name=
         self.assertTrue(gw.converged)
         self.assertAlmostEqual(np.max(qp_energy[gw.mo_occ > 0]), ip, 6, msg=name)
         self.assertAlmostEqual(np.min(qp_energy[gw.mo_occ == 0]), ea, 6, msg=name)
-        self.assertAlmostEqual(gf.occupied().energies[-1], ip_full, 6, msg=name)
-        self.assertAlmostEqual(gf.virtual().energies[0], ea_full, 6, msg=name)
+        self.assertAlmostEqual(gf.occupied().energies[-1], ip_full, 5, msg=name)
+        self.assertAlmostEqual(gf.virtual().energies[0], ea_full, 5, msg=name)
 
     def test_regression_simple(self):
         # Quasiparticle energies:
@@ -52,7 +52,7 @@ def test_regression_pbe_srg(self):
         ip = -0.271017148443
         ea = 0.075779760938
         # GF poles:
-        ip_full = -0.393500679528
+        ip_full = -0.39350150932022504
         ea_full = 0.170103953696
         self._test_regression("pbe", dict(srg=1000), 3, ip, ea, ip_full, ea_full, "pbe srg")
 
diff --git a/tests/test_thcgw.py b/tests/test_thcgw.py
index e81c89ca..210eed19 100644
--- a/tests/test_thcgw.py
+++ b/tests/test_thcgw.py
@@ -18,7 +18,7 @@
 class Test_THCTDA(unittest.TestCase):
     @classmethod
     def setUpClass(cls):
-        cell = gto.M()
+        cell = gto.Cell()
         cell.a = np.eye(3) * 3
         cell.atom = """He 0 0 0; He 1 1 1"""
         cell.basis = "6-31g"
@@ -122,7 +122,7 @@ def test_moments_vs_cderi(self):
             self.assertAlmostEqual(dif, 0, 8)
 
     def _test_regression(self, xc, kwargs, nmom_max, ip, ea, name=""):
-        cell = gto.M()
+        cell = gto.Cell()
         cell.a = np.eye(3) * 3
         cell.atom = """He 0 0 0; He 1 1 1"""
         cell.basis = "6-31g"
@@ -147,12 +147,12 @@ def _test_regression(self, xc, kwargs, nmom_max, ip, ea, name=""):
         gw.kernel(nmom_max)
         gf = gw.gf.physical(weight=0.1)
         self.assertTrue(gw.converged)
-        self.assertAlmostEqual(gf.occupied().energies[-1], ip, 7, msg=name)
-        self.assertAlmostEqual(gf.virtual().energies[0], ea, 7, msg=name)
+        self.assertAlmostEqual(gf.occupied().energies[-1], ip, 6, msg=name)
+        self.assertAlmostEqual(gf.virtual().energies[0], ea, 6, msg=name)
 
     def test_regression_pbe_fock_loop(self):
-        ip = -0.2786188906832294
-        ea = 1.0822831284078982
+        ip = -0.2786185073019116
+        ea = 1.082283812661612
         self._test_regression("pbe", dict(), 1, ip, ea, "pbe")