【编程】超快速计算平方根的倒数

C/C++代码技巧，那肯定是鼎鼎大名的 InvSqrt了，用来 快速计算平方根的倒数。

float InvSqrt(float x)
{
    float xhalf = 0.5f*x;
    int i = *(int*)&x;       // get bits for floating value
    i = 0x5f3759df - (i>>1); // gives initial guess y0
    x = *(float*)&i;         // convert bits back to float
    x = x*(1.5f-xhalf*x*x);  // Newton step, repeating increases accuracy
    return x;
}

- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

（以下摘自 0x5f375a86_百度百科）
这是一个传奇算法，此算法最早被认为是由 约翰·卡马克所发明，
但后来的调查显示，该算法在这之前就于计算机图形学的硬件与软件领域有所应用。
它的速度要比标准的牛顿迭代法快上 4 倍，
而其中的关键是一行神秘的代码和一个莫名其妙的数字 0x5f3759df。

普度大学的数学家Lomont觉得很好玩，
决定要研究一下卡马克弄出来的这个猜测值有什么奥秘。
Lomont也是个牛人，在精心研究之后从理论上也推导出一个最佳猜测值，
和卡马克的数字非常接近, 0x5f37642f。
卡马克真牛，他是外星人吗？

传奇并没有在这里结束。
Lomont计算出结果以后非常满意，于是拿自己计算出的起始值和卡马克的神秘数字做比赛，
看看谁的数字能够更快更精确的求得平方根。
结果是卡马克赢了...

谁也不知道卡马克是怎么找到这个数字的。
最后Lomont怒了，采用暴力方法一个数字一个数字试过来，
终于找到一个比卡马克数字要好上那么一丁点的数字，
虽然实际上这两个数字所产生的结果非常近似，
这个暴力得出的数字是 0x5f375a86。

Lomont为此写下一篇论文， Fast inverse square root （见wikipedia）

本文链接：http://www.zhihu.com/question/37692782/answer/73323805