MeshfunctionDFEM Class Reference

Inheritance diagram for MeshfunctionDFEM:

Public Member Functions
	LIBMESH_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
std::unique_ptr< 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.

References libMesh::Elem::build(), libMesh::QUAD4, and libMesh::Elem::set_node().

  {
    _mesh = std::make_unique<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(Elem::build(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(Elem::build(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));
    }

    _mesh->allow_renumbering(true);
    _mesh->prepare_for_use();

    // get a point locator
    _point_locator = _mesh->sub_point_locator();
  }

CPPUNIT_TEST() [1/3]

MeshfunctionDFEM::CPPUNIT_TEST

(

test_point_locator_dfem

)

CPPUNIT_TEST() [2/3]

MeshfunctionDFEM::CPPUNIT_TEST

(

test_mesh_function_dfem

)

CPPUNIT_TEST() [3/3]

MeshfunctionDFEM::CPPUNIT_TEST

(

test_mesh_function_dfem_grad

)

CPPUNIT_TEST_SUITE_END()

MeshfunctionDFEM::CPPUNIT_TEST_SUITE_END

(

)

LIBMESH_CPPUNIT_TEST_SUITE()

MeshfunctionDFEM::LIBMESH_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.

setUp()

void MeshfunctionDFEM::setUp ( )

inline

Definition at line 116 of file mesh_function_dfem.C.

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

tearDown()

void MeshfunctionDFEM::tearDown ( )

inline

Definition at line 123 of file mesh_function_dfem.C.

{}

test_mesh_function_dfem()

void MeshfunctionDFEM::test_mesh_function_dfem ( )

inline

Definition at line 159 of file mesh_function_dfem.C.

References 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.

  {
    LOG_UNIT_TEST;

    // 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);

  }

test_mesh_function_dfem_grad()

void MeshfunctionDFEM::test_mesh_function_dfem_grad ( )

inline

Definition at line 229 of file mesh_function_dfem.C.

References 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.

  {
    LOG_UNIT_TEST;

    // 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);
  }

test_point_locator_dfem()

void MeshfunctionDFEM::test_point_locator_dfem ( )

inline

Definition at line 126 of file mesh_function_dfem.C.

  {
    LOG_UNIT_TEST;

    // 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

std::unique_ptr<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