源自 Quake3 的快速求 InvSqrt() 函数

ID Software 的 Quake3 代码里面有这么一段，可以快速求解 1/sqrt(x)：


float InvSqrt(float x){
float xhalf=0.5f*x;
long i=*(long*)&x;
i=0x5f3759df - (i>>1);
x=*(float *)&i;
x=x*(1.5f-xhalf*x*x);
return x;
} float InvSqrt(float x){
float xhalf=0.5f*x;
long i=*(long*)&x;
i=0x5f3759df - (i>>1);
x=*(float *)&i;
x=x*(1.5f-xhalf*x*x);
return x;
}


Beyond3D.com 的 Ryszard Sommefeldt 一直在想到底是哪个家伙写了这些神奇的代码？2003年Chris Lomont还写了一篇文章对这些代码进行了分析。

对于编译，OS 这些我还有很多要看，CSDN 上有篇帖子对上面的问题进行了讨论。