计算机与代数---如何计算sqrt---优化

0.简介

根据上一篇的内容，今天又看到了比较的好的一个优化方法，在这里记录一下。

1.判断条件优化

之前的二分法和牛顿迭代法的判断条件都有乘法，今天发现一个比较的好的判断方法。

二分法

long double l = 0, r = (n>=1 ? n:1), mid=0;
while (fabs(l-r)>1e-8)
{
	mid = (l + r) / 2;
	if (mid * mid > n) r = mid;
	else l = mid;
}

牛顿迭代法

long double x = n,y = -1,k,b,lx = n/2;
while (fabs(lx - x)>1e-6)//计算两次x的相近程度
{
	lx = x;//保存x
	y = x * x - n;
	k = 2 * x;
	b = y - k * x;
	x = -(b / k);//新的x
}

这样就没有了乘法在判断条件中。

牛顿迭代法的代码再简单整理下。

long double x = n,y = -1,k,b,lx = n/2;
while (fabs(lx - (x = ((lx = x) * x + n) / (2 * x))) > 1e-6);