/test/fast_inv_sqrt.c

https://github.com/cptjazz/rheo · C · 17 lines · 13 code · 2 blank · 2 comment · 0 complexity · 6941306fecab9e409c60cc821ba04606 MD5 · raw file

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