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