src/smpi/smpi_f77.cpp

   1 /* Copyright (c) 2010-2015. The SimGrid Team.
   2  * All rights reserved.                                                     */
   3
   4 /* This program is free software; you can redistribute it and/or modify it
   5  * under the terms of the license (GNU LGPL) which comes with this package. */
   6
   7 #include <limits.h>
   8 #include <stdio.h>
   9 #include "smpi_mpi_dt_private.h"
  10 #include "private.h"
  11 #include "xbt.h"
  12
  13 //TODO : remove all this.
  14 static xbt_dict_t comm_lookup = nullptr;
  15 static xbt_dict_t group_lookup = nullptr;
  16 static xbt_dict_t request_lookup = nullptr;
  17 static xbt_dict_t type_lookup = nullptr;
  18 static xbt_dict_t op_lookup = nullptr;
  19 static xbt_dict_t win_lookup = nullptr;
  20 static xbt_dict_t info_lookup = nullptr;
  21 static int comm_id=0;
  22 static int group_id=0;
  23 static int request_id=0;
  24 static int type_id=0;
  25 static int op_id=0;
  26 static int win_id=0;
  27 static int info_id=0;
  28
  29
  30 static int running_processes = 0;
  31
  32 #if defined(__alpha__) || defined(__sparc64__) || defined(__x86_64__) || defined(__ia64__)
  33 typedef int integer;
  34 typedef unsigned int uinteger;
  35 #else
  36 typedef long int integer;
  37 typedef unsigned long int uinteger;
  38 #endif
  39
  40 /* Bindings for MPI special values */
  41   struct s_smpi_common {
  42     integer _MPI_IN_PLACE;
  43     integer _MPI_BOTTOM;
  44     integer _MPI_STATUS_IGNORE;
  45     integer _MPI_STATUSES_IGNORE;
  46   } smpi_;
  47
  48 /* Convert between Fortran and C */
  49 #define FORT_ADDR(addr, val)                                         \
  50   (((void *)(addr) == (void*) &(smpi_._ ## val))                  \
  51    ? (val) : (void *)(addr))
  52 #define FORT_BOTTOM(addr)          FORT_ADDR(addr, MPI_BOTTOM)
  53 #define FORT_IN_PLACE(addr)        FORT_ADDR(addr, MPI_IN_PLACE)
  54 #define FORT_STATUS_IGNORE(addr)   static_cast<MPI_Status*>(FORT_ADDR(addr, MPI_STATUS_IGNORE))
  55 #define FORT_STATUSES_IGNORE(addr) static_cast<MPI_Status*>(FORT_ADDR(addr, MPI_STATUSES_IGNORE))
  56
  57 #define KEY_SIZE (sizeof(int) * 2 + 1)
  58
  59 static char* get_key(char* key, int id) {
  60   snprintf(key, KEY_SIZE, "%x",id);
  61   return key;
  62 }
  63
  64 static char* get_key_id(char* key, int id) {
  65   snprintf(key, KEY_SIZE, "%x_%d",id, smpi_process_index());
  66   return key;
  67 }
  68
  69 static void smpi_init_fortran_types(){
  70    if(comm_lookup == nullptr){
  71      comm_lookup = xbt_dict_new_homogeneous(nullptr);
  72      smpi_comm_add_f(MPI_COMM_WORLD);
  73      group_lookup = xbt_dict_new_homogeneous(nullptr);
  74      request_lookup = xbt_dict_new_homogeneous(nullptr);
  75      type_lookup = xbt_dict_new_homogeneous(nullptr);
  76      win_lookup = xbt_dict_new_homogeneous(nullptr);
  77      info_lookup = xbt_dict_new_homogeneous(nullptr);
  78      smpi_type_add_f(MPI_BYTE);//MPI_BYTE
  79      smpi_type_add_f(MPI_CHAR);//MPI_CHARACTER
  80 #if defined(__alpha__) || defined(__sparc64__) || defined(__x86_64__) || defined(__ia64__)
  81      smpi_type_add_f(MPI_C_BOOL);//MPI_LOGICAL
  82      smpi_type_add_f(MPI_INT);//MPI_INTEGER
  83 #else
  84      smpi_type_add_f(MPI_C_BOOL);//MPI_LOGICAL
  85      smpi_type_add_f(MPI_LONG);//MPI_INTEGER
  86 #endif
  87      smpi_type_add_f(MPI_INT8_T);//MPI_INTEGER1
  88      smpi_type_add_f(MPI_INT16_T);//MPI_INTEGER2
  89      smpi_type_add_f(MPI_INT32_T);//MPI_INTEGER4
  90      smpi_type_add_f(MPI_INT64_T);//MPI_INTEGER8
  91      smpi_type_add_f(MPI_REAL);//MPI_REAL
  92      smpi_type_add_f(MPI_REAL4);//MPI_REAL4
  93      smpi_type_add_f(MPI_REAL8);//MPI_REAL8
  94      smpi_type_add_f(MPI_DOUBLE);//MPI_DOUBLE_PRECISION
  95      smpi_type_add_f(MPI_C_FLOAT_COMPLEX);//MPI_COMPLEX
  96      smpi_type_add_f(MPI_C_DOUBLE_COMPLEX);//MPI_DOUBLE_COMPLEX
  97 #if defined(__alpha__) || defined(__sparc64__) || defined(__x86_64__) || defined(__ia64__)
  98      smpi_type_add_f(MPI_2INT);//MPI_2INTEGER
  99 #else
 100      smpi_type_add_f(MPI_2LONG);//MPI_2INTEGER
 101 #endif
 102      smpi_type_add_f(MPI_UINT8_T);//MPI_LOGICAL1
 103      smpi_type_add_f(MPI_UINT16_T);//MPI_LOGICAL2
 104      smpi_type_add_f(MPI_UINT32_T);//MPI_LOGICAL4
 105      smpi_type_add_f(MPI_UINT64_T);//MPI_LOGICAL8
 106      smpi_type_add_f(MPI_2FLOAT);//MPI_2REAL
 107      smpi_type_add_f(MPI_2DOUBLE);//MPI_2DOUBLE_PRECISION
 108      smpi_type_add_f(MPI_PTR);//MPI_AINT
 109      smpi_type_add_f(MPI_OFFSET);//MPI_OFFSET
 110      smpi_type_add_f(MPI_AINT);//MPI_COUNT
 111      smpi_type_add_f(MPI_DATATYPE_NULL);//MPI_REAL16
 112      smpi_type_add_f(MPI_DATATYPE_NULL);//MPI_COMPLEX8
 113      smpi_type_add_f(MPI_DATATYPE_NULL);//MPI_COMPLEX16
 114      smpi_type_add_f(MPI_DATATYPE_NULL);//MPI_COMPLEX32
 115      smpi_type_add_f(MPI_DATATYPE_NULL);//MPI_PACKED
 116      op_lookup = xbt_dict_new_homogeneous(nullptr);
 117      smpi_op_add_f(MPI_MAX);
 118      smpi_op_add_f(MPI_MIN);
 119      smpi_op_add_f(MPI_MAXLOC);
 120      smpi_op_add_f(MPI_MINLOC);
 121      smpi_op_add_f(MPI_SUM);
 122      smpi_op_add_f(MPI_PROD);
 123      smpi_op_add_f(MPI_LAND);
 124      smpi_op_add_f(MPI_LOR);
 125      smpi_op_add_f(MPI_LXOR);
 126      smpi_op_add_f(MPI_BAND);
 127      smpi_op_add_f(MPI_BOR);
 128      smpi_op_add_f(MPI_BXOR);
 129    }
 130 }
 131
 132 template <typename T> int smpi_add_f(T t, xbt_dict_t dict, int* id) {
 133   char key[KEY_SIZE];
 134   xbt_dict_set(dict, get_key(key, *id), t, nullptr);
 135   (*id)++;
 136   return *id-1;
 137 }
 138
 139 template <> int smpi_add_f<MPI_Comm>(MPI_Comm comm, xbt_dict_t dict, int* id) {
 140   char key[KEY_SIZE];
 141   xbt_dict_set(dict, comm==MPI_COMM_WORLD? get_key(key, *id) : get_key_id(key, *id), comm, nullptr);
 142   (*id)++;
 143   return *id-1;
 144 }
 145
 146 template <> int smpi_add_f<MPI_Request>(MPI_Request request, xbt_dict_t dict, int* id) {
 147   char key[KEY_SIZE];
 148   xbt_dict_set(dict, get_key_id(key, *id), request, nullptr);
 149   (*id)++;
 150   return *id-1;
 151 }
 152
 153 template <typename T> int smpi_c2f(T t, xbt_dict_t dict, int* id) {
 154   char* existing_key = xbt_dict_get_key(dict, t);
 155   if(existing_key!=nullptr)
 156     return atoi(existing_key);
 157   else
 158     return smpi_add_f<T>(t,dict,id);
 159 }
 160
 161 template <typename T> void free_f(int id, xbt_dict_t dict) {
 162   char key[KEY_SIZE];
 163   xbt_dict_remove(dict, get_key(key, id));
 164 }
 165
 166 template <> void free_f<MPI_Comm>(int id, xbt_dict_t dict) {
 167   char key[KEY_SIZE];
 168   xbt_dict_remove(dict, id==0? get_key(key, id) : get_key_id(key, id));
 169 }
 170
 171 template <> void free_f<MPI_Request>(int id, xbt_dict_t dict) {
 172   char key[KEY_SIZE];
 173   if(id!=MPI_FORTRAN_REQUEST_NULL)
 174   xbt_dict_remove(dict, get_key_id(key, id));
 175 }
 176
 177 template <typename T> T smpi_f2c(int , xbt_dict_t ) {
 178   return NULL;
 179 }
 180
 181 template <> MPI_Datatype smpi_f2c<MPI_Datatype>(int id, xbt_dict_t dict) {
 182   smpi_init_fortran_types();
 183   char key[KEY_SIZE];
 184   return id >= 0 ? static_cast<MPI_Datatype>(xbt_dict_get_or_null(dict, get_key(key, id))): MPI_DATATYPE_NULL;
 185 }
 186
 187 template <> MPI_Op smpi_f2c<MPI_Op> (int id, xbt_dict_t dict) {
 188   smpi_init_fortran_types();
 189   char key[KEY_SIZE];
 190   return id >= 0 ? static_cast<MPI_Op>(xbt_dict_get_or_null(dict, get_key(key, id))): MPI_OP_NULL;
 191 }
 192
 193 template <> MPI_Info smpi_f2c<MPI_Info>(int id, xbt_dict_t dict) {
 194   smpi_init_fortran_types();
 195   char key[KEY_SIZE];
 196   return id >= 0 ? static_cast<MPI_Info>(xbt_dict_get_or_null(dict, get_key(key, id))): MPI_INFO_NULL;
 197 }
 198
 199 template <> MPI_Win smpi_f2c<MPI_Win>(int id, xbt_dict_t dict) {
 200   smpi_init_fortran_types();
 201   char key[KEY_SIZE];
 202   return id >= 0 ? static_cast<MPI_Win>(xbt_dict_get_or_null(dict, get_key(key, id))): MPI_WIN_NULL;
 203 }
 204
 205 template <> MPI_Comm smpi_f2c<MPI_Comm>(int comm, xbt_dict_t dict) {
 206   smpi_init_fortran_types();
 207   if(comm == -2) {
 208     return MPI_COMM_SELF;
 209   } else if(comm==0){
 210     return MPI_COMM_WORLD;
 211   } else if(dict != nullptr && comm >= 0) {
 212       char key[KEY_SIZE];
 213       MPI_Comm tmp =  static_cast<MPI_Comm>(xbt_dict_get_or_null(dict,get_key_id(key, comm)));
 214       return tmp != nullptr ? tmp : MPI_COMM_NULL ;
 215   } else {
 216     return MPI_COMM_NULL;
 217   }
 218 }
 219
 220 template <> MPI_Group smpi_f2c<MPI_Group>(int group, xbt_dict_t dict) {
 221   smpi_init_fortran_types();
 222   if(group == -2) {
 223     return MPI_GROUP_EMPTY;
 224   } else if(dict != nullptr && group >= 0) {
 225     char key[KEY_SIZE];
 226     return static_cast<MPI_Group>(xbt_dict_get_or_null(dict, get_key(key, group)));
 227   } else {
 228     return MPI_GROUP_NULL;
 229   }
 230 }
 231
 232 template <> MPI_Request smpi_f2c<MPI_Request>(int request, xbt_dict_t dict) {
 233   smpi_init_fortran_types();
 234   char key[KEY_SIZE];
 235   if(request==MPI_FORTRAN_REQUEST_NULL)
 236     return MPI_REQUEST_NULL;
 237   return static_cast<MPI_Request>(xbt_dict_get(dict, get_key_id(key, request)));
 238 }
 239
 240 #define SMPI_F2C_C2F(type, name)\
 241 int smpi_##name##_add_f(type name){\
 242   return smpi_add_f<type>(name, name##_lookup, &name##_id);\
 243 }\
 244 int smpi_##name##_c2f(type name){\
 245   return smpi_c2f<type>(name, name##_lookup, &name##_id);\
 246 }\
 247 static void free_##name(int id) {\
 248   free_f<type>(id, name##_lookup);\
 249 }\
 250 type smpi_##name##_f2c(int id){\
 251   return smpi_f2c<type>(id, name##_lookup);\
 252 }
 253
 254 extern "C" { // This should really use the C linkage to be usable from Fortran
 255
 256 SMPI_F2C_C2F(MPI_Comm, comm)
 257 SMPI_F2C_C2F(MPI_Group, group)
 258 SMPI_F2C_C2F(MPI_Request, request)
 259 SMPI_F2C_C2F(MPI_Datatype, type)
 260 SMPI_F2C_C2F(MPI_Win, win)
 261 SMPI_F2C_C2F(MPI_Op, op)
 262 SMPI_F2C_C2F(MPI_Info, info)
 263
 264 void mpi_init_(int* ierr) {
 265     smpi_init_fortran_types();
 266    *ierr = MPI_Init(nullptr, nullptr);
 267    running_processes++;
 268 }
 269
 270 void mpi_finalize_(int* ierr) {
 271    *ierr = MPI_Finalize();
 272    running_processes--;
 273    if(running_processes==0){
 274      xbt_dict_free(&op_lookup);
 275      xbt_dict_free(&type_lookup);
 276      xbt_dict_free(&request_lookup);
 277      xbt_dict_free(&group_lookup);
 278      xbt_dict_free(&comm_lookup);
 279      xbt_dict_free(&win_lookup);
 280      xbt_dict_free(&info_lookup);
 281    }
 282 }
 283
 284 void mpi_abort_(int* comm, int* errorcode, int* ierr) {
 285   *ierr = MPI_Abort(smpi_comm_f2c(*comm), *errorcode);
 286 }
 287
 288 void mpi_comm_rank_(int* comm, int* rank, int* ierr) {
 289    *ierr = MPI_Comm_rank(smpi_comm_f2c(*comm), rank);
 290 }
 291
 292 void mpi_comm_size_(int* comm, int* size, int* ierr) {
 293    *ierr = MPI_Comm_size(smpi_comm_f2c(*comm), size);
 294 }
 295
 296 double mpi_wtime_() {
 297    return MPI_Wtime();
 298 }
 299
 300 double mpi_wtick_() {
 301   return MPI_Wtick();
 302 }
 303
 304 void mpi_comm_dup_(int* comm, int* newcomm, int* ierr) {
 305   MPI_Comm tmp;
 306
 307   *ierr = MPI_Comm_dup(smpi_comm_f2c(*comm), &tmp);
 308   if(*ierr == MPI_SUCCESS) {
 309     *newcomm = smpi_comm_add_f(tmp);
 310   }
 311 }
 312
 313 void mpi_comm_create_(int* comm, int* group, int* newcomm, int* ierr) {
 314   MPI_Comm tmp;
 315
 316   *ierr = MPI_Comm_create(smpi_comm_f2c(*comm),smpi_group_f2c(*group), &tmp);
 317   if(*ierr == MPI_SUCCESS) {
 318     *newcomm = smpi_comm_add_f(tmp);
 319   }
 320 }
 321
 322 void mpi_comm_free_(int* comm, int* ierr) {
 323   MPI_Comm tmp = smpi_comm_f2c(*comm);
 324
 325   *ierr = MPI_Comm_free(&tmp);
 326
 327   if(*ierr == MPI_SUCCESS) {
 328     free_comm(*comm);
 329   }
 330 }
 331
 332 void mpi_comm_split_(int* comm, int* color, int* key, int* comm_out, int* ierr) {
 333   MPI_Comm tmp;
 334
 335   *ierr = MPI_Comm_split(smpi_comm_f2c(*comm), *color, *key, &tmp);
 336   if(*ierr == MPI_SUCCESS) {
 337     *comm_out = smpi_comm_add_f(tmp);
 338   }
 339 }
 340
 341 void mpi_group_incl_(int* group, int* n, int* ranks, int* group_out, int* ierr) {
 342   MPI_Group tmp;
 343
 344   *ierr = MPI_Group_incl(smpi_group_f2c(*group), *n, ranks, &tmp);
 345   if(*ierr == MPI_SUCCESS) {
 346     *group_out = smpi_group_add_f(tmp);
 347   }
 348 }
 349
 350 void mpi_comm_group_(int* comm, int* group_out,  int* ierr) {
 351   MPI_Group tmp;
 352
 353   *ierr = MPI_Comm_group(smpi_comm_f2c(*comm), &tmp);
 354   if(*ierr == MPI_SUCCESS) {
 355     *group_out = smpi_group_c2f(tmp);
 356   }
 357 }
 358
 359 void mpi_initialized_(int* flag, int* ierr){
 360   *ierr = MPI_Initialized(flag);
 361 }
 362
 363 void mpi_send_init_(void *buf, int* count, int* datatype, int* dst, int* tag, int* comm, int* request, int* ierr) {
 364   MPI_Request req;
 365   buf = static_cast<char *>(FORT_BOTTOM(buf));
 366   *ierr = MPI_Send_init(buf, *count, smpi_type_f2c(*datatype), *dst, *tag, smpi_comm_f2c(*comm), &req);
 367   if(*ierr == MPI_SUCCESS) {
 368     *request = smpi_request_add_f(req);
 369   }
 370 }
 371
 372 void mpi_isend_(void *buf, int* count, int* datatype, int* dst, int* tag, int* comm, int* request, int* ierr) {
 373   MPI_Request req;
 374   buf = static_cast<char *>(FORT_BOTTOM(buf));
 375   *ierr = MPI_Isend(buf, *count, smpi_type_f2c(*datatype), *dst, *tag, smpi_comm_f2c(*comm), &req);
 376   if(*ierr == MPI_SUCCESS) {
 377     *request = smpi_request_add_f(req);
 378   }
 379 }
 380
 381 void mpi_irsend_(void *buf, int* count, int* datatype, int* dst, int* tag, int* comm, int* request, int* ierr) {
 382   MPI_Request req;
 383   buf = static_cast<char *>(FORT_BOTTOM(buf));
 384   *ierr = MPI_Irsend(buf, *count, smpi_type_f2c(*datatype), *dst, *tag, smpi_comm_f2c(*comm), &req);
 385   if(*ierr == MPI_SUCCESS) {
 386     *request = smpi_request_add_f(req);
 387   }
 388 }
 389
 390 void mpi_send_(void* buf, int* count, int* datatype, int* dst, int* tag, int* comm, int* ierr) {
 391   buf = static_cast<char *>(FORT_BOTTOM(buf));
 392    *ierr = MPI_Send(buf, *count, smpi_type_f2c(*datatype), *dst, *tag, smpi_comm_f2c(*comm));
 393 }
 394
 395 void mpi_rsend_(void* buf, int* count, int* datatype, int* dst, int* tag, int* comm, int* ierr) {
 396   buf = static_cast<char *>(FORT_BOTTOM(buf));
 397    *ierr = MPI_Rsend(buf, *count, smpi_type_f2c(*datatype), *dst, *tag, smpi_comm_f2c(*comm));
 398 }
 399
 400 void mpi_sendrecv_(void* sendbuf, int* sendcount, int* sendtype, int* dst, int* sendtag, void *recvbuf, int* recvcount,
 401                    int* recvtype, int* src, int* recvtag, int* comm, MPI_Status* status, int* ierr) {
 402   sendbuf = static_cast<char *>( FORT_BOTTOM(sendbuf));
 403   recvbuf = static_cast<char *>( FORT_BOTTOM(recvbuf));
 404    *ierr = MPI_Sendrecv(sendbuf, *sendcount, smpi_type_f2c(*sendtype), *dst, *sendtag, recvbuf, *recvcount,
 405                         smpi_type_f2c(*recvtype), *src, *recvtag, smpi_comm_f2c(*comm), FORT_STATUS_IGNORE(status));
 406 }
 407
 408 void mpi_recv_init_(void *buf, int* count, int* datatype, int* src, int* tag, int* comm, int* request, int* ierr) {
 409   MPI_Request req;
 410   buf = static_cast<char *>( FORT_BOTTOM(buf));
 411   *ierr = MPI_Recv_init(buf, *count, smpi_type_f2c(*datatype), *src, *tag, smpi_comm_f2c(*comm), &req);
 412   if(*ierr == MPI_SUCCESS) {
 413     *request = smpi_request_add_f(req);
 414   }
 415 }
 416
 417 void mpi_irecv_(void *buf, int* count, int* datatype, int* src, int* tag, int* comm, int* request, int* ierr) {
 418   MPI_Request req;
 419   buf = static_cast<char *>( FORT_BOTTOM(buf));
 420   *ierr = MPI_Irecv(buf, *count, smpi_type_f2c(*datatype), *src, *tag, smpi_comm_f2c(*comm), &req);
 421   if(*ierr == MPI_SUCCESS) {
 422     *request = smpi_request_add_f(req);
 423   }
 424 }
 425
 426 void mpi_recv_(void* buf, int* count, int* datatype, int* src, int* tag, int* comm, MPI_Status* status, int* ierr) {
 427   buf = static_cast<char *>( FORT_BOTTOM(buf));
 428   *ierr = MPI_Recv(buf, *count, smpi_type_f2c(*datatype), *src, *tag, smpi_comm_f2c(*comm), status);
 429 }
 430
 431 void mpi_start_(int* request, int* ierr) {
 432   MPI_Request req = smpi_request_f2c(*request);
 433
 434   *ierr = MPI_Start(&req);
 435 }
 436
 437 void mpi_startall_(int* count, int* requests, int* ierr) {
 438   MPI_Request* reqs;
 439   int i;
 440
 441   reqs = xbt_new(MPI_Request, *count);
 442   for(i = 0; i < *count; i++) {
 443     reqs[i] = smpi_request_f2c(requests[i]);
 444   }
 445   *ierr = MPI_Startall(*count, reqs);
 446   xbt_free(reqs);
 447 }
 448
 449 void mpi_wait_(int* request, MPI_Status* status, int* ierr) {
 450    MPI_Request req = smpi_request_f2c(*request);
 451
 452    *ierr = MPI_Wait(&req, FORT_STATUS_IGNORE(status));
 453    if(req==MPI_REQUEST_NULL){
 454      free_request(*request);
 455      *request=MPI_FORTRAN_REQUEST_NULL;
 456    }
 457 }
 458
 459 void mpi_waitany_(int* count, int* requests, int* index, MPI_Status* status, int* ierr) {
 460   MPI_Request* reqs;
 461   int i;
 462
 463   reqs = xbt_new(MPI_Request, *count);
 464   for(i = 0; i < *count; i++) {
 465     reqs[i] = smpi_request_f2c(requests[i]);
 466   }
 467   *ierr = MPI_Waitany(*count, reqs, index, status);
 468   if(reqs[*index]==MPI_REQUEST_NULL){
 469       free_request(requests[*index]);
 470       requests[*index]=MPI_FORTRAN_REQUEST_NULL;
 471   }
 472   xbt_free(reqs);
 473 }
 474
 475 void mpi_waitall_(int* count, int* requests, MPI_Status* status, int* ierr) {
 476   MPI_Request* reqs;
 477   int i;
 478
 479   reqs = xbt_new(MPI_Request, *count);
 480   for(i = 0; i < *count; i++) {
 481     reqs[i] = smpi_request_f2c(requests[i]);
 482   }
 483   *ierr = MPI_Waitall(*count, reqs, FORT_STATUSES_IGNORE(status));
 484   for(i = 0; i < *count; i++) {
 485       if(reqs[i]==MPI_REQUEST_NULL){
 486           free_request(requests[i]);
 487           requests[i]=MPI_FORTRAN_REQUEST_NULL;
 488       }
 489   }
 490
 491   xbt_free(reqs);
 492 }
 493
 494 void mpi_barrier_(int* comm, int* ierr) {
 495   *ierr = MPI_Barrier(smpi_comm_f2c(*comm));
 496 }
 497
 498 void mpi_bcast_(void *buf, int* count, int* datatype, int* root, int* comm, int* ierr) {
 499   *ierr = MPI_Bcast(buf, *count, smpi_type_f2c(*datatype), *root, smpi_comm_f2c(*comm));
 500 }
 501
 502 void mpi_reduce_(void* sendbuf, void* recvbuf, int* count, int* datatype, int* op, int* root, int* comm, int* ierr) {
 503   sendbuf = static_cast<char *>( FORT_IN_PLACE(sendbuf));
 504   sendbuf = static_cast<char *>( FORT_BOTTOM(sendbuf));
 505   recvbuf = static_cast<char *>( FORT_BOTTOM(recvbuf));
 506   *ierr = MPI_Reduce(sendbuf, recvbuf, *count, smpi_type_f2c(*datatype), smpi_op_f2c(*op), *root, smpi_comm_f2c(*comm));
 507 }
 508
 509 void mpi_allreduce_(void* sendbuf, void* recvbuf, int* count, int* datatype, int* op, int* comm, int* ierr) {
 510   sendbuf = static_cast<char *>( FORT_IN_PLACE(sendbuf));
 511   *ierr = MPI_Allreduce(sendbuf, recvbuf, *count, smpi_type_f2c(*datatype), smpi_op_f2c(*op), smpi_comm_f2c(*comm));
 512 }
 513
 514 void mpi_reduce_scatter_(void* sendbuf, void* recvbuf, int* recvcounts, int* datatype, int* op, int* comm, int* ierr) {
 515   sendbuf = static_cast<char *>( FORT_IN_PLACE(sendbuf));
 516   *ierr = MPI_Reduce_scatter(sendbuf, recvbuf, recvcounts, smpi_type_f2c(*datatype),
 517                         smpi_op_f2c(*op), smpi_comm_f2c(*comm));
 518 }
 519
 520 void mpi_scatter_(void* sendbuf, int* sendcount, int* sendtype, void* recvbuf, int* recvcount, int* recvtype,
 521                    int* root, int* comm, int* ierr) {
 522   recvbuf = static_cast<char *>( FORT_IN_PLACE(recvbuf));
 523   *ierr = MPI_Scatter(sendbuf, *sendcount, smpi_type_f2c(*sendtype),
 524                       recvbuf, *recvcount, smpi_type_f2c(*recvtype), *root, smpi_comm_f2c(*comm));
 525 }
 526
 527 void mpi_scatterv_(void* sendbuf, int* sendcounts, int* displs, int* sendtype,
 528                    void* recvbuf, int* recvcount, int* recvtype, int* root, int* comm, int* ierr) {
 529   recvbuf = static_cast<char *>( FORT_IN_PLACE(recvbuf));
 530   *ierr = MPI_Scatterv(sendbuf, sendcounts, displs, smpi_type_f2c(*sendtype),
 531                       recvbuf, *recvcount, smpi_type_f2c(*recvtype), *root, smpi_comm_f2c(*comm));
 532 }
 533
 534 void mpi_gather_(void* sendbuf, int* sendcount, int* sendtype, void* recvbuf, int* recvcount, int* recvtype,
 535                   int* root, int* comm, int* ierr) {
 536   sendbuf = static_cast<char *>( FORT_IN_PLACE(sendbuf));
 537   sendbuf = static_cast<char *>( FORT_BOTTOM(sendbuf));
 538   recvbuf = static_cast<char *>( FORT_BOTTOM(recvbuf));
 539   *ierr = MPI_Gather(sendbuf, *sendcount, smpi_type_f2c(*sendtype),
 540                      recvbuf, *recvcount, smpi_type_f2c(*recvtype), *root, smpi_comm_f2c(*comm));
 541 }
 542
 543 void mpi_gatherv_(void* sendbuf, int* sendcount, int* sendtype,
 544                   void* recvbuf, int* recvcounts, int* displs, int* recvtype, int* root, int* comm, int* ierr) {
 545   sendbuf = static_cast<char *>( FORT_IN_PLACE(sendbuf));
 546   sendbuf = static_cast<char *>( FORT_BOTTOM(sendbuf));
 547   recvbuf = static_cast<char *>( FORT_BOTTOM(recvbuf));
 548   *ierr = MPI_Gatherv(sendbuf, *sendcount, smpi_type_f2c(*sendtype),
 549                      recvbuf, recvcounts, displs, smpi_type_f2c(*recvtype), *root, smpi_comm_f2c(*comm));
 550 }
 551
 552 void mpi_allgather_(void* sendbuf, int* sendcount, int* sendtype, void* recvbuf, int* recvcount, int* recvtype,
 553                      int* comm, int* ierr) {
 554   sendbuf = static_cast<char *>( FORT_IN_PLACE(sendbuf));
 555   *ierr = MPI_Allgather(sendbuf, *sendcount, smpi_type_f2c(*sendtype),
 556                         recvbuf, *recvcount, smpi_type_f2c(*recvtype), smpi_comm_f2c(*comm));
 557 }
 558
 559 void mpi_allgatherv_(void* sendbuf, int* sendcount, int* sendtype,
 560                      void* recvbuf, int* recvcounts,int* displs, int* recvtype, int* comm, int* ierr) {
 561   sendbuf = static_cast<char *>( FORT_IN_PLACE(sendbuf));
 562   *ierr = MPI_Allgatherv(sendbuf, *sendcount, smpi_type_f2c(*sendtype),
 563                         recvbuf, recvcounts, displs, smpi_type_f2c(*recvtype), smpi_comm_f2c(*comm));
 564 }
 565
 566 void mpi_scan_(void* sendbuf, void* recvbuf, int* count, int* datatype, int* op, int* comm, int* ierr) {
 567   *ierr = MPI_Scan(sendbuf, recvbuf, *count, smpi_type_f2c(*datatype),
 568                    smpi_op_f2c(*op), smpi_comm_f2c(*comm));
 569 }
 570
 571 void mpi_alltoall_(void* sendbuf, int* sendcount, int* sendtype,
 572                     void* recvbuf, int* recvcount, int* recvtype, int* comm, int* ierr) {
 573   *ierr = MPI_Alltoall(sendbuf, *sendcount, smpi_type_f2c(*sendtype),
 574                        recvbuf, *recvcount, smpi_type_f2c(*recvtype), smpi_comm_f2c(*comm));
 575 }
 576
 577 void mpi_alltoallv_(void* sendbuf, int* sendcounts, int* senddisps, int* sendtype,
 578                     void* recvbuf, int* recvcounts, int* recvdisps, int* recvtype, int* comm, int* ierr) {
 579   *ierr = MPI_Alltoallv(sendbuf, sendcounts, senddisps, smpi_type_f2c(*sendtype),
 580                        recvbuf, recvcounts, recvdisps, smpi_type_f2c(*recvtype), smpi_comm_f2c(*comm));
 581 }
 582
 583 void mpi_test_ (int * request, int *flag, MPI_Status * status, int* ierr){
 584   MPI_Request req = smpi_request_f2c(*request);
 585   *ierr= MPI_Test(&req, flag, FORT_STATUS_IGNORE(status));
 586   if(req==MPI_REQUEST_NULL){
 587       free_request(*request);
 588       *request=MPI_FORTRAN_REQUEST_NULL;
 589   }
 590 }
 591
 592 void mpi_testall_ (int* count, int * requests,  int *flag, MPI_Status * statuses, int* ierr){
 593   int i;
 594   MPI_Request* reqs = xbt_new(MPI_Request, *count);
 595   for(i = 0; i < *count; i++) {
 596     reqs[i] = smpi_request_f2c(requests[i]);
 597   }
 598   *ierr= MPI_Testall(*count, reqs, flag, FORT_STATUSES_IGNORE(statuses));
 599   for(i = 0; i < *count; i++) {
 600     if(reqs[i]==MPI_REQUEST_NULL){
 601         free_request(requests[i]);
 602         requests[i]=MPI_FORTRAN_REQUEST_NULL;
 603     }
 604   }
 605   xbt_free(reqs);
 606 }
 607
 608 void mpi_get_processor_name_(char *name, int *resultlen, int* ierr){
 609   *ierr = MPI_Get_processor_name(name, resultlen);
 610 }
 611
 612 void mpi_get_count_(MPI_Status * status, int* datatype, int *count, int* ierr){
 613   *ierr = MPI_Get_count(FORT_STATUS_IGNORE(status), smpi_type_f2c(*datatype), count);
 614 }
 615
 616 void mpi_attr_get_(int* comm, int* keyval, void* attr_value, int* flag, int* ierr ){
 617   *ierr = MPI_Attr_get(smpi_comm_f2c(*comm), *keyval, attr_value, flag);
 618 }
 619
 620 void mpi_type_extent_(int* datatype, MPI_Aint * extent, int* ierr){
 621   *ierr= MPI_Type_extent(smpi_type_f2c(*datatype),  extent);
 622 }
 623
 624 void mpi_type_commit_(int* datatype,  int* ierr){
 625   MPI_Datatype tmp= smpi_type_f2c(*datatype);
 626   *ierr= MPI_Type_commit(&tmp);
 627 }
 628
 629 void mpi_type_vector_(int* count, int* blocklen, int* stride, int* old_type, int* newtype,  int* ierr){
 630   MPI_Datatype tmp;
 631   *ierr= MPI_Type_vector(*count, *blocklen, *stride, smpi_type_f2c(*old_type), &tmp);
 632   if(*ierr == MPI_SUCCESS) {
 633     *newtype = smpi_type_add_f(tmp);
 634   }
 635 }
 636
 637 void mpi_type_create_vector_(int* count, int* blocklen, int* stride, int* old_type, int* newtype,  int* ierr){
 638   MPI_Datatype tmp;
 639   *ierr= MPI_Type_vector(*count, *blocklen, *stride, smpi_type_f2c(*old_type), &tmp);
 640   if(*ierr == MPI_SUCCESS) {
 641     *newtype = smpi_type_add_f(tmp);
 642   }
 643 }
 644
 645 void mpi_type_hvector_(int* count, int* blocklen, MPI_Aint* stride, int* old_type, int* newtype,  int* ierr){
 646   MPI_Datatype tmp;
 647   *ierr= MPI_Type_hvector (*count, *blocklen, *stride, smpi_type_f2c(*old_type), &tmp);
 648   if(*ierr == MPI_SUCCESS) {
 649     *newtype = smpi_type_add_f(tmp);
 650   }
 651 }
 652
 653 void mpi_type_create_hvector_(int* count, int* blocklen, MPI_Aint* stride, int* old_type, int* newtype,  int* ierr){
 654   MPI_Datatype tmp;
 655   *ierr= MPI_Type_hvector(*count, *blocklen, *stride, smpi_type_f2c(*old_type), &tmp);
 656   if(*ierr == MPI_SUCCESS) {
 657     *newtype = smpi_type_add_f(tmp);
 658   }
 659 }
 660
 661 void mpi_type_free_(int* datatype, int* ierr){
 662   MPI_Datatype tmp= smpi_type_f2c(*datatype);
 663   *ierr= MPI_Type_free (&tmp);
 664   if(*ierr == MPI_SUCCESS) {
 665     free_type(*datatype);
 666   }
 667 }
 668
 669 void mpi_type_ub_(int* datatype, MPI_Aint * disp, int* ierr){
 670   *ierr= MPI_Type_ub(smpi_type_f2c(*datatype), disp);
 671 }
 672
 673 void mpi_type_lb_(int* datatype, MPI_Aint * extent, int* ierr){
 674   *ierr= MPI_Type_extent(smpi_type_f2c(*datatype), extent);
 675 }
 676
 677 void mpi_type_size_(int* datatype, int *size, int* ierr)
 678 {
 679   *ierr = MPI_Type_size(smpi_type_f2c(*datatype), size);
 680 }
 681
 682 void mpi_error_string_(int* errorcode, char* string, int* resultlen, int* ierr){
 683   *ierr = MPI_Error_string(*errorcode, string, resultlen);
 684 }
 685
 686 void mpi_win_fence_( int* assert,  int* win, int* ierr){
 687   *ierr =  MPI_Win_fence(* assert, smpi_win_f2c(*win));
 688 }
 689
 690 void mpi_win_free_( int* win, int* ierr){
 691   MPI_Win tmp = smpi_win_f2c(*win);
 692   *ierr =  MPI_Win_free(&tmp);
 693   if(*ierr == MPI_SUCCESS) {
 694     free_win(*win);
 695   }
 696 }
 697
 698 void mpi_win_create_( int *base, MPI_Aint* size, int* disp_unit, int* info, int* comm, int *win, int* ierr){
 699   MPI_Win tmp;
 700   *ierr =  MPI_Win_create( static_cast<void*>(base), *size, *disp_unit, smpi_info_f2c(*info), smpi_comm_f2c(*comm),&tmp);
 701  if(*ierr == MPI_SUCCESS) {
 702    *win = smpi_win_add_f(tmp);
 703  }
 704 }
 705
 706 void mpi_win_post_(int* group, int assert, int* win, int* ierr){
 707   *ierr =  MPI_Win_post(smpi_group_f2c(*group), assert, smpi_win_f2c(*win));
 708 }
 709
 710 void mpi_win_start_(int* group, int assert, int* win, int* ierr){
 711   *ierr =  MPI_Win_start(smpi_group_f2c(*group), assert, smpi_win_f2c(*win));
 712 }
 713
 714 void mpi_win_complete_(int* win, int* ierr){
 715   *ierr =  MPI_Win_complete(smpi_win_f2c(*win));
 716 }
 717
 718 void mpi_win_wait_(int* win, int* ierr){
 719   *ierr =  MPI_Win_wait(smpi_win_f2c(*win));
 720 }
 721
 722 void mpi_win_set_name_ (int*  win, char * name, int* ierr, int size){
 723   //handle trailing blanks
 724   while(name[size-1]==' ')
 725     size--;
 726   while(*name==' '){//handle leading blanks
 727     size--;
 728     name++;
 729   }
 730   char* tname = xbt_new(char,size+1);
 731   strncpy(tname, name, size);
 732   tname[size]='\0';
 733   *ierr = MPI_Win_set_name(smpi_win_f2c(*win), tname);
 734   xbt_free(tname);
 735 }
 736
 737 void mpi_win_get_name_ (int*  win, char * name, int* len, int* ierr){
 738   *ierr = MPI_Win_get_name(smpi_win_f2c(*win),name,len);
 739   if(*len>0)
 740     name[*len]=' ';//blank padding, not \0
 741 }
 742
 743 void mpi_info_create_( int *info, int* ierr){
 744   MPI_Info tmp;
 745   *ierr =  MPI_Info_create(&tmp);
 746   if(*ierr == MPI_SUCCESS) {
 747     *info = smpi_info_add_f(tmp);
 748   }
 749 }
 750
 751 void mpi_info_set_( int *info, char *key, char *value, int* ierr, unsigned int keylen, unsigned int valuelen){
 752   //handle trailing blanks
 753   while(key[keylen-1]==' ')
 754     keylen--;
 755   while(*key==' '){//handle leading blanks
 756     keylen--;
 757     key++;
 758   }
 759   char* tkey = xbt_new(char,keylen+1);
 760   strncpy(tkey, key, keylen);
 761   tkey[keylen]='\0';
 762
 763   while(value[valuelen-1]==' ')
 764     valuelen--;
 765   while(*value==' '){//handle leading blanks
 766     valuelen--;
 767     value++;
 768   }
 769   char* tvalue = xbt_new(char,valuelen+1);
 770   strncpy(tvalue, value, valuelen);
 771   tvalue[valuelen]='\0';
 772
 773   *ierr =  MPI_Info_set( smpi_info_f2c(*info), tkey, tvalue);
 774   xbt_free(tkey);
 775   xbt_free(tvalue);
 776 }
 777
 778 void mpi_info_get_ (int* info,char *key,int* valuelen, char *value, int *flag, int* ierr, unsigned int keylen ){
 779   while(key[keylen-1]==' ')
 780     keylen--;
 781   while(*key==' '){//handle leading blanks
 782     keylen--;
 783     key++;
 784   }
 785   char* tkey = xbt_new(char,keylen+1);
 786   strncpy(tkey, key, keylen);
 787   tkey[keylen]='\0';
 788   *ierr = MPI_Info_get(smpi_info_f2c(*info),tkey,*valuelen, value, flag);
 789   xbt_free(tkey);
 790   if(*flag!=0){
 791     int replace=0;
 792     int i=0;
 793     for (i=0; i<*valuelen; i++){
 794       if(value[i]=='\0')
 795         replace=1;
 796       if(replace)
 797         value[i]=' ';
 798     }
 799   }
 800 }
 801
 802 void mpi_info_free_(int* info, int* ierr){
 803   MPI_Info tmp = smpi_info_f2c(*info);
 804   *ierr =  MPI_Info_free(&tmp);
 805   if(*ierr == MPI_SUCCESS) {
 806     free_info(*info);
 807   }
 808 }
 809
 810 void mpi_get_( int *origin_addr, int* origin_count, int* origin_datatype, int *target_rank,
 811     MPI_Aint* target_disp, int *target_count, int* tarsmpi_type_f2c, int* win, int* ierr){
 812   *ierr =  MPI_Get( static_cast<void*>(origin_addr),*origin_count, smpi_type_f2c(*origin_datatype),*target_rank,
 813       *target_disp, *target_count,smpi_type_f2c(*tarsmpi_type_f2c), smpi_win_f2c(*win));
 814 }
 815
 816 void mpi_accumulate_( int *origin_addr, int* origin_count, int* origin_datatype, int *target_rank,
 817     MPI_Aint* target_disp, int *target_count, int* tarsmpi_type_f2c, int* op, int* win, int* ierr){
 818   *ierr =  MPI_Accumulate( static_cast<void *>(origin_addr),*origin_count, smpi_type_f2c(*origin_datatype),*target_rank,
 819       *target_disp, *target_count,smpi_type_f2c(*tarsmpi_type_f2c), smpi_op_f2c(*op), smpi_win_f2c(*win));
 820 }
 821
 822 void mpi_put_( int *origin_addr, int* origin_count, int* origin_datatype, int *target_rank,
 823     MPI_Aint* target_disp, int *target_count, int* tarsmpi_type_f2c, int* win, int* ierr){
 824   *ierr =  MPI_Put( static_cast<void *>(origin_addr),*origin_count, smpi_type_f2c(*origin_datatype),*target_rank,
 825       *target_disp, *target_count,smpi_type_f2c(*tarsmpi_type_f2c), smpi_win_f2c(*win));
 826 }
 827
 828 //following are automatically generated, and have to be checked
 829 void mpi_finalized_ (int * flag, int* ierr){
 830
 831  *ierr = MPI_Finalized(flag);
 832 }
 833
 834 void mpi_init_thread_ (int* required, int *provided, int* ierr){
 835   smpi_init_fortran_types();
 836   *ierr = MPI_Init_thread(nullptr, nullptr,*required, provided);
 837   running_processes++;
 838 }
 839
 840 void mpi_query_thread_ (int *provided, int* ierr){
 841
 842  *ierr = MPI_Query_thread(provided);
 843 }
 844
 845 void mpi_is_thread_main_ (int *flag, int* ierr){
 846
 847  *ierr = MPI_Is_thread_main(flag);
 848 }
 849
 850 void mpi_address_ (void *location, MPI_Aint * address, int* ierr){
 851
 852  *ierr = MPI_Address(location, address);
 853 }
 854
 855 void mpi_get_address_ (void *location, MPI_Aint * address, int* ierr){
 856
 857  *ierr = MPI_Get_address(location, address);
 858 }
 859
 860 void mpi_type_dup_ (int*  datatype, int* newdatatype, int* ierr){
 861  MPI_Datatype tmp;
 862  *ierr = MPI_Type_dup(smpi_type_f2c(*datatype), &tmp);
 863  if(*ierr == MPI_SUCCESS) {
 864    *newdatatype = smpi_type_add_f(tmp);
 865  }
 866 }
 867
 868 void mpi_type_set_name_ (int*  datatype, char * name, int* ierr, int size){
 869  char* tname = xbt_new(char, size+1);
 870  strncpy(tname, name, size);
 871  tname[size]='\0';
 872  *ierr = MPI_Type_set_name(smpi_type_f2c(*datatype), tname);
 873  xbt_free(tname);
 874 }
 875
 876 void mpi_type_get_name_ (int*  datatype, char * name, int* len, int* ierr){
 877  *ierr = MPI_Type_get_name(smpi_type_f2c(*datatype),name,len);
 878   if(*len>0)
 879     name[*len]=' ';
 880 }
 881
 882 void mpi_type_get_attr_ (int* type, int* type_keyval, void *attribute_val, int* flag, int* ierr){
 883
 884  *ierr = MPI_Type_get_attr ( smpi_type_f2c(*type), *type_keyval, attribute_val,flag);
 885 }
 886
 887 void mpi_type_set_attr_ (int* type, int* type_keyval, void *attribute_val, int* ierr){
 888
 889  *ierr = MPI_Type_set_attr ( smpi_type_f2c(*type), *type_keyval, attribute_val);
 890 }
 891
 892 void mpi_type_delete_attr_ (int* type, int* type_keyval, int* ierr){
 893
 894  *ierr = MPI_Type_delete_attr ( smpi_type_f2c(*type),  *type_keyval);
 895 }
 896
 897 void mpi_type_create_keyval_ (void* copy_fn, void*  delete_fn, int* keyval, void* extra_state, int* ierr){
 898
 899  *ierr = MPI_Type_create_keyval(reinterpret_cast<MPI_Type_copy_attr_function*>(copy_fn), reinterpret_cast<MPI_Type_delete_attr_function*>(delete_fn),
 900                                 keyval,  extra_state) ;
 901 }
 902
 903 void mpi_type_free_keyval_ (int* keyval, int* ierr) {
 904  *ierr = MPI_Type_free_keyval( keyval);
 905 }
 906
 907 void mpi_pcontrol_ (int* level , int* ierr){
 908  *ierr = MPI_Pcontrol(*static_cast<const int*>(level));
 909 }
 910
 911 void mpi_type_get_extent_ (int* datatype, MPI_Aint * lb, MPI_Aint * extent, int* ierr){
 912
 913  *ierr = MPI_Type_get_extent(smpi_type_f2c(*datatype), lb, extent);
 914 }
 915
 916 void mpi_type_get_true_extent_ (int* datatype, MPI_Aint * lb, MPI_Aint * extent, int* ierr){
 917
 918  *ierr = MPI_Type_get_true_extent(smpi_type_f2c(*datatype), lb, extent);
 919 }
 920
 921 void mpi_op_create_ (void * function, int* commute, int* op, int* ierr){
 922   MPI_Op tmp;
 923  *ierr = MPI_Op_create(reinterpret_cast<MPI_User_function*>(function),*commute, &tmp);
 924  if(*ierr == MPI_SUCCESS) {
 925    smpi_op_set_fortran(tmp);
 926    *op = smpi_op_add_f(tmp);
 927  }
 928 }
 929
 930 void mpi_op_free_ (int* op, int* ierr){
 931   MPI_Op tmp=smpi_op_f2c(*op);
 932   *ierr = MPI_Op_free(& tmp);
 933   if(*ierr == MPI_SUCCESS) {
 934     free_op(*op);
 935   }
 936 }
 937
 938 void mpi_group_free_ (int* group, int* ierr){
 939  MPI_Group tmp=smpi_group_f2c(*group);
 940  *ierr = MPI_Group_free(&tmp);
 941  if(*ierr == MPI_SUCCESS) {
 942    free_group(*group);
 943  }
 944 }
 945
 946 void mpi_group_size_ (int* group, int *size, int* ierr){
 947
 948  *ierr = MPI_Group_size(smpi_group_f2c(*group), size);
 949 }
 950
 951 void mpi_group_rank_ (int* group, int *rank, int* ierr){
 952
 953  *ierr = MPI_Group_rank(smpi_group_f2c(*group), rank);
 954 }
 955
 956 void mpi_group_translate_ranks_ (int* group1, int* n, int *ranks1, int* group2, int *ranks2, int* ierr)
 957 {
 958
 959  *ierr = MPI_Group_translate_ranks(smpi_group_f2c(*group1), *n, ranks1, smpi_group_f2c(*group2), ranks2);
 960 }
 961
 962 void mpi_group_compare_ (int* group1, int* group2, int *result, int* ierr){
 963
 964  *ierr = MPI_Group_compare(smpi_group_f2c(*group1), smpi_group_f2c(*group2), result);
 965 }
 966
 967 void mpi_group_union_ (int* group1, int* group2, int* newgroup, int* ierr){
 968  MPI_Group tmp;
 969  *ierr = MPI_Group_union(smpi_group_f2c(*group1), smpi_group_f2c(*group2), &tmp);
 970  if(*ierr == MPI_SUCCESS) {
 971    *newgroup = smpi_group_add_f(tmp);
 972  }
 973 }
 974
 975 void mpi_group_intersection_ (int* group1, int* group2, int* newgroup, int* ierr){
 976  MPI_Group tmp;
 977  *ierr = MPI_Group_intersection(smpi_group_f2c(*group1), smpi_group_f2c(*group2), &tmp);
 978  if(*ierr == MPI_SUCCESS) {
 979    *newgroup = smpi_group_add_f(tmp);
 980  }
 981 }
 982
 983 void mpi_group_difference_ (int* group1, int* group2, int* newgroup, int* ierr){
 984  MPI_Group tmp;
 985  *ierr = MPI_Group_difference(smpi_group_f2c(*group1), smpi_group_f2c(*group2), &tmp);
 986  if(*ierr == MPI_SUCCESS) {
 987    *newgroup = smpi_group_add_f(tmp);
 988  }
 989 }
 990
 991 void mpi_group_excl_ (int* group, int* n, int *ranks, int* newgroup, int* ierr){
 992   MPI_Group tmp;
 993  *ierr = MPI_Group_excl(smpi_group_f2c(*group), *n, ranks, &tmp);
 994  if(*ierr == MPI_SUCCESS) {
 995    *newgroup = smpi_group_add_f(tmp);
 996  }
 997 }
 998
 999 void mpi_group_range_incl_ (int* group, int* n, int ranges[][3], int* newgroup, int* ierr)
1000 {
1001   MPI_Group tmp;
1002  *ierr = MPI_Group_range_incl(smpi_group_f2c(*group), *n, ranges, &tmp);
1003  if(*ierr == MPI_SUCCESS) {
1004    *newgroup = smpi_group_add_f(tmp);
1005  }
1006 }
1007
1008 void mpi_group_range_excl_ (int* group, int* n, int ranges[][3], int* newgroup, int* ierr)
1009 {
1010  MPI_Group tmp;
1011  *ierr = MPI_Group_range_excl(smpi_group_f2c(*group), *n, ranges, &tmp);
1012  if(*ierr == MPI_SUCCESS) {
1013    *newgroup = smpi_group_add_f(tmp);
1014  }
1015 }
1016
1017 void mpi_comm_get_attr_ (int* comm, int* comm_keyval, void *attribute_val, int *flag, int* ierr){
1018
1019  *ierr = MPI_Comm_get_attr (smpi_comm_f2c(*comm), *comm_keyval, attribute_val, flag);
1020 }
1021
1022 void mpi_comm_set_attr_ (int* comm, int* comm_keyval, void *attribute_val, int* ierr){
1023
1024  *ierr = MPI_Comm_set_attr ( smpi_comm_f2c(*comm), *comm_keyval, attribute_val);
1025 }
1026
1027 void mpi_comm_delete_attr_ (int* comm, int* comm_keyval, int* ierr){
1028
1029  *ierr = MPI_Comm_delete_attr (smpi_comm_f2c(*comm),  *comm_keyval);
1030 }
1031
1032 void mpi_comm_create_keyval_ (void* copy_fn, void* delete_fn, int* keyval, void* extra_state, int* ierr){
1033
1034  *ierr = MPI_Comm_create_keyval(reinterpret_cast<MPI_Comm_copy_attr_function*>(copy_fn),  reinterpret_cast<MPI_Comm_delete_attr_function*>(delete_fn),
1035          keyval,  extra_state) ;
1036 }
1037
1038 void mpi_comm_free_keyval_ (int* keyval, int* ierr) {
1039  *ierr = MPI_Comm_free_keyval( keyval);
1040 }
1041
1042 void mpi_comm_get_name_ (int* comm, char* name, int* len, int* ierr){
1043  *ierr = MPI_Comm_get_name(smpi_comm_f2c(*comm), name, len);
1044   if(*len>0)
1045     name[*len]=' ';
1046 }
1047
1048 void mpi_comm_compare_ (int* comm1, int* comm2, int *result, int* ierr){
1049
1050  *ierr = MPI_Comm_compare(smpi_comm_f2c(*comm1), smpi_comm_f2c(*comm2), result);
1051 }
1052
1053 void mpi_comm_disconnect_ (int* comm, int* ierr){
1054  MPI_Comm tmp=smpi_comm_f2c(*comm);
1055  *ierr = MPI_Comm_disconnect(&tmp);
1056  if(*ierr == MPI_SUCCESS) {
1057    free_comm(*comm);
1058  }
1059 }
1060
1061 void mpi_request_free_ (int* request, int* ierr){
1062   MPI_Request tmp=smpi_request_f2c(*request);
1063  *ierr = MPI_Request_free(&tmp);
1064  if(*ierr == MPI_SUCCESS) {
1065    free_request(*request);
1066  }
1067 }
1068
1069 void mpi_sendrecv_replace_ (void *buf, int* count, int* datatype, int* dst, int* sendtag, int* src, int* recvtag,
1070                             int* comm, MPI_Status* status, int* ierr)
1071 {
1072  *ierr = MPI_Sendrecv_replace(buf, *count, smpi_type_f2c(*datatype), *dst, *sendtag, *src,
1073  *recvtag, smpi_comm_f2c(*comm), FORT_STATUS_IGNORE(status));
1074 }
1075
1076 void mpi_testany_ (int* count, int* requests, int *index, int *flag, MPI_Status* status, int* ierr)
1077 {
1078   MPI_Request* reqs;
1079   int i;
1080
1081   reqs = xbt_new(MPI_Request, *count);
1082   for(i = 0; i < *count; i++) {
1083     reqs[i] = smpi_request_f2c(requests[i]);
1084   }
1085   *ierr = MPI_Testany(*count, reqs, index, flag, FORT_STATUS_IGNORE(status));
1086   if(*index!=MPI_UNDEFINED && reqs[*index]==MPI_REQUEST_NULL){
1087     free_request(requests[*index]);
1088     requests[*index]=MPI_FORTRAN_REQUEST_NULL;
1089   }
1090   xbt_free(reqs);
1091 }
1092
1093 void mpi_waitsome_ (int* incount, int* requests, int *outcount, int *indices, MPI_Status* status, int* ierr)
1094 {
1095   MPI_Request* reqs;
1096   int i;
1097
1098   reqs = xbt_new(MPI_Request, *incount);
1099   for(i = 0; i < *incount; i++) {
1100     reqs[i] = smpi_request_f2c(requests[i]);
1101   }
1102   *ierr = MPI_Waitsome(*incount, reqs, outcount, indices, status);
1103   for(i=0;i<*outcount;i++){
1104     if(reqs[indices[i]]==MPI_REQUEST_NULL){
1105         free_request(requests[indices[i]]);
1106         requests[indices[i]]=MPI_FORTRAN_REQUEST_NULL;
1107     }
1108   }
1109   xbt_free(reqs);
1110 }
1111
1112 void mpi_reduce_local_ (void *inbuf, void *inoutbuf, int* count, int* datatype, int* op, int* ierr){
1113
1114  *ierr = MPI_Reduce_local(inbuf, inoutbuf, *count, smpi_type_f2c(*datatype), smpi_op_f2c(*op));
1115 }
1116
1117 void mpi_reduce_scatter_block_ (void *sendbuf, void *recvbuf, int* recvcount, int* datatype, int* op, int* comm,
1118                                 int* ierr)
1119 {
1120   sendbuf = static_cast<char *>( FORT_IN_PLACE(sendbuf));
1121  *ierr = MPI_Reduce_scatter_block(sendbuf, recvbuf, *recvcount, smpi_type_f2c(*datatype), smpi_op_f2c(*op),
1122                                   smpi_comm_f2c(*comm));
1123 }
1124
1125 void mpi_pack_size_ (int* incount, int* datatype, int* comm, int* size, int* ierr) {
1126  *ierr = MPI_Pack_size(*incount, smpi_type_f2c(*datatype), smpi_comm_f2c(*comm), size);
1127 }
1128
1129 void mpi_cart_coords_ (int* comm, int* rank, int* maxdims, int* coords, int* ierr) {
1130  *ierr = MPI_Cart_coords(smpi_comm_f2c(*comm), *rank, *maxdims, coords);
1131 }
1132
1133 void mpi_cart_create_ (int* comm_old, int* ndims, int* dims, int* periods, int* reorder, int*  comm_cart, int* ierr) {
1134   MPI_Comm tmp;
1135  *ierr = MPI_Cart_create(smpi_comm_f2c(*comm_old), *ndims, dims, periods, *reorder, &tmp);
1136  if(*ierr == MPI_SUCCESS) {
1137    *comm_cart = smpi_comm_add_f(tmp);
1138  }
1139 }
1140
1141 void mpi_cart_get_ (int* comm, int* maxdims, int* dims, int* periods, int* coords, int* ierr) {
1142  *ierr = MPI_Cart_get(smpi_comm_f2c(*comm), *maxdims, dims, periods, coords);
1143 }
1144
1145 void mpi_cart_map_ (int* comm_old, int* ndims, int* dims, int* periods, int* newrank, int* ierr) {
1146  *ierr = MPI_Cart_map(smpi_comm_f2c(*comm_old), *ndims, dims, periods, newrank);
1147 }
1148
1149 void mpi_cart_rank_ (int* comm, int* coords, int* rank, int* ierr) {
1150  *ierr = MPI_Cart_rank(smpi_comm_f2c(*comm), coords, rank);
1151 }
1152
1153 void mpi_cart_shift_ (int* comm, int* direction, int* displ, int* source, int* dest, int* ierr) {
1154  *ierr = MPI_Cart_shift(smpi_comm_f2c(*comm), *direction, *displ, source, dest);
1155 }
1156
1157 void mpi_cart_sub_ (int* comm, int* remain_dims, int*  comm_new, int* ierr) {
1158  MPI_Comm tmp;
1159  *ierr = MPI_Cart_sub(smpi_comm_f2c(*comm), remain_dims, &tmp);
1160  if(*ierr == MPI_SUCCESS) {
1161    *comm_new = smpi_comm_add_f(tmp);
1162  }
1163 }
1164
1165 void mpi_cartdim_get_ (int* comm, int* ndims, int* ierr) {
1166  *ierr = MPI_Cartdim_get(smpi_comm_f2c(*comm), ndims);
1167 }
1168
1169 void mpi_graph_create_ (int* comm_old, int* nnodes, int* index, int* edges, int* reorder, int*  comm_graph, int* ierr) {
1170   MPI_Comm tmp;
1171  *ierr = MPI_Graph_create(smpi_comm_f2c(*comm_old), *nnodes, index, edges, *reorder, &tmp);
1172  if(*ierr == MPI_SUCCESS) {
1173    *comm_graph = smpi_comm_add_f(tmp);
1174  }
1175 }
1176
1177 void mpi_graph_get_ (int* comm, int* maxindex, int* maxedges, int* index, int* edges, int* ierr) {
1178  *ierr = MPI_Graph_get(smpi_comm_f2c(*comm), *maxindex, *maxedges, index, edges);
1179 }
1180
1181 void mpi_graph_map_ (int* comm_old, int* nnodes, int* index, int* edges, int* newrank, int* ierr) {
1182  *ierr = MPI_Graph_map(smpi_comm_f2c(*comm_old), *nnodes, index, edges, newrank);
1183 }
1184
1185 void mpi_graph_neighbors_ (int* comm, int* rank, int* maxneighbors, int* neighbors, int* ierr) {
1186  *ierr = MPI_Graph_neighbors(smpi_comm_f2c(*comm), *rank, *maxneighbors, neighbors);
1187 }
1188
1189 void mpi_graph_neighbors_count_ (int* comm, int* rank, int* nneighbors, int* ierr) {
1190  *ierr = MPI_Graph_neighbors_count(smpi_comm_f2c(*comm), *rank, nneighbors);
1191 }
1192
1193 void mpi_graphdims_get_ (int* comm, int* nnodes, int* nedges, int* ierr) {
1194  *ierr = MPI_Graphdims_get(smpi_comm_f2c(*comm), nnodes, nedges);
1195 }
1196
1197 void mpi_topo_test_ (int* comm, int* top_type, int* ierr) {
1198  *ierr = MPI_Topo_test(smpi_comm_f2c(*comm), top_type);
1199 }
1200
1201 void mpi_error_class_ (int* errorcode, int* errorclass, int* ierr) {
1202  *ierr = MPI_Error_class(*errorcode, errorclass);
1203 }
1204
1205 void mpi_errhandler_create_ (void* function, void* errhandler, int* ierr) {
1206  *ierr = MPI_Errhandler_create(reinterpret_cast<MPI_Handler_function*>(function), static_cast<MPI_Errhandler*>(errhandler));
1207 }
1208
1209 void mpi_errhandler_free_ (void* errhandler, int* ierr) {
1210  *ierr = MPI_Errhandler_free(static_cast<MPI_Errhandler*>(errhandler));
1211 }
1212
1213 void mpi_errhandler_get_ (int* comm, void* errhandler, int* ierr) {
1214  *ierr = MPI_Errhandler_get(smpi_comm_f2c(*comm), static_cast<MPI_Errhandler*>(errhandler));
1215 }
1216
1217 void mpi_errhandler_set_ (int* comm, void* errhandler, int* ierr) {
1218  *ierr = MPI_Errhandler_set(smpi_comm_f2c(*comm), *static_cast<MPI_Errhandler*>(errhandler));
1219 }
1220
1221 void mpi_comm_set_errhandler_ (int* comm, void* errhandler, int* ierr) {
1222  *ierr = MPI_Errhandler_set(smpi_comm_f2c(*comm), *static_cast<MPI_Errhandler*>(errhandler));
1223 }
1224
1225 void mpi_comm_get_errhandler_ (int* comm, void* errhandler, int* ierr) {
1226  *ierr = MPI_Errhandler_set(smpi_comm_f2c(*comm), static_cast<MPI_Errhandler*>(errhandler));
1227 }
1228
1229 void mpi_type_contiguous_ (int* count, int* old_type, int*  newtype, int* ierr) {
1230   MPI_Datatype tmp;
1231   *ierr = MPI_Type_contiguous(*count, smpi_type_f2c(*old_type), &tmp);
1232   if(*ierr == MPI_SUCCESS) {
1233     *newtype = smpi_type_add_f(tmp);
1234   }
1235 }
1236
1237 void mpi_cancel_ (int* request, int* ierr) {
1238   MPI_Request tmp=smpi_request_f2c(*request);
1239  *ierr = MPI_Cancel(&tmp);
1240 }
1241
1242 void mpi_buffer_attach_ (void* buffer, int* size, int* ierr) {
1243  *ierr = MPI_Buffer_attach(buffer, *size);
1244 }
1245
1246 void mpi_buffer_detach_ (void* buffer, int* size, int* ierr) {
1247  *ierr = MPI_Buffer_detach(buffer, size);
1248 }
1249
1250 void mpi_testsome_ (int* incount, int*  requests, int* outcount, int* indices, MPI_Status*  statuses, int* ierr) {
1251   MPI_Request* reqs;
1252   int i;
1253
1254   reqs = xbt_new(MPI_Request, *incount);
1255   for(i = 0; i < *incount; i++) {
1256     reqs[i] = smpi_request_f2c(requests[i]);
1257     indices[i]=0;
1258   }
1259   *ierr = MPI_Testsome(*incount, reqs, outcount, indices, FORT_STATUSES_IGNORE(statuses));
1260   for(i=0;i<*incount;i++){
1261     if(indices[i] && reqs[indices[i]]==MPI_REQUEST_NULL){
1262       free_request(requests[indices[i]]);
1263       requests[indices[i]]=MPI_FORTRAN_REQUEST_NULL;
1264     }
1265   }
1266   xbt_free(reqs);
1267 }
1268
1269 void mpi_comm_test_inter_ (int* comm, int* flag, int* ierr) {
1270  *ierr = MPI_Comm_test_inter(smpi_comm_f2c(*comm), flag);
1271 }
1272
1273 void mpi_unpack_ (void* inbuf, int* insize, int* position, void* outbuf, int* outcount, int* type, int* comm,
1274                   int* ierr) {
1275  *ierr = MPI_Unpack(inbuf, *insize, position, outbuf, *outcount, smpi_type_f2c(*type), smpi_comm_f2c(*comm));
1276 }
1277
1278 void mpi_pack_external_size_ (char *datarep, int* incount, int* datatype, MPI_Aint *size, int* ierr){
1279  *ierr = MPI_Pack_external_size(datarep, *incount, smpi_type_f2c(*datatype), size);
1280 }
1281
1282 void mpi_pack_external_ (char *datarep, void *inbuf, int* incount, int* datatype, void *outbuf, MPI_Aint* outcount,
1283                          MPI_Aint *position, int* ierr){
1284  *ierr = MPI_Pack_external(datarep, inbuf, *incount, smpi_type_f2c(*datatype), outbuf, *outcount, position);
1285 }
1286
1287 void mpi_unpack_external_ ( char *datarep, void *inbuf, MPI_Aint* insize, MPI_Aint *position, void *outbuf,
1288                             int* outcount, int* datatype, int* ierr){
1289  *ierr = MPI_Unpack_external( datarep, inbuf, *insize, position, outbuf, *outcount, smpi_type_f2c(*datatype));
1290 }
1291
1292 void mpi_type_hindexed_ (int* count, int* blocklens, MPI_Aint* indices, int* old_type, int*  newtype, int* ierr) {
1293   MPI_Datatype tmp;
1294   *ierr = MPI_Type_hindexed(*count, blocklens, indices, smpi_type_f2c(*old_type), &tmp);
1295   if(*ierr == MPI_SUCCESS) {
1296     *newtype = smpi_type_add_f(tmp);
1297   }
1298 }
1299
1300 void mpi_type_create_hindexed_(int* count, int* blocklens, MPI_Aint* indices, int* old_type, int*  newtype, int* ierr){
1301   MPI_Datatype tmp;
1302   *ierr = MPI_Type_create_hindexed(*count, blocklens, indices, smpi_type_f2c(*old_type), &tmp);
1303   if(*ierr == MPI_SUCCESS) {
1304     *newtype = smpi_type_add_f(tmp);
1305   }
1306 }
1307
1308 void mpi_type_create_hindexed_block_ (int* count, int* blocklength, MPI_Aint* indices, int* old_type, int*  newtype,
1309                                       int* ierr) {
1310   MPI_Datatype tmp;
1311   *ierr = MPI_Type_create_hindexed_block(*count, *blocklength, indices, smpi_type_f2c(*old_type), &tmp);
1312   if(*ierr == MPI_SUCCESS) {
1313     *newtype = smpi_type_add_f(tmp);
1314   }
1315 }
1316
1317 void mpi_type_indexed_ (int* count, int* blocklens, int* indices, int* old_type, int*  newtype, int* ierr) {
1318   MPI_Datatype tmp;
1319   *ierr = MPI_Type_indexed(*count, blocklens, indices, smpi_type_f2c(*old_type), &tmp);
1320   if(*ierr == MPI_SUCCESS) {
1321     *newtype = smpi_type_add_f(tmp);
1322   }
1323 }
1324
1325 void mpi_type_create_indexed_block_ (int* count, int* blocklength, int* indices,  int* old_type,  int*newtype,
1326                                      int* ierr){
1327   MPI_Datatype tmp;
1328   *ierr = MPI_Type_create_indexed_block(*count, *blocklength, indices, smpi_type_f2c(*old_type), &tmp);
1329   if(*ierr == MPI_SUCCESS) {
1330     *newtype = smpi_type_add_f(tmp);
1331   }
1332 }
1333
1334 void mpi_type_struct_ (int* count, int* blocklens, MPI_Aint* indices, int* old_types, int*  newtype, int* ierr) {
1335   MPI_Datatype tmp;
1336   int i=0;
1337   MPI_Datatype* types = static_cast<MPI_Datatype*>(xbt_malloc(*count*sizeof(MPI_Datatype)));
1338   for(i=0; i< *count; i++){
1339     types[i] = smpi_type_f2c(old_types[i]);
1340   }
1341   *ierr = MPI_Type_struct(*count, blocklens, indices, types, &tmp);
1342   if(*ierr == MPI_SUCCESS) {
1343     *newtype = smpi_type_add_f(tmp);
1344   }
1345   xbt_free(types);
1346 }
1347
1348 void mpi_type_create_struct_(int* count, int* blocklens, MPI_Aint* indices, int*  old_types, int*  newtype, int* ierr){
1349   MPI_Datatype tmp;
1350   int i=0;
1351   MPI_Datatype* types = static_cast<MPI_Datatype*>(xbt_malloc(*count*sizeof(MPI_Datatype)));
1352   for(i=0; i< *count; i++){
1353     types[i] = smpi_type_f2c(old_types[i]);
1354   }
1355   *ierr = MPI_Type_create_struct(*count, blocklens, indices, types, &tmp);
1356   if(*ierr == MPI_SUCCESS) {
1357     *newtype = smpi_type_add_f(tmp);
1358   }
1359   xbt_free(types);
1360 }
1361
1362 void mpi_ssend_ (void* buf, int* count, int* datatype, int* dest, int* tag, int* comm, int* ierr) {
1363  *ierr = MPI_Ssend(buf, *count, smpi_type_f2c(*datatype), *dest, *tag, smpi_comm_f2c(*comm));
1364 }
1365
1366 void mpi_ssend_init_ (void* buf, int* count, int* datatype, int* dest, int* tag, int* comm, int* request, int* ierr) {
1367   MPI_Request tmp;
1368  *ierr = MPI_Ssend_init(buf, *count, smpi_type_f2c(*datatype), *dest, *tag, smpi_comm_f2c(*comm), &tmp);
1369  if(*ierr == MPI_SUCCESS) {
1370    *request = smpi_request_add_f(tmp);
1371  }
1372 }
1373
1374 void mpi_intercomm_create_ (int* local_comm, int *local_leader, int* peer_comm, int* remote_leader, int* tag,
1375                             int* comm_out, int* ierr) {
1376   MPI_Comm tmp;
1377   *ierr = MPI_Intercomm_create(smpi_comm_f2c(*local_comm), *local_leader,smpi_comm_f2c(*peer_comm), *remote_leader,
1378                                *tag, &tmp);
1379   if(*ierr == MPI_SUCCESS) {
1380     *comm_out = smpi_comm_add_f(tmp);
1381   }
1382 }
1383
1384 void mpi_intercomm_merge_ (int* comm, int* high, int*  comm_out, int* ierr) {
1385  MPI_Comm tmp;
1386  *ierr = MPI_Intercomm_merge(smpi_comm_f2c(*comm), *high, &tmp);
1387  if(*ierr == MPI_SUCCESS) {
1388    *comm_out = smpi_comm_add_f(tmp);
1389  }
1390 }
1391
1392 void mpi_bsend_ (void* buf, int* count, int* datatype, int *dest, int* tag, int* comm, int* ierr) {
1393  *ierr = MPI_Bsend(buf, *count, smpi_type_f2c(*datatype), *dest, *tag, smpi_comm_f2c(*comm));
1394 }
1395
1396 void mpi_bsend_init_ (void* buf, int* count, int* datatype, int *dest, int* tag, int* comm, int*  request, int* ierr) {
1397   MPI_Request tmp;
1398   *ierr = MPI_Bsend_init(buf, *count, smpi_type_f2c(*datatype), *dest, *tag, smpi_comm_f2c(*comm), &tmp);
1399  if(*ierr == MPI_SUCCESS) {
1400    *request = smpi_request_add_f(tmp);
1401  }
1402 }
1403
1404 void mpi_ibsend_ (void* buf, int* count, int* datatype, int *dest, int* tag, int* comm, int*  request, int* ierr) {
1405   MPI_Request tmp;
1406   *ierr = MPI_Ibsend(buf, *count, smpi_type_f2c(*datatype), *dest, *tag, smpi_comm_f2c(*comm), &tmp);
1407  if(*ierr == MPI_SUCCESS) {
1408    *request = smpi_request_add_f(tmp);
1409  }
1410 }
1411
1412 void mpi_comm_remote_group_ (int* comm, int*  group, int* ierr) {
1413   MPI_Group tmp;
1414  *ierr = MPI_Comm_remote_group(smpi_comm_f2c(*comm), &tmp);
1415  if(*ierr == MPI_SUCCESS) {
1416    *group = smpi_group_c2f(tmp);
1417  }
1418 }
1419
1420 void mpi_comm_remote_size_ (int* comm, int* size, int* ierr) {
1421  *ierr = MPI_Comm_remote_size(smpi_comm_f2c(*comm), size);
1422 }
1423
1424 void mpi_issend_ (void* buf, int* count, int* datatype, int *dest, int* tag, int* comm, int*  request, int* ierr) {
1425   MPI_Request tmp;
1426   *ierr = MPI_Issend(buf, *count, smpi_type_f2c(*datatype), *dest, *tag, smpi_comm_f2c(*comm), &tmp);
1427  if(*ierr == MPI_SUCCESS) {
1428    *request = smpi_request_add_f(tmp);
1429  }
1430 }
1431
1432 void mpi_probe_ (int* source, int* tag, int* comm, MPI_Status*  status, int* ierr) {
1433  *ierr = MPI_Probe(*source, *tag, smpi_comm_f2c(*comm), FORT_STATUS_IGNORE(status));
1434 }
1435
1436 void mpi_attr_delete_ (int* comm, int* keyval, int* ierr) {
1437  *ierr = MPI_Attr_delete(smpi_comm_f2c(*comm), *keyval);
1438 }
1439
1440 void mpi_attr_put_ (int* comm, int* keyval, void* attr_value, int* ierr) {
1441  *ierr = MPI_Attr_put(smpi_comm_f2c(*comm), *keyval, attr_value);
1442 }
1443
1444 void mpi_rsend_init_ (void* buf, int* count, int* datatype, int *dest, int* tag, int* comm, int*  request, int* ierr) {
1445   MPI_Request tmp;
1446   *ierr = MPI_Rsend_init(buf, *count, smpi_type_f2c(*datatype), *dest, *tag, smpi_comm_f2c(*comm), &tmp);
1447  if(*ierr == MPI_SUCCESS) {
1448    *request = smpi_request_add_f(tmp);
1449  }
1450 }
1451
1452 void mpi_keyval_create_ (void* copy_fn, void* delete_fn, int* keyval, void* extra_state, int* ierr) {
1453  *ierr = MPI_Keyval_create(reinterpret_cast<MPI_Copy_function*>(copy_fn),reinterpret_cast<MPI_Delete_function*>(delete_fn), keyval, extra_state);
1454 }
1455
1456 void mpi_keyval_free_ (int* keyval, int* ierr) {
1457  *ierr = MPI_Keyval_free(keyval);
1458 }
1459
1460 void mpi_test_cancelled_ (MPI_Status*  status, int* flag, int* ierr) {
1461  *ierr = MPI_Test_cancelled(status, flag);
1462 }
1463
1464 void mpi_pack_ (void* inbuf, int* incount, int* type, void* outbuf, int* outcount, int* position, int* comm, int* ierr) {
1465  *ierr = MPI_Pack(inbuf, *incount, smpi_type_f2c(*type), outbuf, *outcount, position, smpi_comm_f2c(*comm));
1466 }
1467
1468 void mpi_get_elements_ (MPI_Status*  status, int* datatype, int* elements, int* ierr) {
1469  *ierr = MPI_Get_elements(status, smpi_type_f2c(*datatype), elements);
1470 }
1471
1472 void mpi_dims_create_ (int* nnodes, int* ndims, int* dims, int* ierr) {
1473  *ierr = MPI_Dims_create(*nnodes, *ndims, dims);
1474 }
1475
1476 void mpi_iprobe_ (int* source, int* tag, int* comm, int* flag, MPI_Status*  status, int* ierr) {
1477  *ierr = MPI_Iprobe(*source, *tag, smpi_comm_f2c(*comm), flag, status);
1478 }
1479
1480 void mpi_type_get_envelope_ ( int* datatype, int *num_integers, int *num_addresses, int *num_datatypes, int *combiner,
1481                               int* ierr){
1482  *ierr = MPI_Type_get_envelope(  smpi_type_f2c(*datatype), num_integers,
1483  num_addresses, num_datatypes, combiner);
1484 }
1485
1486 void mpi_type_get_contents_ (int* datatype, int* max_integers, int* max_addresses, int* max_datatypes,
1487                              int* array_of_integers, MPI_Aint* array_of_addresses,
1488  int* array_of_datatypes, int* ierr){
1489  *ierr = MPI_Type_get_contents(smpi_type_f2c(*datatype), *max_integers, *max_addresses,*max_datatypes,
1490                                array_of_integers, array_of_addresses, reinterpret_cast<MPI_Datatype*>(array_of_datatypes));
1491 }
1492
1493 void mpi_type_create_darray_ (int* size, int* rank, int* ndims, int* array_of_gsizes, int* array_of_distribs,
1494                               int* array_of_dargs, int* array_of_psizes,
1495  int* order, int* oldtype, int*newtype, int* ierr) {
1496   MPI_Datatype tmp;
1497   *ierr = MPI_Type_create_darray(*size, *rank, *ndims,  array_of_gsizes,
1498   array_of_distribs,  array_of_dargs,  array_of_psizes,
1499   *order,  smpi_type_f2c(*oldtype), &tmp) ;
1500   if(*ierr == MPI_SUCCESS) {
1501     *newtype = smpi_type_add_f(tmp);
1502   }
1503 }
1504
1505 void mpi_type_create_resized_ (int* oldtype,MPI_Aint* lb, MPI_Aint* extent, int*newtype, int* ierr){
1506   MPI_Datatype tmp;
1507   *ierr = MPI_Type_create_resized(smpi_type_f2c(*oldtype),*lb, *extent, &tmp);
1508   if(*ierr == MPI_SUCCESS) {
1509     *newtype = smpi_type_add_f(tmp);
1510   }
1511 }
1512
1513 void mpi_type_create_subarray_ (int* ndims,int *array_of_sizes, int *array_of_subsizes, int *array_of_starts,
1514                                 int* order, int* oldtype, int*newtype, int* ierr){
1515   MPI_Datatype tmp;
1516   *ierr = MPI_Type_create_subarray(*ndims,array_of_sizes, array_of_subsizes, array_of_starts, *order,
1517                                    smpi_type_f2c(*oldtype), &tmp);
1518   if(*ierr == MPI_SUCCESS) {
1519     *newtype = smpi_type_add_f(tmp);
1520   }
1521 }
1522
1523 void mpi_type_match_size_ (int* typeclass,int* size,int* datatype, int* ierr){
1524   MPI_Datatype tmp;
1525   *ierr = MPI_Type_match_size(*typeclass,*size,&tmp);
1526   if(*ierr == MPI_SUCCESS) {
1527     *datatype = smpi_type_c2f(tmp);
1528   }
1529 }
1530
1531 void mpi_alltoallw_ ( void *sendbuf, int *sendcnts, int *sdispls, int* sendtypes, void *recvbuf, int *recvcnts,
1532                       int *rdispls, int* recvtypes, int* comm, int* ierr){
1533  *ierr = MPI_Alltoallw( sendbuf, sendcnts, sdispls, reinterpret_cast<MPI_Datatype*>(sendtypes), recvbuf, recvcnts, rdispls,
1534                         reinterpret_cast<MPI_Datatype*>(recvtypes), smpi_comm_f2c(*comm));
1535 }
1536
1537 void mpi_exscan_ (void *sendbuf, void *recvbuf, int* count, int* datatype, int* op, int* comm, int* ierr){
1538  *ierr = MPI_Exscan(sendbuf, recvbuf, *count, smpi_type_f2c(*datatype), smpi_op_f2c(*op), smpi_comm_f2c(*comm));
1539 }
1540
1541 void mpi_comm_set_name_ (int* comm, char* name, int* ierr, int size){
1542  char* tname = xbt_new(char, size+1);
1543  strncpy(tname, name, size);
1544  tname[size]='\0';
1545  *ierr = MPI_Comm_set_name (smpi_comm_f2c(*comm), tname);
1546  xbt_free(tname);
1547 }
1548
1549 void mpi_comm_dup_with_info_ (int* comm, int* info, int* newcomm, int* ierr){
1550   MPI_Comm tmp;
1551   *ierr = MPI_Comm_dup_with_info(smpi_comm_f2c(*comm),smpi_info_f2c(*info),&tmp);
1552   if(*ierr == MPI_SUCCESS) {
1553     *newcomm = smpi_comm_add_f(tmp);
1554   }
1555 }
1556
1557 void mpi_comm_split_type_ (int* comm, int* split_type, int* key, int* info, int* newcomm, int* ierr){
1558   MPI_Comm tmp;
1559   *ierr = MPI_Comm_split_type(smpi_comm_f2c(*comm), *split_type, *key, smpi_info_f2c(*info), &tmp);
1560   if(*ierr == MPI_SUCCESS) {
1561     *newcomm = smpi_comm_add_f(tmp);
1562   }
1563 }
1564
1565 void mpi_comm_set_info_ (int* comm, int* info, int* ierr){
1566  *ierr = MPI_Comm_set_info (smpi_comm_f2c(*comm), smpi_info_f2c(*info));
1567 }
1568
1569 void mpi_comm_get_info_ (int* comm, int* info, int* ierr){
1570  MPI_Info tmp;
1571  *ierr = MPI_Comm_get_info (smpi_comm_f2c(*comm), &tmp);
1572  if(*ierr==MPI_SUCCESS){
1573    *info = smpi_info_c2f(tmp);
1574  }
1575 }
1576
1577 void mpi_comm_create_errhandler_ ( void *function, void *errhandler, int* ierr){
1578  *ierr = MPI_Comm_create_errhandler( reinterpret_cast<MPI_Comm_errhandler_fn*>(function), static_cast<MPI_Errhandler*>(errhandler));
1579 }
1580
1581 void mpi_add_error_class_ ( int *errorclass, int* ierr){
1582  *ierr = MPI_Add_error_class( errorclass);
1583 }
1584
1585 void mpi_add_error_code_ (  int* errorclass, int *errorcode, int* ierr){
1586  *ierr = MPI_Add_error_code(*errorclass, errorcode);
1587 }
1588
1589 void mpi_add_error_string_ ( int* errorcode, char *string, int* ierr){
1590  *ierr = MPI_Add_error_string(*errorcode, string);
1591 }
1592
1593 void mpi_comm_call_errhandler_ (int* comm,int* errorcode, int* ierr){
1594  *ierr = MPI_Comm_call_errhandler(smpi_comm_f2c(*comm), *errorcode);
1595 }
1596
1597 void mpi_info_dup_ (int* info, int* newinfo, int* ierr){
1598  MPI_Info tmp;
1599  *ierr = MPI_Info_dup(smpi_info_f2c(*info), &tmp);
1600  if(*ierr==MPI_SUCCESS){
1601    *newinfo= smpi_info_add_f(tmp);
1602  }
1603 }
1604
1605 void mpi_info_get_valuelen_ ( int* info, char *key, int *valuelen, int *flag, int* ierr, unsigned int keylen){
1606   while(key[keylen-1]==' ')
1607     keylen--;
1608   while(*key==' '){//handle leading blanks
1609     keylen--;
1610     key++;
1611   }
1612   char* tkey = xbt_new(char, keylen+1);
1613   strncpy(tkey, key, keylen);
1614   tkey[keylen]='\0';
1615   *ierr = MPI_Info_get_valuelen( smpi_info_f2c(*info), tkey, valuelen, flag);
1616   xbt_free(tkey);
1617 }
1618
1619 void mpi_info_delete_ (int* info, char *key, int* ierr, unsigned int keylen){
1620   while(key[keylen-1]==' ')
1621     keylen--;
1622   while(*key==' '){//handle leading blanks
1623     keylen--;
1624     key++;
1625   }
1626   char* tkey = xbt_new(char, keylen+1);
1627   strncpy(tkey, key, keylen);
1628   tkey[keylen]='\0';
1629   *ierr = MPI_Info_delete(smpi_info_f2c(*info), tkey);
1630   xbt_free(tkey);
1631 }
1632
1633 void mpi_info_get_nkeys_ ( int* info, int *nkeys, int* ierr){
1634  *ierr = MPI_Info_get_nkeys(  smpi_info_f2c(*info), nkeys);
1635 }
1636
1637 void mpi_info_get_nthkey_ ( int* info, int* n, char *key, int* ierr, unsigned int keylen){
1638   *ierr = MPI_Info_get_nthkey( smpi_info_f2c(*info), *n, key);
1639   unsigned int i = 0;
1640   for (i=strlen(key); i<keylen; i++)
1641     key[i]=' ';
1642 }
1643
1644 void mpi_get_version_ (int *version,int *subversion, int* ierr){
1645  *ierr = MPI_Get_version (version,subversion);
1646 }
1647
1648 void mpi_get_library_version_ (char *version,int *len, int* ierr){
1649  *ierr = MPI_Get_library_version (version,len);
1650 }
1651
1652 void mpi_request_get_status_ ( int* request, int *flag, MPI_Status* status, int* ierr){
1653  *ierr = MPI_Request_get_status( smpi_request_f2c(*request), flag, status);
1654 }
1655
1656 void mpi_grequest_start_ ( void *query_fn, void *free_fn, void *cancel_fn, void *extra_state, int*request, int* ierr){
1657   MPI_Request tmp;
1658   *ierr = MPI_Grequest_start( reinterpret_cast<MPI_Grequest_query_function*>(query_fn), reinterpret_cast<MPI_Grequest_free_function*>(free_fn),
1659                               reinterpret_cast<MPI_Grequest_cancel_function*>(cancel_fn), extra_state, &tmp);
1660  if(*ierr == MPI_SUCCESS) {
1661    *request = smpi_request_add_f(tmp);
1662  }
1663 }
1664
1665 void mpi_grequest_complete_ ( int* request, int* ierr){
1666  *ierr = MPI_Grequest_complete( smpi_request_f2c(*request));
1667 }
1668
1669 void mpi_status_set_cancelled_ (MPI_Status* status,int* flag, int* ierr){
1670  *ierr = MPI_Status_set_cancelled(status,*flag);
1671 }
1672
1673 void mpi_status_set_elements_ ( MPI_Status* status, int* datatype, int* count, int* ierr){
1674  *ierr = MPI_Status_set_elements( status, smpi_type_f2c(*datatype), *count);
1675 }
1676
1677 void mpi_comm_connect_ ( char *port_name, int* info, int* root, int* comm, int*newcomm, int* ierr){
1678   MPI_Comm tmp;
1679   *ierr = MPI_Comm_connect( port_name, *reinterpret_cast<MPI_Info*>(info), *root, smpi_comm_f2c(*comm), &tmp);
1680   if(*ierr == MPI_SUCCESS) {
1681     *newcomm = smpi_comm_add_f(tmp);
1682   }
1683 }
1684
1685 void mpi_publish_name_ ( char *service_name, int* info, char *port_name, int* ierr){
1686  *ierr = MPI_Publish_name( service_name, *reinterpret_cast<MPI_Info*>(info), port_name);
1687 }
1688
1689 void mpi_unpublish_name_ ( char *service_name, int* info, char *port_name, int* ierr){
1690  *ierr = MPI_Unpublish_name( service_name, *reinterpret_cast<MPI_Info*>(info), port_name);
1691 }
1692
1693 void mpi_lookup_name_ ( char *service_name, int* info, char *port_name, int* ierr){
1694  *ierr = MPI_Lookup_name( service_name, *reinterpret_cast<MPI_Info*>(info), port_name);
1695 }
1696
1697 void mpi_comm_join_ ( int* fd, int* intercomm, int* ierr){
1698   MPI_Comm tmp;
1699   *ierr = MPI_Comm_join( *fd, &tmp);
1700   if(*ierr == MPI_SUCCESS) {
1701     *intercomm = smpi_comm_add_f(tmp);
1702   }
1703 }
1704
1705 void mpi_open_port_ ( int* info, char *port_name, int* ierr){
1706  *ierr = MPI_Open_port( *reinterpret_cast<MPI_Info*>(info),port_name);
1707 }
1708
1709 void mpi_close_port_ ( char *port_name, int* ierr){
1710  *ierr = MPI_Close_port( port_name);
1711 }
1712
1713 void mpi_comm_accept_ ( char *port_name, int* info, int* root, int* comm, int*newcomm, int* ierr){
1714   MPI_Comm tmp;
1715   *ierr = MPI_Comm_accept( port_name, *reinterpret_cast<MPI_Info*>(info), *root, smpi_comm_f2c(*comm), &tmp);
1716   if(*ierr == MPI_SUCCESS) {
1717     *newcomm = smpi_comm_add_f(tmp);
1718   }
1719 }
1720
1721 void mpi_comm_spawn_ ( char *command, char *argv, int* maxprocs, int* info, int* root, int* comm, int* intercomm,
1722                        int* array_of_errcodes, int* ierr){
1723   MPI_Comm tmp;
1724   *ierr = MPI_Comm_spawn( command, nullptr, *maxprocs, *reinterpret_cast<MPI_Info*>(info), *root, smpi_comm_f2c(*comm), &tmp,
1725                           array_of_errcodes);
1726   if(*ierr == MPI_SUCCESS) {
1727     *intercomm = smpi_comm_add_f(tmp);
1728   }
1729 }
1730
1731 void mpi_comm_spawn_multiple_ ( int* count, char *array_of_commands, char** array_of_argv, int* array_of_maxprocs,
1732                                 int* array_of_info, int* root,
1733  int* comm, int* intercomm, int* array_of_errcodes, int* ierr){
1734  MPI_Comm tmp;
1735  *ierr = MPI_Comm_spawn_multiple(* count, &array_of_commands, &array_of_argv, array_of_maxprocs,
1736  reinterpret_cast<MPI_Info*>(array_of_info), *root, smpi_comm_f2c(*comm), &tmp, array_of_errcodes);
1737  if(*ierr == MPI_SUCCESS) {
1738    *intercomm = smpi_comm_add_f(tmp);
1739  }
1740 }
1741
1742 void mpi_comm_get_parent_ ( int* parent, int* ierr){
1743   MPI_Comm tmp;
1744   *ierr = MPI_Comm_get_parent( &tmp);
1745   if(*ierr == MPI_SUCCESS) {
1746     *parent = smpi_comm_c2f(tmp);
1747   }
1748 }
1749
1750 void mpi_file_close_ ( int* file, int* ierr){
1751   *ierr= MPI_File_close(reinterpret_cast<MPI_File*>(*file));
1752 }
1753
1754 void mpi_file_delete_ ( char* filename, int* info, int* ierr){
1755   *ierr= MPI_File_delete(filename, smpi_info_f2c(*info));
1756 }
1757
1758 void mpi_file_open_ ( int* comm, char* filename, int* amode, int* info, int* fh, int* ierr){
1759   *ierr= MPI_File_open(smpi_comm_f2c(*comm), filename, *amode, smpi_info_f2c(*info), reinterpret_cast<MPI_File*>(*fh));
1760 }
1761
1762 void mpi_file_set_view_ ( int* fh, long long int* offset, int* etype, int* filetype, char* datarep, int* info, int* ierr){
1763   *ierr= MPI_File_set_view(reinterpret_cast<MPI_File>(*fh) , reinterpret_cast<MPI_Offset>(*offset), smpi_type_f2c(*etype), smpi_type_f2c(*filetype), datarep, smpi_info_f2c(*info));
1764 }
1765
1766 void mpi_file_read_ ( int* fh, void* buf, int* count, int* datatype, MPI_Status* status, int* ierr){
1767   *ierr=  MPI_File_read(reinterpret_cast<MPI_File>(*fh), buf, *count, smpi_type_f2c(*datatype), status);
1768 }
1769
1770 void mpi_file_write_ ( int* fh, void* buf, int* count, int* datatype, MPI_Status* status, int* ierr){
1771   *ierr=  MPI_File_write(reinterpret_cast<MPI_File>(*fh), buf, *count, smpi_type_f2c(*datatype), status);
1772 }
1773
1774 } // extern "C"