Fortran 数组无法在函数 : not a DUMMY variable 中返回

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作为内部函数意味着您不必在主程序中声明它。用行话来说:它有一个明确的接口(interface)。

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

函数子程序的定义

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 .通常，虚拟参数是 ( 之间的那些东西。和 ) ,