マジックナンバー0x5f3759df の導き方について

0x5f3759df

は、ビット操作と対数の性質を組み合わせた非常に巧妙な近似手法によって導き出されます。この数値は、計算の誤差を最小化するために、数式による解析と経験的な調整によって最適化されました。 

以下に、導出の主要なステップをまとめます。 

1. 浮動小数点数の表現を利用する 

このアルゴリズムの核となるのは、IEEE 754形式の32ビット浮動小数点数（float）を、そのビットパターンをそのまま32ビット整数（long）として解釈するというアイデアです。 

IEEE 754形式では、32ビット浮動小数点数

$x$

は、符号ビット

$S$

、指数部

$E$

、仮数部

$M$

に分けられ、その値は以下の式で表されます。 

$x=(-1{)}^S\cdot (1+M)\cdot 2^{E-B}$

ここで、

$B$

はバイアス値（32ビット浮動小数点数の場合は127）です。 

2. 対数に変換して近似する 

逆平方根

$y=1/\sqrt{x}$

の両辺の2を底とする対数をとります。 

${\log }_2\left(y\right)={\log }_2\left(x^{-1/2}\right)=-\frac{1}{2}{\log }_2\left(x\right)$

ここで、浮動小数点数を整数にキャストする操作が、この対数変換を近似する役割を果たします。浮動小数点数を整数として見なすと、その値は指数部と仮数部を使って以下のように近似できます。 

$I_x\approx L\cdot {\log }_2\left(x\right)+C$

ここで、

$I_x$

は

$x$

を整数として解釈した値、

$L=2^{23}$

は仮数部のビット数、

$C$

は定数です。
同様に、

$y$

についても以下の式が成り立ちます。 

$I_y\approx L\cdot {\log }_2\left(y\right)+C$

これらの関係式を最初の対数式に代入して整理すると、以下のような関係式が導かれます。 

$I_y\approx \frac{3}{2}L(B-\sigma )-\frac{1}{2}{I}_x$

ここで、

$B=127$

（バイアス値）、

$\sigma$

は、対数近似の誤差を最小化するために導入された補正項です。 

3. マジックナンバーの導出 

導出された関係式は、次のようなビット操作に相当します。 

$I_y\approx R-(I_x\gg 1)$

ここで、>> 1は1ビット右シフトで、

$\frac{1}{2}{I}_x$

を意味します。

$R$

は

$\frac{3}{2}L(B-\sigma )$

の部分であり、これがマジックナンバーの候補となります。 

$L=2^{23}$

と

$B=127$

を代入すると、 

$R\approx \frac{3}{2}\cdot 2^{23}(127-\sigma )$

この式から、

$\sigma$

の値を調整することで、R の値が変わります。 

• 数式から導かれる理論値:
  

  $\sigma$
  の最適値として0.0450466などが導かれます。
• 経験的調整: 最終的なマジックナンバー0x5f3759dfは、すべての浮動小数点数における相対誤差を最小化するように、経験的に微調整されました。 

4. ニュートン法による最適化 

このビット操作による初期近似値は、ニュートン法を1回適用することで、さらに精度が向上します。

$y=y\cdot (1.5f-x_{half}\cdot y\cdot y)$

この一連のステップにより、当時の一般的な方法よりもはるかに高速に逆平方根を計算することが可能になったのです。 

結論 

マジックナンバー0x5f3759dfは、偶然発見されたものではなく、浮動小数点数の内部表現と対数の性質を巧みに利用した数式から、厳密に導出された値です。そして、その精度を最大化するために、多数の浮動小数点数でテストされて最適化された、理論と実践が融合した結果なのです。