diff options
author | Gael Guennebaud <g.gael@free.fr> | 2010-02-25 16:30:58 +0100 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2010-02-25 16:30:58 +0100 |
commit | 959a1b5d6335833e9ad49a088502705bb6967ff5 (patch) | |
tree | 0b41aa2f93a5b009b3699809ccfd4000582c92bf /Eigen/src/LU/PartialPivLU.h | |
parent | d9ca0c0d3643f4b777de686a2c0cddde075aa063 (diff) |
detect and implement inplace permutations
Diffstat (limited to 'Eigen/src/LU/PartialPivLU.h')
-rw-r--r-- | Eigen/src/LU/PartialPivLU.h | 5 |
1 files changed, 2 insertions, 3 deletions
diff --git a/Eigen/src/LU/PartialPivLU.h b/Eigen/src/LU/PartialPivLU.h index fcffc2458..ad0d6b810 100644 --- a/Eigen/src/LU/PartialPivLU.h +++ b/Eigen/src/LU/PartialPivLU.h @@ -412,10 +412,9 @@ MatrixType PartialPivLU<MatrixType>::reconstructedMatrix() const // LU MatrixType res = m_lu.template triangularView<UnitLower>().toDenseMatrix() * m_lu.template triangularView<Upper>(); - + // P^{-1}(LU) - // FIXME implement inplace permutation - res = (m_p.inverse() * res).eval(); + res = m_p.inverse() * res; return res; } |