| 5 | 1/1 | 返回列表 |
| 查看: 1722 | 回復(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賦值?單精度不是最多七位嗎,多了有意義嗎? |
|
奧,當(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)上多的很,樓主可以自己找一找,如果感興趣的話。 |


金蟲(chóng) (正式寫(xiě)手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 300求調(diào)劑,材料科學(xué)英一數(shù)二 +8 | leaflight 2026-03-24 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 071000生物學(xué)求調(diào)劑,初試成績(jī)343 +7 | 小小甜面團(tuán) 2026-03-25 | 7/350 |
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +9 | 丹青奶蓋 2026-03-26 | 10/500 |
|
|
[考研] 331環(huán)境科學(xué)與工程求調(diào)劑 +3 | 熠然好運(yùn)氣 2026-03-27 | 3/150 |
|
|
[考研] 352分 化工與材料 +5 | 海納百川Ly 2026-03-27 | 5/250 |
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
[考研] 266求調(diào)劑 +11 | 陽(yáng)陽(yáng)哇塞 2026-03-27 | 12/600 |
|
|
[考研] 085600,材料與化工321分調(diào)劑 +4 | 大饞小子 2026-03-27 | 6/300 |
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +4 | 柳溪邊 2026-03-26 | 4/200 |
|
|
[考研] 352求調(diào)劑 +4 | 大米飯! 2026-03-22 | 4/200 |
|
|
[考研] 一志愿河工大 081700 276求調(diào)劑 +4 | 地球繞著太陽(yáng)轉(zhuǎn) 2026-03-23 | 4/200 |
|
|
[考研] 考研一志愿蘇州大學(xué)初始315(英一)求調(diào)劑 +3 | sbdksD 2026-03-24 | 4/200 |
|
|
[考研] 0854人工智能方向招收調(diào)劑 +4 | 章小魚(yú)567 2026-03-24 | 4/200 |
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
[考研] 材料專(zhuān)碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
[考研] 284求調(diào)劑 +3 | yanzhixue111 2026-03-23 | 6/300 |
|