PeriodicBCTest Class Reference

Inheritance diagram for PeriodicBCTest:

Public Member Functions
	LIBMESH_CPPUNIT_TEST_SUITE (PeriodicBCTest)
	CPPUNIT_TEST (testPeriodicLagrange2)
	CPPUNIT_TEST_SUITE_END ()

Private Member Functions
void	testPeriodicBC (FEType fe_type)
void	testPeriodicLagrange2 ()

Detailed Description

Definition at line 161 of file periodic_bc_test.C.

Member Function Documentation

CPPUNIT_TEST()

PeriodicBCTest::CPPUNIT_TEST

(

testPeriodicLagrange2

)

CPPUNIT_TEST_SUITE_END()

PeriodicBCTest::CPPUNIT_TEST_SUITE_END

(

)

LIBMESH_CPPUNIT_TEST_SUITE()

PeriodicBCTest::LIBMESH_CPPUNIT_TEST_SUITE

(

PeriodicBCTest

)

testPeriodicBC()

void PeriodicBCTest::testPeriodicBC ( FEType  fe_type )

inlineprivate

Definition at line 175 of file periodic_bc_test.C.

References libMesh::DofMap::add_dirichlet_boundary(), libMesh::DofMap::add_periodic_boundary(), libMesh::BoundaryInfo::add_side(), libMesh::EquationSystems::add_system(), libMesh::System::add_variable(), libMesh::MeshBase::allow_remote_element_removal(), libMesh::System::attach_assemble_function(), bottom_id, libMesh::MeshBase::delete_remote_elements(), libMesh::MeshBase::get_boundary_info(), libMesh::System::get_dof_map(), libMesh::EquationSystems::init(), libMesh::libmesh_real(), mesh, libMesh::PeriodicBoundaryBase::myboundary, libMesh::PeriodicBoundaryBase::pairedboundary, periodic_bc_test_poisson(), libMesh::System::point_value(), quadratic_solution(), libMesh::MeshBase::read(), libMesh::Real, side_id, libMesh::LinearImplicitSystem::solve(), TestCommWorld, libMesh::TOLERANCE, top_id, and libMesh::ExodusII_IO::write_equation_systems().

  {
    Mesh mesh(*TestCommWorld);
    BoundaryInfo & boundary = mesh.get_boundary_info();

    EquationSystems es(mesh);
    LinearImplicitSystem &sys =
      es.add_system<LinearImplicitSystem> ("PBCSys");

    // We need to change around BCIDs and set up the periodic BC first
    mesh.allow_remote_element_removal(false);

    mesh.read("meshes/shark_tooth_tri6.xda.gz");

    // Add some BCIDs on top and bottom for our periodic boundary
    for (const Elem * elem : mesh.active_element_ptr_range())
      {
        for (unsigned int s=0; s != 3; ++s)
          if (!elem->neighbor_ptr(s))
            {
              unsigned v2 = (s+1)%3;
              // Left-running edges are on bottom, because this xda
              // has a flipped mesh
              if (elem->point(s)(0) - elem->point(v2)(0) < -.1)
                boundary.add_side(elem, s, top_id);

              // Right-running edges are on top
              else if (elem->point(s)(0) - elem->point(v2)(0) > .1)
                boundary.add_side(elem, s, bottom_id);

              // Vertical edges are Dirichlet
              else
                boundary.add_side(elem, s, side_id);
            }
      }

    PeriodicBoundary vert(RealVectorValue(0., 4.0));
    vert.myboundary = bottom_id;
    vert.pairedboundary = top_id;
    sys.get_dof_map().add_periodic_boundary(vert);

    // Okay, *now* the mesh knows to save periodic neighbors
    mesh.allow_remote_element_removal(true);
    mesh.delete_remote_elements();

    const unsigned int u_var = sys.add_variable("u", fe_type);
    sys.attach_assemble_function (periodic_bc_test_poisson);

    std::set<boundary_id_type> side_bdy { side_id };
    std::vector<unsigned int> all_vars (1, u_var);
    WrappedFunction<Number> exact_val(sys, quadratic_solution);
    DirichletBoundary exact_bc(side_bdy, all_vars, exact_val);
    sys.get_dof_map().add_dirichlet_boundary(exact_bc);

    es.init();

    sys.solve();

    ExodusII_IO(mesh).write_equation_systems("periodic.e", es);

    Parameters params;

    for (Real x = 0.1; x < 10; x += 0.2)
      for (Real ya = -2.9; ya < 1; ya += 0.2)
        {
          // Integrate the sharktooth shape
          const Real offset = 1-std::abs(std::fmod(x,2)-1);
          const Real y = ya + offset;
          Point p{x,y};
          const Number exact = quadratic_solution(p,params,"","");
          const Number approx = sys.point_value(0,p);
          LIBMESH_ASSERT_FP_EQUAL(libmesh_real(exact),
                                  libmesh_real(approx),
                                  TOLERANCE*TOLERANCE*20);
        }
  }

testPeriodicLagrange2()

void PeriodicBCTest::testPeriodicLagrange2 ( )

inlineprivate

Definition at line 253 of file periodic_bc_test.C.

References libMesh::LAGRANGE, and libMesh::SECOND.

{ LOG_UNIT_TEST; testPeriodicBC(FEType(SECOND, LAGRANGE)); }

---

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

• tests/systems/periodic_bc_test.C