平方根を求めて遊ぶ。(ほぼニュートン法) (0x5f3759df)

参考

注意

• 過去ログをみよ!!!
• M5NanoC6での予定は、未定
• 非営利、研究調査で引用しました。

プログラム

オンラインコンパイラ paiza

//https://en.wikipedia.org/wiki/Fast_inverse_square_root
float q_rsqrt(float number)
{
  long i;
  float x2, y;
  const float threehalfs = 1.5F;

  x2 = number * 0.5F;
  y  = number;
  i  = * ( long * ) &y;                       // evil floating point bit level hacking
  i  = 0x5f3759df - ( i >> 1 );               // what the fuck?
  y  = * ( float * ) &i;
  y  = y * ( threehalfs - ( x2 * y * y ) );   // 1st iteration
  // y  = y * ( threehalfs - ( x2 * y * y ) );   // 2nd iteration, this can be removed

  return y;
}


#include <iostream>
using namespace std;
int main(void){
    // Your code here!
    
    float in;
    float f;
    
    in = 2.0;
    
    f = q_rsqrt(in); //逆平方根
    
    f = in * f;
    
    printf("%f\n",f);
    
}

1.413860