Actual source code: mpi.h
1: /*
2: This is a special set of bindings for uni-processor use of MPI by the PETSc library.
3:
4: NOT ALL THE MPI CALLS ARE IMPLEMENTED CORRECTLY! Only those needed in PETSc.
6: For example,
7: * Does not implement send to self.
8: * Does not implement attributes correctly.
9: */
11: /*
12: The following info is a response to one of the petsc-maint questions
13: regarding MPIUNI.
15: MPIUNI was developed with the aim of getting PETSc compiled, and
16: usable in the absence of a full MPI implementation. With this, we
17: were able to provide PETSc on Windows, Windows64 even before any MPI
18: implementation was available on these platforms. [Or with certain
19: compilers - like borland, that do not have a useable MPI
20: implementation]
22: However - providing a seqential, standards compliant MPI
23: implementation is *not* the goal of MPIUNI. The development strategy
24: was - to make enough changes to it so that PETSc sources, examples
25: compile without errors, and runs in the uni-processor mode. This is
26: the reason each function is not documented.
28: PETSc usage of MPIUNI is primarily from C. However a minimal fortran
29: interface is also provided - to get PETSc fortran examples with a
30: few MPI calls working.
32: One of the optimzation with MPIUNI, is to avoid the function call
33: overhead, when possible. Hence most of the C functions are
34: implemented as macros. However the function calls cannot be avoided
35: with fortran usage.
37: Most PETSc objects have both sequential and parallel
38: implementations, which are separate. For eg: We have two types of
39: sparse matrix storage formats - SeqAIJ, and MPIAIJ. Some MPI
40: routines are used in the Seq part, but most of them are used in the
41: MPI part. The send/receive calls can be found mostly in the MPI
42: part.
44: When MPIUNI is used, only the Seq version of the PETSc objects are
45: used, even though the MPI variant of the objects are compiled. Since
46: there are no send/receive calls in the Seq variant, PETSc works fine
47: with MPIUNI in seq mode.
49: The reason some send/receive functions are defined to abort(), is to
50: detect sections of code that use send/receive functions, and gets
51: executed in the sequential mode. (which shouldn't happen in case of
52: PETSc).
54: Proper implementation of send/receive would involve writing a
55: function for each of them. Inside each of these functions, we have
56: to check if the send is to self or receive is from self, and then
57: doing the buffering accordingly (until the receive is called) - or
58: what if a nonblocking receive is called, do a copy etc.. Handling
59: the buffering aspects might be complicated enough, that in this
60: case, a proper implementation of MPI might as well be used. This is
61: the reason the send to self is not implemented in MPIUNI, and never
62: will be.
63:
64: Proper implementations of MPI [for eg: MPICH & OpenMPI] are
65: available for most machines. When these packages are available, Its
66: generally preferable to use one of them instead of MPIUNI - even if
67: the user is using PETSc sequentially.
69: - MPIUNI does not support all MPI functions [or functionality].
70: Hence it might not work with external packages or user code that
71: might have MPI calls in it.
73: - MPIUNI is not a standards compliant implementation for np=1.
74: For eg: if the user code has send/recv to self, then it will
75: abort. [Similar issues with a number of other MPI functionality]
76: However MPICH & OpenMPI are the correct implementations of MPI
77: standard for np=1.
79: - When user code uses multiple MPI based packages that have their
80: own *internal* stubs equivalent to MPIUNI - in sequential mode,
81: invariably these multiple implementations of MPI for np=1 conflict
82: with each other. The correct thing to do is: make all such
83: packages use the *same* MPI implementation for np=1. MPICH/OpenMPI
84: satisfy this requirement correctly [and hence the correct choice].
86: - Using MPICH/OpenMPI sequentially should have minimal
87: disadvantages. [for eg: these binaries can be run without
88: mpirun/mpiexec as ./executable, without requiring any extra
89: configurations for ssh/rsh/daemons etc..]. This should not be a
90: reason to avoid these packages for sequential use.
92: Instructions for building standalone MPIUNI [for eg: linux/gcc+gfortran]:
93: - extract include/mpiuni/mpi.h,mpif.f, src/sys/mpiuni/mpi.c from PETSc
94: - remove reference to petscconf.h from mpi.h
95: - gcc -c mpi.c -DPETSC_HAVE_STDLIB_H -DPETSC_HAVE_FORTRAN_UNDERSCORE
96: - ar cr libmpiuni.a mpi.o
98: */
103: /* Requred by abort() in mpi.c & for win64 */
104: #include "petscconf.h"
106: #if defined(__cplusplus)
108: #endif
110: /* require an int variable large enough to hold a pointer */
111: #if !defined(MPIUNI_INTPTR)
112: #define MPIUNI_INTPTR long
113: #endif
115: /*
117: MPIUNI_TMP is used in the macros below only to stop various C/C++ compilers
118: from generating warning messages about unused variables while compiling PETSc.
119: */
122: #define MPI_COMM_WORLD 1
123: #define MPI_COMM_SELF MPI_COMM_WORLD
124: #define MPI_COMM_NULL 0
125: #define MPI_SUCCESS 0
126: #define MPI_IDENT 0
127: #define MPI_CONGRUENT 0
128: #define MPI_SIMILAR 0
129: #define MPI_UNEQUAL 3
130: #define MPI_ANY_SOURCE (-2)
131: #define MPI_KEYVAL_INVALID 0
132: #define MPI_ERR_UNKNOWN 18
133: #define MPI_ERR_INTERN 21
134: #define MPI_ERR_OTHER 1
135: #define MPI_TAG_UB 0
136: #define MPI_ERRORS_RETURN 0
137: #define MPI_UNDEFINED (-32766)
139: /* External types */
140: typedef int MPI_Comm;
141: typedef void *MPI_Request;
142: typedef void *MPI_Group;
143: typedef struct {int MPI_TAG,MPI_SOURCE,MPI_ERROR;} MPI_Status;
144: typedef char *MPI_Errhandler;
145: typedef int MPI_Fint;
146: typedef int MPI_File;
147: typedef int MPI_Info;
148: typedef int MPI_Offset;
152: /* In order to handle datatypes, we make them into "sizeof(raw-type)";
153: this allows us to do the MPIUNI_Memcpy's easily */
154: #define MPI_Datatype int
155: #define MPI_FLOAT sizeof(float)
156: #define MPI_DOUBLE sizeof(double)
157: #define MPI_LONG_DOUBLE sizeof(long double)
158: #define MPI_CHAR sizeof(char)
159: #define MPI_BYTE sizeof(char)
160: #define MPI_INT sizeof(int)
161: #define MPI_LONG sizeof(long)
162: #define MPI_LONG_LONG_INT sizeof(long long)
163: #define MPI_SHORT sizeof(short)
164: #define MPI_UNSIGNED_SHORT sizeof(unsigned short)
165: #define MPI_UNSIGNED sizeof(unsigned)
166: #define MPI_UNSIGNED_CHAR sizeof(unsigned char)
167: #define MPI_UNSIGNED_LONG sizeof(unsigned long)
168: #define MPI_COMPLEX 2*sizeof(float)
169: #define MPI_C_DOUBLE_COMPLEX 2*sizeof(double)
170: #define MPI_FLOAT_INT (sizeof(float) + sizeof(int))
171: #define MPI_DOUBLE_INT (sizeof(double) + sizeof(int))
172: #define MPI_LONG_INT (sizeof(long) + sizeof(int))
173: #define MPI_SHORT_INT (sizeof(short) + sizeof(int))
174: #define MPI_2INT (2* sizeof(int))
176: #define MPI_REQUEST_NULL ((MPI_Request)0)
177: #define MPI_GROUP_NULL ((MPI_Group)0)
178: #define MPI_INFO_NULL ((MPI_Info)0)
179: #define MPI_BOTTOM (void *)0
180: typedef int MPI_Op;
182: #define MPI_MODE_RDONLY 0
183: #define MPI_MODE_WRONLY 0
184: #define MPI_MODE_CREATE 0
186: #define MPI_SUM 0
187: #define MPI_MAX 0
188: #define MPI_MIN 0
189: #define MPI_ANY_TAG (-1)
190: #define MPI_DATATYPE_NULL 0
191: #define MPI_PACKED 0
192: #define MPI_MAX_ERROR_STRING 2056
193: #define MPI_STATUS_IGNORE (MPI_Status *)1
194: #define MPI_ORDER_FORTRAN 57
195: #define MPI_IN_PLACE (void *) -1
197: /*
198: Prototypes of some functions which are implemented in mpi.c
199: */
200: typedef int (MPI_Copy_function)(MPI_Comm,int,void *,void *,void *,int *);
201: typedef int (MPI_Delete_function)(MPI_Comm,int,void *,void *);
202: typedef void (MPI_User_function)(void*, void *, int *, MPI_Datatype *);
204: /*
205: In order that the PETSc MPIUNI can be used with another package that has its
206: own MPIUni we map the following function names to a unique PETSc name. Those functions
207: are defined in mpi.c and put into the libpetscsys.a or libpetsc.a library.
209: Note that this does not work for the MPIUni Fortran symbols which are explicitly in the
210: PETSc libraries unless the flag MPIUNI_AVOID_MPI_NAMESPACE is set.
211: */
212: #define MPI_Abort Petsc_MPI_Abort
213: #define MPI_Attr_get Petsc_MPI_Attr_get
214: #define MPI_Keyval_free Petsc_MPI_Keyval_free
215: #define MPI_Attr_put Petsc_MPI_Attr_put
216: #define MPI_Attr_delete Petsc_MPI_Attr_delete
217: #define MPI_Keyval_create Petsc_MPI_Keyval_create
218: #define MPI_Comm_free Petsc_MPI_Comm_free
219: #define MPI_Comm_dup Petsc_MPI_Comm_dup
220: #define MPI_Comm_create Petsc_MPI_Comm_create
221: #define MPI_Init Petsc_MPI_Init
222: #define MPI_Finalize Petsc_MPI_Finalize
223: #define MPI_Initialized Petsc_MPI_Initialized
224: #define MPI_Finalized Petsc_MPI_Finalized
241: #define MPI_Aint MPIUNI_INTPTR
242: /*
243: Routines we have replace with macros that do nothing
244: Some return error codes others return success
245: */
247: #define MPI_Comm_f2c(comm) (MPI_Comm)(comm)
248: #define MPI_Comm_c2f(comm) (MPI_Fint)(comm)
250: #define MPI_Send(buf,count,datatype,dest,tag,comm) \
251: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (buf),\
252: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (count),\
253: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (datatype),\
254: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (dest),\
255: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (tag),\
256: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
257: MPI_Abort(MPI_COMM_WORLD,0))
258: #define MPI_Recv(buf,count,datatype,source,tag,comm,status) \
259: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (buf),\
260: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (count),\
261: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (datatype),\
262: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (source),\
263: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (tag),\
264: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
265: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (status),\
266: MPI_Abort(MPI_COMM_WORLD,0))
267: #define MPI_Get_count(status, datatype,count) \
268: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (status),\
269: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (datatype),\
270: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (count),\
271: MPI_Abort(MPI_COMM_WORLD,0))
272: #define MPI_Bsend(buf,count,datatype,dest,tag,comm) \
273: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (buf),\
274: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (count),\
275: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (datatype),\
276: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (dest),\
277: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (tag),\
278: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
279: MPI_Abort(MPI_COMM_WORLD,0))
280: #define MPI_Ssend(buf,count, datatype,dest,tag,comm) \
281: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (buf),\
282: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (count),\
283: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (datatype),\
284: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (dest),\
285: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (tag),\
286: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
287: MPI_Abort(MPI_COMM_WORLD,0))
288: #define MPI_Rsend(buf,count, datatype,dest,tag,comm) \
289: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (buf),\
290: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (count),\
291: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (datatype),\
292: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (dest),\
293: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (tag),\
294: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
295: MPI_Abort(MPI_COMM_WORLD,0))
296: #define MPI_Buffer_attach(buffer,size) \
297: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (buffer),\
298: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (size),\
299: MPI_SUCCESS)
300: #define MPI_Buffer_detach(buffer,size)\
301: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (buffer),\
302: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (size),\
303: MPI_SUCCESS)
304: #define MPI_Ibsend(buf,count, datatype,dest,tag,comm,request) \
305: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (buf),\
306: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (count),\
307: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (datatype),\
308: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (dest),\
309: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (tag),\
310: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
311: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (request),\
312: MPI_Abort(MPI_COMM_WORLD,0))
313: #define MPI_Issend(buf,count, datatype,dest,tag,comm,request) \
314: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (buf),\
315: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (count),\
316: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (datatype),\
317: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (dest),\
318: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (tag),\
319: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
320: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (request),\
321: MPI_Abort(MPI_COMM_WORLD,0))
322: #define MPI_Irsend(buf,count, datatype,dest,tag,comm,request) \
323: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (buf),\
324: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (count),\
325: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (datatype),\
326: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (dest),\
327: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (tag),\
328: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
329: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (request),\
330: MPI_Abort(MPI_COMM_WORLD,0))
331: #define MPI_Irecv(buf,count, datatype,source,tag,comm,request) \
332: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (buf),\
333: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (count),\
334: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (datatype),\
335: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (source),\
336: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (tag),\
337: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
338: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (request),\
339: MPI_Abort(MPI_COMM_WORLD,0))
340: #define MPI_Isend(buf,count, datatype,dest,tag,comm,request) \
341: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (buf),\
342: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (count),\
343: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (datatype),\
344: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (dest),\
345: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (tag),\
346: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
347: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (request),\
348: MPI_Abort(MPI_COMM_WORLD,0))
349: #define MPI_Wait(request,status) \
350: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (request),\
351: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (status),\
352: MPI_SUCCESS)
353: #define MPI_Test(request,flag,status) \
354: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (request),\
355: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (status),\
356: *(flag) = 0, \
357: MPI_SUCCESS)
358: #define MPI_Request_free(request) \
359: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (request),\
360: MPI_SUCCESS)
361: #define MPI_Waitany(a,b,c,d) \
362: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (a),\
363: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (b),\
364: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (c),\
365: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (d),\
366: MPI_SUCCESS)
367: #define MPI_Testany(a,b,c,d,e) \
368: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (a),\
369: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (b),\
370: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (c),\
371: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (d),\
372: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (e),\
373: MPI_SUCCESS)
374: #define MPI_Waitall(count,array_of_requests,array_of_statuses) \
375: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (count),\
376: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (array_of_requests),\
377: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (array_of_statuses),\
378: MPI_SUCCESS)
379: #define MPI_Testall(count,array_of_requests,flag,array_of_statuses) \
380: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (count),\
381: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (array_of_requests),\
382: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (flag),\
383: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (array_of_statuses),\
384: MPI_SUCCESS)
385: #define MPI_Waitsome(incount,array_of_requests,outcount,\
386: array_of_indices,array_of_statuses) \
387: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (incount),\
388: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (array_of_requests),\
389: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (outcount),\
390: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (array_of_indices),\
391: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (array_of_statuses),\
392: MPI_SUCCESS)
393: #define MPI_Comm_group(comm,group) \
394: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
395: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (group),\
396: MPI_SUCCESS)
397: #define MPI_Group_incl(group,n,ranks,newgroup) \
398: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (group),\
399: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (n),\
400: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (ranks),\
401: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (newgroup),\
402: MPI_SUCCESS)
403: #define MPI_Testsome(incount,array_of_requests,outcount,\
404: array_of_indices,array_of_statuses) MPI_SUCCESS
405: #define MPI_Iprobe(source,tag,comm,flag,status) (*(flag)=0, MPI_SUCCESS)
406: #define MPI_Probe(source,tag,comm,status) MPI_SUCCESS
407: #define MPI_Cancel(request) (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (request),MPI_SUCCESS)
408: #define MPI_Test_cancelled(status,flag) (*(flag)=0,MPI_SUCCESS)
409: #define MPI_Send_init(buf,count, datatype,dest,tag,comm,request) \
410: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (buf),\
411: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (count),\
412: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (datatype),\
413: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (dest),\
414: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (tag),\
415: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
416: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (request),\
417: MPI_SUCCESS)
418: #define MPI_Bsend_init(buf,count, datatype,dest,tag,comm,request) \
419: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (buf),\
420: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (count),\
421: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (datatype),\
422: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (dest),\
423: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (tag),\
424: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
425: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (request),\
426: MPI_SUCCESS)
427: #define MPI_Ssend_init(buf,count, datatype,dest,tag,comm,request) \
428: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (buf),\
429: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (count),\
430: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (datatype),\
431: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (dest),\
432: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (tag),\
433: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
434: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (request),\
435: MPI_SUCCESS)
436: #define MPI_Bsend_init(buf,count, datatype,dest,tag,comm,request) \
437: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (buf),\
438: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (count),\
439: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (datatype),\
440: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (dest),\
441: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (tag),\
442: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
443: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (request),\
444: MPI_SUCCESS)
445: #define MPI_Rsend_init(buf,count, datatype,dest,tag,comm,request) \
446: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (buf),\
447: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (count),\
448: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (datatype),\
449: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (dest),\
450: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (tag),\
451: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
452: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (request),\
453: MPI_SUCCESS)
454: #define MPI_Recv_init(buf,count, datatype,source,tag,comm,request) \
455: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (buf),\
456: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (count),\
457: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (datatype),\
458: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (source),\
459: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (tag),\
460: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
461: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (request),\
462: MPI_SUCCESS)
463: #define MPI_Start(request) (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (request),MPI_SUCCESS)
464: #define MPI_Startall(count,array_of_requests) \
465: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (count),\
466: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (array_of_requests),\
467: MPI_SUCCESS)
468: #define MPI_Op_create(function,commute,op) \
469: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (function),\
470: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (commute),\
471: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (op),\
472: MPI_SUCCESS)
473: #define MPI_Op_free(op) \
474: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (op),\
475: MPI_SUCCESS)
476: /* Need to determine sizeof "sendtype" */
477: #define MPI_Sendrecv(sendbuf,sendcount, sendtype,\
478: dest,sendtag,recvbuf,recvcount,\
479: recvtype,source,recvtag,\
480: comm,status) \
481: MPIUNI_Memcpy(recvbuf,sendbuf,(sendcount) * (sendtype))
482: #define MPI_Sendrecv_replace(buf,count, datatype,dest,sendtag,\
483: source,recvtag,comm,status) MPI_SUCCESS
484: #define MPI_Type_contiguous(count, oldtype,newtype) \
485: (*(newtype) = (count)*(oldtype),MPI_SUCCESS)
486: #define MPI_Type_vector(count,blocklength,stride,oldtype, newtype) MPI_SUCCESS
487: #define MPI_Type_hvector(count,blocklength,stride,oldtype, newtype) MPI_SUCCESS
488: #define MPI_Type_indexed(count,array_of_blocklengths,\
489: array_of_displacements, oldtype,\
490: newtype) MPI_SUCCESS
491: #define MPI_Type_hindexed(count,array_of_blocklengths,\
492: array_of_displacements, oldtype,\
493: newtype) MPI_SUCCESS
494: #define MPI_Type_struct(count,array_of_blocklengths,\
495: array_of_displacements,\
496: array_of_types, newtype) MPI_SUCCESS
497: #define MPI_Address(location,address) \
498: (*(address) = (MPIUNI_INTPTR)(char *)(location),MPI_SUCCESS)
499: #define MPI_Type_extent(datatype,extent) *(extent) = datatype
500: #define MPI_Type_size(datatype,size) *(size) = datatype
501: #define MPI_Type_lb(datatype,displacement) \
502: MPI_Abort(MPI_COMM_WORLD,0)
503: #define MPI_Type_ub(datatype,displacement) \
504: MPI_Abort(MPI_COMM_WORLD,0)
505: #define MPI_Type_commit(datatype) (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (datatype),\
506: MPI_SUCCESS)
507: #define MPI_Type_free(datatype) MPI_SUCCESS
508: #define MPI_Get_elements(status, datatype,count) \
509: MPI_Abort(MPI_COMM_WORLD,0)
510: #define MPI_Pack(inbuf,incount, datatype,outbuf,\
511: outsize,position, comm) \
512: MPI_Abort(MPI_COMM_WORLD,0)
513: #define MPI_Unpack(inbuf,insize,position,outbuf,\
514: outcount, datatype,comm) \
515: MPI_Abort(MPI_COMM_WORLD,0)
516: #define MPI_Pack_size(incount, datatype,comm,size) \
517: MPI_Abort(MPI_COMM_WORLD,0)
518: #define MPI_Barrier(comm) \
519: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
520: MPI_SUCCESS)
521: #define MPI_Bcast(buffer,count,datatype,root,comm) \
522: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (buffer),\
523: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (count),\
524: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (datatype),\
525: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
526: MPI_SUCCESS)
527: #define MPI_Gather(sendbuf,sendcount, sendtype,\
528: recvbuf,recvcount, recvtype,\
529: root,comm) \
530: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (recvcount),\
531: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (root),\
532: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (recvtype),\
533: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
534: MPIUNI_Memcpy(recvbuf,sendbuf,(sendcount)* (sendtype)),\
535: MPI_SUCCESS)
536: #define MPI_Gatherv(sendbuf,sendcount, sendtype,\
537: recvbuf,recvcounts,displs,\
538: recvtype,root,comm) \
539: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (recvcounts),\
540: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (displs),\
541: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (recvtype),\
542: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (root),\
543: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
544: MPIUNI_Memcpy(recvbuf,sendbuf,(sendcount)* (sendtype)),\
545: MPI_SUCCESS)
546: #define MPI_Scatter(sendbuf,sendcount, sendtype,\
547: recvbuf,recvcount, recvtype,\
548: root,comm) \
549: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (sendbuf),\
550: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (sendcount),\
551: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (sendtype),\
552: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (recvbuf),\
553: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (recvcount),\
554: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (recvtype),\
555: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (root),\
556: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),MPI_Abort(MPI_COMM_WORLD,0))
557: #define MPI_Scatterv(sendbuf,sendcounts,displs,\
558: sendtype, recvbuf,recvcount,\
559: recvtype,root,comm) \
560: (MPIUNI_Memcpy(recvbuf,sendbuf,(recvcount)*(recvtype)),\
561: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (displs),\
562: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (sendtype),\
563: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (sendcounts),\
564: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (root),\
565: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
566: MPI_SUCCESS)
567: #define MPI_Allgather(sendbuf,sendcount, sendtype,\
568: recvbuf,recvcount, recvtype,comm) \
569: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (recvcount),\
570: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (recvtype),\
571: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
572: MPIUNI_Memcpy(recvbuf,sendbuf,(sendcount)* (sendtype)),\
573: MPI_SUCCESS)
574: #define MPI_Allgatherv(sendbuf,sendcount, sendtype,\
575: recvbuf,recvcounts,displs,recvtype,comm) \
576: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (recvcounts),\
577: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (displs),\
578: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (recvtype),\
579: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
580: (sendbuf != MPI_IN_PLACE) ? MPIUNI_Memcpy((recvbuf),(sendbuf),(sendcount)* (sendtype)) : 0, \
581: MPI_SUCCESS)
582: #define MPI_Alltoall(sendbuf,sendcount, sendtype,\
583: recvbuf,recvcount, recvtype,comm) \
584: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (recvcount),\
585: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (recvtype),\
586: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
587: MPIUNI_Memcpy(recvbuf,sendbuf,(sendcount)* (sendtype)),\
588: MPI_SUCCESS)
589: #define MPI_Alltoallv(sendbuf,sendcounts,sdispls,\
590: sendtype, recvbuf,recvcounts,\
591: rdispls, recvtype,comm) MPI_Abort(MPI_COMM_WORLD,0)
592: #define MPI_Alltoallw(sendbuf,sendcounts,sdispls,\
593: sendtypes, recvbuf,recvcounts,\
594: rdispls, recvtypes,comm) MPI_Abort(MPI_COMM_WORLD,0)
595: #define MPI_Reduce(sendbuf, recvbuf,count,\
596: datatype,op,root,comm) \
597: (MPIUNI_Memcpy(recvbuf,sendbuf,(count)*(datatype)),\
598: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),MPI_SUCCESS)
599: #define MPI_Allreduce(sendbuf, recvbuf,count,datatype,op,comm) \
600: (MPIUNI_Memcpy(recvbuf,sendbuf,(count)*(datatype)),\
601: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),MPI_SUCCESS)
602: #define MPI_Scan(sendbuf, recvbuf,count,datatype,op,comm) \
603: (MPIUNI_Memcpy(recvbuf,sendbuf,(count)*(datatype)),\
604: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),MPI_SUCCESS)
605: #define MPI_Reduce_scatter(sendbuf, recvbuf,recvcounts,\
606: datatype,op,comm) \
607: MPI_Abort(MPI_COMM_WORLD,0)
608: #define MPI_Group_size(group,size) (*(size)=1,MPI_SUCCESS)
609: #define MPI_Group_rank(group,rank) (*(rank)=0,MPI_SUCCESS)
610: #define MPI_Group_translate_ranks (group1,n,ranks1,\
611: group2,ranks2) MPI_Abort(MPI_COMM_WORLD,0)
612: #define MPI_Group_compare(group1,group2,result) \
613: (*(result)=1,MPI_SUCCESS)
614: #define MPI_Group_union(group1,group2,newgroup) MPI_SUCCESS
615: #define MPI_Group_intersection(group1,group2,newgroup) MPI_SUCCESS
616: #define MPI_Group_difference(group1,group2,newgroup) MPI_SUCCESS
617: #define MPI_Group_excl(group,n,ranks,newgroup) MPI_SUCCESS
618: #define MPI_Group_range_incl(group,n,ranges,newgroup) MPI_SUCCESS
619: #define MPI_Group_range_excl(group,n,ranges,newgroup) MPI_SUCCESS
620: #define MPI_Group_free(group) \
621: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (group),\
622: MPI_SUCCESS)
623: #define MPI_Comm_size(comm,size) \
624: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
625: *(size)=1,\
626: MPI_SUCCESS)
627: #define MPI_Comm_rank(comm,rank) \
628: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
629: *(rank)=0,\
630: MPI_SUCCESS)
631: #define MPI_Comm_compare(comm1,comm2,result) \
632: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm1),\
633: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm2),\
634: *(result)=MPI_IDENT,\
635: MPI_SUCCESS)
636: #define MPI_Comm_split(comm,color,key,newcomm) MPI_SUCCESS
637: #define MPI_Comm_test_inter(comm,flag) (*(flag)=1,MPI_SUCCESS)
638: #define MPI_Comm_remote_size(comm,size) (*(size)=1,MPI_SUCCESS)
639: #define MPI_Comm_remote_group(comm,group) MPI_SUCCESS
640: #define MPI_Intercomm_create(local_comm,local_leader,peer_comm,\
641: remote_leader,tag,newintercomm) MPI_SUCCESS
642: #define MPI_Intercomm_merge(intercomm,high,newintracomm) MPI_SUCCESS
644: #define MPI_Topo_test(comm,status) MPI_SUCCESS
645: #define MPI_Cart_create(comm_old,ndims,dims,periods,\
646: reorder,comm_cart) MPI_SUCCESS
647: #define MPI_Dims_create(nnodes,ndims,dims) MPI_Abort(MPI_COMM_WORLD,0)
648: #define MPI_Graph_create(comm,a,b,c,d,e) MPI_SUCCESS
649: #define MPI_Graphdims_Get(comm,nnodes,nedges) MPI_Abort(MPI_COMM_WORLD,0)
650: #define MPI_Graph_get(comm,a,b,c,d) MPI_Abort(MPI_COMM_WORLD,0)
651: #define MPI_Cartdim_get(comm,ndims) MPI_Abort(MPI_COMM_WORLD,0)
652: #define MPI_Cart_get(comm,maxdims,dims,periods,coords) \
653: MPI_Abort(MPI_COMM_WORLD,0)
654: #define MPI_Cart_rank(comm,coords,rank) MPI_Abort(MPI_COMM_WORLD,0)
655: #define MPI_Cart_coords(comm,rank,maxdims,coords) \
656: MPI_Abort(MPI_COMM_WORLD,0)
657: #define MPI_Graph_neighbors_count(comm,rank,nneighbors) \
658: MPI_Abort(MPI_COMM_WORLD,0)
659: #define MPI_Graph_neighbors(comm,rank,maxneighbors,neighbors) \
660: MPI_Abort(MPI_COMM_WORLD,0)
661: #define MPI_Cart_shift(comm,direction,disp,rank_source,rank_dest) \
662: MPI_Abort(MPI_COMM_WORLD,0)
663: #define MPI_Cart_sub(comm,remain_dims,newcomm) MPI_Abort(MPI_COMM_WORLD,0)
664: #define MPI_Cart_map(comm,ndims,dims,periods,newrank) MPI_Abort(MPI_COMM_WORLD,0)
665: #define MPI_Graph_map(comm,a,b,c,d) MPI_Abort(MPI_COMM_WORLD,0)
666: #define MPI_Get_processor_name(name,result_len) \
667: (MPIUNI_Memcpy(name,"localhost",9*sizeof(char)),name[10] = 0,*(result_len) = 10)
668: #define MPI_Errhandler_create(function,errhandler) \
669: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (errhandler),\
670: MPI_SUCCESS)
671: #define MPI_Errhandler_set(comm,errhandler) \
672: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm),\
673: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (errhandler),\
674: MPI_SUCCESS)
675: #define MPI_Errhandler_get(comm,errhandler) MPI_SUCCESS
676: #define MPI_Errhandler_free(errhandler) MPI_SUCCESS
677: #define MPI_Error_string(errorcode,string,result_len) MPI_SUCCESS
678: #define MPI_Error_class(errorcode,errorclass) MPI_SUCCESS
679: #define MPI_Wtick() 1.0
680: #define MPI_Wtime() 0.0
681: #define MPI_Pcontrol(level) MPI_SUCCESS
683: #define MPI_NULL_COPY_FN 0
684: #define MPI_NULL_DELETE_FN 0
686: /* MPI-IO additions */
688: #define MPI_File_open(comm,filename,amode,info,mpi_fh) \
689: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (comm), \
690: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (filename), \
691: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (amode), \
692: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (info), \
693: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (mpi_fh), \
694: MPI_Abort(MPI_COMM_WORLD,0))
696: #define MPI_File_close(mpi_fh) \
697: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (mpi_fh), \
698: MPI_Abort(MPI_COMM_WORLD,0))
700: #define MPI_File_set_view(mpi_fh,disp,etype,filetype,datarep,info) \
701: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (mpi_fh), \
702: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (disp), \
703: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (etype), \
704: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (filetype), \
705: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (datarep), \
706: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (info), \
707: MPI_Abort(MPI_COMM_WORLD,0))
709: #define MPI_Type_get_extent(datatype,lb,extent) \
710: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (datatype), \
711: *(lb) = 0, *(extent) = datatype,0)
713: #define MPI_File_write_all(mpi_fh,buf,count,datatype,status) \
714: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (mpi_fh), \
715: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (buf), \
716: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (count), \
717: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (datatype), \
718: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (status), \
719: MPI_Abort(MPI_COMM_WORLD,0))
721: #define MPI_File_read_all(mpi_fh,buf,count,datatype,status) \
722: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (mpi_fh), \
723: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (buf), \
724: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (count), \
725: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (datatype), \
726: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (status), \
727: MPI_Abort(MPI_COMM_WORLD,0))
729: /* called from PetscInitialize() - so return success */
730: #define MPI_Register_datarep(name,read_conv_fn,write_conv_fn,extent_fn,state) \
731: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (name), \
732: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (read_conv_fn), \
733: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (write_conv_fn), \
734: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (extent_fn), \
735: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (state), \
736: MPI_SUCCESS)
738: #define MPI_Type_create_subarray(ndims,array_of_sizes,array_of_subsizes,array_of_starts,order,oldtype,newtype) \
739: (MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (ndims), \
740: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (array_of_sizes), \
741: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (array_of_subsizes), \
742: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (array_of_starts), \
743: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (order), \
744: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (oldtype), \
745: MPIUNI_TMP = (void*)(MPIUNI_INTPTR) (newtype), \
746: MPI_Abort(MPI_COMM_WORLD,0))
748: #if defined(__cplusplus)
749: }
750: #endif
751: #endif