Re: [問題] 關於fotran數值運算的誤差問題

看板Programming作者jinming.時間18年前 (2006/11/07 02:01)推噓1(1推 0噓 5→)

留言6則, 1人參與討論串2/5 (看更多)

我不懂 fortran 但可能有兩種原因只顯示 0.00000 1.fortran 儲存 real 只用不到 8 bytes(應該不太可能?) 2.write要設定顯示科學記號表示法 ※ 引述《sjgau.bbs@ptt.cc (sjgau)》之銘言： > #include "stdafx.h" > #include <math.h> > int main(int argc, char* argv[]) > { > double a, b, c; > printf("\nsizeof(c)= %d bytes\n", sizeof(c)); > a=0.100000000001e300; > b=0.100000000000e300; > c= a - b; > printf("\nc= %25.22le\n", c); > return 0; > } > /* > sizeof(c)= 8 bytes > c= 1.0000002832399531000000e+288 > Press any key to continue > */ > 手邊暫時沒有 fortran 的 compiler > 先使用對等的 C 幫你看看 > 你看看那個 c 的答案是否為你所要的， > 然後，私下用 e-mail 討論吧 > sjgau4311@gmail.com > 這種簡單的問題討論， > 是會被公幹的 > ※ 引述《jas2902 (..)》之銘言： > : 我想要請教板上的高手關於fortran數值計算上面的誤差問題 > : 參考下面這個簡單的程式 > : -------------------------------- > : program test > : real(kind=8) a, b > : a=0.100000000001*1D300 > : b=0.100000000000*1D300 > : write(*,*) a-b > : end test > : -------------------------------- > : 我發現在螢幕上顯示出來的結果是: > : 0.0000000000000000 > : 但實際的結果應該是: > : 0.000000000001*10^300=10^288 > : 誤差差了 10^288 ! > : 因為類似這樣的計算誤差 > : 讓我在計算量子力學裡面穿透機率的問題 > : 在某些點發生了機率大於1的狀況(機率一定會小於1) > : 想了很久我不知道遇到這種狀況要怎麼樣處理 > : 不知道板上的高手有沒有什麼好辦法可以解決 > : 希望可以幫我個忙謝謝 -- 信言不美，美言不信善者不辯，辯者不善知者不博，博者不知 -- 夫兵者不祥之器物或惡之故有道者不處君子居則貴左用兵則貴右兵者不祥之器非君子之器不得已而用之恬淡為上勝而不美而美之者是樂殺人夫樂殺人者則不可得志於天下矣吉事尚左凶事尚右偏將軍居左上將軍居右言以喪禮處之殺人之眾以哀悲泣之戰勝以喪禮處之道常無名樸雖小天下莫能臣侯王若能守之萬物將自賓天地相合以降甘露民莫之令而自均始制有名名亦既有夫亦將知止知止218-163-120-51.dynamic.hinet.net海