From 13f31b8daf4d67d6310c567de36a34afa7c6e18f Mon Sep 17 00:00:00 2001 From: Benoit Jacob Date: Tue, 20 Oct 2009 00:36:07 -0400 Subject: * make PartialLU avoid to generate inf/nan when given a singular matrix (result undefined, but at least it won't take forever on intel 387) * add lots of comments, especially to LU.h * fix stuff I had broken in Inverse.h * split inverse test --- test/inverse.cpp | 14 ++++++++------ 1 file changed, 8 insertions(+), 6 deletions(-) (limited to 'test/inverse.cpp') diff --git a/test/inverse.cpp b/test/inverse.cpp index 65dfbc73e..6fc65786c 100644 --- a/test/inverse.cpp +++ b/test/inverse.cpp @@ -83,19 +83,21 @@ void test_inverse() { int s; for(int i = 0; i < g_repeat; i++) { - CALL_SUBTEST( inverse(Matrix()) ); - CALL_SUBTEST( inverse(Matrix2d()) ); - CALL_SUBTEST( inverse(Matrix3f()) ); - CALL_SUBTEST( inverse(Matrix4f()) ); + CALL_SUBTEST1( inverse(Matrix()) ); + CALL_SUBTEST2( inverse(Matrix2d()) ); + CALL_SUBTEST3( inverse(Matrix3f()) ); + CALL_SUBTEST4( inverse(Matrix4f()) ); s = ei_random(50,320); - CALL_SUBTEST( inverse(MatrixXf(s,s)) ); + CALL_SUBTEST5( inverse(MatrixXf(s,s)) ); s = ei_random(25,100); - CALL_SUBTEST( inverse(MatrixXcd(s,s)) ); + CALL_SUBTEST6( inverse(MatrixXcd(s,s)) ); } +#ifdef EIGEN_TEST_PART_4 // test some tricky cases for 4x4 matrices VERIFY_IS_APPROX((Matrix4f() << 0,0,1,0, 1,0,0,0, 0,1,0,0, 0,0,0,1).finished().inverse(), (Matrix4f() << 0,1,0,0, 0,0,1,0, 1,0,0,0, 0,0,0,1).finished()); VERIFY_IS_APPROX((Matrix4f() << 1,0,0,0, 0,0,1,0, 0,0,0,1, 0,1,0,0).finished().inverse(), (Matrix4f() << 1,0,0,0, 0,0,0,1, 0,1,0,0, 0,0,1,0).finished()); +#endif } -- cgit v1.2.3