From b85bcd91bf815fee8d42ac6f2add45e246f0ed51 Mon Sep 17 00:00:00 2001
From: Gael Guennebaud <g.gael@free.fr>
Date: Thu, 1 Dec 2011 18:17:19 +0100
Subject: remove GSL dependency in the unit tests

---
 test/cholesky.cpp | 32 --------------------------------
 1 file changed, 32 deletions(-)

(limited to 'test/cholesky.cpp')

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());
-- 
cgit v1.2.3