我试图计算没有数学库的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;
}
解决方法:
指令e = e *((a * a)/((a)*(a)));由于undefined behavior,可以产生不同的结果
您需要以这种方式更改代码:
e *= a * a / ((a + 1) * (a + 2));
a += 2;
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。