fixed how initial condition is stored

PMEAL · Aug 3, 2023 · 7e2d622 · 7e2d622
1 parent 7e75cef
commit 7e2d622
Showing 1 changed file with 1 addition and 2 deletions.
diff --git a/openpnm/contrib/_generic_multiphysics.py b/openpnm/contrib/_generic_multiphysics.py
@@ -104,8 +104,7 @@ def run(self, x0, tspan, interval=None, integrator=None, g_tol=None,
         self.soln = SolutionContainer()
         for alg in algs:
             # transient solution for each alg
-            x0 = alg.x.reshape((alg.Np, 1))
-            print(t0)
+            x0 = alg.x.reshape((alg.Np, 1)).copy()
             soln = TransientSolution(t0, x0)
             self.soln[alg.settings['quantity']] = soln
         phase = network.project.phases[0]