Set correct options in MUMPS to better detect null pivots (singular m…

…atrices)
cvanaret · Nov 1, 2024 · 1316296 · 1316296
1 parent bcea64b
commit 1316296
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/uno/solvers/MUMPS/MUMPSSolver.cpp b/uno/solvers/MUMPS/MUMPSSolver.cpp
@@ -28,6 +28,8 @@ namespace uno {
       this->mumps_structure.icntl[1] = -1;
       this->mumps_structure.icntl[2] = -1;
       this->mumps_structure.icntl[3] = 0;
+      this->mumps_structure.icntl[12] = 1;
+      this->mumps_structure.icntl[23] = 1; // ICNTL(24) controls the detection of “null pivot rows”
    }
 
    MUMPSSolver::~MUMPSSolver() {