C语言之linux内核实现平方根计算算法

关于平方根的计算，在linux内核中也有实现，就像math.h数学库里的sqrt这个函数一样。

平方根的公式定义:

如果一个非负数x的平方等于a，即

  

，

  

，那么这个非负数x叫做a的算术平方根。a的算术平方根记为

  

，读作“根号a”，a叫做被开方数（radicand）。求一个非负数a的平方根的运算叫做开平方。结论：被开方数越大，对应的算术平方根也越大（对所有正数都成立）。

一个正数如果有平方根，那么必定有两个，它们互为相反数。显然，如果我们知道了这两个平方根的一个，那么就可以及时的根据相反数的概念得到它的另一个平方根。 

哈哈，小学生都懂，不解释不解释，直接来看代码：

一样的，从内核里把代码取出来：

#include <stdio.h>
#ifdef CONFIG_64BIT
#define BITS_PER_LONG 64
#else
#define BITS_PER_LONG 32
#endif
/**
 * int_sqrt - rough approximation to sqrt
 * @x: integer of which to calculate the sqrt
 *
 * A very rough approximation to the sqrt() function.
 */
unsigned long int_sqrt(unsigned long x)
{
	unsigned long op, res, one;
	op = x;
	res = 0;
	one = 1UL << (BITS_PER_LONG - 2);
	while (one > op)
		one >>= 2;

	while (one != 0) {
		if (op >= res + one) {
			op = op - (res + one);
			res = res +  2 * one;
		}
		res /= 2;
		one /= 4;
	}
	return res;
}

int main(void)
{
	printf("%d\n",int_sqrt(16))	;
	return 0 ;
}

运行结果：