Geometry.hpp

Go to the documentation of this file.

/*******************************************************************************
 Causal Dynamical Triangulations in C++ using CGAL
 
 Copyright © 2018 Adam Getchell
 ******************************************************************************/
 
 
#ifndef CDT_PLUSPLUS_GEOMETRY_HPP
#define CDT_PLUSPLUS_GEOMETRY_HPP
 
#include "Foliated_triangulation.hpp"
 
template <int dimension>
struct Geometry;
 
template <>
struct [[nodiscard("This contains data!")]] Geometry<3>
{
  Int_precision N3{0};  // NOLINT
 
  Int_precision N3_31{0};  // NOLINT
 
  Int_precision N3_13{0};  // NOLINT
 
  Int_precision N3_31_13{0};  // NOLINT
 
  Int_precision N3_22{0};  // NOLINT
 
  Int_precision N2{0};  // NOLINT
 
  Int_precision N1{0};  // NOLINT
 
  Int_precision N1_TL{0};  // NOLINT
 
  Int_precision N1_SL{0};  // NOLINT
 
  Int_precision N0{0};  // NOLINT
 
  Geometry() = default;
 
  explicit Geometry(
      foliated_triangulations::FoliatedTriangulation_3 const& triangulation)
 
      : N3{static_cast<Int_precision>(triangulation.number_of_finite_cells())}
      , N3_31{static_cast<Int_precision>(triangulation.get_three_one().size())}
      , N3_13{static_cast<Int_precision>(triangulation.get_one_three().size())}
      , N3_31_13{N3_31 + N3_13}
      , N3_22{static_cast<Int_precision>(triangulation.get_two_two().size())}
      , N2{static_cast<Int_precision>(triangulation.number_of_finite_facets())}
      , N1{static_cast<Int_precision>(triangulation.number_of_finite_edges())}
      , N1_TL{triangulation.N1_TL()}
      , N1_SL{triangulation.N1_SL()}
      , N0{static_cast<Int_precision>(triangulation.number_of_vertices())}
 
  {}
 
  friend void swap(Geometry& swap_from, Geometry& swap_into) noexcept
  {
#ifndef NDEBUG
    spdlog::debug("{} called.\n", __PRETTY_FUNCTION__);
#endif
    using std::swap;
    swap(swap_from.N3, swap_into.N3);
    swap(swap_from.N3_31, swap_into.N3_31);
    swap(swap_from.N3_13, swap_into.N3_13);
    swap(swap_from.N3_31_13, swap_into.N3_31_13);
    swap(swap_from.N3_22, swap_into.N3_22);
    swap(swap_from.N2, swap_into.N2);
    swap(swap_from.N1, swap_into.N1);
    swap(swap_from.N1_TL, swap_into.N1_TL);
    swap(swap_from.N1_SL, swap_into.N1_SL);
    swap(swap_from.N0, swap_into.N0);
  }  // swap
};  // struct Geometry<3>
 
using Geometry_3 = Geometry<3>;
 
template <>
struct [[nodiscard("This contains data!")]] Geometry<4>
{
  Int_precision N4{0};
  Int_precision N3{0};
  Int_precision N2{0};
  Int_precision N1{0};
  Int_precision N0{0};
};  // struct Geometry<4>
 
using Geometry_4 = Geometry<4>;
 
#endif  // CDT_PLUSPLUS_GEOMETRY_HPP