Tweaked mesh and solver for fracture problem.

alucantonio · Jul 14, 2023 · 9b57134 · 9b57134
1 parent 5226e20
commit 9b57134
Show file tree

Hide file tree

Showing 2 changed files with 150 additions and 57 deletions.
diff --git a/examples/pacman_fracture.ipynb b/examples/pacman_fracture.ipynb
diff --git a/src/dctkit/math/opt/optctrl.py b/src/dctkit/math/opt/optctrl.py
@@ -140,7 +140,7 @@ def __init__(self, prb: OptimizationProblem):
         init()
         # create default solver and settings
         self.tao = PETSc.TAO().create()
-        self.tao.setType(PETSc.TAO.Type.LMVM)  # Specify the solver type
+        self.tao.setType(PETSc.TAO.Type.BQNLS)  # Specify the solver type
         # create variable to store the gradient of the objective function
         self.g = PETSc.Vec().createSeq(self.prb.dim)