当前位置：首页 > news >正文

Fast Inverse Square Root（快速平方根倒数算法）

news 2026/3/27 6:04:42

Fast Inverse Square Root（快速平方根倒数算法）被誉为计算机图形学史上最“臭名昭著”也最天才的代码。

这段代码最早出现在《雷神之锤 III 竞技场》（Quake III Arena）的源码中。

在那个 3D 算力贫瘠的年代，它用一种“近乎巫术”的位操作，将原本极慢的浮点数除法和开方运算，变成了一次简单的整数减法和位移。

1. 为什么要计算“平方根倒数”？

在 3D 引擎中，计算光的反射、阴影或物理碰撞，核心操作是向量归一化。
归一化的公式是：

这意味着每一帧画面，显卡/CPU 都要进行成千上万次的运算。在 90 年代，浮点数开方和除法指令极其昂贵（消耗几十个时钟周期），而这正是游戏掉帧的元凶。

2. 传奇代码解析（C 语言版）

这是《雷神之锤 III》源码中原封不动的注释，甚至保留了那个著名的脏话：

float Q_rsqrt( float number )
{long i;float x2, y;const float threehalfs = 1.5F;x2 = number * 0.5F;y  = number;i  = * ( long * ) &y;                       // 邪恶的位平面浮点转换i  = 0x5f3759df - ( i >> 1 );               // 这他妈是怎么回事？ (what the xxxx?)y  = * ( float * ) &i;y  = y * ( threehalfs - ( x2 * y * y ) );   // 第一轮牛顿迭代// y  = y * ( threehalfs - ( x2 * y * y ) );   // 第二轮迭代（可以省略）return y;
}

核心步骤拆解：

邪恶的类型转换 (Type Punning)：
i = * ( long * ) &y;
这行代码在内存层面强制把一个 float 当作 long（整数）来看待。这正是你之前看到的 union 或 FieldOffset(0) 的作用！
那个“魔数” (The Magic Number)：
0x5f3759df
这一行利用了 IEEE 754 浮点数表示法的数学特性。通过把浮点数当整数右移一位（相当于除以 2），再用魔数减去它，竟然能在空间内近似得到的结果！
牛顿迭代法：
最后一步利用牛顿迭代法进行精度补偿。只需一次迭代，误差就能控制在 1% 以内，速度却是标准算法的 4 倍。

3. C# 深度实现：复刻这一奇迹

在 C# 中，我们不能像 C 语言那样随意玩弄指针（除非在 unsafe 模块下），但我们可以利用你刚学的 StructLayout 来实现这个“邪恶转换”。

using System.Runtime.InteropServices;public static class FastMath
{// 定义一个共用体，模拟 C 语言的指针强转[StructLayout(LayoutKind.Explicit)]struct FloatIntUnion{[FieldOffset(0)] public float FloatValue;[FieldOffset(0)] public int IntValue;}public static float InverseSqrt(float number){FloatIntUnion u;float x2 = number * 0.5f;float y = number;u.IntValue = 0; // 初始化u.FloatValue = y;// 1. 类型转换：把 float 看作 intint i = u.IntValue;// 2. 核心魔数操作：位移并相减// 这一步在对数空间内完成了对平方根倒数的近似估计i = 0x5f3759df - (i >> 1);// 3. 转换回 floatu.IntValue = i;y = u.FloatValue;// 4. 牛顿迭代法 (Newton's method)，提高精度// y = y * (1.5 - (x2 * y * y))y = y * (1.5f - (x2 * y * y));return y;}
}