f2py支持在fortran语言中调用其他Fortran函数或C代码或Python代码

f2py支持在fortran语言中调用其他Fortran函数或C代码或Python代码

分类： Python Numpy_Scipy fortran MinGW_GCC_Boost f2py 2012-03-21 13:20 722人阅读 评论(0) 收藏 举报

fortran python 语言 c subroutine

请看示例F77文件:

[plain] view plain copy

1. C FILE: CALLBACK.F  
2.       SUBROUTINE FOO(FUN,R)  
3.       EXTERNAL FUN  
4.       INTEGER I  
5.       REAL*8 R  
6. Cf2py intent(out) r  
7.       R = 0D0  
8.       DO I=-5,5  
9.          R = R + FUN(I)  
10.       ENDDO  
11.       END  
12. C END OF FILE CALLBACK.F  

在DOS窗口输入:

[plain] view plain copy

1. f2py -c -m callback callback.f  

下面是测试:

[python] view plain copy

1. >>> import callback  
2. >>> print callback.foo.__doc__  
3. foo - Function signature:  
4.   r = foo(fun,[fun_extra_args])  
5. Required arguments:  
6.   fun : call-back function  
7. Optional arguments:  
8.   fun_extra_args := () input tuple  
9. Return objects:  
10.   r : float  
11. Call-back functions:  
12.   def fun(i): return r  
13.   Required arguments:  
14.     i : input int  
15.   Return objects:  
16.     r : float  
17.   
18. >>> def f(i): return i*i  
19. ...   
20. >>> print callback.foo(f)       
21. 110.0  
22. >>> print callback.foo(lambda i:1)  
23. 11.0  

更纤细的说明,请看官方网站:

http://cens.ioc.ee/projects/f2py2e/

 

可以再看下面的例子:

[plain] view plain copy

1.       subroutine calculate(x,n)  
2. cf2py intent(callback) func  
3.       external func  
4. c     The following lines define the signature of func for F2PY:  
5. cf2py real*8 y  
6. cf2py y = func(y)  
7. c  
8. cf2py intent(in,out,copy) x  
9.       integer n,i  
10.       real*8 x(n)  
11.       do i=1,n  
12.          x(i) = func(x(i))  
13.       end do  
14.       end  

在DOS窗口输入:

[python] view plain copy

1. f2py -c -m foo calculate.f  

下面是在Python中测试:

 

[python] view plain copy

1. >>> import foo  
2. >>> foo.calculate(range(5), lambda x: x*x)  
3. array([  0.,   1.,   4.,   9.,  16.])  
4. >>> import math  
5. >>> foo.calculate(range(5), math.exp)  
6. array([  1.        ,   2.71828175,   7.38905621,  20.08553696,  54.59814835])  

以上都是官方的示例