评论

#曾登高 发表于2006-03-20 13:12:00  IP: 211.100.21.*

TrackBack来自《上周技术关注：第16届JOLT大奖获奖名单公布》：

BOOKS GENERAL
Jolt Winner: Prefactoring by Ken Pugh (O'Reilly)

Productivity Winners:
* Innovation Happens Elsewhere: Open Source as Business Strategy by Ron Goldman, Richard P. Gabriel (Morgan Kaufmann)
* Producing Open Source Software: How to Run a...

#太假了吧 发表于2006-03-31 14:23:00  IP: 210.22.82.*

用到这地步，服了

#蜗牛 发表于2006-04-01 22:06:00  IP: 218.1.134.*

就是按照上边的函数算吗？我按照这样算出来不对呀，是不是我漏了什么呢？

#Zhe 发表于2006-04-01 23:04:00  IP: 218.80.248.*

这个标题有误导性，这个函数是求平方根的倒数的
return y; 改成 return number * y; 求的就是平方根了

#wangjun 发表于2006-04-06 17:37:00  IP: 61.172.57.*

一个比一个厉害！　服了！

#yanx8844 发表于2006-04-07 22:09:00  IP: 211.143.181.*

国人有没有这种程度的？

#咋回事？ 发表于2006-04-10 09:54:00  IP: 202.204.57.*

结果根本就不对！

#yyxxh 发表于2006-05-10 22:56:00  IP: 58.51.105.*

没人用数值计算方法分析一下吗？

#确实是计算出了平方根的倒数。 发表于2006-05-13 23:54:00  IP: 221.232.32.*

//测试代码如下：
#include<iostream>
using namespace std;

float Q_rsqrt( float number );

void main()
{
float y;
y = Q_rsqrt(4);
printf( "y = %f", y);
}

float Q_rsqrt( float number )
{
long i;
float x2, y;
const float threehalfs = 1.5F;

x2 = number * 0.5F;
y = number;
i = * ( long * ) &y; // evil floating point bit level hacking
i = 0x5f3759df - ( i >> 1 ); // what the fuck?
y = * ( float * ) &i;
y = y * ( threehalfs - ( x2 * y * y ) ); // 1st iteration
// y = y * ( threehalfs - ( x2 * y * y ) ); // 2nd iteration, this can be removed

#ifndef Q3_VM
#ifdef __linux__
assert( !isnan(y) ); // bk010122 - FPE?
#endif
#endif

// add by huanliang89@hotmail.com :watch out here!
return number*y;
}

#yl 发表于2006-06-23 14:17:00  IP: 218.249.66.*

服了。
都是超强人。

#chinadev 发表于2006-06-27 10:10:00  IP: 210.83.202.*

这个函数不是卡玛克写的，他曾经承认过~~~~

#jianburufei9133 发表于2007-12-03 21:23:10  IP: 218.67.135.*

各位大侠，请问“const float threehalfs = 1.5F;“中的"F"代表什么？谢谢。

2007-12-04 10:11:35作者回复

1.5F表示一个类型为const float的浮点数1.5。C里1.5这类常数的类型是const double。如果你不加F，编译器会认为你把一个const double赋给一个const float，多半会生成警告信息。