Inheritance diagram for NonManifoldCouplingTestBase:

[legend]

Protected Member Functions
void	read_mesh (const std::string &mesh_filename)

void	init_es ()

Protected Attributes
std::unique_ptr< MeshBase >	_mesh

std::unique_ptr< EquationSystems >	_es

Detailed Description

Definition at line 115 of file nonmanifold_coupling_test.C.

Member Function Documentation

◆ init_es()

void NonManifoldCouplingTestBase::init_es ( )

inlineprotected

Definition at line 142 of file nonmanifold_coupling_test.C.

References libMesh::DofMap::add_coupling_functor(), libMesh::System::add_variable(), libMesh::FIRST, libMesh::System::get_dof_map(), and libMesh::LAGRANGE.

   {
     _es = std::make_unique<EquationSystems>(*_mesh);
 
     // Add System with a single linear Lagrange variable on it.
     LinearImplicitSystem & sys = _es->add_system<LinearImplicitSystem> ("SimpleSystem");
     sys.add_variable("u", FIRST, LAGRANGE);
 
     // Attach ghosting functor to the System. We use the std::shared_ptr interface
     // so that the DofMap takes ownership of the object. Note: if you comment out
     // this line, then the test _can_ fail when run in parallel, but it depends on
     // the number of processors used. For small numbers of procs, sometimes we can
     // get "lucky" and wind up with all the non-manifold neighbor DOFs either local
     // to or ghosted on all procs where they are needed. In practice, this test should
     // be run on 3 or more procs to make it fail when the NonManifoldGhostingFunctor
     // is not used.
     auto ghosting_functor = std::make_shared<NonManifoldGhostingFunctor>(*_mesh);
     sys.get_dof_map().add_coupling_functor(ghosting_functor, /*to_mesh=*/true);
 
     // Initialize the DofMap, etc. for all Systems
     _es->init();
   }

◆ read_mesh()

void NonManifoldCouplingTestBase::read_mesh ( const std::string & mesh_filename )

inlineprotected

Definition at line 123 of file nonmanifold_coupling_test.C.

   {
     // We are making assumptions in various places about the presence
     // of the elements on the current processor so we're restricting to
     // ReplicatedMesh for now.
     _mesh = std::make_unique<ReplicatedMesh>(*TestCommWorld);
     _mesh->read(mesh_filename);
 
     // Use our custom partitioner that assigns non-manifold edge
     // neighbors to different processors.
     _mesh->partitioner() = std::make_unique<NonManifoldTestPartitioner>();
 
     _mesh->prepare_for_use();
 
     // Debugging: print processor ids determined by partitioner
     // for (const auto & elem : _mesh->element_ptr_range())
     //   libMesh::out << "Elem " << elem->id() << " is on processor " << elem->processor_id() << std::endl;
   }