WriteSidesetData Class Reference

Inheritance diagram for WriteSidesetData:

Public Member Functions
	CPPUNIT_TEST_SUITE (WriteSidesetData)
	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_sideset_data.C.

Member Function Documentation

CPPUNIT_TEST()

WriteSidesetData::CPPUNIT_TEST

(

testWrite

)

CPPUNIT_TEST_SUITE()

WriteSidesetData::CPPUNIT_TEST_SUITE

(

WriteSidesetData

)

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

CPPUNIT_TEST_SUITE_END()

WriteSidesetData::CPPUNIT_TEST_SUITE_END

(

)

testWrite()

void WriteSidesetData::testWrite ( )

inline

Definition at line 31 of file write_sideset_data.C.

  {
    Mesh mesh(*TestCommWorld);
 
    // We set our initial conditions based on build_square node ids
    mesh.allow_renumbering(false);
 
    MeshTools::Generation::build_square(mesh,
                                        /*nx=*/5, /*ny=*/5,
                                        -1., 1.,
                                        -1., 1.,
                                        QUAD4);
 
    // Get list of all (elem, side, id) tuples
    std::vector<BoundaryInfo::BCTuple> all_bc_tuples =
      mesh.get_boundary_info().build_side_list();
 
    // Data structures to be passed to ExodusII_IO::write_sideset_data
    std::vector<std::string> var_names = {"var1", "var2"};
    std::vector<std::set<boundary_id_type>> side_ids =
      {
        {0, 2}, // var1 is defined on sidesets 0 and 2
        {1, 3}  // var2 is defined on sidesets 1 and 3
      };
 
    // Data structure mapping (elem, side, id) tuples to Real values that
    // will be passed to Exodus.
    std::vector<std::map<BoundaryInfo::BCTuple, Real>> bc_vals(var_names.size());
 
    // For each var_names[i], construct bc_vals[i]
    for (unsigned int i=0; i<var_names.size(); ++i)
      {
        // const auto & var_name = var_names[i];
        auto & vals = bc_vals[i];
 
        for (const auto & t : all_bc_tuples)
          {
            // dof_id_type elem_id = std::get<0>(t);
            // unsigned int side_id = std::get<1>(t);
            boundary_id_type b_id = std::get<2>(t);
 
            if (side_ids[i].count(b_id))
              {
                // Compute a value. This could in theory depend on
                // var_name, elem_id, side_id, and/or b_id.
                Real val = static_cast<Real>(b_id);
 
                // Insert into the vals map.
                vals.insert(std::make_pair(t, val));
              }
          }
 
        // If we have a distributed mesh, write_sideset_data wants our
        // ghost data too; we'll just serialize everything here.
        if (!mesh.is_serial())
          TestCommWorld->set_union(vals);
 
      } // done constructing bc_vals
 
#ifdef LIBMESH_HAVE_EXODUS_API
 
    // We write the file in the ExodusII format.
    {
      ExodusII_IO writer(mesh);
      writer.write("write_sideset_data.e");
      writer.write_sideset_data (/*timestep=*/1, var_names, side_ids, bc_vals);
    }
 
    // Make sure that the writing is done before the reading starts.
    TestCommWorld->barrier();
 
    // Now read it back in
    Mesh read_mesh(*TestCommWorld);
    ExodusII_IO reader(read_mesh);
    reader.read("write_sideset_data.e");
 
    std::vector<std::string> read_in_var_names;
    std::vector<std::set<boundary_id_type>> read_in_side_ids;
    std::vector<std::map<BoundaryInfo::BCTuple, Real>> read_in_bc_vals;
    reader.read_sideset_data
      (/*timestep=*/1, read_in_var_names, read_in_side_ids, read_in_bc_vals);
 
    // Assert that we got back out what we put in.
    CPPUNIT_ASSERT(read_in_var_names == var_names);
    CPPUNIT_ASSERT(read_in_side_ids == side_ids);
    CPPUNIT_ASSERT(read_in_bc_vals == bc_vals);
 
#endif // #ifdef LIBMESH_HAVE_EXODUS_API
  }

References libMesh::MeshBase::allow_renumbering(), libMesh::BoundaryInfo::build_side_list(), libMesh::MeshTools::Generation::build_square(), libMesh::MeshBase::get_boundary_info(), libMesh::MeshBase::is_serial(), mesh, libMesh::QUAD4, libMesh::ExodusII_IO::read(), libMesh::ExodusII_IO::read_sideset_data(), libMesh::Real, TestCommWorld, libMesh::ExodusII_IO::write(), and libMesh::ExodusII_IO::write_sideset_data().

---

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

• tests/mesh/write_sideset_data.C