推 ledia:整數->浮點數會有誤差, 浮點數->整數不會有誤差吧? @@ 09/16 11:54
→ ledia:浮點數->整數 指的是 cast (round) 成 int 09/16 11:55
推 VictorTom:小弟我倒是忘了考慮表示法vs整數的特性會不會剛好在這段 09/16 12:09
→ VictorTom:不產生誤差了Orz 回l大, 我說的誤差是浮點數->浮點數後 09/16 12:10
→ VictorTom:面那個浮點數, 因為floor/ceil的回傳值都是浮點數@_@" 09/16 12:10
推 wowtiger:float(4294967295) = ? XD, 另外 x87 有 fisttp 可以用 09/17 00:16
推 phterry:浮點數->整數也有可能發生誤差的情形,請看下例: 09/17 11:32
→ phterry:float x = 1.801; int y = (int)(x * 1000.0); 09/17 11:32
→ phterry:這樣在 x86 上, y的結果會變成1800,而不是1801 09/17 11:33
→ phterry:但是在其他的processer,就能得到正確的1801 09/17 11:34
推 zerodevil:樓上那個是乘法的誤差 不是cast的誤差吧 09/17 11:48
推 phterry:這不單純是乘法的誤差,可以看看下例: 09/17 11:57
→ phterry:int y = (int)(1.801 * 1000.0); 09/17 11:58
→ phterry:以及 float x = 1.801; int y; x *= 1000.0; y = (int)x; 09/17 11:59
→ phterry:以上這兩例都可以在x86上取得正確的1801 09/17 11:59
推 zerodevil:這變因又更多了 剛剛忘了還有10進位轉二進位的誤差 09/17 12:01
→ zerodevil:你要就直接給一個浮點數的binary表示法, 證明它轉成int 09/17 12:01
→ zerodevil:會錯. 要不然就證明這些多餘的動作不會影響結果 09/17 12:02
推 zerodevil:順便一題 你的第一個例子和後兩個的差別在 09/17 12:35
→ zerodevil:第一個是double->int, 後面是double->float->int 09/17 12:35
推 zerodevil: 提^ 09/17 12:38
推 phterry:不知道你有沒有親身驗證過這個問題... 上述的code你可以改 09/17 12:54
→ phterry:成double試看看,結果我就不說明了. 09/17 12:55
推 zerodevil:測過了 你想表達什麼 09/17 12:58
→ phterry:float改成double後,結果有什麼差異否? 09/17 13:00
推 zerodevil:全部是1801 所以呢 09/17 13:01
→ zerodevil:你又舉了三個*完全不同*的case 這有又什麼影響 09/17 13:01
→ phterry:1801是哪一段啊? 09/17 13:05
→ phterry:我實在不想和你戰,覺得沒什麼意義,我只不過想提醒大家可能 09/17 13:07
→ phterry:會被乎略的錯誤而已. x86 的浮點有一些特性會產生和其他 09/17 13:08
→ phterry:CPU不同的結果,如此而已. 09/17 13:10
推 zerodevil:不要越扯越遠 你一開始是說浮點數->整數可能會有誤差 09/17 13:10
→ zerodevil:你是你舉的例子誤差是在double->float的四捨五入 就這樣 09/17 13:11
→ zerodevil:但^ 09/17 13:11
→ phterry:其實我不是很懂我說的這個跟你說的有什麼關係耶... 09/17 13:16
→ phterry:或許你用你的想法解釋為何會有1800的結果會比較理想 09/17 13:17
推 zerodevil:1000.0是double, 1000.0f才是float 09/17 13:18
→ zerodevil:所以你所有運算會被先cast成double算 09/17 13:18
→ zerodevil:(int)(x*1000.0)是double直接變成int 09/17 13:19
→ zerodevil:x*=1000.0; (int)x是做完double乘法後轉float再轉int 09/17 13:19
推 phterry:double x; int y = (int)(x * 1000.0); 請問你的結果? 09/17 13:19
推 zerodevil:1800啊= = 有什麼就直說好嗎 這樣很煩 09/17 13:21
→ zerodevil:啊看錯 1801.. 09/17 13:21
推 phterry:不好意思,我剛x忘了給初值,應該要double x = 1.801; 09/17 13:22
→ phterry:我用x86 32bit Linux執行的結果是1800 09/17 13:23
→ phterry:在ARM Linux上跑,結果卻是1801,而1801也是我們所希望的 09/17 13:23
推 zerodevil:x*1000.0在x是float跟double時算出來會有一點誤差 09/17 13:24
→ zerodevil:在cast成int前就有差了 09/17 13:24
→ phterry:再來 float x = 1.801; int y = (int)(x * 1000.0f); 09/17 13:26
→ phterry:上面這樣寫,x86 -> 1800, ARM->1801 09/17 13:26
推 zerodevil:乘法的實作方式不同吧.. 09/17 13:31
→ phterry:你的意思是compiler還是CPU? 09/17 13:32
推 zerodevil:cpu 09/17 13:32
推 zerodevil:可以用 cout << (unsigned int&)(x * 1000.0f); 09/17 13:33
→ zerodevil:看有沒有不一樣 09/17 13:33
→ phterry:那原先的double->float四捨五入之說應該.... 09/17 13:34
推 zerodevil:啊? 那是同機器上不同code的比較啊... 09/17 13:35
→ phterry:有一篇文章是一個日本人寫的,可能有點久了,我暫時找不到 09/17 13:36
→ phterry:大意是說x86的浮點有一個擴張精度,會讓指數部份變成15bit 09/17 13:37
→ phterry:而不是一般所認知的11bit 09/17 13:37