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