f2py封装Fortran子程序非纯特性及状态提取技术问询

我来帮你拆解这个f2py和Fortran状态保留的问题——这其实是Fortran本身的特性结合f2py封装机制导致的常见情况，下面分点给你说明：

1. 状态如何以及在何处被保留？ #

这些跨调用的状态，本质上存储在Fortran代码里的持久化变量中，f2py封装后，这些变量会依附于Python的扩展模块内存空间，属于进程级的全局状态：

• 如果Fortran代码用了模块变量（比如module里定义的变量），这些变量的生命周期和整个模块绑定，只要Python的扩展模块没被卸载，变量值就会一直保留；
• 另一种情况是COMMON块（Fortran的老式全局变量机制），COMMON块里的变量会被分配到固定的内存区域，同样在模块加载期间持续存在；
• 还有带**save属性**的局部变量，即使是函数内定义的变量，加了save后也会在函数调用之间保留值，不会在函数退出时被销毁。

这些变量的内存由f2py生成的扩展模块管理，每次调用封装后的函数，都会直接访问这些内存里的现有值。

2. 现象产生的原因是什么？ #

核心原因是Fortran本身支持非纯函数的状态持久化，而f2py在封装时完整保留了Fortran代码的变量生命周期逻辑。

举个简单的例子，假设你的Fortran代码是这样的：

module accumulator
  real :: total = 0.0
contains
  function add_val(val) result(res)
    real, intent(in) :: val
    real :: res
    total = total + val
    res = total
  end function add_val
end module accumulator

用f2py封装后，每次在Python里调用add_val(1.0)，返回值会依次是1.0、2.0、3.0……因为模块里的total变量在每次调用后都会保留更新后的值，下一次调用会基于这个值计算。

简单说，就是这个Fortran函数没有做到“输入相同则输出相同”，而是依赖了函数外部的、会被修改的持久化变量，所以变成了非纯函数。

3. 是否可以提取函数调用前已存在的状态？ #

可以，但需要你手动给Fortran代码添加状态访问接口，再用f2py重新封装。

比如上面的累加器例子，你可以在模块里新增一个读取total的函数：

function get_total() result(current_val)
  real :: current_val
  current_val = total
end function get_total

重新编译后，在Python里就能调用get_total()来获取当前的状态值。如果是COMMON块里的变量，同样可以写一个专门的函数来读取块内的变量；如果是带save属性的局部变量，也需要添加对应的访问函数来暴露这个值——毕竟f2py默认不会把这些内部变量直接暴露给Python，必须显式提供接口才能访问。

注意：如果你的代码是多线程环境，这种全局状态会有线程安全问题，多个线程同时调用函数可能会导致状态混乱，这时候可能需要额外的同步机制，但这是另一个场景的问题了。

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