From 8b6c1caa3e41d5659e35b1852b791a0f9bd880d2 Mon Sep 17 00:00:00 2001
From: Gael Guennebaud <g.gael@free.fr>
Date: Fri, 14 Jan 2011 20:29:55 +0100
Subject: [PATCH] fix compilation of rowmajor sparse time diagonal

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

diff --git a/Eigen/src/Sparse/SparseDiagonalProduct.h b/Eigen/src/Sparse/SparseDiagonalProduct.h
index 82b2dff25..994bf163e 100644
--- a/Eigen/src/Sparse/SparseDiagonalProduct.h
+++ b/Eigen/src/Sparse/SparseDiagonalProduct.h
@@ -166,18 +166,17 @@ class sparse_diagonal_product_inner_iterator_selector
   : public CwiseBinaryOp<
       scalar_product_op<typename Rhs::Scalar>,
       SparseInnerVectorSet<Lhs,1>,
-      Transpose<typename Rhs::DiagonalVectorType> >::InnerIterator
+      Transpose<const typename Rhs::DiagonalVectorType> >::InnerIterator
 {
     typedef typename CwiseBinaryOp<
       scalar_product_op<typename Rhs::Scalar>,
       SparseInnerVectorSet<Lhs,1>,
-      Transpose<typename Rhs::DiagonalVectorType> >::InnerIterator Base;
+      Transpose<const typename Rhs::DiagonalVectorType> >::InnerIterator Base;
     typedef typename Lhs::Index Index;
   public:
     inline sparse_diagonal_product_inner_iterator_selector(
               const SparseDiagonalProductType& expr, Index outer)
-      : Base(expr.lhs().innerVector(outer) .cwiseProduct(expr.rhs().const_cast_derived().diagonal().transpose()), 0)
-       // the const_cast_derived above is to get it to compile. once Sparse is const correct, that shouldn't be needed anymore.
+      : Base(expr.lhs().innerVector(outer) .cwiseProduct(expr.rhs().diagonal().transpose()), 0)
     {}
 };