This class is for unit testing dual coefficient and shape function values. More...

Inheritance diagram for DualShapeTest:

Public Member Functions
	LIBMESH_CPPUNIT_TEST_SUITE (DualShapeTest)

	CPPUNIT_TEST (testEdge2Lagrange)

	CPPUNIT_TEST_SUITE_END ()

void	testEdge2Lagrange ()

void	setUp ()

void	tearDown ()

Private Attributes
std::unique_ptr< Mesh >	_mesh

std::unique_ptr< QGauss >	_qrule

std::unique_ptr< FEBase >	_fe

Elem *	_elem

Detailed Description

This class is for unit testing dual coefficient and shape function values.

Definition at line 19 of file dual_shape_verification_test.C.

Member Function Documentation

◆ CPPUNIT_TEST()

DualShapeTest::CPPUNIT_TEST ( testEdge2Lagrange )

◆ CPPUNIT_TEST_SUITE_END()

DualShapeTest::CPPUNIT_TEST_SUITE_END ( )

◆ LIBMESH_CPPUNIT_TEST_SUITE()

DualShapeTest::LIBMESH_CPPUNIT_TEST_SUITE ( DualShapeTest )

◆ setUp()

void DualShapeTest::setUp ( )

inline

Definition at line 74 of file dual_shape_verification_test.C.

References libMesh::FEGenericBase< OutputType >::build(), libMesh::MeshTools::Generation::build_line(), libMesh::FEType::default_quadrature_order(), libMesh::EDGE2, libMesh::FIRST, and libMesh::LAGRANGE.

     {
       FEType fe_type(FIRST, LAGRANGE);
       _fe = FEBase::build(1, fe_type);
       _fe->get_phi();
       _fe->get_dual_phi();
 
       _mesh = std::make_unique<Mesh>(*TestCommWorld);
 
       MeshTools::Generation::build_line(*_mesh, 1, -1, 1, EDGE2);
 
       auto rng = _mesh->active_local_element_ptr_range();
       _elem = rng.begin() == rng.end() ? nullptr : *(rng.begin());
 
       _qrule = std::make_unique<QGauss>(1, fe_type.default_quadrature_order());
       _fe->attach_quadrature_rule(_qrule.get());
     }

◆ tearDown()

void DualShapeTest::tearDown ( )

inline

Definition at line 92 of file dual_shape_verification_test.C.

92 {}

◆ testEdge2Lagrange()

void DualShapeTest::testEdge2Lagrange ( )

inline

Definition at line 34 of file dual_shape_verification_test.C.

References libMesh::index_range(), libMesh::DenseMatrixBase< T >::m(), libMesh::DenseMatrixBase< T >::n(), libMesh::Real, and libMesh::TOLERANCE.

   {
     LOG_UNIT_TEST;
     if (!_elem)
       return;
 
     _fe->reinit(_elem);
 
     const DenseMatrix<Real> & dual_coeff = _fe->get_dual_coeff();
 
     CPPUNIT_ASSERT_EQUAL(dual_coeff.m(), unsigned(2));
     CPPUNIT_ASSERT_EQUAL(dual_coeff.n(), unsigned(2));
 
     // TOLERANCE*TOLERANCE works with double but not float128
     Real my_tol = TOLERANCE*std::sqrt(TOLERANCE);
 
     LIBMESH_ASSERT_FP_EQUAL(2, dual_coeff(0,0), my_tol);
     LIBMESH_ASSERT_FP_EQUAL(-1, dual_coeff(0,1), my_tol);
     LIBMESH_ASSERT_FP_EQUAL(-1, dual_coeff(1,0), my_tol);
     LIBMESH_ASSERT_FP_EQUAL(2, dual_coeff(1,1), my_tol);
 
     const auto & dual_phi = _fe->get_dual_phi();
 
     CPPUNIT_ASSERT_EQUAL(std::size_t(2), dual_phi.size());
 
     const auto & qpoints = _qrule->get_points();
 
     CPPUNIT_ASSERT_EQUAL(qpoints.size(), dual_phi[0].size());
 
     for (auto qp : index_range(dual_phi[0]))
       LIBMESH_ASSERT_FP_EQUAL(1./2. * (1. - 3.*qpoints[qp](0)), dual_phi[0][qp],
         my_tol);
 
     CPPUNIT_ASSERT_EQUAL(qpoints.size(), dual_phi[1].size());
 
     for (auto qp : index_range(dual_phi[1]))
       LIBMESH_ASSERT_FP_EQUAL(1./2. * (1. + 3.*qpoints[qp](0)), dual_phi[1][qp],
         my_tol);
   }