src/smpi/smpi_f77.cpp

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