Temporarily restore analysis at each iteration

cvanaret · Dec 5, 2024 · 30c7304 · 30c7304
1 parent 5f08f59
commit 30c7304
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 1 deletion.
diff --git a/uno/linear_algebra/SymmetricIndefiniteLinearSystem.hpp b/uno/linear_algebra/SymmetricIndefiniteLinearSystem.hpp
@@ -91,7 +91,7 @@ namespace uno {
    template <typename ElementType>
    void SymmetricIndefiniteLinearSystem<ElementType>::factorize_matrix(DirectSymmetricIndefiniteLinearSolver<size_t, ElementType>& linear_solver,
          const WarmstartInformation& warmstart_information) {
-      if (warmstart_information.problem_structure_changed) {
+      if (true) {
          DEBUG << "Performing symbolic analysis of the indefinite system\n";
          linear_solver.do_symbolic_analysis(this->matrix);
       }

diff --git a/uno/solvers/MUMPS/MUMPSSolver.cpp b/uno/solvers/MUMPS/MUMPSSolver.cpp
@@ -35,6 +35,13 @@ namespace uno {
       this->mumps_structure.icntl[12] = 1;
       this->mumps_structure.icntl[23] = 1; // ICNTL(24) controls the detection of “null pivot rows”
 
+      /*
+      // debug for MUMPS team
+      this->mumps_structure.icntl[1] = 6; // ICNTL(2)=6
+      this->mumps_structure.icntl[2] = 6; // ICNTL(3)=6
+      this->mumps_structure.icntl[3] = 6; // ICNTL(4)=2
+       */
+
       this->row_indices.reserve(number_nonzeros);
       this->column_indices.reserve(number_nonzeros);
    }