X-Git-Url: https://git.ucc.asn.au/?p=ipdf%2Fcode.git;a=blobdiff_plain;f=src%2Freal.h;h=57036713cca488cef07bcca4a833720250363295;hp=1d3fb8b0ecc2b717bcbab282de65e83633bbbb36;hb=b0c05b83db47aa91cb2da4244f1401aaf2c41ba0;hpb=32cd67ec48271c176eba218e082fcf5824aa8e6c

diff --git a/src/real.h b/src/real.h
index 1d3fb8b..5703671 100644
--- a/src/real.h
+++ b/src/real.h
@@ -6,8 +6,8 @@
 namespace IPDF
 {
 
-//#define REAL_FLOAT
-#define REAL_DOUBLE
+#define REAL_SINGLE
+//#define REAL_DOUBLE
 //#define REAL_HALF
 
 #ifdef REAL_SINGLE
@@ -20,24 +20,29 @@ namespace IPDF
 	struct Real
 	{
 		Real() = default;
-		Real(float r) : value(r) 
+		Real(double r) : value(r) 
 		{
 			int & a = *((int*)(&value)); // um...
 			// mask out extra bits in exponent
 			 //1000 1111 1000 0000 0000 0011 1111 1111
 			// Endianness matters
-			a &= 0xFF3008F8;//0x8F8003FF;
+			a &= 0xFFC001F1; //0x8F8003FF;
 
 		}	
 	
-		Real operator+(float f) {return Real(value+f);}
-		Real operator-(float f) {return Real(value+f);}
-		Real operator/(float f) {return Real(value/f);}
-		Real operator*(float f) {return Real(value*f);}
-		Real operator+(const Real & r) {return this->operator+(r.value);}
-		Real operator-(const Real & r) {return this->operator-(r.value);}
-		Real operator*(const Real & r) {return this->operator*(r.value);}
-		Real operator/(const Real & r) {return this->operator/(r.value);}
+		Real operator+(float f) const {return Real(value+f);}
+		Real operator-(float f) const {return Real(value+f);}
+		Real operator/(float f) const {return Real(value/f);}
+		Real operator*(float f) const {return Real(value*f);}
+		Real operator+(const Real & r) const {return Real(this->value + r.value);}
+		Real operator-(const Real & r) const {return Real(this->value - r.value);}
+		Real operator*(const Real & r) const {return Real(this->value * r.value);}
+		Real operator/(const Real & r) const {return Real(this->value / r.value);}
+		Real & operator+=(const Real & r) {this->value += r.value; return *this;}
+		Real & operator-=(const Real & r) {this->value -= r.value; return *this;}
+		Real & operator/=(const Real & r) {this->value /= r.value; return *this;}
+		Real & operator*=(const Real & r) {this->value *= r.value; return *this;}
+
 		float value;
 	};
 	inline float Float(Real r) {return r.value;}