Debian如何集成Fortran与Python

在Debian系统中集成Fortran和Python：三种实战路径解析

将历史悠久的Fortran与灵活现代的Python结合，是科学计算领域一个经典且高效的策略。在Debian系统上，有几种成熟的技术路线可以实现这一目标，每种方法各有侧重，适用于不同的场景和需求。

方法一：使用f2py——最直接的“翻译官”

对于Python生态的用户而言，f2py通常是集成Fortran代码的首选工具。它本质上是NumPy套件中的一个“翻译官”，能够将Fortran子程序无缝封装成可以直接导入的Python模块。整个过程相当直观。

1. 第一步，自然是搭建好基础环境。确保系统已安装Fortran编译器和NumPy：

   sudo apt-get update
   sudo apt-get install gfortran
   pip install numpy

2. 接着，准备好你的Fortran源代码。这里以一个简单的加法子程序为例，保存为example.f90。注意，为了与C/Python交互，我们使用了iso_c_binding和bind(c)属性：

   subroutine add(a, b, c) bind(c, name="add")
       use, intrinsic :: iso_c_binding
       real(c_double), intent(in) :: a, b
       real(c_double), intent(out) :: c
       c = a + b
   end subroutine add

3. 关键一步来了。使用f2py命令进行编译和封装：

   f2py -c example.f90 -m example

   这个命令会生成一个名为example的Python扩展模块。

4. 最后，在Python中调用它，就像调用普通函数一样自然：

   import example
   result = example.add(1.0, 2.0)
   print(result)  # 输出 3.0

看，整个过程是不是行云流水？f2py自动处理了类型转换和接口生成，让Fortran代码几乎“零成本”地融入Python世界。

方法二：使用cython——更灵活的“中间层”

如果你需要更精细的控制，或者项目本身就在使用Cython，那么通过Cython来包装Fortran代码是一个极具灵活性的方案。它相当于在Fortran和Python之间搭建了一个可定制的“中间层”。

1. 首先，安装必要的工具链：

   sudo apt-get update
   sudo apt-get install cython gfortran

2. Fortran源代码部分与之前类似，这里不再赘述。

3. 接下来是核心：编写Cython接口文件（example.pyx）。这个文件负责声明Fortran函数的C接口，并定义一个Python友好的包装函数：

   cdef extern from "example.f90":
       void add_(double *a, double *b, double *c)
   
   def add(double a, double b):
       cdef double c
       add_(&a, &b, &c)
       return c

4. 然后，创建一个setup.py文件来指导构建过程。注意这里需要显式链接Fortran运行时库：

   from setuptools import setup
   from Cython.Build import cythonize
   
   setup(
       ext_modules=cythonize("example.pyx"),
       extra_link_args=['-lgfortran']
   )

5. 执行构建命令，生成模块：

   python setup.py build_ext --inplace

6. 使用方式与第一种方法完全一致：

   import example
   result = example.add(1.0, 2.0)
   print(result)  # 输出 3.0

这种方法虽然步骤稍多，但赋予了开发者对内存布局和调用细节的完全掌控权，适合处理复杂的、高性能的交互场景。

方法三：使用iso_c_binding与ctypes——最底层的“握手”

当你追求极致的简洁，或者不想引入额外的工具依赖时，直接利用Fortran 2003标准引入的iso_c_binding模块与Python的ctypes库进行交互，是最轻量、最底层的方案。这相当于让两种语言在ABI（应用二进制接口）层面直接“握手”。

1. 编写Fortran代码时，确保使用bind(c)来创建标准的C接口。可以将其封装在一个模块中：

   module example_module
       use iso_c_binding
       implicit none
   
       interface
           subroutine add(a, b, c) bind(c, name="add")
               import :: c_double
               real(c_double), intent(in) :: a, b
               real(c_double), intent(out) :: c
           end subroutine add
       end interface
   
   end module example_module

2. 将其编译为共享库（.so文件）：

   gfortran -fPIC -c example.f90 -o example.o
   gfortran -shared example.o -o libexample.so

3. 在Python端，使用ctypes来加载这个共享库并调用函数。这一步需要精确匹配C接口的定义：

   import ctypes
   
   # 加载共享库
   lib = ctypes.CDLL('./libexample.so')
   
   # 定义函数原型：参数类型和返回类型
   lib.add.argtypes = [ctypes.c_double, ctypes.c_double, ctypes.POINTER(ctypes.c_double)]
   lib.add.restype = None
   
   # 准备参数并调用
   a = ctypes.c_double(1.0)
   b = ctypes.c_double(2.0)
   c = ctypes.c_double()  # 用于接收输出
   
   lib.add(a, b, ctypes.byref(c))
   print(c.value)  # 输出 3.0

这种方法绕过了所有中间工具，直接进行二进制交互，虽然代码稍显繁琐，但依赖最少，理解起来也最直接，是理解跨语言调用本质的绝佳途径。

总结来说，三种方法各有千秋：f2py胜在便捷，cython强在灵活，而iso_c_binding+ctypes则赢在轻量和透明。根据你的项目对便捷性、控制力和依赖管理的不同要求，选择最适合的那把“钥匙”即可。