[[201284]]

昨天阿楠發(fā)現(xiàn)了項目中的一個 bug ，是因為浮點運算的前后不一致導致的。明明是完全相同的 C 代碼，參數(shù)也嚴格一致，但是計算出了不相同的結果。我對這個現(xiàn)象非常感興趣，仔細研究了一下成因。

原始代碼比較繁雜。在弄清楚原理后，我簡化了出問題的代碼，重現(xiàn)了這個問題：

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