/*$Id: mpi.c,v 1.1.6.1 2002/04/24 03:25:59 erik Exp $*//*      This provides a few of the MPI-uni functions that cannot be implemented    with C macros*/#include "mpi.h"#if defined (MPIUNI_USE_STDCALL)#define MPIUNI_STDCALL __stdcall#else#define MPIUNI_STDCALL#endif/*     The system call exit() is not properly prototyped on all systems   hence we fake it by creating our own prototype*//* This file does not include petsc.h hence cannot use EXTERN_C_XXXX*/#if defined(__cplusplus)extern "C" {  void exit(int);}#else  void exit(int);#endif#define MPI_SUCCESS 0#define MPI_FAILURE 1void    *MPIUNI_TMP        = 0;int     MPIUNI_DATASIZE[5] = { sizeof(int),sizeof(float),sizeof(double),2*sizeof(double),sizeof(char)};/*       With MPI Uni there is only one communicator, which is called 1.*/#define MAX_ATTR 128typedef struct {  void                *extra_state;  void                *attribute_val;  int                 active;  MPI_Delete_function *del;} MPI_Attr;static MPI_Attr attr[MAX_ATTR];static int      num_attr = 1,mpi_tag_ub = 100000000;/*    To avoid problems with prototypes to the system memcpy() it is duplicated here*/int MPIUNI_Memcpy(void *a,void* b,int n) {  int  i;  char *aa= (char*)a;  char *bb= (char*)b;  for (i=0; i<n; i++) aa[i] = bb[i];  return 0;}/*   Used to set the built-in MPI_TAG_UB attribute*/static int Keyval_setup(void){  attr[0].active        = 1;  attr[0].attribute_val = &mpi_tag_ub;  return 0;}/*         These functions are mapped to the Petsc_ name by ./mpi.h*/int Petsc_MPI_Keyval_create(MPI_Copy_function *copy_fn,MPI_Delete_function *delete_fn,int *keyval,void *extra_state){  if (num_attr >= MAX_ATTR) MPI_Abort(MPI_COMM_WORLD,1);  attr[num_attr].extra_state = extra_state;  attr[num_attr].del         = delete_fn;  *keyval                    = num_attr++;  return 0;}int Petsc_MPI_Keyval_free(int *keyval){  attr[*keyval].active = 0;  return MPI_SUCCESS;}int Petsc_MPI_Attr_put(MPI_Comm comm,int keyval,void *attribute_val){  attr[keyval].active        = 1;  attr[keyval].attribute_val = attribute_val;  return MPI_SUCCESS;}  int Petsc_MPI_Attr_delete(MPI_Comm comm,int keyval){  if (attr[keyval].active && attr[keyval].del) {    (*(attr[keyval].del))(comm,keyval,attr[keyval].attribute_val,attr[keyval].extra_state);  }  attr[keyval].active        = 0;  attr[keyval].attribute_val = 0;  return MPI_SUCCESS;}int Petsc_MPI_Attr_get(MPI_Comm comm,int keyval,void *attribute_val,int *flag){  if (!keyval) Keyval_setup();  *flag                  = attr[keyval].active;  *(int **)attribute_val = (int *)attr[keyval].attribute_val;  return MPI_SUCCESS;}static int dups = 0;int Petsc_MPI_Comm_dup(MPI_Comm comm,MPI_Comm *out){  *out = comm;  dups++;  return 0;}int Petsc_MPI_Comm_free(MPI_Comm *comm){  int i;  if (--dups) return MPI_SUCCESS;  for (i=0; i<num_attr; i++) {    if (attr[i].active && attr[i].del) {      (*attr[i].del)(*comm,i,attr[i].attribute_val,attr[i].extra_state);    }    attr[i].active = 0;  }  return MPI_SUCCESS;}/* --------------------------------------------------------------------------*/int Petsc_MPI_Abort(MPI_Comm comm,int errorcode) {  exit(errorcode);   return MPI_SUCCESS;}static int MPI_was_initialized = 0;int Petsc_MPI_Initialized(int *flag){  *flag = MPI_was_initialized;  return 0;}int Petsc_MPI_Finalize(void){  MPI_was_initialized = 0;  return 0;}/* -------------------     Fortran versions of several routines ------------------ */#if defined(__cplusplus)extern "C" {#endif/******mpi_init*******/void MPIUNI_STDCALL  mpi_init(int *ierr){  MPI_was_initialized = 1;  *ierr = MPI_SUCCESS;}void MPIUNI_STDCALL  mpi_init_(int *ierr){  MPI_was_initialized = 1;  *ierr = MPI_SUCCESS;}void MPIUNI_STDCALL  mpi_init__(int *ierr){  MPI_was_initialized = 1;  *ierr = MPI_SUCCESS;}void MPIUNI_STDCALL  MPI_INIT(int *ierr){  MPI_was_initialized = 1;  *ierr = MPI_SUCCESS;}/******mpi_finalize*******/void MPIUNI_STDCALL  mpi_finalize(int *ierr){  *ierr = MPI_SUCCESS;}void MPIUNI_STDCALL  mpi_finalize_(int *ierr){  *ierr = MPI_SUCCESS;}void MPIUNI_STDCALL  mpi_finalize__(int *ierr){  *ierr = MPI_SUCCESS;}void MPIUNI_STDCALL  MPI_FINALIZE(int *ierr){  *ierr = MPI_SUCCESS;}/******mpi_comm_size*******/void MPIUNI_STDCALL mpi_comm_size(MPI_Comm *comm,int *size,int *ierr) {  *size = 1;  *ierr = 0;}void MPIUNI_STDCALL mpi_comm_size_(MPI_Comm *comm,int *size,int *ierr) {  *size = 1;  *ierr = 0;}void MPIUNI_STDCALL mpi_comm_size__(MPI_Comm *comm,int *size,int *ierr) {  *size = 1;  *ierr = 0;}void MPIUNI_STDCALL MPI_COMM_SIZE(MPI_Comm *comm,int *size,int *ierr) {  *size = 1;  *ierr = 0;}/******mpi_comm_rank*******/void MPIUNI_STDCALL mpi_comm_rank(MPI_Comm *comm,int *rank,int *ierr){  *rank=0;  *ierr=MPI_SUCCESS;}void MPIUNI_STDCALL mpi_comm_rank_(MPI_Comm *comm,int *rank,int *ierr){  *rank=0;  *ierr=MPI_SUCCESS;}void MPIUNI_STDCALL mpi_comm_rank__(MPI_Comm *comm,int *rank,int *ierr){  *rank=0;  *ierr=MPI_SUCCESS;}void MPIUNI_STDCALL MPI_COMM_RANK(MPI_Comm *comm,int *rank,int *ierr){  *rank=0;  *ierr=MPI_SUCCESS;}/*******mpi_abort******/void MPIUNI_STDCALL mpi_abort(MPI_Comm *comm,int *errorcode,int *ierr) {  exit(*errorcode);   *ierr = MPI_SUCCESS;}void MPIUNI_STDCALL mpi_abort_(MPI_Comm *comm,int *errorcode,int *ierr) {  exit(*errorcode);  *ierr = MPI_SUCCESS;}void MPIUNI_STDCALL mpi_abort__(MPI_Comm *comm,int *errorcode,int *ierr) {  exit(*errorcode);  *ierr = MPI_SUCCESS;}void MPIUNI_STDCALL MPI_ABORT(MPI_Comm *comm,int *errorcode,int *ierr) {  exit(*errorcode);  *ierr = MPI_SUCCESS;}/*******mpi_allreduce******/void MPIUNI_STDCALL mpi_allreduce(void *sendbuf,void *recvbuf,int *count,int *datatype,int *op,int *comm,int *ierr) {  MPIUNI_Memcpy(recvbuf,sendbuf,(*count)*MPIUNI_DATASIZE[*datatype]);  *ierr = MPI_SUCCESS;} void MPIUNI_STDCALL mpi_allreduce_(void *sendbuf,void *recvbuf,int *count,int *datatype,int *op,int *comm,int *ierr) {  MPIUNI_Memcpy(recvbuf,sendbuf,(*count)*MPIUNI_DATASIZE[*datatype]);  *ierr = MPI_SUCCESS;} void MPIUNI_STDCALL mpi_allreduce__(void *sendbuf,void *recvbuf,int *count,int *datatype,int *op,int *comm,int *ierr) {  MPIUNI_Memcpy(recvbuf,sendbuf,(*count)*MPIUNI_DATASIZE[*datatype]);  *ierr = MPI_SUCCESS;} void MPIUNI_STDCALL MPI_ALLREDUCE(void *sendbuf,void *recvbuf,int *count,int *datatype,int *op,int *comm,int *ierr) {  MPIUNI_Memcpy(recvbuf,sendbuf,(*count)*MPIUNI_DATASIZE[*datatype]);  *ierr = MPI_SUCCESS;} #if defined(__cplusplus)}#endif