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

时间:2013-08-09 20:21:12

标签: c++ fortran dynamic-arrays language-interoperability

在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声明和设置正确的值?

1 个答案:

答案 0 :(得分:3)

你可以:

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引用中对其进行硬编码。

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