src/smpi/smpi_f77.c

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