MergeMesh.hpp

Go to the documentation of this file.

#ifndef MERGEMESH_HPP
#define MERGEMESH_HPP

#include "moab/Interface.hpp"
#include "moab/Range.hpp"

namespace moab {

class AdaptiveKDTree;
    
class MergeMesh 
{
public:
    /* \brief Constructor
     */
  MergeMesh(moab::Interface *mbImpl, bool printErrorIn = true);
  
    /* \brief Destructor
     */
  virtual ~MergeMesh();

    /* \brief Merge vertices in elements passed in
     */
   moab::ErrorCode merge_entities(moab::EntityHandle *elems,
                      int elems_size,
                      const double merge_tol,
                      const int do_merge = true,
                      const int update_sets = false,
              moab::Tag merge_tag = 0,
              bool do_higher_dim = true);

  moab::ErrorCode merge_entities(moab::Range &elems,
                                 const double merge_tol,
                                 const int do_merge = true,
                                 const int update_sets = false,
                                 moab::Tag merge_tag = 0,
                 bool do_higher_dim = true);
  
  //Identify higher dimension to be merged
  moab::ErrorCode merge_higher_dimensions(moab::Range &elems);

      //- perform the actual merge
  moab::ErrorCode perform_merge(moab::Tag merged_to);
private:
  //iMesh_Instance imeshImpl;

  double mergeTol, mergeTolSq;

  moab::Tag mergeTag;

    //- given a kdtree, set tag on vertices in leaf nodes with vertices
    //- to which they should be merged
  moab::ErrorCode find_merged_to(moab::EntityHandle &tree_root,
                                 moab::AdaptiveKDTree &tree,
                 moab::Tag merged_to);
  
  moab::Interface *mbImpl;

    //- the tag pointing to the entity to which an entity will be merged
  moab::Tag mbMergeTag;

    //- entities which will go away after the merge
  moab::Range deadEnts;

  //Allow a warning to be suppressed when no merging is done
  bool printError;
};

  inline MergeMesh::MergeMesh(Interface *impl, bool printErrorIn) 
    : mbImpl(impl), printError(printErrorIn)
{
}

inline MergeMesh::~MergeMesh() 
{
  if (mbMergeTag) mbImpl->tag_delete(mbMergeTag);
}

}

#endif