DefaultCouplingTest Class Reference

Inheritance diagram for DefaultCouplingTest:

Public Member Functions
	CPPUNIT_TEST_SUITE (DefaultCouplingTest)
	CPPUNIT_TEST (testCouplingOnEdge3)
	CPPUNIT_TEST (testCouplingOnQuad9)
	CPPUNIT_TEST (testCouplingOnTri6)
	CPPUNIT_TEST (testCouplingOnHex27)
	CPPUNIT_TEST_SUITE_END ()
void	setUp ()
void	tearDown ()
void	testCoupling (const ElemType elem_type)
void	testCouplingOnEdge3 ()
void	testCouplingOnQuad9 ()
void	testCouplingOnTri6 ()
void	testCouplingOnHex27 ()

Detailed Description

Definition at line 31 of file default_coupling_test.C.

Member Function Documentation

CPPUNIT_TEST() [1/4]

DefaultCouplingTest::CPPUNIT_TEST

(

testCouplingOnEdge3

)

CPPUNIT_TEST() [2/4]

DefaultCouplingTest::CPPUNIT_TEST

(

testCouplingOnHex27

)

CPPUNIT_TEST() [3/4]

DefaultCouplingTest::CPPUNIT_TEST

(

testCouplingOnQuad9

)

CPPUNIT_TEST() [4/4]

DefaultCouplingTest::CPPUNIT_TEST

(

testCouplingOnTri6

)

CPPUNIT_TEST_SUITE()

DefaultCouplingTest::CPPUNIT_TEST_SUITE

(

DefaultCouplingTest

)

CPPUNIT_TEST_SUITE_END()

DefaultCouplingTest::CPPUNIT_TEST_SUITE_END

(

)

setUp()

void DefaultCouplingTest::setUp ( )

inline

Definition at line 49 of file default_coupling_test.C.

  {}

tearDown()

void DefaultCouplingTest::tearDown ( )

inline

Definition at line 52 of file default_coupling_test.C.

  {}

testCoupling()

void DefaultCouplingTest::testCoupling ( const ElemType  elem_type )

inline

Definition at line 55 of file default_coupling_test.C.

  {
    Mesh mesh(*TestCommWorld);
 
    EquationSystems es(mesh);
    System &sys = es.add_system<System> ("SimpleSystem");
    sys.add_variable("u", THIRD, HIERARCHIC);
    sys.get_dof_map().default_algebraic_ghosting().set_n_levels(3);
 
    const unsigned int n_elem_per_side = 5;
    const std::unique_ptr<Elem> test_elem = Elem::build(elem_type);
    const unsigned int ymax = test_elem->dim() > 1;
    const unsigned int zmax = test_elem->dim() > 2;
    const unsigned int ny = ymax * n_elem_per_side;
    const unsigned int nz = zmax * n_elem_per_side;
 
    MeshTools::Generation::build_cube (mesh,
                                       n_elem_per_side,
                                       ny,
                                       nz,
                                       0., 1.,
                                       0., ymax,
                                       0., zmax,
                                       elem_type);
 
    es.init();
    sys.project_solution(cubic_default_coupling_test, nullptr, es.parameters);
 
    std::set<dof_id_type> evaluable_elements;
 
    for (const auto & elem : mesh.active_local_element_ptr_range())
      {
        CPPUNIT_ASSERT(sys.get_dof_map().is_evaluable(*elem));
        evaluable_elements.insert(elem->id());
 
        for (unsigned int s1=0; s1 != elem->n_neighbors(); ++s1)
          {
            const Elem * n1 = elem->neighbor_ptr(s1);
            // Let's speed up this test by only checking the ghosted
            // elements which are most likely to break.
            if (!n1 ||
                n1->processor_id() == mesh.processor_id())
              continue;
 
            if (!evaluable_elements.count(n1->id()))
              {
                CPPUNIT_ASSERT(sys.get_dof_map().is_evaluable(*n1));
                evaluable_elements.insert(n1->id());
              }
 
            for (unsigned int s2=0; s2 != elem->n_neighbors(); ++s2)
              {
                const Elem * n2 = n1->neighbor_ptr(s2);
                if (!n2 ||
                    n2->processor_id() == mesh.processor_id())
                  continue;
 
                if (!evaluable_elements.count(n2->id()))
                  {
                    CPPUNIT_ASSERT(sys.get_dof_map().is_evaluable(*n2));
                    evaluable_elements.insert(n2->id());
                  }
 
                for (unsigned int s3=0; s3 != elem->n_neighbors(); ++s3)
                  {
                    const Elem * n3 = n2->neighbor_ptr(s3);
                    if (!n3 ||
                        n3->processor_id() == mesh.processor_id() ||
                        evaluable_elements.count(n3->id()))
                      continue;
 
                    CPPUNIT_ASSERT(sys.get_dof_map().is_evaluable(*n3));
                    evaluable_elements.insert(n3->id());
 
                    Point p = n3->centroid();
 
                    LIBMESH_ASSERT_FP_EQUAL(libmesh_real(sys.point_value(0,p,n3)),
                                            libmesh_real(cubic_default_coupling_test(p,es.parameters,"","")),
                                            TOLERANCE*TOLERANCE);
                  }
              }
          }
      }
 
    const std::size_t n_evaluable =
      std::distance(mesh.evaluable_elements_begin(sys.get_dof_map()),
                    mesh.evaluable_elements_end(sys.get_dof_map()));
 
    CPPUNIT_ASSERT_EQUAL(evaluable_elements.size(), n_evaluable);
  }

References libMesh::MeshBase::active_local_element_ptr_range(), libMesh::EquationSystems::add_system(), libMesh::System::add_variable(), libMesh::Elem::build(), libMesh::MeshTools::Generation::build_cube(), libMesh::Elem::centroid(), cubic_default_coupling_test(), libMesh::DofMap::default_algebraic_ghosting(), distance(), libMesh::MeshBase::evaluable_elements_begin(), libMesh::MeshBase::evaluable_elements_end(), libMesh::System::get_dof_map(), libMesh::HIERARCHIC, libMesh::DofObject::id(), libMesh::EquationSystems::init(), libMesh::DofMap::is_evaluable(), libMesh::libmesh_real(), mesh, libMesh::Elem::neighbor_ptr(), libMesh::EquationSystems::parameters, libMesh::System::point_value(), libMesh::ParallelObject::processor_id(), libMesh::DofObject::processor_id(), libMesh::System::project_solution(), libMesh::DefaultCoupling::set_n_levels(), TestCommWorld, libMesh::THIRD, and libMesh::TOLERANCE.

testCouplingOnEdge3()

void DefaultCouplingTest::testCouplingOnEdge3 ( )

inline

Definition at line 148 of file default_coupling_test.C.

{ testCoupling(EDGE3); }

References libMesh::EDGE3.

testCouplingOnHex27()

void DefaultCouplingTest::testCouplingOnHex27 ( )

inline

Definition at line 151 of file default_coupling_test.C.

{ testCoupling(HEX27); }

References libMesh::HEX27.

testCouplingOnQuad9()

void DefaultCouplingTest::testCouplingOnQuad9 ( )

inline

Definition at line 149 of file default_coupling_test.C.

{ testCoupling(QUAD9); }

References libMesh::QUAD9.

testCouplingOnTri6()

void DefaultCouplingTest::testCouplingOnTri6 ( )

inline

Definition at line 150 of file default_coupling_test.C.

{ testCoupling(TRI6); }

References libMesh::TRI6.

---

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

• tests/base/default_coupling_test.C