Really use zero guess in ConjugateGradients::solve as documented

and expected for consistency with other methods.
This commit is contained in:
Jan Blechta 2015-02-18 14:26:10 +01:00
parent bb3e5c29cc
commit a72bf09e6d

View File

@ -224,7 +224,7 @@ public:
template<typename Rhs,typename Dest> template<typename Rhs,typename Dest>
void _solve(const Rhs& b, Dest& x) const void _solve(const Rhs& b, Dest& x) const
{ {
x.setOnes(); x.setZero();
_solveWithGuess(b,x); _solveWithGuess(b,x);
} }