为什么case语句中，没有的例子，值会很大？

芥子纳须弥 · 发表于 2022-5-26 12:13:31

本帖最后由芥子纳须弥于 2022-5-26 12:14 编辑

[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode

program select
   INTEGER medium,isal
   real ew
   
   medium=1
   isal=0
   call DIEL_WAT(medium, isal, ew)
   write(*,*) ew
   write(*,*) isal
  
  end 
  


  subroutine DIEL_WAT(medium, isal, ew)
  INTEGER medium,isal,a,b
  real ew

  SELECT CASE (isal)
  CASE ( 1 )
    a =20
  CASE ( 2 )
    a =40
  endselect
  
    SELECT CASE (medium)
  CASE ( 0 )
    b =0
  CASE ( 1 )
   b=1
  CASE ( 2 )
   b=2
  endselect
 
  ew = a+b
  end

得到的结果是

vvt · 发表于 2022-5-26 15:24:19

因为变量 a， b 没有初始化，其值不确定。

vvt · 发表于 2022-5-26 15:29:17

如果你的 isal , medium 可能的取值是连续的。可以尝试用数组表来代替 select

[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode

program select
  integer :: medium = 1 , isal = 1
  real :: ew
  call DIEL_WAT(medium, isal, ew)
  write(*,*) ew, isal
end program select

subroutine DIEL_WAT(medium, isal, ew)
  Integer , parameter :: AList(0:*) = [0,20,40]
  Integer , parameter :: BList(0:*) = [0,1,2]
  iNTEGER medium,isal
  real ew
  ew = AList(isal)+BList(medium)
end subroutine DIEL_WAT

芥子纳须弥 · 发表于 2022-5-26 21:16:25

vvt 发表于 2022-5-26 15:29
如果你的 isal , medium 可能的取值是连续的。可以尝试用数组表来代替 select
[mw_shl_code=fortran,true]p ...

感谢，解决了我的问题

zjk0112 · 发表于 2022-5-26 23:01:45

我用vs2022+oneAPI的Fortran编辑器跑mkl里的例子就是遇到过这样的情况。没有转置时，索引的值是0，一转置就变成了几千万。最后发现就是没有初始化，没给赋初值的原因。

芥子纳须弥 · 发表于 2022-7-18 22:25:11

zjk0112 发表于 2022-5-26 23:01
我用vs2022+oneAPI的Fortran编辑器跑mkl里的例子就是遇到过这样的情况。没有转置时，索引的值是0，一转置就 ...

是的,没错

[流程控制] 为什么case语句中，没有的例子，值会很大？

规矩勋章

元老勋章

新人勋章

水王勋章

热心勋章