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