Logo AND Algorithmique Numérique Distribuée

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