src/smpi/smpi_f77.c

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