Fortran读取可变阶数(n≤50)矩阵出现EOF运行时错误的解决方法

看起来你遇到的是程序预期读取的数据量和输入文件实际内容不匹配导致的EOF（文件结束）错误，这在Fortran的顺序文件读取里非常常见，我来帮你一步步排查解决：

1. 最可能的原因：输入文件缺少矩阵维度n的定义 #

从你的测试用例来看，你准备的是9个元素（对应3×3矩阵），但如果你的程序开头是先读取矩阵的维度n（比如代码里有read(9,*) n），那你的输入文件gaussj.inp必须先写出n的值，再跟上n²个矩阵元素。

举个例子，你的测试输入文件应该写成这样：

3
5.0 6.0 7.0 4.0 3.0 8.0 1.0 4.0 2.0

如果没写开头的3，程序会把第一个元素5.0当成整数n（自动截断为5），然后尝试读取5×5=25个元素，但文件里只剩8个，自然会触发EOF错误。

2. 核对你的读取代码逻辑 #

看你给出的代码片段，假设你的读取部分是类似这样的嵌套循环：

integer::i,j,n
integer,parameter::nmax=50
real,dimension(nmax,nmax)::A

open(unit=9, file='gaussj.inp', status='old')
read(9,*) n  ! 先读取维度n
do i=1,n
  do j=1,n
    read(9,*) A(i,j)  ! 逐元素读取矩阵
  end do
end do

这里的关键是：

• 必须确保输入文件开头的n值和后续元素数量一致（n²个）
• Fortran的列表导向读取（*格式）会自动跳过空格、换行，所以元素写在一行或多行都没问题，只要总数对。

如果你的程序没有先读n，而是直接用了固定值，那也会导致元素数量不匹配——比如硬编码n=5，但输入只有9个元素。

3. 快速调试：打印读取到的n值 #

在程序里加一行打印n的代码，就能快速确认问题：

read(9,*) n
print *, "程序读取到的矩阵维度n = ", n

如果输出的n是5，那就能100%确认是输入文件没写n的问题。

4. 可选方案：让程序自动计算矩阵维度 #

如果你不想在输入文件里写n，可以先读取所有元素到临时数组，再通过元素总数计算n（要求元素个数是完全平方数），代码示例：

integer::i,j,n,count,ios
integer,parameter::nmax=50
real,dimension(nmax,nmax)::A
real,dimension(nmax*nmax)::temp

open(unit=9, file='gaussj.inp', status='old')
count = 0
! 先读取所有元素到临时数组
do
  count = count + 1
  read(9,*,iostat=ios) temp(count)
  if(ios /= 0) exit  ! 遇到文件结束或错误就停止
end do
count = count - 1  ! 减去最后一次失败的读取

! 计算矩阵维度n
n = int(sqrt(real(count)))
! 检查元素个数是否为完全平方数
if(n*n /= count) then
  print *, "错误：输入的元素个数不是完全平方数，无法构成方阵"
  stop
end if

! 将临时数组转成n×n矩阵
do i=1,n
  do j=1,n
    A(i,j) = temp((i-1)*n + j)
  end do
end do

这种方式不需要输入n，程序自动适配，但要保证输入的元素个数刚好是某个整数的平方。

内容的提问来源于stack exchange，提问作者Rasheed Aadeez