错误: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