-march native对x86_64 Linux平台浮点计算精度有什么影响？

问题解答 #

1. -march=native对浮点误差的影响 #

两种情况都存在。开启-march=native后，编译器会自动适配当前CPU支持的所有指令集，同时开启对应平台的性能优化，这些优化对浮点精度的影响不是单向的：

• 降低误差的场景：你提到的融合乘加FMA指令就是典型例子，普通a*b + c运算需要先做乘法舍入一次、再加法舍入一次，FMA将两步运算合并，全程只做一次舍入，理论误差更低。你使用的Intel Xeon Gold 6152 CPU的flags中已经标注了fma特性，原生支持该指令。
• 增大误差的场景：一是编译器为了利用AVX/AVX512等SIMD向量指令做并行运算，会重排浮点运算顺序，而浮点加法/乘法不满足严格结合律，运算顺序变化会导致累积误差变化；二是x86平台默认的x87浮点单元使用80位扩展精度做中间运算，而SSE/AVX等SIMD指令使用固定的32位单精度/64位双精度，编译器切换运算单元后中间精度下降，也会导致最终结果误差增大。

2. 精度优化是不是新浮点指令的普遍特性？ #

不是。绝大多数新推出的浮点指令都是在严格符合IEEE 754浮点标准的前提下做性能优化，FMA这类附带精度提升效果的属于少数特例。你当前CPU支持的AVX、AVX2、AVX512系列的基础浮点运算指令，精度和更早的SSE指令完全一致，只是并行度更高、延迟更低。

3. 是否存在效率更高但精度更低的新增指令？ #

存在，你使用的CPU就支持这类指令，典型场景如下：

• 近似运算指令：比如AVX指令集中的VRSQRTPS（单精度倒数平方根近似）、VRCPPS（单精度倒数近似），这类指令的运算延迟只有精确指令的1/3不到，但精度仅为12位左右，远低于单精度浮点数的23位尾数精度，编译器开启-ffast-math类激进优化时会用这类指令替代精确运算，导致误差明显变大。
• 归约运算指令：比如AVX512的向量归约求和指令，会改变原本的串行求和顺序，不同的累加顺序会带来不同的浮点累积误差，大概率会比串行求和的误差更大。

场景适配建议 #

如果需要保持原有测试用例的结果一致性，同时尽可能享受新指令集的性能提升，可以在开启-march=native的同时添加以下编译选项：

• -ffp-contract=off：关闭FMA自动融合，避免乘加运算的精度变化
• -mfpmath=sse：强制所有浮点运算使用SSE/SIMD单元，避免x87和SIMD单元的精度差异
• -fno-approx-func：禁止编译器用近似指令替代精确浮点运算

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