SlitMeshTest Class Reference

Inheritance diagram for SlitMeshTest:

Public Member Functions
	CPPUNIT_TEST_SUITE (SlitMeshTest)
	The goal of this test is to ensure that a 2D mesh with nodes overlapping on opposite sides of an internal, "slit" edge is usable. More...
	CPPUNIT_TEST (testMesh)
	CPPUNIT_TEST_SUITE_END ()
void	setUp ()
void	tearDown ()
void	testMesh ()

Protected Member Functions
void	build_mesh ()

Protected Attributes
Mesh *	_mesh

Detailed Description

Definition at line 61 of file slit_mesh_test.C.

Member Function Documentation

build_mesh()

void SlitMeshTest::build_mesh ( )

inlineprotected

Definition at line 79 of file slit_mesh_test.C.

  {
    _mesh = new Mesh(*TestCommWorld);
 
    // (0,1)           (1,1)           (2,1)
    // x---------------x---------------x
    // |               |               |
    // |               |               |
    // |               |               |
    // |               |               |
    // |               |               |
    // x---------------x---------------x
    // (0,0)           (1,0)          (2,0)
    // |               |               |
    // |               |               |
    // |               |               |
    // |               |               |
    // x---------------x---------------x
    // (0,-1)          (1,-1)         (2,-1)
 
    _mesh->set_mesh_dimension(2);
 
    _mesh->add_point( Point(0.0, 0.0), 0 );
    _mesh->add_point( Point(1.0, 0.0), 1 );
    _mesh->add_point( Point(1.0, 1.0), 2 );
    _mesh->add_point( Point(0.0, 1.0), 3 );
    _mesh->add_point( Point(0.0,-1.0), 4 );
    _mesh->add_point( Point(1.0,-1.0), 5 );
    _mesh->add_point( Point(1.0, 0.0), 6 );
    _mesh->add_point( Point(2.0, 0.0), 7 );
    _mesh->add_point( Point(2.0, 1.0), 8 );
    _mesh->add_point( Point(2.0,-1.0), 9 );
 
    {
      Elem* elem_top_left = new Quad4;
      elem_top_left->set_node(0) = _mesh->node_ptr(0);
      elem_top_left->set_node(1) = _mesh->node_ptr(1);
      elem_top_left->set_node(2) = _mesh->node_ptr(2);
      elem_top_left->set_node(3) = _mesh->node_ptr(3);
      elem_top_left->set_id() = 0;
      _mesh->add_elem(elem_top_left);
 
      Elem* elem_bottom_left = new Quad4;
      elem_bottom_left->set_node(0) = _mesh->node_ptr(4);
      elem_bottom_left->set_node(1) = _mesh->node_ptr(5);
      elem_bottom_left->set_node(2) = _mesh->node_ptr(6);
      elem_bottom_left->set_node(3) = _mesh->node_ptr(0);
      elem_bottom_left->set_id() = 1;
      _mesh->add_elem(elem_bottom_left);
 
      Elem* elem_top_right = new Quad4;
      elem_top_right->set_node(0) = _mesh->node_ptr(1);
      elem_top_right->set_node(1) = _mesh->node_ptr(7);
      elem_top_right->set_node(2) = _mesh->node_ptr(8);
      elem_top_right->set_node(3) = _mesh->node_ptr(2);
      elem_top_right->set_id() = 2;
      _mesh->add_elem(elem_top_right);
 
      Elem* elem_bottom_right = new Quad4;
      elem_bottom_right->set_node(0) = _mesh->node_ptr(5);
      elem_bottom_right->set_node(1) = _mesh->node_ptr(9);
      elem_bottom_right->set_node(2) = _mesh->node_ptr(7);
      elem_bottom_right->set_node(3) = _mesh->node_ptr(6);
      elem_bottom_right->set_id() = 3;
      _mesh->add_elem(elem_bottom_right);
    }
 
    // libMesh shouldn't renumber, or our based-on-initial-id
    // assertions later may fail.
    _mesh->allow_renumbering(false);
 
    _mesh->prepare_for_use();
  }

References libMesh::DistributedMesh::add_elem(), libMesh::DistributedMesh::add_point(), libMesh::MeshBase::allow_renumbering(), libMesh::DistributedMesh::node_ptr(), libMesh::MeshBase::prepare_for_use(), libMesh::DofObject::set_id(), libMesh::MeshBase::set_mesh_dimension(), libMesh::Elem::set_node(), and TestCommWorld.

CPPUNIT_TEST()

SlitMeshTest::CPPUNIT_TEST

(

testMesh

)

CPPUNIT_TEST_SUITE()

SlitMeshTest::CPPUNIT_TEST_SUITE

(

SlitMeshTest

)

The goal of this test is to ensure that a 2D mesh with nodes overlapping on opposite sides of an internal, "slit" edge is usable.

CPPUNIT_TEST_SUITE_END()

SlitMeshTest::CPPUNIT_TEST_SUITE_END

(

)

setUp()

void SlitMeshTest::setUp ( )

inline

Definition at line 154 of file slit_mesh_test.C.

  {
#if LIBMESH_DIM > 1
    this->build_mesh();
#endif
  }

tearDown()

void SlitMeshTest::tearDown ( )

inline

Definition at line 161 of file slit_mesh_test.C.

  {
    delete _mesh;
  }

testMesh()

void SlitMeshTest::testMesh ( )

inline

Definition at line 166 of file slit_mesh_test.C.

  {
    // There'd better be 4 elements
    CPPUNIT_ASSERT_EQUAL( (dof_id_type)4, _mesh->n_elem() );
 
    // There'd better still be a full 10 nodes
    CPPUNIT_ASSERT_EQUAL( (dof_id_type)10, _mesh->n_nodes() );
 
    /* The middle nodes should still be distinct between the top and
     * bottom elements */
    if (_mesh->query_elem_ptr(0) && _mesh->query_elem_ptr(1))
      CPPUNIT_ASSERT( _mesh->elem_ref(0).node_id(1) != _mesh->elem_ref(1).node_id(2) );
    if (_mesh->query_elem_ptr(2) && _mesh->query_elem_ptr(3))
      CPPUNIT_ASSERT( _mesh->elem_ref(2).node_id(0) != _mesh->elem_ref(3).node_id(3) );
 
    /* The middle nodes should still be shared between left and right
     * elements on top and bottom */
    if (_mesh->query_elem_ptr(0) && _mesh->query_elem_ptr(2))
      CPPUNIT_ASSERT_EQUAL( _mesh->elem_ref(0).node_id(1),
                            _mesh->elem_ref(2).node_id(0) );
    if (_mesh->query_elem_ptr(1) && _mesh->query_elem_ptr(3))
      CPPUNIT_ASSERT_EQUAL( _mesh->elem_ref(1).node_id(2),
                            _mesh->elem_ref(3).node_id(3) );
  }

References libMesh::MeshBase::elem_ref(), libMesh::DistributedMesh::n_elem(), libMesh::DistributedMesh::n_nodes(), libMesh::Elem::node_id(), and libMesh::DistributedMesh::query_elem_ptr().

Member Data Documentation

_mesh

Mesh* SlitMeshTest::_mesh

protected

Definition at line 77 of file slit_mesh_test.C.

---

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

• tests/mesh/slit_mesh_test.C