24小时热门版块排行榜

返回列表

guojunhang

金虫 (正式写手)

应助: 7 (幼儿园)
金币: 3820.6
帖子: 516
在线: 174.2小时
虫号: 1438713

[交流] 问一个关于函数调用的小问题

请教个问题，看看我这个函数，fortran77写的：
   FUNCTION HEAVI(PRIS)
         REAL*8 PRIS
   IF(PRIS.GT.0.0) THEN
      HEAVI=1.0
   ENDIF
   RETURN
C
      IF(PRIS.LT.0.0) THEN
      HEAVI=0.0
      ENDIF
   RETURN
      END

调用的时候出现下列情况：
HEAVI(-3)=1.00000000000000，
HEAVI(-3.0D0)=0.000000000000000E+000，
   HEAVI(-3.0)=1.00000000000000
HEAVI(3)=1.00000000000000，
HEAVI(3.0D0)=1.000000000000000E+000，
HEAVI(3.0)=1.00000000000000
不知道什么原因？

改成这样怎么就好了呢？

   Program ex005
   real HEAVI
   write(*,*) HEAVI(3.0d0)
   end

   FUNCTION HEAVI(PRIS)
   REAL PRIS
   Real HEAVI
   IF(PRIS.GT.0.0) THEN
      HEAVI=1.0
   else
      HEAVI=0.0
   ENDIF
   RETURN
   END
语言Fortran77，gfortran环境。
求解释。

回复此楼

» 猜你喜欢

一篇论文同时出现在两个期刊，一模一样，这算不算学术不端，请各位老师斧正。已经有12人回复
希望面上有个好结果已经有7人回复
有谁可曾问过你过的还好吗？已经有15人回复
E0414, 我的本子有没有希望？已经有5人回复
今年也是没消息就是没中么已经有16人回复
三区计算机方向期刊推荐已经有5人回复
申博已经有5人回复
sci论文二审求助已经有5人回复
函评已经有7人回复
以CTAB等为模板水热合成的产物，如何除去模板？已经有6人回复

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2011-11-11 21:53:59

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
dubo(金币+1): 欢迎讨论 2011-11-11 23:31:16
guojunhang(金币+1): 这个问题我也注意到了，所以改正后了程序流程 2011-11-13 18:27:43

你第一个程序中第一个 return 直接就结束了那个函数，后面的那个条件永远不会被执行啊……

赞一下(2人)

回复此楼

2楼2011-11-11 22:53:36

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

maomao1210

金虫 (正式写手)

程序强帖: 5
应助: 2 (幼儿园)
金币: 1432.8
帖子: 991
在线: 442.1小时
虫号: 253215

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
jjdg(金币+1): 感谢参与 2011-11-12 00:41:19
余泽成: 专家考核存档 2011-11-13 21:35:24

引用回帖:

1楼: Originally posted by guojunhang at 2011-11-11 21:53:59:
请教个问题，看看我这个函数，fortran77写的：
   FUNCTION HEAVI(PRIS)
         REAL*8 PRIS
   IF(PRIS.GT.0.0) THEN
      HEAVI=1.0
   ENDIF
   RETURN
C
      IF(PRIS. ...

貌似变量定义也有问题，第一个函数中用的是real*8,主程序给的有时候是real*8 ,有时候是real*4

[ Last edited by maomao1210 on 2011-11-11 at 23:00 ]

赞一下(3人)

回复此楼

3楼2011-11-11 22:55:13

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖