From 1316296199b0f10c9865237662e63b620246071d Mon Sep 17 00:00:00 2001
From: Charlie Vanaret <vanaret@zib.de>
Date: Fri, 1 Nov 2024 11:47:05 +0100
Subject: [PATCH] Set correct options in MUMPS to better detect null pivots
 (singular matrices)

---
 uno/solvers/MUMPS/MUMPSSolver.cpp | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/uno/solvers/MUMPS/MUMPSSolver.cpp b/uno/solvers/MUMPS/MUMPSSolver.cpp
index 1a7fc5ed..7c3b901c 100644
--- 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() {