diff options
author | Gael Guennebaud <g.gael@free.fr> | 2009-08-04 11:28:02 +0200 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2009-08-04 11:28:02 +0200 |
commit | 4bec101470091aae27e7469025c80e31b889f566 (patch) | |
tree | 2f611c04043575272e7f1a24cda4e1199b7196c4 /test/inverse.cpp | |
parent | 912da9fade7a02086747d086e0634cab0d6ff4b6 (diff) |
implement two levels of blocking in PartialLU => high speedup
Diffstat (limited to 'test/inverse.cpp')
-rw-r--r-- | test/inverse.cpp | 7 |
1 files changed, 5 insertions, 2 deletions
diff --git a/test/inverse.cpp b/test/inverse.cpp index b4eef73b6..65dfbc73e 100644 --- a/test/inverse.cpp +++ b/test/inverse.cpp @@ -81,13 +81,16 @@ template<typename MatrixType> void inverse(const MatrixType& m) void test_inverse() { + int s; for(int i = 0; i < g_repeat; i++) { CALL_SUBTEST( inverse(Matrix<double,1,1>()) ); CALL_SUBTEST( inverse(Matrix2d()) ); CALL_SUBTEST( inverse(Matrix3f()) ); CALL_SUBTEST( inverse(Matrix4f()) ); - CALL_SUBTEST( inverse(MatrixXf(72,72)) ); - CALL_SUBTEST( inverse(MatrixXcd(56,56)) ); + s = ei_random<int>(50,320); + CALL_SUBTEST( inverse(MatrixXf(s,s)) ); + s = ei_random<int>(25,100); + CALL_SUBTEST( inverse(MatrixXcd(s,s)) ); } // test some tricky cases for 4x4 matrices |