建站
咨詢
售后
建站咨詢
新聞中心
昨天阿楠發(fā)現(xiàn)了項目中的一個 bug ，是因為浮點運算的前后不一致導致的。明明是完全相同的 C 代碼，參數(shù)也嚴格一致，但是計算出了不相同的結果。我對這個現(xiàn)象非常感興趣，仔細研究了一下成因。
專注于為中小企業(yè)提供做網(wǎng)站、網(wǎng)站設計服務,電腦端+手機端+微信端的三站合一,更高效的管理,為中小企業(yè)清江浦免費做網(wǎng)站提供優(yōu)質(zhì)的服務。我們立足成都，凝聚了一批互聯(lián)網(wǎng)行業(yè)人才，有力地推動了上千家企業(yè)的穩(wěn)健成長，幫助中小企業(yè)通過網(wǎng)站建設實現(xiàn)規(guī)模擴充和轉(zhuǎ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 4.9.2 ，強制使用 x87 浮點運算編譯運行，你會發(fā)現(xiàn)令人詫異的結果。
- gcc a.c -mfpmath=387
- 19
- 20
前一次的輸出是 19 ，后一次是 20 。
這是為什么呢?讓我們來看看 gcc 生成的代碼，我截取了相關的段落：
- flds 16(%rbp)
- flds .LC0(%rip)
- fmulp %st, %st(1)
- fstps -4(%rbp) ; 1. x * 0.01f 結果保存到內(nèi)存中的 float 變量中
- flds 16(%rbp)
- flds .LC0(%rip)
- fmulp %st, %st(1)
- fisttpl -20(%rbp) ; 2. x * 0.01f 結果直接轉(zhuǎn)換為整型
- 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
這里我做了三行注釋。
首先，0.01 是無法精確表示成 2 進制的，所以 * 0.01 這個操作一定會存在誤差。
兩次運算都是 x * 0.01f ，雖然按 C 語言的轉(zhuǎn)換規(guī)則，表達式中都是 float 時，按 float 精度運算。但這里 gcc 生成的代碼并沒有嚴格設置 FPU 的精度控制，在注釋 2 這個地方，乘法結果是直接從浮點寄存器轉(zhuǎn)換為整數(shù)的。而在注釋 1 這個地方，把乘法結果通過 fstps 以低精度形式保存到內(nèi)存，再在注釋 3 的地方 flds 讀回。
所以在注釋 2 和注釋 3 的地方，浮點寄存器 st 內(nèi)的值其實是有差別的，這導致了 fisttpl 轉(zhuǎn)換為整數(shù)后結果不同。
文章名稱：浮點運算潛在的結果不一致問題
URL標題：http://www.fisionsoft.com.cn/article/dpohjhc.html
