我正在尝试计算没有数学库的arkus sin()。它工作,但在不同的平台上产生不同的结果。
Windows和Mac(正确):5.2359877560e-001 Linux:5.1532736669e-01
问题在哪里?
double my_asin(double k)
{
double x = k;
double a = 1;
double s = x;
const double hodnota = 0.00000000001f;
double s1 = 0;
double d = x;
double e = 1.0;
double y;
y = x * x;
while(my_abs(s1 - s) >= hodnota) {
s1 = s;
d = d * y;
e = e * ((a * a) / ((++a) * (++a)) );
s = s + (d * e);
}
return s;
}
答案 0 :(得分:9)
e = e * ((a * a) / ((++a) * (++a)) );
可能会产生不同的结果
您需要以这种方式更改代码:
e *= a * a / ((a + 1) * (a + 2));
a += 2;
答案 1 :(得分:2)
此行具有未定义的效果。
e = e * ((a * a) / ((++a) * (++a)) );
你需要移动一个,理想情况下两个都可以移动到不同的行。例如:
e = e * ((a * a) / ((a+2) * (a+1)) );
a+=2
虽然你需要玩替换线,直到它实际上做你希望的。
顺便说一句,这段代码很容易因编译器版本而异,而不仅仅是编译器品牌和操作系统。