f2py数组值函数

Question

最新版本的f2py是否支持包含数组值的fortran函数？在一些古代文献中,这不受支持.现在怎么样？

让我们将以下函数保存为func.f95.

function func(x)
    implicit none
    double precision :: x(:),func(size(x))
    integer :: i
    do i=1,size(x)
        func(i) = i*x(i)
    end do
end function

我用它编译 f2py --fcompiler=gnu95 -c -m func func.f95

然后让以下python代码为test_func.py

import func
from numpy import array

x = array(xrange(1,10),dtype='float64')
print 'x=',x

y = func.func(x)
print 'func(x)=',y

输出
python test_func.py是

x= [ 1.  2.  3.  4.  5.  6.  7.  8.  9.]
Segmentation fault

Answer 1

f2py 的机制将 Fortran子例程转换为 python 函数。它不明白如何将 Fortran函数转换为 python 函数。我发现我需要用子例程包装所有 Fortran 函数，或者更好的是，将它们重写为子例程。