从Fortran访问C ++结构数组？

Question

在C ++中，我分配了一个S数组。在Fortran中，我想访问这个数组的元素。 我怎样才能做到这一点？

C ++：

struct S {double a; double b;};
S *arrayOfS;
arrayOfS = (S *)new S[123]; // allocate

Fortran 2003：

USE ISO_C_Binding
TYPE, BIND(C) :: SFortran
REAL(c_double) :: a,b
END TYPE SFortran

S和SFortran现在应该是可互操作的，但我还需要有一种方法来访问在C ++中声明的数组元素。 我希望Fortran中的SC（5）％a对应于C ++中的arrayOfS [4] .a。 如何为具有此访问权限的Fortran阵列SC声明和设置正确的值？

Answer 1

你可以：

1）将C ++数组传递给Fortran BIND（C）过程，该过程采用适当的数组参数。

SUBROUTINE proc(array) BIND(C, NAME='proc')
  ...
  TYPE(SFortran) :: array(*)

使用这种方法，您可能还希望传递数组的大小并使数组参数显式化。

b）将数组指针作为C ++端的外部“C”全局，然后通过Fortran模块变量与BIND（C）进行互操作。

MODULE some_module
  USE, INTRINSIC :: ISO_C_BINDING, ONLY: C_PTR, C_F_POINTER
  ...
  TYPE(C_PTR), BIND(C, NAME='arrayOfS') :: array_ptr
  ...
  ! In a procedure in the module...
  TYPE(SFortran), POINTER :: array(:)
  CALL C_F_POINTER(array_ptr, array, [123])

同样，它可能适合单独提供阵列的大小，以避免在C_F_POINTER引用中对其进行硬编码。

哪种方法最好取决于您的要求！