Quake III中的一段小代码

以下是最近热议的新开源的Quake III源码中的一段tricky code，我稍作修改后如下：

 

float CarmackInvSqrt(float x, unsigned long magicNumber) { float xhalf = 0.5f*x; int i = *(int*)&x; // get bits for floating VALUE i = magicNumber - (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; }

对于其中魔法数magic number，代码作者Carmack的版本是0x5f375a86，Lomont 的版本是0x5f3759df。

发现那篇中文分析文章做得不够严谨。虽然仔细分析这个东西不是像人们想象那么神秘深奥，但也不能草草了事。在做以下之前还是要先膜拜一下作者Carmack的睿智，倒不是仅仅因为这一段代码，而是因为能写出n多的代码。

 

首先对于最后一段，尽管注释和分析都说是牛顿法，或者牛顿迭代法（Newton Iteration</