Logo AND Algorithmique Numérique Distribuée

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