开启-ffast-math时如何检测double类型的NaN值？

为什么常规isnan会失效 #

开启-ffast-math后，编译器会启用一系列浮点数优化，包括假设浮点数不存在NaN/Inf，这会导致标准库的isnan()宏或函数被优化为恒假，或者行为不符合预期。但正如你所说，基础运算生成的NaN依然会保留IEEE 754定义的比特位结构，所以可以通过其他可靠方法检测。

可靠的检测方案 #

1. 利用NaN的核心特性：自不等性 #

IEEE 754标准规定，NaN是唯一不等于自身的浮点数。即使开启-ffast-math，主流编译器（GCC/Clang/MSVC）通常不会优化掉x != x这个比较逻辑——不过为了防止极端情况下的优化，可以给变量加上volatile修饰：

int is_nan(double x) {
    volatile double tmp = x;
    return tmp != tmp;
}

这个方法完全跨平台，不需要依赖任何特定硬件或编译器特性，是最通用的 fallback 方案。

2. 基于IEEE 754比特位的安全位运算 #

绝大多数现代CPU都遵循IEEE 754双精度标准：double类型占8字节，其中最高1位是符号位，接下来11位是指数位，最后52位是尾数位。NaN的判定规则是：指数位全为1，且尾数位不全为0。

为了避免严格别名问题（直接用指针转换可能触发未定义行为），可以用union来访问double的比特位：

#include <stdint.h>

int is_nan(double x) {
    union {
        double d_val;
        uint64_t u_val;
    } converter = {.d_val = x};
    // 检查指数位全1，且尾数位非零
    const uint64_t exponent_mask = 0x7FF0000000000000ULL;
    const uint64_t mantissa_mask = 0x000FFFFFFFFFFFFFULL;
    return ((converter.u_val & exponent_mask) == exponent_mask) &&
           ((converter.u_val & mantissa_mask) != 0);
}

这个方法的兼容性极强，只要目标平台使用IEEE 754双精度格式（现在几乎所有平台都满足），就能稳定工作。

3. 编译器内置函数（推荐） #

主流编译器提供了专门的内置函数来检测NaN，即使开启-ffast-math也能正常工作，且性能最优：

#ifdef __GNUC__  // GCC、Clang、MinGW等
#define IS_NAN(x) __builtin_isnan(x)
#elif _MSC_VER   // MSVC
#define IS_NAN(x) _isnan(x)
#else
// fallback到自不等性检测
#define IS_NAN(x) ((volatile double)(x) != (x))
#endif

使用时直接调用IS_NAN(x)即可，这种方案兼顾了性能和跨平台兼容性。

总结 #

• 优先选择编译器内置函数，性能最高且最可靠；
• 自不等性方法是最通用的跨平台 fallback，无需依赖任何硬件特性；
• 位运算方法需要确认平台遵循IEEE 754，但现代平台几乎都满足，union方式避免了未定义行为。

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