发一个吊的不行的开平方算法

247153481

发一个吊的不行的开平方算法 [复制链接]

 

FAST INVERSE SQUARE ROOT
看到这个名字有没有什么想法

这个可是号称世界上最快的浮点开方算法哦，怎么做到的？0x5f3759df就是它
原算法好像是在雷神之锤里面实现的。
大家看看吧

1. #include <math.h>
   
2. float InvSqrt(float x)
   
3. {
   
4. float xhalf = 0.5f*x;
   
5. int i = *(int*)&x; // get bits for floating VALUE
   
6. i = 0x5f375a86- (i>>1); // gives initial guess y0
   
7. x = *(float*)&i; // convert bits BACK to float
   
8. x = x*(1.5f-xhalf*x*x); // Newton step, repeating increases accuracy
   
9. return x;
   
10. }
    
11. 
    
12. int main()
    
13. {
    
14.   printf("%lf",1/InvSqrt(3));
    
15. 
    
16.    return 0;
    
17. }
    

复制代码

在附上论文
Fast_Inverse_Square_Root.pdf (195.54 KB, 下载次数: 36)

2015-8-25 14:27 上传

点击文件名下载附件

lyzhangxiang

我看过的版本注释是这样的

what's the fuck of 0x5f375a86 ??

lyzhangxiang

特地去找了下，号称比一般的算法快5倍

lcofjp

确实叼！

247153481

lyzhangxiang 发表于 2015-8-25 15:43
我看过的版本注释是这样的

what's the fuck of 0x5f375a86 ??

是啊，然后一个教授专门又去研究这个数，通过理论证明，后面好像有找到一个更好的数字，跟这个数字非常接近

lyzhangxiang

247153481 发表于 2015-8-25 16:39
是啊，然后一个教授专门又去研究这个数，通过理论证明，后面好像有找到一个更好的数字，跟这个数字非常接 ...

但是没这个数准。。

247153481

lyzhangxiang 发表于 2015-8-25 16:39
但是没这个数准。。

反正不知道是怎么想到的

wsmysyn

http://www.cnblogs.com/pkuoliver ... try-about-sqrt.html
好牛啊，好想看看里面都还有哪些NB的算法，可惜链接的网页里提供下载的都失效了

wsmysyn

特意找了一个QUAKE-III的原代码，找到了那个传说中的Q_math.c，里面真的有这个函数，好强大

白丁

看题目以为会是段汇编

chenzhufly

好牛逼！

cj_in

早有耳闻，今天见到真身了

chenzhufly

这个算法非常好！

woody_chen

niubility