c++ - 为什么使用英特尔 C++ 编译器时 NaN - NaN == 0.0？

Question

众所周知，NaN在算术中传播，但是我找不到任何演示，所以我写了一个小测试:

#include <limits>
#include <cstdio>

int main(int argc, char* argv[]) {
    float qNaN = std::numeric_limits<float>::quiet_NaN();

    float neg = -qNaN;

    float sub1 = 6.0f - qNaN;
    float sub2 = qNaN - 6.0f;
    float sub3 = qNaN - qNaN;

    float add1 = 6.0f + qNaN;
    float add2 = qNaN + qNaN;

    float div1 = 6.0f / qNaN;
    float div2 = qNaN / 6.0f;
    float div3 = qNaN / qNaN;

    float mul1 = 6.0f * qNaN;
    float mul2 = qNaN * qNaN;

    printf(
        "neg: %f\nsub: %f %f %f\nadd: %f %f\ndiv: %f %f %f\nmul: %f %f\n",
        neg, sub1,sub2,sub3, add1,add2, div1,div2,div3, mul1,mul2
    );

    return 0;
}

这个例子(running live here)基本上产生了我所期望的(否定有点奇怪，但有点道理):

neg: -nan
sub: nan nan nan
add: nan nan
div: nan nan nan
mul: nan nan

MSVC 2015 产生了类似的东西。但是，英特尔 C++ 15 会产生:

neg: -nan(ind)
sub: nan nan 0.000000
add: nan nan
div: nan nan nan
mul: nan nan

具体来说，qNaN - qNaN == 0.0.

这……不可能吧？相关标准(ISO C、ISO C++、IEEE 754)对此有何规定，为什么编译器之间的行为存在差异？

Answer 1

英特尔 C++ 编译器中的默认浮点处理是 /fp:fast，它不安全地处理 NaN(这也会导致 NaN == NaN 是 true 例如)。尝试指定 /fp:strict 或 /fp:precise 看看是否有帮助。