34 #ifndef EIGEN_QR_MKL_H 35 #define EIGEN_QR_MKL_H 37 #include "Eigen/src/Core/util/MKL_support.h" 45 #define EIGEN_MKL_QR_NOPIV(EIGTYPE, MKLTYPE, MKLPREFIX) \ 46 template<typename MatrixQR, typename HCoeffs> \ 47 void householder_qr_inplace_blocked(MatrixQR& mat, HCoeffs& hCoeffs, \ 48 typename MatrixQR::Index maxBlockSize=32, \ 49 EIGTYPE* tempData = 0) \ 51 lapack_int m = mat.rows(); \ 52 lapack_int n = mat.cols(); \ 53 lapack_int lda = mat.outerStride(); \ 54 lapack_int matrix_order = (MatrixQR::IsRowMajor) ? LAPACK_ROW_MAJOR : LAPACK_COL_MAJOR; \ 55 LAPACKE_##MKLPREFIX##geqrf( matrix_order, m, n, (MKLTYPE*)mat.data(), lda, (MKLTYPE*)hCoeffs.data()); \ 56 hCoeffs.adjointInPlace(); \ 60 EIGEN_MKL_QR_NOPIV(
double,
double, d)
61 EIGEN_MKL_QR_NOPIV(
float,
float, s)
62 EIGEN_MKL_QR_NOPIV(dcomplex, MKL_Complex16, z)
63 EIGEN_MKL_QR_NOPIV(scomplex, MKL_Complex8, c)
69 #endif // EIGEN_QR_MKL_H
Definition: Eigen_Colamd.h:54