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环境。
求解释。

回复此楼

» 猜你喜欢

江苏省自然基金什么时候出结果已经有6人回复
面上有专家说实验设备不是我们单位的已经有3人回复
广西大学-广州大学招聘博士后欢迎广大优秀人才！！！已经有11人回复
闲聊已经有3人回复
如何查询2026年获批项目名称？已经有3人回复
两块石头已经有10人回复
科研人应该花精力去思考如何解决问题，而不是去凝练问题已经有19人回复
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有15人回复
面上没中，邀请各位路过的虫友分析一下分数已经有13人回复
初秋的晨风已经有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的回帖