c-为什么QNAN==QNAN不会引发FE_INVALID异常

Question

代码(t125.c(：

#include <fenv.h>
#include <stdint.h>
#include <stdio.h>
#if _MSC_VER
#pragma fenv_access (on)
#else
#pragma STDC FENV_ACCESS ON
#endif
void show_fe_exceptions(void)
{
printf("exceptions raised: ");
if (fetestexcept(FE_DIVBYZERO))     printf(" FE_DIVBYZERO");
if (fetestexcept(FE_INEXACT))       printf(" FE_INEXACT");
if (fetestexcept(FE_INVALID))       printf(" FE_INVALID");
if (fetestexcept(FE_OVERFLOW))      printf(" FE_OVERFLOW");
if (fetestexcept(FE_UNDERFLOW))     printf(" FE_UNDERFLOW");
if (fetestexcept(FE_ALL_EXCEPT)==0) printf(" none");
printf("n");
}
typedef union { uint32_t u; float f; } ufloat;
int main(void)
{
_Bool b;
ufloat uqnan;
volatile float f;
uqnan.u = 0x7fc00000;
f = uqnan.f;
b = f == f;
show_fe_exceptions();
return b ? 1 : 0;
}

调用：

$ gcc t125.c -Wall -Wextra -pedantic -std=c17 && ./a.exe
t125.c:7: warning: ignoring ‘#pragma STDC FENV_ACCESS’ [-Wunknown-pragmas]
7 | #pragma STDC FENV_ACCESS ON
|
exceptions raised:  none
$ clang t125.c -Wall -Wextra -pedantic -std=c17 && ./a.exe
t125.c:7:14: warning: pragma STDC FENV_ACCESS ON is not supported, ignoring pragma [-Wunknown-pragmas]
#pragma STDC FENV_ACCESS ON
^
1 warning generated.
exceptions raised:  none
$ cl t125.c /fp:strict /std:c17 && t125
exceptions raised:  none

问题：为什么QNAN == QNAN不导致引发FE_INVALID异常？

UPD。问题原因：(false，见下文(假设<any_NAN> == <any_NAN>导致引发FE_INVALID异常。

UPD2.更改代码：从f = *(float*)&qnan更改为f = uqnan.f(通过并集类型punning(。这是为了避免违反C标准的混叠规则。

Answer 1

ISO/IEC 9899:2011(E((增加强调(：

5.2.4.2.2浮动类型的特征<float.h>

3安静的NaN几乎在每个算术运算中传播，而不会引发浮点异常；信号NaN通常在作为算术操作数出现时引发浮点异常。

另请参阅：安静NaN和信令NaN之间的区别是什么？。

UPD。是的，相等似乎不算为此目的进行的算术运算。下面引用IEEE 754-2008的一句话(加上重点(：

5.11比较谓词的详细信息明确考虑安静NaN操作数可能性的程序可以使用表5.3中的无序安静谓词，不会发出此类无效操作异常的信号。

例如，谓词LT EQ不应导致引发无效操作异常。然而，我们看到(在上面的评论中(对于f <= f，gcc和cl (msvc)都提高了FE_INVALID。这是错误/缺陷吗？尽管它们都没有用定义1来定义CCD_ 12。但是，它们不是必需的:

__STDC_IEC_559__整数常数1，旨在表示符合附录F(IEC 60559浮点运算(中的规范。