如何用F2PY包装依赖Fortran模块的有限差分子程序生成fdq.so？

问题场景 #

我有一套用于生成FDq-matrices的遗留Fortran代码，需要将其中的finiteDifferences.f90子程序调用到Python中，该子程序依赖5个Fortran模块：Lagrange_Interpolation_Method.f、Linear_Algebra_Lapack.f、Nonlinear_Algebra.f、Piecewise_Polynomial_Interpolation.f、Optimum_Grid_Methods.f。尝试用F2PY编译生成可导入Python的fdq.so文件时，出现编译错误：

/var/folders/bj/4qqzkp4560d75dw4wr7tpmr00000gn/T/tmp39eelon3/src.macosx-13-arm64-3.11/fdqmodule.c:578:23: error: expected expression nodes_Dims[0]=1 + ; ^

解决方法 #

1. 排查Fortran代码的数组维度语法错误 #

错误提示显示F2PY解析时遇到了不完整的表达式，核心原因是某个与nodes相关的数组维度定义存在语法缺陷：

• 检查所有依赖模块和finiteDifferences.f90中的数组维度声明，确认没有类似1 +这种缺少右侧表达式的写法（正确应为1 + n这类完整表达式）。
• 修正所有语法不完整的维度定义，确保每个数组维度都有合法的计算表达式。

2. 一次性编译所有依赖文件 #

F2PY处理多模块依赖时，必须将所有相关源文件一起传入编译命令，不能仅编译单个子程序。正确的编译命令格式如下：

f2py -c -m fdq finiteDifferences.f90 Lagrange_Interpolation_Method.f Linear_Algebra_Lapack.f Nonlinear_Algebra.f Piecewise_Polynomial_Interpolation.f Optimum_Grid_Methods.f

参数说明：-c表示编译生成扩展模块，-m fdq指定生成的Python模块名为fdq，后续依次列出所有需要编译的Fortran源文件（包括依赖模块和目标子程序）。

3. 确保Fortran语法与F2PY兼容并链接依赖库 #

• 检查所有模块是否使用F2PY支持的标准Fortran语法，避免非标准或过时的特性。
• 如果代码中调用了LAPACK接口，编译时需明确链接LAPACK库，可在命令中添加--f90flags="-llapack"（根据系统环境调整库名）：

f2py -c -m fdq --f90flags="-llapack" finiteDifferences.f90 Lagrange_Interpolation_Method.f Linear_Algebra_Lapack.f Nonlinear_Algebra.f Piecewise_Polynomial_Interpolation.f Optimum_Grid_Methods.f

4. 分步编译模块与子程序 #

若一次性编译仍有问题，可尝试拆分编译步骤：

1. 先将所有依赖模块编译为目标文件：

gfortran -c Lagrange_Interpolation_Method.f Linear_Algebra_Lapack.f Nonlinear_Algebra.f Piecewise_Polynomial_Interpolation.f Optimum_Grid_Methods.f

2. 再用F2PY编译目标子程序并链接已生成的目标文件：

f2py -c -m fdq finiteDifferences.f90 Lagrange_Interpolation_Method.o Linear_Algebra_Lapack.o Nonlinear_Algebra.o Piecewise_Polynomial_Interpolation.o Optimum_Grid_Methods.o

内容的提问来源于stack exchange，提问作者vpezlar