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

(grafted from b49ef99617b11404007db283ce68f5dd5b2eff0e
)
This commit is contained in:
Gael Guennebaud 2014-08-21 22:14:25 +02:00
parent 0321449944
commit c7331ebb06

View File

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