GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-06-04 18:54:00 +08:00
Code Issues Packages Projects Releases Wiki Activity

fix dirty triangular unit test

Browse Source
This commit is contained in:
Gael Guennebaud 2010-01-04 19:02:43 +01:00
parent 8eab9fb87e
commit b6898996d4
1 changed files with 2 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
test/triangular.cpp
View File

@ -110,24 +110,14 @@ template<typename MatrixType> void triangular_square(const MatrixType& m)
// check M * inv(L) using in place API // check M * inv(L) using in place API
m4 = m3; m4 = m3;
<<<<<<< local
m3.transpose().template triangularView<Eigen::UpperTriangular>().solveInPlace(trm4); m3.transpose().template triangularView<Eigen::UpperTriangular>().solveInPlace(trm4);
VERIFY(m4.cwiseAbs().isIdentity(test_precision<RealScalar>())); VERIFY(m4.cwiseAbs().isIdentity(test_precision<RealScalar>()));
=======
m3.transpose().template triangularView<UpperTriangular>().solveInPlace(trm4);
VERIFY(m4.cwise().abs().isIdentity(test_precision<RealScalar>()));
>>>>>>> other
// check M * inv(U) using in place API // check M * inv(U) using in place API
m3 = m1.template triangularView<UpperTriangular>(); m3 = m1.template triangularView<UpperTriangular>();
m4 = m3; m4 = m3;
<<<<<<< local
m3.transpose().template triangularView<Eigen::LowerTriangular>().solveInPlace(trm4); m3.transpose().template triangularView<Eigen::LowerTriangular>().solveInPlace(trm4);
VERIFY(m4.cwiseAbs().isIdentity(test_precision<RealScalar>())); VERIFY(m4.cwiseAbs().isIdentity(test_precision<RealScalar>()));
=======
m3.transpose().template triangularView<LowerTriangular>().solveInPlace(trm4);
VERIFY(m4.cwise().abs().isIdentity(test_precision<RealScalar>()));
>>>>>>> other
// check solve with unit diagonal // check solve with unit diagonal
m3 = m1.template triangularView<UnitUpperTriangular>(); m3 = m1.template triangularView<UnitUpperTriangular>();
@ -220,7 +210,7 @@ template<typename MatrixType> void triangular_rect(const MatrixType& m)
VERIFY(m2.diagonal().isMuchSmallerThan(RealScalar(1))); VERIFY(m2.diagonal().isMuchSmallerThan(RealScalar(1)));
m2 = m1.template triangularView<UnitUpperTriangular>(); m2 = m1.template triangularView<UnitUpperTriangular>();
VERIFY(m2.isUpperTriangular()); VERIFY(m2.isUpperTriangular());
m2.diagonal().cwise() -= Scalar(1); m2.diagonal().array() -= Scalar(1);
VERIFY(m2.diagonal().isMuchSmallerThan(RealScalar(1))); VERIFY(m2.diagonal().isMuchSmallerThan(RealScalar(1)));
m2 = m1.template triangularView<LowerTriangular>(); m2 = m1.template triangularView<LowerTriangular>();
VERIFY(m2.isLowerTriangular()); VERIFY(m2.isLowerTriangular());
@ -230,7 +220,7 @@ template<typename MatrixType> void triangular_rect(const MatrixType& m)
VERIFY(m2.diagonal().isMuchSmallerThan(RealScalar(1))); VERIFY(m2.diagonal().isMuchSmallerThan(RealScalar(1)));
m2 = m1.template triangularView<UnitLowerTriangular>(); m2 = m1.template triangularView<UnitLowerTriangular>();
VERIFY(m2.isLowerTriangular()); VERIFY(m2.isLowerTriangular());
m2.diagonal().cwise() -= Scalar(1); m2.diagonal().array() -= Scalar(1);
VERIFY(m2.diagonal().isMuchSmallerThan(RealScalar(1))); VERIFY(m2.diagonal().isMuchSmallerThan(RealScalar(1)));
// test swap // test swap
m1.setOnes(); m1.setOnes();