SphereDecomp.hpp

Go to the documentation of this file.

#ifndef SPHERE_DECOMP_HPP
#define SPHERE_DECOMP_HPP

#include "moab/Interface.hpp"

class SphereDecomp 
{
public:
  SphereDecomp(moab::Interface *impl);

  moab::ErrorCode build_sphere_mesh(const char *sphere_radii_tag_name, 
                                moab::EntityHandle *hex_set = NULL);
  
private:

  moab::ErrorCode compute_nodes(const int dim);

  moab::ErrorCode build_hexes(std::vector<moab::EntityHandle> &sphere_hexes,
                          std::vector<moab::EntityHandle> &interstic_hexes);
  
  moab::ErrorCode subdivide_tet(moab::EntityHandle tet, 
                            std::vector<moab::EntityHandle> &sphere_hexes,
                            std::vector<moab::EntityHandle> &interstic_hexes);
  
  moab::ErrorCode retrieve_subdiv_verts(moab::EntityHandle tet, moab::EntityHandle this_ent,
                                    const moab::EntityHandle *tet_conn,
                                    const int dim, moab::EntityHandle *subdiv_verts);
  
  moab::Tag sphereRadiiTag;

  moab::Tag subdivVerticesTag;
  
  moab::Interface *mbImpl;
  
};
#endif