33 #ifndef EIGEN_COMPLEX_SCHUR_MKL_H 34 #define EIGEN_COMPLEX_SCHUR_MKL_H 36 #include "Eigen/src/Core/util/MKL_support.h" 42 #define EIGEN_MKL_SCHUR_COMPLEX(EIGTYPE, MKLTYPE, MKLPREFIX, MKLPREFIX_U, EIGCOLROW, MKLCOLROW) \ 43 template<> template<typename InputType> inline \ 44 ComplexSchur<Matrix<EIGTYPE, Dynamic, Dynamic, EIGCOLROW> >& \ 45 ComplexSchur<Matrix<EIGTYPE, Dynamic, Dynamic, EIGCOLROW> >::compute(const EigenBase<InputType>& matrix, bool computeU) \ 47 typedef Matrix<EIGTYPE, Dynamic, Dynamic, EIGCOLROW> MatrixType; \ 48 typedef MatrixType::RealScalar RealScalar; \ 49 typedef std::complex<RealScalar> ComplexScalar; \ 51 eigen_assert(matrix.cols() == matrix.rows()); \ 53 m_matUisUptodate = false; \ 54 if(matrix.cols() == 1) \ 56 m_matT = matrix.derived().template cast<ComplexScalar>(); \ 57 if(computeU) m_matU = ComplexMatrixType::Identity(1,1); \ 59 m_isInitialized = true; \ 60 m_matUisUptodate = computeU; \ 63 lapack_int n = matrix.cols(), sdim, info; \ 64 lapack_int matrix_order = MKLCOLROW; \ 65 char jobvs, sort='N'; \ 66 LAPACK_##MKLPREFIX_U##_SELECT1 select = 0; \ 67 jobvs = (computeU) ? 'V' : 'N'; \ 68 m_matU.resize(n, n); \ 69 lapack_int ldvs = m_matU.outerStride(); \ 71 lapack_int lda = m_matT.outerStride(); \ 72 Matrix<EIGTYPE, Dynamic, Dynamic> w; \ 74 info = LAPACKE_##MKLPREFIX##gees( matrix_order, jobvs, sort, select, n, (MKLTYPE*)m_matT.data(), lda, &sdim, (MKLTYPE*)w.data(), (MKLTYPE*)m_matU.data(), ldvs ); \ 78 m_info = NoConvergence; \ 80 m_isInitialized = true; \ 81 m_matUisUptodate = computeU; \ 86 EIGEN_MKL_SCHUR_COMPLEX(dcomplex, MKL_Complex16, z, Z,
ColMajor, LAPACK_COL_MAJOR)
87 EIGEN_MKL_SCHUR_COMPLEX(scomplex, MKL_Complex8, c, C,
ColMajor, LAPACK_COL_MAJOR)
88 EIGEN_MKL_SCHUR_COMPLEX(dcomplex, MKL_Complex16, z, Z,
RowMajor, LAPACK_ROW_MAJOR)
89 EIGEN_MKL_SCHUR_COMPLEX(scomplex, MKL_Complex8, c, C, RowMajor, LAPACK_ROW_MAJOR)
93 #endif // EIGEN_COMPLEX_SCHUR_MKL_H
Definition: Constants.h:320
Definition: Constants.h:322