ParallelMesher.hpp

Go to the documentation of this file.

// It does parallel meshing
// Read and distribute geometry in parallel by iGeom
// Inserts shared vertex/edge/surface entity mesh operations and
// exchange operations to graph
// shared entities are meshed and exchanged, then volume mesh is
// performed at last
// 

#ifndef PARALLELMESHER_HPP
#define PARALLELMESHER_HPP

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

namespace MeshKit
{
class ModelEnt;

using namespace moab;

enum PARALLEL_OP_TYPE
{
  MESH_VERTEX = 0,
  MESH_EDGE,
  MESH_INTER_SURF,
  MESH_NINTER_SURF,
  EXCHANGE_VERTEX,
  EXCHANGE_EDGE,
  SEND_POST_SURF_MESH,
  RECV_SURF_MESH,
  MESH_VOLUME
};

class ParallelMesher : public MeshScheme
{
public:

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

  virtual ~ParallelMesher();

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

  static bool can_mesh(iBase_EntityType dim)
    { return iBase_REGION == dim; }
   
  static bool can_mesh(ModelEnt *me)
    { return canmesh_region(me); }
  
  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();

  void set_mesh_size(double size, int interval);

  // debug function
  void print_mesh();

  // debug function
  void print_geom_info(ModelEnt* me, const int dim,
                       const bool local);
  
private:

  MeshOp* get_mesher(PARALLEL_OP_TYPE type);

  void add_parallel_mesh_op(PARALLEL_OP_TYPE type, bool after = false);

  void check_partition(TDParallel* td_par, ModelEnt* me, int dim);

  ParallelComm* m_mpcomm; // mesh parallel communication

  unsigned int m_rank;

  iBase_TagHandle m_mPuniqueIDTag;
  
  std::vector< MEntSet > m_sEntity;
};
}

#endif