LS: properly set up KSP debug options for the split system

optimad · May 29, 2024 · ee7a287 · ee7a287
1 parent ef2cf6a
commit ee7a287
Show file tree

Hide file tree

Showing 2 changed files with 29 additions and 0 deletions.
diff --git a/src/LA/system_solvers_split.cpp b/src/LA/system_solvers_split.cpp
@@ -1572,6 +1572,33 @@ void SplitSystemSolver::setupSplitKrylovs()
     }
 }
 
+/*!
+ * Perform actions before Krylov subspace method setup.
+ */
+void SplitSystemSolver::preKrylovSetupActions()
+{
+    // Execute base actions
+    SystemSolver::preKrylovSetupActions();
+
+    // Set KSP debug options
+    if (m_debug) {
+        int nSplits = getSplitCount();
+        for (int split = 0; split < nSplits; ++split) {
+            std::string prefix = "fieldsplit_" + std::to_string(split) + "_";
+
+#if (PETSC_VERSION_MAJOR >= 3 && PETSC_VERSION_MINOR >= 7)
+            PetscOptionsSetValue(PETSC_NULLPTR, ("-" + prefix + "ksp_monitor_true_residual").c_str(), "");
+            PetscOptionsSetValue(PETSC_NULLPTR, ("-" + prefix + "ksp_converged_reason").c_str(), "");
+            PetscOptionsSetValue(PETSC_NULLPTR, ("-" + prefix + "ksp_monitor_singular_value").c_str(), "");
+#else
+            PetscOptionsSetValue(("-" + prefix + "ksp_monitor_true_residual").c_str(), "");
+            PetscOptionsSetValue(("-" + prefix + "ksp_converged_reason").c_str(), "");
+            PetscOptionsSetValue(("-" + prefix + "ksp_monitor_singular_value").c_str(), "");
+#endif
+        }
+    }
+}
+
 /*!
  * Generate the split permutation.
  *

diff --git a/src/LA/system_solvers_split.hpp b/src/LA/system_solvers_split.hpp
@@ -142,6 +142,8 @@ class SplitSystemSolver : public SystemSolver {
     using SystemSolver::setupKrylov;
     virtual void setupSplitKrylovs();
 
+    void preKrylovSetupActions() override;
+
     void postKSPSolveActions() override;
 
     void initializeKSPOptions() override;