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