In the Eigen backend, the preconditioners and direct solvers might reuse the...

In the Eigen backend, the preconditioners and direct solvers might reuse the symbolic pattern of factorization

src/amdis/linear_algebra/eigen/DirectRunner.hpp

@@ -27,12 +27,17 @@ namespace AMDiS
       : solver_{}
     {
       SolverConfig<LU>::init(prefix, solver_);
+      Parameters::get(prefix + "->reuse pattern", reusePattern_);
     }
 
     /// Implementes \ref RunnerInterface::init()
     virtual void init(Matrix const& A) override
     {
-      solver_.compute(A);
+      if (!reusePattern_ || !initialized_) {
+        solver_.analyzePattern(A);
+        initialized_ = true;
+      }
+      solver_.factorize(A);
 
       test_exit(solver_.info() == Eigen::Success,
         "Error in solver.compute(matrix)");
@@ -60,5 +65,7 @@ namespace AMDiS
 
   private:
     LU solver_;
+    bool reusePattern_ = false;
+    bool initialized_ = false;
   };
 }

src/amdis/linear_algebra/eigen/IterativeRunner.hpp

@@ -22,13 +22,18 @@ namespace AMDiS
     {
       SolverCfg::init(prefix, solver_);
       PreconCfg::init(prefix + "->precon", solver_.preconditioner());
+      Parameters::get(prefix + "->reuse pattern", reusePattern_);
     }
 
 
     /// Implementes \ref RunnerInterface::init()
     virtual void init(Matrix const& A) override
     {
-      solver_.compute(A);
+      if (!reusePattern_ || !initialized_) {
+        solver_.analyzePattern(A);
+        initialized_ = true;
+      }
+      solver_.factorize(A);
 
       test_exit(solver_.info() == Eigen::Success,
         "Error in solver.compute(matrix)");
@@ -58,6 +63,8 @@ namespace AMDiS
 
   private:
     IterativeSolver solver_;
+    bool reusePattern_ = false;
+    bool initialized_ = false;
   };
 
 } // end namespace AMDiS