GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-04-21 17:19:36 +08:00
Code Issues Packages Projects Releases Wiki Activity

forward port slight changes in the 4x4 inverse test

Browse Source
This commit is contained in:
Benoit Jacob 2009-11-25 09:08:03 -05:00
parent 684d76eba3
commit 1e2bcba5e4
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
test/prec_inverse_4x4.cpp
View File

@ -79,7 +79,7 @@ template<typename MatrixType> void inverse_general_4x4(int repeat)
do { do {
m = MatrixType::Random(); m = MatrixType::Random();
absdet = ei_abs(m.determinant()); absdet = ei_abs(m.determinant());
} while(absdet < 2*epsilon<Scalar>() ); } while(absdet < 10 * epsilon<Scalar>());
MatrixType inv = m.inverse(); MatrixType inv = m.inverse();
double error = double( (m*inv-MatrixType::Identity()).norm() * absdet / epsilon<Scalar>() ); double error = double( (m*inv-MatrixType::Identity()).norm() * absdet / epsilon<Scalar>() );
error_sum += error; error_sum += error;
@ -90,7 +90,7 @@ template<typename MatrixType> void inverse_general_4x4(int repeat)
EIGEN_DEBUG_VAR(error_avg); EIGEN_DEBUG_VAR(error_avg);
EIGEN_DEBUG_VAR(error_max); EIGEN_DEBUG_VAR(error_max);
VERIFY(error_avg < (NumTraits<Scalar>::IsComplex ? 8.4 : 1.4) ); VERIFY(error_avg < (NumTraits<Scalar>::IsComplex ? 8.4 : 1.4) );
VERIFY(error_max < (NumTraits<Scalar>::IsComplex ? 150.0 : 75.) ); VERIFY(error_max < (NumTraits<Scalar>::IsComplex ? 160.0 : 75.) );
} }
void test_prec_inverse_4x4() void test_prec_inverse_4x4()