Optimize matrix functions: m_fT is triangular and trmm is faster than gemm

This commit is contained in:
Chen-Pang He 2012-09-16 14:42:42 +08:00
parent 04f315d692
commit d5d99dd1f0

View File

@ -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 */