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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » numerical recipes in fortan 代码的疑惑
81/1返回列表
查看: 1380  |  回复: 7
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者liming9398将赠送您 2 个金币

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返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫 SCI英文润色翻译-EnPapers美团队-先修改后付款
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 豫北虫友互识 +10 xuhongli903 2024-06-18 11/550 2024-06-20 00:25 by xuhongli903
[基金申请] 江南大学到瑞士招聘，称取消非升即走，改预聘+长聘 +21 babu2015 2024-06-18 22/1100 2024-06-19 23:03 by feng6531
[基金申请] F口信息学部拿面上，大概需要什么样的成果 +7 _奋黎_ 2024-06-16 18/900 2024-06-19 20:02 by _奋黎_
[考博] 34岁读博士晚吗 +39 emitdne 2024-06-13 39/1950 2024-06-19 19:03 by Dqtjs
[基金申请] 面上博士后 +3 jsqy 2024-06-19 3/150 2024-06-19 17:58 by bxwisdom_886
[论文投稿] 审稿人含糊拒稿，还需要回复吗？如何回复？ 20+4 BruceChum 2024-06-15 22/1100 2024-06-19 08:00 by kanyechris
[精细化工] 试剂价格 +4 旭必上岸 2024-06-17 4/200 2024-06-18 23:58 by 又幂青松顶
[找工作] 初始合伙人来啦！（生物试剂耗材标准品） +14 欢快的小科研人 2024-06-15 25/1250 2024-06-18 20:35 by 小飞来虫
[教师之家] 请问事业编制和年薪制冲突吗？ +12 ZHONGWU_U 2024-06-14 12/600 2024-06-18 19:31 by fangyl2005
[论文投稿] 投稿求助 10+3 2022_灵魂工程师 2024-06-16 8/400 2024-06-18 12:12 by 投必得科研顾问
[高分子] 烧瓶内合成聚酯 +3 大帝国乐 2024-06-17 7/350 2024-06-18 11:10 by 大帝国乐
[论文投稿] 论文已接收，但发现修改稿传的是旧版该怎么办？ +3 zsq54321 2024-06-15 6/300 2024-06-17 19:45 by wjykycg
[催化] 镍负载氧化铝的保存问题 8+3 lwn0130 2024-06-15 4/200 2024-06-17 10:48 by adaihao
[基金申请] 博士后创新人才支持计划公示 +9 aishida144 2024-06-14 15/750 2024-06-16 09:52 by msjy
[基金申请] 关于博后基金的bug问题 +6 lxr1991 2024-06-14 9/450 2024-06-15 21:17 by since—2010
[论文投稿] 投稿时忘记修改一作 +7 gll123456 2024-06-13 11/550 2024-06-15 11:49 by gll123456
[论文投稿] 投了一篇4区的SCI，审稿人一个拒稿，一个小修，编辑给了大修。 +9 安稳22123 2024-06-13 10/500 2024-06-14 23:45 by jurkat.1640
[有机交流] ππ堆积会发生在有机溶剂中吗 5+3 zibuyu0420 2024-06-13 4/200 2024-06-14 14:17 by 小肉干
[基金申请] 工材E10口函评结束了吗 10+3 我1的飞翔 2024-06-13 5/250 2024-06-14 06:35 by nono2009
[基金申请] 博士后面上项目状态还是专家评审吗 10+9 Thatcheremu 2024-06-13 55/2750 2024-06-13 21:23 by 乌合麒麟
信息提示
关闭
请填处理意见
关闭