From b0adbfbae716bbeac362aae1d931de24b8bb83f5 Mon Sep 17 00:00:00 2001
From: Gael Guennebaud <g.gael@free.fr>
Date: Wed, 21 Sep 2011 18:08:08 +0200
Subject: [PATCH] BiCGSTAB does not like starting from 0...

---
 unsupported/Eigen/src/IterativeSolvers/BiCGSTAB.h | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/unsupported/Eigen/src/IterativeSolvers/BiCGSTAB.h b/unsupported/Eigen/src/IterativeSolvers/BiCGSTAB.h
index 4f365bd96..d04158bb6 100644
--- a/unsupported/Eigen/src/IterativeSolvers/BiCGSTAB.h
+++ b/unsupported/Eigen/src/IterativeSolvers/BiCGSTAB.h
@@ -243,7 +243,7 @@ struct solve_retval<BiCGSTAB<_MatrixType, _Preconditioner>, Rhs>
 
   template<typename Dest> void evalTo(Dest& dst) const
   {
-    dst.setZero();
+    dst.setOnes();
     dec()._solve(rhs(),dst);
   }
 };