MinGW GCC 4.9.1和浮点确定性

作者: 枫
发布时间: 2024-12-31 09:39:12 (1月前)
转自：

3 条回复

0#
回复此人
时间 | 2019-08-31 10-32

<div class =“post-text”itemprop =“text”> <P> 浮点数的精度存在差异，具体取决于存储位置。如果编译器将一个变量保存在寄存器中，则它作为存储在内存中的变量具有更高的精度。您可以尝试强制将变量存储在内存中，例如： </p> <pre> <code> int main() { std::array<double, 3u> arr = { 4.0, -2.0, 6.0 }; volatile double v1 = norm(arr); volatile double v2 = norm(arr); std::cout << v1 - v2 << '\n'; } </code> </pre> <P> 这样可以得到0的预期结果。 </p> </DIV>

编辑
1#
回复此人
猫南北 | 2019-08-31 10-32

<div class =“post-text”itemprop =“text”> <P> 正如@MatthiasB指出的那样，这似乎是gcc暂时将80位浮点值存储到64位寄存器/存储器位置的问题。请考虑以下简化程序仍然可以重现该问题： </p> <pre> <code> #include <cmath> #include <iostream> double norm() { double res = 4.0 * 4.0 + (-2.0 * -2.0) + (6.0 * 6.0); return std::sqrt(res); } int main() { std::cout << norm() - norm() << '\n'; return 0; } </code> </pre> <P> 必不可少的部分的汇编代码 <code> norm() - norm() </code> 看起来像这样（使用32位mingw 4.8.0编译器） </p> <pre> <code> ... call __Z4normv ; call norm() fstpl -16(%ebp) ; store result (80 bit) in temporary (64 bit!) call __Z4normv ; call norm() again fsubrl -16(%ebp) ; subtract result (80 bit) from temporary (64 bit!) ... </code> </pre> <P> 从本质上讲，我会认为这是一个gcc bug，但它似乎是一个 <a href="https://gcc.gnu.org/bugzilla/show_bug.cgi?id=12331" rel="noreferrer"> 复杂的话题 </A> ... </p> </DIV>

编辑

登录后才能参与评论