Merge pull request #61 from cvanaret/cvanaret-mumps-fix

Set correct options in MUMPS to better detect null pivots (singular matrices)
cvanaret · Nov 1, 2024 · c56cbbd · c56cbbd
2 parents bcea64b + 1316296
commit c56cbbd
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() {