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

苏坡，man

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

苏坡，man

我在想，这个函数指针变量不需要定义也可以。直接传递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);
  }

fcode

根据你的描述，沙发的代码即可实现。我看不清楚你代码里 __acemesh_task_generator_with_neighbors 和  stencil_core 是什么？

苏坡，man

苏坡，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);
}

苏坡，man

哦，我是想建立一个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);

fcode

你这个代码好绕啊，错误也很多。你想实现什么？我改了一下，但不知是否和你意？

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