浮点运算潜在的结果不一致问题

[[201284]]

昨天阿楠发现了项目中的一个 bug ，是因为浮点运算的前后不一致导致的。明明是完全相同的 C 代码，参数也严格一致，但是计算出了不相同的结果。我对这个现象非常感兴趣，仔细研究了一下成因。

原始代码比较繁杂。在弄清楚原理后，我简化了出问题的代码，重现了这个问题：

static void 
foo(float x) { 
    float xxx = x * 0.01f; 
    printf("%d\n", (int)(x * 0.01f)); 
    printf("%d\n", (int)xx); 
} 
 
int 
main() { 
    foo(2000.0f); 
    return 0; 
} 

gcc a.c -mfpmath=387 
 
19 
20 

这是为什么呢?让我们来看看 gcc 生成的代码，我截取了相关的段落：

flds    16(%rbp) 
 flds    .LC0(%rip) 
 fmulp   %st, %st(1) 
 fstps   -4(%rbp)          ; 1. x * 0.01f 结果保存到内存中的 float 变量中 
 flds    16(%rbp) 
 flds    .LC0(%rip) 
 fmulp   %st, %st(1) 
 fisttpl -20(%rbp)        ; 2. x * 0.01f 结果直接转换为整型 
 movl    -20(%rbp), %eax 
 movl    %eax, %edx 
 leaq    .LC1(%rip), %rcx 
 call    printf 
 flds    -4(%rbp)                 ; 3. 读出 1. 保存的乘法结果 
 fisttpl -20(%rbp) 
 movl    -20(%rbp), %eax 
 movl    %eax, %edx 
 leaq    .LC1(%rip), %rcx 
 call    printf