麻烦大家审查下我的代码。提出宝贵意见
Testf.f:module tts
Use ISO_C_BINDING
Integer(c_int)pt,loc
Implicit none
contains
subroutinets(ats)
Use ISO_C_BINDING
class(c_ptr),pointer,value ::ats
pt = loc(ats)
print *,'pt=',pt
! return
End subroutine ts
End module tts
!********************************************
!********************************************
program main
use ISO_C_BINDING
use tts
Implicit none
Interface
!******************************************
subroutine ss(as) bind(C,name = 'ss_' )
class(c_ptr),value ::as
End subroutine ss
!******************************************
subroutine TASK_FUNCPTR(args)
class(c_ptr),value :: args
End subroutine TASK_FUNCPTR
Endinterface
!******************************************
procedure(TASK_FUNCPTR),pointer :: task => null() !指针的初始化
task => ts;
call ss(c_funptr(task))!c_funcptr将函数指针转换为地址传递给c
End
Testc.c:
#include <stdio.h>
void ss_(void *as){
printf(“as=%p\n”,as);
} 你这个代码好绕啊，错误也很多。你想实现什么？我改了一下，但不知是否和你意？module tts
Use ISO_C_BINDING
Implicit none
Integer(c_int) , pointer ::pt
contains
subroutine ts( ats ) Bind( C )
type(c_ptr),value ::ats
call c_f_pointer( ats , pt ) !// 把 C 语言的指针 ats 转换成 fortran 的指针 pt
print *,'pt=',pt
End subroutine ts
End module tts
!********************************************
!********************************************
program main
use ISO_C_BINDING
use tts
Implicit none
Interface
subroutine ss(as) bind(C,name = 'ss_' )
import
type(c_funptr),value ::as
End subroutine ss
End Interface
call ss( c_funloc(ts) )
End Program Main#include <stdio.h>
void ss_(void *as){
printf("s=%p\n",as);
}
哦，我是想建立一个fortran的函数指针变量a，其中有一个fortran的子程序b(args),args是type(c_ptr),value::x,然后将这个函数指针变量a=>b.再将a的地址传递给c函数即可。a.h:
typedef void (*TASK_FUNCPTR)(void *args);
void __acemesh_task_generator_with_neighbors(void* taskfptr,void* args,unsigned int args_size,
void* cxx_this_pointer, void* get_neighbors_funcptr, void *neighbor_args);
b.c:
__acemesh_task_generator_with_neighbors((TASK_FUNCPTR)stencil_core,
(void*)(&each), sizeof(targs),
NULL, (NEIGHBOR_FUNCPTR)my_neighbors, &each); 苏坡，man 发表于 2016-1-13 18:00
哦，我是想建立一个fortran的函数指针变量a，其中有一个fortran的子程序b(args),args是type(c_ptr),value:: ...
void stencil_core(targs* args)
{
int starti=args->starti;
int endi=args->endi;
int startj=args->startj;
int endj=args->endj;
int iter=args->iter;
int i,j,k;
int dest = iter% 2 ;
int src = 1 - dest;
for(i = starti; i < endi; i++)
for(j = startj; j < endj; j++)
for(k = 1; k < args->k_N - 1; k++)
Do3d7pkernel(dest,src,i,j,k);
} 根据你的描述，沙发的代码即可实现。我看不清楚你代码里 __acemesh_task_generator_with_neighbors 和stencil_core 是什么？ 我在想，这个函数指针变量不需要定义也可以。直接传递fortran子程序的地址即可,就是按照你写的代码意思。
在原先的c代码中，实现的方式如下：
c代码：
void stencil_core(targs* args)
{
int starti=args->starti;
int endi=args->endi;
int startj=args->startj;
int endj=args->endj;
int iter=args->iter;
int i,j,k;
int dest = iter% 2 ;
int src = 1 - dest;
for(i = starti; i < endi; i++)
for(j = startj; j < endj; j++)
for(k = 1; k < args->k_N - 1; k++)
Do3d7pkernel(dest,src,i,j,k);
}
int Do3d7p(const int bkn, const int iter, const int blocsizex, const int blocsizey)
{
__acemesh_task_generator_with_neighbors((TASK_FUNCPTR)stencil_core,
(void*)(&each), sizeof(targs),
NULL, (NEIGHBOR_FUNCPTR)my_neighbors, &each);
｝
头文件，b.h：
typedef void (*NEIGHBOR_FUNCPTR)(void** neighbor_addrs, int* num_neigh, void* args);
void __acemesh_task_generator_with_neighbors(void* taskfptr,void* args,unsigned int args_size,
void* cxx_this_pointer, void* get_neighbors_funcptr, void *neighbor_args);
b.cpp:
__acemesh_task_generator_with_neighbors(void* taskfptr,void* args,unsigned int args_size,
void *cxx_this_pointer, void* get_neighbors_funcptr, void *neighbor_args)
{
TASK_FUNCPTR funcptr = (TASK_FUNCPTR)taskfptr;
ci_task* _task = ACEMESH_NEW ci_task(funcptr);
_task->define(args, args_size);
_task->define_neighbors(cxx_this_pointer, get_neighbors_funcptr, neighbor_args);
}
页:
[1]