你想快速計算某數的平方根倒數。
在 Quake Ⅲ Arena 的源代碼中，這個奇怪的算法對一個幻數進行整數運算，來計算平方根倒數的浮點近似值。
在 CoffeeScript 中，他使用經典原始的變量，以及由 Chris Lomont 發(fā)現的新的最優(yōu) 32 位幻數。除此之外，還使用 64 位大小的幻數。
另一特征是可以通過控制牛頓迭代法的迭代次數來改變其精確度。
相比于傳統(tǒng)的，該算法在性能上更勝一籌，這歸功于使用的機器及其精確度。
運行的時候使用 coffee -c script.coffee 來編譯 script：
然后復制粘貼編譯的 JS 代碼到瀏覽器的 JavaScript 控制臺。
注意：你需要一個支持類型數組的瀏覽器
參考：
以下的代碼來源于：https://gist.github.com/1036533
###
Author: Jason Giedymin <jasong _a_t_ apache -dot- org>
http://www.jasongiedymin.com
https://github.com/JasonGiedymin
在 Quake Ⅲ Arena 的源代碼 [1](ftp://ftp.idsoftware.com/idstuff/source/quake3-1.32b-source.zip) 中，這個奇怪的算法對一個幻數進行整數運算，來計算平方根倒數的浮點近似值 [5](http://en.wikipedia.org/wiki/Fast_inverse_square_root)。
在 CoffeeScript 中，我使用經典原始的變量，以及由 Chris Lomont [2](http://www.lomont.org/Math/Papers/2003/InvSqrt.pdf) 發(fā)現的新的最優(yōu) 32 位幻數。除此之外，還使用 64 位大小的幻數。
另一特征是可以通過控制牛頓迭代法 [3](http://en.wikipedia.org/wiki/Newton%27s_method) 的迭代次數來改變其精確度。
相比于傳統(tǒng)的，該算法在性能上更勝一籌，歸功于使用的機器及其精確度。
運行的時候使用 coffee -c script.coffee 來編譯 script：
然后復制粘貼編譯的 JS 代碼到瀏覽器的 JavaScript 控制臺。
注意：你需要一個支持類型數組 [4](https://developer.mozilla.org/en/JavaScript_typed_arrays) 的瀏覽器
###
approx_const_quake_32 = 0x5f3759df # See [1]
approx_const_32 = 0x5f375a86 # See [2]
approx_const_64 = 0x5fe6eb50c7aa19f9 # See [2]
fastInvSqrt_typed = (n, precision=1) ->
# 使用類型數組。現在只能在瀏覽器中操作。
# Node.JS 的版本即將推出。
y = new Float32Array(1)
i = new Int32Array(y.buffer)
y[0] = n
i[0] = 0x5f375a86 - (i[0] >> 1)
for iter in [1...precision]
y[0] = y[0] * (1.5 - ((n * 0.5) * y[0] * y[0]))
return y[0]
### 單次運行示例
testSingle = () ->
example_n = 10
console.log("Fast InvSqrt of 10, precision 1: #{fastInvSqrt_typed(example_n)}")
console.log("Fast InvSqrt of 10, precision 5: #{fastInvSqrt_typed(example_n, 5)}")
console.log("Fast InvSqrt of 10, precision 10: #{fastInvSqrt_typed(example_n, 10)}")
console.log("Fast InvSqrt of 10, precision 20: #{fastInvSqrt_typed(example_n, 20)}")
console.log("Classic of 10: #{1.0 / Math.sqrt(example_n)}")
testSingle()
掃碼關注鍍金池微信號
手機訂課，優(yōu)惠多多