| 10 | 1/1 | 返回列表 |
| 查看: 4137 | 回復(fù): 9 | |||
fish8900金蟲(chóng) (小有名氣)
|
[求助]
fortran中運(yùn)行中出現(xiàn)無(wú)窮大的數(shù)怎么處理
|
|
在計(jì)算的過(guò)程中出現(xiàn)了無(wú)窮大的數(shù)(然后程序終止),我想進(jìn)行以下幾種嘗試,不知如何進(jìn)行 1、當(dāng)A變得無(wú)窮大后,我想對(duì)它進(jìn)行賦值,令A(yù)=1.0E+10 2、當(dāng)A變得無(wú)窮大后,想用一個(gè)判據(jù)來(lái)進(jìn)行后來(lái)的其他計(jì)算,這個(gè)判據(jù)該怎么設(shè)置 3、如何判斷A變得無(wú)窮大了 |
北大科研分享 |
木蟲(chóng) (小有名氣)
金蟲(chóng) (小有名氣)
|
您好,這是根據(jù)要求寫(xiě)的一點(diǎn)代碼,希望對(duì)您有幫助。 LOGICAL ISLARGE !邏輯變量,是否是大數(shù) REAL LARGENUM, A PARAMETER(LARGENUM=1.0E10) !也可按需求設(shè)為其他大數(shù) ISLARGE=.FALSE. !判斷A是否是大數(shù),只考慮了正數(shù)情況, IF(A.GE.LARGENUM) !如果還要考慮負(fù)數(shù)情況則改為IF(ABS(A).GE.LARGENUM) A=1.0E10 !是大數(shù)則設(shè)定值為1.0E10 ISLARGE=.TRUE. ENDIF IF(ISLARGE) PRINT 'A is very large' !輸出信息,判斷A為大數(shù)之后的后續(xù)操作 ENDIF |
金蟲(chóng) (小有名氣)
金蟲(chóng) (小有名氣)
金蟲(chóng) (小有名氣)
送紅花一朵 |
1、您說(shuō)的誤差隨步長(zhǎng)增加而增長(zhǎng)了是什么意思呢?我在寫(xiě)程序過(guò)程中是規(guī)定了當(dāng)前后的計(jì)算值(T-T1)/T<=1.0e-03時(shí),認(rèn)定計(jì)算結(jié)果達(dá)到收斂。而且我在調(diào)試過(guò)程中是將某些數(shù)值調(diào)大一點(diǎn)后才出現(xiàn)無(wú)窮大的情況。 2、當(dāng)A的相關(guān)變量變得不合理時(shí),我試著將它變成一個(gè)可以接受的值,但是發(fā)現(xiàn)還是會(huì)出現(xiàn)計(jì)算停止的現(xiàn)象(結(jié)果顯示也還是出現(xiàn)sqrt的問(wèn)題),想問(wèn)一下您猜測(cè)原因可能是什么? (可能有點(diǎn)憑空說(shuō)了,不好意思哈,沒(méi)具體程序給你看) |
金蟲(chóng) (小有名氣)
|
之前我猜測(cè)可能的原因之一是算法數(shù)值不穩(wěn)定導(dǎo)致誤差隨步長(zhǎng)增加而增長(zhǎng)(有類(lèi)似大數(shù)吃小數(shù),相近數(shù)相減等情況導(dǎo)致數(shù)值誤差),但不一定就是你的情況。 既然不能貼出具體代碼,就得靠你自己調(diào)試找錯(cuò)了。 調(diào)試的時(shí)候可以看看相關(guān)變量的值在調(diào)函數(shù)和子程序之前和之后的變化情況是否合理,還可以試試用Excel按程序里的步驟代入算一遍,看是否得到相同的結(jié)果。 |
金蟲(chóng) (小有名氣)
新蟲(chóng) (初入文壇)
銀蟲(chóng) (正式寫(xiě)手)
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 303求調(diào)劑 +5 | DLkz1314. 2026-03-30 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué)求調(diào)劑 +6 | 丹青奶蓋 2026-03-26 | 8/400 |
|
|
[考研] 293求調(diào)劑 +3 | 末未mm 2026-03-30 | 5/250 |
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
[考研] 求調(diào)劑 +7 | 研研,接電話 2026-03-24 | 8/400 |
|
|
[考研] 求調(diào)劑 +7 | 青春裁為三截 2026-03-29 | 7/350 |
|
|
[考研] 085701求調(diào)劑初試286分 +5 | secret0328 2026-03-28 | 5/250 |
|
|
[考研] 085602 化學(xué)工程專(zhuān)碩 340分求調(diào)劑 +4 | qianbai11 2026-03-29 | 4/200 |
|
|
[考研] 一志愿雙一流機(jī)械285分求調(diào)劑 +4 | 幸運(yùn)的三木 2026-03-29 | 5/250 |
|
|
[考研] 343求調(diào)劑 +6 | 愛(ài)羈絆 2026-03-29 | 6/300 |
|
|
[考研] 332求92調(diào)劑 +8 | 蕉蕉123 2026-03-28 | 8/400 |
|
|
[考研] 學(xué)碩274求調(diào)劑 +9 | Li李魚(yú) 2026-03-26 | 9/450 |
|
|
[考研] 一志愿北化085600材料專(zhuān)碩275|有文章專(zhuān)利|求調(diào)劑 +7 | Micky11223 2026-03-25 | 7/350 |
|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
[材料工程] 一志愿C9材料與化工專(zhuān)業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 0703化學(xué)338求調(diào)劑! +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
[考研] 309求調(diào)劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|