X-Git-Url: https://git.ucc.asn.au/?p=ipdf%2Fcode.git;a=blobdiff_plain;f=src%2Freal.h;h=0068472c040b3b4a4b3e682af2ca9436f3957bfe;hp=57036713cca488cef07bcca4a833720250363295;hb=c9ce50e17952b0feeec68758444a50302a0aafba;hpb=a269b3e29535918a390f448829a3459e4853425b diff --git a/src/real.h b/src/real.h index 5703671..0068472 100644 --- a/src/real.h +++ b/src/real.h @@ -3,53 +3,30 @@ #include "common.h" -namespace IPDF -{ +#define REAL_SINGLE 0 +#define REAL_DOUBLE 1 +#define REAL_LONG_DOUBLE 2 + +#ifndef REAL + #error "REAL was not defined!" +#endif -#define REAL_SINGLE -//#define REAL_DOUBLE -//#define REAL_HALF +namespace IPDF +{ + extern const char * g_real_name[]; -#ifdef REAL_SINGLE +#if REAL == REAL_SINGLE typedef float Real; inline float Float(Real r) {return r;} -#elif defined REAL_DOUBLE +#elif REAL == REAL_DOUBLE typedef double Real; inline double Float(Real r) {return r;} -#elif defined REAL_HALF - struct Real - { - Real() = default; - 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 &= 0xFFC001F1; //0x8F8003FF; - - } - - 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;} - - inline std::ostream & operator<<(std::ostream & os, Real & r) {return os << r.value;} // yuk - -#endif //REAL_HALF +#elif REAL == REAL_LONG_DOUBLE + typedef long double Real; + inline long double Float(Real r) {return r;} +#else + #error "Type of Real unspecified." +#endif //REAL }