WriteEdgesetData Class Reference

Inheritance diagram for WriteEdgesetData:

Public Member Functions
	CPPUNIT_TEST_SUITE (WriteEdgesetData)
	This test ensures you can write both vector and scalar variables. More...
	CPPUNIT_TEST (testWrite)
	CPPUNIT_TEST_SUITE_END ()
void	testWrite ()

Detailed Description

Definition at line 17 of file write_edgeset_data.C.

Member Function Documentation

CPPUNIT_TEST()

WriteEdgesetData::CPPUNIT_TEST

(

testWrite

)

CPPUNIT_TEST_SUITE()

WriteEdgesetData::CPPUNIT_TEST_SUITE

(

WriteEdgesetData

)

This test ensures you can write both vector and scalar variables.

CPPUNIT_TEST_SUITE_END()

WriteEdgesetData::CPPUNIT_TEST_SUITE_END

(

)

testWrite()

void WriteEdgesetData::testWrite ( )

inline

Definition at line 31 of file write_edgeset_data.C.

  {
    ReplicatedMesh mesh(*TestCommWorld);
 
    MeshTools::Generation::build_cube (mesh,
                                       5, 5, 5,
                                       0., 1.,
                                       0., 1.,
                                       0., 1.,
                                       HEX8);
 
    // Add named edgesets containing the "front" and "back" edges of
    // the cube.
    BoundaryInfo & bi = mesh.get_boundary_info();
 
    std::vector<BoundaryInfo::BCTuple> bc_tuples =
      bi.build_side_list();
 
    for (const auto & t : bc_tuples)
      {
        dof_id_type elem_id = std::get<0>(t);
        boundary_id_type b_id = std::get<2>(t);
 
        // Side 0 corresponds to the "back" sideset, with edges {0,1,2,3}
        // Side 5 corresponds to the "front" sideset, with edges {4,5,6,7}
        // If we are on either of those sides, add all the corresponding edges.
        if (b_id == 0)
          for (unsigned short edge=0; edge<4; ++edge)
            bi.add_edge(elem_id, edge, b_id);
 
        if (b_id == 5)
          for (unsigned short edge=8; edge<12; ++edge)
            bi.add_edge(elem_id, edge, b_id);
      }
 
    // Name the edgesets
    bi.edgeset_name(0) = "back_edgeset";
    bi.edgeset_name(5) = "front_edgeset";
 
#ifdef LIBMESH_HAVE_EXODUS_API
 
    // We write the file in the ExodusII format.
    {
      ExodusII_IO writer(mesh);
      writer.write("write_edgeset_data.e");
    }
 
    // Make sure that the writing is done before the reading starts.
    TestCommWorld->barrier();
 
    // Now read it back in
    ReplicatedMesh read_mesh(*TestCommWorld);
    ExodusII_IO reader(read_mesh);
    reader.read("write_edgeset_data.e");
 
    // Assert that we got back out what we put in.
    BoundaryInfo & read_bi = read_mesh.get_boundary_info();
    CPPUNIT_ASSERT(read_bi.edgeset_name(0) == "back_edgeset");
    CPPUNIT_ASSERT(read_bi.edgeset_name(5) == "front_edgeset");
    CPPUNIT_ASSERT(read_bi.n_edge_conds() == 200);
    CPPUNIT_ASSERT(read_bi.get_edge_boundary_ids().size() == 2);
 
    // Make sure that we have the expected amount of edges in each edgeset.
    std::map<boundary_id_type, unsigned int> counts;
    std::vector<BoundaryInfo::BCTuple> edge_tuples = bi.build_edge_list();
    for (const auto & t : edge_tuples)
      counts[std::get<2>(t)]++;
    CPPUNIT_ASSERT(counts[0] == 100);
    CPPUNIT_ASSERT(counts[5] == 100);
 
#endif // #ifdef LIBMESH_HAVE_EXODUS_API
  }

References libMesh::BoundaryInfo::add_edge(), libMesh::MeshTools::Generation::build_cube(), libMesh::BoundaryInfo::build_edge_list(), libMesh::BoundaryInfo::build_side_list(), libMesh::BoundaryInfo::edgeset_name(), libMesh::MeshBase::get_boundary_info(), libMesh::BoundaryInfo::get_edge_boundary_ids(), libMesh::HEX8, mesh, libMesh::BoundaryInfo::n_edge_conds(), libMesh::ExodusII_IO::read(), TestCommWorld, and libMesh::ExodusII_IO::write().

---

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

• tests/mesh/write_edgeset_data.C