From: Sam Moore <matches@ucc.asn.au>
Date: Wed, 26 Mar 2014 17:42:04 +0000 (+0800)
Subject: Tester to compare custom Real to IEEE float
X-Git-Url: https://git.ucc.asn.au/?a=commitdiff_plain;h=95dd47f77239fd252e0292dca21baac29513bc7b;p=ipdf%2Fcode.git

Tester to compare custom Real to IEEE float

Also got the endianness right this time. I think.
---

diff --git a/src/main.cpp b/src/main.cpp
index 339ea8d..1b2073e 100644
--- a/src/main.cpp
+++ b/src/main.cpp
@@ -15,7 +15,8 @@ int main(int argc, char ** argv)
 	else
 	{
 		Debug("Add random object");
-		doc.Add(RECT_FILLED, Rect(Random()*0.5, Random()*0.5, Random()*0.5, Random()*0.5));
+		//doc.Add(RECT_FILLED, Rect(Random()*0.5, Random()*0.5, Random()*0.5, Random()*0.5));
+		doc.Add(RECT_FILLED, Rect(0.25,0.25, 0.5, 0.5));
 	}
 	MainLoop(doc);
 	return 0;
diff --git a/src/real.h b/src/real.h
index 1d3fb8b..7b30949 100644
--- a/src/real.h
+++ b/src/real.h
@@ -7,8 +7,8 @@ namespace IPDF
 {
 
 //#define REAL_FLOAT
-#define REAL_DOUBLE
-//#define REAL_HALF
+//#define REAL_DOUBLE
+#define REAL_HALF
 
 #ifdef REAL_SINGLE
 	typedef float Real;
@@ -26,7 +26,7 @@ namespace IPDF
 			// mask out extra bits in exponent
 			 //1000 1111 1000 0000 0000 0011 1111 1111
 			// Endianness matters
-			a &= 0xFF3008F8;//0x8F8003FF;
+			a &= 0xFFC001F1; //0x8F8003FF;
 
 		}	
 	
diff --git a/src/tests/realtofloat.cpp b/src/tests/realtofloat.cpp
index b09bd49..32ee20b 100644
--- a/src/tests/realtofloat.cpp
+++ b/src/tests/realtofloat.cpp
@@ -1,9 +1,26 @@
 #include "main.h"
-float test = 1e4;
+#include <cmath>
+float test_min = 0;
+float test_max = 1e-6;
+unsigned test_count = 100;
+float error_thresh = (test_max - test_min)/1e1;
 int main(int argc, char ** argv)
 {
-	Real r(test);
-	Debug("test float %.20f", test);
-	Debug("test real %.20f", Float(r));
+	Debug("TEST STARTING - Comparing Float(Real(test)) for %u trials between %.20f and %.20f", test_count, test_min, test_max);
+	float error;
+	for (unsigned i = 0; i < test_count; ++i)
+	{
+		float test = test_min + (test_max-test_min)*((float)(rand() % (int)1e6)/1e6);
+		Real real(test);
+		float thiserror = abs(test - real.value);
+		error += thiserror;
+		Debug("#%u: |test %.20f - real %.20f| = %.20f [mean %f]", i, test, real.value, thiserror, error/(i+1));
+	}
+
+	if (error/test_count > error_thresh)
+	{
+		Fatal("TEST FAILED - Average error %.20f exceeds threshold %.20f", error/test_count, error_thresh);
+	}
+	Debug("TEST SUCCEEDED - Average error %.20f is %f percent of range", error/test_count, 1e2*(error/test_count) / (test_max - test_min));
 	return 0;
 }