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

Question

我在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

为了确定，流程是

array_cptr 与具有 C_LOC 功能的 array 关联。（也许 'associate' 在这里不是一个合适的术语，因为它是 Fortran 指针相关的术语。我应该怎么称呼它？）
array 在我调用 get_value_array 时被修改，因为子例程（或 C 函数）修改了 array_cptr 指向的位置。

C_F_POINTER

array_fptr 通过 array_cptr 关联到 array

最后 array_cptr 和 array_fptr 是彼此不同的对象 (?)，同时指向同一个目标 array，对吗？

Answer 1

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

用Fortran语句array_cptr = C_LOC(array)，array_cptr被赋值为array.
的C地址
调用get_value_array时，array的值最终被C函数修改。 Fortran 代码还允许子例程（或 C 函数）修改存储在 array_cptr 中的 C 地址 - 伪参数声明为 INTENT(INOUT) 并且没有 VALUE 属性，并且与C 函数通过引用获取指向数组的指针（即 C 形式参数被声明为指向双精度指针的指针）。但是C函数不会对array_cptr中存储的C地址做这样的修改。对于显示的代码，存在不必要的间接级别。
Fortran 代码将 array_fptr 与 array_cptr 中的 C 地址引用的事物相关联。因为array_cptr中的C地址没有被C函数改变，所以array_fptr继续与array关联。在这种情况下执行 C_F_POINTER 不会改变程序的状态。

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

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

Check for ISO_C_BINDING usage in Fortran calling C function

c

fortran

pointers

fortran-iso-c-binding