Saddle point distributed in work

Wells-Group · Nov 14, 2024 · d416df5 · d416df5
1 parent 15db039
commit d416df5
Showing 1 changed file with 5 additions and 0 deletions.
diff --git a/demo/saddle_point_dlrbnicsx_distributed/dlrbnicsx_saddle_point_distributed.py b/demo/saddle_point_dlrbnicsx_distributed/dlrbnicsx_saddle_point_distributed.py
@@ -319,11 +319,16 @@ def solve(self, mu):
         ksp = PETSc.KSP()
         ksp.create(mesh.comm)
         ksp.setOperators(A)
+        '''
         # Set GMRES solver
         ksp.setType("gmres")
         ksp.setGMRESRestart(100)
         pc = ksp.getPC()
         pc.setType("bjacobi")
+        '''
+        ksp.setType("preonly")
+        pc = ksp.getPC()
+        pc.setType("lu")
         # Convergence criteria based on residual tolerance
         ksp.rtol = 1.e-8
         # Solve and see convergence details