fork download
copy 
  1. #include <stdio.h>
  2. #include <math.h>
  3.  
  4. // https://e...content-available-to-author-only...a.org/wiki/Fast_inverse_square_root
  5. float Q_rsqrt( float number )
  6. {
  7.   long i;
  8.   float x2, y;
  9.   const float threehalfs = 1.5F;
  10.  
  11.   x2 = number * 0.5F;
  12.   y  = number;
  13.   i  = * ( long * ) &y;                       // evil floating point bit level hacking
  14.   i  = 0x5f3759df - ( i >> 1 );               // what the fuck?
  15.   y  = * ( float * ) &i;
  16.   y  = y * ( threehalfs - ( x2 * y * y ) );   // 1st iteration
  17. //y  = y * ( threehalfs - ( x2 * y * y ) );   // 2nd iteration, this can be removed
  18.  
  19.   return y;
  20. }
  21.  
  22. void print(float n) {
  23.   float result_1 = 1 / n;
  24.   float result_2 = Q_rsqrt(n) * Q_rsqrt(n);
  25.   printf("%f %f\n", result_1, result_2);
  26. }
  27.  
  28. int main() {
  29.   print(1000.0f);
  30.   print(100.0f);
  31.   print(10.0f);
  32.   print(1.0f);
  33.   print(0.1f);
  34.   print(0.01f);
  35.   return 0;
  36. }
  37.  
Success #stdin #stdout 0s 4268KB
comments (?)
stdin
copy 
Standard input is empty
stdout
copy 
0.001000 0.000997
0.010000 0.009969
0.100000 0.099658
1.000000 0.996617
10.000000 9.968114
100.000000 99.650749
https://ideone.com/n8krf8
language:
C (clang 8.0)
created:
5 years ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!