collectivecommunication.hh

Go to the documentation of this file.

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#ifndef DUNE_COLLECTIVECOMMUNICATION_HH
#define DUNE_COLLECTIVECOMMUNICATION_HH
#include <iostream>
#include <complex>
#include <algorithm>
 
#include <dune/common/binaryfunctions.hh>
#include <dune/common/exceptions.hh>
#include <dune/common/unused.hh>
 
namespace Dune
{
 
  /* define some type that definitely differs from MPI_Comm */
  struct No_Comm {};
 
 
  template<typename Communicator>
  class CollectiveCommunication
  {
  public:
    CollectiveCommunication()
    {}
 
    CollectiveCommunication (const Communicator&)
    {}
 
    int rank () const
    {
      return 0;
    }
 
    int size () const
    {
      return 1;
    }
 
    template<typename T>
    T sum (const T& in) const
    {
      return in;
    }
 
    template<typename T>
    int sum (T* inout, int len) const
    {
      DUNE_UNUSED_PARAMETER(inout);
      DUNE_UNUSED_PARAMETER(len);
      return 0;
    }
 
    template<typename T>
    T prod (const T& in) const
    {
      return in;
    }
 
    template<typename T>
    int prod (T* inout, int len) const
    {
      DUNE_UNUSED_PARAMETER(inout);
      DUNE_UNUSED_PARAMETER(len);
      return 0;
    }
 
    template<typename T>
    T min (const T& in) const
    {
      return in;
    }
 
    template<typename T>
    int min (T* inout, int len) const
    {
      DUNE_UNUSED_PARAMETER(inout);
      DUNE_UNUSED_PARAMETER(len);
      return 0;
    }
 
    template<typename T>
    T max (const T& in) const
    {
      return in;
    }
 
    template<typename T>
    int max (T* inout, int len) const
    {
      DUNE_UNUSED_PARAMETER(inout);
      DUNE_UNUSED_PARAMETER(len);
      return 0;
    }
 
    int barrier () const
    {
      return 0;
    }
 
    template<typename T>
    int broadcast (T* inout, int len, int root) const
    {
      DUNE_UNUSED_PARAMETER(inout);
      DUNE_UNUSED_PARAMETER(len);
      DUNE_UNUSED_PARAMETER(root);
      return 0;
    }
 
    template<typename T>
    int gather (const T* in, T* out, int len, int root) const     // note out must have same size as in
    {
      DUNE_UNUSED_PARAMETER(root);
      for (int i=0; i<len; i++)
        out[i] = in[i];
      return 0;
    }
 
    template<typename T>
    int gatherv (const T* in, int sendlen, T* out, int* recvlen, int* displ, int root) const
    {
      DUNE_UNUSED_PARAMETER(recvlen);
      DUNE_UNUSED_PARAMETER(root);
      for (int i=*displ; i<sendlen; i++)
        out[i] = in[i];
      return 0;
    }
 
    template<typename T>
    int scatter (const T* send, T* recv, int len, int root) const // note out must have same size as in
    {
      DUNE_UNUSED_PARAMETER(root);
      for (int i=0; i<len; i++)
        recv[i] = send[i];
      return 0;
    }
 
    template<typename T>
    int scatterv (const T* send, int* sendlen, int* displ, T* recv, int recvlen, int root) const
    {
      DUNE_UNUSED_PARAMETER(recvlen);
      DUNE_UNUSED_PARAMETER(root);
      for (int i=*displ; i<*sendlen; i++)
        recv[i] = send[i];
      return 0;
    }
 
    template<typename T>
    int allgather(const T* sbuf, int count, T* rbuf) const
    {
      for(const T* end=sbuf+count; sbuf < end; ++sbuf, ++rbuf)
        *rbuf=*sbuf;
      return 0;
    }
 
    template<typename T>
    int allgatherv (const T* in, int sendlen, T* out, int* recvlen, int* displ) const
    {
      DUNE_UNUSED_PARAMETER(recvlen);
      for (int i=*displ; i<sendlen; i++)
        out[i] = in[i];
      return 0;
    }
 
    template<typename BinaryFunction, typename Type>
    int allreduce(Type* inout, int len) const
    {
      DUNE_UNUSED_PARAMETER(inout);
      DUNE_UNUSED_PARAMETER(len);
      return 0;
    }
 
    template<typename BinaryFunction, typename Type>
    void allreduce(const Type* in, Type* out, int len) const
    {
      std::copy(in, in+len, out);
      return;
    }
 
  };
}
 
#endif