gfortran优化标志下嵌套隐式do循环数据读取异常问询

问题现象 #

使用gfortran 12.2.0编译Fortran程序时，遇到如下差异表现：

• 无优化编译（gfortran Read_fct.f90 -o Read）：程序正常读取Input.txt中的8个浮点值，输出符合预期结果。
• 带-O1及以上优化编译（如gfortran -O2 Read_fct.f90 -o Read）：输出全为0.00000000，无错误或警告；添加-fcheck=bounds调试后，触发数组越界错误：

  At line 18 of file Read_fct.f90
  Fortran runtime error: Index '0' of dimension 1 of array 'noheight' below lower bound of 1

而使用Cray Fortran 16.0.1或ifort 2021.7.1以-O3编译时，程序可正常输出正确结果。

最小复现示例 #

程序代码（Read_fct.f90） #

program Read_fct
implicit none
integer ::  ntraj, mpuff
integer :: jheight, jt, jpuff
real,allocatable :: sttheight(:,:,:)

integer, dimension(:), allocatable :: noheight

open(1, file='Input.txt')
ntraj = 8
mpuff = 1

allocate(noheight(mpuff))
noheight(1) = 1

allocate(sttheight(noheight(1), ntraj, mpuff))

read (1,*) (((sttheight(jheight,jt,jpuff),jheight=1,noheight(jpuff)),jt=1,ntraj),jpuff=1,mpuff)
write(*,*) sttheight
close(1)
end program

输入文件（Input.txt） #

0.00000000,14370.01149861,14370.01149861,14370.01149861,14370.01149861,14370.01149861,14370.01149861,14370.01149861     ! sttheight

核心疑问与解答 #

1. 代码是否存在问题？ #

代码存在未定义行为。问题出在嵌套隐式DO循环的上限表达式noheight(jpuff)：该表达式引用了外层循环变量jpuff，但Fortran标准并未强制规定隐式DO循环上限的求值时机。gfortran在优化时会提前计算所有循环的上限值，此时jpuff尚未进入循环迭代（值未定义，通常为0），导致访问noheight(0)触发越界；而其他编译器选择在每次外层循环迭代时计算内层上限，因此能正常运行。这种依赖编译器特定行为的写法，本质上是可移植性缺陷。

2. 这是gfortran优化标志下的预期行为吗？ #

这是gfortran优化策略导致的结果，但不属于"预期行为"范畴——它是代码未定义行为被优化策略暴露后的表现。gfortran的优化器为提升性能，会提前计算循环上限等可预判的表达式，但这也会暴露代码中依赖未初始化变量的问题。

3. 这是gfortran的bug吗？ #

不算gfortran的bug。Fortran标准允许编译器在未定义行为场景下采用不同处理逻辑，gfortran的优化行为符合标准的灵活性规定。问题根源在于代码写法依赖了非标准的编译器行为，而非编译器本身的错误。

修复建议 #

将嵌套隐式DO循环改为显式DO循环，确保内层循环上限在每次外层迭代时正确求值，避免未定义行为：

program Read_fct
implicit none
integer ::  ntraj, mpuff
integer :: jheight, jt, jpuff
real,allocatable :: sttheight(:,:,:)

integer, dimension(:), allocatable :: noheight

open(1, file='Input.txt')
ntraj = 8
mpuff = 1

allocate(noheight(mpuff))
noheight(1) = 1

allocate(sttheight(noheight(1), ntraj, mpuff))

! 替换隐式DO为显式循环
do jpuff = 1, mpuff
  do jt = 1, ntraj
    read(1,*) (sttheight(jheight, jt, jpuff), jheight=1, noheight(jpuff))
  end do
end do

write(*,*) sttheight
close(1)
end program

内容的提问来源于stack exchange，提问作者Kristian HM