使用Clang -Ofast编译时isinf()与isfinite()失效的解决办法及疑问

一、std::isinf()无法识别inf的原因 #

核心问题出在**-Ofast编译选项带来的编译器优化行为**：

• -Ofast是Clang的最高级优化选项，它会启用-ffast-math集合，其中包含-fno-math-errno、-fno-trapping-math等子选项。编译器会基于“浮点运算不会产生NaN、inf等异常值”的假设做激进优化。
• 基于这个假设，编译器会直接将std::isinf(tmp)的调用优化为常量0——哪怕实际运行时生成了inf，它也认为这种情况不可能发生。
• 位操作能识别inf，是因为它直接读取浮点值的二进制表示，不依赖标准库的数学函数，不会被编译器的优化逻辑干扰。

可行解决方法 #

如果你不想用位操作，可以通过以下方式规避：

• 调整优化选项：将-Ofast替换为-O3 -ffast-math -fmath-errno，保留大部分优化效果的同时，让编译器保留浮点错误的处理逻辑，这样std::isinf就能正常工作。
• 局部禁用优化：用编译器指令包裹检测代码块，强制编译器不对这部分代码做激进优化：

  #pragma clang optimize off
  std::cout << "??? " << tmp << " " << std::isinf(tmp) << std::endl;
  #pragma clang optimize on
  

二、负指数std::pow()生成inf的原因 #

你可能混淆了运算逻辑：std::pow(U(2), -I * wbits)等价于1.0 / std::pow(U(2), I * wbits)。
当I * wbits的值足够大时：

• std::pow(U(2), I * wbits)会因指数过大，超出浮点类型U的最小可表示非零值范围，直接下溢为0。
• 此时1.0 / 0在浮点运算中会生成正inf，再乘以非0、非NaN的f，结果自然还是inf。

比如对于float类型，当I*wbits > 127时，2^(I*wbits)会下溢为0，最终导致pow(2, -I*wbits)变成inf。

内容的提问来源于stack exchange，提问作者user1095108