MeshfunctionDFEM Class Reference

Inheritance diagram for MeshfunctionDFEM:

Public Member Functions
	CPPUNIT_TEST_SUITE (MeshfunctionDFEM)
	The goal of this test is to verify proper operation of mesh function for discontinuous shape functions and points close to the boundary. More...
	CPPUNIT_TEST (test_point_locator_dfem)
	CPPUNIT_TEST (test_mesh_function_dfem)
	CPPUNIT_TEST (test_mesh_function_dfem_grad)
	CPPUNIT_TEST_SUITE_END ()
void	setUp ()
void	tearDown ()
void	test_point_locator_dfem ()
void	test_mesh_function_dfem ()
void	test_mesh_function_dfem_grad ()

Protected Member Functions
void	build_mesh ()

Protected Attributes
ReplicatedMesh *	_mesh
std::unique_ptr< PointLocatorBase >	_point_locator

Detailed Description

Definition at line 42 of file mesh_function_dfem.C.

Member Function Documentation

build_mesh()

void MeshfunctionDFEM::build_mesh ( )

inlineprotected

Definition at line 65 of file mesh_function_dfem.C.

  {
    _mesh = new ReplicatedMesh(*TestCommWorld);
 
    // (0,1)           (1,1)
    // x---------------x
    // |               |
    // |               |
    // |               |
    // |               |
    // |               |
    // x---------------x
    // (0,0)           (1,0)
    // |               |
    // |               |
    // |               |
    // |               |
    // x---------------x
    // (0,-1)          (1,-1)
 
    _mesh->set_mesh_dimension(2);
 
    _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), 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, 0.0), 0 );
 
    {
      Elem* elem_top = _mesh->add_elem( new Quad4 );
      elem_top->set_node(0) = _mesh->node_ptr(0);
      elem_top->set_node(1) = _mesh->node_ptr(1);
      elem_top->set_node(2) = _mesh->node_ptr(2);
      elem_top->set_node(3) = _mesh->node_ptr(3);
 
      Elem* elem_bottom = _mesh->add_elem( new Quad4 );
      elem_bottom->set_node(0) = _mesh->node_ptr(4);
      elem_bottom->set_node(1) = _mesh->node_ptr(5);
      elem_bottom->set_node(2) = _mesh->node_ptr(1);
      elem_bottom->set_node(3) = _mesh->node_ptr(0);
    }
 
    // libMesh will renumber, but we numbered according to its scheme
    // anyway. We do this because when we call uniformly_refine subsequently,
    // it's going use skip_renumber=false.
    _mesh->prepare_for_use(false /*skip_renumber*/);
 
    // get a point locator
    _point_locator = _mesh->sub_point_locator();
  }

References libMesh::ReplicatedMesh::add_elem(), libMesh::ReplicatedMesh::add_point(), libMesh::ReplicatedMesh::node_ptr(), libMesh::MeshBase::prepare_for_use(), libMesh::MeshBase::set_mesh_dimension(), libMesh::Elem::set_node(), libMesh::MeshBase::sub_point_locator(), and TestCommWorld.

CPPUNIT_TEST() [1/3]

MeshfunctionDFEM::CPPUNIT_TEST

(

test_mesh_function_dfem

)

CPPUNIT_TEST() [2/3]

MeshfunctionDFEM::CPPUNIT_TEST

(

test_mesh_function_dfem_grad

)

CPPUNIT_TEST() [3/3]

MeshfunctionDFEM::CPPUNIT_TEST

(

test_point_locator_dfem

)

CPPUNIT_TEST_SUITE()

MeshfunctionDFEM::CPPUNIT_TEST_SUITE

(

MeshfunctionDFEM

)

The goal of this test is to verify proper operation of mesh function for discontinuous shape functions and points close to the boundary.

CPPUNIT_TEST_SUITE_END()

MeshfunctionDFEM::CPPUNIT_TEST_SUITE_END

(

)

setUp()

void MeshfunctionDFEM::setUp ( )

inline

Definition at line 118 of file mesh_function_dfem.C.

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

tearDown()

void MeshfunctionDFEM::tearDown ( )

inline

Definition at line 125 of file mesh_function_dfem.C.

  {
#if LIBMESH_DIM > 1
    delete _mesh;
#endif
  }

test_mesh_function_dfem()

void MeshfunctionDFEM::test_mesh_function_dfem ( )

inline

Definition at line 164 of file mesh_function_dfem.C.

  {
    // this point is in top element only
    Point top(0.5, 0.5, 0.0);
 
    // this point is in bottom element only
    Point bottom(0.5, -0.5, 0.0);
 
    // this point is on the face between bottom and top
    Point face(0.5, 0.0, 0.0);
 
    // set up some necessary objects
    EquationSystems es(*_mesh);
    System & sys = es.add_system<System> ("SimpleSystem");
    sys.add_variable("u", CONSTANT, MONOMIAL);
 
    es.init();
    sys.project_solution(position_function, nullptr, es.parameters);
 
    std::vector<unsigned int> variables;
    sys.get_all_variable_numbers(variables);
    std::unique_ptr<NumericVector<Number>> mesh_function_vector =
      NumericVector<Number>::build(es.comm());
    mesh_function_vector->init(sys.n_dofs(), false, SERIAL);
    sys.solution->localize(*mesh_function_vector);
 
    MeshFunction mesh_function (es, *mesh_function_vector,
                                sys.get_dof_map(), variables);
    mesh_function.init();
 
    // test mesh function in top
    std::map<const Elem *, Number> top_val = mesh_function.discontinuous_value(top);
 
    // check that there is only one value
    CPPUNIT_ASSERT (top_val.size() == 1);
 
    // check that this one value is the right one
    for (std::map<const Elem *, Number>::const_iterator it = top_val.begin(); it != top_val.end(); ++it)
      CPPUNIT_ASSERT (it->first->id() == 0 && std::abs(it->second) < 1.0e-10);
 
    // test mesh function in bottom
    std::map<const Elem *, Number> bottom_val = mesh_function.discontinuous_value(bottom);
 
    // check that there is only one value
    CPPUNIT_ASSERT (bottom_val.size() == 1);
 
    // check that this one value is the right one
    for (std::map<const Elem *, Number>::const_iterator it = bottom_val.begin(); it != bottom_val.end(); ++it)
      CPPUNIT_ASSERT (it->first->id() == 1 && std::abs(it->second - 1.0) < 1.0e-10);
 
    // test mesh function in face
    std::map<const Elem *, Number> face_val = mesh_function.discontinuous_value(face);
 
    // check that there are two values
    CPPUNIT_ASSERT (face_val.size() == 2);
 
    // check that the two values are attached to the right element
    for (std::map<const Elem *, Number>::const_iterator it = face_val.begin(); it != face_val.end(); ++it)
      if (it->first->id() == 0)
        CPPUNIT_ASSERT (std::abs(it->second) < 1.0e-10);
      else if (it->first->id() == 1)
        CPPUNIT_ASSERT (std::abs(it->second - 1.0) < 1.0e-10);
      else
        CPPUNIT_ASSERT (false);
 
  }

References std::abs(), libMesh::EquationSystems::add_system(), libMesh::System::add_variable(), libMesh::NumericVector< T >::build(), libMesh::ParallelObject::comm(), libMesh::CONSTANT, libMesh::System::get_all_variable_numbers(), libMesh::System::get_dof_map(), libMesh::MeshFunction::init(), libMesh::EquationSystems::init(), libMesh::MONOMIAL, libMesh::System::n_dofs(), libMesh::EquationSystems::parameters, position_function(), libMesh::System::project_solution(), libMesh::SERIAL, and libMesh::System::solution.

test_mesh_function_dfem_grad()

void MeshfunctionDFEM::test_mesh_function_dfem_grad ( )

inline

Definition at line 232 of file mesh_function_dfem.C.

  {
    // this point is in top element only
    Point top(0.5, 0.5, 0.0);
 
    // this point is in bottom element only
    Point bottom(0.5, -0.5, 0.0);
 
    // this point is on the face between bottom and top
    Point face(0.5, 0.0, 0.0);
 
    // set up some necessary objects
    EquationSystems es(*_mesh);
    System & sys = es.add_system<System> ("SimpleSystem");
    sys.add_variable("u", FIRST, LAGRANGE);
 
    es.init();
    sys.project_solution(position_function2, nullptr, es.parameters);
 
    std::vector<unsigned int> variables;
    sys.get_all_variable_numbers(variables);
    std::unique_ptr<NumericVector<Number>> mesh_function_vector =
      NumericVector<Number>::build(es.comm());
    mesh_function_vector->init(sys.n_dofs(), false, SERIAL);
    sys.solution->localize(*mesh_function_vector);
 
    MeshFunction mesh_function (es, *mesh_function_vector,
                                sys.get_dof_map(), variables);
    mesh_function.init();
 
    // test mesh function in top
    std::map<const Elem *, Gradient> top_val = mesh_function.discontinuous_gradient(top);
 
    // check that there is only one value
    CPPUNIT_ASSERT (top_val.size() == 1);
 
    // check that this one value is the right one
    for (std::map<const Elem *, Gradient>::const_iterator it = top_val.begin(); it != top_val.end(); ++it)
      CPPUNIT_ASSERT (it->first->id() == 0 && std::abs(it->second(1) - 2.0) < 1.0e-10);
 
    // test mesh function in bottom
    std::map<const Elem *, Gradient> bottom_val = mesh_function.discontinuous_gradient(bottom);
 
    // check that there is only one value
    CPPUNIT_ASSERT (bottom_val.size() == 1);
 
    // check that this one value is the right one
    for (std::map<const Elem *, Gradient>::const_iterator it = bottom_val.begin(); it != bottom_val.end(); ++it)
      CPPUNIT_ASSERT (it->first->id() == 1 && std::abs(it->second(1) - 1.0) < 1.0e-10);
 
    // test mesh function in face
    std::map<const Elem *, Gradient> face_val = mesh_function.discontinuous_gradient(face);
 
    // check that there are two values
    CPPUNIT_ASSERT (face_val.size() == 2);
 
    // check that the two values are attached to the right element
    for (std::map<const Elem *, Gradient>::const_iterator it = face_val.begin(); it != face_val.end(); ++it)
      if (it->first->id() == 0)
        CPPUNIT_ASSERT (std::abs(it->second(1) - 2.0) < 1.0e-10);
      else if (it->first->id() == 1)
        CPPUNIT_ASSERT (std::abs(it->second(1) - 1.0) < 1.0e-10);
      else
        CPPUNIT_ASSERT (false);
  }

References std::abs(), libMesh::EquationSystems::add_system(), libMesh::System::add_variable(), libMesh::NumericVector< T >::build(), libMesh::ParallelObject::comm(), libMesh::FIRST, libMesh::System::get_all_variable_numbers(), libMesh::System::get_dof_map(), libMesh::MeshFunction::init(), libMesh::EquationSystems::init(), libMesh::LAGRANGE, libMesh::System::n_dofs(), libMesh::EquationSystems::parameters, position_function2(), libMesh::System::project_solution(), libMesh::SERIAL, and libMesh::System::solution.

test_point_locator_dfem()

void MeshfunctionDFEM::test_point_locator_dfem ( )

inline

Definition at line 133 of file mesh_function_dfem.C.

  {
    // this point is in bottom element only
    Point interior(0.5, -0.5, 0.0);
 
    // this point is on the face between bottom and top
    Point face(0.5, 0.0, 0.0);
 
    // test interior point
    std::set<const Elem *> int_cand;
    (*_point_locator)(interior, int_cand, nullptr);
    CPPUNIT_ASSERT (int_cand.size() == 1);
    for (std::set<const Elem *>::iterator it = int_cand.begin(); it != int_cand.end(); ++it)
      CPPUNIT_ASSERT ((*it)->id() == 1);
 
    // test interior point
    std::set<const Elem *> face_cand;
    (*_point_locator)(face, face_cand, nullptr);
    CPPUNIT_ASSERT (face_cand.size() == 2);
    int array[2] = {0, 0};
    for (std::set<const Elem *>::iterator it = face_cand.begin(); it != face_cand.end(); ++it)
      {
        // first test that array entry hasn't been set before
        CPPUNIT_ASSERT (array[(*it)->id()] == 0);
        array[(*it)->id()] = 1;
      }
    CPPUNIT_ASSERT (array[0] == 1);
    CPPUNIT_ASSERT (array[1] == 1);
  }

Member Data Documentation

_mesh

ReplicatedMesh* MeshfunctionDFEM::_mesh

protected

Definition at line 62 of file mesh_function_dfem.C.

_point_locator

std::unique_ptr<PointLocatorBase> MeshfunctionDFEM::_point_locator

protected

Definition at line 63 of file mesh_function_dfem.C.

---

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

• tests/mesh/mesh_function_dfem.C