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