AllTriTest Class Reference

Inheritance diagram for AllTriTest:

Public Member Functions
	LIBMESH_CPPUNIT_TEST_SUITE (AllTriTest)
	The goal of this test is to verify proper operation of the Mesh Extruder with the optional object callback for setting custom subdomain IDs. More...
	CPPUNIT_TEST (testAllTriTri)
	CPPUNIT_TEST (testAllTriQuad)
	CPPUNIT_TEST (testAllTriQuad8)
	CPPUNIT_TEST (testAllTriQuad9)
	CPPUNIT_TEST (testAllTriC0Polygon)
	CPPUNIT_TEST (testAllTriC0PolygonOctagon)
	CPPUNIT_TEST (testAllTriPrism6)
	CPPUNIT_TEST (testAllTriPrism18)
	CPPUNIT_TEST (testAllTriPrism20)
	CPPUNIT_TEST (testAllTriPrism21)
	CPPUNIT_TEST (testAllTriC0PolyhedronCube)
	CPPUNIT_TEST (testAllTriC0PolyhedronHexagonalPrism)
	CPPUNIT_TEST_SUITE_END ()
void	setUp ()
void	tearDown ()
void	testAllTriTri ()
void	testAllTriQuad ()
void	testAllTriQuad8 ()
void	testAllTriQuad9 ()
void	testAllTriPrism6 ()
void	testAllTriPrism18 ()
void	testAllTriPrism20 ()
void	testAllTriPrism21 ()
void	testAllTriC0Polygon ()
void	testAllTriC0PolygonOctagon ()
void	testAllTriC0PolyhedronCube ()
void	testAllTriC0PolyhedronHexagonalPrism ()

Protected Member Functions
void	test_helper_2D (ElemType elem_type, dof_id_type n_elem_expected, std::size_t n_boundary_conds_expected)
void	test_helper_3D (ElemType elem_type, dof_id_type n_elem_expected, std::size_t n_boundary_conds_expected)
void	test_helper_c0polyhedron (const std::vector< Point > &points, const std::vector< std::vector< unsigned int >> &nodes_on_side)

Detailed Description

Definition at line 21 of file all_tri.C.

Member Function Documentation

CPPUNIT_TEST() [1/12]

AllTriTest::CPPUNIT_TEST

(

testAllTriTri

)

CPPUNIT_TEST() [2/12]

AllTriTest::CPPUNIT_TEST

(

testAllTriQuad

)

CPPUNIT_TEST() [3/12]

AllTriTest::CPPUNIT_TEST

(

testAllTriQuad8

)

CPPUNIT_TEST() [4/12]

AllTriTest::CPPUNIT_TEST

(

testAllTriQuad9

)

CPPUNIT_TEST() [5/12]

AllTriTest::CPPUNIT_TEST

(

testAllTriC0Polygon

)

CPPUNIT_TEST() [6/12]

AllTriTest::CPPUNIT_TEST

(

testAllTriC0PolygonOctagon

)

CPPUNIT_TEST() [7/12]

AllTriTest::CPPUNIT_TEST

(

testAllTriPrism6

)

CPPUNIT_TEST() [8/12]

AllTriTest::CPPUNIT_TEST

(

testAllTriPrism18

)

CPPUNIT_TEST() [9/12]

AllTriTest::CPPUNIT_TEST

(

testAllTriPrism20

)

CPPUNIT_TEST() [10/12]

AllTriTest::CPPUNIT_TEST

(

testAllTriPrism21

)

CPPUNIT_TEST() [11/12]

AllTriTest::CPPUNIT_TEST

(

testAllTriC0PolyhedronCube

)

CPPUNIT_TEST() [12/12]

AllTriTest::CPPUNIT_TEST

(

testAllTriC0PolyhedronHexagonalPrism

)

CPPUNIT_TEST_SUITE_END()

AllTriTest::CPPUNIT_TEST_SUITE_END

(

)

LIBMESH_CPPUNIT_TEST_SUITE()

AllTriTest::LIBMESH_CPPUNIT_TEST_SUITE

(

AllTriTest

)

The goal of this test is to verify proper operation of the Mesh Extruder with the optional object callback for setting custom subdomain IDs.

We pass a custom object for generating subdomains based on the old element ID and the current layer and assert the proper values.

setUp()

void AllTriTest::setUp ( )

inline

Definition at line 105 of file all_tri.C.

{}

tearDown()

void AllTriTest::tearDown ( )

inline

Definition at line 107 of file all_tri.C.

{}

test_helper_2D()

void AllTriTest::test_helper_2D ( ElemType  elem_type, dof_id_type  n_elem_expected, std::size_t  n_boundary_conds_expected  )

inlineprotected

Definition at line 56 of file all_tri.C.

References libMesh::MeshTools::Modification::all_tri(), libMesh::MeshTools::Generation::build_square(), libMesh::MeshBase::get_boundary_info(), mesh, libMesh::BoundaryInfo::n_boundary_conds(), libMesh::MeshBase::n_elem(), and TestCommWorld.

  {
    ReplicatedMesh mesh(*TestCommWorld, /*dim=*/2);

    // Build a 2x1 TRI3 mesh and ask to split it into triangles.
    // Should be a no-op
    MeshTools::Generation::build_square(mesh,
                                        /*nx=*/2, /*ny=*/1,
                                        /*xmin=*/0., /*xmax=*/1.,
                                        /*ymin=*/0., /*ymax=*/1.,
                                        elem_type);

    MeshTools::Modification::all_tri(mesh);

    // Make sure that the expected number of elements is found.
    CPPUNIT_ASSERT_EQUAL(n_elem_expected, mesh.n_elem());

    // Make sure the expected number of BCs is found.
    CPPUNIT_ASSERT_EQUAL(n_boundary_conds_expected, mesh.get_boundary_info().n_boundary_conds());
  }

test_helper_3D()

void AllTriTest::test_helper_3D ( ElemType  elem_type, dof_id_type  n_elem_expected, std::size_t  n_boundary_conds_expected  )

inlineprotected

Definition at line 80 of file all_tri.C.

References libMesh::MeshTools::Modification::all_tri(), libMesh::MeshTools::Generation::build_cube(), libMesh::MeshBase::get_boundary_info(), mesh, libMesh::BoundaryInfo::n_boundary_conds(), libMesh::MeshBase::n_elem(), and TestCommWorld.

  {
    ReplicatedMesh mesh(*TestCommWorld, /*dim=*/3);

    // Build a 2x1 TRI3 mesh and ask to split it into triangles.
    // Should be a no-op
    MeshTools::Generation::build_cube(mesh,
                                      /*nx=*/1, /*ny=*/1, /*nz=*/1,
                                      /*xmin=*/0., /*xmax=*/1.,
                                      /*ymin=*/0., /*ymax=*/1.,
                                      /*zmin=*/0., /*zmax=*/1.,
                                      elem_type);

    MeshTools::Modification::all_tri(mesh);

    // Make sure that the expected number of elements is found.
    CPPUNIT_ASSERT_EQUAL(n_elem_expected, mesh.n_elem());

    // Make sure the expected number of BCs is found.
    CPPUNIT_ASSERT_EQUAL(n_boundary_conds_expected, mesh.get_boundary_info().n_boundary_conds());
  }

test_helper_c0polyhedron()

void AllTriTest::test_helper_c0polyhedron ( const std::vector< Point > &  points, const std::vector< std::vector< unsigned int >> &  nodes_on_side  )

inlineprotected

Definition at line 215 of file all_tri.C.

References libMesh::MeshBase::add_elem(), libMesh::MeshBase::add_node(), libMesh::MeshBase::add_point(), libMesh::BoundaryInfo::add_side(), libMesh::MeshTools::Modification::all_tri(), b, libMesh::MeshBase::get_boundary_info(), libMesh::BoundaryInfo::has_boundary_id(), libMesh::index_range(), libMesh::invalid_int, libMesh::make_range(), mesh, libMesh::BoundaryInfo::n_boundary_conds(), libMesh::MeshBase::n_elem(), libMesh::Elem::n_sides(), libMesh::Elem::n_sub_elem(), libMesh::MeshBase::node_ptr(), libMesh::MeshBase::prepare_for_use(), libMesh::Real, libMesh::DofObject::set_id(), TestCommWorld, libMesh::TET4, libMesh::TOLERANCE, and libMesh::MeshTools::volume().

  {
    ReplicatedMesh mesh(*TestCommWorld, /*dim=*/3);

    for (auto p : index_range(points))
      mesh.add_point(points[p], /*id=*/p);

    std::vector<std::shared_ptr<Polygon>> sides(nodes_on_side.size());
    for (auto s : index_range(nodes_on_side))
      {
        const auto & nodes_on_s = nodes_on_side[s];
        sides[s] = std::make_shared<C0Polygon>(nodes_on_s.size());
        for (auto i : index_range(nodes_on_s))
          sides[s]->set_node(i, mesh.node_ptr(nodes_on_s[i]));
      }

    std::unique_ptr<Node> mid_elem_node;
    std::unique_ptr<Elem> polyhedron =
      std::make_unique<C0Polyhedron>(sides, mid_elem_node);
    if (mid_elem_node)
      mesh.add_node(std::move(mid_elem_node));
    polyhedron->set_id() = 0;
    Elem * elem = mesh.add_elem(std::move(polyhedron));

    const auto * poly = cast_ptr<const C0Polyhedron *>(elem);
    const dof_id_type n_elem_expected = poly->n_subelements();

    // Mark every external face with a boundary id so we can verify
    // boundary information is transferred to the new tets.
    for (unsigned int s = 0; s < elem->n_sides(); ++s)
      mesh.get_boundary_info().add_side(elem, s, /*bnd_id=*/s);

    // The number of boundary triangles produced is the total number of
    // subtriangles across the polyhedron's polygonal faces.
    std::size_t n_bcs_expected = 0;
    for (unsigned int s = 0; s < elem->n_sides(); ++s)
      n_bcs_expected += sides[s]->n_subtriangles();

    // Keep this map for future checks
    std::vector<std::array<int, 4>> sub_elem_sides_to_parent_side(elem->n_sub_elem());
    for (unsigned int b = 0; b < elem->n_sub_elem(); ++b)
      sub_elem_sides_to_parent_side[b] = poly->subelement_sides_to_poly_sides(b);

    Real poly_volume = poly->volume();
    mesh.prepare_for_use();

    MeshTools::Modification::all_tri(mesh);

    CPPUNIT_ASSERT_EQUAL(n_elem_expected, mesh.n_elem());
    CPPUNIT_ASSERT_EQUAL(n_bcs_expected,
                         mesh.get_boundary_info().n_boundary_conds());

    for (const Elem * e : mesh.element_ptr_range())
      CPPUNIT_ASSERT_EQUAL(ElemType(TET4), e->type());

    // Check that the boundary info and numbering is as expected
    for (const Elem * e : mesh.element_ptr_range())
      for (const auto s : make_range(e->n_sides()))
      {
        // Get parent side
        const auto side = sub_elem_sides_to_parent_side[e->id()][s];
        // Check boundary exists on the tet side if the tet side is on a poly side
        if (side != invalid_int)
          CPPUNIT_ASSERT_EQUAL(mesh.get_boundary_info().has_boundary_id(e, s, side), true);
      }
    LIBMESH_ASSERT_NUMBERS_EQUAL(poly_volume, MeshTools::volume(mesh), TOLERANCE);
  }

testAllTriC0Polygon()

void AllTriTest::testAllTriC0Polygon ( )

inline

Definition at line 129 of file all_tri.C.

References libMesh::MeshTools::Modification::all_tri(), libMesh::MeshTools::Generation::build_square(), libMesh::C0POLYGON, libMesh::MeshBase::get_boundary_info(), mesh, libMesh::BoundaryInfo::n_boundary_conds(), libMesh::MeshBase::n_elem(), libMesh::Polygon::n_subtriangles(), TestCommWorld, and libMesh::TRI3.

  {
    LOG_UNIT_TEST;

    ReplicatedMesh mesh(*TestCommWorld, /*dim=*/2);

    MeshTools::Generation::build_square(mesh,
                                        /*nx=*/2, /*ny=*/2,
                                        /*xmin=*/0., /*xmax=*/1.,
                                        /*ymin=*/0., /*ymax=*/1.,
                                        C0POLYGON);

    // The post-all_tri element count is the sum of the per-polygon
    // subtriangle counts, and external sides should be preserved one
    // for one as TRI3 sides.
    dof_id_type n_elem_expected = 0;
    for (const Elem * elem : mesh.element_ptr_range())
      {
        const Polygon * poly = dynamic_cast<const Polygon *>(elem);
        CPPUNIT_ASSERT(poly != nullptr);
        n_elem_expected += poly->n_subtriangles();
      }

    const std::size_t n_bcs_before =
      mesh.get_boundary_info().n_boundary_conds();

    MeshTools::Modification::all_tri(mesh);

    CPPUNIT_ASSERT_EQUAL(n_elem_expected, mesh.n_elem());
    CPPUNIT_ASSERT_EQUAL(n_bcs_before,
                         mesh.get_boundary_info().n_boundary_conds());

    for (const Elem * elem : mesh.element_ptr_range())
      CPPUNIT_ASSERT_EQUAL(ElemType(TRI3), elem->type());
  }

testAllTriC0PolygonOctagon()

void AllTriTest::testAllTriC0PolygonOctagon ( )

inline

Definition at line 168 of file all_tri.C.

References libMesh::MeshBase::add_elem(), libMesh::MeshBase::add_point(), libMesh::BoundaryInfo::add_side(), libMesh::MeshTools::Modification::all_tri(), libMesh::MeshBase::get_boundary_info(), mesh, libMesh::BoundaryInfo::n_boundary_conds(), libMesh::MeshBase::n_elem(), libMesh::pi, libMesh::MeshBase::prepare_for_use(), libMesh::Real, libMesh::DofObject::set_id(), TestCommWorld, libMesh::TOLERANCE, libMesh::TRI3, and libMesh::MeshTools::volume().

  {
    LOG_UNIT_TEST;

    constexpr unsigned int n_sides = 8;

    ReplicatedMesh mesh(*TestCommWorld, /*dim=*/2);

    std::unique_ptr<Elem> octagon = std::make_unique<C0Polygon>(n_sides);
    for (unsigned int i = 0; i < n_sides; ++i)
      {
        const Real angle = 2 * libMesh::pi * i / n_sides;
        Node * node = mesh.add_point(Point(std::cos(angle), std::sin(angle), 0.),
                                     /*id=*/i);
        octagon->set_node(i, node);
      }
    octagon->set_id() = 0;
    Real poly_volume = octagon->volume();
    Elem * elem = mesh.add_elem(std::move(octagon));

    // Mark every external side with a boundary id so we can verify
    // boundary information is transferred to the new triangles.
    for (unsigned int s = 0; s < n_sides; ++s)
      mesh.get_boundary_info().add_side(elem, s, /*bnd_id=*/0);

    mesh.prepare_for_use();

    MeshTools::Modification::all_tri(mesh);

    // n_sides - 2 = 6 subtriangles
    CPPUNIT_ASSERT_EQUAL(dof_id_type(n_sides - 2), mesh.n_elem());
    CPPUNIT_ASSERT_EQUAL(std::size_t(n_sides),
                         mesh.get_boundary_info().n_boundary_conds());

    for (const Elem * e : mesh.element_ptr_range())
      CPPUNIT_ASSERT_EQUAL(ElemType(TRI3), e->type());

    LIBMESH_ASSERT_NUMBERS_EQUAL(poly_volume, MeshTools::volume(mesh), TOLERANCE);
  }

testAllTriC0PolyhedronCube()

void AllTriTest::testAllTriC0PolyhedronCube ( )

inline

Definition at line 288 of file all_tri.C.

  {
    LOG_UNIT_TEST;

    const std::vector<Point> points =
      { {0,0,0}, {1,0,0}, {1,1,0}, {0,1,0},
        {0,0,1}, {1,0,1}, {1,1,1}, {0,1,1} };

    const std::vector<std::vector<unsigned int>> nodes_on_side =
      { {0, 1, 2, 3},   // min z
        {0, 1, 5, 4},   // min y
        {2, 6, 5, 1},   // max x
        {2, 3, 7, 6},   // max y
        {0, 4, 7, 3},   // min x
        {5, 6, 7, 4} }; // max z

    test_helper_c0polyhedron(points, nodes_on_side);
  }

testAllTriC0PolyhedronHexagonalPrism()

void AllTriTest::testAllTriC0PolyhedronHexagonalPrism ( )

inline

Definition at line 309 of file all_tri.C.

  {
    LOG_UNIT_TEST;

    const std::vector<Point> points =
      { { 0, -2, 0}, {-1, -1, 0}, {-1, 1, 0},
        { 0,  2, 0}, { 1,  1, 0}, { 1, -1, 0},
        { 0, -2, 1}, {-1, -1, 1}, {-1, 1, 1},
        { 0,  2, 1}, { 1,  1, 1}, { 1, -1, 1} };

    const std::vector<std::vector<unsigned int>> nodes_on_side =
      { {0, 1, 2, 3, 4, 5},
        {0, 1,  7,  6},
        {1, 2,  8,  7},
        {2, 3,  9,  8},
        {3, 4, 10,  9},
        {4, 5, 11, 10},
        {5, 0,  6, 11},
        {6, 7,  8,  9, 10, 11} };

    test_helper_c0polyhedron(points, nodes_on_side);
  }

testAllTriPrism18()

void AllTriTest::testAllTriPrism18 ( )

inline

Definition at line 123 of file all_tri.C.

References libMesh::PRISM18.

{ LOG_UNIT_TEST; test_helper_3D(PRISM18, /*nelem=*/6, /*nbcs=*/12); }

testAllTriPrism20()

void AllTriTest::testAllTriPrism20 ( )

inline

Definition at line 124 of file all_tri.C.

References libMesh::PRISM20.

{ LOG_UNIT_TEST; test_helper_3D(PRISM20, /*nelem=*/6, /*nbcs=*/12); }

testAllTriPrism21()

void AllTriTest::testAllTriPrism21 ( )

inline

Definition at line 125 of file all_tri.C.

References libMesh::PRISM21.

{ LOG_UNIT_TEST; test_helper_3D(PRISM21, /*nelem=*/6, /*nbcs=*/12); }

testAllTriPrism6()

void AllTriTest::testAllTriPrism6 ( )

inline

Definition at line 122 of file all_tri.C.

References libMesh::PRISM6.

{ LOG_UNIT_TEST; test_helper_3D(PRISM6, /*nelem=*/6, /*nbcs=*/12); }

testAllTriQuad()

void AllTriTest::testAllTriQuad ( )

inline

Definition at line 113 of file all_tri.C.

References libMesh::QUAD4.

{ LOG_UNIT_TEST; test_helper_2D(QUAD4, /*nelem=*/4, /*nbcs=*/6); }

testAllTriQuad8()

void AllTriTest::testAllTriQuad8 ( )

inline

Definition at line 116 of file all_tri.C.

References libMesh::QUAD8.

{ LOG_UNIT_TEST; test_helper_2D(QUAD8, /*nelem=*/4, /*nbcs=*/6); }

testAllTriQuad9()

void AllTriTest::testAllTriQuad9 ( )

inline

Definition at line 119 of file all_tri.C.

References libMesh::QUAD9.

{ LOG_UNIT_TEST; test_helper_2D(QUAD9, /*nelem=*/4, /*nbcs=*/6); }

testAllTriTri()

void AllTriTest::testAllTriTri ( )

inline

Definition at line 110 of file all_tri.C.

References libMesh::TRI3.

{ LOG_UNIT_TEST; test_helper_2D(TRI3, /*nelem=*/4, /*nbcs=*/6); }

---

The documentation for this class was generated from the following file:

• tests/mesh/all_tri.C