t8_gtest_boundary_extrude.cxx

/*
  This file is part of t8code.
  t8code is a C library to manage a collection (a forest) of multiple
  connected adaptive space-trees of general element classes in parallel.

  Copyright (C) 2015 the developers

  t8code is free software; you can redistribute it and/or modify
  it under the terms of the GNU General Public License as published by
  the Free Software Foundation; either version 2 of the License, or
  (at your option) any later version.

  t8code is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  GNU General Public License for more details.

  You should have received a copy of the GNU General Public License
  along with t8code; if not, write to the Free Software Foundation, Inc.,
  51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
*/

#include <gtest/gtest.h>
#include <t8_eclass.h>
#include <t8_schemes/t8_default/t8_default.hxx>
#include <test/t8_gtest_custom_assertion.hxx>
#include <test/t8_gtest_macros.hxx>
#include "t8_gtest_dfs_base.hxx"
#include <test/t8_gtest_macros.hxx>

struct class_test_boundary_extrude: public TestDFS
{
 private:
  /* For elements that are on the face of the root element, check that creating the boundary element
   * and extruding it results in the original element
    */
  void
  check_element () override
  {
    const int num_faces = scheme->element_get_num_faces (eclass, element);
    for (int iface = 0; iface < num_faces; iface++) {
      /* Iterate over all faces that are also root faces and determine the face element */
      if (scheme->element_is_root_boundary (eclass, element, iface)) {
        /* Get face scheme */
        const int tree_face = scheme->element_get_tree_face (eclass, element, iface);

        /* Note: This wont work with non-default schemes, where the order of schemes is not the same as
         * in the default scheme. */
        const t8_eclass_t face_eclass = (t8_eclass_t) t8_eclass_face_types[tree_class][tree_face];

        t8_element_t *boundary;
        scheme->element_new (face_eclass, 1, &boundary);

        scheme->element_get_boundary_face (eclass, element, iface, boundary);

        scheme->element_extrude_face (eclass, boundary, check, tree_face);

        EXPECT_ELEM_EQ (scheme, eclass, element, check);

        scheme->element_destroy (face_eclass, 1, &boundary);
      }
    }
  }

 protected:
  void
  SetUp () override
  {
    dfs_test_setup ();
    /* Get element and initialize it */
    scheme->element_new (eclass, 1, &check);
    tree_class = scheme->get_eclass_scheme_eclass (eclass);
  }
  void
  TearDown () override
  {
    /* Destroy element */
    scheme->element_destroy (eclass, 1, &check);

    /* Destroy DFS test */
    dfs_test_teardown ();
  }
  t8_eclass_t tree_class;
  t8_element_t *check;
};

TEST_P (class_test_boundary_extrude, test_boundary_extrude_dfs)
{
#if T8_TEST_LEVEL_INT >= 1
  const int maxlvl = 4;
#else
  const int maxlvl = 6;
#endif
  check_recursive_dfs_to_max_lvl (maxlvl);
}

INSTANTIATE_TEST_SUITE_P (t8_gtest_test_all_imps, class_test_boundary_extrude, AllSchemes, print_all_schemes);