不能传递子程序名字

pilipili2021

env：ubuntu 20.04+ifort
CODE:

[Fortran] 纯文本查看 复制代码

查看代码

打印?

01	program  main

02

CALL B()

03	CALL FKB(B)

04

05

06

contains

07	subroutine B()

08	IMPLICIT NONE

09	write (*,*) 'fk'

10	end subroutine B

11

12	subroutine FKB(fcn)

13	implicit none

14

INTERFACE

15	SUBROUTINE fcn()

16	IMPLICIT NONE

17	END SUBROUTINE fcn

18	END INTERFACE

19	write (*,*) '1'

20

call B()

21	write (*,*) '2'

22	call fcn()

23	write (*,*) '3'

24	end subroutine FKB

25

26

end

很简单的程序，把子程序名字‘B' 通过'fcn'传递给子程序'FKB'里，编译也不报错，但是运行以后发现

ifort -g main.F90
./a.out

fk
1
fk
2
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image              PC                Routine            Line        Source
a.out              000000000040490A  Unknown               Unknown  Unknown
libpthread-2.31.s  00007F86D9BE53C0  Unknown               Unknown  Unknown
Unknown            00007FFFF6983F50  Unknown               Unknown  Unknown


显然调用'fcn'的时候出现了问题。
请问这到底是为什么？如何解决？

pilipili2021

或者有人能帮我试着跑一下可以么，看看报不报错

kyra

gfortran处理主程序contains下的函数有特殊的方式。
你把 subroutine B 放到 program 外面去就行了，或者放到一个module里。

另外的问题，建议用 Procedure 代替子程序内的 interface

[Fortran] 纯文本查看 复制代码

查看代码

打印?

01	Module fcnMod

02

Interface

03	Subroutine IF_fcn()

04	End Subroutine IF_fcn

05	End Interface

06

contains

07	subroutine B()

08	write (*,*) 'fk'

09	end subroutine B

10	End Module fcnMod

11

12	program  main

13	use fcnMod

14	IMPLICIT NONE

15

16

CALL B()

17	CALL FKB(B)

18

19

contains

20

21	subroutine FKB(fcn)

22	Procedure(IF_fcn) :: fcn

23	write (*,*) '1'

24

call B()

25	write (*,*) '2'

26	call fcn()

27	write (*,*) '3'

28	end subroutine FKB

29

30

end

或者

[Fortran] 纯文本查看 复制代码

查看代码

打印?

01	Subroutine C()

02	write (*,*) 'fc'

03	End Subroutine C

04

05	program  main

06	IMPLICIT NONE

07

Interface

08	Subroutine IF_fcn()

09	End Subroutine IF_fcn

10	End Interface

11	Procedure(IF_fcn) :: C

12

CALL C()

13	CALL FKB(C)

14

15

contains

16

17	subroutine FKB(fcn)

18	Procedure(IF_fcn) :: fcn

19	write (*,*) '1'

20

call C()

21	write (*,*) '2'

22	call fcn()

23	write (*,*) '3'

24	end subroutine FKB

25

26

end

pilipili2021

非常有启发性的建议。
实际上我写的东西是一个鲍威尔寻路解非线性方程组的MODULE，所以必须要传递函数名字。
MODULE里面再加MODULE，这个我不确定，但是的确提供了一种思路。
至于用Procedure去替代interface 我都没想过。。。。
太感谢你了，我好好鼓捣鼓捣

kyra

不需要module里再加module，在module里（而不是program）就行。可以在同一个module，或者另一个module。

Fortran 应该是不允许module里再写module的。但是module允许有SubModule