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