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 的高效协同，就不再是难题了。