+ inline Real RealFromStr(const char * str) {return strtold(str, NULL);}
+#elif REALTYPE == REAL_VFPU
+ typedef VFPU::VFloat Real;
+ inline float Float(const Real & r) {return r.m_value;}
+ inline double Double(const Real & r) {return r.m_value;}
+ inline Real RealFromStr(const char * str) {return Real(strtod(str, NULL));}
+#elif REALTYPE == REAL_RATIONAL
+ typedef Rational<int64_t> Real;
+ inline float Float(const Real & r) {return (float)r.ToDouble();}
+ 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 Arbint // Set to Gmpint or Arbint here
+
+ typedef Rational<ARBINT> 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<ARBINT> Sqrt(const Rational<ARBINT> & 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();}
+ inline float Float(const Real & r) {return r.toDouble();}
+ inline int64_t Int64(const Real & r) {return r.toLong();}
+ 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));}
+#elif REALTYPE == REAL_IRRAM
+ typedef iRRAM::REAL Real;
+ inline double Double(const Real & r) {return r.as_double(53);}
+ inline float Float(const Real & r) {return r.as_double(53);}
+ 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;}
+