CentOS中Fortran与Python如何协同工作

在 CentOS 系统里，让 Fortran 和 Python 协同工作，其实有很多成熟的套路可走。下面这几种方法，基本上涵盖了主流的实践路径。 首先，最省心的方案是借助现成的接口工具。比如 f2py，这是 NumPy 社区提供的利器，它能把 Fortran 代码直接“翻译”成 Python 模块，之后你就可以像调用普通 Python 函数一样调用 Fortran 的子程序。另一个工具叫 SWIG，它原本是为 C/C++ 设计的，但 Fortran 代码很容易先编译成 C 接口，然后通过 SWIG 生成 Python 包装，同样很顺手。 如果你不想引入额外工具，直接调用编译好的共享库也是一种经典做法。把 Fortran 代码编译成 `.so` 文件，然后在 Python 里用 `ctypes` 或者 `cdll` 库加载它，按约定传入参数即可。这种方式虽然需要手动处理数据类型对齐，但灵活度很高。 还有些场景更简单——直接用系统调用。在 Python 脚本里用 `subprocess` 模块，调起 gfortran 编译 Fortran 源码，再运行生成的可执行文件。适合那些计算逻辑独立、不需要频繁交互的场景。 下面用一个典型的 f2py 例子来说明流程。假设你有一段 Fortran 子程序，专门做加法运算： **Fortran 代码（example.f90）** ```fortran 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 ``` 然后在 Python 里这样调用： ```python from example import add import numpy as np a = np.array([1.0, 2.0, 3.0], dtype=np.float64) b = np.array([4.0, 5.0, 6.0], dtype=np.float64) c = np.zeros_like(a) add(a, b, c) print(c) # 输出: [5. 7. 9.] ``` 编译只需要一行命令： ```bash f2py -c example.f90 -m example ``` 执行完后会生成 `example.so` 和 `example.pyf` 接口文件。 有几个细节值得留意： - Fortran 代码里的数据类型和调用约定必须与 C 兼容，尤其是 `bind(c)` 声明不能少。 - 编译时可能需要额外指定编译器标志，确保生成的共享库与 Python 解释器使用的 ABI 一致。 - 如果 Fortran 代码依赖了第三方库（比如 LAPACK、BLAS），记得要确保这些库在系统库路径中，或者在编译时用 `-L` 和 `-l` 明确指定。 掌握了这些方法，在 CentOS 上让 Fortran 和 Python 高效协同工作，基本就不是什么难事了。