将数组传递给具有假定形状的函数

Question

我正在使用 Python 调用由 f2py 绑定的 Fortran 函数。当我显式标注数组时，我的 sum 调用返回所需的结果，但是当我使用假定形状时，它返回 0

! foo.f95
function sum_test(arr)
  IMPLICIT NONE
  integer(8), dimension(:), intent(in) :: arr
  integer(8) :: sum_test
  sum_test = sum(arr)
end function sum_test

Python 方面：

import foo
foo.sum_test([1,2,3])
0L

如果我在 Fortran 中明确标注：

! foo.f95
function sum_test(arr)
  IMPLICIT NONE
  integer(8), dimension(3), intent(in) :: arr
  integer(8) :: sum_test
  sum_test = sum(arr)
end function sum_test

Python 方面：

import foo
foo.sum_test([1,2,3])
6L

请注意，如果我像这样在假定形状的版本上打印出我的值：

write(*,*) arr

我可以看到数组中的值。

我显然在这里遗漏了一些关键部分！

Answer 1

假定形状的数组要求调用者传递有关数组边界的信息以及数据地址。这样做的机制是依赖于实现的，并不是所有的实现都记录了它们的方法。您的 Python 代码只是传递数据地址，但 Fortran 需要（通常）一个“描述符”数据结构。 （Fortran 2015 指定了执行此操作的标准方法，作为“进一步的 C 互操作性”的一部分。）

我建议改为 Fortran 所谓的“可调整数组”，您可以在其中将范围作为单独的参数传递。