diff options
author | Gael Guennebaud <g.gael@free.fr> | 2019-11-13 21:16:53 +0100 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2019-11-13 21:16:53 +0100 |
commit | 8496f86f84a78e1d8b2424ad5048c5fdab37f58f (patch) | |
tree | 36590fde0fd271fcccc3c8a31a605dc2cdc44939 /test/permutationmatrices.cpp | |
parent | 002e5b6db6b89d3e2543a17f031121c7d97a2f33 (diff) |
Enable CompleteOrthogonalDecomposition::pseudoInverse with non-square fixed-size matrices.
Diffstat (limited to 'test/permutationmatrices.cpp')
-rw-r--r-- | test/permutationmatrices.cpp | 6 |
1 files changed, 5 insertions, 1 deletions
diff --git a/test/permutationmatrices.cpp b/test/permutationmatrices.cpp index a26b4ad45..d4b68b2d4 100644 --- a/test/permutationmatrices.cpp +++ b/test/permutationmatrices.cpp @@ -54,7 +54,11 @@ template<typename MatrixType> void permutationmatrices(const MatrixType& m) m_permuted = m_original; VERIFY_EVALUATION_COUNT(m_permuted = lp * m_permuted * rp, 1); VERIFY_IS_APPROX(m_permuted, lm*m_original*rm); - + + LeftPermutationType lpi; + lpi = lp.inverse(); + VERIFY_IS_APPROX(lpi*m_permuted,lp.inverse()*m_permuted); + VERIFY_IS_APPROX(lp.inverse()*m_permuted*rp.inverse(), m_original); VERIFY_IS_APPROX(lv.asPermutation().inverse()*m_permuted*rv.asPermutation().inverse(), m_original); VERIFY_IS_APPROX(MapLeftPerm(lv.data(),lv.size()).inverse()*m_permuted*MapRightPerm(rv.data(),rv.size()).inverse(), m_original); |