检查 Fortran 调用 C 函数中的 ISO_C_BINDING 用法

Check for ISO_C_BINDING usage in Fortran calling C function

我在ISO_C_BINDING上问了几个问题,我做了一个例程。由于我在处理 Fortran POINTER 时总是犯了很多错误,所以我想确保没有错误或一些不推荐的奇怪点。

(Fortran 部分)

program test
  use iso_c_binding
  implicit none
  interface
    subroutine get_value_array(in, num) bind(C, name='get_value_Array')
      use iso_c_binding
      implicit none
      type(C_PTR),           intent(inout)  :: in
      integer(C_INT), value, intent(in)     :: num
    end subroutine
  end interface

  real(C_DOUBLE), allocatable, target  :: array(:)
  real(C_DOUBLE),              pointer :: array_fptr(:)
  type(C_PTR)                          :: array_cptr

  integer :: array_len

  allocate(array(12))
  array_len = size(array,1)

  array_cptr = C_LOC(array)

  call get_value_array (array_cptr, array_len)

  call C_F_POINTER(array_cptr, array_fptr, [array_len])

  print *, 'array_fptr'
  print *, array_fptr

  print *, 'array'
  print *, array

end program

(C部分)

void get_value_Array(double **in, int num) {
  int i;
  for (i = 0; i < num; i++) {
    (*in)[i] = i+1;
  } 
} 

(输出)

 array_fptr
   1.00000000000000        2.00000000000000        3.00000000000000
   4.00000000000000        5.00000000000000        6.00000000000000
   7.00000000000000        8.00000000000000        9.00000000000000
   10.0000000000000        11.0000000000000        12.0000000000000
 array
   1.00000000000000        2.00000000000000        3.00000000000000 
   4.00000000000000        5.00000000000000        6.00000000000000 
   7.00000000000000        8.00000000000000        9.00000000000000 
   10.0000000000000        11.0000000000000        12.0000000000000 

为了确定,流程是

  1. array_cptr 与具有 C_LOC 功能的 array 关联。 (也许 'associate' 在这里不是一个合适的术语,因为它是 Fortran 指针相关的术语。我应该怎么称呼它?)
  2. array 在我调用 get_value_array 时被修改,因为子例程(或 C 函数)修改了 array_cptr 指向的位置。
  3. 当我调用 C_F_POINTER 时,
  4. array_fptr 通过 array_cptr 关联到 array

最后 array_cptrarray_fptr 是彼此不同的对象 (?),同时指向同一个目标 array,对吗?

问题中的陈述含糊不清。准确地说。

  1. 用Fortran语句array_cptr = C_LOC(array)array_cptr被赋值为array.

    的C地址
  2. 调用get_value_array时,array的值最终被C函数修改。 Fortran 代码还允许子例程(或 C 函数)修改存储在 array_cptr 中的 C 地址 - 伪参数声明为 INTENT(INOUT) 并且没有 VALUE 属性,并且与C 函数通过引用获取指向数组的指针(即 C 形式参数被声明为指向双精度指针的指针)。但是C函数不会对array_cptr中存储的C地址做这样的修改。对于显示的代码,存在不必要的间接级别。

  3. Fortran 代码将 array_fptrarray_cptr 中的 C 地址引用的事物相关联。因为array_cptr中的C地址没有被C函数改变,所以array_fptr继续与array关联。在这种情况下执行 C_F_POINTER 不会改变程序的状态。

array_cptrarray_fptr 不同 - 但它们具有相同的目标(array_cptr 是一个对象,array_fptr 与对象相关联 - 中的指针与 C 不同,Fortran 本身不是对象。