【Beautiful Algorithm】Fast Inverse Square Root（快速求平方根倒数）

来自Quick ||| Algorithm, Fast Inverse Square Root $\frac{1}{\sqrt{x}}$

代码

float Q_rsqrt(float number) {
	long i;
	float x, y;
	const float threehalfs = 1.5F; // <1>

	x = number * 0.5F; // <2>
	y = number;
	i = *(long *) &y; // <3>
	i = 0x5f3759df - (i >> 1); // <4>
	y = *(float *) &i; // <5>
	y = y * (threehalfs - (x * y * y)); // 1st iteration <6>
	y = y * (threehalfs - (x * y * y)); // 2st iteration, can be removed // <7>
	return y;
}

浮点数的表示 - IEEE 754

对于float, 其浮点数表示如下图所示。

其他方面

• normalised numbers
• denormalised numbers
• NaN
• infinities
• 0 & -0

分析¹

• <1>用在<6>和<7>处，放到<6>和<7>处分析。

• <3>将y的内容变为long的形式，如果使用i = (long)y,只会得到取整后的值，所以需要改变指针的意义，使得某个指向long的地址指向y。<3>和<5>是配对的，<3>把y中的内容复制出来，<5>把新的内容放到y中。

• <4> 对于浮点数0:E:M,其值为
  $$(1+\frac{M}{2^{23}})\times 2^{E-127}$$
  对值取$lo{g}_2$,得
  $$\begin{array}{rl}lo{g}_2((1+& \frac{M}{2^{23}})\times 2^{E-127})\\ & =lo{g}_2(1+\frac{M}{2^{23}})+E-127\\ & =\frac{M}{2^{23}}+\mu +E-127(duetolo{g}_2(1+x)\approx x+\mu ,0<x<1）\\ & =\frac{1}{2^{23}}(M+2^{23}\times E)+\mu -127\end{array}$$
  其中$(M+2^{23}\times E)$就是float的数字表示结果，即$i$得值。
  $$i=(M+2^{23}\times E)$$
  令$\Gamma =\frac{1}{\sqrt{x}}$，其中$x=(1+\frac{M}{2^{23}})\times 2^{E-127}$，则
  $$\begin{array}{rl}{\log }_2(\Gamma )& ={\log }_2\frac{1}{\sqrt{x}}\\ forleft:{\log }_2\frac{1}{\sqrt{x}}& ={\log }_2(x^{-\frac{1}{2}})\\ & =-\frac{1}{2}{\log }_2(x)\\ & =-\frac{1}{2}(\frac{1}{2^{23}}(M_x+2^{23}\times E_x)+\mu -127)\\ forright:{\log }_2(\Gamma )& =\frac{1}{2^{23}}(M_{\Gamma }+2^{23}\times E_{\Gamma })+\mu -127\\ combineleftandright:& \\ (M_{\Gamma }+2^{23}\times E_{\Gamma })& =\frac{2}{3}{2}^{23}(127-\mu )-\frac{1}{2}(M_x+2^{23}\times E_x)\\ \Rightarrow & \\ (M_{\Gamma }+2^{23}\times E_{\Gamma })& =0\times 5f3759df-(i>>1)\end{array}$$
  $(M_{\Gamma }+2^{23}\times E_{\Gamma })$是$\Gamma =\frac{1}{\sqrt{x}}$的整数值表示，步骤<5>则将其转化为浮点数。

• <6>所示的是迭代公式，不断逼近目标值
  Newton Method 用于解决$f(x)=0$的零点查找问题。迭代公式为：
  $$x_{new}=x-\frac{f(x)}{f^{{}^{\prime }}(x)}$$
  对于本题，$x$为输入的number，我们需要求的值是$\frac{1}{\sqrt{(}x)}=y$，故而目标函数是：
  $$f(y)=\frac{1}{y^2}-x$$
  故而迭代公式是
  $$y_{new}=y(\frac{3}{2}-\frac{1}{2}x{y}^2)$$

评估本函数的误差

TODO

---

1. https://www.youtube.com/watch?v=p8u_k2LIZyo ↩︎