Triangulation_traits.hpp

Go to the documentation of this file.

/*******************************************************************************
 Causal Dynamical Triangulations in C++ using CGAL
 
 Copyright © 2021 Adam Getchell
 ******************************************************************************/

 
#ifndef CDT_PLUSPLUS_TRIANGULATION_TRAITS_HPP
#define CDT_PLUSPLUS_TRIANGULATION_TRAITS_HPP
 
#include <CGAL/Delaunay_triangulation_3.h>
#include <CGAL/Delaunay_triangulation_cell_base_3.h>
#include <CGAL/Exact_predicates_inexact_constructions_kernel.h>
#include <CGAL/point_generators_3.h>
#include <CGAL/tags.h>
#include <CGAL/Triangulation_cell_base_with_info_3.h>
#include <CGAL/Triangulation_data_structure_3.h>
#include <CGAL/Triangulation_vertex_base_with_info_3.h>
 
#include <utility>
#include <vector>
 
#include "Settings.hpp"
 
namespace cdt::detail
{
  template <int dimension>
  struct TriangulationTraits;
 
  template <>
  struct TriangulationTraits<3>
  {
    using Kernel = CGAL::Exact_predicates_inexact_constructions_kernel;
    using Vertex_base =
        CGAL::Triangulation_vertex_base_with_info_3<Int_precision, Kernel>;
    using Delaunay_cell_base = CGAL::Delaunay_triangulation_cell_base_3<Kernel>;
    using Cell_base =
        CGAL::Triangulation_cell_base_with_info_3<Int_precision, Kernel,
                                                  Delaunay_cell_base>;
#if defined(CDT_ENABLE_PARALLEL_TRIANGULATION) && \
    CDT_ENABLE_PARALLEL_TRIANGULATION
#ifndef CGAL_LINKED_WITH_TBB
#error "Parallel CDT++ triangulations require CGAL::TBB_support."
#endif
    using Concurrency_tag = CGAL::Parallel_tag;
#else
    using Concurrency_tag = CGAL::Sequential_tag;
#endif
    using Tds = CGAL::Triangulation_data_structure_3<Vertex_base, Cell_base,
                                                     Concurrency_tag>;
    using Delaunay         = CGAL::Delaunay_triangulation_3<Kernel, Tds>;
 
    using Cell_handle      = Delaunay::Cell_handle;
    using Facet            = Delaunay::Facet;
    using Edge_handle      = Delaunay::Edge;
    using Vertex_handle    = Delaunay::Vertex_handle;
    using Point            = Delaunay::Point;
    using Causal_vertices  = std::vector<std::pair<Point, Int_precision>>;

    using squared_distance = Kernel::Compute_squared_distance_3;
 
    using Spherical_points_generator = CGAL::Random_points_on_sphere_3<Point>;
 
    static inline Point const ORIGIN_POINT = Point{0, 0, 0};
  };  // TriangulationTraits<3>
}  // namespace cdt::detail
 
#endif  // CDT_PLUSPLUS_TRIANGULATION_TRAITS_HPP