X-Git-Url: https://git.ucc.asn.au/?p=ipdf%2Fcode.git;a=blobdiff_plain;f=src%2Freal.h;h=20ec68101b70611a6b76a7172b4c9ffc932bf301;hp=c5b16ac82100853677dd84b6b6a2502ce44e457a;hb=6c0dfe752994312ee58d307b383948bfeb2d6e2e;hpb=ea748154f1bc7dbc81cb52611a52865e63109439

diff --git a/src/real.h b/src/real.h
index c5b16ac..20ec681 100644
--- a/src/real.h
+++ b/src/real.h
@@ -3,6 +3,7 @@
 
 #include "common.h"
 #include <cmath>
+#include <cfloat>
 
 
 #define REAL_SINGLE 0
@@ -20,6 +21,10 @@
 	#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
@@ -36,6 +41,7 @@
 
 #if REALTYPE == REAL_MPFRCPP
 	#include <mpreal.h>
+
 #endif //REALTYPE
 
 #if REALTYPE == REAL_IRRAM
@@ -57,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;}
@@ -90,6 +99,11 @@ 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);}	
+	inline mpfr::mpreal Exp(const mpfr::mpreal & a) {return mpfr::pow(2.817, a);}
+	
 #elif REALTYPE == REAL_IRRAM
 	typedef iRRAM::REAL Real;
 	inline double Double(const Real & r) {return r.as_double(53);}
@@ -113,6 +127,7 @@ namespace IPDF
 	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;}
+	inline std::string Str(const Real & r) {return r.Str();}
 	
 #else
 	#error "Type of Real unspecified."
@@ -128,7 +143,27 @@ 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(const Real & a)
+	{
+		double d = Double(a);
+		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)
 	{
@@ -178,7 +213,108 @@ namespace IPDF
 	inline Real RealFromStr(const std::string & str) {return RealFromStr(str.c_str());}
 
 
+	// things stolen from wikipedia and googling
+	inline const char * HumanScale(const Real & r)
+	{
+		double f = Double(r);
+		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 < 0.5)
+			return "iPhone";
+		if (f < 5)
+			return "Person";
+		if (f < 1e2)
+			return "Building";
+		if (f < 1e3)
+			return "Football Field";
+		if (f < 2e3)
+			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
+		#if REALTYPE == REAL_MPFRCPP
+			Debug("Precision of MPFR is %d", mpfr_get_default_prec());
+		#endif
+		
+		#ifdef TRANSFORM_BEZIERS_TO_PATH
+		Debug("PathReal = %d => \"%s\"", PATHREAL, g_real_name[PATHREAL]);
+		#endif
+	}
 }
 
 #endif //_REAL_H