change precision for tests : put it to 3e-16 for double and 3e-6 for float

1 files changed, 25 insertions, 21 deletions

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);
 	 }
 
 }