aboutsummaryrefslogtreecommitdiffhomepage
path: root/test/sparse_solvers.cpp
diff options
context:
space:
mode:
authorGravatar Gael Guennebaud <g.gael@free.fr>2009-05-18 17:55:50 +0200
committerGravatar Gael Guennebaud <g.gael@free.fr>2009-05-18 17:55:50 +0200
commite1867288672b7f8acff34fbdee59e7bb05ab9ead (patch)
treea484d649ed4188ed941ef22293ec9083ca939c83 /test/sparse_solvers.cpp
parente0832d5d9312d722f488704f1571b7b9780d463e (diff)
fix #1 : need to nest by value the affine part in homogeneous product
Diffstat (limited to 'test/sparse_solvers.cpp')
-rw-r--r--test/sparse_solvers.cpp10
1 files changed, 5 insertions, 5 deletions
diff --git a/test/sparse_solvers.cpp b/test/sparse_solvers.cpp
index 1f159801b..de308e018 100644
--- a/test/sparse_solvers.cpp
+++ b/test/sparse_solvers.cpp
@@ -72,12 +72,12 @@ template<typename Scalar> void sparse_solvers(int rows, int cols)
initSparse<Scalar>(density, refMat2, m2, ForceNonZeroDiag|MakeUpperTriangular, &zeroCoords, &nonzeroCoords);
VERIFY_IS_APPROX(refMat2.template marked<UpperTriangular>().solveTriangular(vec2),
m2.template triangular<UpperTriangular>().solve(vec3));
-
+
// TODO test row major
-
+
SparseMatrix<Scalar> matB(rows, rows);
DenseMatrix refMatB = DenseMatrix::Zero(rows, rows);
-
+
// lower - sparse
initSparse<Scalar>(density, refMat2, m2, ForceNonZeroDiag|MakeLowerTriangular);
initSparse<Scalar>(density, refMatB, matB);
@@ -91,7 +91,7 @@ template<typename Scalar> void sparse_solvers(int rows, int cols)
refMat2.template marked<UpperTriangular>().solveTriangularInPlace(refMatB);
m2.template triangular<UpperTriangular>().solveInPlace(matB);
VERIFY_IS_APPROX(matB, refMatB);
-
+
// test deprecated API
initSparse<Scalar>(density, refMat2, m2, ForceNonZeroDiag|MakeLowerTriangular, &zeroCoords, &nonzeroCoords);
VERIFY_IS_APPROX(refMat2.template marked<LowerTriangular>().solveTriangular(vec2),
@@ -122,7 +122,7 @@ template<typename Scalar> void sparse_solvers(int rows, int cols)
SparseLLT<SparseSelfAdjointMatrix,Cholmod>(m2).solveInPlace(x);
VERIFY(refX.isApprox(x,test_precision<Scalar>()) && "LLT: cholmod");
#endif
-
+
#ifdef EIGEN_TAUCS_SUPPORT
x = b;
SparseLLT<SparseSelfAdjointMatrix,Taucs>(m2,IncompleteFactorization).solveInPlace(x);