#define FLT_MASK   0xBE6F0000
#define DBL_MASK   0xBFCD4600

union FTOS
{
   FTOS( const float X )
      : f( X )
   {
i = ( FLT_MASK - i ) >> 1;   
   }
   float  f;
   size_t i;
};

inline float fast_sqrt( const float X )
{
   float f1, f2;
   FTOS Tmp = X;

   f1 = Tmp.f;
   f2 = X * 0.5f;

   f1 *= 1.5f - f1 * f1 * f2; /* Dup this line to increase precision. */
   f1 *= 1.5f - f1 * f1 * f2;
   
   return f1 * X;
}