using different direct solver

pmgbergen · Aug 22, 2024 · d58e565 · d58e565
1 parent f17b453
commit d58e565
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/src/darsia/measure/wasserstein.py b/src/darsia/measure/wasserstein.py
@@ -114,6 +114,7 @@ def __init__(
                     - "direct": direct solver
                     - "amg": algebraic multigrid solver
                     - "cg": conjugate gradient solver preconditioned with AMG
+                    - "ksp": PETSc KSP solver
                 - formulation (str): formulation of the linear system. Defaults to
                     "pressure". Supported formulations are:
                     - "full": full system
@@ -462,6 +463,7 @@ def setup_ksp_solver(self, matrix: sps.csc_matrix) -> None:
         kernel = kernel / np.linalg.norm(kernel)
         nullspace=[kernel]
         self.linear_solver = darsia.linalg.KSP(matrix)#, nullspace=nullspace)
+
 
         # Define solver options
         linear_solver_options = self.options.get("linear_solver_options", {})
@@ -472,7 +474,8 @@ def setup_ksp_solver(self, matrix: sps.csc_matrix) -> None:
         if approach == "direct":
             self.solver_options = {
                 "ksp_type": "preonly",
-                'pc_type': 'lu',   
+                'pc_type': 'lu',
+                "pc_factor_mat_solver_type": "mumps",  
             }
         else:
             self.solver_options = {