diff options
| author | torset | 2009-02-27 13:59:12 +0000 |
|---|---|---|
| committer | torset | 2009-02-27 13:59:12 +0000 |
| commit | 6e7f75a3a916ae7b7b6145c0badad7b649265c0f (patch) | |
| tree | d6150181ae1536cf65d4a1a873345a7e7e8e8f50 /src/matrixOperations/expm/testMatrixExponential.c | |
| parent | 05fb6a62c83ae7cea37aa2795e61697b3560f135 (diff) | |
| download | scilab2c-6e7f75a3a916ae7b7b6145c0badad7b649265c0f.tar.gz scilab2c-6e7f75a3a916ae7b7b6145c0badad7b649265c0f.tar.bz2 scilab2c-6e7f75a3a916ae7b7b6145c0badad7b649265c0f.zip | |
change precision for tests : put it to 3e-16 for double and 3e-6 for float
Diffstat (limited to 'src/matrixOperations/expm/testMatrixExponential.c')
| -rw-r--r-- | src/matrixOperations/expm/testMatrixExponential.c | 22 |
1 files changed, 11 insertions, 11 deletions
diff --git a/src/matrixOperations/expm/testMatrixExponential.c b/src/matrixOperations/expm/testMatrixExponential.c index 13bf82d5..72586c53 100644 --- a/src/matrixOperations/expm/testMatrixExponential.c +++ b/src/matrixOperations/expm/testMatrixExponential.c @@ -421,14 +421,14 @@ static void sexpmaTest (void ) { for ( i = 0 ; i < LEADDIM*LEADDIM ; i++ ) { printf ( "\t\t %d out : %e\tresult : %e\tassert : %e \n" , i , out[i] , result[i] , fabs ( out[i] - result[i] ) / fabs( out[i]) ) ; - assert ( fabs ( out[i] - result[i] ) / fabs( out[i]) < 1e-5 ) ; + assert ( fabs ( out[i] - result[i] ) / fabs( out[i]) < 3e-6 ) ; } sexpma(mon_test,out_mon_test,3); for ( i = 0 ; i < 9 ; i++ ) { printf ( "\t\t %d out : %e\tresult : %e\tassert : %e \n" , i , out[i] , result[i] , fabs ( out[i] - result[i] ) / fabs( out[i]) ) ; - assert ( fabs ( out_mon_test[i] - result_mon_test[i] ) / fabs( out_mon_test[i]) < 1e-5 ) ; + assert ( fabs ( out_mon_test[i] - result_mon_test[i] ) / fabs( out_mon_test[i]) < 3e-6 ) ; } } @@ -450,7 +450,7 @@ static void dexpmaTest (void ) { for ( i = 0 ; i < LEADDIM*LEADDIM ; i++ ) { printf ( "\t\t %d out : %e\tresult : %e\tassert : %e \n" , i , out[i] , result[i] , fabs ( out[i] - result[i] ) / fabs( out[i]) ) ; - assert ( fabs ( out[i] - result[i] ) / fabs( out[i]) < 1e-14 ) ; + assert ( fabs ( out[i] - result[i] ) / fabs( out[i]) < 3e-16 ) ; } } @@ -483,16 +483,16 @@ static void cexpmaTest ( void) { fabs( creals(out[i]) - creals (Result[i]) ) / fabs (creals (out[i])) , fabs( cimags(out[i]) - cimags (Result[i]) ) / fabs (cimags (out[i]))); - if ( creals(out[i]) < 1e-14 && creals (Result[i]) < 1e-18 ) + if ( creals(out[i]) < 3e-16 && creals (Result[i]) < 1e-18 ) assert ( 1 ) ; else - assert ( fabs( creals(out[i]) - creals (Result[i]) ) / fabs (creals (out[i])) < 1e-4 ); + assert ( fabs( creals(out[i]) - creals (Result[i]) ) / fabs (creals (out[i])) < 3e-6 ); - if ( cimags(out[i]) < 1e-14 && cimags (Result[i]) < 1e-18 ) + if ( cimags(out[i]) < 3e-16 && cimags (Result[i]) < 1e-18 ) assert ( 1 ) ; else - assert ( fabs( cimags(out[i]) - cimags (Result[i]) ) / fabs (cimags (out[i])) < 1e-4 ) ; + assert ( fabs( cimags(out[i]) - cimags (Result[i]) ) / fabs (cimags (out[i])) < 3e-6 ) ; } } @@ -527,16 +527,16 @@ static void zexpmaTest ( void) { fabs( zreals(out[i]) - zreals (Result[i]) ) / fabs (zreals (out[i])) , fabs( zimags(out[i]) - zimags (Result[i]) ) / fabs (zimags (out[i]))); - if ( zreals(out[i]) < 1e-14 && zreals (Result[i]) < 1e-18 ) + if ( zreals(out[i]) < 3e-16 && zreals (Result[i]) < 1e-18 ) assert ( 1 ) ; else - assert ( fabs( zreals(out[i]) - zreals (Result[i]) ) / fabs (zreals (out[i])) < 1e-12 ); + assert ( fabs( zreals(out[i]) - zreals (Result[i]) ) / fabs (zreals (out[i])) < 3e-16 ); - if ( zimags(out[i]) < 1e-14 && zimags (Result[i]) < 1e-18 ) + if ( zimags(out[i]) < 3e-16 && zimags (Result[i]) < 1e-18 ) assert ( 1 ) ; else - assert ( fabs( zimags(out[i]) - zimags (Result[i]) ) / fabs (zimags (out[i])) < 1e-12 ) ; + assert ( fabs( zimags(out[i]) - zimags (Result[i]) ) / fabs (zimags (out[i])) < 3e-16 ) ; } } |