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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » numerical recipes in fortan 代码的疑惑
81/1返回列表
查看: 1993  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:37 by etmYJ6d2rquH
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:03 by etmYJ6d2rquH
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:46 by etmYJ6d2rquH
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-29 14:24 by etmYJ6d2rquH
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +6 IDs3scOF0tjC 2026-09-28 6/300 2026-09-29 14:13 by etmYJ6d2rquH
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意！ +8 zju2000 2026-09-26 9/450 2026-09-29 13:47 by ahsoarli
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[公派出国] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[有机交流] 同一个分子，一条来自文献，一条来自AI——不告诉你答案，你会选哪条？ +3 tianxiaoxian 2026-09-28 6/300 2026-09-28 22:44 by tianxiaoxian
[找工作] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
关闭
请填处理意见
关闭