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) \
41template<> template<typename InputType> inline \
42RealSchur<Matrix<EIGTYPE, Dynamic, Dynamic, EIGCOLROW> >& \
43RealSchur<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; \
70EIGEN_LAPACKE_SCHUR_REAL(
double,
double, d, D, ColMajor, LAPACK_COL_MAJOR)
71EIGEN_LAPACKE_SCHUR_REAL(
float,
float, s, S, ColMajor, LAPACK_COL_MAJOR)
72EIGEN_LAPACKE_SCHUR_REAL(
double,
double, d, D, RowMajor, LAPACK_ROW_MAJOR)
73EIGEN_LAPACKE_SCHUR_REAL(
float,
float, s, S, RowMajor, LAPACK_ROW_MAJOR)
Namespace containing all symbols from the Eigen library.
Definition: Core:141