为何开启编译器优化会改变double转unsigned int的标量转换结果？

测试代码 #

#include <stdio.h>
#include <xmmintrin.h>

int main ()
{
    double d[2] = { 161785254.0, -215713672.0 };
    int i[4];
    unsigned u[4];

    __m128d vd = _mm_loadu_pd(d);
    __m128i vi = _mm_cvtpd_epi32(vd);
    _mm_storeu_si128((void *)i, vi);
    _mm_storeu_si128((void *)u, vi);

    printf("Doubles to convert            : %.1f %.1f\n", d[0], d[1]);
    printf("\n");
    printf("SIMD double to signed int     : %i   %i\n", i[0], i[1]);
    printf("SIMD double stored to unsigned: %u   %u\n", u[0], u[1]);
    printf("\n");
    printf("Scalar double to signed int   : %i   %i\n", (int)d[0], (int)d[1]);
    printf("Scalar double to unsigned     : %u   %u\n", (unsigned)d[0], (unsigned)d[1]);
}

不同编译选项的输出 #

-O0编译输出（符合预期）

Doubles to convert            : 161785254.0 -215713672.0

SIMD double to signed int     : 161785254   -215713672
SIMD double stored to unsigned: 161785254   4079253624

Scalar double to signed int   : 161785254   -215713672
Scalar double to unsigned     : 161785254   4079253624

-O1编译输出（标量转unsigned结果异常）

Scalar double to unsigned     : 161785254   0

核心原因：未定义行为导致的优化差异 #

• C语言标准的明确规定：将负数的浮点数转换为无符号整数属于未定义行为。标准仅保证当浮点数的值在目标无符号类型范围内（≥0且≤unsigned类型最大值）时，转换结果确定；超出范围的转换，编译器可自由选择处理方式。
• 无优化（-O0）时的行为：编译器直接遵循硬件指令的默认逻辑：先把double转为signed int（得到-215713672），再将该有符号整数转为unsigned int——这一步会执行模2^32运算，最终得到4079253624（-215713672 + 2^32 = 4079253624）。
• 开启-O1优化后的行为：编译器会利用未定义行为做优化。既然标准允许任意处理这种非法转换，编译器选择直接返回0——这是合理的优化选择，因为它可以假设程序不会触发这类未定义行为，或用最简洁的方式处理。
• SIMD转换不受影响的原因：_mm_cvtpd_epi32指令行为明确：直接将double截断为signed int并存储到内存；后续读取为unsigned int只是对内存二进制位做无符号解释，全程未触发“负数浮点数转unsigned int”的未定义行为，因此结果始终一致。

内容的提问来源于stack exchange，提问作者aganm