From 05fb6a62c83ae7cea37aa2795e61697b3560f135 Mon Sep 17 00:00:00 2001 From: torset Date: Fri, 27 Feb 2009 13:30:30 +0000 Subject: change precision for tests : put it to 3e-16 for double and 3e-6 for float --- src/elementaryFunctions/fix/testDoubleFix.c | 46 ++++++++++++++++------------- 1 file changed, 25 insertions(+), 21 deletions(-) (limited to 'src/elementaryFunctions/fix/testDoubleFix.c') diff --git a/src/elementaryFunctions/fix/testDoubleFix.c b/src/elementaryFunctions/fix/testDoubleFix.c index 95c2e0b3..2d7c6cb5 100644 --- a/src/elementaryFunctions/fix/testDoubleFix.c +++ b/src/elementaryFunctions/fix/testDoubleFix.c @@ -18,25 +18,25 @@ static void dfixsTest(void) { assert ((fabs( dfixs(-2) - (-2))/fabs(dfixs(-2)))<1e-16); assert ((fabs( dfixs(-1.9) - (-1))/fabs(dfixs(-1.9)))<1e-16); - assert ((fabs( dfixs(-1.8) - (-1))/fabs(dfixs(-1.9)))<1e-16); - assert ((fabs( dfixs(-1.7) - (-1))/fabs(dfixs(-1.9)))<1e-16); - assert ((fabs( dfixs(-1.6) - (-1))/fabs(dfixs(-1.9)))<1e-16); - assert ((fabs( dfixs(-1.5) - (-1))/fabs(dfixs(-1.9)))<1e-16); - assert ((fabs( dfixs(-1.4) - (-1))/fabs(dfixs(-1.9)))<1e-16); - assert ((fabs( dfixs(-1.3) - (-1))/fabs(dfixs(-1.9)))<1e-16); - assert ((fabs( dfixs(-1.2) - (-1))/fabs(dfixs(-1.9)))<1e-16); - assert ((fabs( dfixs(-1.1) - (-1))/fabs(dfixs(-1.9)))<1e-16); - assert ((fabs( dfixs(-1.0) - (-1))/fabs(dfixs(-1.9)))<1e-16); + assert ((fabs( dfixs(-1.8) - (-1))/fabs(dfixs(-1.8)))<1e-16); + assert ((fabs( dfixs(-1.7) - (-1))/fabs(dfixs(-1.7)))<1e-16); + assert ((fabs( dfixs(-1.6) - (-1))/fabs(dfixs(-1.6)))<1e-16); + assert ((fabs( dfixs(-1.5) - (-1))/fabs(dfixs(-1.5)))<1e-16); + assert ((fabs( dfixs(-1.4) - (-1))/fabs(dfixs(-1.4)))<1e-16); + assert ((fabs( dfixs(-1.3) - (-1))/fabs(dfixs(-1.3)))<1e-16); + assert ((fabs( dfixs(-1.2) - (-1))/fabs(dfixs(-1.2)))<1e-16); + assert ((fabs( dfixs(-1.1) - (-1))/fabs(dfixs(-1.1)))<1e-16); + assert ((fabs( dfixs(-1.0) - (-1))/fabs(dfixs(-1.0)))<1e-16); assert (fabs( dfixs(-.9))<1e-16); - assert ((fabs( dfixs(1.0) - (1))/fabs(dfixs(1.9)))<1e-16); - assert ((fabs( dfixs(1.1) - (1))/fabs(dfixs(1.9)))<1e-16); - assert ((fabs( dfixs(1.2) - (1))/fabs(dfixs(1.9)))<1e-16); - assert ((fabs( dfixs(1.3) - (1))/fabs(dfixs(1.9)))<1e-16); - assert ((fabs( dfixs(1.4) - (1))/fabs(dfixs(1.9)))<1e-16); - assert ((fabs( dfixs(1.5) - (1))/fabs(dfixs(1.9)))<1e-16); - assert ((fabs( dfixs(1.6) - (1))/fabs(dfixs(1.9)))<1e-16); - assert ((fabs( dfixs(1.7) - (1))/fabs(dfixs(1.9)))<1e-16); - assert ((fabs( dfixs(1.8) - (1))/fabs(dfixs(1.9)))<1e-16); + assert ((fabs( dfixs(1.0) - (1))/fabs(dfixs(1.0)))<1e-16); + assert ((fabs( dfixs(1.1) - (1))/fabs(dfixs(1.1)))<1e-16); + assert ((fabs( dfixs(1.2) - (1))/fabs(dfixs(1.2)))<1e-16); + assert ((fabs( dfixs(1.3) - (1))/fabs(dfixs(1.3)))<1e-16); + assert ((fabs( dfixs(1.4) - (1))/fabs(dfixs(1.4)))<1e-16); + assert ((fabs( dfixs(1.5) - (1))/fabs(dfixs(1.5)))<1e-16); + assert ((fabs( dfixs(1.6) - (1))/fabs(dfixs(1.6)))<1e-16); + assert ((fabs( dfixs(1.7) - (1))/fabs(dfixs(1.7)))<1e-16); + assert ((fabs( dfixs(1.8) - (1))/fabs(dfixs(1.8)))<1e-16); assert ((fabs( dfixs(1.9) - (1))/fabs(dfixs(1.9)))<1e-16); assert ((fabs( dfixs(2.0) - (2))/fabs(dfixs(2.0)))<1e-16); @@ -107,7 +107,9 @@ static void dfixaTest(void) { dfixa(in,10,out); - for (i=0;i<10;i++) assert(fabs(out[i]-res[i]) <1e-16); + for (i=0;i<10;i++) + if (out[i]!=0) assert (fabs(out[i]-res[i])/fabs(out[i]) <1e-16); + else assert(fabs(out[i]-res[i]) <1e-16); } @@ -123,8 +125,10 @@ static void zfixaTest(void) { zfixa(in,10,out); for (i=0;i<10;i++){ - assert(fabs(zreals(out[i])-resR[i]) <1e-16); - assert(fabs(zimags(out[i])-resI[i]) <1e-16); + if (zreals(out[i])!=0) assert(fabs(zreals(out[i])-resR[i])/fabs(zreals(out[i])) <1e-16); + else assert(fabs(zreals(out[i])-resR[i]) <1e-16); + if (zimags(out[i])!=0) assert(fabs(zimags(out[i])-resI[i])/fabs(zimags(out[i])) <1e-16); + else assert(fabs(zimags(out[i])-resI[i]) <1e-16); } } -- cgit