| 5 | 1/1 | 返回列表 |
| 查看: 1725 | 回復(fù): 7 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
liming9398金蟲 (正式寫手)
|
[求助]
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賦值?單精度不是最多七位嗎,多了有意義嗎? |
金蟲 (小有名氣)
|
c語(yǔ)言中,如果最多有效位的數(shù)是這段函數(shù)的有效位,所以較少有效數(shù)位的數(shù)其實(shí)省去了最后的0而已,不知道是不是可以類比 發(fā)自小木蟲Android客戶端 |


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

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工調(diào)劑一志愿大連海事085600,349 +6 | 吃的不少 2026-03-30 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 抱歉 +4 | 田洪有 2026-03-30 | 4/200 |
|
|
[考研] 281求調(diào)劑 +5 | 亞克西good 2026-03-26 | 7/350 |
|
|
[考研] 一志愿鄭州大學(xué),080500學(xué)碩,總分317分求調(diào)劑 +10 | 舉個(gè)栗子oi 2026-03-24 | 11/550 |
|
|
[考研] 0703化學(xué)321分求調(diào)劑 +10 | 三dd. 2026-03-30 | 11/550 |
|
|
[考研] 295材料工程專碩求調(diào)劑 +10 | 1428151015 2026-03-27 | 10/500 |
|
|
[考研] 292求調(diào)劑 +13 | 是妍子也是研子 2026-03-30 | 13/650 |
|
|
[考研] 070300化學(xué)354求調(diào)劑 +7 | 101次希望 2026-03-28 | 7/350 |
|
|
[考研] 求調(diào)劑 +7 | 青春裁為三截 2026-03-29 | 7/350 |
|
|
[考研] 318求調(diào)劑 +7 | 陳晨79 2026-03-30 | 7/350 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-26 | 7/350 |
|
|
[考研] 【求調(diào)劑】085601材料工程專碩 | 總分272 | +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
[考研] 343求調(diào)劑 +5 | 愛(ài)羈絆 2026-03-28 | 5/250 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
[考研] 283求調(diào)劑 +7 | A child 2026-03-28 | 7/350 |
|
|
[考研] 085600,材料與化工321分調(diào)劑 +4 | 大饞小子 2026-03-27 | 6/300 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|