循环展开后的汇编差异导致浮点计算结果不一致

问题根源 #

核心原因是浮点数加法不满足结合律：浮点数是有限精度表示，不同的累加顺序会产生不同的舍入误差。两个逻辑一致的循环经GCC 4.9.2展开后，编译器选择了不同的累加顺序，导致最终结果出现细微差异；而未展开时，两者均按i=0到7的顺序依次累加，因此结果完全匹配。

解决办法 #

1. 禁止测试循环的展开优化 #

让测试循环保持和原循环未展开时的执行逻辑一致，通过编译器属性实现：

__attribute__((optimize("no-unroll-loops")))
float TestFunc(OuterStruct *b)
{
    float a = 0.0f;
    for (int i = 0; i < 8; i++)
        a += b->c[i]->d;
    return a;
}

也可以在编译测试代码时添加-fno-unroll-loops选项（针对整个编译单元）。这种方式无需修改循环逻辑，就能保证结果和原循环未展开时完全一致。

2. 手动按固定顺序展开循环 #

观察原循环展开后的汇编，它严格按i=0到7的顺序逐个累加。手动将测试循环展开成相同顺序，强制编译器生成一致的执行逻辑：

float TestFunc(OuterStruct *b)
{
    float a = 0.0f;
    a += b->c[0]->d;
    a += b->c[1]->d;
    a += b->c[2]->d;
    a += b->c[3]->d;
    a += b->c[4]->d;
    a += b->c[5]->d;
    a += b->c[6]->d;
    a += b->c[7]->d;
    return a;
}

这种方式通过固定累加顺序，确保编译器生成的汇编和原循环展开后的逻辑对齐，彻底消除结果差异。

3. 用编译选项强制浮点运算一致性 #

若需要保留循环展开的性能优势，可尝试添加-ffloat-store编译选项。该选项会强制将浮点运算结果写回内存，避免因寄存器精度（部分架构的浮点寄存器精度高于float）导致的额外误差，同时可能让编译器的展开逻辑更接近原循环。不过这个选项会带来轻微的性能损耗。

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