Fortran Coder

查看: 10340|回复: 5
打印 上一主题 下一主题

[混编] 麻烦大家审查下我的代码。提出宝贵意见

[复制链接]

13

帖子

8

主题

0

精华

入门

F 币
83 元
贡献
51 点
跳转到指定楼层
楼主
发表于 2016-1-12 20:45:58 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
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);
        }
分享到:  微信微信
收藏收藏 点赞点赞 点踩点踩

2015

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
1574 元
贡献
676 点

美女勋章热心勋章星光勋章新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

沙发
发表于 2016-1-13 10:20:41 | 只看该作者
你这个代码好绕啊,错误也很多。你想实现什么?我改了一下,但不知是否和你意?
[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);
        }


13

帖子

8

主题

0

精华

入门

F 币
83 元
贡献
51 点
板凳
 楼主| 发表于 2016-1-13 18:00:11 | 只看该作者
哦,我是想建立一个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);

13

帖子

8

主题

0

精华

入门

F 币
83 元
贡献
51 点
地板
 楼主| 发表于 2016-1-13 18:01:11 | 只看该作者
苏坡,man 发表于 2016-1-13 18:00
哦,我是想建立一个fortran的函数指针变量a,其中有一个fortran的子程序b(args),args是type(c_ptr),value:: ...

[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);
}

2015

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
1574 元
贡献
676 点

美女勋章热心勋章星光勋章新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

5#
发表于 2016-1-13 18:14:16 | 只看该作者
根据你的描述,沙发的代码即可实现。我看不清楚你代码里 __acemesh_task_generator_with_neighbors 和  stencil_core 是什么?

13

帖子

8

主题

0

精华

入门

F 币
83 元
贡献
51 点
6#
 楼主| 发表于 2016-1-13 19:08:09 | 只看该作者
我在想,这个函数指针变量不需要定义也可以。直接传递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);
  }
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-11-14 14:39

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表