Add non-zero initial guess in test

Author: finsberg

.cspell_dict.txt

@@ -92,6 +92,7 @@ lmbda
 Luca
 MATIS
 Matteo
+maxstep
 mechanicsproblem
 ment
 meshtags
@@ -146,8 +147,8 @@ Taras
 tomek
 Tomek
 TRAME
-unassembled
 transmural
+unassembled
 ureg
 Usyk
 varepsilon

src/pulse/problem.py

@@ -176,12 +176,14 @@ def default_parameters():
                 "snes_type": "newtonls",
                 "snes_atol": 1e-6,
                 "snes_rtol": 1e-10,
-                # "snes_stol": 1e-8,
+                "snes_stol": 1e-8,
                 "snes_max_it": 50,
+                # "snes_linesearch_maxstep": 50,
+                # "snes_view": None,
                 # "snes_type": "newtontr",
                 # "snes_type": "vinewtonrsls",
                 # "snes_linesearch_type": "none",
-                "snes_linesearch_type": "secant",
+                "snes_linesearch_type": "l2",
                 # "mat_mumps_icntl_24": 1,  # Zero pivot detection
                 # "mat_mumps_icntl_25": 0,  # Which nullspace to extract
                 # "mat_mumps_icntl_4": 1,  # Verbosity

tests/test_static_problem.py

@@ -150,7 +150,7 @@ def handle_control_lv(
         def gen(problem):
             dvlp = (target_lvp - initial_lvp) / 3.0
 
-            for lvp in np.arange(0.0, target_lvp + 1e-9, dvlp):
+            for lvp in np.arange(initial_lvp + dvlp, target_lvp + 1e-9, dvlp):
                 print(f"Solve Pressure: {lvp}")
                 traction.assign(lvp)
                 problem.solve()
@@ -166,7 +166,8 @@ def gen(problem):
         def gen(problem):
             dlv = (target_lvv - initial_volume) / 4.0
 
-            for lvv in np.arange(initial_volume, target_lvv + 1e-9, dlv):
+            for lvv in np.arange(initial_volume + dlv, target_lvv + 1e-9, dlv):
+                problem.u.x.array[:] = 0.00001
                 print("Solve Volume: ", lvv)
                 volume.value = lvv
                 problem.solve()
@@ -247,7 +248,7 @@ def test_static_problem_lv(
         },
         cavities=cavities,
     )
-    problem.solve()
+
     for lvp, lvv in gen(problem):
         print(f"LVP: {lvp}, LVV: {lvv}")