f2py处理inout参数未生成正确包装器的问题咨询

问题核心原因 #

1. inout参数的返回逻辑：f2py对intent(inout)的默认处理是原地修改传入对象，而非将其作为返回值返回。这是因为Fortran的inout语义本身就是修改实参，f2py映射时严格遵循这一底层逻辑。
2. 标量被转为数组的原因：Python原生数值类型（如float、int）是不可变的，无法实现原地修改。f2py会自动将这些标量转换为rank-0 numpy数组（即0维数组），这样Fortran端才能修改其内容，保证inout语义生效。

解决方案 #

方案1：正确使用原地修改（推荐） #

不要传入Python原生标量，改用numpy数组承载inout参数：

import numpy as np
import fortran_file

io = np.array(7.5, dtype=np.float32)  # 用numpy数组替代原生float
out1, out2, out3 = fortran_file.func(5, 2.5, False, io)
print(io[0])  # 直接读取修改后的数组值，输出15.0

方案2：让inout参数作为返回值返回 #

如果需要将inout参数作为返回值的一部分，可以在签名文件（.pyf）中给参数同时添加intent(inout, out)：
修改后的sig_example.pyf：

!    -*- f90 -*-
python module fortran_file
    interface
        subroutine func(inputA, inputB, inputC, ioD, outE, outF, outG)
            integer, intent(in) :: inputA
            real, intent(in) :: inputB
            logical, intent(in) :: inputC
            real, intent(inout, out) :: ioD  # 同时指定inout和out属性
            real, intent(out) :: outE, outF, outG
        end subroutine func
    end interface
end python module fortran_file

重新编译后，调用函数即可获取4个返回值：

io = 7.5
io, out1, out2, out3 = fortran_file.func(5, 2.5, False, io)

批量处理（无需手动编辑pyf文件） #

针对大量Fortran文件，有两种高效方式批量指定参数意图：

方式1：命令行参数指定 #

编译时通过--intent-inout、--intent-in、--intent-out参数批量设置：

python.exe -m numpy.f2py -m fortran_file ^
  --intent-in=inputA,inputB,inputC ^
  --intent-inout=ioD ^
  --intent-out=outE,outF,outG ^
  -c example.f

如果需要让inout参数作为返回值，可额外添加--intent-out=ioD（与--intent-inout同时使用）。

方式2：Fortran代码内添加注释指令 #

在Fortran代码中通过!f2py注释直接指定参数意图，f2py编译时会自动识别：

subroutine func(inputA, inputB, inputC, ioD, outE, outF, outG)

integer inputA
!f2py intent(in) :: inputA
real inputB
!f2py intent(in) :: inputB
logical inputC
!f2py intent(in) :: inputC

real ioD
!f2py intent(inout, out) :: ioD  # 如需返回则添加out属性
real outE, outF, outG
!f2py intent(out) :: outE, outF, outG
real localH

if(.not.inputC) then
   localH= ioD
else
   localH= inputB
endif

ioD= ioD + localH
outE= inputA + 10.5
outF= inputA + 5.5
outG= inputA + 1.5
return
end

直接编译即可，无需生成和修改pyf文件：

python.exe -m numpy.f2py -m fortran_file -c example.f

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