Fortran中传递至外部函数后二维数组为空的问题排查

问题根源 #

你遇到的数组大小异常问题，核心原因是外部函数缺少显式接口。

当你在主程序中用real(real64), external :: test_function声明外部函数时，编译器只能依赖隐式接口处理函数调用，但你的函数参数是假定形状数组（u(:,:)）——这种类型的参数需要编译器知道数组的维度、形状等元信息才能正确传递，隐式接口无法提供这些信息，导致数组的形状数据在传递时丢失，所以函数内部读取到的数组大小变成了错误的1x1。

解决方案 #

方案1：使用模块（推荐的现代Fortran写法） #

将函数封装到模块中，主程序通过use语句引入模块，编译器会自动生成并识别显式接口，这是Fortran中管理跨代码单元接口的标准方式。

步骤1：创建模块文件func_mod.f90

module func_mod
    use iso_fortran_env
    implicit none
contains
    function test_function(u) result(res)
        real(real64), intent(inout):: u(:,:)
        real(real64) res

        print *, "Entering test_function"
        print *, "Size of u: ", size(u, 1), size(u, 2)
        print *, "Contents of u (first few values): ", u(1:5, 1:5)

        res = 0.
    end function test_function
end module func_mod

步骤2：修改主程序program_test.f90

program test
    use iso_fortran_env
    use func_mod  ! 引入模块，自动获取函数的显式接口
    implicit none
    real(real64), allocatable :: u(:, :)
    real(real64) :: res

    allocate(u(17, 17))
    if (.not. allocated(u)) then
        print *, "Memory allocation failed for u or f."
        stop
    end if
    print *, "Before entering test_function"
    print *, "Size of u: ", size(u, 1), size(u, 2)
    print *, "Contents of u (first few values): ", u(1:2, 1:2)
    u = 0.

    res = test_function(u)

end program test

编译运行

保持原编译命令不变：

gfortran program_test.f90 func_mod.f90 -o program_test.out -Wall && ./program_test.out

方案2：手动编写显式接口块 #

如果不想使用模块，可以在主程序中为外部函数手动编写显式接口块，明确告诉编译器函数的参数类型、形状和返回值：

修改主程序program_test.f90：

program test
    use iso_fortran_env
    implicit none
    real(real64), allocatable :: u(:, :)
    real(real64) :: res

    ! 为外部函数定义显式接口
    interface
        function test_function(u) result(res)
            use iso_fortran_env
            implicit none
            real(real64), intent(inout):: u(:,:)
            real(real64) res
        end function test_function
    end interface

    allocate(u(17, 17))
    if (.not. allocated(u)) then
        print *, "Memory allocation failed for u or f."
        stop
    end if
    print *, "Before entering test_function"
    print *, "Size of u: ", size(u, 1), size(u, 2)
    print *, "Contents of u (first few values): ", u(1:2, 1:2)
    u = 0.

    res = test_function(u)

end program test

编译运行命令不变，即可正常传递数组信息。

关键说明 #

Fortran中，以下特性必须依赖显式接口才能正常工作：

• 假定形状数组（u(:,:)）
• 可选参数
• 关键字参数
• 数组指针或分配数组作为参数
• 返回值为数组或派生类型的函数

使用external声明只能告知编译器这是一个外部函数，但无法提供上述特性所需的元信息，因此会导致参数传递异常。

内容的提问来源于stack exchange，提问作者andreas De Keersmaeker