Logo AND Algorithmique Numérique Distribuée

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