ParExchangeMesh.hpp

Go to the documentation of this file.

// It exchanges entities and sets
// calls moab::ParallelComm::exchange_owned_meshs
// 

#ifndef PAR_EXCHANGE_MESH_HPP
#define PAR_EXCHANGE_MESH_HPP

#include "meshkit/MeshScheme.hpp"
#include "meshkit/ModelEnt.hpp"
#include "moab/ParallelComm.hpp"

namespace MeshKit
{
using namespace moab;

class ParExchangeMesh : public MeshScheme
{
public:

  ParExchangeMesh(MKCore *mkcore, const MEntVector &me_vec);

  virtual ~ParExchangeMesh();

  static const char* name() 
    { return "ParExchangeMesh"; }

  static bool can_mesh(iBase_EntityType dim)
    { return iBase_ALL_TYPES == dim; }
   
  static bool can_mesh(ModelEnt *me)
    { return (me->dimension() == 4); }
  
  static const moab::EntityType* output_types();

  virtual const moab::EntityType* mesh_types_arr() const
    { return output_types(); }

  virtual void setup_this();

  virtual void execute_this();

  // debug function
  void print_mesh();

private:

  ParallelComm* m_mpcomm; // mesh parallel communication

  unsigned int m_rank, m_proc_size;

  iBase_TagHandle m_mPuniqueIDTag;

  std::vector< Range* > m_shared_entities;

  std::vector< unsigned int > m_shared_procs;

  std::vector< MPI_Request > m_recv_reqs, m_recv_remoteh_reqs;

  int get_shared_list(const int to_proc);
};
}

#endif