Delphi和MSVC不会以相同的方式将+NAN与零进行比较

我正在将 C 代码移植到 Delphi，并发现编译器（Delphi 10.4.1 和 MSVC2019，均针对 x32 平台）处理 +NAN 与零的比较的方式存在问题。两种编译器都使用 IEEE754 表示双浮点值。我发现了这个问题，因为我移植到 Delphi 的 C 代码附带了一堆数据来验证代码的正确性。

原始源代码很复杂，但我能够在 Delphi 和 C 中生成最小的可重现应用程序。

C代码：

#include <stdio.h>
#include <math.h>

double AngRound(double x) {
    const double z = 1 / (double)(16);
    volatile double y;
    if (x == 0) 
        return 0;
    y = fabs(x);
    /* The compiler mustn't "simplify" z - (z - y) to y */
    if (y < z)
        y = z - (z - y);      // <= This line is *NOT* executed
    if (x < 0)
        return -y;
    else
        return y;             // <= This line is executed
}

union {
    double d;
    int bits[2];
} u;


int main()
{
    double lon12;
    double ar;
    int    lonsign;

    // Create a +NAN number IEEE754
    u.bits[0] = 0;
    u.bits[1] = 0x7ff80000;

    lon12    = u.d;                // Debugger shows lon12 is +nan
    if (lon12 >= 0)
        lonsign = 1;
    else
        lonsign = -1;              // <= This line is executed
    // Now lonsign is -1

    ar = AngRound(lon12);
    // Now ar is +nan

    lon12 = lonsign * ar;
    // Now lon12 is +nan
}

德尔福代码：

program NotANumberTest;

{$APPTYPE CONSOLE}

{$R *.res}

uses
  System.SysUtils;

type
    TRec = record
       case t : Boolean of
       TRUE:  (d    : Double);
       FALSE: (bits : array [0..1] of Cardinal);
    end;

function AngRound(x : Double) : Double;
const
    z : Double = 1 / Double(16);
var
    y : Double;
begin
    if x = 0 then
        Result := 0
    else begin
        y := abs(x);
        if y < z then
            // The compiler mustn't "simplify" z - (z - y) to y
            y := z - (z - y);           // <= This line is executed
        if x < 0 then
            Result := -y                // <= This line is executed
        else
            Result := y;
    end;
end;

var
    u       : TRec;
    lon12   : Double;
    lonsign : Integer;
    ar      : Double;
begin
    // Create a +NAN number IEEE754
    u.bits[0] := 0;
    u.bits[1] := $7ff80000;

    lon12 := u.d;                       // Debugger shows lon12 is +NAN
    if lon12 >= 0 then
        lonsign := 1                    // <= This line is executed
    else
        lonsign := -1;
    // Now lonsign is +1

    ar := AngRound(lon12);
    // Now ar is -NAN

    lon12 := lonsign * ar;
    // Now lon12 is -NAN
end.

我已经标记了比较后执行的行。当 lon12 变量等于 +NAN 时，Delphi 将 (lon12 >= 0) 评估为 TRUE。当 lon12 变量等于 +NAN 时，MSVC 将 (lon12 >= 0) 评估为 FALSE。

lonsign 在 C 和 Delphi 中有不同的值。

AngRound 接收 +NAN 作为参数返回不同的值。

lon12 的最终值是（致命的）不同的。

编译器生成的机器码是不同的：

Delphi生成的机器码：

MSVC2019生成的机器码：

Delphi 中的比较结果似乎更合乎逻辑：当 lon12 为 +NAN 时，(lon12 >= 0) 为 TRUE。这是否意味着错误在 MSVC2019 编译器中？我应该考虑原始 C-Code 的测试数据集进位错误吗？

回答

{$APPTYPE CONSOLE}

uses
  System.Math;

var
  d: Double;
begin
  SetFPUExceptionMask(exAllArithmeticExceptions);
  SetSSEExceptionMask(exAllArithmeticExceptions);
  d := NaN;
  Writeln(d > 0);
  Writeln(d >= 0);
  Writeln(d < 0);
  Writeln(d <= 0);
  Writeln(d = d);
  Writeln(d <> d);
end.

真的
真的
错误的
错误的
真的
错误的

真的
真的
真的
真的
错误的
真的

错误的
错误的
错误的
错误的
错误的
真的