/Fast_Inv_Sqrt/C/jcla1/fast_inv_sqrt.c

https://github.com/andrewcyu/Algorithm-Implementations · C · 16 lines · 13 code · 2 blank · 1 comment · 0 complexity · 0068b7d8e4f9f1ba2ea3b4a43d63aaa7 MD5 · raw file

  1. float Q_rsqrt( float number )
  2. {
  3. long i;
  4. float x2, y;
  5. const float threehalfs = 1.5F;
  6. x2 = number * 0.5F;
  7. y = number;
  8. i = * ( long * ) &y; // evil floating point bit level hacking
  9. i = 0x5f3759df - ( i >> 1 ); // what the fuck?
  10. y = * ( float * ) &i;
  11. y = y * ( threehalfs - ( x2 * y * y ) ); // 1st iteration
  12. // y = y * ( threehalfs - ( x2 * y * y ) ); // 2nd iteration, this can be removed
  13. return y;
  14. }