From 9bd902ed9cd55198683f4248cd685b05851acaee Mon Sep 17 00:00:00 2001 From: Gael Guennebaud Date: Sun, 4 Dec 2011 14:57:43 +0100 Subject: fix bug #341: trisove on MappedSparseMatrix --- test/sparse_solvers.cpp | 7 +++++++ 1 file changed, 7 insertions(+) (limited to 'test/sparse_solvers.cpp') 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 void sparse_solvers(int rows, int cols) m2.template triangularView().solve(vec3)); VERIFY_IS_APPROX(refMat2.conjugate().template triangularView().solve(vec2), m2.conjugate().template triangularView().solve(vec3)); + { + SparseMatrix cm2(m2); + //Index rows, Index cols, Index nnz, Index* outerIndexPtr, Index* innerIndexPtr, Scalar* valuePtr + MappedSparseMatrix mm2(rows, cols, cm2.nonZeros(), cm2.outerIndexPtr(), cm2.innerIndexPtr(), cm2.valuePtr()); + VERIFY_IS_APPROX(refMat2.conjugate().template triangularView().solve(vec2), + mm2.conjugate().template triangularView().solve(vec3)); + } // lower - transpose initSparse(density, refMat2, m2, ForceNonZeroDiag|MakeLowerTriangular, &zeroCoords, &nonzeroCoords); -- cgit v1.2.3