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的chanism將Fortran 子程序轉換爲python函數。它不知道如何將Fortran 函數轉換爲python函數。我發現我需要用子程序包裝所有Fortran函數，或者甚至更好，將它們重寫爲子程序。