33 #ifndef EIGEN_REAL_SCHUR_LAPACKE_H 34 #define EIGEN_REAL_SCHUR_LAPACKE_H 40 #define EIGEN_LAPACKE_SCHUR_REAL(EIGTYPE, LAPACKE_TYPE, LAPACKE_PREFIX, LAPACKE_PREFIX_U, EIGCOLROW, LAPACKE_COLROW) \ 41 template<> template<typename InputType> inline \ 42 RealSchur<Matrix<EIGTYPE, Dynamic, Dynamic, EIGCOLROW> >& \ 43 RealSchur<Matrix<EIGTYPE, Dynamic, Dynamic, EIGCOLROW> >::compute(const EigenBase<InputType>& matrix, bool computeU) \ 45 eigen_assert(matrix.cols() == matrix.rows()); \ 47 lapack_int n = internal::convert_index<lapack_int>(matrix.cols()), sdim, info; \ 48 lapack_int matrix_order = LAPACKE_COLROW; \ 49 char jobvs, sort='N'; \ 50 LAPACK_##LAPACKE_PREFIX_U##_SELECT2 select = 0; \ 51 jobvs = (computeU) ? 'V' : 'N'; \ 52 m_matU.resize(n, n); \ 53 lapack_int ldvs = internal::convert_index<lapack_int>(m_matU.outerStride()); \ 55 lapack_int lda = internal::convert_index<lapack_int>(m_matT.outerStride()); \ 56 Matrix<EIGTYPE, Dynamic, Dynamic> wr, wi; \ 57 wr.resize(n, 1); wi.resize(n, 1); \ 58 info = LAPACKE_##LAPACKE_PREFIX##gees( matrix_order, jobvs, sort, select, n, (LAPACKE_TYPE*)m_matT.data(), lda, &sdim, (LAPACKE_TYPE*)wr.data(), (LAPACKE_TYPE*)wi.data(), (LAPACKE_TYPE*)m_matU.data(), ldvs ); \ 62 m_info = NoConvergence; \ 64 m_isInitialized = true; \ 65 m_matUisUptodate = computeU; \ 70 EIGEN_LAPACKE_SCHUR_REAL(
double,
double, d, D,
ColMajor, LAPACK_COL_MAJOR)
71 EIGEN_LAPACKE_SCHUR_REAL(
float,
float, s, S,
ColMajor, LAPACK_COL_MAJOR)
72 EIGEN_LAPACKE_SCHUR_REAL(
double,
double, d, D,
RowMajor, LAPACK_ROW_MAJOR)
73 EIGEN_LAPACKE_SCHUR_REAL(
float,
float, s, S,
RowMajor, LAPACK_ROW_MAJOR)
77 #endif // EIGEN_REAL_SCHUR_LAPACKE_H Definition: Constants.h:319
Namespace containing all symbols from the Eigen library.
Definition: Core:141
Definition: Constants.h:321