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