我遇到了一个错误，我的变量与无穷大不等于无穷大。我都不确定是哪种情况。
下面是最小的例子：

#include <iostream>
#include <limits>
#include <cmath>

using namespace std;

constexpr auto NT_INF = std::numeric_limits<double>::infinity();

// prints a number byte-by-byte
void printN(double n){
    auto tmp = (uint8_t*)&n;
    for (int i = 0; i < 8; i++) {
        wprintf(L"%4d ", tmp[i]);
    }
    wprintf(L"\n");
};

double * genNumbers() {
    double * numbers = new double[6];
    
    numbers[0] = 100;
    numbers[1] = 1;
    numbers[2] = 198263783302;
    numbers[3] = 198263783302;
    numbers[4] = 100;
    numbers[5] = 0;
    
    return numbers;
}

int main()
{
    // leave no chance for the compiler to optimize variables away
    double * numbers = genNumbers();
    
    // floor((100. * 1.) / (198263783302. - 198263783302.)) / 100.;
    numbers[5] = floor((numbers[0] * numbers[1]) / (numbers[2] - numbers[3])) / numbers[4];
    double number = numbers[5];

    printN(number);
    printN(NT_INF);
    
    wprintf(L"%f, %f, %d, %d, %d, %d\n", number, NT_INF, number == 1.0 / 0.0, number == NT_INF, isinf(NT_INF), isinf(number));
    
    delete numbers;

    return 0;
}

我的GCC编译器标志是：-fno-reorder-blocks -ffast-math -march=native -fno-exceptions -Ofast -fmodulo-sched -fgcse-sm -fgcse-las -fno-inline-small-functions -flto=2 -fgcse-lm -fira-region=all
如果您尝试在onlinegdb中运行它，那么您将得到类似于

0    0    0    0    0    0  240  127 
   0    0    0    0    0    0  240  127 
inf, inf, 1, 1, 1, 0

所以NT_INF和number是相同的，从printN输出判断。同时，isinf说number不是无穷大。但printf说它是？？？
如果我尝试在本地编译，情况会更糟（我想，这是因为-march=native）。

0    0    0    0    0    0  240  127 
   0    0    0    0    0    0  240  127 
inf, inf, 0, 0, 1, 0

因此，这里的number变量甚至不等于无穷大，使用简单的==
我完全不知道这个例子（我的真实的代码）有什么问题。也许这是由于一些旗帜或他们的组合，但有太多的组合，以尝试他们。