X-Git-Url: https://git.ucc.asn.au/?p=ipdf%2Fcode.git;a=blobdiff_plain;f=src%2Freal.h;h=13395fb75cac4457db0db17777c4c981ae8668d5;hp=5d4e1333ee14bf0204696c614c5a53119fc4c235;hb=f28146cf72fc68c8d3690814b1f89d47b8c4e7b9;hpb=6ce000e7212d9f5db6e5998c41df15bcad2022c8 diff --git a/src/real.h b/src/real.h index 5d4e133..13395fb 100644 --- a/src/real.h +++ b/src/real.h @@ -3,6 +3,7 @@ #include "common.h" #include +#include #define REAL_SINGLE 0 @@ -13,11 +14,17 @@ #define REAL_RATIONAL_ARBINT 5 #define REAL_MPFRCPP 6 #define REAL_IRRAM 7 +#define REAL_PARANOIDNUMBER 8 +#define REAL_GMPRAT 9 #ifndef REALTYPE #error "REALTYPE was not defined!" #endif +#define XSTR(x) STR(x) +#define STR(x) #x +//#pragma message "REALTYPE = " XSTR(REALTYPE) + #if REALTYPE == REAL_VFPU #include "vfpu.h" #endif @@ -34,12 +41,21 @@ #if REALTYPE == REAL_MPFRCPP #include + #endif //REALTYPE #if REALTYPE == REAL_IRRAM #include "../contrib/iRRAM/include/iRRAM/lib.h" #endif +#if REALTYPE == REAL_PARANOIDNUMBER + #include "paranoidnumber.h" +#endif + +#if REALTYPE == REAL_GMPRAT + #include "gmprat.h" +#endif + namespace IPDF { extern const char * g_real_name[]; @@ -47,12 +63,15 @@ namespace IPDF #if REALTYPE == REAL_SINGLE typedef float Real; inline Real RealFromStr(const char * str) {return strtof(str, NULL);} + inline std::string Str(const Real & a) {std::stringstream s; s << a; return s.str();} #elif REALTYPE == REAL_DOUBLE typedef double Real; inline Real RealFromStr(const char * str) {return strtod(str, NULL);} + inline std::string Str(const Real & a) {std::stringstream s; s << a; return s.str();} #elif REALTYPE == REAL_LONG_DOUBLE typedef long double Real; inline Real RealFromStr(const char * str) {return strtold(str, NULL);} + inline std::string Str(const Real & a) {std::stringstream s; s << a; return s.str();} #elif REALTYPE == REAL_VFPU typedef VFPU::VFloat Real; inline float Float(const Real & r) {return r.m_value;} @@ -64,13 +83,14 @@ namespace IPDF inline double Double(const Real & r) {return r.ToDouble();} inline Real RealFromStr(const char * str) {return Real(strtod(str, NULL));} #elif REALTYPE == REAL_RATIONAL_ARBINT - #define ARBINT Gmpint // Set to Gmpint or Arbint here + #define ARBINT Arbint // Set to Gmpint or Arbint here typedef Rational Real; inline float Float(const Real & r) {return (float)r.ToDouble();} inline double Double(const Real & r) {return r.ToDouble();} inline int64_t Int64(const Real & r) {return r.ToInt64();} inline Rational Sqrt(const Rational & r) {return r.Sqrt();} + inline Real RealFromStr(const char * str) {return Real(strtod(str, NULL));} #elif REALTYPE == REAL_MPFRCPP typedef mpfr::mpreal Real; inline double Double(const Real & r) {return r.toDouble();} @@ -79,6 +99,10 @@ namespace IPDF inline Real Sqrt(const Real & r) {return mpfr::sqrt(r, mpfr::mpreal::get_default_rnd());} inline Real Abs(const Real & r) {return mpfr::abs(r, mpfr::mpreal::get_default_rnd());} inline Real RealFromStr(const char * str) {return Real(strtod(str, NULL));} + inline std::string Str(const mpfr::mpreal & a) {std::stringstream s; s << a; return s.str();} + inline size_t Size(mpfr::mpreal & a) {return a.get_prec();} + inline mpfr::mpreal Log10(const mpfr::mpreal & a) {return mpfr::log10(a);} + #elif REALTYPE == REAL_IRRAM typedef iRRAM::REAL Real; inline double Double(const Real & r) {return r.as_double(53);} @@ -86,6 +110,24 @@ namespace IPDF inline int64_t Int64(const Real & r) {return (int64_t)r.as_double(53);} inline Real Sqrt(const Real & r) {return iRRAM::sqrt(r);} inline Real RealFromStr(const char * str) {return Real(strtod(str, NULL));} +#elif REALTYPE == REAL_PARANOIDNUMBER + typedef ParanoidNumber Real; + inline double Double(const Real & r) {return r.Digit();} + inline float Float(const Real & r) {return r.Digit();} + inline int64_t Int64(const Real & r) {return (int64_t)r.Digit();} + inline Real Sqrt(const Real & r) {return Real(sqrt(r.Digit()));} + inline Real RealFromStr(const char * str) {return Real(str);} + inline Real Abs(const Real & a) {return Real(fabs(a.Digit()));} +#elif REALTYPE == REAL_GMPRAT + typedef Gmprat Real; + inline double Double(const Real & r) {return r.ToDouble();} + inline float Float(const Real & r) {return (float)(r.ToDouble());} + inline int64_t Int64(const Real & r) {return (int64_t)r.ToDouble();} + inline Real Sqrt(const Real & r) {return Real(sqrt(r.ToDouble()));} + inline Real RealFromStr(const char * str) {return Real(strtod(str, NULL));} + inline Real Abs(const Real & a) {return (a > Real(0)) ? a : Real(0)-a;} + + #else #error "Type of Real unspecified." #endif //REALTYPE @@ -100,7 +142,26 @@ namespace IPDF inline double Double(long double f) {return (double)(f);} inline double Sqrt(double f) {return sqrt(f);} inline double Abs(double a) {return fabs(a);} - inline int64_t Int64(double a){return (int64_t)a;} + inline double Log10(double a) {return log(a)/log(10.0);} + inline size_t Size(float a) {return sizeof(a);} + inline size_t Size(double a) {return sizeof(a);} + + + // Don't cause an exception + inline float ClampFloat(double d) + { + float f = (fabs(d) < FLT_MAX) ? ((fabs(d) > FLT_MIN) ? (float)d : FLT_MIN) : FLT_MAX; + return copysign(f, d); + } + + inline int64_t Int64(double a) + { + if (a < INT64_MIN) + return INT64_MIN; + if (a > INT64_MAX) + return INT64_MAX; + return (int64_t)(a); + } inline Real Power(const Real & a, int n) { @@ -117,7 +178,7 @@ namespace IPDF { Real x; Real y; - Vec2() : x(0), y(0) {} + Vec2() : x(0.0), y(0.0) {} Vec2(Real _x, Real _y) : x(_x), y(_y) {} Vec2(const std::pair & p) : x(p.first), y(p.second) {} #if REALTYPE != REAL_IRRAM @@ -144,10 +205,106 @@ namespace IPDF }; + //TODO: Make sure there is actually a RealFromStr(const char * str) function + // Or this will recurse infinitely + // (If you remove this it will also break). inline Real RealFromStr(const std::string & str) {return RealFromStr(str.c_str());} + // things stolen from wikipedia and googling + inline const char * HumanScale(double f) + { + if (f < 1e-36) + return "RATHER SMALL"; + if (f < 1e-35) + return "Plank Length"; + if (f < 1e-25) + return "Turtles all the way"; + if (f < 1e-24) + return "More turtles"; + if (f < 1e-23) + return "Turtles"; + if (f < 1e-22) + return "This small"; + if (f < 1e-21) + return "To find things"; + if (f < 1e-20) + return "It is pretty difficult"; + if (f < 1e-19) + return "Not much"; + if (f < 1e-17) + return "Weak Force"; + if (f < 1e-16) + return "Proton"; + if (f < 1e-15) + return "(Classical) Electron"; + if (f < 1e-11) + return "Inter atomic (still)"; + if (f < 1e-10) + return "Inter atomic"; + if (f < 1e-9) + return "Atom"; + if (f < 1e-8) + return "DNA"; + if (f < 1e-7) + return "Virus"; + if (f < 1e-6) + return "Light"; + if (f < 1e-5) + return "Bacteria"; + if (f < 1e-4) + return "4004 Transistor"; + if (f < 1e-3) + return "Ant"; + if (f < 1e-2) + return "Coin"; + if (f < 1e-1) + return "iPhone"; + if (f < 1e0) + return "Person"; + if (f < 1e1) + return "Building"; + if (f < 1e2) + return "Football Field"; + if (f < 1e3) + return "Mountain"; + if (f < 1e4) + return "Clouds"; + if (f < 1e5) + return "Countries"; + if (f < 1e6) + return "Earth"; + if (f < 1e8) + return "Between Earth and Moon"; + if (f < 1e9) + return "Solar System"; + if (f < 1e13) + return "Distance to nearest Star"; + if (f < 1e21) + return "Milky Way"; + if (f < 1e26) + return "Universe"; + if (f < 1e27) + return "A bigger Universe"; + if (f < 1e28) + return "Really big things"; + if (f < 1e29) + return "Almost as big as..."; + if (f < 1e30) + return "Wolfram's Magestic Ego"; + return "QUITE BIG"; + + } + inline void DebugRealInfo() + { + Debug("Compiled with REAL = %d => \"%s\" sizeof(Real) == %d bytes", REALTYPE, g_real_name[REALTYPE], sizeof(Real)); + #if REALTYPE == REAL_PARANOIDNUMBER + #ifdef PARANOID_SIZE_LIMIT + Debug("Size limit of %d is being enforced", PARANOID_SIZE_LIMIT); + #endif + #endif + } } #endif //_REAL_H