diff --git a/Eigen/src/SparseCore/SparseBlock.h b/Eigen/src/SparseCore/SparseBlock.h index a6cde9df3..2b092dc21 100644 --- a/Eigen/src/SparseCore/SparseBlock.h +++ b/Eigen/src/SparseCore/SparseBlock.h @@ -243,7 +243,7 @@ class SparseInnerVectorSet, Size> else if(m_outerSize.value()==0) return 0; else - return Map >(m_matrix.innerNonZeroPtr(), m_outerSize.value()).sum(); + return Map >(m_matrix.innerNonZeroPtr()+m_outerStart, m_outerSize.value()).sum(); } const Scalar& lastCoeff() const diff --git a/test/sparse_basic.cpp b/test/sparse_basic.cpp index 276db65aa..f36fcb36d 100644 --- a/test/sparse_basic.cpp +++ b/test/sparse_basic.cpp @@ -237,6 +237,26 @@ template void sparse_basic(const SparseMatrixType& re VERIFY_IS_APPROX(m2.innerVector(j0)+m2.innerVector(j1), refMat2.row(j0)+refMat2.row(j1)); else VERIFY_IS_APPROX(m2.innerVector(j0)+m2.innerVector(j1), refMat2.col(j0)+refMat2.col(j1)); + + SparseMatrixType m3(rows,rows); + m3.reserve(VectorXi::Constant(rows,rows/2)); + for(int j=0; j0) + VERIFY(j==m3.innerVector(j).lastCoeff()); + } + m3.makeCompressed(); + for(int j=0; j0) + VERIFY(j==m3.innerVector(j).lastCoeff()); + } + //m2.innerVector(j0) = 2*m2.innerVector(j1); //refMat2.col(j0) = 2*refMat2.col(j1); //VERIFY_IS_APPROX(m2, refMat2);