Compute symbolic factorization only when the problem structure change…

…d (once at the beginning and at every phase switch)
cvanaret · Nov 19, 2024 · 8166937 · 8166937
1 parent a1b3b9f
commit 8166937
Showing 1 changed file with 2 additions and 3 deletions.
diff --git a/uno/linear_algebra/SymmetricIndefiniteLinearSystem.hpp b/uno/linear_algebra/SymmetricIndefiniteLinearSystem.hpp
@@ -90,9 +90,8 @@ namespace uno {
 
    template <typename ElementType>
    void SymmetricIndefiniteLinearSystem<ElementType>::factorize_matrix(DirectSymmetricIndefiniteLinearSolver<size_t, ElementType>& linear_solver,
-         const WarmstartInformation& /*warmstart_information*/) {
-      // TODO compute the symbolic factorization only when the problem definition changes
-      if (true) {
+         const WarmstartInformation& warmstart_information) {
+      if (warmstart_information.problem_structure_changed) {
          DEBUG << "Performing symbolic factorization of the indefinite system\n";
          linear_solver.do_symbolic_factorization(this->matrix);
       }