Windows10下用f2py编译Fortran90遇__imp_PyFloat_Type未定义引用错误

嘿，我之前也踩过类似的坑——gfortran单独编译Fortran代码完全没问题，但用f2py做跨语言绑定时就卡壳，尤其在Win10+Python(x,y)的环境里，大概率是调用方式或环境配置的小细节出了问题，咱们一步步来排查：

1. 先确认f2py的调用命令是否正确 #

很多新手容易忽略-m参数（指定生成的Python模块名），或者没明确指定Fortran标准。试试用这个标准命令：

f2py -c -m ljlib ljlib3.f90 --f90flags="-std=f90"

这里-m ljlib会让生成的Python模块名为ljlib，之后你就能在Python里通过import ljlib调用子程序；--f90flags="-std=f90"强制gfortran用Fortran 90标准编译，避免和新版本编译器的默认标准产生兼容性冲突。

2. 检查环境变量与编译器指定 #

Win10的环境变量优先级很容易踩坑，确保Mingw-w64的bin目录（比如C:\mingw-w64\mingw64\bin）放在系统PATH的最顶端，避免和其他可能存在的编译器（比如MSVC）冲突。

另外，有时候f2py会默认找错编译器，你可以手动指定使用gnu的Fortran编译器：

f2py --fcompiler=gnu95 -c -m ljlib ljlib3.f90

3. 确保Fortran代码的接口符合f2py要求 #

虽然gfortran能正常编译，但f2py对参数的传递方向（intent）很敏感。如果你的子程序没有明确声明参数的intent(in)/intent(out)，或者是单独的子程序（不在模块内），可以加上f2py专属注释来提示：
比如你的子程序可以改成这样：

subroutine lj_calc(r, energy)
    real*8 :: r, energy
    !f2py intent(in) r
    !f2py intent(out) energy
    energy = 4.0d0 * ((1.0d0/r)**12 - (1.0d0/r)**6)
end subroutine lj_calc

这样f2py能更精准地生成Python接口，避免绑定过程中出现参数解析错误。

4. 分步编译定位问题 #

如果还是报错，试试拆分编译步骤，先生成接口文件再编译，这样能快速定位是接口解析问题还是链接问题：

1. 生成接口文件：

f2py -h ljlib.pyf ljlib3.f90

打开生成的ljlib.pyf文件，检查里面的子程序声明是否和你的Fortran代码一致，有没有遗漏或错误。
2. 用接口文件编译：

f2py -c ljlib.pyf ljlib3.f90

如果这一步出错，就能更精准地看到是接口定义的问题，还是编译链接环节的问题。

要是能提供具体的报错信息（比如是链接错误还是语法解析错误），还能更快锁定问题根源~

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