如何利用c_ptr实现支持不同类型参数的Fortran子程序

Fortran 语法强制要求所有局部变量的声明必须放在子程序的声明段（即第一条可执行语句之前），在if/select case等执行分支内声明变量、用block块声明跨块使用的变量都是不符合规范的。可以通过以下两种合规方式实现需求：

• 方案1：预声明所有类型指针，分支内完成指针绑定
  这是对现有代码改动最小的实现，在子程序开头的声明区提前定义所有可能传入的类型指针，进入分支后仅对匹配类型标记的指针执行c_f_pointer绑定，重复的公共逻辑抽离为内部子程序复用即可，参考实现：

subroutine myfoo(i, params) 
  use, intrinsic :: iso_c_binding
  implicit none
  integer, intent(in) :: i
  type(c_ptr), intent(in) :: params
  ! 声明段提前定义所有支持的类型指针
  type(type1_t), pointer :: pars1
  type(type2_t), pointer :: pars2

  select case(i)
  case(1)
    call c_f_pointer(params, pars1)
    ! 编写type1_t类型的专属处理逻辑
    ! 调用公共逻辑时传入对应成员
    call shared_process(pars1%a, pars1%b, pars1%c)
  case(2)
    call c_f_pointer(params, pars2)
    ! 编写type2_t类型的专属处理逻辑
    call shared_process(pars2%a, pars2%e, pars2%b, pars2%c, pars2%d)
  case default
    error stop "传入的类型标记无效"
  end select

contains
  ! 抽离两类参数的公共处理逻辑，避免重复代码
  subroutine shared_process(base_a, val1, c, val2, val3)
    real, intent(in) :: base_a
    integer, intent(in) :: c
    real, optional, intent(in) :: val1, val2
    integer, optional, intent(in) :: val3
    ! 此处编写两类参数共用的业务代码
  end subroutine shared_process
end subroutine myfoo

注意：不要尝试通过内存共用、偏移量计算等方式跨类型直接访问内存，这类写法违反Fortran的类型别名规则，会触发未定义行为，编译器开启优化后大概率出现逻辑错误。

• 方案2：使用无限多态指针实现类型安全分发（推荐）
  如果不需要严格兼容C侧的回调接口，更推荐用Fortran原生的无限多态特性替代手动维护类型标记+c_ptr的写法，编译器会自动完成类型检查，避免标记传错导致的内存越界问题。如果必须保留c_ptr接口，只需要多做一次指针转换即可：

subroutine myfoo(i, params)
  use, intrinsic :: iso_c_binding
  implicit none
  integer, intent(in) :: i
  type(c_ptr), intent(in) :: params
  class(*), pointer :: pars ! 无限多态指针，可以指向任意类型

  call c_f_pointer(params, pars)
  select type(pars)
  type is(type1_t)
    ! 分支内可直接访问type1_t的所有成员
    print *, pars%a, pars%b, pars%c
  type is(type2_t)
    ! 分支内可直接访问type2_t的所有成员
    print *, pars%a, pars%e, pars%b, pars%c, pars%d
  class default
    error stop "不支持的参数类型"
  end select
end subroutine myfoo

这种写法不需要手动维护类型标记和实际类型的对应关系，公共逻辑可以直接在select type的对应分支里调用，代码可维护性更高。

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