Apply_move_test.cpp

Go to the documentation of this file.

/*******************************************************************************
 Causal Dynamical Triangulations in C++ using CGAL
 
 Copyright © 2019 Adam Getchell
 ******************************************************************************/
 
 
#include "Apply_move.hpp"
 
#include <doctest/doctest.h>
 
#include "Ergodic_moves_3.hpp"
 
using namespace std;
 
SCENARIO("Apply an ergodic move to 2+1 manifolds" *
         doctest::test_suite("apply"))
{
  GIVEN("A 2+1 dimensional spherical manifold.")
  {
    auto constexpr desired_simplices  = 9600;
    auto constexpr desired_timeslices = 7;
    manifolds::Manifold_3 manifold(desired_simplices, desired_timeslices);
    REQUIRE(manifold.is_correct());
    // Copy of manifold
    auto manifold_before = manifold;
    WHEN("A null move is applied to the manifold.")
    {
      spdlog::debug("Applying null move to manifold.\n");
      if (auto result = apply_move(manifold, ergodic_moves::null_move); result)
      {
        manifold = result.value();
        // Update Geometry and Foliated_triangulation with new info
        manifold.update();
      }
      else
      {
        spdlog::debug("{}", result.error());
        REQUIRE(result.has_value());
      }
      THEN("The resulting manifold is valid and unchanged.")
      {
        CHECK(manifold.is_valid());
        CHECK_EQ(manifold_before.simplices(), manifold.simplices());
        CHECK_EQ(manifold_before.faces(), manifold.faces());
        CHECK_EQ(manifold_before.edges(), manifold.edges());
        CHECK_EQ(manifold_before.vertices(), manifold.vertices());
        // Human verification
        fmt::print("Old manifold.\n");
        manifold_before.print_details();
        fmt::print("New manifold after null move:\n");
        manifold.print_details();
      }
    }
    WHEN("A (2,3) move is applied to the manifold.")
    {
      spdlog::debug("Applying (2,3) move to manifold.\n");
      if (auto result = apply_move(manifold, ergodic_moves::do_23_move); result)
      {
        manifold = result.value();
        // Update Geometry and Foliated_triangulation with new info
        manifold.update();
      }
      else
      {
        spdlog::debug("{}", result.error());
        REQUIRE(result.has_value());
      }
      THEN("The resulting manifold has the applied move.")
      {
        CHECK(ergodic_moves::check_move(manifold_before, manifold,
                                        move_tracker::move_type::TWO_THREE));
        // Human verification
        fmt::print("Old manifold.\n");
        manifold_before.print_details();
        fmt::print("New manifold after (2,3) move:\n");
        manifold.print_details();
      }
    }
    WHEN("A (3,2) move is applied to the manifold.")
    {
      spdlog::debug("Applying (3,2) move to manifold.\n");
      if (auto result = apply_move(manifold, ergodic_moves::do_32_move); result)
      {
        manifold = result.value();
        // Update Geometry and Foliated_triangulation with new info
        manifold.update();
      }
      else
      {
        spdlog::debug("{}", result.error());
        // Stop further tests
        REQUIRE(result.has_value());
      }
      THEN("The resulting manifold has the applied move.")
      {
        CHECK(ergodic_moves::check_move(manifold_before, manifold,
                                        move_tracker::move_type::THREE_TWO));
        // Human verification
        fmt::print("Old manifold.\n");
        manifold_before.print_details();
        fmt::print("New manifold after (3,2) move:\n");
        manifold.print_details();
      }
    }
    WHEN("A (2,6) move is applied to the manifold.")
    {
      spdlog::debug("Applying (2,6) move to manifold.\n");
      if (auto result = apply_move(manifold, ergodic_moves::do_26_move); result)
      {
        manifold = result.value();
        // Update Geometry and Foliated_triangulation with new info
        manifold.update();
      }
      else
      {
        spdlog::debug("{}", result.error());
        // Stop further tests
        REQUIRE(result.has_value());
      }
      THEN("The resulting manifold has the applied move.")
      {
        CHECK(ergodic_moves::check_move(manifold_before, manifold,
                                        move_tracker::move_type::TWO_SIX));
        // Human verification
        fmt::print("Old manifold.\n");
        manifold_before.print_details();
        fmt::print("New manifold after (2,6) move:\n");
        manifold.print_details();
      }
    }
    WHEN("A (6,2) move is applied to the manifold.")
    {
      spdlog::debug("Applying (6,2) move to manifold.\n");
      auto result = apply_move(manifold, ergodic_moves::do_62_move);
      if (result)
      {
        manifold = result.value();
        // Update Geometry and Foliated_triangulation with new info
        manifold.update();
        THEN("The resulting manifold has the applied move.")
        {
          CHECK(ergodic_moves::check_move(manifold_before, manifold,
                                          move_tracker::move_type::SIX_TWO));
          // Human verification
          fmt::print("Old manifold.\n");
          manifold_before.print_details();
          fmt::print("New manifold after (6,2) move:\n");
          manifold.print_details();
        }
      }
      else
      {
        spdlog::warn("Cannot apply (6,2) move: {}", result.error());
        doctest::skip("No valid (6,2) move exists in this manifold configuration.");
      }
    }
    WHEN("A (4,4) move is applied to the manifold.")
    {
      spdlog::debug("Applying (4,4) move to manifold.\n");
      if (auto result = apply_move(manifold, ergodic_moves::do_44_move); result)
      {
        manifold = result.value();
        // Update Geometry and Foliated_triangulation with new info
        manifold.update();
      }
      else
      {
        spdlog::debug("{}", result.error());
        // Stop further tests
        REQUIRE(result.has_value());
      }
      THEN("The resulting manifold has the applied move.")
      {
        CHECK(ergodic_moves::check_move(manifold_before, manifold,
                                        move_tracker::move_type::FOUR_FOUR));
        // Human verification
        fmt::print("Old manifold.\n");
        manifold_before.print_details();
        fmt::print("New manifold after (4,4) move:\n");
        manifold.print_details();
      }
    }
  }
}