无限循环heisenbug：如果我添加打印输出，它就会退出

这是我的源代码：

#include <iostream>
#include <cmath>

using namespace std;

double up = 19.0 + (61.0/125.0);
double down = -32.0 - (2.0/3.0);
double rectangle = (up - down) * 8.0;

double f(double x) {
    return (pow(x, 4.0)/500.0) - (pow(x, 2.0)/200.0) - 0.012;
}

double g(double x) {
    return -(pow(x, 3.0)/30.0) + (x/20.0) + (1.0/6.0);
}

double area_upper(double x, double step) {
    return (((up - f(x)) + (up - f(x + step))) * step) / 2.0;
}

double area_lower(double x, double step) {
    return (((g(x) - down) + (g(x + step) - down)) * step) / 2.0;
}

double area(double x, double step) {
    return area_upper(x, step) + area_lower(x, step);
}

int main() {
    double current = 0, last = 0, step = 1.0;
    
    do {
        last = current;
        step /= 10.0;
        current = 0;
        
        for(double x = 2.0; x < 10.0; x += step) current += area(x, step);
        
        current = rectangle - current;
        current = round(current * 1000.0) / 1000.0;
        //cout << current << endl;
    } while(current != last);
    
    cout << current << endl;
    return 0;
}

它的作用是计算曲线之间的面积。 main() 中有一个循环 - 其目的是计算尽可能精确到小数点后 3 位的值。

它不起作用。为了调试的目的，我添加了唯一被注释的行。我想知道循环内部发生了什么。

//cout << current << endl;

当该行存在时 - 当我取消注释时 - 一切正常。如果不是的话 - 循环似乎是无限的。

神圣编译器，为什么？

据我所知，这不是不精确的浮点数的问题。循环内容在 4 次重复内完成when我正在输出其中的当前值。

@Skizz 的评论给出了可能的问题，但要详细说明：

浮点数学很棘手，特别是经常会出现舍入错误。诸如 1/1000.0 之类的数字（您的结果roundcall) 无法用浮点数精确表示。

更复杂的是，一方面是速度，另一方面是一致、直观的结果。例如，Intel 处理器的 FPU 以 80 位扩展精度格式存储值，而 C/C++double通常为 64 位。为了提高性能，编译器可能会将值保留在 FPU 中，作为 80 位临时值，尽管这可能会产生与将它们截断为 64 位时所得到的结果不同的结果。

启用调试语句后，current可能会存储到内存中，将其截断为 64 位，这样就可以与last.

禁用调试语句后，current可能是存储在 FPU 寄存器中的 80 位值，因此它可以never equal last， 只要last是一个 64 位值，它们都试图存储不精确的浮点表示x/1000.0.

解决方案是使用一个浮点比较与一些允许的误差（因为直接检查与浮点数的相等性几乎从来都不是一个好主意）。

进一步说明：我还没有查看汇编输出来验证情况是否如此；如果你愿意的话，你可以自己做。只有启用优化才能重现该问题。您也许可以通过调整编译器标志以选择一致性而不是速度来“修复”该错误，但正确的解决方案是使用不精确的比较而不是直接检查相等性。