1

代码(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。问题的原因:(错误,见下文)<any_NAN> == <any_NAN>导致引发FE_INVALID异常的假设。

UPD2。更改代码:从f = *(float*)&qnanf = uqnan.f(通过联合键入双关语)。这是为了避免违反 C 标准的别名规则。

4

1 回答 1

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和raise 。这是一个错误/缺陷吗?尽管它们都没有定义为 1。但是,它们不是必需的:gcccl (msvc)FE_INVALID__STDC_IEC_559__

__STDC_IEC_559__整数常量 1,旨在表明符合附件 F(IEC 60559 浮点运算)中的规范。

于 2021-03-24T19:17:20.537 回答