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