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