diff --git a/uno/linear_algebra/SymmetricIndefiniteLinearSystem.hpp b/uno/linear_algebra/SymmetricIndefiniteLinearSystem.hpp index 37aa9bf8..b53d072c 100644 --- a/uno/linear_algebra/SymmetricIndefiniteLinearSystem.hpp +++ b/uno/linear_algebra/SymmetricIndefiniteLinearSystem.hpp @@ -91,7 +91,7 @@ namespace uno { template void SymmetricIndefiniteLinearSystem::factorize_matrix(DirectSymmetricIndefiniteLinearSolver& 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 index 713f473a..62d6482f 100644 --- 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); }