变量定义上下文中的 INTENT(IN) 虚拟参数

Question

当我运行我的代码时，我在我的子例程中收到一条错误消息。

此代码来自 Kincaid 和 Cheney 关于高斯赛德尔方法求解椭圆偏微分方程的书中的练习。

错误信息是：

dummy argument 'u' with INTENT(IN) in variable definition context (assignment) at (1).

我参考下面代码中的（1）。我怎样才能修复子例程，使错误消息不出现？

subroutine seidel(ax,ay,nx,ny,h,itmax,u)     
  real, dimension(0:nx,0:ny), intent(in) :: u        
  real, intent(in) :: ax,ay, h                 
  integer, intent(in) :: nx, ny, itmax         
  integer:: i,j,k

  do  k = 1,itmax      
    do  j = 1,ny-1     
      y = ay + real(j)*h
      do  i = 1,nx-1   
        x = ax + real(i)*h      
        v = u(i+1,j) + u(i-1,j) + u(i,j+1) + u(i,j-1) 
        u(i,j) = (v - h*h*g(x,y))/(4.0 - h*h*f(x,y))           (1)
      end do
    end do
  end do
end subroutine seidel

Answer 1

intent(in) 是对编译器的承诺，即例程不会尝试更新参数的值。您显示的代码违反了该行中的承诺：

u(i,j) = (v - h*h*g(x,y))/(4.0 - h*h*f(x,y))

通过不违背你的承诺来解决这个问题，或者在这种情况下可能更合适，像这样inout这样

real, dimension(0:nx,0:ny), intent(inout) :: u

inout 告诉编译器例程将被传递参数并可能对其进行修改。

（我认为这可能是重复的，但还找不到。）

变量定义上下文中的 INTENT(IN) 虚拟参数

INTENT(IN) dummy argument in variable definition context

fortran

differential-equations

discretization