如何访问用 Python 编写的 FORTRAN 函数？

Question

有没有办法在 FORTRAN 中使用 python 函数？我得到了一个包含一些函数的 python 脚本，我需要从 FORTRAN 代码访问这个函数。

我见过“f2py”，它允许从 Python 访问 FORTRAN 子例程，以及 py2exe，它将 python 脚本编译为可执行文件。有“py2f”的东西吗？ Python 脚本可以编译为目标文件吗？然后我可以将它与 FORTRAN 代码链接起来。

例如，将“mypython_func.py”视为包含函数的 Python 脚本，将“mainfortran.f”视为调用 Python 函数的主程序 FORTRAN 代码。我想： 从 'mypython_func.py' 编译为 'mypython_func.o'，从 'mainfortran.f' 编译为 'mainfortran.o' (>> gfortran -c mainfortran.f)，然后链接这些文件 (>> gfortran -c mainfortran.o mypython_func.o -o myexec.exe）。这样的事情可能吗？

感谢您的时间和帮助。

文斯

Answer 1

不要浪费大量时间编译和翻译。这样做。

1. Fortran 第 1 部分程序编写了一个供 Python 执行的内容的文件。写入标准输出。将此称为 F1

2. Python 读取文件，进行 Python 计算，将响应写入 Fortran 文件。将此称为 P。

3. Fortran 第 2 部分程序从标准输入读取内容文件。这些是Python计算的结果。

连接它们

F1 | python p.py | F2

您无需重新编译任何内容。另请注意，所有三个同时运行，这可能会带来相当大的加速。

Python的中间部分应该是这样的。

import sys
import my_python_module
for line in sys.stdin:
    x, y, p, q = map( float, line.split() )
    print ("%8.3f"*6) % ( x, y, z, p, q, my_python_module.some_function( x, y, p, q ) )

函数的简单包装器，以 Fortran 友好的格式读取 stdin 并写入 stdout。

Answer 2

@S.Lott 的解决方案是可行的方法。但要回答你的问题——是的，可以从 Fortran 调用 Python。我首先使用 Cython 的 C API 公开 Python，然后使用 iso_c_binding 模块创建这些 C 函数的 Fortran 接口，最后从 Fortran 调用它们。它非常重量级，在大多数实际问题中不值得（使用管道方法），但它是可行的（它对我有用）。