版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

liming9398

金蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 780.9
散金: 89
紅花: 1
帖子: 361
在線: 1009.2小時(shí)
蟲號(hào): 610070
注冊(cè): 2008-09-23
專業(yè): 代數(shù)學(xué)

[求助] 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賦值？單精度不是最多七位嗎，多了有意義嗎？

回復(fù)此樓

» 猜你喜歡

焦慮已經(jīng)有13人回復(fù)
一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑已經(jīng)有8人回復(fù)
0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有9人回復(fù)
346求調(diào)劑[0856] 已經(jīng)有5人回復(fù)
285化工學(xué)碩求調(diào)劑（081700）已經(jīng)有11人回復(fù)
267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑已經(jīng)有9人回復(fù)
274求調(diào)劑已經(jīng)有6人回復(fù)
求材料調(diào)劑已經(jīng)有10人回復(fù)
本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑已經(jīng)有4人回復(fù)
國(guó)自科面上基金字體已經(jīng)有8人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

求 Risk, Hazards & Crisis in Public Policy期刊縮寫已經(jīng)有6人回復(fù)
應(yīng)蟲友要求，共享《Catalysis in C1 Chemistry》已經(jīng)有292人回復(fù)
【轉(zhuǎn)載】物理著作 The Theory of Transport Phenomena in Solids 已經(jīng)有86人回復(fù)
Nuclear-Power-in-China-Market-Outlook-to-2020 已經(jīng)有19人回復(fù)
Postdoctoral Position in the MSE in NC state university 已經(jīng)有15人回復(fù)
求分子動(dòng)力學(xué)入門fortran程序代碼已經(jīng)有1人回復(fù)
急問(wèn)：連續(xù)資助項(xiàng)目的申請(qǐng)代碼1必須與原青年科學(xué)基金項(xiàng)目保持一致嗎？已經(jīng)有8人回復(fù)
International Journal for Numerical Methods in Engineering二審時(shí)間已經(jīng)有13人回復(fù)
International Journal for Numerical Methods in Engineering 審稿周期已經(jīng)有10人回復(fù)
Numerical Recipes in Fortran 已經(jīng)有181人回復(fù)
【原創(chuàng)+分享】書寫金Au(111)面塊體電極的代碼已經(jīng)有4人回復(fù)

1樓 2012-05-14 11:05:53

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hklulu

木蟲 (正式寫手)

000

應(yīng)助: 18 (小學(xué)生)
金幣: 2645
散金: 61
紅花: 3
帖子: 341
在線: 154.9小時(shí)
蟲號(hào): 108850
注冊(cè): 2005-11-18
專業(yè): 原子和分子物理

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

哈哈，新手吧。

八位不是指的數(shù)字的位數(shù)，而是二進(jìn)制的8B，就是通常用64位二進(jìn)制數(shù)字來(lái)表示雙精度數(shù)字。

此位非彼位。

贊一下

回復(fù)此樓

2樓2012-05-14 11:26:38

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

liming9398

金蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 780.9
散金: 89
紅花: 1
帖子: 361
在線: 1009.2小時(shí)
蟲號(hào): 610070
注冊(cè): 2008-09-23
專業(yè): 代數(shù)學(xué)

可能我沒(méi)說(shuō)清楚3.14159265不就是九位的嗎，可RAD是real類型的也就是單精度類型的，最多七位有效數(shù)字啊，賦值時(shí)寫這么多位還有意義嗎，同樣的問(wèn)題出現(xiàn)在給real類型數(shù)據(jù)as賦值時(shí)表達(dá)式中出現(xiàn)29.105356，也超過(guò)7位有效數(shù)字啊

贊一下

回復(fù)此樓

3樓2012-05-14 11:51:31

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hklulu

木蟲 (正式寫手)

000

應(yīng)助: 18 (小學(xué)生)
金幣: 2645
散金: 61
紅花: 3
帖子: 341
在線: 154.9小時(shí)
蟲號(hào): 108850
注冊(cè): 2005-11-18
專業(yè): 原子和分子物理

【答案】應(yīng)助回帖

奧，當(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ù)此樓

4樓2012-05-14 13:47:26

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Carpathia

銀蟲 (小有名氣)

應(yīng)助: 35 (小學(xué)生)
金幣: 5.6
散金: 12
紅花: 2
帖子: 140
在線: 24.9小時(shí)
蟲號(hào): 1070633
注冊(cè): 2010-08-06
性別: GG
專業(yè): 動(dòng)力學(xué)與控制

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

建議ls看清楚lz的問(wèn)題
lz說(shuō)得沒(méi)錯(cuò)，給單精度賦值超過(guò)位數(shù)是沒(méi)有意義的，同樣給雙精度數(shù)賦值也要加上D0才有意義。Fortran默認(rèn)的是單精度數(shù)。

贊一下

回復(fù)此樓

追求理想的道路上東張西望

5樓2012-05-15 15:02:48

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

liming9398

金蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 780.9
散金: 89
紅花: 1
帖子: 361
在線: 1009.2小時(shí)
蟲號(hào): 610070
注冊(cè): 2008-09-23
專業(yè): 代數(shù)學(xué)

引用回帖:

4樓: Originally posted by hklulu at 2012-05-14 13:47:26:
奧，當(dāng)然有意義：
32位機(jī)的單雙精度的數(shù)據(jù)：?jiǎn)尉鹊挠行?shù)字：6位；雙精度的有效數(shù)字：15位
64位機(jī)器上有最大18位有效數(shù)字和4932為指數(shù)的類型。

有些計(jì)算需要有效位數(shù)1000多位，比如Pi的計(jì)算等，這時(shí)候就需 ...

不太明白你的回復(fù)，你解釋一下嗎？我只是覺(jué)得這是本有名的書，似乎不應(yīng)該犯這種低級(jí)錯(cuò)誤，或許是作者另有初衷，我沒(méi)理解到嗎？

贊一下

回復(fù)此樓

6樓2012-05-15 15:16:59

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jackplay

金蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 1120.8
紅花: 1
帖子: 128
在線: 87.6小時(shí)
蟲號(hào): 2945895
注冊(cè): 2014-01-24
性別: GG
專業(yè): 分析化學(xué)

c語(yǔ)言中，如果最多有效位的數(shù)是這段函數(shù)的有效位，所以較少有效數(shù)位的數(shù)其實(shí)省去了最后的0而已，不知道是不是可以類比

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

maths physics chemistry

7樓2016-01-09 08:21:30

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

43fd6ys

木蟲 (正式寫手)

應(yīng)助: 37 (小學(xué)生)
金幣: 2050.4
紅花: 7
帖子: 379
在線: 147小時(shí)
蟲號(hào): 3935320
注冊(cè): 2015-06-22
專業(yè): 零件成形制造

【答案】應(yīng)助回帖

應(yīng)該沒(méi)錯(cuò)，可能是
當(dāng)表達(dá)式的數(shù)字長(zhǎng)度超過(guò)單精度后，會(huì)采用雙精度進(jìn)行運(yùn)算，最后賦值再轉(zhuǎn)為單精度數(shù)。

贊一下(1人)

回復(fù)此樓

8樓2016-01-09 21:20:09

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 liming9398 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 求材料調(diào)劑 +10	隔壁陳先生 2026-03-12	10/500	2026-03-19 09:15 by 無(wú)際的草原
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +10	秋有木北 2026-03-14	10/500	2026-03-19 05:52 by anny19840123
[考研] 354求調(diào)劑 +4	Tyoumou 2026-03-18	7/350	2026-03-18 21:45 by Tyoumou
[考研] 266求調(diào)劑 +5	陽(yáng)陽(yáng)哇塞 2026-03-14	9/450	2026-03-18 15:05 by stone_128
[考研] 08工科 320總分求調(diào)劑 +5	梨花珞晚風(fēng) 2026-03-17	5/250	2026-03-18 14:49 by haxia
[考研] 311求調(diào)劑 +11	冬十三 2026-03-15	12/600	2026-03-18 14:36 by 星空星月
[考研] 303求調(diào)劑 +4	睿08 2026-03-17	6/300	2026-03-18 11:01 by Iveryant
[考研] 0703化學(xué)336分求調(diào)劑 +6	zbzihdhd 2026-03-15	7/350	2026-03-18 09:53 by zhukairuo
[考研] 268求調(diào)劑 +6	簡(jiǎn)單點(diǎn)0 2026-03-17	6/300	2026-03-18 09:04 by 無(wú)際的草原
[考研] 293求調(diào)劑 +11	zjl的號(hào) 2026-03-16	16/800	2026-03-18 08:10 by zhukairuo
[考研] 本人考085602 化學(xué)工程專碩 +16	不知道叫什么！ 2026-03-15	18/900	2026-03-17 17:05 by ruiyingmiao
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無(wú)懈可擊111
[考研] 211本，11408一志愿中科院277分，曾在中科院自動(dòng)化所實(shí)習(xí) +6	Losir 2026-03-12	7/350	2026-03-17 12:09 by danranxie
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運(yùn)氣yunqi
[考研] 機(jī)械專碩調(diào)劑 +3	笨笨兔子 2026-03-12	3/150	2026-03-15 20:02 by 栗子粥?
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3	麥田耕者 2026-03-14	3/150	2026-03-14 20:04 by 外星文明
[考研] 297一志愿上交085600求調(diào)劑 +5	指尖八千里 2026-03-14	5/250	2026-03-14 17:26 by a不易
[考研] 復(fù)試調(diào)劑 +4	z1z2z3879 2026-03-14	5/250	2026-03-14 16:30 by JourneyLucky
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠

24小時(shí)熱門版塊排行榜

[求助] numerical recipes in fortan 代碼的疑惑 已有1人參與

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

[求助] numerical recipes in fortan 代碼的疑惑已有1人參與

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助: