苏坡,man 发表于 2016-4-11 20:22:21

关于c_f_pointer的使用

我下面有一段代码,是想c指针转换为fortran指针,然后在c代码那边可以调用fortran的变量的值。代码如下: 我认为我的老师使用的c_f_pointer不对。 我想请教下大家应该如何将数组地址的值传递给 c中的void *addr使用?
SUBROUTINE calc1_neighbor1(neighbor_addrs,args)
    use ISO_C_BINDING
    IMPLICIT NONE
#include "aceMesh_runtime_f.h"
type,bind(c) :: ptr_array
       integer :: len
      type(c_ptr) :: addr!在c中是定义的 void *addr类型
   end
    real(c_double) src,src1
    type(calc1_type1),target :: args
    integer(c_int) ii,jj,size,num
type(ptr_array) ::neighbor_addrs!派生类型变量
    type(c_ptr),pointer::array(:)!定义c指针
    size=MAX_NEIGHBORS

call c_f_pointer(neighbor_addrs%arr,array,) !使用该转换函数使得c指针转换为f指针,然后将array数组指针传递给c代码使用
    ii=args%ii
    jj=args%jj
    num=0
    if(ii+iblock<=M)then
      num=num+1
    array(num)=c_loc(p(ii+iblock,jj))   !获得fortran数组p的地址
    endif

neighbor_addrs%len=num !传递数量num给c使用
    RETURN
    END SUBROUTINE calc1_neighbor1

我应该如何获得数组的地址,传递给c代码中的void *addr变量使用?并且需要传递计数的num变量给c代码使用?

vvt 发表于 2016-4-11 21:24:03

看不懂你的代码啊。。。。
你是 C 调用 fortran,还是 fortran 调用 C?

calc1_neighbor1(neighbor_addrs,args)
这里是给 C 调用的?两个参数分别什么含义?

苏坡,man 发表于 2016-4-11 21:44:13

不想使用c_f_pointer。


我下面有一段代码,是想c指针转换为fortran指针,然后在c代码那边可以调用fortran的变量的值。代码如下: 我认为我的老师使用的c_f_pointer不对。 我想请教下大家应该如何将数组地址的值传递给 c中的void *addr使用。

简单的说是: c调用fortran,在fortran中有数组p需要传递数组每个地址给c代码的结构ptr_array中的void *addr,当addr得到地址后再使用地址进行相应的计算。
我应该如何获得数组p的地址,传递给c代码中的void *addr变量使用?并且需要传递计数的num变量给c代码使用?

neighbor_addrs: type(ptr_array) ::neighbor_addrs
args: type(calc1_type1),target :: args
calc1_type:
type calc1_type1   
integer(c_int) ii,jj   end type calc1_type1摘去的部分代码:
SUBROUTINE calc1_neighbor1(neighbor_addrs,args)
    use ISO_C_BINDING
    IMPLICIT NONE
#include "aceMesh_runtime_f.h"
   type,bind(c) :: ptr_array
       integer :: len
      type(c_ptr) :: addr!在c中是定义的 void *addr类型
   end
    real(c_double) src,src1
    type(calc1_type1),target :: args
    integer(c_int) ii,jj,size,num
type(ptr_array) ::neighbor_addrs!派生类型变量
    type(c_ptr),pointer::array(:)!定义c指针
    size=MAX_NEIGHBORS

call c_f_pointer(neighbor_addrs%arr,array,) !使用该转换函数使得c指针转换为f指针,然后将array数组指针传递给c代码使用
    ii=args%ii
    jj=args%jj
    num=0
    if(ii+iblock<=M)then
      num=num+1
      array(num)=c_loc(p(ii+iblock,jj))   !获得fortran数组p的地址
    endif

neighbor_addrs%len=num !传递数量num给c使用
    RETURN
    END SUBROUTINE calc1_neighbor1


我的想法是将fortran代码中的type(c_ptr) :: addr(该代码在c中是被定义为void *addr,注释:不能修改c代码),改为 integer(c_int) :: addr(32),
然后    array(num)=c_loc(p(ii+iblock,jj)) 改为 addr(num)=loc(p(ii+iblock,jj)),不知我的想法如何?

pasuka 发表于 2016-4-11 22:51:24

请问lz编译出错还是运行报错?
获得Fortran数组p的地址这行没有看出语法错误啊~数组是连续分配内存的话,只要指针指向某个元素即可,是否数组越界那是另外一个检查
lz老师的用法没发现有问题,gofrtran的手册说得很明白:If a pointer is a dummy-argument of an interoperable procedure, it usually has to be declared using the VALUE attribute. void* matches TYPE(C_PTR), VALUE, while TYPE(C_PTR) alone matches void**.
Name        Type
C_NULL_PTR        C_PTR
C_NULL_FUNPTR        C_FUNPTR
Both are equivalent to the value NULL in C.
页: [1]
查看完整版本: 关于c_f_pointer的使用