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 <type_traits>
 
#include "Foliated_triangulation.hpp"
 
namespace cdt
{
  template <int dimension>
  struct Geometry;

  template <>
  struct [[nodiscard("This contains data!")]] Geometry<3>
  {
    static_assert(std::is_nothrow_swappable_v<Int_precision>,
                  "Geometry swap requires non-throwing scalar swaps.");

    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.number_of_three_one_cells())}
        , N3_13{static_cast<Int_precision>(
              triangulation.number_of_one_three_cells())}
        , N3_31_13{N3_31 + N3_13}
        , N3_22{static_cast<Int_precision>(
              triangulation.number_of_two_two_cells())}
        , 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
    {
      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>;
}  // namespace cdt
 
#endif  // CDT_PLUSPLUS_GEOMETRY_HPP