From 6e7f75a3a916ae7b7b6145c0badad7b649265c0f Mon Sep 17 00:00:00 2001
From: torset
Date: Fri, 27 Feb 2009 13:59:12 +0000
Subject: change precision for tests : put it to 3e-16 for double and 3e-6 for
 float

---
 src/matrixOperations/powm/testDoublePowm.c | 19 ++++++++-----------
 1 file changed, 8 insertions(+), 11 deletions(-)

(limited to 'src/matrixOperations/powm/testDoublePowm.c')

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