From 8fa7e92e77f550efb6d67635f19dfd7672d3259c Mon Sep 17 00:00:00 2001
From: Gael Guennebaud <g.gael@free.fr>
Date: Tue, 26 Jul 2011 09:02:41 +0200
Subject: [PATCH] fix sparse selfadjoint time dense such that the other
 triangular part is not used at all

---
 Eigen/src/Sparse/SparseSelfAdjointView.h | 10 +++++++---
 1 file changed, 7 insertions(+), 3 deletions(-)

diff --git a/Eigen/src/Sparse/SparseSelfAdjointView.h b/Eigen/src/Sparse/SparseSelfAdjointView.h
index bdee96452..866c315ba 100644
--- a/Eigen/src/Sparse/SparseSelfAdjointView.h
+++ b/Eigen/src/Sparse/SparseSelfAdjointView.h
@@ -227,10 +227,14 @@ class SparseSelfAdjointTimeDenseProduct
       for (Index j=0; j<m_lhs.outerSize(); ++j)
       {
         LhsInnerIterator i(m_lhs,j);
-        if (ProcessSecondHalf && i && (i.index()==j))
+        if (ProcessSecondHalf)
         {
-          dest.row(j) += i.value() * m_rhs.row(j);
-          ++i;
+          while (i.index()<j) ++i;
+          if(i && i.index()==j)
+          {
+            dest.row(j) += i.value() * m_rhs.row(j);
+            ++i;
+          }
         }
         for(; (ProcessFirstHalf ? i && i.index() < j : i) ; ++i)
         {