Logo AND Algorithmique Numérique Distribuée

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