应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » numerical recipes in fortan 代码的疑惑
北京石油化工学院2026年研究生招生接收调剂公告
81/1返回列表
查看: 1744  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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标准谢的，为什么要用八位或者九位的常数给real数据RAD，as,am赋值？单精度不是最多七位吗，多了有意义吗？
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-05-14 11:05:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hklulu

木虫 (正式写手)

000

【答案】应助回帖

感谢参与，应助指数 +1
哈哈，新手吧。

八位不是指的数字的位数，而是二进制的8B，就是通常用64位二进制数字来表示双精度数字。

此位非彼位。
赞一下 回复此楼
00
2楼2012-05-14 11:26:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liming9398

金虫 (正式写手)
可能我没说清楚3.14159265不就是九位的吗，可RAD是real类型的也就是单精度类型的，最多七位有效数字啊，赋值时写这么多位还有意义吗，同样的问题出现在给real类型数据as赋值时表达式中出现29.105356，也超过7位有效数字啊
赞一下 回复此楼
3楼2012-05-14 11:51:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hklulu

木虫 (正式写手)

000

【答案】应助回帖

奥，当然有意义：
32位机的单双精度的数据：单精度的有效数字：6位；双精度的有效数字：15位
64位机器上有最大18位有效数字和4932为指数的类型。

有些计算需要有效位数1000多位，比如Pi的计算等，这时候就需要自己编写程序来控制有效位数的输出。我记着10年前，我为了在PC上输出20位有效位数，亲自编写了一个多有效位（名为doubledoubledboule类型）的程序，不过现在网上多的很，楼主可以自己找一找，如果感兴趣的话。
赞一下 回复此楼
00
4楼2012-05-14 13:47:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Carpathia

银虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
建议ls看清楚lz的问题
lz说得没错，给单精度赋值超过位数是没有意义的，同样给双精度数赋值也要加上D0才有意义。Fortran默认的是单精度数。
赞一下 回复此楼
追求理想的道路上东张西望
5楼2012-05-15 15:02:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liming9398

金虫 (正式写手)
引用回帖:
4楼: Originally posted by hklulu at 2012-05-14 13:47:26:
奥，当然有意义：
32位机的单双精度的数据：单精度的有效数字：6位；双精度的有效数字：15位
64位机器上有最大18位有效数字和4932为指数的类型。

有些计算需要有效位数1000多位，比如Pi的计算等，这时候就需 ...

不太明白你的回复，你解释一下吗？我只是觉得这是本有名的书，似乎不应该犯这种低级错误，或许是作者另有初衷，我没理解到吗？
赞一下 回复此楼
6楼2012-05-15 15:16:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jackplay

金虫 (小有名气)
c语言中，如果最多有效位的数是这段函数的有效位，所以较少有效数位的数其实省去了最后的0而已，不知道是不是可以类比

发自小木虫Android客户端
赞一下 回复此楼
maths physics chemistry
7楼2016-01-09 08:21:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

43fd6ys

木虫 (正式写手)

【答案】应助回帖

应该没错，可能是
当表达式的数字长度超过单精度后，会采用双精度进行运算，最后赋值再转为单精度数。
赞一下(1人) 回复此楼
8楼2016-01-09 21:20:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liming9398 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 319求调剂 +3 handrui 2026-04-05 3/150 2026-04-06 09:33 by jp9609
[考研] 一志愿C9的化学工程（085602） 340分，感觉校内调剂无望，求调剂 +12 万事宜臻 2026-04-04 12/600 2026-04-06 07:46 by 无际的草原
[考研] 本科211 分数293请求调剂 +4 莲菜就是藕吧 2026-04-01 4/200 2026-04-04 22:32 by hemengdong
[考研] 材料调剂 +10 懒羊羊轻置玉臀 2026-04-02 11/550 2026-04-04 21:56 by laoshidan
[考研] 一志愿上海大学生物学346 +3 上海大学346调剂 2026-04-03 3/150 2026-04-04 20:20 by dongzh2009
[考研] 一志愿华北电力大学（北京），材料科学与工程学硕265，求调剂 +11 yelck 2026-04-03 12/600 2026-04-04 19:52 by dongzh2009
[考研] 333求调剂 +9 阿科逸 2026-03-31 9/450 2026-04-04 18:25 by macy2011
[考研] 一志愿武理材料工程302调剂环化或化工 +19 Doleres 2026-03-31 20/1000 2026-04-04 16:44 by 啊俊！
[考研] 22408求调剂 354分 可跨专业 +3 hannnnnnn 2026-04-04 3/150 2026-04-04 14:35 by 土木硕士招生
[考研] 复试调剂 +6 范根培 2026-04-04 6/300 2026-04-04 14:27 by 土木硕士招生
[考研] 348分环境工程·调剂 +10 吴彦祖24k 2026-04-03 11/550 2026-04-04 14:19 by 无际的草原
[考研] 本9一志愿2 0854低分专硕286求调剂 +9 芒种111 2026-04-04 9/450 2026-04-04 11:01 by tangruihua
[考研] 282求调剂 +20 ycy1201 2026-04-01 22/1100 2026-04-04 00:42 by userper
[考研] 调剂0855-288 +5 x熊二a 2026-04-03 5/250 2026-04-04 00:19 by 猪会飞
[考研] 320求调剂 +5 振—TZ 2026-04-02 5/250 2026-04-03 14:42 by fxue1114
[考研] 一志愿安徽大学0817化学工程与技术，求调剂 +14 我不是只因 2026-04-02 15/750 2026-04-03 09:49 by 蓝云思雨
[考研] 0856材料与化工调剂，339 +14 10213207 2026-03-31 14/700 2026-04-02 21:01 by 1104338198
[考研] 283求调剂 +3 jiouuu 2026-04-02 4/200 2026-04-02 14:08 by 哒哒哒呱呱呱
[考研] 材料与化工（0856）304求B区调剂 +8 邱gl 2026-03-30 16/800 2026-04-01 17:58 by 邱gl
[考研] 求调剂 生物学 377分 +6 zzll03 2026-03-31 6/300 2026-03-31 17:33 by 唐沐儿
信息提示
关闭
请填处理意见
关闭