macOS不同架构下clang编译C程序结果不一致问题排查

问题背景 #

将复杂程序移植到macOS/arm64平台时，发现程序行为与macOS/x86_64存在差异，用以下简化C代码复现问题（使用conda-forge提供的原生osx/arm64 clang 14.0.6，交叉编译x86_64版本）：

#include "assert.h"
#include "stdio.h"
int main()
{
    double y[2] = {-0.01,0.9};
    double r;
    r = y[0]+0.03*y[1];
    printf("r = %24.26e\n",r);
    assert(r == 0.017);
}

运行结果 #

• arm64平台：

$ clang -arch arm64 test.c -o test; ./test
Assertion failed: (r == 0.017), function main, file test.c, line 9.
r = 1.69999999999999977517983751e-02
zsh: abort      ./test

• x86_64平台（交叉编译后通过Rosetta运行或原生运行）：

$ clang -arch x86_64 test.c -o test; ./test
r = 1.70000000000000012212453271e-02
$       

已完成的排查 #

1. 验证两个平台的数值约定一致（舍入模式、精度等均符合IEEE754标准）：

#include <iostream>
#include <limits>

#define LOG(x) std::cout << #x " = " << x << '\n'

int main()
{
    using l = std::numeric_limits<double>;
    LOG(l::digits);
    LOG(l::round_style);
    LOG(l::epsilon());
    LOG(l::min());

    return 0;
}

输出：

l::digits = 53
l::round_style = 1
l::epsilon() = 2.22045e-16
l::min() = 2.22507e-308

2. arm64平台下表达式计算结果与调用refBLAS ddot（向量{1,0.03}与y）的结果一致。
3. 使用macOS 11.6.1默认Apple clang 13.0.0编译时，两个架构结果一致，锁定问题源于conda包cxx-compiler 1.5.2工具链。

可能导致差异的编译标志 #

1. 浮点运算优化类标志 #

• -ffast-math：该标志会启用一系列违反IEEE754标准的优化，包括允许重新排列浮点运算顺序、忽略NaN/无穷大检查等。conda-forge的clang可能默认开启了此标志，而Apple clang默认未启用。不同架构下运算重排的逻辑不同，会直接导致最终计算结果差异。
• -funsafe-math-optimizations：作为-ffast-math的子集，允许编译器进行可能改变浮点结果的优化（如合并运算项、重新关联表达式），arm64和x86_64的优化路径差异会引发数值偏差。
• -fassociative-math：允许编译器修改浮点运算的结合顺序（如a+(b+c)改为(a+b)+c），由于浮点运算不满足严格结合律，不同架构的优化策略会产生不同的中间结果，最终导致输出差异。

2. 架构特定的指令集标志 #

• -march=native：conda-forge的clang可能默认针对arm64启用了本地指令集优化（如使用NEON高级浮点指令），而x86_64交叉编译时的指令集配置不同。这些指令在处理中间结果的精度和舍入行为上存在细微差异。
• -mfpu/-mfloat-abi：针对arm架构的浮点单元和ABI配置，不同设置会影响浮点运算的精度和计算方式，与x86_64的对应参数不匹配时会引发结果差异。

3. 编译优化级别 #

• -O2/-O3：较高优化级别下，不同架构的编译器对浮点运算的优化策略可能不同。比如arm64版本可能将y[0]+0.03*y[1]优化为特定指令序列，而x86_64版本的优化路径不同，最终导致数值偏差。conda-forge工具链可能默认使用了比Apple clang更高或不同的优化级别。

4. 浮点环境控制标志 #

• -frounding-math：该标志会让编译器避免依赖默认舍入模式进行优化，若conda-forge的clang在不同架构下对该标志的处理不一致，可能导致舍入行为差异。不过已验证舍入模式一致，此可能性较低，但仍可排查。

排查建议 #

• 显式禁用浮点优化标志：编译时添加-fno-fast-math、-fno-associative-math，观察两个架构的结果是否一致。
• 对比默认编译参数：用clang -arch arm64 -v test.c查看conda-forge clang的完整编译参数，再与Apple clang的参数对比，找出差异项。
• 调整优化级别：从-O0（无优化）开始逐步测试，确认是否是优化级别导致的差异。

内容的提问来源于stack exchange，提问作者Stéphane Mottelet