Fortran 切片矩阵输出错误

我试图为n*n矩阵编写一个非常简单的代码，以生成切片矩阵。然而，当使用Geany和Gfortran运行时，我得到了错误的切片矩阵

代码是：

program mat_slice

implicit none
integer, parameter:: n = 3
integer, dimension(n,n) :: mat
integer, dimension(n-1,n-1) :: slice
logical, dimension(n,n) :: mask
integer :: i,j
mat = reshape((/1,2,3,4,5,6,7,8,9/),(/n,n/))

print*, mat

mask = .true.
mask(i,:) = .false.
mask(:,j) = .false.

slice = reshape((pack(mat,mask)),(/n-1,n-1/))

do i  =  1,n-1
    do  j = 1,n-1
        print*, slice(i,j)
    end do
end do

端程序块

输出是 0 4. 0 5

然而，当我在Jdoodle上运行相同的代码时，我得到了正确的答案，即

输出是

1 4. 2. 5

---

如果您有任何想法，我们将不胜感激。

您的

i

和

j

值在首次使用时未定义

mask(i,:) = .false. 
mask(:,j) = .false.

因此，您的程序是无效的

当程序访问数组时，可能会使用值0。始终首先编译程序，并进行所有调试检查和警告。在gfortran中，使用

-g-fbacktrace-Wall-fcheck=all

启动。有关更多信息，请参阅手册