MPISum.c

Go to the documentation of this file.

#ifndef serial
#include <mpi.h>
#endif

int MPISum_integer(
                    int   *global, 
                    int   *var,    
                    int   size,    
                    void  *comm    
                  )
{
#ifdef serial
  int i;
  for (i = 0; i < size; i++)  global[i] = var[i];
#else
  MPI_Allreduce((var==global?MPI_IN_PLACE:var),global,size,MPI_INT,MPI_SUM,*((MPI_Comm*)comm));
#endif
  return(0);
}

int MPISum_double(
                    double  *global, 
                    double  *var,    
                    int     size,    
                    void    *comm    
                 )
{
#ifdef serial
  int i;
  for (i = 0; i < size; i++)  global[i] = var[i];
#else
  MPI_Allreduce((var==global?MPI_IN_PLACE:var),global,size,MPI_DOUBLE,MPI_SUM,*((MPI_Comm*)comm));
#endif
  return(0);
}