C++ exp10与pow不同（10）
首先，我意识到大多数以10为基数的数字不能精确地以2为基数表示，因此我的问题并不是关于浮点运算的缺陷 我正在尝试编写一个函数，通过检查最后6个有意义的数字是否在某个公差范围内，并将其更改为高于某个假定精确值的下一个可表示的值（仅用于显示目的-除非它是整数或二的幂），尝试更正受累积舍入错误影响的双精度 我的函数中让我吃惊的一个部分是exp10的输出；就我所知，只要两个double之间的间距小于2，那么存储为double的整数值应该是精确的——尽管10^14正在推它，但这应该是一个精确的整数（因为10^14=~2^46.507<2^53）。然而，这不是我的测试所显示的 我的调试工作摘录（没有什么比这更明显）和输出如下：C++ exp10与pow不同（10）,c++,floating-point,math.h,C++,Floating Point,Math.h,首先，我意识到大多数以10为基数的数字不能精确地以2为基数表示，因此我的问题并不是关于浮点运算的缺陷 我正在尝试编写一个函数，通过检查最后6个有意义的数字是否在某个公差范围内，并将其更改为高于某个假定精确值的下一个可表示的值（仅用于显示目的-除非它是整数或二的幂），尝试更正受累积舍入错误影响的双精度 我的函数中让我吃惊的一个部分是exp10的输出；就我所知，只要两个double之间的间距小于2，那么存储为double的整数值应该是精确的——尽管10^14正在推它，但这应该是一个精确的整数（因为1
double test = 0.000699;
double tmp = fabs(test);
double exp = 10.0 - floor(log10(tmp));
double powTen = exp10(10.0 - floor(log10(tmp)));
double powTen2 = exp10(exp);
double powTen3 = exp10((int)exp);
double powTen4 = exp10(exp);
double powTen5 = pow(10, exp);
printf("exp: %.16lf\n", exp);
printf("powTen: %.16lf\n", powTen);
printf("powTen2: %.16lf\n", powTen2);
printf("powTen3: %.16lf\n", powTen3);
printf("powTen4: %.16lf\n", powTen4);
//these two are exact
printf("10^14: %.16lf\n", exp10(14));
printf("powTen5: %.16lf\n", powTen5);
printf("exp == 14.0: %d\n", exp == 14.0);
exp: 14.0000000000000000
powTen: 100000000000000.1250000000000000
powTen2: 100000000000000.1250000000000000
powTen3: 100000000000000.1250000000000000
powTen4: 100000000000000.1250000000000000
10^14: 100000000000000.0000000000000000
powTen5: 100000000000000.0000000000000000
exp == 14.0: 1
编辑-通过O3、O2、O1优化，我将获得预期的输出-除非运行时才能知道数据。此时exp10仍有问题。您的
exp10doubleexp10-fno-builtinexp10eglibcexp10double
__ieee754_exp10 (double arg)
{
/* This is a very stupid and inprecise implementation. It'll get
replaced sometime (soon?). */
return __ieee754_exp (M_LN10 * arg);
}
不要期望它能很好地工作。谢谢，我已经关注了那篇文章，但是exp10的这种行为是不正确的-除非我对它的使用不正确-我不是在问为什么0.6看起来像0.599999999….+垃圾或为什么0.3-0.2-0.1！=0.0等等。。。我想问的是，为什么可以精确表示为整数的东西不是用exp10表示的，而是用pow
exp10（14）exp==14.0exp10exp10l（）doubleexp10exp10[1,2]double-doubledouble-doubledouble-double-doubleexp10