Fortran 数组不能在函数中返回：不是 DUMMY 变量

Question

作为 Fortran 90 自由格式的新手，我真的很想知道为什么以下代码片段不起作用：

program test2
    implicit none
    !!! A program to practice f90 writing.
    ! Define double precision data
    integer, parameter :: dp = kind(1.d0)
    real(dp) :: a(3), b(3)
    integer :: i
    a = (/(i, i=1, 3)/)
    b = (/(i, i=1, 3)/)
    write (*, *) m31tensorprod(a, b)

contains
    function m31tensorprod(a, b)
        real(dp), dimension(3), intent(in) :: a, b
        real(dp), intent(out) :: m31tensorprod(3, 3)
        integer :: k1, k2
        forall(k1=1:3, k2=1:3)
            m31tensorprod(k1, k2) = a(k1) * b(k2)
        end forall
        return
    end function m31tensorprod
end program test2

当我尝试通过 gfortran test2.f90 编译它时，它说：

test2.f90:13.4:

function m31tensorprod(a, b)
1 Error: Symbol at (1) is not a DUMMY variable

我想因为m31tensorprod是一个内部函数，所以不需要声明。我哪里做错了？

谢谢，

Answer 1

你是对的，m31tensorprod是一个内部函数意味着你不必在主程序中声明它。用行话来说：它有一个显式接口。

但是，这不是您的代码的问题。函数定义本身出了问题。 [诚然，编译器消息并不太有用。]

函数子程序的定义

function m31tensorprod(a, b)

用结果变量定义一个函数m31tensorprod。此结果变量以您的声明为准

    real(dp), intent(out) :: m31tensorprod(3, 3)

这个声明是不正确的。您可以声明类型 (real(dp)) 和维度 ((3,3))，但 intent(out) 是错误的。

intent属性，用Fortran standard的话来说，就是受约束（C538）

An entity with the INTENT attribute shall be a dummy data object or a dummy procedure pointer.

回到编译器消息，m31tensorprod 不是虚拟变量。在这种情况下，伪参数是 a 和 b。一般来说，虚拟参数是 ( 和 )、

之间的那些东西

Fortran 数组不能在函数中返回：不是 DUMMY 变量

Fortran array cannot be returned in function: not a DUMMY variable

fortran

gfortran

fortran90