GhostLowerDElems.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html
 
// App includes
#include "GhostLowerDElems.h"
#include "Executioner.h"
#include "FEProblemBase.h"
#include "MooseApp.h"
 
// libMesh includes
#include "libmesh/elem.h"
#include "libmesh/mesh_base.h"
#include "libmesh/boundary_info.h"
 
registerMooseObject("MooseApp", GhostLowerDElems);
 
using namespace libMesh;
 
InputParameters
GhostLowerDElems::validParams()
{
  InputParameters params = RelationshipManager::validParams();
  params.set<bool>("attach_geometric_early") = false;
  return params;
}
 
GhostLowerDElems::GhostLowerDElems(const InputParameters & params) : RelationshipManager(params) {}
 
GhostLowerDElems::GhostLowerDElems(const GhostLowerDElems & other) : RelationshipManager(other) {}
 
std::string
GhostLowerDElems::getInfo() const
{
  std::ostringstream oss;
  oss << "GhostLowerDElems";
  return oss.str();
}
 
void
GhostLowerDElems::operator()(const MeshBase::const_element_iterator & range_begin,
                             const MeshBase::const_element_iterator & range_end,
                             const processor_id_type p,
                             map_type & coupled_elements)
{
  mooseAssert(_moose_mesh,
              "The MOOSE mesh must be non-null in order for this relationship manager to work.");
  if (!_moose_mesh->hasLowerD())
    return;
 
  static const CouplingMatrix * const null_mat = nullptr;
 
  for (const Elem * const elem : as_range(range_begin, range_end))
    for (const auto s : elem->side_index_range())
    {
      const Elem * const neighbor = elem->neighbor_ptr(s);
      const bool elem_owns_lowerd = !neighbor || elem->id() < neighbor->id();
 
      const Elem * const lower_d_elem =
          elem_owns_lowerd
              ? _moose_mesh->getLowerDElem(elem, s)
              : _moose_mesh->getLowerDElem(neighbor, neighbor->which_neighbor_am_i(elem));
 
      if (lower_d_elem && lower_d_elem->processor_id() != p)
        coupled_elements.emplace(lower_d_elem, null_mat);
    }
}
 
bool
GhostLowerDElems::operator>=(const RelationshipManager & other) const
{
  return dynamic_cast<const GhostLowerDElems *>(&other);
}
 
std::unique_ptr<GhostingFunctor>
GhostLowerDElems::clone() const
{
  return _app.getFactory().copyConstruct(*this);
}