diff options
| author |  Chen-Pang He <jdh8@ms63.hinet.net> | 2012-09-16 14:42:42 +0800 |
|---|---|---|
| committer | Chen-Pang He <jdh8@ms63.hinet.net> | 2012-09-16 14:42:42 +0800 |
| commit | d5d99dd1f001b6bb05d530ab40e455a69703dcc6 (patch) | |
| tree | 235087811cc139adaba1c57d54a81b3621fa4e8c /unsupported | |
| parent | 04f315d692d4b2b01635981d34b44743ba63571c (diff) | |
Optimize matrix functions: m_fT is triangular and trmm is faster than gemm
Diffstat (limited to 'unsupported')
| -rw-r--r-- | unsupported/Eigen/src/MatrixFunctions/MatrixFunction.h | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/unsupported/Eigen/src/MatrixFunctions/MatrixFunction.h b/unsupported/Eigen/src/MatrixFunctions/MatrixFunction.h index c57ca87ed..e87a28f6c 100644 --- a/unsupported/Eigen/src/MatrixFunctions/MatrixFunction.h +++ b/unsupported/Eigen/src/MatrixFunctions/MatrixFunction.h @@ -209,7 +209,7 @@ void MatrixFunction<MatrixType,AtomicType,1>::compute(ResultType& result) permuteSchur(); computeBlockAtomic(); computeOffDiagonal(); - result = m_U * m_fT * m_U.adjoint(); + result = m_U * (m_fT.template triangularView<Upper>() * m_U.adjoint()); } /** \brief Store the Schur decomposition of #m_A in #m_T and #m_U */ |