求助Float Invalid Operation

2332222 · 发表于 2017-4-12 17:41:02

程序最后没有出现错误和警告，但运行之后出现：“Float Invalid Operation” 然后黑色窗口闪了一下，就关了

2332222 · 发表于 2017-4-12 17:43:38

这是这个程序

fcode · 发表于 2017-4-12 21:15:45

1. 你的各种变量没有初始化。比如主程序 uu , ud , du , dd
2. 学会调试，你会自己解决各种问题。详解 http://v.fcode.cn/video-debugger.html
3. 根据我的调试，你的主要问题大概在
SUBROUTINE tsc(theta,m,k,t)
IMPLICIT NONE
REAL                :: d,PI
INTEGER,INTENT(in) :: m
REAL,INTENT(in)    :: theta,k
complex,INTENT(out)    :: t
  d=0.15
  PI=4.0*ATAN(1.0)
  t=SQRT(2.0*(COS(theta))/(k*d))*(SIN((k*SIN(theta)+m*PI/d)*d/2.0)/(SIN(theta)+m*PI/(k*d))&
&+SIN((k*SIN(theta)-m*PI/d)*d/2.0)/(SIN(theta)-m*PI/(k*d)))
END SUBROUTINE

当 theta 大于 PI/2 时，cos 为负数。开 sqrt 出错。

此外，117 行的  DO y=x,y+10,1
这是一个很奇怪的用法。

2332222 · 发表于 2017-4-13 22:17:59

fcode 发表于 2017-4-12 21:15
1. 你的各种变量没有初始化。比如主程序 uu , ud , du , dd
2. 学会调试，你会自己解决各种问题。详解 http ...

谢谢你的建议，我刚接触fortran，对它还不是太熟悉，非常感谢你能给我帮助

[数值问题] 求助Float Invalid Operation

美女勋章

热心勋章

星光勋章

新人勋章

贡献勋章

管理勋章

帅哥勋章

爱心勋章

规矩勋章

元老勋章

水王勋章