Date: 2022-03-31T04:45:54+00:00

Actual source code: mpi.h

      1: /*
      2:    This is a special set of bindings for uni-processor use of MPI by the PETSc library.

      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 usable MPI
     20:   implementation]

     22:   However - providing a sequential, 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.

     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: */

     94: #if !defined(MPIUNI_H)
     95: #define MPIUNI_H

     97: /* Required by abort() in mpi.c & for win64 */
     98: #include <petscconf.h>
     99: #include <stddef.h>

    101: /*  This is reproduced from petscsys.h so that mpi.h can be used standalone without first including petscsys.h */
    102: #if defined(_WIN32) && defined(PETSC_USE_SHARED_LIBRARIES)
    103: #  define MPIUni_ __declspec(dllexport)
    104: #  define MPIUni_PETSC_DLLIMPORT __declspec(dllimport)
    105: #elif defined(PETSC_USE_VISIBILITY_CXX) && defined(__cplusplus)
    106: #  define MPIUni_ __attribute__((visibility ("default")))
    107: #  define MPIUni_PETSC_DLLIMPORT __attribute__((visibility ("default")))
    108: #elif defined(PETSC_USE_VISIBILITY_C) && !defined(__cplusplus)
    109: #  define MPIUni_ __attribute__((visibility ("default")))
    110: #  define MPIUni_PETSC_DLLIMPORT __attribute__((visibility ("default")))
    111: #else
    112: #  define MPIUni_
    113: #  define MPIUni_PETSC_DLLIMPORT
    114: #endif

    116: #if defined(petsc_EXPORTS)
    117: #  define MPIUni_PETSC_VISIBILITY_PUBLIC MPIUni_
    118: #else  /* Win32 users need this to import symbols from petsc.dll */
    119: #  define MPIUni_PETSC_VISIBILITY_PUBLIC MPIUni_PETSC_DLLIMPORT
    120: #endif

    122: #if defined(__cplusplus)
    123: #define MPIUni_PETSC_EXTERN extern "C" MPIUni_PETSC_VISIBILITY_PUBLIC
    124: #else
    125: #define MPIUni_PETSC_EXTERN extern MPIUni_PETSC_VISIBILITY_PUBLIC
    126: #endif

    128: #if defined(__cplusplus)
    129: extern "C" {
    130: #endif

    132: /* MPI_Aint has to be a signed integral type large enough to hold a pointer */
    133: typedef ptrdiff_t MPI_Aint;

    135: /* old 32bit MS compiler does not support long long */
    136: #if defined(PETSC_SIZEOF_LONG_LONG)
    137: typedef long long MPIUNI_INT64;
    138: typedef unsigned long long MPIUNI_UINT64;
    139: #elif defined(PETSC_HAVE___INT64)
    140: typedef _int64 MPIUNI_INT64;
    141: typedef unsigned _int64 MPIUNI_UINT64;
    142: #else
    143: #error "cannot determine MPIUNI_INT64, MPIUNI_UINT64 types"
    144: #endif

    146: /*

    148:  MPIUNI_ARG is used in the macros below only to stop various C/C++ compilers
    149:  from generating warning messages about unused variables while compiling PETSc.
    150: */
    151: MPIUni_PETSC_EXTERN void *MPIUNI_TMP;
    152: #define MPIUNI_ARG(arg) (MPIUNI_TMP = (void *)(MPI_Aint) (arg))

    154: #define MPI_IDENT            0
    155: #define MPI_CONGRUENT        1
    156: #define MPI_SIMILAR          2
    157: #define MPI_UNEQUAL          3

    159: #define MPI_BOTTOM   ((void *) 0)
    160: #define MPI_IN_PLACE ((void *)-1)

    162: #define MPI_PROC_NULL      (-1)
    163: #define MPI_ANY_SOURCE     (-2)
    164: #define MPI_ANY_TAG        (-1)
    165: #define MPI_UNDEFINED  (-32766)

    167: #define MPI_SUCCESS          0
    168: #define MPI_ERR_COUNT        2
    169: #define MPI_ERR_OTHER       17
    170: #define MPI_ERR_UNKNOWN     18
    171: #define MPI_ERR_INTERN      21
    172: #define MPI_ERR_NOSUPPORT   22

    174: #define MPI_KEYVAL_INVALID   0
    175: #define MPI_TAG_UB           0

    177: #define MPI_MAX_PROCESSOR_NAME 1024
    178: #define MPI_MAX_ERROR_STRING   2056
    179: #define MPI_MAX_OBJECT_NAME    1024

    181: typedef int MPI_Comm;
    182: #define MPI_COMM_NULL  0
    183: #define MPI_COMM_SELF  1
    184: #define MPI_COMM_WORLD 2
    185: #define MPI_COMM_TYPE_SHARED 1

    187: typedef int MPI_Info;
    188: #define MPI_INFO_NULL 0

    190: typedef struct {int MPI_SOURCE,MPI_TAG,MPI_ERROR;} MPI_Status;
    191: #define MPI_STATUS_IGNORE   (MPI_Status *)0
    192: #define MPI_STATUSES_IGNORE (MPI_Status *)0

    194: /* 32-bit packing scheme: [combiner:4 | type-index:8 | count:12 | base-bytes:8] */
    195: /* Any changes here must also be reflected in mpif.h */
    196: typedef int MPI_Datatype;
    197: #define MPI_DATATYPE_NULL      0
    198: #define MPI_PACKED             0

    200: #define MPI_FLOAT              (1 << 20 | 1 << 8 | (int)sizeof(float))
    201: #define MPI_DOUBLE             (1 << 20 | 1 << 8 | (int)sizeof(double))
    202: #define MPI_LONG_DOUBLE        (1 << 20 | 1 << 8 | (int)sizeof(long double))

    204: #define MPI_COMPLEX            (2 << 20 | 1 << 8 | 2*(int)sizeof(float))
    205: #define MPI_C_COMPLEX          (2 << 20 | 1 << 8 | 2*(int)sizeof(float))
    206: #define MPI_C_FLOAT_COMPLEX    (2 << 20 | 1 << 8 | 2*(int)sizeof(float))
    207: #define MPI_DOUBLE_COMPLEX     (2 << 20 | 1 << 8 | 2*(int)sizeof(double))
    208: #define MPI_C_DOUBLE_COMPLEX   (2 << 20 | 1 << 8 | 2*(int)sizeof(double))

    210: #define MPI_CHAR               (3 << 20 | 1 << 8 | (int)sizeof(char))
    211: #define MPI_BYTE               (3 << 20 | 1 << 8 | (int)sizeof(char))
    212: #define MPI_SIGNED_CHAR        (3 << 20 | 1 << 8 | (int)sizeof(signed char))
    213: #define MPI_UNSIGNED_CHAR      (3 << 20 | 1 << 8 | (int)sizeof(unsigned char))

    215: #define MPI_SHORT              (4 << 20 | 1 << 8 | (int)sizeof(short))
    216: #define MPI_INT                (4 << 20 | 1 << 8 | (int)sizeof(int))
    217: #define MPI_LONG               (4 << 20 | 1 << 8 | (int)sizeof(long))
    218: #define MPI_LONG_LONG          (4 << 20 | 1 << 8 | (int)sizeof(MPIUNI_INT64))
    219: #define MPI_LONG_LONG_INT      MPI_LONG_LONG
    220: #define MPI_INTEGER8           MPI_LONG_LONG
    221: #define MPI_INT8_T             (5 << 20 | 1 << 8 | (int)sizeof(int8_t))
    222: #define MPI_INT16_T            (5 << 20 | 1 << 8 | (int)sizeof(int16_t))
    223: #define MPI_INT32_T            (5 << 20 | 1 << 8 | (int)sizeof(int32_t))

    225: #define MPI_UNSIGNED_SHORT     (5 << 20 | 1 << 8 | (int)sizeof(unsigned short))
    226: #define MPI_UNSIGNED           (5 << 20 | 1 << 8 | (int)sizeof(unsigned))
    227: #define MPI_UNSIGNED_LONG      (5 << 20 | 1 << 8 | (int)sizeof(unsigned long))
    228: #define MPI_UNSIGNED_LONG_LONG (5 << 20 | 1 << 8 | (int)sizeof(MPIUNI_UINT64))

    230: #define MPI_FLOAT_INT          (10 << 20 | 1 << 8 | (int)(sizeof(float) + sizeof(int)))
    231: #define MPI_DOUBLE_INT         (11 << 20 | 1 << 8 | (int)(sizeof(double) + sizeof(int)))
    232: #define MPI_LONG_INT           (12 << 20 | 1 << 8 | (int)(sizeof(long) + sizeof(int)))
    233: #define MPI_SHORT_INT          (13 << 20 | 1 << 8 | (int)(sizeof(short) + sizeof(int)))
    234: #define MPI_2INT               (14 << 20 | 1 << 8 | (int)(2*sizeof(int)))
    235: #define MPI_2DOUBLE            (15 << 20 | 1 << 8 | (int)(2*sizeof(double)))

    237: /* Fortran datatypes; Jed Brown says they should be defined here */
    238: #define MPI_INTEGER MPI_INT
    239: #define MPI_DOUBLE_PRECISION MPI_DOUBLE
    240: #define MPI_COMPLEX16 MPI_C_DOUBLE_COMPLEX
    241: #define MPI_2DOUBLE_PRECISION MPI_2DOUBLE

    243: #define MPI_ORDER_C            0
    244: #define MPI_ORDER_FORTRAN      1

    246: #define MPI_sizeof_default(datatype) ((((datatype) >> 8) & 0xfff) * ((datatype) & 0xff))
    247: #if defined(PETSC_USE_REAL___FP16)
    248: MPIUni_PETSC_EXTERN MPI_Datatype MPIU___FP16;
    249: #define MPI_sizeof(datatype) ((datatype == MPIU___FP16) ? (int)(2*sizeof(char)) : MPI_sizeof_default(datatype))
    250: #elif defined(PETSC_USE_REAL___FLOAT128)
    251: MPIUni_PETSC_EXTERN MPI_Datatype MPIU___FLOAT128;
    252: #define MPI_sizeof(datatype) ((datatype == MPIU___FLOAT128) ? (int)(2*sizeof(double)) : MPI_sizeof_default(datatype))
    253: #else
    254: #define MPI_sizeof(datatype) (MPI_sizeof_default(datatype))
    255: #endif

    257: MPIUni_PETSC_EXTERN int MPIUNI_Memcpy(void*,const void*,int);

    259: typedef int MPI_Request;
    260: #define MPI_REQUEST_NULL 0

    262: typedef int MPI_Group;
    263: #define MPI_GROUP_NULL  0
    264: #define MPI_GROUP_EMPTY 0

    266: typedef int MPI_Op;
    267: #define MPI_OP_NULL    0
    268: #define MPI_SUM        1
    269: #define MPI_MAX        2
    270: #define MPI_MIN        3
    271: #define MPI_REPLACE    4
    272: #define MPI_PROD       5
    273: #define MPI_LAND       6
    274: #define MPI_BAND       7
    275: #define MPI_LOR        8
    276: #define MPI_BOR        9
    277: #define MPI_LXOR       10
    278: #define MPI_BXOR       11
    279: #define MPI_MAXLOC     12
    280: #define MPI_MINLOC     13

    282: typedef void (MPI_User_function)(void*, void *, int *, MPI_Datatype *);

    284: typedef int MPI_Errhandler;
    285: #define MPI_ERRHANDLER_NULL  0
    286: #define MPI_ERRORS_RETURN    0
    287: #define MPI_ERRORS_ARE_FATAL 0
    288: #define MPI_ERR_LASTCODE     0x3fffffff
    289: typedef void (MPI_Handler_function)(MPI_Comm *, int *, ...);

    291: /*
    292:   Prototypes of some functions which are implemented in mpi.c
    293: */
    294: typedef int (MPI_Copy_function)(MPI_Comm,int,void *,void *,void *,int *);
    295: typedef int (MPI_Delete_function)(MPI_Comm,int,void *,void *);
    296: #define MPI_NULL_COPY_FN   (MPI_Copy_function*)0
    297: #define MPI_NULL_DELETE_FN (MPI_Delete_function*)0

    299: #define MPI_THREAD_SINGLE 0
    300: #define MPI_THREAD_FUNNELED 1
    301: #define MPI_THREAD_SERIALIZED 2
    302: #define MPI_THREAD_MULTIPLE 3

    304: /*
    305:   To enable linking PETSc+MPIUNI with any other package that might have its
    306:   own MPIUNI (equivalent implementation) we need to avoid using 'MPI'
    307:   namespace for MPIUNI functions that go into the petsc library.

    309:   For C functions below (that get compiled into petsc library) - we map
    310:   the 'MPI' functions to use 'Petsc_MPI' namespace.

    312:   With fortran we use similar mapping - thus requiring the use of
    313:   c-preprocessor with mpif.h
    314: */
    315: #define MPI_Abort         Petsc_MPI_Abort
    316: #define MPIUni_Abort      Petsc_MPIUni_Abort
    317: #define MPI_Attr_get      Petsc_MPI_Attr_get
    318: #define MPI_Keyval_free   Petsc_MPI_Keyval_free
    319: #define MPI_Attr_put      Petsc_MPI_Attr_put
    320: #define MPI_Attr_delete   Petsc_MPI_Attr_delete
    321: #define MPI_Keyval_create Petsc_MPI_Keyval_create
    322: #define MPI_Comm_free     Petsc_MPI_Comm_free
    323: #define MPI_Comm_dup      Petsc_MPI_Comm_dup
    324: #define MPI_Comm_create   Petsc_MPI_Comm_create
    325: #define MPI_Init          Petsc_MPI_Init
    326: #define MPI_Init_thread   Petsc_MPI_Init_thread
    327: #define MPI_Query_thread  Petsc_MPI_Query_thread
    328: #define MPI_Finalize      Petsc_MPI_Finalize
    329: #define MPI_Initialized   Petsc_MPI_Initialized
    330: #define MPI_Finalized     Petsc_MPI_Finalized
    331: #define MPI_Comm_size     Petsc_MPI_Comm_size
    332: #define MPI_Comm_rank     Petsc_MPI_Comm_rank
    333: #define MPI_Wtime         Petsc_MPI_Wtime
    334: #define MPI_Type_get_envelope Petsc_MPI_Type_get_envelope
    335: #define MPI_Type_get_contents Petsc_MPI_Type_get_contents
    336: #define MPI_Add_error_class   Petsc_MPI_Add_error_class
    337: #define MPI_Add_error_code    Petsc_MPI_Add_error_code

    339: /* identical C bindings */
    340: #define MPI_Comm_copy_attr_function   MPI_Copy_function
    341: #define MPI_Comm_delete_attr_function MPI_Delete_function
    342: #define MPI_COMM_NULL_COPY_FN         MPI_NULL_COPY_FN
    343: #define MPI_COMM_NULL_DELETE_FN       MPI_NULL_DELETE_FN
    344: #define MPI_Comm_create_keyval        Petsc_MPI_Keyval_create
    345: #define MPI_Comm_free_keyval          Petsc_MPI_Keyval_free
    346: #define MPI_Comm_get_attr             Petsc_MPI_Attr_get
    347: #define MPI_Comm_set_attr             Petsc_MPI_Attr_put
    348: #define MPI_Comm_delete_attr          Petsc_MPI_Attr_delete
    349: #define MPI_Comm_get_name             Petsc_MPI_Comm_get_name
    350: #define MPI_Comm_set_name             Petsc_MPI_Comm_set_name

    352: MPIUni_PETSC_EXTERN int    MPIUni_Abort(MPI_Comm,int);
    353: MPIUni_PETSC_EXTERN int    MPI_Abort(MPI_Comm,int);
    354: MPIUni_PETSC_EXTERN int    MPI_Attr_get(MPI_Comm comm,int keyval,void *attribute_val,int *flag);
    355: MPIUni_PETSC_EXTERN int    MPI_Keyval_free(int*);
    356: MPIUni_PETSC_EXTERN int    MPI_Attr_put(MPI_Comm,int,void *);
    357: MPIUni_PETSC_EXTERN int    MPI_Attr_delete(MPI_Comm,int);
    358: MPIUni_PETSC_EXTERN int    MPI_Keyval_create(MPI_Copy_function *,MPI_Delete_function *,int *,void *);
    359: MPIUni_PETSC_EXTERN int    MPI_Comm_free(MPI_Comm*);
    360: MPIUni_PETSC_EXTERN int    MPI_Comm_dup(MPI_Comm,MPI_Comm *);
    361: MPIUni_PETSC_EXTERN int    MPI_Comm_create(MPI_Comm,MPI_Group,MPI_Comm *);
    362: MPIUni_PETSC_EXTERN int    MPI_Init(int *, char ***);
    363: MPIUni_PETSC_EXTERN int    MPI_Init_thread(int *, char ***, int, int *);
    364: MPIUni_PETSC_EXTERN int    MPI_Query_thread(int *);
    365: MPIUni_PETSC_EXTERN int    MPI_Finalize(void);
    366: MPIUni_PETSC_EXTERN int    MPI_Initialized(int*);
    367: MPIUni_PETSC_EXTERN int    MPI_Finalized(int*);
    368: MPIUni_PETSC_EXTERN int    MPI_Comm_size(MPI_Comm,int*);
    369: MPIUni_PETSC_EXTERN int    MPI_Comm_rank(MPI_Comm,int*);
    370: MPIUni_PETSC_EXTERN double MPI_Wtime(void);

    372: MPIUni_PETSC_EXTERN int MPI_Type_get_envelope(MPI_Datatype,int*,int*,int*,int*);
    373: MPIUni_PETSC_EXTERN int MPI_Type_get_contents(MPI_Datatype,int,int,int,int*,MPI_Aint*,MPI_Datatype*);
    374: MPIUni_PETSC_EXTERN int MPI_Add_error_class(int*);
    375: MPIUni_PETSC_EXTERN int MPI_Add_error_code(int,int*);
    376: MPIUni_PETSC_EXTERN int MPI_Comm_get_name(MPI_Comm,char*,int*);
    377: MPIUni_PETSC_EXTERN int MPI_Comm_set_name(MPI_Comm,const char*);
    378: /*
    379:     Routines we have replace with macros that do nothing
    380:     Some return error codes others return success
    381: */

    383: typedef int MPI_Fint;
    384: #define MPI_Comm_f2c(comm) (MPI_Comm)(comm)
    385: #define MPI_Comm_c2f(comm) (MPI_Fint)(comm)
    386: #define MPI_Type_f2c(type) (MPI_Datatype)(type)
    387: #define MPI_Type_c2f(type) (MPI_Fint)(type)
    388: #define MPI_Op_f2c(op)     (MPI_Op)(op)
    389: #define MPI_Op_c2f(op)     (MPI_Fint)(op)

    391: #define MPI_Send(buf,count,datatype,dest,tag,comm)  \
    392:      (MPIUNI_ARG(buf),\
    393:       MPIUNI_ARG(count),\
    394:       MPIUNI_ARG(datatype),\
    395:       MPIUNI_ARG(dest),\
    396:       MPIUNI_ARG(tag),\
    397:       MPIUNI_ARG(comm),\
    398:       MPIUni_Abort(MPI_COMM_WORLD,0))
    399: #define MPI_Recv(buf,count,datatype,source,tag,comm,status) \
    400:      (MPIUNI_ARG(buf),\
    401:       MPIUNI_ARG(count),\
    402:       MPIUNI_ARG(datatype),\
    403:       MPIUNI_ARG(source),\
    404:       MPIUNI_ARG(tag),\
    405:       MPIUNI_ARG(comm),\
    406:       MPIUNI_ARG(status),\
    407:       MPIUni_Abort(MPI_COMM_WORLD,0))
    408: #define MPI_Get_count(status,datatype,count) \
    409:      (MPIUNI_ARG(status),\
    410:       MPIUNI_ARG(datatype),\
    411:       MPIUNI_ARG(count),\
    412:       MPIUni_Abort(MPI_COMM_WORLD,0))
    413: #define MPI_Bsend(buf,count,datatype,dest,tag,comm)  \
    414:      (MPIUNI_ARG(buf),\
    415:       MPIUNI_ARG(count),\
    416:       MPIUNI_ARG(datatype),\
    417:       MPIUNI_ARG(dest),\
    418:       MPIUNI_ARG(tag),\
    419:       MPIUNI_ARG(comm),\
    420:       MPIUni_Abort(MPI_COMM_WORLD,0))
    421: #define MPI_Ssend(buf,count,datatype,dest,tag,comm) \
    422:      (MPIUNI_ARG(buf),\
    423:       MPIUNI_ARG(count),\
    424:       MPIUNI_ARG(datatype),\
    425:       MPIUNI_ARG(dest),\
    426:       MPIUNI_ARG(tag),\
    427:       MPIUNI_ARG(comm),\
    428:       MPIUni_Abort(MPI_COMM_WORLD,0))
    429: #define MPI_Rsend(buf,count,datatype,dest,tag,comm) \
    430:      (MPIUNI_ARG(buf),\
    431:       MPIUNI_ARG(count),\
    432:       MPIUNI_ARG(datatype),\
    433:       MPIUNI_ARG(dest),\
    434:       MPIUNI_ARG(tag),\
    435:       MPIUNI_ARG(comm),\
    436:       MPIUni_Abort(MPI_COMM_WORLD,0))
    437: #define MPI_Buffer_attach(buffer,size) \
    438:      (MPIUNI_ARG(buffer),\
    439:       MPIUNI_ARG(size),\
    440:       MPI_SUCCESS)
    441: #define MPI_Buffer_detach(buffer,size)\
    442:      (MPIUNI_ARG(buffer),\
    443:       MPIUNI_ARG(size),\
    444:       MPI_SUCCESS)
    445: #define MPI_Ibsend(buf,count,datatype,dest,tag,comm,request) \
    446:      (MPIUNI_ARG(buf),\
    447:       MPIUNI_ARG(count),\
    448:       MPIUNI_ARG(datatype),\
    449:       MPIUNI_ARG(dest),\
    450:       MPIUNI_ARG(tag),\
    451:       MPIUNI_ARG(comm),\
    452:       MPIUNI_ARG(request),\
    453:       MPIUni_Abort(MPI_COMM_WORLD,0))
    454: #define MPI_Issend(buf,count,datatype,dest,tag,comm,request) \
    455:      (MPIUNI_ARG(buf),\
    456:       MPIUNI_ARG(count),\
    457:       MPIUNI_ARG(datatype),\
    458:       MPIUNI_ARG(dest),\
    459:       MPIUNI_ARG(tag),\
    460:       MPIUNI_ARG(comm),\
    461:       MPIUNI_ARG(request),\
    462:       MPIUni_Abort(MPI_COMM_WORLD,0))
    463: #define MPI_Irsend(buf,count,datatype,dest,tag,comm,request) \
    464:      (MPIUNI_ARG(buf),\
    465:       MPIUNI_ARG(count),\
    466:       MPIUNI_ARG(datatype),\
    467:       MPIUNI_ARG(dest),\
    468:       MPIUNI_ARG(tag),\
    469:       MPIUNI_ARG(comm),\
    470:       MPIUNI_ARG(request),\
    471:       MPIUni_Abort(MPI_COMM_WORLD,0))
    472: #define MPI_Irecv(buf,count,datatype,source,tag,comm,request) \
    473:      (MPIUNI_ARG(buf),\
    474:       MPIUNI_ARG(count),\
    475:       MPIUNI_ARG(datatype),\
    476:       MPIUNI_ARG(source),\
    477:       MPIUNI_ARG(tag),\
    478:       MPIUNI_ARG(comm),\
    479:       MPIUNI_ARG(request),\
    480:       MPIUni_Abort(MPI_COMM_WORLD,0))
    481: #define MPI_Isend(buf,count,datatype,dest,tag,comm,request) \
    482:      (MPIUNI_ARG(buf),\
    483:       MPIUNI_ARG(count),\
    484:       MPIUNI_ARG(datatype),\
    485:       MPIUNI_ARG(dest),\
    486:       MPIUNI_ARG(tag),\
    487:       MPIUNI_ARG(comm),\
    488:       MPIUNI_ARG(request),\
    489:       MPIUni_Abort(MPI_COMM_WORLD,0))
    490: #define MPI_Wait(request,status) \
    491:      (MPIUNI_ARG(request),\
    492:       MPIUNI_ARG(status),\
    493:       MPI_SUCCESS)
    494: #define MPI_Test(request,flag,status) \
    495:      (MPIUNI_ARG(request),\
    496:       MPIUNI_ARG(status),\
    497:       *(flag) = 0,\
    498:       MPI_SUCCESS)
    499: #define MPI_Request_free(request) \
    500:      (MPIUNI_ARG(request),\
    501:       MPI_SUCCESS)
    502: #define MPI_Waitany(count,array_of_requests,index,status) \
    503:      (MPIUNI_ARG(count),\
    504:       MPIUNI_ARG(array_of_requests),\
    505:       MPIUNI_ARG(status),\
    506:       (*(status)).MPI_SOURCE = 0,               \
    507:       *(index) = 0,\
    508:       MPI_SUCCESS)
    509: #define MPI_Testany(a,b,c,d,e) \
    510:      (MPIUNI_ARG(a),\
    511:       MPIUNI_ARG(b),\
    512:       MPIUNI_ARG(c),\
    513:       MPIUNI_ARG(d),\
    514:       MPIUNI_ARG(e),\
    515:       MPI_SUCCESS)
    516: #define MPI_Waitall(count,array_of_requests,array_of_statuses) \
    517:      (MPIUNI_ARG(count),\
    518:       MPIUNI_ARG(array_of_requests),\
    519:       MPIUNI_ARG(array_of_statuses),\
    520:       MPI_SUCCESS)
    521: #define MPI_Testall(count,array_of_requests,flag,array_of_statuses) \
    522:      (MPIUNI_ARG(count),\
    523:       MPIUNI_ARG(array_of_requests),\
    524:       MPIUNI_ARG(flag),\
    525:       MPIUNI_ARG(array_of_statuses),\
    526:       MPI_SUCCESS)
    527: #define MPI_Waitsome(incount,array_of_requests,outcount,\
    528:                      array_of_indices,array_of_statuses)        \
    529:      (MPIUNI_ARG(incount),\
    530:       MPIUNI_ARG(array_of_requests),\
    531:       MPIUNI_ARG(outcount),\
    532:       MPIUNI_ARG(array_of_indices),\
    533:       MPIUNI_ARG(array_of_statuses),\
    534:       MPI_SUCCESS)
    535: #define MPI_Comm_group(comm,group) \
    536:      (MPIUNI_ARG(comm),\
    537:       *group = 1,\
    538:       MPI_SUCCESS)
    539: #define MPI_Group_excl(group,n,ranks,newgroup) \
    540:      (MPIUNI_ARG(group),\
    541:       MPIUNI_ARG(n),\
    542:       MPIUNI_ARG(ranks),\
    543:       MPIUNI_ARG(newgroup),\
    544:       MPI_SUCCESS)
    545: #define MPI_Group_incl(group,n,ranks,newgroup) \
    546:      (MPIUNI_ARG(group),\
    547:       MPIUNI_ARG(n),\
    548:       MPIUNI_ARG(ranks),\
    549:       MPIUNI_ARG(newgroup),\
    550:       MPI_SUCCESS)
    551: #define MPI_Testsome(incount,array_of_requests,outcount,\
    552:                      array_of_indices,array_of_statuses) \
    553:      (MPIUNI_ARG(incount),\
    554:       MPIUNI_ARG(array_of_requests),\
    555:       MPIUNI_ARG(outcount),\
    556:       MPIUNI_ARG(array_of_indices),\
    557:       MPIUNI_ARG(array_of_statuses),\
    558:       MPI_SUCCESS)
    559: #define MPI_Iprobe(source,tag,comm,flag,status) \
    560:      (MPIUNI_ARG(source),\
    561:       MPIUNI_ARG(tag),\
    562:       MPIUNI_ARG(comm),\
    563:       *(flag)=0,\
    564:       MPIUNI_ARG(status),\
    565:       MPI_SUCCESS)
    566: #define MPI_Probe(source,tag,comm,status) \
    567:      (MPIUNI_ARG(source),\
    568:       MPIUNI_ARG(tag),\
    569:       MPIUNI_ARG(comm),\
    570:       MPIUNI_ARG(status),\
    571:       MPI_SUCCESS)
    572: #define MPI_Cancel(request) \
    573:      (MPIUNI_ARG(request),\
    574:       MPI_SUCCESS)
    575: #define MPI_Test_cancelled(status,flag) \
    576:      (MPIUNI_ARG(status),\
    577:       *(flag)=0,\
    578:       MPI_SUCCESS)
    579: #define MPI_Send_init(buf,count,datatype,dest,tag,comm,request) \
    580:      (MPIUNI_ARG(buf),\
    581:       MPIUNI_ARG(count),\
    582:       MPIUNI_ARG(datatype),\
    583:       MPIUNI_ARG(dest),\
    584:       MPIUNI_ARG(tag),\
    585:       MPIUNI_ARG(comm),\
    586:       MPIUNI_ARG(request),\
    587:       MPI_SUCCESS)
    588: #define MPI_Bsend_init(buf,count,datatype,dest,tag,comm,request) \
    589:      (MPIUNI_ARG(buf),\
    590:       MPIUNI_ARG(count),\
    591:       MPIUNI_ARG(datatype),\
    592:       MPIUNI_ARG(dest),\
    593:       MPIUNI_ARG(tag),\
    594:       MPIUNI_ARG(comm),\
    595:       MPIUNI_ARG(request),\
    596:       MPI_SUCCESS)
    597: #define MPI_Ssend_init(buf,count,datatype,dest,tag,comm,request) \
    598:      (MPIUNI_ARG(buf),\
    599:       MPIUNI_ARG(count),\
    600:       MPIUNI_ARG(datatype),\
    601:       MPIUNI_ARG(dest),\
    602:       MPIUNI_ARG(tag),\
    603:       MPIUNI_ARG(comm),\
    604:       MPIUNI_ARG(request),\
    605:       MPI_SUCCESS)
    606: #define MPI_Bsend_init(buf,count,datatype,dest,tag,comm,request) \
    607:      (MPIUNI_ARG(buf),\
    608:       MPIUNI_ARG(count),\
    609:       MPIUNI_ARG(datatype),\
    610:       MPIUNI_ARG(dest),\
    611:       MPIUNI_ARG(tag),\
    612:       MPIUNI_ARG(comm),\
    613:       MPIUNI_ARG(request),\
    614:       MPI_SUCCESS)
    615: #define MPI_Rsend_init(buf,count,datatype,dest,tag,comm,request) \
    616:      (MPIUNI_ARG(buf),\
    617:       MPIUNI_ARG(count),\
    618:       MPIUNI_ARG(datatype),\
    619:       MPIUNI_ARG(dest),\
    620:       MPIUNI_ARG(tag),\
    621:       MPIUNI_ARG(comm),\
    622:       MPIUNI_ARG(request),\
    623:       MPI_SUCCESS)
    624: #define MPI_Recv_init(buf,count,datatype,source,tag,comm,request) \
    625:      (MPIUNI_ARG(buf),\
    626:       MPIUNI_ARG(count),\
    627:       MPIUNI_ARG(datatype),\
    628:       MPIUNI_ARG(source),\
    629:       MPIUNI_ARG(tag),\
    630:       MPIUNI_ARG(comm),\
    631:       MPIUNI_ARG(request),\
    632:       MPI_SUCCESS)
    633: #define MPI_Start(request) \
    634:      (MPIUNI_ARG(request),\
    635:       MPI_SUCCESS)
    636: #define MPI_Startall(count,array_of_requests) \
    637:      (MPIUNI_ARG(count),\
    638:       MPIUNI_ARG(array_of_requests),\
    639:       MPI_SUCCESS)
    640: #define MPI_Sendrecv(sendbuf,sendcount,sendtype,\
    641:                      dest,sendtag,recvbuf,recvcount,\
    642:                      recvtype,source,recvtag,\
    643:                      comm,status) \
    644:      (MPIUNI_ARG(dest),\
    645:       MPIUNI_ARG(sendtag),\
    646:       MPIUNI_ARG(recvcount),\
    647:       MPIUNI_ARG(recvtype),\
    648:       MPIUNI_ARG(source),\
    649:       MPIUNI_ARG(recvtag),\
    650:       MPIUNI_ARG(comm),\
    651:       MPIUNI_ARG(status),\
    652:       MPIUNI_Memcpy(recvbuf,sendbuf,(sendcount)*MPI_sizeof(sendtype)))
    653: #define MPI_Sendrecv_replace(buf,count,datatype,dest,sendtag,\
    654:                              source,recvtag,comm,status) \
    655:      (MPIUNI_ARG(buf),\
    656:       MPIUNI_ARG(count),\
    657:       MPIUNI_ARG(datatype),\
    658:       MPIUNI_ARG(dest),\
    659:       MPIUNI_ARG(sendtag),\
    660:       MPIUNI_ARG(source),\
    661:       MPIUNI_ARG(recvtag),\
    662:       MPIUNI_ARG(comm),\
    663:       MPIUNI_ARG(status),\
    664:       MPI_SUCCESS)

    666: #define MPI_COMBINER_NAMED      0
    667: #define MPI_COMBINER_DUP        1
    668: #define MPI_COMBINER_CONTIGUOUS 2
    669:   /* 32-bit packing scheme: [combiner:4 | type-index:8 | count:12 | base-bytes:8] */
    670: #define MPI_Type_dup(oldtype,newtype) \
    671:      (*(newtype) = oldtype, MPI_SUCCESS)
    672: #define MPI_Type_contiguous(count,oldtype,newtype) \
    673:      (*(newtype) = (MPI_COMBINER_CONTIGUOUS<<28)|((oldtype)&0x0ff00000)|(((oldtype)>>8&0xfff)*(count))<<8|((oldtype)&0xff), MPI_SUCCESS)
    674: #define MPI_Type_vector(count,blocklength,stride,oldtype,newtype) \
    675:      (MPIUNI_ARG(count),\
    676:       MPIUNI_ARG(blocklength),\
    677:       MPIUNI_ARG(stride),\
    678:       MPIUNI_ARG(oldtype),\
    679:       MPIUNI_ARG(newtype),\
    680:       MPIUni_Abort(MPI_COMM_WORLD,0))
    681: #define MPI_Type_hvector(count,blocklength,stride,oldtype,newtype) \
    682:      (MPIUNI_ARG(count),\
    683:       MPIUNI_ARG(blocklength),\
    684:       MPIUNI_ARG(stride),\
    685:       MPIUNI_ARG(oldtype),\
    686:       MPIUNI_ARG(newtype),\
    687:       MPIUni_Abort(MPI_COMM_WORLD,0))
    688: #define MPI_Type_indexed(count,array_of_blocklengths,array_of_displacements,oldtype,newtype) \
    689:      (MPIUNI_ARG(count),\
    690:       MPIUNI_ARG(array_of_blocklengths),\
    691:       MPIUNI_ARG(array_of_displacements),\
    692:       MPIUNI_ARG(oldtype),\
    693:       MPIUNI_ARG(newtype),\
    694:       MPIUni_Abort(MPI_COMM_WORLD,0))
    695: #define MPI_Type_hindexed(count,array_of_blocklengths,array_of_displacements,oldtype,newtype) \
    696:      (MPIUNI_ARG(count),\
    697:       MPIUNI_ARG(array_of_blocklengths),\
    698:       MPIUNI_ARG(array_of_displacements),\
    699:       MPIUNI_ARG(oldtype),\
    700:       MPIUNI_ARG(newtype),\
    701:       MPIUni_Abort(MPI_COMM_WORLD,0))
    702: #define MPI_Type_struct(count,array_of_blocklengths,array_of_displacements,array_of_types,newtype) \
    703:      (MPIUNI_ARG(count),\
    704:       MPIUNI_ARG(array_of_blocklengths),\
    705:       MPIUNI_ARG(array_of_displacements),\
    706:       MPIUNI_ARG(array_of_types),\
    707:       MPIUNI_ARG(newtype),\
    708:       MPIUni_Abort(MPI_COMM_WORLD,0))
    709: #define MPI_Address(location,address) \
    710:      (*(address) = (MPI_Aint)((char *)(location)), MPI_SUCCESS)
    711: #define MPI_Type_size(datatype,size) (*(size) = MPI_sizeof((datatype)), MPI_SUCCESS)
    712: #define MPI_Type_lb(datatype,lb) (MPIUNI_ARG(datatype), *(lb) = 0, MPI_SUCCESS)
    713: #define MPI_Type_ub(datatype,ub) (*(ub) = MPI_sizeof((datatype)), MPI_SUCCESS)
    714: #define MPI_Type_extent(datatype,extent) \
    715:      (*(extent) = MPI_sizeof((datatype)), MPI_SUCCESS)
    716: #define MPI_Type_get_extent(datatype,lb,extent) \
    717:      (*(lb) = 0, *(extent) = MPI_sizeof((datatype)), MPI_SUCCESS)
    718: #define MPI_Type_commit(datatype) (MPIUNI_ARG(datatype), MPI_SUCCESS)
    719: #define MPI_Type_free(datatype) (*(datatype) = MPI_DATATYPE_NULL, MPI_SUCCESS)
    720: #define MPI_Get_elements(status,datatype,count) \
    721:      (MPIUNI_ARG(status),\
    722:       MPIUNI_ARG(datatype),\
    723:       MPIUNI_ARG(count),\
    724:       MPIUni_Abort(MPI_COMM_WORLD,0))
    725: #define MPI_Pack(inbuf,incount,datatype,outbuf,outsize,position,comm) \
    726:      (MPIUNI_ARG(inbuf),\
    727:       MPIUNI_ARG(incount),\
    728:       MPIUNI_ARG(datatype),\
    729:       MPIUNI_ARG(outbuf),\
    730:       MPIUNI_ARG(outsize),\
    731:       MPIUNI_ARG(position),\
    732:       MPIUNI_ARG(comm),\
    733:       MPIUni_Abort(MPI_COMM_WORLD,0))
    734: #define MPI_Unpack(inbuf,insize,position,outbuf,outcount,datatype,comm) \
    735:      (MPIUNI_ARG(inbuf),\
    736:       MPIUNI_ARG(insize),\
    737:       MPIUNI_ARG(position),\
    738:       MPIUNI_ARG(outbuf),\
    739:       MPIUNI_ARG(outcount),\
    740:       MPIUNI_ARG(datatype),\
    741:       MPIUNI_ARG(comm),\
    742:       MPIUni_Abort(MPI_COMM_WORLD,0))
    743: #define MPI_Pack_size(incount,datatype,comm,size) \
    744:      (MPIUNI_ARG(incount),\
    745:       MPIUNI_ARG(datatype),\
    746:       MPIUNI_ARG(comm),\
    747:       MPIUNI_ARG(size),\
    748:       MPIUni_Abort(MPI_COMM_WORLD,0))
    749: #define MPI_Barrier(comm) \
    750:      (MPIUNI_ARG(comm),\
    751:       MPI_SUCCESS)
    752: #define MPI_Bcast(buffer,count,datatype,root,comm) \
    753:      (MPIUNI_ARG(buffer),\
    754:       MPIUNI_ARG(count),\
    755:       MPIUNI_ARG(datatype),\
    756:       MPIUNI_ARG(root),\
    757:       MPIUNI_ARG(comm),\
    758:       MPI_SUCCESS)
    759: #define MPI_Gather(sendbuf,sendcount,sendtype,\
    760:                    recvbuf,recvcount, recvtype,\
    761:                    root,comm) \
    762:      (MPIUNI_ARG(recvcount),\
    763:       MPIUNI_ARG(root),\
    764:       MPIUNI_ARG(recvtype),\
    765:       MPIUNI_ARG(comm),\
    766:       MPIUNI_Memcpy(recvbuf,sendbuf,(sendcount)*MPI_sizeof(sendtype)))
    767: #define MPI_Gatherv(sendbuf,sendcount,sendtype,\
    768:                     recvbuf,recvcounts,displs,\
    769:                     recvtype,root,comm) \
    770:      (MPIUNI_ARG(recvcounts),\
    771:       MPIUNI_ARG(displs),\
    772:       MPIUNI_ARG(recvtype),\
    773:       MPIUNI_ARG(root),\
    774:       MPIUNI_ARG(comm),\
    775:       MPIUNI_Memcpy(recvbuf,sendbuf,(sendcount)*MPI_sizeof(sendtype)))
    776: #define MPI_Scatter(sendbuf,sendcount,sendtype,\
    777:                     recvbuf,recvcount,recvtype,\
    778:                     root,comm) \
    779:      (MPIUNI_ARG(sendcount),\
    780:       MPIUNI_ARG(sendtype),\
    781:       MPIUNI_ARG(recvbuf),\
    782:       MPIUNI_ARG(recvtype),\
    783:       MPIUNI_ARG(root),\
    784:       MPIUNI_ARG(comm),\
    785:       MPIUNI_Memcpy(recvbuf,sendbuf,(recvcount)*MPI_sizeof(recvtype)))
    786: #define MPI_Scatterv(sendbuf,sendcounts,displs,\
    787:                      sendtype,recvbuf,recvcount,\
    788:                      recvtype,root,comm) \
    789:      (MPIUNI_ARG(displs),\
    790:       MPIUNI_ARG(sendtype),\
    791:       MPIUNI_ARG(sendcounts),\
    792:       MPIUNI_ARG(root),\
    793:       MPIUNI_ARG(comm),\
    794:       MPIUNI_Memcpy(recvbuf,sendbuf,(recvcount)*MPI_sizeof(recvtype)))
    795: #define MPI_Allgather(sendbuf,sendcount,sendtype,\
    796:                      recvbuf,recvcount,recvtype,comm) \
    797:      (MPIUNI_ARG(recvcount),\
    798:       MPIUNI_ARG(recvtype),\
    799:       MPIUNI_ARG(comm),\
    800:       MPIUNI_Memcpy(recvbuf,sendbuf,(sendcount)*MPI_sizeof(sendtype)))
    801: #define MPI_Allgatherv(sendbuf,sendcount,sendtype,\
    802:      recvbuf,recvcounts,displs,recvtype,comm) \
    803:      (MPIUNI_ARG(recvcounts),\
    804:       MPIUNI_ARG(displs),\
    805:       MPIUNI_ARG(recvtype),\
    806:       MPIUNI_ARG(comm),\
    807:       MPIUNI_Memcpy(recvbuf,sendbuf,(sendcount)*MPI_sizeof(sendtype)))
    808: #define MPI_Alltoall(sendbuf,sendcount,sendtype,\
    809:                      recvbuf,recvcount,recvtype,comm) \
    810:      (MPIUNI_ARG(recvcount),\
    811:       MPIUNI_ARG(recvtype),\
    812:       MPIUNI_ARG(comm),\
    813:       MPIUNI_Memcpy(recvbuf,sendbuf,(sendcount)*MPI_sizeof(sendtype)))
    814: #define MPI_Alltoallv(sendbuf,sendcounts,sdispls,sendtype,\
    815:                       recvbuf,recvcounts,rdispls,recvtype,comm) \
    816:      (MPIUNI_ARG(sendbuf),\
    817:       MPIUNI_ARG(sendcounts),\
    818:       MPIUNI_ARG(sdispls),\
    819:       MPIUNI_ARG(sendtype),\
    820:       MPIUNI_ARG(recvbuf),\
    821:       MPIUNI_ARG(recvcounts),\
    822:       MPIUNI_ARG(rdispls),\
    823:       MPIUNI_ARG(recvtype),\
    824:       MPIUNI_ARG(comm),\
    825:       MPIUni_Abort(MPI_COMM_WORLD,0))
    826: #define MPI_Alltoallw(sendbuf,sendcounts,sdispls,sendtypes,\
    827:                       recvbuf,recvcounts,rdispls,recvtypes,comm) \
    828:      (MPIUNI_ARG(sendbuf),\
    829:       MPIUNI_ARG(sendcounts),\
    830:       MPIUNI_ARG(sdispls),\
    831:       MPIUNI_ARG(sendtypes),\
    832:       MPIUNI_ARG(recvbuf),\
    833:       MPIUNI_ARG(recvcount),\
    834:       MPIUNI_ARG(rdispls),\
    835:       MPIUNI_ARG(recvtypes),\
    836:       MPIUNI_ARG(comm),\
    837:       MPIUni_Abort(MPI_COMM_WORLD,0))
    838: #define MPI_Reduce(sendbuf,recvbuf,count,datatype,op,root,comm) \
    839:      (MPIUNI_ARG(op),\
    840:       MPIUNI_ARG(root),\
    841:       MPIUNI_ARG(comm),\
    842:       MPIUNI_Memcpy(recvbuf,sendbuf,(count)*MPI_sizeof(datatype)))
    843: #define MPI_Allreduce(sendbuf,recvbuf,count,datatype,op,comm) \
    844:      (MPIUNI_ARG(op),\
    845:       MPIUNI_ARG(comm),\
    846:       MPIUNI_Memcpy(recvbuf,sendbuf,(count)*MPI_sizeof(datatype)))
    847: #define MPI_Iallreduce(sendbuf,recvbuf,count,datatype,op,comm,request) \
    848:      (MPIUNI_ARG(op),\
    849:       MPIUNI_ARG(comm),\
    850:       MPIUNI_ARG(request),\
    851:       MPIUNI_Memcpy(recvbuf,sendbuf,(count)*MPI_sizeof(datatype)))
    852: #define MPI_Scan(sendbuf, recvbuf,count,datatype,op,comm) \
    853:      (MPIUNI_ARG(op),\
    854:       MPIUNI_ARG(comm),\
    855:       MPIUNI_Memcpy(recvbuf,sendbuf,(count)*MPI_sizeof(datatype)))
    856: #define MPI_Exscan(sendbuf,recvbuf,count,datatype,op,comm) \
    857:      (MPIUNI_ARG(sendbuf),\
    858:       MPIUNI_ARG(recvbuf),\
    859:       MPIUNI_ARG(count),\
    860:       MPIUNI_ARG(datatype),\
    861:       MPIUNI_ARG(op),\
    862:       MPIUNI_ARG(comm),\
    863:       MPI_SUCCESS)
    864: #define MPI_Reduce_scatter(sendbuf,recvbuf,recvcounts,datatype,op,comm) \
    865:      (MPIUNI_ARG(op),\
    866:       MPIUNI_ARG(comm),\
    867:       MPIUNI_Memcpy(recvbuf,sendbuf,(*recvcounts)*MPI_sizeof(datatype)))
    868: #define MPI_Op_create(function,commute,op) \
    869:      (MPIUNI_ARG(function),\
    870:       MPIUNI_ARG(commute),\
    871:       MPIUNI_ARG(op),\
    872:       MPI_SUCCESS)
    873: #define MPI_Op_free(op) \
    874:      (*(op) = MPI_OP_NULL, MPI_SUCCESS)

    876: #define MPI_Group_size(group,size) \
    877:   (MPIUNI_ARG(group),\
    878:    *(size)=1,\
    879:    MPI_SUCCESS)
    880: #define MPI_Group_rank(group,rank) \
    881:   (MPIUNI_ARG(group),\
    882:    *(rank)=0,\
    883:    MPI_SUCCESS)
    884: #define MPI_Group_translate_ranks(group1,n,ranks1,group2,ranks2) \
    885:      (MPIUNI_ARG(group1),\
    886:       MPIUNI_ARG(group2),\
    887:       MPIUNI_Memcpy((ranks2),(ranks1),(n)*MPI_sizeof(MPI_INT)))
    888: #define MPI_Group_compare(group1,group2,result) \
    889:     (MPIUNI_ARG(group1),\
    890:      MPIUNI_ARG(group2),\
    891:      *(result)=1,\
    892:      MPI_SUCCESS)
    893: #define MPI_Group_union(group1,group2,newgroup) \
    894:     (MPIUNI_ARG(group1),\
    895:      MPIUNI_ARG(group2),\
    896:      *(newgroup)=1,\
    897:      MPI_SUCCESS)
    898: #define MPI_Group_intersection(group1,group2,newgroup) \
    899:     (MPIUNI_ARG(group1),\
    900:      MPIUNI_ARG(group2),\
    901:      *(newgroup)=1,\
    902:      MPI_SUCCESS)
    903: #define MPI_Group_difference(group1,group2,newgroup) \
    904:     (MPIUNI_ARG(group1),\
    905:      MPIUNI_ARG(group2),\
    906:      *(newgroup)=MPI_GROUP_EMPTY,\
    907:      MPI_SUCCESS)
    908: #define MPI_Group_range_incl(group,n,ranges,newgroup) \
    909:     (MPIUNI_ARG(group),\
    910:      MPIUNI_ARG(n),\
    911:      MPIUNI_ARG(ranges),\
    912:      *(newgroup)=1,\
    913:      MPI_SUCCESS)
    914: #define MPI_Group_range_excl(group,n,ranges,newgroup) \
    915:     (MPIUNI_ARG(group),\
    916:      MPIUNI_ARG(n),\
    917:      MPIUNI_ARG(ranges),\
    918:      *(newgroup)=MPI_GROUP_EMPTY,\
    919:      MPI_SUCCESS)
    920: #define MPI_Group_free(group) \
    921:      (*(group) = MPI_GROUP_NULL, MPI_SUCCESS)

    923: #define MPI_Comm_compare(comm1,comm2,result) \
    924:      (MPIUNI_ARG(comm1),\
    925:       MPIUNI_ARG(comm2),\
    926:       *(result)=MPI_IDENT,\
    927:       MPI_SUCCESS)
    928: #define MPI_Comm_split(comm,color,key,newcomm) \
    929:      (MPIUNI_ARG(color),\
    930:       MPIUNI_ARG(key),\
    931:       MPI_Comm_dup(comm,newcomm))
    932: #define MPI_Comm_split_type(comm,color,key,info,newcomm) \
    933:      (MPIUNI_ARG(color),\
    934:       MPIUNI_ARG(key),\
    935:       MPIUNI_ARG(info),\
    936:       MPI_Comm_dup(comm,newcomm))
    937: #define MPI_Comm_test_inter(comm,flag) (*(flag)=1, MPI_SUCCESS)
    938: #define MPI_Comm_remote_size(comm,size) (*(size)=1 ,MPI_SUCCESS)
    939: #define MPI_Comm_remote_group(comm,group) MPI_SUCCESS
    940: #define MPI_Intercomm_create(local_comm,local_leader,peer_comm,\
    941:      remote_leader,tag,newintercomm) MPI_SUCCESS
    942: #define MPI_Intercomm_merge(intercomm,high,newintracomm) MPI_SUCCESS
    943: #define MPI_Topo_test(comm,flag) MPI_SUCCESS
    944: #define MPI_Cart_create(comm_old,ndims,dims,periods,\
    945:      reorder,comm_cart) MPIUni_Abort(MPI_COMM_WORLD,0)
    946: #define MPI_Dims_create(nnodes,ndims,dims) MPIUni_Abort(MPI_COMM_WORLD,0)
    947: #define MPI_Graph_create(comm,a,b,c,d,e) MPIUni_Abort(MPI_COMM_WORLD,0)
    948: #define MPI_Graphdims_Get(comm,nnodes,nedges) MPIUni_Abort(MPI_COMM_WORLD,0)
    949: #define MPI_Graph_get(comm,a,b,c,d) MPIUni_Abort(MPI_COMM_WORLD,0)
    950: #define MPI_Cartdim_get(comm,ndims) MPIUni_Abort(MPI_COMM_WORLD,0)
    951: #define MPI_Cart_get(comm,maxdims,dims,periods,coords) \
    952:      MPIUni_Abort(MPI_COMM_WORLD,0)
    953: #define MPI_Cart_rank(comm,coords,rank) MPIUni_Abort(MPI_COMM_WORLD,0)
    954: #define MPI_Cart_coords(comm,rank,maxdims,coords) \
    955:      MPIUni_Abort(MPI_COMM_WORLD,0)
    956: #define MPI_Graph_neighbors_count(comm,rank,nneighbors) \
    957:      MPIUni_Abort(MPI_COMM_WORLD,0)
    958: #define MPI_Graph_neighbors(comm,rank,maxneighbors,neighbors) \
    959:      MPIUni_Abort(MPI_COMM_WORLD,0)
    960: #define MPI_Cart_shift(comm,direction,disp,rank_source,rank_dest) \
    961:      MPIUni_Abort(MPI_COMM_WORLD,0)
    962: #define MPI_Cart_sub(comm,remain_dims,newcomm) MPIUni_Abort(MPI_COMM_WORLD,0)
    963: #define MPI_Cart_map(comm,ndims,dims,periods,newrank) MPIUni_Abort(MPI_COMM_WORLD,0)
    964: #define MPI_Graph_map(comm,a,b,c,d) MPIUni_Abort(MPI_COMM_WORLD,0)

    966: #define MPI_Get_processor_name(name,result_len)                         \
    967:      (*(result_len) = 9,MPIUNI_Memcpy(name,"localhost",10*MPI_sizeof(MPI_CHAR)))
    968: #define MPI_Errhandler_create(function,errhandler) \
    969:      (MPIUNI_ARG(function),\
    970:       *(errhandler) = MPI_ERRORS_RETURN,\
    971:       MPI_SUCCESS)
    972: #define MPI_Errhandler_set(comm,errhandler) \
    973:      (MPIUNI_ARG(comm),\
    974:       MPIUNI_ARG(errhandler),\
    975:       MPI_SUCCESS)
    976: #define MPI_Errhandler_get(comm,errhandler) \
    977:      (MPIUNI_ARG(comm),\
    978:       (*errhandler) = MPI_ERRORS_RETURN,\
    979:       MPI_SUCCESS)
    980: #define MPI_Errhandler_free(errhandler) \
    981:      (*(errhandler) = MPI_ERRHANDLER_NULL,\
    982:       MPI_SUCCESS)
    983: #define MPI_Error_string(errorcode,string,result_len) \
    984:      (MPIUNI_ARG(errorcode),\
    985:      (errorcode == MPI_ERR_NOSUPPORT) ? \
    986:        (*(result_len) = 35, MPIUNI_Memcpy(string,"MPI error, not supported by MPI-uni",35*MPI_sizeof(MPI_CHAR))) : \
    987:       (*(result_len) = 9, MPIUNI_Memcpy(string,"MPI error",9*MPI_sizeof(MPI_CHAR))))
    988: #define MPI_Error_class(errorcode,errorclass) \
    989:      (*(errorclass) = errorcode, MPI_SUCCESS)
    990: #define MPI_Wtick() 1.0
    991: #define MPI_Pcontrol(level) MPI_SUCCESS

    993: /* MPI-IO additions */

    995: typedef int MPI_File;
    996: #define MPI_FILE_NULL 0

    998: typedef int MPI_Offset;

    1000: #define MPI_MODE_RDONLY  0
    1001: #define MPI_MODE_WRONLY  0
    1002: #define MPI_MODE_CREATE  0

    1004: #define MPI_File_open(comm,filename,amode,info,mpi_fh) \
    1005:   (MPIUNI_ARG(comm),\
    1006:    MPIUNI_ARG(filename),\
    1007:    MPIUNI_ARG(amode),\
    1008:    MPIUNI_ARG(info),\
    1009:    MPIUNI_ARG(mpi_fh),\
    1010:    MPIUni_Abort(MPI_COMM_WORLD,0))

    1012: #define MPI_File_close(mpi_fh) \
    1013:   (MPIUNI_ARG(mpi_fh),\
    1014:    MPIUni_Abort(MPI_COMM_WORLD,0))

    1016: #define MPI_File_set_view(mpi_fh,disp,etype,filetype,datarep,info) \
    1017:   (MPIUNI_ARG(mpi_fh),\
    1018:    MPIUNI_ARG(disp),\
    1019:    MPIUNI_ARG(etype),\
    1020:    MPIUNI_ARG(filetype),\
    1021:    MPIUNI_ARG(datarep),\
    1022:    MPIUNI_ARG(info),\
    1023:    MPIUni_Abort(MPI_COMM_WORLD,0))

    1025: #define MPI_File_write_all(mpi_fh,buf,count,datatype,status) \
    1026:   (MPIUNI_ARG(mpi_fh),\
    1027:    MPIUNI_ARG(buf),\
    1028:    MPIUNI_ARG(count),\
    1029:    MPIUNI_ARG(datatype),\
    1030:    MPIUNI_ARG(status),\
    1031:    MPIUni_Abort(MPI_COMM_WORLD,0))

    1033: #define MPI_File_read_all(mpi_fh,buf,count,datatype,status) \
    1034:   (MPIUNI_ARG(mpi_fh),\
    1035:    MPIUNI_ARG(buf),\
    1036:    MPIUNI_ARG(count),\
    1037:    MPIUNI_ARG(datatype),\
    1038:    MPIUNI_ARG(status),\
    1039:    MPIUni_Abort(MPI_COMM_WORLD,0))

    1041: #define MPI_File_write_at(mpi_fh,off,buf,count,datatype,status) \
    1042:   (MPIUNI_ARG(mpi_fh),\
    1043:    MPIUNI_ARG(off),\
    1044:    MPIUNI_ARG(buf),\
    1045:    MPIUNI_ARG(count),\
    1046:    MPIUNI_ARG(datatype),\
    1047:    MPIUNI_ARG(status),\
    1048:    MPIUni_Abort(MPI_COMM_WORLD,0))

    1050: #define MPI_File_read_at(mpi_fh,off,buf,count,datatype,status) \
    1051:   (MPIUNI_ARG(mpi_fh),\
    1052:    MPIUNI_ARG(off),\
    1053:    MPIUNI_ARG(buf),\
    1054:    MPIUNI_ARG(count),\
    1055:    MPIUNI_ARG(datatype),\
    1056:    MPIUNI_ARG(status),\
    1057:    MPIUni_Abort(MPI_COMM_WORLD,0))

    1059: #define MPI_File_write_at_all(mpi_fh,off,buf,count,datatype,status) \
    1060:   (MPIUNI_ARG(mpi_fh),\
    1061:    MPIUNI_ARG(off),\
    1062:    MPIUNI_ARG(buf),\
    1063:    MPIUNI_ARG(count),\
    1064:    MPIUNI_ARG(datatype),\
    1065:    MPIUNI_ARG(status),\
    1066:    MPIUni_Abort(MPI_COMM_WORLD,0))

    1068: #define MPI_File_read_at_all(mpi_fh,off,buf,count,datatype,status) \
    1069:   (MPIUNI_ARG(mpi_fh),\
    1070:    MPIUNI_ARG(off),\
    1071:    MPIUNI_ARG(buf),\
    1072:    MPIUNI_ARG(count),\
    1073:    MPIUNI_ARG(datatype),\
    1074:    MPIUNI_ARG(status),\
    1075:    MPIUni_Abort(MPI_COMM_WORLD,0))

    1077:   /* called from PetscInitialize() - so return success */
    1078: #define MPI_Register_datarep(name,read_conv_fn,write_conv_fn,extent_fn,state) \
    1079:   (MPIUNI_ARG(name),\
    1080:    MPIUNI_ARG(read_conv_fn),\
    1081:    MPIUNI_ARG(write_conv_fn),\
    1082:    MPIUNI_ARG(extent_fn),\
    1083:    MPIUNI_ARG(state),\
    1084:    MPI_SUCCESS)

    1086: #define MPI_Type_create_subarray(ndims,array_of_sizes,array_of_subsizes,array_of_starts,order,oldtype,newtype) \
    1087:   (MPIUNI_ARG(ndims),\
    1088:    MPIUNI_ARG(array_of_sizes),\
    1089:    MPIUNI_ARG(array_of_subsizes),\
    1090:    MPIUNI_ARG(array_of_starts),\
    1091:    MPIUNI_ARG(order),\
    1092:    MPIUNI_ARG(oldtype),\
    1093:    MPIUNI_ARG(newtype),\
    1094:    MPIUni_Abort(MPI_COMM_WORLD,0))

    1096: #define MPI_Type_create_resized(oldtype,lb,extent,newtype) \
    1097:   (MPIUNI_ARG(oldtype),\
    1098:    MPIUNI_ARG(lb),\
    1099:    MPIUNI_ARG(extent),\
    1100:    MPIUNI_ARG(newtype),\
    1101:    MPIUni_Abort(MPI_COMM_WORLD,0))

    1103: #define MPI_Type_create_indexed_block(count,blocklength,array_of_displacements,oldtype,newtype) \
    1104:   (MPIUNI_ARG(count),\
    1105:    MPIUNI_ARG(blocklength),\
    1106:    MPIUNI_ARG(array_of_displacements),\
    1107:    MPIUNI_ARG(oldtype),\
    1108:    MPIUNI_ARG(newtype),\
    1109:    MPIUni_Abort(MPI_COMM_WORLD,0))

    1111: #if defined(__cplusplus)
    1112: }
    1113: #endif
    1114: #endif