C/C++代码技巧,那肯定是鼎鼎大名的
InvSqrt了,用来
快速计算平方根的倒数。
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
(以下摘自 0x5f375a86_百度百科)
这是一个传奇算法,此算法最早被认为是由 约翰·卡马克所发明,
但后来的调查显示,该算法在这之前就于计算机图形学的硬件与软件领域有所应用。
它的速度要比标准的牛顿迭代法快上 4 倍,
而其中的关键是一行神秘的代码和一个莫名其妙的数字 0x5f3759df。
普度大学的数学家Lomont觉得很好玩,
决定要研究一下卡马克弄出来的这个猜测值有什么奥秘。
Lomont也是个牛人,在精心研究之后从理论上也推导出一个最佳猜测值,
和卡马克的数字非常接近, 0x5f37642f。
卡马克真牛,他是外星人吗?
传奇并没有在这里结束。
Lomont计算出结果以后非常满意,于是拿自己计算出的起始值和卡马克的神秘数字做比赛,
看看谁的数字能够更快更精确的求得平方根。
结果是卡马克赢了...
谁也不知道卡马克是怎么找到这个数字的。
最后Lomont怒了,采用暴力方法一个数字一个数字试过来,
终于找到一个比卡马克数字要好上那么一丁点的数字,
虽然实际上这两个数字所产生的结果非常近似,
这个暴力得出的数字是 0x5f375a86。
Lomont为此写下一篇论文, Fast inverse square root (见wikipedia)
float InvSqrt(float x)
{
float xhalf = 0.5f*x;
int i = *(int*)&x; // get bits for floating value
i = 0x5f3759df - (i>>1); // gives initial guess y0
x = *(float*)&i; // convert bits back to float
x = x*(1.5f-xhalf*x*x); // Newton step, repeating increases accuracy
return x;
}
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
(以下摘自 0x5f375a86_百度百科)
这是一个传奇算法,此算法最早被认为是由 约翰·卡马克所发明,
但后来的调查显示,该算法在这之前就于计算机图形学的硬件与软件领域有所应用。
它的速度要比标准的牛顿迭代法快上 4 倍,
而其中的关键是一行神秘的代码和一个莫名其妙的数字 0x5f3759df。
普度大学的数学家Lomont觉得很好玩,
决定要研究一下卡马克弄出来的这个猜测值有什么奥秘。
Lomont也是个牛人,在精心研究之后从理论上也推导出一个最佳猜测值,
和卡马克的数字非常接近, 0x5f37642f。
卡马克真牛,他是外星人吗?
传奇并没有在这里结束。
Lomont计算出结果以后非常满意,于是拿自己计算出的起始值和卡马克的神秘数字做比赛,
看看谁的数字能够更快更精确的求得平方根。
结果是卡马克赢了...
谁也不知道卡马克是怎么找到这个数字的。
最后Lomont怒了,采用暴力方法一个数字一个数字试过来,
终于找到一个比卡马克数字要好上那么一丁点的数字,
虽然实际上这两个数字所产生的结果非常近似,
这个暴力得出的数字是 0x5f375a86。
Lomont为此写下一篇论文, Fast inverse square root (见wikipedia)
本文链接:http://www.zhihu.com/question/37692782/answer/73323805