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)

这个声明是不正确的.您可以声明type(real(dp))和dimension((3,3)),但这intent(out)是错误的.

Fortran标准中的intent属性受约束(C538)

具有INTENT属性的实体应为虚拟数据对象或虚拟过程指针.

回到编译器消息,m31tensorprod不是虚拟变量.在这种情况下,伪参数是a和b.一般来说,伪参数是(和之间的事物),