src/smpi/smpi_f77.c

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