一维数组的分割错误

Question

在这里，我知道 A 只是一个长度为 170 的一维数组。

subroutine gonewrong(q,e,A)

implicit none

integer, dimension(:) :: q

integer, dimension(:,:), intent(in) :: e

integer, dimension(size(e,1)*2), intent(out) :: A
print *, A(2)
end subroutine gonewrong

例如，当我尝试查找 A(2) 时，它给了我一个分段错误！ e 具有维度 (85,2)。

会不会是因为我在用于编译它的程序中将 A 声明为可分配数组？

我使用的程序：

program prog_1
use module_1
implicit none
integer::qm,a,b,c
integer,allocatable,dimension(:)::q
integer,allocatable,dimension(:,:)::e
!integer, dimension(size(e,1)*2) :: A
a = 5
b = 3
c = 10
allocate(q(a+c))
allocate(e(a+b*c,2))

call subr1(a,b,c,qm,q,e) !Outputs are qm,q and e.
call gonewrong(q,e,A) !gonewrong takes q and e as arguments
end program prog_1

Answer 1

另一个需要考虑的问题是在声明中：

integer,allocatable,dimension(:,:)::e
integer, dimension(size(e,1)*2) :: A

e 的大小在声明 A 时是未定义的（零？）。如果A 的大小取决于e 的大小，您需要在编译时修复e 的大小，或者将A 设置为可分配并在@987654328 的大小后单独分配它@ 已经设置好了。要点是A 不会随着e 的大小变化而动态调整其大小。在这种情况下，看起来A 声明的长度为零，因为size(e) 返回零，因为只定义了e 的形状，而不是它的大小。引用 A(2) 可能会引发分段违规（读取数组末尾），这与您观察到的内容一致。

除此之外：许多从动态语言进入 Fortran 的人没有意识到 Fortran 例程至少分为三个“区域”；包含用于导入外部资源的use 语句的初始部分，声明参数和局部变量的声明部分，以及放置可执行语句的例程主体（鼓励学究阅读语言规范；我在这里概括为了实用主义）。 implicit 语句将导入与声明分开，按照惯例，我用一个简单的continue 语句将声明与可执行语句分开。 continue 被认为是可执行的，因此如果声明出现在第一个 continue 之后，编译器应该会报错。

Answer 2

错误很明显：

integer::qm,a

使a 成为一个整数（标量）。

子程序需要一个数组：

integer, dimension(size(e,1)*2), intent(out) :: A

所以你在这里传递了错误的东西：

call gonewrong(q,e,A) !gonewrong takes q and e as arguments

因为 Fortran 不区分大小写！您将标量整数 a/A 传递给需要数组的子例程。

最好在A 的子例程中使用相同的(:) 并让编译器告诉您您正在出错。

参数是一个数组，所以你必须在那里传递一个数组。

注意：在您的原始代码中，您已将这一行注释掉

! integer, dimension(size(e,1)*2) :: A

当你有两条线时

integer::qm,a
integer, dimension(size(e,1)*2) :: A

程序将无法编译，因为a/A 被声明了两次。

---

对我来说不太清楚的是你的实际意图。你的评论

!gonewrong takes q and e as arguments

真的很混乱。这三件事都是gonewrong 的参数。 q、e 和 A。