当我运行以下简单程序时
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
显示
错误:无法分类的陈述
答案 0 :(得分:10)
您尚未声明子例程中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
答案 1 :(得分:2)
就像ifort提示那样(变量z)这个名字还没有被声明为数组或者函数.u需要将变量x,y,z声明为子程序multul中的数组。