【Python 線上課程試看搶先體驗 Part1】
📍課程募資倒數1週，55折敬請把握▶︎ https://pse.is/BSPZA
📍真正適合新手的Python課程▶︎計算機概論章節 #已解鎖送給學員！
📍達到350人還會解鎖程式人都該學的『git & github』章節喔！
#影片主題：為什麼在程式裡， 0.1 + 0.2 不等於 0.3？
#重點筆記：
✏️程式語言中無法表示「無限位數」的小數，因此改稱為「浮點數」
✏️浮點數本身其實具有誤差，為什麼呢？
✏️有經驗的工程師，在進行計算時又是如何驗證與除錯的？
這些問題的答案，看影片都有解答喔！
我總是跟同學說：
沒有一套程式教學適合所有人學習，但你可以透過我的教學影片體驗我的教學風格。再決定你要不要跟著紀老師學Python，我能為入門同學做到的是：
1. 注意專有名詞的解說
2. 運用生活例子說明難懂的程式邏輯
3. 課程中你的程式碼問題一定設法幫您解答到您滿意
希望同學能透過我的線上課程效率學習、收穫滿滿
⋯⋯⋯⋯
其他免費教學影片傳送門，這邊走▶︎ https://pse.is/BSPZA
55 折早鳥優惠只到11/28，2019 新技能趁現在搶先投資！
浮點數計算機 在 [問題] IEEE 754 浮點數運算以及rounding請益- 看板C_and_CPP 的推薦與評價
大家好
小弟目前在修系上計算機方面的課程
老師希望我們寫出一個浮點數的運算單元,
然後以IEEE 754 單精度浮點數為標準,rounding方式是round to nearest even
而老師給的測資都是標準二進位浮點數的格式了,不用自己再轉為IEEE 754浮點數
即輸入會像是這樣:
0 10010011 0000 0000 0000 0000 1111 111
| |------| |--------------------------|
sign exponent mantissa
而我在寫之前想要先確定一下自己觀念對不對
不然到時候觀念不對怕浪費太多時間改程式
所以想要問一個運算的問題:
(以IEEE 754 單精度浮點數為例
,即1個sign bit,8個 exponent bit,23個mantissabit)
兩個浮點數在算加減法的時候，exponent小的mantissa要對齊exponent大的mantissa
也就是要看兩個浮點數的exponent差距多少來看mantissa要移位多少
那如果exponent小的那個的mantissa在移位過後超過mantissa所能表示的範圍
要把超過範圍的那幾個bit一起算，還是要捨去呢?
舉例來說
我要算兩個浮點數相減
第一個數:
0 10010011 0000 0000 0000 0000 1111 111
| |------| |--------------------------|
sign exponent mantissa
第二個數:
1 10001110 0000 0000 0000 0111 1111 111
| |------| |--------------------------|
sign exponent mantissa
第一個數的exponent換成十進位是147，第二個數的exponent換成十進位是142
而147-127(bias)=20,142-127=15
所以事實上上面兩個數可以變為:
第一個數:
1.0000 0000 0000 0000 1111 111 * 2^20
第二個數:
-1.0000 0000 0000 0111 1111 111 * 2^15
因為第二個數比第一個數的次方少五,所以要右移5個bit
那麼問題來了，移完之後是會變成
(一)所有bit都保留,因此共要28bit表示mantissa
-0.0000 1000 0000 0000 0011 1111 1111 *(2^20)
|----|
這五個bit超過23bit
(二)超過23bit之後直接砍掉,因此滿足23bit表示mantissa
-0.0000 1000 0000 0000 0011 111 *(2^20)
(三)加入round,guard,sticky三個bit去考慮,因此用25bit表示mantissa
-0.0000 1000 0000 0000 0011 1111 1 且設S=1(因為砍掉後面三個1)
| |
G R
是上面(一)、(二)、(三)的哪一種呢?
因為這三種不同的移位方式會造成最後運算答案都不一樣，
所以我想IEEE 754應該會有明確的規範。
我個人是比較傾向於第(三)種，
因為如果是第(一)種的話，兩個浮點數若exponent差太多
那就要保存一大堆數字，像是兩數的exponent如果差了一百
那小的exponent很可能就要保存一百個0外加原本的23個mantissa
等於要保存123個bit，以硬體的角度而言應該是不會這樣設計?
還請各位替我解答一下
老師上課浮點數就只是帶過 然後就要我們一周寫出來
實在是有點頭痛Orz
感謝各位!
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 1.160.137.74 (臺灣)
※ 文章網址: https://www.ptt.cc/bbs/C_and_CPP/M.1590905575.A.34F.html
那可以請問一下如果要符合要求,小exponent的那個mantissa到底該怎麼移嗎?
因為要對齊大exponent的話,一定會有幾個bit超過mantissa能保存的範圍(23bit)
我有去翻Computer Organization and Design: The Hardware/Software Interface
這本書,但是他舉的例子好像比較接近我提的第(三)種QQ
那請問我的G,R,S有取錯嗎?
(G=第24bit,R=第25bit,S=第25bit後所有bit做OR,所以G=1,R=1,S=1)
那您的意思是:
加入round,guard,sticky三個bit去考慮,
-0.0000 1000 0000 0000 0011 1111 1 且設S=1(因為砍掉後面三個1)
| |
G R
而因為G,R = 1
=>進位
=> -0.0000 1000 0000 0000 0011 111
+) 1
----------------------------------------------------
-0.0000 1000 0000 0000 0100 000
=>第二個數在移完之後變成:-0.0000 1000 0000 0000 0100 000 * 2^20
=>和第一個數運算
是這樣子嗎?
在麻煩您解答一下了,感激不盡
f大不好意思,
如果照您說的"rounding是運算完再做"的話,
那我所說的第二個浮點數:
-1.0000 0000 0000 0111 1111 111 * 2^15
要右移五次方的話,超過mantissa 23bit的那五個bit
到底該怎麼處理,去和第一個浮點數:
1.0000 0000 0000 0000 1111 111 * 2^20
運算完後才有正確的結果呢?
我現在用: https://weitz.de/ieee/ 這個網站去算我這兩個數字
答案會是 1.1111 0000 0000 0001 0111 110 * 2^19
但是我怎麼算都算不出這個結果
不好意思問題很多 但是我真的被這問題困擾很久了QQ
我能理解您所說的"做完運算(減完後)"還要再rounding,
但是我的問題比較像是,
在中間過程中,如果有一個數字在移位後,
他的mantissa超過了原本格式所能表示的範圍,
那這些超過的bit要如何處理它們,
才會使我可以正確的運算並且得到正確的結果
然後最後再做rounding呢?
因為如果這幾個bit不正確處理的話,最後運算結果再rounding應該也不會對
所以我目前是卡在中間過程中不知道該如何處理這些超過mantissa的bits QQ
您說的"超過的部分在進位前做保留"，意思是不管超過幾位，
運算的時候就是把這些超過的存起來照算嗎？
如果是這樣的話我再算一次的過程是：
第一個數:
0 10010011 0000 0000 0000 0000 1111 111
| |------| |--------------------------|
sign exponent mantissa
第二個數:
1 10001110 0000 0000 0000 0111 1111 111
| |------| |--------------------------|
sign exponent mantissa
化成指數形式:
第一個數:
1.0000 0000 0000 0000 1111 111 * 2^20
第二個數:
-1.0000 0000 0000 0111 1111 111 * 2^15
第二個數右移5BIT
=>-0.0000 1000 0000 0000 0011 1111 1111 *(2^20)
|----|
這五個bit超過23bit,但是保留
然後第一個數的1.mantissa減第二個數的1.mantissa:
1.0000 0000 0000 0000 1111 1110 0000
- 0.0000 1000 0000 0000 0011 1111 1111
---------------------------------------------
0.1111 1000 0000 0000 1011 1110 0001
因此目前算出的答案是:0.1111 1000 0000 0000 1011 1110 0001 * 2^20
接著normalized:1.111 1000 0000 0000 1011 1110 0001 * 2^19
但是1.111 1000 0000 0000 1011 1110 0001 * 2^19
||
GR
而有4bit超過23bit,因此要rounding,而G=0,R=0>捨去
其中我G取第24bit,R取第25bit,S=1
因此最後答案為:1.111 1000 0000 0000 1011 1110 * 2^19
這樣子我的運算觀念是對的嗎?
答案和 https://weitz.de/ieee/ 算出來的一樣
但是我怕我會不會是觀念錯誤只是剛好矇對(像是GR取錯之類的)
再請f大幫我看一下了,謝謝您
※ 編輯: ayn775437403 (1.160.137.74 臺灣), 05/31/2020 16:13:13
太好了!真的非常感謝您不厭其煩的替我解答這麼久!
希望如果有問題的話還可以再請教您一下!謝謝!
※ 編輯: ayn775437403 (1.160.137.74 臺灣), 05/31/2020 16:14:38
有，我有去翻IEEE Standard for Floating-Point Arithmetic這個規格書
但是似乎都是在講"運算過後的結果如何rounding"
好像沒有提到如果"運算過程中移位後超過mantissa範圍的bit該如何處理"
當然可能是我沒看很仔細 QQ
※ 編輯: ayn775437403 (123.195.195.29 臺灣), 05/31/2020 22:54:52
了解,謝謝您的回復
所以不管中間運算過程中的浮點數超過了mantissa所能表示範圍的多少bit
仍然要把它一起加入運算就對了
然後最後算完的結果如果有多出mantissa所能表示的範圍再rounding掉
應該是這樣理解沒錯吧@@?
謝謝
※ 編輯: ayn775437403 (123.195.195.29 臺灣), 06/01/2020 15:43:18
... <看更多>