From 4d54e3dd332b7f07eccb428df2e1e3074abc3aeb Mon Sep 17 00:00:00 2001 From: Gael Guennebaud Date: Mon, 25 Jul 2016 17:55:07 +0200 Subject: bug #173: remove dependency to MKL for LAPACKe backend. --- Eigen/src/Cholesky/LLT_MKL.h | 16 ++++++++-------- 1 file changed, 8 insertions(+), 8 deletions(-) (limited to 'Eigen/src/Cholesky') diff --git a/Eigen/src/Cholesky/LLT_MKL.h b/Eigen/src/Cholesky/LLT_MKL.h index f5be3b240..f0f8600e4 100644 --- a/Eigen/src/Cholesky/LLT_MKL.h +++ b/Eigen/src/Cholesky/LLT_MKL.h @@ -40,10 +40,10 @@ namespace Eigen { namespace internal { -template struct mkl_llt; +template struct lapacke_llt; -#define EIGEN_MKL_LLT(EIGTYPE, MKLTYPE, MKLPREFIX) \ -template<> struct mkl_llt \ +#define EIGEN_MKL_LLT(EIGTYPE, BLASTYPE, MKLPREFIX) \ +template<> struct lapacke_llt \ { \ template \ static inline Index potrf(MatrixType& m, char uplo) \ @@ -59,7 +59,7 @@ template<> struct mkl_llt \ a = &(m.coeffRef(0,0)); \ lda = convert_index(m.outerStride()); \ \ - info = LAPACKE_##MKLPREFIX##potrf( matrix_order, uplo, size, (MKLTYPE*)a, lda ); \ + info = LAPACKE_##MKLPREFIX##potrf( matrix_order, uplo, size, (BLASTYPE*)a, lda ); \ info = (info==0) ? -1 : info>0 ? info-1 : size; \ return info; \ } \ @@ -69,7 +69,7 @@ template<> struct llt_inplace \ template \ static Index blocked(MatrixType& m) \ { \ - return mkl_llt::potrf(m, 'L'); \ + return lapacke_llt::potrf(m, 'L'); \ } \ template \ static Index rankUpdate(MatrixType& mat, const VectorType& vec, const typename MatrixType::RealScalar& sigma) \ @@ -80,7 +80,7 @@ template<> struct llt_inplace \ template \ static Index blocked(MatrixType& m) \ { \ - return mkl_llt::potrf(m, 'U'); \ + return lapacke_llt::potrf(m, 'U'); \ } \ template \ static Index rankUpdate(MatrixType& mat, const VectorType& vec, const typename MatrixType::RealScalar& sigma) \ @@ -92,8 +92,8 @@ template<> struct llt_inplace \ EIGEN_MKL_LLT(double, double, d) EIGEN_MKL_LLT(float, float, s) -EIGEN_MKL_LLT(dcomplex, MKL_Complex16, z) -EIGEN_MKL_LLT(scomplex, MKL_Complex8, c) +EIGEN_MKL_LLT(dcomplex, lapack_complex_double, z) +EIGEN_MKL_LLT(scomplex, lapack_complex_float, c) } // end namespace internal -- cgit v1.2.3