float Sqrt(float x);
float InvSqrt(float x);
}
//实现平方根的倒数
float InvSqrt(float x)
{
float xhalf = 0.5f*x;
int i = *(int*)&x; // get bits for floating value
i = 0x5f375a86- (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 牛顿迭代法,次数越多越准确
x = x*(1.5f-xhalf*x*x); // Newton step, repeating increases accuracy
x = x*(1.5f-xhalf*x*x); // Newton step, repeating increases accuracy
return x;
}
原文地址:http://blog.csdn.net/u013467442/article/details/41147815