基于ISO_C_BINDING，Fortran调用C代码的结构体映射问题咨询

给定C语言定义的嵌套结构体：

struct Vector
{
  int size;
  double data[3];
};

struct VerticesVect
{
  int size;
  struct Vector* data;
};

需要在Fortran中通过ISO_C_BINDING读取该结构体（仅读不修改），但尝试声明VerticesVect时遇到错误：

• 使用type(Vector), dimension(:) :: data时，报错“data必须具有显式形状”
• 使用type(Vector), pointer :: data时，报错“组件'data'不能具有POINTER属性，因为它是BIND(C)派生类型'verticesvect'的成员”

由于后续嵌套结构体尺寸过大，无法指定静态数组覆盖所有情况，需要找到编译器可接受且映射正确的声明方式。

1. 正确的结构体声明：用c_ptr对应C指针 #

在绑定到C的Fortran派生类型中，不能直接使用Fortran的pointer或allocatable组件（这类组件带有Fortran特有的元数据，和C指针内存布局不兼容），必须用ISO_C_BINDING提供的c_ptr类型来对应C的指针。

正确的声明代码：

use, intrinsic :: ISO_C_BINDING

! 正确映射C的struct Vector
type, bind(c) :: Vector
    integer(kind=C_INT32_T) :: size
    real(kind=C_DOUBLE) :: data(3)
end type Vector

! 正确映射C的struct VerticesVect
type, bind(c) :: VerticesVect
    integer(kind=C_INT) :: size
    type(c_ptr) :: data  ! 对应C的struct Vector*
end type VerticesVect

2. 转换C指针为Fortran可访问的指针数组 #

读取数据时，通过c_f_pointer函数将c_ptr类型转换为Fortran的指针数组，这样就能像访问普通Fortran数组一样操作数据：

subroutine process_vertices(vertices)
    use, intrinsic :: ISO_C_BINDING
    type(VerticesVect), intent(in) :: vertices
    type(Vector), pointer :: vec_array(:)
    integer :: i

    ! 将C指针转换为Fortran指针数组，指定数组长度为vertices%size
    call c_f_pointer(vertices%data, vec_array, [vertices%size])

    ! 示例：遍历并打印数据
    do i = 1, vertices%size
        print *, "Vector ", i, " size: ", vec_array(i)%size
        print *, "Data values: ", vec_array(i)%data
    end do

    ! 无需手动释放内存，原内存由C代码管理
end subroutine process_vertices

关键说明 #

• 不能移除bind(c)：带bind(c)的派生类型会严格遵循C的内存布局（成员顺序、字节填充规则），如果移除，Fortran编译器可能会调整内存布局，导致和C结构体不匹配，引发内存访问错误。
• 必须使用c_ptr：这是ISO_C_BINDING标准中定义的、唯一能和C指针直接对应的Fortran类型，是跨语言交互的规范做法。

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