当我运行以下简单程序时
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它显示了
错误:不可分类的语句
发布于 2013-06-05 01:42:36
您还没有声明子例程中有什么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发布于 2013-06-06 03:13:12
正如ifort提示(变量z)这个名称没有被声明为数组或函数一样,u需要在子例程polymul中声明变量x、y、z为数组。
https://stackoverflow.com/questions/16930096
复制相似问题