sqrt函数实现分析

看到一篇博文讲述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