Logo AND Algorithmique Numérique Distribuée

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