diff options
author | Gael Guennebaud <g.gael@free.fr> | 2011-12-01 18:17:19 +0100 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2011-12-01 18:17:19 +0100 |
commit | b85bcd91bf815fee8d42ac6f2add45e246f0ed51 (patch) | |
tree | 1ada8572d606a3382366e5910d13a702ac0b9a81 /test/cholesky.cpp | |
parent | 7aaae9d6dfbb57390a4ab9370a3553e7e1501fdd (diff) |
remove GSL dependency in the unit tests
Diffstat (limited to 'test/cholesky.cpp')
-rw-r--r-- | test/cholesky.cpp | 32 |
1 files changed, 0 insertions, 32 deletions
diff --git a/test/cholesky.cpp b/test/cholesky.cpp index 2b646e329..66eae547f 100644 --- a/test/cholesky.cpp +++ b/test/cholesky.cpp @@ -41,10 +41,6 @@ static int nb_temporaries; VERIFY( (#XPR) && nb_temporaries==N ); \ } -#ifdef HAS_GSL -#include "gsl_helper.h" -#endif - template<typename MatrixType> void cholesky(const MatrixType& m) { typedef typename MatrixType::Index Index; @@ -77,34 +73,6 @@ template<typename MatrixType> void cholesky(const MatrixType& m) // FIXME: currently that fails !! //symm.template part<StrictlyLower>().setZero(); - #ifdef HAS_GSL -// if (internal::is_same<RealScalar,double>::value) -// { -// typedef GslTraits<Scalar> Gsl; -// typename Gsl::Matrix gMatA=0, gSymm=0; -// typename Gsl::Vector gVecB=0, gVecX=0; -// convert<MatrixType>(symm, gSymm); -// convert<MatrixType>(symm, gMatA); -// convert<VectorType>(vecB, gVecB); -// convert<VectorType>(vecB, gVecX); -// Gsl::cholesky(gMatA); -// Gsl::cholesky_solve(gMatA, gVecB, gVecX); -// VectorType vecX(rows), _vecX, _vecB; -// convert(gVecX, _vecX); -// symm.llt().solve(vecB, &vecX); -// Gsl::prod(gSymm, gVecX, gVecB); -// convert(gVecB, _vecB); -// // test gsl itself ! -// VERIFY_IS_APPROX(vecB, _vecB); -// VERIFY_IS_APPROX(vecX, _vecX); -// -// Gsl::free(gMatA); -// Gsl::free(gSymm); -// Gsl::free(gVecB); -// Gsl::free(gVecX); -// } - #endif - { LLT<SquareMatrixType,Lower> chollo(symmLo); VERIFY_IS_APPROX(symm, chollo.reconstructedMatrix()); |