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中的m31tensorprod属性受限制（C538）

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

回到编译器消息，a不是虚拟变量。在这种情况下，伪参数是b和(。通常，伪参数是)和{{1}}之间的事物，

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

1 个答案: