Inheritance diagram for PeriodicBCTest:

[legend]

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.

  {
    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.prepare_for_use();
 
    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_NUMBERS_EQUAL(exact, approx,
                                      TOLERANCE*TOLERANCE*20);
        }
  }