计算机与数学 —— 雷神之锤3源码中的快速逆平方根算法

这篇博客介绍了在雷神之锤3源代码中快速求逆平方根的算法。

---

源码

雷神之锤3中的逆平方根算法如下：

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 ;
}

函数接受一个浮点数作为输入，输出的结果是平方根的倒数。

单精度浮点数

在计算机中，单精度浮点数使用32位来储存表示：

在这32位中，最高位为符号位，后面的8位为整数ex，代表浮点数的指数，而最后的23位表示的是小数部分mx，小数部分第一位表示$2^{-1}$，第二位表示$2^{-2}$……

所以如果x是一个正浮点数，则有：

$$x=2^{e_x}(1+m_x)$$

另外，如果需要把一个浮点数转化为整数形式，则需要做如下的运算：

$$Ix=E_{X}L+M_x=L(e_x+B+m_x)$$

在上面的公式中，L指的是指数部分需要的唯一次数，这里是$2^{23}$，B是127，而M是小数部分对应的整数版本。

牛顿法

牛顿迭代法（简称牛顿法），是用于计算机求解任意连续函数的根值的一种方法。

假设对于如下的函数，我们想要求这个函数的根：

该如何求得这个根呢？首先，我们先猜一个解 x 0 x_{0}