fortran - 错误：fortran 中的不可分类语句

Question

当我运行以下简单程序时

program test
! integer m,n,r,i
double precision x(2),y(3),z(4)
x=(/2.0,1.0/)
y=(/1.0,2.0,1.0/)
call polymul(x,2,y,3,z,4)
print *,z
 end

subroutine polymul(x,m,y,n,z,r)
! polynominal multipy
integer i,j,k
do i=1,r
z(i)=0.0
end do
do i=1,m
  do j=1,n
    k=i+j-1
    z(k)=z(k)+x(i)*y(j)
  end do
end do
end

这显示了

错误：无法分类的陈述

Answer 1

您尚未在子例程中声明什么x、、y和z。Fortran 不知道这些变量是函数（尚未定义）还是数组。修复很简单：在子例程中显式声明数组：

    subroutine polymul(x, m, y, n, z, r)
       implicit none
       integer m, n, r
       double precision x(m), y(n), z(r)
       integer i, j, k
       do i=1,r
          z(i)=0.0
       enddo
       do i=1,m
          do j=1,n
             k=i+j-1
             z(k)=z(k)+x(i)*y(j)
          enddo
       enddo
    end subroutine

Answer 2

就像ifort提示的那样（变量z）这个名字没有被声明为数组或者函数。u需要在子程序polymul中声明变量x,y,z为数组。