平方根倒数的一个快速近似计算算法

已知一个浮点数$x$，现在需要求$\frac{1}{\sqrt{x}}$，这是一个非线性方程组的求解问题，可以采用不动点迭代法、牛顿迭代法之类的方法来解决。但是早在上个世纪，在游戏《Quake III》中，开发者就采用了一种非常快速求解这个问题的算法，这个算法非常有意思，而且现在也被一些GPU在硬件上实现，下面来看一下它的思路。

首先，设$x$按照IEEE浮点数表示规则对应的无符号整数为$I_x$，$y=\frac{1}{\sqrt{x}}$，y对应的无符号整数为$I_y$，在IEEE的浮点数表示规则中，一个单精度浮点数由32位二进制位表示，最高位代表符号位$S$，23-31这8位表示指数$E$，0-22这23位表示尾数$F$，如果$x=sign\ast 2^{e_x}\ast (1+f)$，那么对应的二进制表示中：$S=0(sign=1),1(sign=-1)$， $E=e_x+127$，$F=2^{23}\ast f$

因此如果对x取对数，就有$$lo{g}_2(x)=lo{g}_2(sign\ast 2^{e_x}$$