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
",int_sqrt(16))	;
	return 0 ;
}

执行结果：