看到一篇博文讲述sqrt的实现方法,其中二分法没有区分区间,当x位于(0,1)和(1, 正无穷) 时初始值设置应该有所区别
const float eps=0.000001; // eps的值可能影响最后计算精度,甚至导致无限循环
// 二分法,注意区分x的取值区间
float SqrtByBisection(float x)
{
if(x<0) // 负数
return x;
if(x<=eps) // 正数0
return 0.0f;
if(fabs(x-1)<=eps) // 正数1
return 1.0f;
float left, right;
float mid;
if(x>eps&&x<1.0f-eps) // (0,1)区间
{
left=x;
right=1.0f;
}
else
{
left=1.0f;
right=x;
}
while(right-left>eps) // (1,)区间
{
mid=(left+right)/2;
if(mid*mid>x+eps)
right=mid;
else if(mid*mid<x-eps)
left=mid;
else
return mid;
}
return mid;
}
上述方法虽然考虑了区间,但是验证时还是出现细节问题,如sqrt(0.25)=0.499999等
与二分法相比,牛顿迭代法速度更快,计算也更准确
// 牛顿迭代法
float SqrtByNewton(float x)
{
float val=x;
float last;
while(fabs(val-last)>eps)
{
last=val;
val=(val+x/val)/2;
}
return val;
}
更快的方法,原来用在一款游戏中,参加:http://zh.wikipedia.org/wiki/%E5%B9%B3%E6%96%B9%E6%A0%B9%E5%80%92%E6%95%B0%E9%80%9F%E7%AE%97%E6%B3%95