Python 中使用 FORTRAN 怎么搭

发布网友 发布时间：2022-05-11 01:43

我来回答

共3个回答

热心网友 时间：2022-05-11 03:12

f2py是numpy自带的一个工具,只要安装了numpy,一般都会安装f2py.exe和f2py.py两个文件.

f2py.exe适合在命令行中使用,而f2py.py则一般在Python代码中使用.

废话少说, 先建立一个testfortran.f90的文件如下:

[plain] view plain copy
!SUBROUTINE
SUBROUTINE ADDSUB(A,B,C,D)
IMPLICIT NONE
DOUBLE PRECISION A,B,C,D
!f2py intent(in) :: A,B
!f2py intent(out) :: C,D
C = A + B
D = A - B
print*, "ADDSUB From Fortran!"
print*, "ADD=",C
print*, "SUB=",D
RETURN
END

注意这两行的代码:

[plain] view plain copy
!f2py intent(in) :: A,B
!f2py intent(out) :: C,D

对于Fortran只是注释,但对于f2py却很重要,相当于"签名".

注意签名的注释前面不能有空格!

当然也可以写成如下的形式:

[plain] view plain copy
Cf2py intent(in) :: A,B
Cf2py intent(out) :: C,D
下面开始编译Fortran代码为python模块,打开CMD窗口,输入如下命令:

[plain] view plain copy
f2py -m testfortran -c testfortran.f90
会在当前目录下生成testfortran.pyd的文件.

下面就可以再python中使用这个模块了:

[python] view plain copy
In [1]: import testfortran

In [2]: print testfortran.__doc__
This mole 'testfortran' is auto-generated with f2py (version:2).
Functions:
c,d = addsub(a,b)
.

In [3]: x=testfortran.addsub(4,9)
ADDSUB From Fortran!
ADD= 13.000
SUB= -5.000
In [4]: x
Out[4]: (13.0, -5.0)

In [5]:

热心网友 时间：2022-05-11 04:30

只能是fortran程序编译为可执行文件，python去调用执行。