diff options
author | Gael Guennebaud <g.gael@free.fr> | 2016-07-04 17:19:38 +0200 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2016-07-04 17:19:38 +0200 |
commit | fbcfc2f86225891dd21ccdefa4a01b7f9409a992 (patch) | |
tree | daf1ce97fd7cbc3faf25f6c7810ef6fcf8bd3ffd /test/sparse_solver.h | |
parent | 7f7839c12f9bbd88bee1f3c9d7f5dbbe556045d4 (diff) |
Add unit test for solveWithGuess, and fix template resolution.
Diffstat (limited to 'test/sparse_solver.h')
-rw-r--r-- | test/sparse_solver.h | 21 |
1 files changed, 21 insertions, 0 deletions
diff --git a/test/sparse_solver.h b/test/sparse_solver.h index b67653496..8cf35bef7 100644 --- a/test/sparse_solver.h +++ b/test/sparse_solver.h @@ -11,6 +11,21 @@ #include <Eigen/SparseCore> #include <sstream> +template<typename Solver, typename Rhs, typename Guess,typename Result> +void solve_with_guess(IterativeSolverBase<Solver>& solver, const MatrixBase<Rhs>& b, const Guess& g, Result &x) { + x = solver.derived().solveWithGuess(b.derived(),g); +} + +template<typename Solver, typename Rhs, typename Guess,typename Result> +void solve_with_guess(SparseSolverBase<Solver>& solver, const MatrixBase<Rhs>& b, const Guess& , Result& x) { + x = solver.derived().solve(b); +} + +template<typename Solver, typename Rhs, typename Guess,typename Result> +void solve_with_guess(SparseSolverBase<Solver>& solver, const SparseMatrixBase<Rhs>& b, const Guess& , Result& x) { + x = solver.derived().solve(b); +} + template<typename Solver, typename Rhs, typename DenseMat, typename DenseRhs> void check_sparse_solving(Solver& solver, const typename Solver::MatrixType& A, const Rhs& b, const DenseMat& dA, const DenseRhs& db) { @@ -37,6 +52,12 @@ void check_sparse_solving(Solver& solver, const typename Solver::MatrixType& A, } VERIFY(oldb.isApprox(b) && "sparse solver testing: the rhs should not be modified!"); VERIFY(x.isApprox(refX,test_precision<Scalar>())); + + x.setZero(); + solve_with_guess(solver, b, x, x); + VERIFY(solver.info() == Success && "solving failed when using analyzePattern/factorize API"); + VERIFY(oldb.isApprox(b) && "sparse solver testing: the rhs should not be modified!"); + VERIFY(x.isApprox(refX,test_precision<Scalar>())); x.setZero(); // test the analyze/factorize API |