diff --git a/tests/test_meanfield.py b/tests/test_meanfield.py
index 3c670469780556a465a465a1409bcbbb33b05943..bf0a694ef946c1d995f210fd41dab7c7f257496f 100644
--- a/tests/test_meanfield.py
+++ b/tests/test_meanfield.py
@@ -10,4 +10,4 @@ def test_meanfield():
     network_params = {}
     theory_params = {}
     M0 = MultiAreaModel(network_params, theory=True, theory_spec=theory_params)
-    p, r0 = M0.theory.integrate_siegert()
+    p, r0 = M0.theory.integrate_siegert_nest()
diff --git a/tests/test_replace_cc.py b/tests/test_replace_cc.py
index 2b2791d5bd6ebea87e7a3856cbedc0d7ea096600..421361112c3abaa8a9fb9b181eb25704248a5fdf 100644
--- a/tests/test_replace_cc.py
+++ b/tests/test_replace_cc.py
@@ -16,7 +16,7 @@ def test_het_poisson_stat_mf():
     network_params = {}
     theory_params = {}
     M0 = MultiAreaModel(network_params, theory=True, theory_spec=theory_params)
-    p, r0 = M0.theory.integrate_siegert()
+    p, r0 = M0.theory.integrate_siegert_nest()
 
     rates = vector_to_dict(r0[:, -1], M0.area_list, M0.structure)
     with open('mf_rates.json', 'w') as f:
@@ -26,7 +26,7 @@ def test_het_poisson_stat_mf():
                                             'replace_cc_input_source': 'mf_rates.json'}}
     theory_params = {}
     M = MultiAreaModel(network_params, theory=True, theory_spec=theory_params)
-    p, r = M.theory.integrate_siegert()
+    p, r = M.theory.integrate_siegert_nest()
 
     assert(np.allclose(r0[:, -1], r[:, -1]))
 
@@ -35,7 +35,7 @@ def test_hom_poisson_stat_mf():
     network_params = {'connection_params': {'replace_cc': 'hom_poisson_stat'}}
     theory_params = {}
     M = MultiAreaModel(network_params, theory=True, theory_spec=theory_params)
-    p, r = M.theory.integrate_siegert()
+    p, r = M.theory.integrate_siegert_nest()
 
     mu, sigma = M.theory.replace_cc_input()
     # Test for V1