+ // Allow us to call Float on the primative types
+ // Useful so I can template some things that could be either (a more complicated) Real or a primitive type
+ // Mostly in the testers.
+ inline float Float(float f) {return (float)f;}
+ inline float Float(double f) {return (float)f;}
+ inline float Float(long double f) {return (float)(f);}
+ inline double Double(float f) {return (double)f;}
+ inline double Double(double f) {return (double)f;}
+ inline double Double(long double f) {return (double)(f);}
+ inline double Sqrt(double f) {return sqrt(f);}
+
+ inline Real Power(const Real & a, int n)
+ {
+ if (n < 0)
+ {
+ return Power(Real(1)/a, -n);
+ }
+ Real r(1);
+ for (int i = 0; i < n; ++i)
+ r *= a;
+ return r;
+ }
+
+
+