| 查看: 4659 | 回復(fù): 14 | |||
hshx9291銀蟲 (初入文壇)
|
[求助]
Fortran計算10的十次方以上運算結(jié)果不正確,是精度問題嗎? 已有1人參與
|
| 最近用Fortran進行模擬計算,計算結(jié)果顯示錯誤,往回找某一參數(shù)的計算結(jié)果發(fā)現(xiàn)為0.0000000E+000,后來檢查式子發(fā)現(xiàn)有個10的12次方計算。然后我就另建新項目計算了10的幾次方計算,發(fā)現(xiàn)10的9次方計算結(jié)果是1000000000正常,但是10的十次方計算結(jié)果就成了1410065408,11次方為1215752192,12次方就成了-727379968。一開始懷疑是精度問題,但是設(shè)置成精度為real(8),real(16)結(jié)果計算都是這樣。求問大神是為什么? |
至尊木蟲 (正式寫手)
巾幗將軍
|
是因為整數(shù),fortran對整數(shù)的存儲是有極限的,你把做計算的變量定義成雙精度實型,最后它會以科學(xué)計數(shù)法的形式給出結(jié)果。 發(fā)自小木蟲IOS客戶端 |

木蟲 (著名寫手)
|
你用的四字節(jié)整數(shù),就只能到這個數(shù)量級。如果需要更多,可以告訴編譯器使用八字節(jié)整數(shù)(ifort 的flag 是-i8),可以到二十位數(shù)。再多一般就沒有使用精確數(shù)值的必要了,可以用浮點數(shù)近似代替 發(fā)自小木蟲IOS客戶端 |
銀蟲 (初入文壇)
銀蟲 (初入文壇)
木蟲 (著名寫手)
|
你這個變量是整數(shù),不是實型。八字節(jié)整數(shù)定義為integer*8 如果不想改代碼可以用flag讓編譯器自動把所有整形視為長整數(shù),windows沒用過, 但猜測應(yīng)該是/i8 發(fā)自小木蟲IOS客戶端 |
至尊木蟲 (正式寫手)
巾幗將軍

銀蟲 (初入文壇)
銀蟲 (初入文壇)
木蟲 (著名寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 297求調(diào)劑 +11 | 戲精丹丹丹 2026-03-17 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-21 | 3/150 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂29 2026-03-18 | 5/250 |
|
|
[考研]
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +9 | 手機用戶 2026-03-17 | 9/450 |
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
[考研] 0854可跨調(diào)劑,一作一項核心論文五項專利,省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨; 2026-03-16 | 3/150 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|