defaults.hh

Go to the documentation of this file.

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#ifndef DUNE_MMESH_CGAL_DEFAULTS_HH
#define DUNE_MMESH_CGAL_DEFAULTS_HH
 
// CGAL includes
#include "includes.hh"
 
#define CGAL_NO_POSTCONDITIONS
 
namespace Dune
{
  namespace MMeshDefaults
  {
 
    template< int dim >
    struct ElementInfo {
      std::size_t insertionIndex;
      std::size_t index;
      std::array< std::size_t, dim+1 > cgalIndex;
      size_t domainMarker = 0;
      int mark = 0;
      bool isNew = false;
      bool mightVanish = false;
      std::size_t componentNumber = 0;
      int rank = 0;
      int partition = 0;
      std::unordered_set<int> connectivity;
    };
 
    struct VertexInfo {
      std::size_t id;
      bool idWasSet = false;
      std::size_t index;
      std::size_t insertionLevel = 0;
      bool isInterface = false;
      int boundaryFlag = -1;
      int partition = 0;
    };
 
    template<>
    class Triangulation<2>
    {
    private:
      typedef CGAL::Exact_predicates_inexact_constructions_kernel K;
 
      typedef CGAL::Triangulation_vertex_base_2<K> Vbbb;
      typedef CGAL::Triangulation_vertex_base_with_info_2<VertexInfo, K, Vbbb> Vbb;
      typedef CGAL::Triangulation_hierarchy_vertex_base_2<Vbb> Vb;
 
      typedef CGAL::Triangulation_face_base_2<K> Fbbb;
      typedef CGAL::Triangulation_face_base_with_info_2<ElementInfo<2>, K, Fbbb> Fbb;
      typedef CGAL::Triangulation_face_base_2<K,Fbb> Fb;
 
      typedef CGAL::Triangulation_data_structure_2<Vb, Fb> Tds;
 
    public:
      typedef CGAL::Triangulation_2<K, Tds> type;
    };
 
 
    template<>
    class Triangulation<3>
    {
    private:
      typedef CGAL::Exact_predicates_inexact_constructions_kernel K;
 
      typedef CGAL::Triangulation_vertex_base_with_info_3<VertexInfo, K> Vb;
 
      typedef CGAL::Triangulation_cell_base_3<K> Fbb;
      typedef CGAL::Triangulation_cell_base_with_info_3<ElementInfo<3>, K, Fbb> Fb;
 
      typedef CGAL::Triangulation_data_structure_3<Vb, Fb> Tds;
 
    public:
      typedef CGAL::Triangulation_3<K, Tds> type;
    };
 
 
    template<>
    class Delaunay<2>
    {
    private:
      typedef CGAL::Exact_predicates_inexact_constructions_kernel K;
 
      typedef CGAL::Delaunay_mesh_vertex_base_2<K> Vbbb;
      typedef CGAL::Triangulation_vertex_base_with_info_2<VertexInfo, K, Vbbb> Vbb;
      typedef CGAL::Triangulation_hierarchy_vertex_base_2<Vbb> Vb;
 
      typedef CGAL::Delaunay_mesh_face_base_2<K> Fbbb;
      typedef CGAL::Triangulation_face_base_with_info_2<ElementInfo<2>, K, Fbbb> Fbb;
      typedef CGAL::Triangulation_face_base_2<K,Fbb> Fb;
 
      typedef CGAL::Triangulation_data_structure_2<Vb, Fb> Tds;
 
    public:
      typedef CGAL::Delaunay_triangulation_2<K, Tds> type;
    };
 
 
    template<>
    class Delaunay<3>
    {
    private:
      typedef CGAL::Exact_predicates_inexact_constructions_kernel K;
 
      typedef CGAL::Triangulation_vertex_base_with_info_3<VertexInfo, K> Vb;
 
      typedef CGAL::Delaunay_triangulation_cell_base_3<K> Fbb;
      typedef CGAL::Triangulation_cell_base_with_info_3<ElementInfo<3>, K, Fbb> Fb;
 
      typedef CGAL::Triangulation_data_structure_3<Vb, Fb> Tds;
 
    public:
      typedef CGAL::Delaunay_triangulation_3<K, Tds> type;
    };
 
 
  } // end namespace MMeshDefaults
 
} // end namespace Dune
 
#endif