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

added a tough test to check the determinant that currently fails

Browse Source
This commit is contained in:
Gael Guennebaud
2008-04-25 23:13:20 +00:00
parent 6f2c72fb53
commit 173e582e3c
5 changed files with 85 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

7
test/linearstructure.cpp
View File

@@ -89,6 +89,12 @@ template<typename MatrixType> void linearStructure(const MatrixType& m)
VERIFY_IS_APPROX((m1*s1)(r,c), (m1(r,c))*s1);
if(NumTraits<Scalar>::HasFloatingPoint)
VERIFY_IS_APPROX((m1/s1)(r,c), (m1(r,c))/s1);
// use .block to disable vectorization and compare to the vectorized version
VERIFY_IS_APPROX(m1+m1.block(0,0,rows,cols), m1+m1);
VERIFY_IS_APPROX(m1.cwiseProduct(m1.block(0,0,rows,cols)), m1.cwiseProduct(m1));
VERIFY_IS_APPROX(m1 - m1.block(0,0,rows,cols), m1 - m1);
VERIFY_IS_APPROX(m1.block(0,0,rows,cols) * s1, m1 * s1);
}
void EigenTest::testLinearStructure()
@@ -97,6 +103,7 @@ void EigenTest::testLinearStructure()
linearStructure(Matrix<float, 1, 1>());
linearStructure(Matrix4d());
linearStructure(MatrixXcf(3, 3));
linearStructure(MatrixXf(8, 12));
linearStructure(MatrixXi(8, 12));
linearStructure(MatrixXcd(20, 20));
}