Ubuntu中Fortran与Python如何交互

在Ubuntu系统中实现Fortran与Python交互的三种高效方法

免费影视、动漫、音乐、游戏、小说资源长期稳定更新！ 👉 点此立即查看 👈

1. 使用f2py工具进行集成

对于希望在Python中调用高性能Fortran计算模块的开发者，f2py是最为成熟和广泛采用的解决方案。作为NumPy生态系统的重要组成部分，它能够将Fortran子程序自动封装为可直接导入的Python扩展模块，实现两种语言的无缝衔接。

安装f2py环境

首先需要配置基础开发环境，确保系统中已安装Python科学计算栈和Fortran编译器。在Ubuntu终端中执行以下命令：

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

编写兼容性Fortran代码

创建Fortran源文件时，必须使用iso_c_binding模块确保与C语言的ABI兼容，这是跨语言调用的技术基础。以下是一个标准示例：

! example.f90
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可调用模块

通过f2py命令行工具将Fortran源代码编译为Python扩展：

f2py -c example.f90 -m example

成功执行后，系统将生成两个关键文件：example.so动态链接库和example.py接口封装文件，完成Fortran功能的Python化封装。

在Python环境中调用Fortran模块

现在可以在Python脚本中像使用原生模块一样调用Fortran函数：

import example

a = 1.0
b = 2.0
c = example.add(a, b)
print(f"The result is {c}")

这种方法完美结合了Fortran的数值计算性能与Python的灵活性和丰富的生态系统。

2. 通过ctypes进行底层调用

对于需要更精细控制内存管理和调用过程的开发者，Python标准库中的ctypes模块提供了直接调用C兼容共享库的能力。这种方法的核心是将Fortran代码编译为标准共享库，再通过Python进行动态加载。

编译C兼容的Fortran共享库

首先创建C语言头文件，明确定义函数接口规范：

// example.h
#ifndef EXAMPLE_H
#define EXAMPLE_H

#ifdef __cplusplus
extern "C" {
#endif

void add_(double *a, double *b, double *c);

#ifdef __cplusplus
}
#endif

#endif // EXAMPLE_H

调整Fortran源代码，确保绑定名称与C接口完全匹配：

! example.f90
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

执行两阶段编译过程，生成标准共享库文件：

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

使用ctypes实现Python调用

在Python中通过ctypes精确声明函数签名并调用共享库：

import ctypes

# 加载共享库
libexample = ctypes.CDLL('./libexample.so')

# 定义函数的参数和返回类型
libexample.add_.argtypes = [ctypes.c_double, ctypes.c_double, ctypes.POINTER(ctypes.c_double)]
libexample.add_.restype = None

# 调用函数
a = ctypes.c_double(1.0)
b = ctypes.c_double(2.0)
c = ctypes.c_double()
libexample.add_(ctypes.byref(a), ctypes.byref(b), ctypes.byref(c))

print(f"The result is {c.value}")

这种方法虽然需要更多手动配置，但提供了对内存布局和调用约定的完全控制，适合高性能计算和系统级编程场景。

3. 利用pyfortran进行代码转换

pyfortran提供了一种创新性的思路——将Fortran源代码直接转换为等效的Python实现。这种方法特别适合希望快速迁移遗留Fortran代码或理解算法逻辑的开发者。

安装pyfortran转换工具

通过Python包管理器直接安装：

pip install pyfortran

准备Fortran源代码

使用标准Fortran语法编写功能模块：

! example.f90
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

执行Fortran到Python的转换

运行转换命令，生成可直接执行的Python模块：

pyfortran example.f90 -o example.py

使用转换后的Python代码

导入生成的模块并调用转换后的函数：

import example

a = 1.0
b = 2.0
c = example.add(a, b)
print(f"The result is {c}")

需要注意的是，pyfortran对现代Fortran语法的支持仍在不断完善中，对于复杂的数值计算程序，建议进行充分的测试验证。

综上所述，在Ubuntu平台上实现Fortran与Python的交互主要有三种技术路径：成熟稳定的f2py方案、底层灵活的ctypes方案以及创新性的pyfortran转换方案。开发者应根据项目的具体需求、性能要求、维护成本和技术栈兼容性等因素，选择最适合的集成策略。对于科学计算和工程应用，合理利用Fortran的计算性能与Python的生态系统，能够显著提升开发效率和程序性能。