函数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++

叶子的小屋

好记性不如多写点笔记

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

搜索

天气

分类

标签

存档

链接