24小时热门版块排行榜

>论坛更新日志 (4681)
>虫友互识 (599)
>导师招生 (527)
>文献求助 (339)
>找工作 (291)
>论文投稿 (273)
>职场人生 (245)
>休闲灌水 (172)
>考博 (165)
>考研 (108)
>硕博家园 (102)
>博后之家 (91)
>学术会议 (90)
>招聘信息布告栏 (74)
>医学 (64)
>绿色求助(高悬赏) (41)

返回列表

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

回复此楼

» 猜你喜欢

关于核磁制样已经有3人回复
郑州大学申博已经有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的回帖