include/crbn/basic/fastmath.h

#ifndef __fast_math_h__
#define __fast_math_h__

#define __FASTMATH_EPSILON__ 2

#define fm_eps_eq( x, y ) (abs( (*(unsigned*)(&x)) - (*(unsigned*)(&y)) ) < __FASTMATH_EPSILON__)

inline int fm_neg( float f )
{
/*   return ((*(unsigned*)(&f)) & 0x80000000); */
  union {
    float f;
    unsigned int u;
  } conv;
  conv.f = f;
  return conv.u & 0x80000000;
}

inline float fm_abs( float f )
{
/*   unsigned u = (*(unsigned*)&f) & 0x7FFFFFFF; */
/*   return (*(float*)&u); */
  union {
    float f;
    unsigned int u;
  } conv;
 
  conv.f = f;
  conv.u &= 0x7FFFFFFF;
  return conv.f;  
}

#undef __FASTMATH_EPSILON__

#endif // __fast_math_h__

---