c - 尝试处理Fortran指针-调用c_f_pointer

我在使用后尝试处理Fortran指针时遇到问题：
call c_f_pointer(dataPtr, dataF, [natoms,countf])，其中dataPtr是带有intent(in)的C指针，而dataF是带有intent(out)的Fortran指针。

natoms和countf是整数，用于定义输出指针的大小。

变量先前被声明为：

use, intrinsic :: iso_c_binding
integer natoms, countf
type(c_ptr) dataPtr
real(c_float), pointer :: dataF(:,:) => NULL()

dataPtr来自另一个函数。

如果打印shape(dataF)，我将得到正确的结果。但是我无法使用dataF打印或复印或执行任何操作。我得到的错误是forrtl：严重（174）：SIGSEGV，发生分段错误。

我按照示例https://software.intel.com/en-us/node/678433进行操作，可以毫无问题地进行编译和操作。



编辑：添加更多的代码。

program test

use, intrinsic :: iso_c_binding
implicit none

interface
    subroutine open(ptr) bind(c, name='open')
        use, intrinsic :: iso_c_binding
        implicit none
        type(c_ptr), intent(out) :: ptr
    end subroutine open
    subroutine gather(ptr, name, typ, countf, data) bind(c, name='gather')
        use, intrinsic :: iso_c_binding
        implicit none
        type(c_ptr), value, intent(in) :: ptr
        type(c_ptr), value, intent(in) :: name
        integer(c_int), value, intent(in) :: typ
        integer(c_int), value, intent(in) :: countf
        type(c_ptr), intent(out) :: data
    end subroutine gather
    subroutine scatter(ptr, name, typ, countf, data) bind(c, name='scatter')
        use, intrinsic :: iso_c_binding
        implicit none
        type(c_ptr), value, intent(in) :: ptr
        type(c_ptr), value, intent(in) :: name
        integer(c_int), value, intent(in) :: typ
        integer(c_int), value, intent(in) :: countf
        type(c_ptr), intent(in) :: data
    end subroutine scatter
end interface

integer natoms, countf, typ
character(len=2,kind=c_char) :: var = c_char_'x'//c_null_char
type(c_ptr) ptr, dataPtr
real(c_float), pointer :: dataF(:,:) => NULL()

!opens an instance of the program (loaded as library into my Fortran)
!ptr is a pointer used to reference to it in other calls
call open(ptr)

!Here comes definitions for the integer variables
!...

!get data
call gather(ptr, c_loc(var), typ, countf, dataPtr)

!need to do stuff with the data here
call c_f_pointer(dataPtr, dataF, [natoms,countf])

!return modified data
call scatter(ptr, c_loc(var), typ, countf, dataPtr)

如果我只是获取数据call gather(...)，然后将其返回call scatter(...)，则程序可以正常工作。但是我需要在两者之间做一些事情。



EDIT2：添加了C函数

void open(void **ptr) {
// ...
}

void gather(void *ptr, const char *name,
                     int type, int count, void *data) {
// ...
}

void scatter(void *ptr, const char *name,
                     int type, int count, void *data) {
// ...
}

最佳答案

显然我明白了。

首先，数据指针必须在子程序intent(in)中输入为value和gather：

subroutine gather(ptr, name, typ, countf, data) bind(c, name='gather')
    ! ...
    type(c_ptr), value, intent(in) :: data
    ! instead of:
    ! type(c_ptr), intent(out) :: data
end subroutine gather

我更改了dataF的声明，不再需要从c_f_pointer进行调用：

real(c_float), target, allocatable :: dataF(:,:)

所以我删除了call c_f_pointer(dataPtr, dataF, [natoms,countf])并添加了：

allocate(dataF(natoms,countf))
dataPtr = c_loc(dataCFloat)

现在，dataF直接包含更新后的数组。

关于c - 尝试处理Fortran指针-调用c_f_pointer，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/50373861/