diff --git a/MathLib/LinAlg/Eigen/EigenLinearSolver.cpp b/MathLib/LinAlg/Eigen/EigenLinearSolver.cpp index e1bbcbc83fe3269becd46064a10933e6e508aba4..3c66f5cef42357099dde496e6035a0761fc1b7fd 100644 --- a/MathLib/LinAlg/Eigen/EigenLinearSolver.cpp +++ b/MathLib/LinAlg/Eigen/EigenLinearSolver.cpp @@ -32,16 +32,17 @@ public: explicit EigenDirectLinearSolver(EigenMatrix::RawMatrixType &A) : _A(A) { INFO("-> initialize with the coefficient matrix"); - _solver.compute(A); - if(_solver.info()!=Eigen::Success) { - ERR("Failed during Eigen linear solver initialization"); - return; - } } void solve(EigenVector::RawVectorType &b, EigenVector::RawVectorType &x, EigenOption &/*opt*/) override { INFO("-> solve"); + _solver.compute(_A); + if(_solver.info()!=Eigen::Success) { + ERR("Failed during Eigen linear solver initialization"); + return; + } + x = _solver.solve(b); if(_solver.info()!=Eigen::Success) { ERR("Failed during Eigen linear solve"); @@ -62,11 +63,6 @@ public: explicit EigenIterativeLinearSolver(EigenMatrix::RawMatrixType &A) : _A(A) { INFO("-> initialize with the coefficient matrix"); - _solver.compute(A); - if(_solver.info()!=Eigen::Success) { - ERR("Failed during Eigen linear solver initialization"); - return; - } } void solve(EigenVector::RawVectorType &b, EigenVector::RawVectorType &x, EigenOption &opt) override @@ -74,6 +70,11 @@ public: INFO("-> solve"); _solver.setTolerance(opt.error_tolerance); _solver.setMaxIterations(opt.max_iterations); + _solver.compute(_A); + if(_solver.info()!=Eigen::Success) { + ERR("Failed during Eigen linear solver initialization"); + return; + } x = _solver.solveWithGuess(b, x); if(_solver.info()!=Eigen::Success) { ERR("Failed during Eigen linear solve");