SideHierarchicNodalSolnTest Class Reference

Inheritance diagram for SideHierarchicNodalSolnTest:

Public Member Functions
	LIBMESH_CPPUNIT_TEST_SUITE (SideHierarchicNodalSolnTest)
	CPPUNIT_TEST (testConstantQuad9)
	CPPUNIT_TEST_SUITE_END ()
void	testConstantQuad9 ()

Detailed Description

Definition at line 539 of file fe_side_test.C.

Member Function Documentation

CPPUNIT_TEST()

SideHierarchicNodalSolnTest::CPPUNIT_TEST

(

testConstantQuad9

)

CPPUNIT_TEST_SUITE_END()

SideHierarchicNodalSolnTest::CPPUNIT_TEST_SUITE_END

(

)

LIBMESH_CPPUNIT_TEST_SUITE()

SideHierarchicNodalSolnTest::LIBMESH_CPPUNIT_TEST_SUITE

(

SideHierarchicNodalSolnTest

)

testConstantQuad9()

void SideHierarchicNodalSolnTest::testConstantQuad9 ( )

inline

Definition at line 546 of file fe_side_test.C.

  {
    LOG_UNIT_TEST;
 
#if LIBMESH_DIM > 1
    Mesh mesh(*TestCommWorld);
    MeshTools::Generation::build_square(mesh, 1, 1, 0., 1., 0., 1., QUAD9);
 
    auto range = mesh.active_local_element_ptr_range();
    if (range.begin() == range.end())
      return;
 
    const Elem * elem = *range.begin();
    CPPUNIT_ASSERT_EQUAL(static_cast<int>(QUAD9),
                         static_cast<int>(elem->type()));
 
    const FEType fe_type(CONSTANT, SIDE_HIERARCHIC);
    const std::vector<Number> elem_soln = {1, 3, 5, 7};
 
    std::vector<Number> nodal_soln;
    FEInterface::nodal_soln(/*dim=*/2,
                            fe_type,
                            elem,
                            elem_soln,
                            nodal_soln,
                            /*add_p_level=*/false);
 
    CPPUNIT_ASSERT_EQUAL(std::size_t(9), nodal_soln.size());
 
    std::vector<Number> expected(nodal_soln.size(), 0);
    std::vector<unsigned int> expected_count(nodal_soln.size(), 0);
    for (const auto side : elem->side_index_range())
      for (const auto n : elem->nodes_on_side(side))
        {
          expected[n] += elem_soln[side];
          ++expected_count[n];
        }
 
    unsigned int n_nodes_without_side_value = 0;
    for (const auto n : index_range(expected))
      {
        if (expected_count[n])
          expected[n] /= expected_count[n];
        else
          ++n_nodes_without_side_value;
 
        LIBMESH_ASSERT_NUMBERS_EQUAL(expected[n],
                                     nodal_soln[n],
                                     TOLERANCE * TOLERANCE);
      }
 
    CPPUNIT_ASSERT_EQUAL(1u, n_nodes_without_side_value);
#endif
  }

References mesh, and TestCommWorld.

---

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

• tests/fe/fe_side_test.C