Ubuntu不同版本下g++编译sin函数的浮点结果差异及对齐需求

我在Ubuntu 18.04 Bionic Docker镜像中开发的代码，部署到Ubuntu 20.04 Focal Docker镜像后单元测试失败，已定位到以下极简代码的差异：

#include <iostream>
#include <iomanip>
#include <math.h>
int main()
{
    const float DEG_TO_RAD_FLOAT = float(M_PI / 180.);
    float theta = 22.0f;
    theta = theta * DEG_TO_RAD_FLOAT;
    std::cout << std::setprecision(20) << theta << ' ' << sin(theta) << std::endl;
    return 0;
}

在Bionic镜像中，我通过以下命令升级了g++版本：

sudo apt-get install -y software-properties-common
sudo add-apt-repository ppa:ubuntu-toolchain-r/test
sudo apt install -y gcc-9 g++-9
sudo update-alternatives --install /usr/bin/gcc gcc /usr/bin/gcc-9 90 --slave /usr/bin/g++ g++ /usr/bin/g++-9 --slave /usr/bin/gcov gcov /usr/bin/gcov-9
sudo update-alternatives --install /usr/bin/gcc gcc /usr/bin/gcc-7 70 --slave /usr/bin/g++ g++ /usr/bin/g++-7 --slave /usr/bin/gcov gcov /usr/bin/gcov-7

两个环境的g++版本均为9.4.0，但输出结果存在差异：

• Ubuntu 18.04输出：0.38397243618965148926 0.37460657954216003418
• Ubuntu 20.04输出：0.38397243618965148926 0.37460660934448242188

可见sin(theta)的结果在第7位小数处不同，推测差异源于libc版本（18.04为2.27，20.04为2.31）。我尝试过-mfpmath=sse、-fPIC、-ffloat-store、-msse、-msse2等g++编译选项，但均无效果。

另外，使用/fp:precise编译的Windows版本代码，结果与Ubuntu 18.04一致。请问如何让g++编译结果与Windows编译器保持一致？

1. 替换libc的数学库实现 #

glibc不同版本的sinf（float版本的sin函数）实现存在精度优化差异，可通过以下方式规避：

• 切换到musl libc：musl的数学函数实现更注重跨平台一致性，可基于alpine等musl镜像编译代码，或在Ubuntu环境中静态链接musl库。
• 使用第三方数学库：替换系统sin函数为GSL（GNU Scientific Library）的gsl_sf_sin，或MPFR高精度浮点库函数，这类库的跨平台计算行为更稳定。

2. 对齐g++与Windows的浮点编译规则 #

Windows的/fp:precise会禁用激进浮点优化，严格遵循IEEE标准，g++可通过以下选项组合对齐：

• -frounding-math：禁止编译器假设默认舍入模式，严格按当前舍入规则计算
• -ffp-contract=off：禁用融合乘法加法（FMA）优化，避免不同平台因指令支持差异导致结果偏差
• -mfpmath=sse：强制使用SSE指令集，与Windows浮点计算指令集对齐

编译命令示例：

g++ -std=c++11 -frounding-math -ffp-contract=off -mfpmath=sse main.cpp -o main

3. 代码层统一精度处理 #

若无法修改编译环境，可在代码中对结果做显式精度截断，比如保留到指定小数位，或用std::round/std::trunc统一处理输出，消除微小精度差异。

4. 静态链接旧版本glibc（临时方案） #

从Ubuntu 18.04系统中复制glibc 2.27库文件，编译时静态链接该版本库。此方法存在兼容性风险，仅适合临时验证场景。

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