aboutsummaryrefslogtreecommitdiffhomepage
path: root/test/cholesky.cpp
diff options
context:
space:
mode:
authorGravatar Gael Guennebaud <g.gael@free.fr>2011-12-01 18:17:19 +0100
committerGravatar Gael Guennebaud <g.gael@free.fr>2011-12-01 18:17:19 +0100
commitb85bcd91bf815fee8d42ac6f2add45e246f0ed51 (patch)
tree1ada8572d606a3382366e5910d13a702ac0b9a81 /test/cholesky.cpp
parent7aaae9d6dfbb57390a4ab9370a3553e7e1501fdd (diff)
remove GSL dependency in the unit tests
Diffstat (limited to 'test/cholesky.cpp')
-rw-r--r--test/cholesky.cpp32
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());