netmpi.h

Go to the documentation of this file.

#ifndef  _MY_NetMPI_
#define  _MY_NetMPI_
#include "globals.h"
#include "vecdouble.h"
/* #if OMPI_MINOR_VERSION < 4 */
/* #ifndef GPU_SRV //the gpu srv has a newer version of open mpi and doesnot need this crap. */
/* # ifdef SEEK_SET */
/* # define MB_SEEK_SET SEEK_SET */
/* # define MB_SEEK_CUR SEEK_CUR */
/* # define MB_SEEK_END SEEK_END */
/* # undef SEEK_SET */
/* # undef SEEK_CUR */
/* # undef SEEK_END */
/* # endif */
/* #include <mpi.h> */
/* # ifdef MB_SEEK_SET */
/* # define SEEK_SET MB_SEEK_SET */
/* # define SEEK_CUR MB_SEEK_CUR */
/* # define SEEK_END MB_SEEK_END */
/* # undef MB_SEEK_SET */
/* # undef MB_SEEK_CUR */
/* # undef MB_SEEK_END */
/* #endif */
/* #endif */
/* #endif */

/* Ok we are having a lot of warnings about not used parameters when we use
mpi with C++. This is annoying. I noted that this problem happens if we use
openmpi or mpich. To avoid this problem, we use a pragma asking the compiler
will ignore the warnings. If you are trying to compile this code in a compiler other
than gcc, please comment the pragma lines above and get used with the warning messages that will pop up. The MPI developers should have included this lines in their code.
This can be considered a bug of the new MPI implementation. For me it is. 
*/
//#pragma GCC diagnostic push
#pragma GCC diagnostic ignored "-Wunused-parameter"
#include <mpi.h>
#pragma GCC diagnostic warning "-Wunused-parameter"
#include <iostream>

typedef std::vector<MPI_Request> VecReq;
class NetMPI 
{
 private:
  static bool m_loop;
  static bool b1,b2,m_print;
 protected:
  static int m_nP,m_rank;
 public:
  static int rank();
  static int nProcess(); 
  static bool isLastProcess(); 
  static std::ostream *m_out;
  static void setTrace(bool b){m_print=b;}
  static std::ostream& getLog();
  static void setLog(std::ostream &out){m_out=&out;}

  static void Isend(VecDouble &v,int dst,int tag,MPI_Request &req)
  {
    if (dst >= 0 && dst < nProcess())
      MPI_Isend(&(v(0)),v.size(),MPI_DOUBLE,dst,tag,MPI_COMM_WORLD,&req);
    
  } 
  static void Barrier()
  {
    MPI_Barrier(MPI_COMM_WORLD);
 }
  static void Irecv(VecDouble &v,int src,int tag,MPI_Request &req)
  {
    if (src >= 0 && src < nProcess())
      MPI_Irecv(&(v(0)),v.size(),MPI_DOUBLE,src,tag,MPI_COMM_WORLD,&req);
  } 
  static void exchangeDataRedBlack(VecIndex &vLSndMap,VecDouble &vLSnd,VecIndex &vLRcvMap,VecDouble &vLRcv,VecIndex &vRSndMap,VecDouble &vRSnd,VecIndex &vRRcvMap,VecDouble &vRRcv,VecDouble &data,int tag); 
  static void exchangeDataRedBlack(VecIndex &vLSndMap,VecDouble &vLSnd,VecDouble &vLRcv,VecIndex &vRSndMap,VecDouble &vRSnd,VecDouble &vRRcv,VecDouble &data,int tag,VecReq &req); 
  static void exchangeDataRedBlack(VecIndex &vLSndMap,VecDouble &vLSnd,VecIndex &vLRcvMap,VecDouble &vLRcv,VecIndex &vRSndMap,VecDouble &vRSnd,VecIndex &vRRcvMap,VecDouble &vRRcv,VecDouble &dataSnd, VecDouble &dataRcv, int  tag); 

  static void trace(const char *format,...); 
  static void debugPoint();
  static double getMinValue(double localValue);
};




#endif

---