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.

  {
    _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();
  }

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

Referenced by NonManifoldGhostingFunctorTest::verify_send_list_entries_helper().

◆ 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;
  }