Delphi와 MSVC는 동일한 방식으로 + NAN을 0과 비교하지 않습니다.

C 코드를 Delphi로 이식 중이며 컴파일러 (Delphi 10.4.1 및 MSVC2019, 둘 다 x32 플랫폼을 대상으로 함)가 + NAN과 0의 비교를 처리하는 방식에 문제가 있습니다. 두 컴파일러 모두 이중 부동 소수점 값에 IEEE754 표현을 사용합니다. 델파이로의 C-Code I 이식이 코드 정확성을 검증하기 위해 많은 데이터와 함께 제공되기 때문에 문제를 발견했습니다.

원본 소스 코드는 복잡하지만 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과 같을 때 (lon12> = 0)을 TRUE로 평가합니다. MSVC는 lon12 변수가 + NAN과 같을 때 (lon12> = 0)을 FALSE로 평가합니다.

lonsign은 C와 Delphi에서 다른 값을 가지고 있습니다.

+ NAN을 인수로받는 AngRound는 다른 값을 반환합니다.

lon12의 최종 값은 (치명적으로) 다릅니다.

컴파일러에 의해 생성 된 기계 코드는 다릅니다.

델파이가 생성 한 기계 코드 :

MSVC2019 생성 기계 코드 :

비교 결과는 Delphi에서 더 논리적으로 보입니다. (lon12> = 0) lon12가 + NAN 일 때 TRUE입니다. 이것은 버그가 MSVC2019 컴파일러에 있음을 의미합니까? 원래 C-Code의 테스트 데이터 세트가 오류가있는 것으로 간주해야합니까?