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

implement two levels of blocking in PartialLU => high speedup

Browse Source
This commit is contained in:
Gael Guennebaud
2009-08-04 11:28:02 +02:00
parent 912da9fade
commit 4bec101470
2 changed files with 142 additions and 84 deletions
Download Patch File Download Diff File Expand all files Collapse all files

7
test/inverse.cpp
View File

@@ -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