From a72bf09e6d4f9bcbe4c39a85617a6a8ccf34a36d Mon Sep 17 00:00:00 2001
From: Jan Blechta <blechta@karlin.mff.cuni.cz>
Date: Wed, 18 Feb 2015 14:26:10 +0100
Subject: [PATCH] Really use zero guess in ConjugateGradients::solve as
 documented and expected for consistency with other methods.

---
 Eigen/src/IterativeLinearSolvers/ConjugateGradient.h | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/Eigen/src/IterativeLinearSolvers/ConjugateGradient.h b/Eigen/src/IterativeLinearSolvers/ConjugateGradient.h
index 35ec189ff..8ba4a8dbe 100644
--- a/Eigen/src/IterativeLinearSolvers/ConjugateGradient.h
+++ b/Eigen/src/IterativeLinearSolvers/ConjugateGradient.h
@@ -224,7 +224,7 @@ public:
   template<typename Rhs,typename Dest>
   void _solve(const Rhs& b, Dest& x) const
   {
-    x.setOnes();
+    x.setZero();
     _solveWithGuess(b,x);
   }