diff options
| author |  Jan Blechta <blechta@karlin.mff.cuni.cz> | 2015-02-18 14:26:10 +0100 |
|---|---|---|
| committer | Jan Blechta <blechta@karlin.mff.cuni.cz> | 2015-02-18 14:26:10 +0100 |
| commit | 168ceb271e631c17da1c722119746ef5931fdca5 (patch) | |
| tree | c4a70738bb5e6313bd055ef265e004e7efdfc0af /Eigen/src/IterativeLinearSolvers | |
| parent | 8fdcaded5ece5691e5166a04f4255da217ae5d30 (diff) | |
Really use zero guess in ConjugateGradients::solve as documented
and expected for consistency with other methods.
Diffstat (limited to 'Eigen/src/IterativeLinearSolvers')
| -rw-r--r-- | Eigen/src/IterativeLinearSolvers/ConjugateGradient.h | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/Eigen/src/IterativeLinearSolvers/ConjugateGradient.h b/Eigen/src/IterativeLinearSolvers/ConjugateGradient.h index 10cd94783..fab5fdb1f 100644 --- a/Eigen/src/IterativeLinearSolvers/ConjugateGradient.h +++ b/Eigen/src/IterativeLinearSolvers/ConjugateGradient.h @@ -208,7 +208,7 @@ public: template<typename Rhs,typename Dest> void _solve_impl(const MatrixBase<Rhs>& b, Dest& x) const { - x.setOnes(); + x.setZero(); _solve_with_guess_impl(b.derived(),x); } |