函数sqrt,asin,acos源代码[备忘]

2015-12-25 admin 备忘

/*
 *sqrt源代码
 */
double my_sqrt(double x)
{
 register double ret;
 __asm__(
  "fsqrt"
  : "=t" (ret) 
  : "0" (x)
  );
 return ret;
}//计算x的平方根
/*
 *asin源代码
 */
double atan2 (double x, double y)
{
 register double ret;
 __asm__(
  "fpatan\n\t"
  "fld %%st(0)"
  : "=t" (ret) 
  : "0" (y), "u" (x)
  );
 return ret;
}//求x / y的反正切值。
double my_asin(double x)
{
 return atan2 (x, sqrt (1.0 - x * x));
}//求x的反正弦值。
/*
*acos源代码
*/
double atan2 (double x, double y)
{
 register double ret;
 __asm__(
  "fpatan\n\t"
  "fld %%st(0)"
  : "=t" (ret) 
  : "0" (y), "u" (x)
  );
 return ret;
}//求x / y的反正切值。
double my_acos(double x)
{
 return atan2 (sqrt (1.0 - x * x), x);
}//求x的反余弦值。

 

atan的泰勒公式 arctanx=∫(0→x)1/(1+x^2)dx=x-x^3/3+x^5/5+……+(-1)^n*x^(2n+1)/(2n+1) 


/***********************/ 以上是你要求的源代码,如果你要自己来实现,可以参考数值计算或数值逼近 sqrt:可以使用牛顿迭代公式,f(x)=x^2 - a 求零点(即方程f(x)=0的根),就是开方了。 而对于后面的两个,sin cos exp 是用泰勒公式和麦克劳林公式来计算,通过这里演化迭代就可以了,对于任何的c还是c++的编译器,在库函数里面都应该有这些基本的数学函数的,你在使用的时候加上#include<math.h>就可以了

标签: 算法 C++

发表评论:

Powered by emlog | |

浙ICP备17021512号 |浙公网安备 33010602008237号