BoundaryMeshSubdomainTest Class Reference

Inheritance diagram for BoundaryMeshSubdomainTest:

Public Member Functions
	LIBMESH_CPPUNIT_TEST_SUITE (BoundaryMeshSubdomainTest)
	Tests the new_subdomain_ids parameter of add_elements(), which assigns subdomain IDs to newly-created lower-D elements.
	CPPUNIT_TEST (testDefaultSubdomain)
	CPPUNIT_TEST (testSingleSubdomain)
	CPPUNIT_TEST (testPerBoundarySubdomain)
	CPPUNIT_TEST_SUITE_END ()
void	testDefaultSubdomain ()
void	testSingleSubdomain ()
void	testPerBoundarySubdomain ()

Detailed Description

Definition at line 510 of file boundary_mesh.C.

Member Function Documentation

CPPUNIT_TEST() [1/3]

BoundaryMeshSubdomainTest::CPPUNIT_TEST

(

testDefaultSubdomain

)

CPPUNIT_TEST() [2/3]

BoundaryMeshSubdomainTest::CPPUNIT_TEST

(

testPerBoundarySubdomain

)

CPPUNIT_TEST() [3/3]

BoundaryMeshSubdomainTest::CPPUNIT_TEST

(

testSingleSubdomain

)

CPPUNIT_TEST_SUITE_END()

BoundaryMeshSubdomainTest::CPPUNIT_TEST_SUITE_END

(

)

LIBMESH_CPPUNIT_TEST_SUITE()

BoundaryMeshSubdomainTest::LIBMESH_CPPUNIT_TEST_SUITE

(

BoundaryMeshSubdomainTest

)

Tests the new_subdomain_ids parameter of add_elements(), which assigns subdomain IDs to newly-created lower-D elements.

Three cases:

1. Empty (default) - elements retain subdomain 0.
2. Single ID - all new elements get that ID.
3. Per-boundary - one-to-one with requested_boundary_ids (sorted).

Uses a 3x3 QUAD4 mesh. build_square() assigns: bottom=0, right=1, top=2, left=3.

testDefaultSubdomain()

void BoundaryMeshSubdomainTest::testDefaultSubdomain ( )

inline

Definition at line 533 of file boundary_mesh.C.

  {
    LOG_UNIT_TEST;
 
    Mesh mesh(*TestCommWorld);
    MeshTools::Generation::build_square(mesh, 3, 3, 0., 1., 0., 1., QUAD4);
 
    const BoundaryInfo & bi = mesh.get_boundary_info();
    const auto & side_boundary_ids = bi.get_side_boundary_ids();
 
    mesh.get_boundary_info().add_elements(side_boundary_ids, mesh);
    mesh.complete_preparation();
 
    for (const auto & elem : mesh.active_element_ptr_range())
      if (elem->dim() < mesh.mesh_dimension())
        CPPUNIT_ASSERT_EQUAL(static_cast<subdomain_id_type>(0), elem->subdomain_id());
  }

References libMesh::BoundaryInfo::add_elements(), libMesh::MeshTools::Generation::build_square(), libMesh::MeshBase::complete_preparation(), libMesh::MeshBase::get_boundary_info(), libMesh::BoundaryInfo::get_side_boundary_ids(), mesh, libMesh::MeshBase::mesh_dimension(), libMesh::QUAD4, and TestCommWorld.

testPerBoundarySubdomain()

void BoundaryMeshSubdomainTest::testPerBoundarySubdomain ( )

inline

Definition at line 579 of file boundary_mesh.C.

  {
    LOG_UNIT_TEST;
 
    Mesh mesh(*TestCommWorld);
    MeshTools::Generation::build_square(mesh, 3, 3, 0., 1., 0., 1., QUAD4);
 
    const BoundaryInfo & bi = mesh.get_boundary_info();
    const boundary_id_type right_id = bi.get_id_by_name("right");
    const boundary_id_type top_id   = bi.get_id_by_name("top");
 
    // The set is iterated in sorted boundary_id order, so the subdomain_ids
    // vector must be ordered accordingly.
    const subdomain_id_type right_sid = 10;
    const subdomain_id_type top_sid   = 20;
    std::set<boundary_id_type> ids = {right_id, top_id};
    std::vector<subdomain_id_type> sids;
    for (const auto id : ids)
      sids.push_back(id == right_id ? right_sid : top_sid);
 
    mesh.get_boundary_info().add_elements(ids, mesh, false, sids);
    mesh.complete_preparation();
 
    unsigned int n_right = 0, n_top = 0;
    for (const auto & elem : mesh.active_element_ptr_range())
      if (elem->dim() < mesh.mesh_dimension())
        {
          const Real cx = elem->vertex_average()(0);
          if (cx > 0.9)
            {
              CPPUNIT_ASSERT_EQUAL(right_sid, elem->subdomain_id());
              if (elem->processor_id() == mesh.processor_id())
                ++n_right;
            }
          else
            {
              CPPUNIT_ASSERT_EQUAL(top_sid, elem->subdomain_id());
              if (elem->processor_id() == mesh.processor_id())
                ++n_top;
            }
        }
 
    mesh.comm().sum(n_right);
    mesh.comm().sum(n_top);
    CPPUNIT_ASSERT_EQUAL(3u, n_right);
    CPPUNIT_ASSERT_EQUAL(3u, n_top);
  }

References libMesh::BoundaryInfo::add_elements(), libMesh::MeshTools::Generation::build_square(), libMesh::ParallelObject::comm(), libMesh::MeshBase::complete_preparation(), libMesh::MeshBase::get_boundary_info(), libMesh::BoundaryInfo::get_id_by_name(), mesh, libMesh::MeshBase::mesh_dimension(), libMesh::ParallelObject::processor_id(), libMesh::QUAD4, libMesh::Real, right_id, libMesh::Parallel::Communicator::sum(), TestCommWorld, and top_id.

testSingleSubdomain()

void BoundaryMeshSubdomainTest::testSingleSubdomain ( )

inline

Definition at line 551 of file boundary_mesh.C.

  {
    LOG_UNIT_TEST;
 
    Mesh mesh(*TestCommWorld);
    MeshTools::Generation::build_square(mesh, 3, 3, 0., 1., 0., 1., QUAD4);
 
    const BoundaryInfo & bi = mesh.get_boundary_info();
    const auto & side_boundary_ids = bi.get_side_boundary_ids();
 
    const subdomain_id_type sid = 42;
    mesh.get_boundary_info().add_elements(side_boundary_ids, mesh, false, {sid});
    mesh.complete_preparation();
 
    unsigned int n_lower = 0;
    for (const auto & elem : mesh.active_element_ptr_range())
      if (elem->dim() < mesh.mesh_dimension())
        {
          CPPUNIT_ASSERT_EQUAL(sid, elem->subdomain_id());
          if (elem->processor_id() == mesh.processor_id())
            ++n_lower;
        }
 
    // 3x3 quad mesh has 12 edges across its sides
    mesh.comm().sum(n_lower);
    CPPUNIT_ASSERT_EQUAL(12u, n_lower);
  }

References libMesh::BoundaryInfo::add_elements(), libMesh::MeshTools::Generation::build_square(), libMesh::ParallelObject::comm(), libMesh::MeshBase::complete_preparation(), libMesh::MeshBase::get_boundary_info(), libMesh::BoundaryInfo::get_side_boundary_ids(), mesh, libMesh::MeshBase::mesh_dimension(), libMesh::ParallelObject::processor_id(), libMesh::QUAD4, libMesh::Parallel::Communicator::sum(), and TestCommWorld.

---

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

• tests/mesh/boundary_mesh.C