要在Python中调用Fortran动态库,可以使用ctypes
库。下面是一个简单的示例:
首先,编写Fortran代码并将其编译为动态库。假设我们有一个名为mylib.f90
的Fortran源文件,其中包含一个名为myfunc
的函数:
! mylib.f90
subroutine myfunc(a, b, c)
implicit none
real*8, intent(in) :: a, b
real*8, intent(out) :: c
c = a + b
end subroutine myfunc
接下来,使用Fortran编译器将源代码编译为动态库。例如,使用gfortran
编译器可以执行以下命令:
gfortran -shared -o mylib.so mylib.f90
这将生成一个名为mylib.so
的动态库文件。
然后,在Python中使用ctypes
库加载并调用这个动态库。例如:
import ctypes
# 加载动态库
mylib = ctypes.CDLL('./mylib.so')
# 声明函数的参数类型
mylib.myfunc.argtypes = (ctypes.c_double, ctypes.c_double, ctypes.POINTER(ctypes.c_double))
mylib.myfunc.restype = None
# 调用函数
a = 3.0
b = 4.0
c = ctypes.c_double()
mylib.myfunc(a, b, ctypes.byref(c))
print(c.value) # 输出结果:7.0
在上述示例中,我们首先使用ctypes.CDLL
函数加载动态库。然后,通过设置argtypes
和restype
属性来指定函数的参数类型和返回类型。最后,使用ctypes.byref
将变量传递给Fortran函数,并使用value
属性获取返回值。
请注意,上述示例假设动态库文件mylib.so
位于当前目录下。如果文件位于不同的位置,需要提供正确的路径。
希望对你有帮助!