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