liming9398

金蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 780.9
散金: 89
紅花: 1
帖子: 361
在線: 1009.2小時
蟲號: 610070
注冊: 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標準謝的，為什么要用八位或者九位的常數(shù)給real數(shù)據(jù)RAD，as,am賦值？單精度不是最多七位嗎，多了有意義嗎？

回復(fù)此樓

» 猜你喜歡

材料專碩英一數(shù)二306 已經(jīng)有5人回復(fù)
085600材料與化工調(diào)劑 324分已經(jīng)有9人回復(fù)
0703化學(xué)調(diào)劑已經(jīng)有12人回復(fù)
0703化學(xué) 305求調(diào)劑已經(jīng)有4人回復(fù)
0703化學(xué)調(diào)劑，求各位老師收留已經(jīng)有10人回復(fù)
271材料工程求調(diào)劑已經(jīng)有5人回復(fù)
281求調(diào)劑（0805）已經(jīng)有16人回復(fù)
304求調(diào)劑已經(jīng)有6人回復(fù)
材料工程專碩調(diào)劑已經(jīng)有6人回復(fù)
一志愿天大材料與化工（085600）總分338 已經(jīng)有4人回復(fù)

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

求 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ù)
求分子動力學(xué)入門fortran程序代碼已經(jīng)有1人回復(fù)
急問：連續(xù)資助項目的申請代碼1必須與原青年科學(xué)基金項目保持一致嗎？已經(jīng)有8人回復(fù)
International Journal for Numerical Methods in Engineering二審時間已經(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小時
蟲號: 108850
注冊: 2005-11-18
專業(yè): 原子和分子物理

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

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

哈哈，新手吧。

八位不是指的數(shù)字的位數(shù)，而是二進制的8B，就是通常用64位二進制數(shù)字來表示雙精度數(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小時
蟲號: 610070
注冊: 2008-09-23
專業(yè): 代數(shù)學(xué)

可能我沒說清楚3.14159265不就是九位的嗎，可RAD是real類型的也就是單精度類型的，最多七位有效數(shù)字啊，賦值時寫這么多位還有意義嗎，同樣的問題出現(xiàn)在給real類型數(shù)據(jù)as賦值時表達式中出現(xiàn)29.105356，也超過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小時
蟲號: 108850
注冊: 2005-11-18
專業(yè): 原子和分子物理

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

奧，當然有意義：
32位機的單雙精度的數(shù)據(jù)：單精度的有效數(shù)字：6位；雙精度的有效數(shù)字：15位
64位機器上有最大18位有效數(shù)字和4932為指數(shù)的類型。

有些計算需要有效位數(shù)1000多位，比如Pi的計算等，這時候就需要自己編寫程序來控制有效位數(shù)的輸出。我記著10年前，我為了在PC上輸出20位有效位數(shù)，親自編寫了一個多有效位（名為doubledoubledboule類型）的程序，不過現(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小時
蟲號: 1070633
注冊: 2010-08-06
性別: GG
專業(yè): 動力學(xué)與控制

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

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

建議ls看清楚lz的問題
lz說得沒錯，給單精度賦值超過位數(shù)是沒有意義的，同樣給雙精度數(shù)賦值也要加上D0才有意義。Fortran默認的是單精度數(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小時
蟲號: 610070
注冊: 2008-09-23
專業(yè): 代數(shù)學(xué)

引用回帖:

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

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

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

贊一下

回復(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小時
蟲號: 2945895
注冊: 2014-01-24
性別: GG
專業(yè): 分析化學(xué)

c語言中，如果最多有效位的數(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小時
蟲號: 3935320
注冊: 2015-06-22
專業(yè): 零件成形制造

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

應(yīng)該沒錯，可能是
當表達式的數(shù)字長度超過單精度后，會采用雙精度進行運算，最后賦值再轉(zhuǎn)為單精度數(shù)。

贊一下(1人)

回復(fù)此樓

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

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

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

[考研] 085600材料與化工調(diào)劑 324分 +9	llllkkkhh 2026-03-18	9/450	2026-03-19 07:40 by BruceLiu320
[考研] 311求調(diào)劑 +4	冬十三 2026-03-18	4/200	2026-03-18 21:47 by 盡舜堯1
[考研] 321求調(diào)劑 +3	何潤采123 2026-03-18	3/150	2026-03-18 21:27 by li123456789.
[考研] 0703化學(xué)調(diào)劑，六級已過，有科研經(jīng)歷 +10	曦熙兮 2026-03-15	10/500	2026-03-18 14:19 by 007_lilei
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 0854，計算機類招收調(diào)劑 +3	胡辣湯放糖 2026-03-15	6/300	2026-03-18 12:09 by 上岸上岸……..
[考研] 303求調(diào)劑 +4	睿08 2026-03-17	6/300	2026-03-18 11:01 by Iveryant
[考研] 265求調(diào)劑 +3	梁梁校校 2026-03-17	3/150	2026-03-18 09:12 by zhukairuo
[考研] 275求調(diào)劑 +4	太陽花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 0703一志愿211 285分求調(diào)劑 +5	ly3471z 2026-03-13	5/250	2026-03-16 16:16 by 哦哦123
[考研] 0856求調(diào)劑 +3	劉夢微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 22408總分284求調(diào)劑 +3	InAspic 2026-03-13	3/150	2026-03-15 11:10 by zhq0425
[考研] 復(fù)試調(diào)劑 +4	z1z2z3879 2026-03-14	5/250	2026-03-14 16:30 by JourneyLucky
[考研] 學(xué)碩285求調(diào)劑 +13	Wisjxn 2026-03-12	46/2300	2026-03-14 10:33 by JourneyLucky
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 304求調(diào)劑 +7	7712b 2026-03-13	7/350	2026-03-13 21:42 by peike
[考研] 工科278分求調(diào)劑 +5	周慢熱啊 2026-03-12	7/350	2026-03-13 15:49 by JourneyLucky
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生，2026.03.20截止 +3	Xiangyu_ou 2026-03-12	3/150	2026-03-13 09:36 by duanwu655