在 FORTRAN-77 中传递二维数组的子数组

Question

我有二维数组

real triangle(0:2, 0:1)

其中“三角形”是一个向量数组（一维数组） 我也有子程序

subroutine vecSub(lhs, rhs, result)
real lhs(0:1), rhs(0:1), result(0:1)
    result(0) = lhs(0) - rhs(0)
    result(1) = lhs(1) - rhs(1)
    return
end

有没有办法将“三角形”变量中的一个向量传递给这个子程序？ Fortran-90 可以做到这一点：triangle(0, :) 它给出了第一个三角形数组，但我只能使用 FORTRAN-77，所以这不行，有什么建议吗？

Answer 1

@Javier Martin 写道“不使用阵列的当前布局”，但错过了提出替代方案的机会。

如果您将变量声明如下：

real triangle(0:1, 0:2)

反转边界的顺序，然后您可以将triangle(0,0)、triangle(0,1) 或triangle(0,2) 传递给子例程，并获得您想要的行为，这要归功于名为“序列关联”的 Fortran 功能。当您将单个数组元素传递给作为数组的虚拟参数时，您将按数组元素顺序隐式传递该数组和后面的元素。这是唯一允许违反正常 Fortran 形状匹配规则的行为，并且是 FORTRAN 77 的一部分。

Answer 2

否，与阵列的当前布局无关，原因有二：

• Fortran 使用数组元素顺序，其中最左边的维度是连续的。也就是说，在大小为(n,m,l) 的数组中，元素之间的距离（步幅）为(1,n,m)，以数组元素为单位（即，不是字节）测量。
• F77 不包括假定形状数组a(:)，它们通常通过传递一个小型描述符结构来实现，该结构传达诸如步幅或元素数量等细节。相反，您只能使用假定长度数组a(*)，它们通常是指向第一个元素的指针，有点像 C 数组。您必须将长度作为单独的参数传递，并且数组元素必须是连续的

这就是为什么您可以“将子数组”“传递给 F77 子例程”的原因，只要该子数组是例如矩阵列：其中的元素是连续的。

一种可能的解决方案（许多当前的 Fortran 编译器实现的解决方案）是，当您尝试将不连续的子数组传递给未知接受它们的函数时，它们会复制该数组，甚至将其写入如果需要，返回内存。这相当于：

! Actual array
integer m(3,5)
integer dummy(5)
dummy = m(2,:)
call myF77sub(dummy, 5)
m(2,:) = dummy

但是，正如其他人所说，您应该尽量不要直接调用 F77 函数，而应使它们适应或至少将它们包装在更新的 Fortran 接口中。然后，您可以在包装器中包含上述代码，并从现代 Fortran 例程中“正常”调用该包装器。然后，您最终可能会在不影响客户端代码的情况下改写现代 Fortran 中的实际实现。