24小时热门版块排行榜    

Znn3bq.jpeg
查看: 532  |  回复: 9
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

gongxd325

至尊木虫 (文坛精英)

[交流] 【求助】一句Fortran语句出错问题(已解决)

编个小程序计算指前因子和速率常数,有个语句出错,向大家请教(kB/h=0..20827*10**11)
      A=0.20827*298.15*EXP(2-13.288/8.314)*10**11
程序计算结果是 0.113E+12,正确的结果应该是 9.28E+12
假如只计算 0.20827*298.15*EXP(2-13.288/8.314) 得到 92.8,乘以10**11结果就正确。
想不明白,不知道为什么?

[ Last edited by cadick on 2009-12-14 at 02:53 ]
回复此楼
一个人做点好事并不难难的是一辈子做好事!!!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by tjpm at 2009-6-1 20:54:
有兴趣的大家可以试试这个

program xxx
real*8 i
i=0.1
write(*,*) "i=",i
end program xxx

使用gfortran和ifort都是
输出i=0.00000001490116

[ Last edited by tjpm on 2009-6-1 at 20 ...

你就是用 i=0.1d0 得到的结果也不可能恰好等于 0.1 啊……
9楼2009-06-01 22:45:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

snoopyzhao

至尊木虫 (职业作家)


小木虫(金币+0.5):给个红包,谢谢回帖交流
CODE:
A=0.20827*298.15*EXP(2-13.288/8.314)*10.0**11

其实改动很小,就是将 10 改成 10.0。如果不改的话,gfortran -Wall 会给出如下提示:
CODE:
      A=0.20827*298.15*EXP(2-13.288/8.314)*10**11
                                            1
Error: Arithmetic overflow at (1)

[ Last edited by snoopyzhao on 2009-6-1 at 11:26 ]
2楼2009-06-01 11:24:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gongxd325

至尊木虫 (文坛精英)

谢谢!
Compaq Visul Fortran不报错,但结果错误,搞了半天。
一个人做点好事并不难难的是一辈子做好事!!!
3楼2009-06-01 11:42:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tjpm

金虫 (正式写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
q68(金币+2,VIP+0):多多指教! 6-1 18:32
你这样写程序会有问题的
建议改一下
A=0.20827d0*298.15d0*EXP(2.0d0-13.288d0/8.314d0)*10.0d0**11

如果是real*16的 就改d0成q0.
引用回帖:
Originally posted by gongxd325 at 2009-6-1 11:13:
编个小程序计算指前因子和速率常数,有个语句出错,向大家请教(kB/h=0..20827*10**11)
      A=0.20827*298.15*EXP(2-13.288/8.314)*10**11
程序计算结果是 0.113E+12,正确的结果应该是 9.28E+12
假如只计算  ...

4楼2009-06-01 16:37:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +6 yexuqing 2026-04-19 7/350 2026-04-23 12:32 by yexuqing
[基金申请] 国自然面上和省基金B类撒花 +18 花田半亩~白 2026-04-21 18/900 2026-04-23 11:31 by 12021227
[考研] 有没有学校收留 +3 蒋昌鹏qtj 2026-04-20 3/150 2026-04-22 20:25 by 学员JpLReM
[考研] 312求调剂 +3 山河似你温柔 2026-04-22 3/150 2026-04-22 20:17 by 学员JpLReM
[考博] 华师大读博 +3 xq83 2026-04-22 5/250 2026-04-22 10:42 by xq83
[论文投稿] 急需审稿人!!! +3 陆小果画大饼 2026-04-21 3/150 2026-04-21 23:54 by jzy_123456
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 8/400 2026-04-21 16:26 by 啃面包的小书虫
[考研] 295分求调剂 +6 ?要上岸? 2026-04-17 6/300 2026-04-21 08:18 by Equinoxhua
[考研] 085600材料与化工调剂 5+3 孜孜不倦2002 2026-04-19 6/300 2026-04-20 21:25 by babero
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
信息提示
请填处理意见