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/powm/testDoublePowm.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/powm/testDoublePowm.c')
| -rw-r--r-- | src/matrixOperations/powm/testDoublePowm.c | 19 |
1 files changed, 8 insertions, 11 deletions
diff --git a/src/matrixOperations/powm/testDoublePowm.c b/src/matrixOperations/powm/testDoublePowm.c index 1eb5a5e6..a81baf68 100644 --- a/src/matrixOperations/powm/testDoublePowm.c +++ b/src/matrixOperations/powm/testDoublePowm.c @@ -40,14 +40,11 @@ static void dpowmaTest(void){ dpowma(in2, 4, expand2, out2); for (i=0;i<4;i++) { - assert( fabs(out1[i]-result1R[i]) / fabs(out1[i]) <1e-15); + assert( fabs(out1[i]-result1R[i]) / fabs(out1[i]) <3e-16); } -/* - FIXME : assert 1e-11 maybe due to spec2 -*/ for (i=0;i<16;i++) { - assert( fabs(out2[i]-result2R[i]) / fabs(out2[i]) <1e-14); + assert( fabs(out2[i]-result2R[i]) / fabs(out2[i]) <3e-16); } } @@ -74,8 +71,8 @@ static void zpowmaTest(void){ for (i=0;i<9;i++) printf("out[%d] = %f+%f*i\n",i,zreals(out[i]),zimags(out[i])); for (i=0;i<9;i++){ - assert( (fabs(zreals(out[i])-resultR[i])/ fabs(zreals(out[i])) ) <1e-14); - assert( (fabs(zimags(out[i])-resultI[i])/ fabs(zimags(out[i])) ) <1e-14); + assert( (fabs(zreals(out[i])-resultR[i])/ fabs(zreals(out[i])) ) <3e-16); + assert( (fabs(zimags(out[i])-resultI[i])/ fabs(zimags(out[i])) ) <3e-16); } } @@ -116,8 +113,8 @@ static void zpowmaTest(void){ zpowma(in2, 4, DoubleComplex(expand2,0), out2); for (i=0;i<4;i++) { - assert( fabs(zreals(out1[i])-result1R[i]) / fabs(zreals(out1[i])) <1e-15); - assert( fabs(zimags(out1[i])-result1I[i]) / fabs(zimags(out1[i])) <1e-15); + assert( fabs(zreals(out1[i])-result1R[i]) / fabs(zreals(out1[i])) <3e-16); + assert( fabs(zimags(out1[i])-result1I[i]) / fabs(zimags(out1[i])) <3e-16); } /* @@ -125,8 +122,8 @@ static void zpowmaTest(void){ */ for (i=0;i<16;i++) { printf("out[%d] = %1.16f+%1.16f*i --- result = %1.16f+%1.16f*i\n",i,zreals(out2[i]),zimags(out2[i]),result2R[i],result2I[i]); - assert( fabs(zreals(out2[i])-result2R[i]) / fabs(zreals(out2[i])) <1e-14); - assert( fabs(zimags(out2[i])-result2I[i]) / fabs(zimags(out2[i])) <1e-11); + assert( fabs(zreals(out2[i])-result2R[i]) / fabs(zreals(out2[i])) <3e-16); + assert( fabs(zimags(out2[i])-result2I[i]) / fabs(zimags(out2[i])) <3e-16); } } |