Fortran Coder
标题: 麻烦大家审查下我的代码。提出宝贵意见 [打印本页]
作者: 苏坡，man 时间: 2016-1-12 20:45
标题: 麻烦大家审查下我的代码。提出宝贵意见
Testf.f:[Fortran] 纯文本查看 复制代码
module tts
Use ISO_C_BINDING
Integer(c_int) pt,loc
Implicit none
contains
subroutine ts(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
End interface
!******************************************
procedure(TASK_FUNCPTR),pointer :: task => null() !指针的初始化
task => ts;
call ss(c_funptr(task)) !c_funcptr将函数指针转换为地址传递给c
End
Testc.c:[C] 纯文本查看 复制代码
#include <stdio.h>
void ss_(void *as){
printf(“as=%p\n”,as);
}
作者: fcode 时间: 2016-1-13 10:20
你这个代码好绕啊，错误也很多。你想实现什么？我改了一下，但不知是否和你意？[Fortran] 纯文本查看 复制代码
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
[C] 纯文本查看 复制代码
#include <stdio.h>
void ss_(void *as){
printf("s=%p\n",as);
}
作者: 苏坡，man 时间: 2016-1-13 18:00
哦，我是想建立一个fortran的函数指针变量a，其中有一个fortran的子程序b(args),args是type(c_ptr),value::x,然后将这个函数指针变量a=>b.再将a的地址传递给c函数即可。[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:01
[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);
}
作者: fcode 时间: 2016-1-13 18:14
根据你的描述，沙发的代码即可实现。我看不清楚你代码里 __acemesh_task_generator_with_neighbors 和 stencil_core 是什么？
作者: 苏坡，man 时间: 2016-1-13 19:08
我在想，这个函数指针变量不需要定义也可以。直接传递fortran子程序的地址即可,就是按照你写的代码意思。
在原先的c代码中，实现的方式如下：
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：[C] 纯文本查看 复制代码
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);
}
| 欢迎光临 Fortran Coder (http://bbs.fcode.cn/) |
Powered by Discuz! X3.2 |