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