Gfortran与Intel Fortran接口混用子程序和函数的差异及解决方法

Fortran 官方标准（从F90到最新的F2018）明确要求：同一个通用接口下绑定的过程必须全部为子例程（subroutine）或全部为函数（function），禁止混合绑定。你遇到的报错是gfortran严格执行标准的正常表现；Intel Fortran默认开启了私有的非标准扩展，允许混合绑定，所以可以直接编译运行。

总改动不超过5行，完全保留核心随机数生成逻辑，兼容所有Fortran编译器：

1. 把原代码中生成单个随机数的函数runif_r1改成子例程，将返回值改为输出参数，和数组版本的子例程类型保持一致
2. 通用接口块无需任何修改，编译器会根据第三个参数是标量还是数组自动匹配对应的子例程，不会出现重载歧义
3. 主程序中原来用赋值方式调用标量版本的代码，改成call语句调用即可

修改后可直接编译运行的完整代码 #

module stats

implicit none
public :: runif
private :: runif_r, runif_r1
interface runif
    module procedure runif_r, runif_r1
end interface 

contains

!get uniform random vector using subroutine
subroutine runif_r(a,b,R) 
implicit none
real, intent(in) :: a,b 
real, intent(out) :: R(:)
call random_number(R)
R = a+(b-a)*R
end subroutine runif_r

! 原函数改为子例程，返回值调整为输出参数，核心逻辑不变
subroutine runif_r1(a,b,r) 
implicit none
real, intent(in) :: a,b 
real, intent(out) :: r
call random_number(r)
r = a+(b-a)*r
end subroutine runif_r1

end module stats


program test

use stats ! 取消原注释，确保正确引用模块过程
integer, parameter :: n=5
real :: x(n), x1

call runif(1.0, 3.0, x)
print*, x

! 原赋值调用改为call语句
call runif(1.0, 3.0, x1) 
print*, x1 

end program test

如果你一定要保留x1 = runif(1.0,3.0)的函数式调用写法，只需要把两个过程都统一为函数即可：数组版本改为返回数组的函数，标量版本保留原函数形式，主程序中数组生成的call语句改为赋值语句，改动量同样很小。

内容的提问来源于stack exchange，提问作者John Stone