3

如果引用数组元素的索引实际上超出了它的假定范围，我很困惑 Fortran 如何处理这种情况。

下面是一个简单的代码来说明这个问题：

PROGRAM test_matrix_out

USE mod_writearray

IMPLICIT NONE
INTEGER :: i,j,m,n
REAL    :: k
REAL, Dimension(:,:),ALLOCATABLE :: A

m = 3
n = 4
ALLOCATE(A(m,n))

k = 1

DO i=1,m
    DO j=1,n
        A(i,j)=k
        k=k+1
    ENDDO
ENDDO

CALL writearray(A)
WRITE(*,*)
WRITE(*,*) A(1,:)
WRITE(*,*)
WRITE(*,*) A(2,:)
WRITE(*,*)
WRITE(*,*) A(0,:)
WRITE(*,*)
WRITE(*,*) A(4,:)
WRITE(*,*)
WRITE(*,*) A(5,:)
WRITE(*,*)
WRITE(*,*) A(100,:)
WRITE(*,*)
WRITE(*,*) A(:,1)
WRITE(*,*)
WRITE(*,*) A(:,2)
WRITE(*,*)
WRITE(*,*) A(:,0)
WRITE(*,*)
WRITE(*,*) A(:,4)
WRITE(*,*)
WRITE(*,*) A(:,5)
WRITE(*,*)
WRITE(*,*) A(:,100)


DEALLOCATE(A)

END PROGRAM test_matrix_out

它给了我以下结果：

   1.000000       2.000000       3.000000       4.000000

   5.000000       6.000000       7.000000       8.000000

  0.0000000E+00   9.000000       10.00000       11.00000

   2.000000       3.000000       4.000000      0.0000000E+00

   6.000000       7.000000       8.000000      0.0000000E+00

  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00

   1.000000       5.000000       9.000000

   2.000000       6.000000       10.00000

 -1.0097448E-28  8.9776148E-39  0.0000000E+00

   4.000000       8.000000       12.00000

  0.0000000E+00  0.0000000E+00  0.0000000E+00

 -3.3631163E-44  1.4293244E-43  0.0000000E+00

为什么会这样？

4

2 回答 2

7

当您编写 A(i,j) 时，编译器会计算该内存位置的地址。例如，参见http://en.wikipedia.org/wiki/Array_data_structure#Multidimensional_arrays。编译器通常不会根据语言规则确定这是否是合法地址。使用超出维度的索引是非法的。程序员有责任不这样做。Fortran 的优点之一是能够为该错误添加运行时检查。传统的传说是运行时下标检查很昂贵，但我在测试时经常发现运行时成本可以忽略不计，有时会在程序的生产版本中保留它。

如果您正在读取内存，索引错误的可能后果将是获得错误的值，除非内存位置太远以至于它超出了属于程序的内存，这将产生错误。如果您正在写入内存，您将破坏数组中其他地方的内存，属于其他一些变量，或者属于程序的内部数据结构。看看没有释放会导致什么样的问题？对于索引错误导致程序运行时问题的示例问题。

于 2012-07-19T06:54:25.530 回答
5

您所看到的是，您用于编译程序的编译器并未在运行时检查数组是否越界。因此，取决于编译器和机器，任何事情都会发生。有时，可能会引用未在内存中显式分配的数组元素，这就是您的示例中发生的情况。在这种情况下，越界元素的值是程序运行时位于该内存地址的任何值。如果请求的内存地址不存在或无法访问，程序将因分段错误而失败。

我只是在查看当前 Fortran 标准的草案，我找不到任何关于越界访问数组元素是否是已定义行为的声明。为了避免这些问题，请使用-C（检查边界）标志编译您的程序。如果可能，该程序将让您知道哪个数组的哪个元素超出了范围。在开发期间使用-C，但不要在生产中使用，因为它会极大地减慢代码速度。

此外，为了将来参考，在询问此类问题时（例如，为什么我的程序会输出这个？），最好包括有关正在使用的编译器（带有版本号）和目标体系结构的信息。

于 2012-07-19T04:42:15.427 回答