Add dist and distSqr functions

[simgear.git] / simgear / math / fastmath.hxx

diff --git a/simgear/math/fastmath.hxx b/simgear/math/fastmath.hxx
index b2fd735b31d583875e3ba04f56c2cc666cafba4c..97851e33e9e88bd2801c2ed3bfc645296197789d 100644 (file)
--- a/simgear/math/fastmath.hxx
+++ b/simgear/math/fastmath.hxx
@@ -2,7 +2,7 @@
  * \file fastmath.hxx
  * fast mathematics routines.
  *
- * Refferences:
+ * References:
  *
  * A Fast, Compact Approximation of the Exponential Function
  * Nicol N. Schraudolph
@@ -51,16 +51,18 @@ void fast_BSR(float &x, register unsigned long shiftAmount);
 
 inline float fast_log2 (float val)
 {
-    int * const  exp_ptr = reinterpret_cast <int *> (&val);
-    int          x = *exp_ptr;
-    const int    log_2 = ((x >> 23) & 255) - 128;
-    x &= ~(255 << 23);
-    x += 127 << 23;
-    *exp_ptr = x;
-
-    val = ((-1.0f/3) * val + 2) * val - 2.0f/3;   // (1)
-
-    return (val + log_2);
+    union {
+        float f;
+        int i;
+    } v;
+    v.f = val;
+    const int log_2 = ((v.i >> 23) & 255) - 128;
+    v.i &= ~(255 << 23);
+    v.i += 127 << 23;
+
+    v.f = ((-1.0f/3) * v.f + 2) * v.f - 2.0f/3;   // (1)
+
+    return (v.f + log_2);
 }