从Fortran回调Python

从Fortran回调Python,python,callback,fortran90,f2py,Python,Callback,Fortran90,F2py,现在我使用f2py从Fortran代码调用Python函数。我试过一个很简单的例子,但没有成功 Fortran90代码: subroutine foo(fun,r) external fun integer ( kind = 4 ) i real ( kind = 8 ) r r=0.0D+00 do i= 1,5 r=r+fun(i) enddo end 使用命令行: f2py-c-m回调callback.f90 Python代码: import callback def f(i)

现在我使用
f2py
从Fortran代码调用Python函数。我试过一个很简单的例子,但没有成功

Fortran90代码:

subroutine foo(fun,r)
external fun
integer ( kind = 4 ) i
real ( kind = 8 ) r
r=0.0D+00
do i= 1,5
    r=r+fun(i)
enddo
end
使用命令行:

f2py-c-m回调callback.f90

Python代码:

import callback

def f(i):
    return i * i
print callback.foo(f)
错误:

Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
TypeError: `Required argument 'r' (pos 2) not found`
回溯(最近一次呼叫最后一次):
文件“”,第1行,在
TypeError:'未找到必需参数'r'(位置2)`

您需要将
r
声明为返回值。。。无论如何,这是一个很好的Fortran 90实践。现在
f2py
假设它是一个输入值

subroutine foo(fun,r)
    external fun
    real ( kind = 8 ), intent(out) :: r
    integer ( kind = 4 )           :: i
    r=0.0D+00
    do i= 1,5
        r=r+fun(i)
    enddo
end
f2py
使用Fortran的intent指令确定传递给函数的内容和返回的内容。

python:

# -*- coding: utf-8 -*-
import MArray

print MArray.__doc__
print MArray.s1dim_array.__doc__
print MArray.s2dim_array.__doc__

print "="*60
print help(MArray)
print "="*60
print help(MArray.s1dim_array)


print "="*60
MArray.s1dim_array([6.,7.,8.])


print "="*60
MArray.s2dim_array([[6.,7.,8.],[1,2,3]])




subroutine S1dim_Array (iN_dim, rArray)
  implicit none  
  integer,  intent(in) :: iN_dim
  real,     intent(in) :: rArray(iN_dim)
  print*, iN_dim
  print*, rArray
  Return
End subroutine


subroutine S2dim_Array (iN_Row, iN_Col, rArray)
  implicit none  
  integer,  intent(in) :: iN_Row, iN_Col
  real,     intent(in) :: rArray(iN_Row, iN_Col)
  integer :: i , j
  print*, iN_Row, iN_Col
  do i = 1 , iN_Row
    write(*,*) (rArray(i,j), j = 1,iN_Col)
  enddo
  Return
End subroutine

我不是python专家,但是这个错误不表明需要两个参数吗?嗨,我已经解决了这个问题。我面临的不是另一个需要的争论。我正在尝试另一个.pyf,使用推荐行自动生成并决定哪个参数应该是输入值,哪个应该是输出值。然后,问题是当我调用函数时,参数的类型不同。所以我把r定为一个整数,问题就解决了。对不起,我不明白这个意思。问题是,我试图首先从python调用fortran代码,然后从fortran调用python函数。fortran中应该有一个外部部分