| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 1726 | 回復(fù): 7 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
liming9398金蟲(chóng) (正式寫(xiě)手)
|
[求助]
numerical recipes in fortan 代碼的疑惑 已有1人參與
|
|
|
SUBROUTINE flmoon(n,nph,jd,frac) INTEGER jd,n,nph REAL frac,RAD PARAMETER (RAD=3.14159265/180.) INTEGER i REAL am,as,c,t,t2,xtra c=n+nph/4. t=c/1236.85 t2=t**2 as=359.2242+29.105356*c am=306.0253+385.816918*c+0.010730*t2 jd=2415020+28*n+7*nph xtra=0.75933+1.53058868*c+(1.178e-4-1.55e-7*t)*t2 if(nph.eq.0.or.nph.eq.2)then xtra=xtra+(0.1734-3.93e-4*t)*sin(RAD*as)-0.4068*sin(RAD*am) else if(nph.eq.1.or.nph.eq.3)then xtra=xtra+(0.1721-4.e-4*t)*sin(RAD*as)-0.6280*sin(RAD*am) else pause 'nph is unknown in flmoon' endif if(xtra.ge.0.)then i=int(xtra) else i=int(xtra-1.) endif jd=jd+i frac=xtra-i return END fotran77b標(biāo)準(zhǔn)謝的,為什么要用八位或者九位的常數(shù)給real數(shù)據(jù)RAD,as,am賦值?單精度不是最多七位嗎,多了有意義嗎? |


金蟲(chóng) (正式寫(xiě)手)
|
奧,當(dāng)然有意義: 32位機(jī)的單雙精度的數(shù)據(jù):?jiǎn)尉鹊挠行?shù)字:6位;雙精度的有效數(shù)字:15位 64位機(jī)器上有最大18位有效數(shù)字和4932為指數(shù)的類(lèi)型。 有些計(jì)算需要有效位數(shù)1000多位,比如Pi的計(jì)算等,這時(shí)候就需要自己編寫(xiě)程序來(lái)控制有效位數(shù)的輸出。我記著10年前,我為了在PC上輸出20位有效位數(shù),親自編寫(xiě)了一個(gè)多有效位(名為doubledoubledboule類(lèi)型)的程序,不過(guò)現(xiàn)在網(wǎng)上多的很,樓主可以自己找一找,如果感興趣的話。 |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 372分材料與化工(085600)一志愿湖南大學(xué)求調(diào)劑 +3 | 藍(lán)箋片 2026-04-02 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 348求調(diào)劑 +3 | 吳彥祖24k 2026-04-02 | 3/150 |
|
|
[考研] 292分,材料與化工,申請(qǐng)調(diào)劑 +19 | 程晴之 2026-04-01 | 22/1100 |
|
|
[考研] 311求調(diào)劑 +16 | zchqwer 2026-04-01 | 18/900 |
|
|
[考研]
|
shzhou_ 2026-04-01 | 6/300 |
|
|
[考研] 372求調(diào)劑 +3 | jj涌77 2026-04-02 | 3/150 |
|
|
[考研] 085601一志愿中山大學(xué)深圳材料工程330求調(diào)劑 +7 | pipiver 2026-03-30 | 7/350 |
|
|
[考研] 求調(diào)劑,一志愿南京師范大學(xué)計(jì)算機(jī)專(zhuān)碩,初試373,六級(jí)通過(guò), +3 | 計(jì)算機(jī)追夢(mèng)人 2026-04-01 | 3/150 |
|
|
[考研] 296求調(diào)劑 +4 | 汪?! 2026-03-31 | 7/350 |
|
|
[考研] 一志愿西安交大材料學(xué)碩(英一數(shù)二)347,求調(diào)劑到高分子/材料相關(guān)專(zhuān)業(yè) +7 | zju51 2026-03-31 | 9/450 |
|
|
[考研] 286求調(diào)劑 +5 | lim0922 2026-03-26 | 5/250 |
|
|
[考研] 材料科學(xué)與工程339求調(diào)劑 +11 | hyz0119 2026-03-31 | 12/600 |
|
|
[考研] 285求調(diào)劑 +7 | AZMK 2026-03-30 | 13/650 |
|
|
[考研] 省雙一流重點(diǎn)一本大學(xué)招收調(diào)劑 +4 | wwwwffffff 2026-03-31 | 7/350 |
|
|
[考研] 材料調(diào)劑 +10 | Eujd1 2026-03-31 | 11/550 |
|
|
[考研] 一志愿:西北大學(xué),英一數(shù)一408-284分求調(diào)劑 +7 | 12.27 2026-03-27 | 7/350 |
|
|
[考研] 一志愿華東師范大學(xué)有機(jī)化學(xué)專(zhuān)業(yè),初試351分,復(fù)試被刷求調(diào)劑! +9 | 真名有冰 2026-03-29 | 10/500 |
|
|
[考研] 材料求調(diào)劑 一志愿哈工大總分298分,前三科223分 +11 | dongfang59 2026-03-27 | 11/550 |
|
|
[考研] 福建理工大學(xué)材料學(xué)院先進(jìn)合金團(tuán)隊(duì)招收考研調(diào)劑學(xué)生 +3 | 大華金商都 2026-03-30 | 4/200 |
|
|
[考研] 296求調(diào)劑 +10 | 彼岸t 2026-03-29 | 10/500 |
|