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

2 files changed, 52 insertions, 46 deletions

diff --git a/src/elementaryFunctions/floor/testDoubleFloor.c b/src/elementaryFunctions/floor/testDoubleFloor.c
index 186aba95..71e0b36f 100644
--- a/src/elementaryFunctions/floor/testDoubleFloor.c
+++ b/src/elementaryFunctions/floor/testDoubleFloor.c
@@ -20,25 +20,25 @@
 static void dfloorsTest(void) {
 	assert ((fabs( dfloors(-2) - (-2))/fabs(dfloors(-2)))<1e-16);	
 	assert ((fabs( dfloors(-1.9) - (-2))/fabs(dfloors(-1.9)))<1e-16);		
-	assert ((fabs( dfloors(-1.8) - (-2))/fabs(dfloors(-1.9)))<1e-16);	
-	assert ((fabs( dfloors(-1.7) - (-2))/fabs(dfloors(-1.9)))<1e-16);	
-	assert ((fabs( dfloors(-1.6) - (-2))/fabs(dfloors(-1.9)))<1e-16);	
-	assert ((fabs( dfloors(-1.5) - (-2))/fabs(dfloors(-1.9)))<1e-16);	
-	assert ((fabs( dfloors(-1.4) - (-2))/fabs(dfloors(-1.9)))<1e-16);	
-	assert ((fabs( dfloors(-1.3) - (-2))/fabs(dfloors(-1.9)))<1e-16);	
-	assert ((fabs( dfloors(-1.2) - (-2))/fabs(dfloors(-1.9)))<1e-16);	
-	assert ((fabs( dfloors(-1.1) - (-2))/fabs(dfloors(-1.9)))<1e-16);	
-	assert ((fabs( dfloors(-1.0) - (-1))/fabs(dfloors(-1.9)))<1e-16);		
-	assert ((fabs( dfloors(-0.9) - (-1))/fabs(dfloors(-1.9)))<1e-16);
-	assert ((fabs( dfloors(1.0) - (1))/fabs(dfloors(1.9)))<1e-16);	
-	assert ((fabs( dfloors(1.1) - (1))/fabs(dfloors(1.9)))<1e-16);	
-	assert ((fabs( dfloors(1.2) - (1))/fabs(dfloors(1.9)))<1e-16);	
-	assert ((fabs( dfloors(1.3) - (1))/fabs(dfloors(1.9)))<1e-16);	
-	assert ((fabs( dfloors(1.4) - (1))/fabs(dfloors(1.9)))<1e-16);	
-	assert ((fabs( dfloors(1.5) - (1))/fabs(dfloors(1.9)))<1e-16);	
-	assert ((fabs( dfloors(1.6) - (1))/fabs(dfloors(1.9)))<1e-16);	
-	assert ((fabs( dfloors(1.7) - (1))/fabs(dfloors(1.9)))<1e-16);	
-	assert ((fabs( dfloors(1.8) - (1))/fabs(dfloors(1.9)))<1e-16);	
+	assert ((fabs( dfloors(-1.8) - (-2))/fabs(dfloors(-1.8)))<1e-16);	
+	assert ((fabs( dfloors(-1.7) - (-2))/fabs(dfloors(-1.7)))<1e-16);	
+	assert ((fabs( dfloors(-1.6) - (-2))/fabs(dfloors(-1.6)))<1e-16);	
+	assert ((fabs( dfloors(-1.5) - (-2))/fabs(dfloors(-1.5)))<1e-16);	
+	assert ((fabs( dfloors(-1.4) - (-2))/fabs(dfloors(-1.4)))<1e-16);	
+	assert ((fabs( dfloors(-1.3) - (-2))/fabs(dfloors(-1.3)))<1e-16);	
+	assert ((fabs( dfloors(-1.2) - (-2))/fabs(dfloors(-1.2)))<1e-16);	
+	assert ((fabs( dfloors(-1.1) - (-2))/fabs(dfloors(-1.1)))<1e-16);	
+	assert ((fabs( dfloors(-1.0) - (-1))/fabs(dfloors(-1.0)))<1e-16);		
+	assert ((fabs( dfloors(-0.9) - (-1))/fabs(dfloors(-0.9)))<1e-16);
+	assert ((fabs( dfloors(1.0) - (1))/fabs(dfloors(1.0)))<1e-16);	
+	assert ((fabs( dfloors(1.1) - (1))/fabs(dfloors(1.1)))<1e-16);	
+	assert ((fabs( dfloors(1.2) - (1))/fabs(dfloors(1.2)))<1e-16);	
+	assert ((fabs( dfloors(1.3) - (1))/fabs(dfloors(1.3)))<1e-16);	
+	assert ((fabs( dfloors(1.4) - (1))/fabs(dfloors(1.4)))<1e-16);	
+	assert ((fabs( dfloors(1.5) - (1))/fabs(dfloors(1.5)))<1e-16);	
+	assert ((fabs( dfloors(1.6) - (1))/fabs(dfloors(1.6)))<1e-16);	
+	assert ((fabs( dfloors(1.7) - (1))/fabs(dfloors(1.7)))<1e-16);	
+	assert ((fabs( dfloors(1.8) - (1))/fabs(dfloors(1.8)))<1e-16);	
 	assert ((fabs( dfloors(1.9) - (1))/fabs(dfloors(1.9)))<1e-16);	
 	assert ((fabs( dfloors(2.0) - (2))/fabs(dfloors(2.0)))<1e-16);
 
@@ -110,9 +110,9 @@ static void dflooraTest(void) {
 	
 	dfloora(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);
 }
 
 static void zflooraTest(void) {
@@ -127,8 +127,10 @@ static void zflooraTest(void) {
 	zfloora(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);
 	 }
 
 }
diff --git a/src/elementaryFunctions/floor/testFloatFloor.c b/src/elementaryFunctions/floor/testFloatFloor.c
index 1b171893..df49ee0d 100644
--- a/src/elementaryFunctions/floor/testFloatFloor.c
+++ b/src/elementaryFunctions/floor/testFloatFloor.c
@@ -18,25 +18,25 @@
 static void sfloorsTest(void) {
 	assert ((fabs( sfloors(-2) - (-2))/fabs(sfloors(-2)))<1e-16);	
 	assert ((fabs( sfloors(-1.9f) - (-2))/fabs(sfloors(-1.9f)))<1e-16);		
-	assert ((fabs( sfloors(-1.8f) - (-2))/fabs(sfloors(-1.9f)))<1e-16);	
-	assert ((fabs( sfloors(-1.7f) - (-2))/fabs(sfloors(-1.9f)))<1e-16);	
-	assert ((fabs( sfloors(-1.6f) - (-2))/fabs(sfloors(-1.9f)))<1e-16);	
-	assert ((fabs( sfloors(-1.5f) - (-2))/fabs(sfloors(-1.9f)))<1e-16);	
-	assert ((fabs( sfloors(-1.4f) - (-2))/fabs(sfloors(-1.9f)))<1e-16);	
-	assert ((fabs( sfloors(-1.3f) - (-2))/fabs(sfloors(-1.9f)))<1e-16);	
-	assert ((fabs( sfloors(-1.2f) - (-2))/fabs(sfloors(-1.9f)))<1e-16);	
-	assert ((fabs( sfloors(-1.1f) - (-2))/fabs(sfloors(-1.9f)))<1e-16);	
-	assert ((fabs( sfloors(-1.0f) - (-1))/fabs(sfloors(-1.9f)))<1e-16);	
-	assert ((fabs( sfloors(-0.9f) - (-1))/fabs(sfloors(-1.9f)))<1e-16);	
-	assert ((fabs( sfloors(1.0f) - (1))/fabs(sfloors(1.9f)))<1e-16);	
-	assert ((fabs( sfloors(1.1f) - (1))/fabs(sfloors(1.9f)))<1e-16);	
-	assert ((fabs( sfloors(1.2f) - (1))/fabs(sfloors(1.9f)))<1e-16);	
-	assert ((fabs( sfloors(1.3f) - (1))/fabs(sfloors(1.9f)))<1e-16);	
-	assert ((fabs( sfloors(1.4f) - (1))/fabs(sfloors(1.9f)))<1e-16);	
-	assert ((fabs( sfloors(1.5f) - (1))/fabs(sfloors(1.9f)))<1e-16);	
-	assert ((fabs( sfloors(1.6f) - (1))/fabs(sfloors(1.9f)))<1e-16);	
-	assert ((fabs( sfloors(1.7f) - (1))/fabs(sfloors(1.9f)))<1e-16);	
-	assert ((fabs( sfloors(1.8f) - (1))/fabs(sfloors(1.9f)))<1e-16);	
+	assert ((fabs( sfloors(-1.8f) - (-2))/fabs(sfloors(-1.8f)))<1e-16);	
+	assert ((fabs( sfloors(-1.7f) - (-2))/fabs(sfloors(-1.7f)))<1e-16);	
+	assert ((fabs( sfloors(-1.6f) - (-2))/fabs(sfloors(-1.6f)))<1e-16);	
+	assert ((fabs( sfloors(-1.5f) - (-2))/fabs(sfloors(-1.5f)))<1e-16);	
+	assert ((fabs( sfloors(-1.4f) - (-2))/fabs(sfloors(-1.4f)))<1e-16);	
+	assert ((fabs( sfloors(-1.3f) - (-2))/fabs(sfloors(-1.3f)))<1e-16);	
+	assert ((fabs( sfloors(-1.2f) - (-2))/fabs(sfloors(-1.2f)))<1e-16);	
+	assert ((fabs( sfloors(-1.1f) - (-2))/fabs(sfloors(-1.1f)))<1e-16);	
+	assert ((fabs( sfloors(-1.0f) - (-1))/fabs(sfloors(-1.0f)))<1e-16);	
+	assert ((fabs( sfloors(-0.9f) - (-1))/fabs(sfloors(-0.9f)))<1e-16);	
+	assert ((fabs( sfloors(1.0f) - (1))/fabs(sfloors(1.0f)))<1e-16);	
+	assert ((fabs( sfloors(1.1f) - (1))/fabs(sfloors(1.1f)))<1e-16);	
+	assert ((fabs( sfloors(1.2f) - (1))/fabs(sfloors(1.2f)))<1e-16);	
+	assert ((fabs( sfloors(1.3f) - (1))/fabs(sfloors(1.3f)))<1e-16);	
+	assert ((fabs( sfloors(1.4f) - (1))/fabs(sfloors(1.4f)))<1e-16);	
+	assert ((fabs( sfloors(1.5f) - (1))/fabs(sfloors(1.5f)))<1e-16);	
+	assert ((fabs( sfloors(1.6f) - (1))/fabs(sfloors(1.6f)))<1e-16);	
+	assert ((fabs( sfloors(1.7f) - (1))/fabs(sfloors(1.7f)))<1e-16);	
+	assert ((fabs( sfloors(1.8f) - (1))/fabs(sfloors(1.8f)))<1e-16);	
 	assert ((fabs( sfloors(1.9f) - (1))/fabs(sfloors(1.9f)))<1e-16);	
 	assert ((fabs( sfloors(2.0f) - (2))/fabs(sfloors(2.0f)))<1e-16);
 }
@@ -107,7 +107,9 @@ static void sflooraTest(void) {
 	
 	sfloora(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);
 }
 
 static void cflooraTest(void) {
@@ -122,8 +124,10 @@ static void cflooraTest(void) {
 	cfloora(in,10,out);
 	
 	for (i=0;i<10;i++){
-	 	assert(fabs(creals(out[i])-resR[i]) <1e-16);
-	 	assert(fabs(cimags(out[i])-resI[i]) <1e-16);
+		if (creals(out[i])!=0) assert(fabs(creals(out[i])-resR[i])/fabs(creals(out[i])) <1e-16);
+		else 	assert(fabs(creals(out[i])-resR[i]) <1e-16);
+		if (cimags(out[i])!=0) assert(fabs(cimags(out[i])-resI[i])/fabs(cimags(out[i])) <1e-16);
+	 	else 	assert(fabs(cimags(out[i])-resI[i]) <1e-16);
 	 }
 }