运行时多态场景下gfortran内存泄漏问题技术问询

一、代码是否符合Fortran标准？ #

你写的这种工厂函数返回class(...) allocatable多态实例的代码，完全符合Fortran 2008及后续版本的标准。Fortran明确允许函数返回可分配的多态对象，并且要求编译器在对象生命周期结束时自动调用析构逻辑释放资源。你遇到的内存未释放问题，是gfortran 7.5.0/9.4.0/10.3.0版本的编译器实现缺陷——ifort能正常处理也佐证了这一点，不是代码本身违反标准。

二、两种泄漏场景的本质 #

• 场景1：结果用于表达式时未立即释放
  当工厂函数的返回值直接用于表达式（比如作为实参传给其他函数），标准要求这个临时对象在表达式执行完后立刻被销毁。gfortran没触发析构，属于内存泄漏的编译器bug。
• 场景2：赋值给class(...), allocatable变量时未销毁旧实例
  给可分配多态变量赋值时，标准要求编译器先自动释放变量已有的实例（如果存在），再分配新实例并拷贝数据。gfortran没正确完成旧实例的销毁步骤，同样是编译器的问题。

三、兼容旧gfortran的正确多态用法 #

如果必须兼容有bug的gfortran版本，可以用以下两种方式规避：

1. 改用指针类型管理多态对象 #

把工厂函数改成返回指针类型，手动控制内存分配与释放：

function create_shape(shape_type) result(res)
  class(Shape), pointer :: res
  integer, intent(in) :: shape_type

  select case(shape_type)
    case(1)
      allocate(Circle :: res)
    case(2)
      allocate(Square :: res)
  end select
end function create_shape

使用时记得手动释放：

class(Shape), pointer :: my_shape
my_shape => create_shape(1)
! 这里使用对象
deallocate(my_shape)

2. 显式控制临时对象生命周期 #

对于直接用表达式的场景，先把工厂函数的返回值存到可分配变量里，用完后显式释放：

! 原泄漏写法：
call print_area(create_shape(1))

! 规避写法：
class(Shape), allocatable :: temp_shape
temp_shape = create_shape(1)
call print_area(temp_shape)
deallocate(temp_shape)

3. 升级gfortran版本 #

最彻底的解决办法是升级到gfortran 11.x及以上版本——这些后续版本已经修复了多态可分配对象的析构bug。

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