devtools/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2026-04-10 11:34:33 +08:00
Code Issues Packages Projects Releases Wiki Activity

Fix more gcc compiler warnings / sort-of bugs

Browse Source
This commit is contained in:
Charles Schlosser
2023-07-14 21:12:45 +00:00
parent 21cd3fe209
commit 81fe2d424f
5 changed files with 24 additions and 16 deletions
Download Patch File Download Diff File Expand all files Collapse all files

25
test/triangular.cpp
View File

@@ -55,6 +55,7 @@ template<typename MatrixType> void triangular_square(const MatrixType& m)
r1(rows, cols),
r2(rows, cols);
VectorType v2 = VectorType::Random(rows);
VectorType v3 = VectorType::Zero(rows);
MatrixType m1up = m1.template triangularView<Upper>();
MatrixType m2up = m2.template triangularView<Upper>();
@@ -96,23 +97,31 @@ template<typename MatrixType> void triangular_square(const MatrixType& m)
Transpose<MatrixType> trm4(m4);
// test back and forward substitution with a vector as the rhs
m3 = m1.template triangularView<Upper>();
VERIFY(v2.isApprox(m3.adjoint() * (m1.adjoint().template triangularView<Lower>().solve(v2)), largerEps));
v3 = m3.adjoint() * (m1.adjoint().template triangularView<Lower>().solve(v2));
VERIFY(v2.isApprox(v3, largerEps));
m3 = m1.template triangularView<Lower>();
VERIFY(v2.isApprox(m3.transpose() * (m1.transpose().template triangularView<Upper>().solve(v2)), largerEps));
v3 = m3.transpose() * (m1.transpose().template triangularView<Upper>().solve(v2));
VERIFY(v2.isApprox(v3, largerEps));
m3 = m1.template triangularView<Upper>();
VERIFY(v2.isApprox(m3 * (m1.template triangularView<Upper>().solve(v2)), largerEps));
v3 = m3 * (m1.template triangularView<Upper>().solve(v2));
VERIFY(v2.isApprox(v3, largerEps));
m3 = m1.template triangularView<Lower>();
VERIFY(v2.isApprox(m3.conjugate() * (m1.conjugate().template triangularView<Lower>().solve(v2)), largerEps));
v3 = m3.conjugate() * (m1.conjugate().template triangularView<Lower>().solve(v2));
VERIFY(v2.isApprox(v3, largerEps));
// test back and forward substitution with a matrix as the rhs
m3 = m1.template triangularView<Upper>();
VERIFY(m2.isApprox(m3.adjoint() * (m1.adjoint().template triangularView<Lower>().solve(m2)), largerEps));
m4 = m3.adjoint() * (m1.adjoint().template triangularView<Lower>().solve(m2));
VERIFY(m2.isApprox(m4, largerEps));
m3 = m1.template triangularView<Lower>();
VERIFY(m2.isApprox(m3.transpose() * (m1.transpose().template triangularView<Upper>().solve(m2)), largerEps));
m4 = m3.transpose() * (m1.transpose().template triangularView<Upper>().solve(m2));
VERIFY(m2.isApprox(m4, largerEps));
m3 = m1.template triangularView<Upper>();
VERIFY(m2.isApprox(m3 * (m1.template triangularView<Upper>().solve(m2)), largerEps));
m4 = m3 * (m1.template triangularView<Upper>().solve(m2));
VERIFY(m2.isApprox(m4, largerEps));
m3 = m1.template triangularView<Lower>();
VERIFY(m2.isApprox(m3.conjugate() * (m1.conjugate().template triangularView<Lower>().solve(m2)), largerEps));
m4 = m3.conjugate() * (m1.conjugate().template triangularView<Lower>().solve(m2));
VERIFY(m2.isApprox(m4, largerEps));
// check M * inv(L) using in place API
m4 = m3;