bug #897: fix regression in BiCGSTAB(mat) ctor (an all other iterative solvers).

Add respective regression unit test.

1 files changed, 9 insertions, 0 deletions

diff --git a/test/sparse_solver.h b/test/sparse_solver.h
index 42b365eaa..45cfdad25 100644
--- a/test/sparse_solver.h
+++ b/test/sparse_solver.h
@@ -83,6 +83,15 @@ void check_sparse_solving(Solver& solver, const typename Solver::MatrixType& A,
     VERIFY(xm.isApprox(refX,test_precision<Scalar>()));
   }
   
+  // test initialization ctor
+  {
+    Rhs x(b.rows(), b.cols());
+    Solver solver2(A);
+    VERIFY(solver2.info() == Success);
+    x = solver2.solve(b);
+    VERIFY(x.isApprox(refX,test_precision<Scalar>()));
+  }
+  
   // test dense Block as the result and rhs:
   {
     DenseRhs x(db.rows(), db.cols());