Do not apply the preconditioner before starting the iterations as this might destroy a very good initial guess.

This commit is contained in:
Gael Guennebaud 2014-08-21 22:14:25 +02:00
parent 9c0aa81fbf
commit b49ef99617

View File

@ -39,7 +39,6 @@ bool bicgstab(const MatrixType& mat, const Rhs& rhs, Dest& x,
int maxIters = iters; int maxIters = iters;
int n = mat.cols(); int n = mat.cols();
x = precond.solve(x);
VectorType r = rhs - mat * x; VectorType r = rhs - mat * x;
VectorType r0 = r; VectorType r0 = r;