fix bug #341: trisove on MappedSparseMatrix

1 files changed, 7 insertions, 0 deletions

diff --git a/test/sparse_solvers.cpp b/test/sparse_solvers.cpp
index 3d1140bd1..a3d79b0ea 100644
--- a/test/sparse_solvers.cpp
+++ b/test/sparse_solvers.cpp
@@ -74,6 +74,13 @@ template<typename Scalar> void sparse_solvers(int rows, int cols)
                      m2.template triangularView<Upper>().solve(vec3));
     VERIFY_IS_APPROX(refMat2.conjugate().template triangularView<Upper>().solve(vec2),
                      m2.conjugate().template triangularView<Upper>().solve(vec3));
+    {
+      SparseMatrix<Scalar> cm2(m2);
+      //Index rows, Index cols, Index nnz, Index* outerIndexPtr, Index* innerIndexPtr, Scalar* valuePtr
+      MappedSparseMatrix<Scalar> mm2(rows, cols, cm2.nonZeros(), cm2.outerIndexPtr(), cm2.innerIndexPtr(), cm2.valuePtr());
+      VERIFY_IS_APPROX(refMat2.conjugate().template triangularView<Upper>().solve(vec2),
+                       mm2.conjugate().template triangularView<Upper>().solve(vec3));
+    }
 
     // lower - transpose
     initSparse<Scalar>(density, refMat2, m2, ForceNonZeroDiag|MakeLowerTriangular, &zeroCoords, &nonzeroCoords);