快速开方之魔法数 0x5f3759df

一、快速开方之魔法数 0x5f3759df

测试代码：

#include 
#include 
#include 
#include 

#define N_NUM 2147483640

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
	return y;
}

int main()
{
	int i, j;
	float res;

	struct timeval start, end;
	int timeuse;

	printf("%f\n", 1.0/Q_rsqrt(3.14159));
	printf("%f\n", sqrt(3.14159));
	printf("%f\n", 1.0/Q_rsqrt(100));
	printf("%f\n", sqrt(100));

	gettimeofday(&start, NULL);
	for(i=0; i

运行结果：

1.773184

1.772453

10.015536

10.000000

time1=0.000s

done

time2=3.439s //耗时差距很大