AugmentSparsityOnInterface.h

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

#pragma once

// App includes
#include "AutomaticMortarGeneration.h"
#include "RelationshipManager.h"

// libMesh includes
#include "libmesh/mesh_base.h"

using libMesh::boundary_id_type;
using libMesh::CouplingMatrix;
using libMesh::Elem;
using libMesh::GhostingFunctor;
using libMesh::MeshBase;
using libMesh::processor_id_type;

class AugmentSparsityOnInterface : public RelationshipManager
{
public:
  AugmentSparsityOnInterface(const InputParameters &);

  AugmentSparsityOnInterface(const AugmentSparsityOnInterface & others);

  static InputParameters validParams();

  virtual void operator()(const MeshBase::const_element_iterator & range_begin,
                          const MeshBase::const_element_iterator & range_end,
                          processor_id_type p,
                          map_type & coupled_elements) override;

  virtual std::unique_ptr<GhostingFunctor> clone() const override;

  virtual void redistribute() override { this->mesh_reinit(); }

  std::string getInfo() const override;

  virtual bool operator>=(const RelationshipManager & other) const override;

protected:
  virtual void internalInitWithMesh(const MeshBase &) override;

  void ghostMortarInterfaceCouplings(const processor_id_type p,
                                     const Elem * const elem,
                                     map_type & coupled_elements,
                                     const AutomaticMortarGeneration & amg) const;

  void ghostLowerDSecondaryElemPointNeighbors(const processor_id_type p,
                                              const Elem * const query_elem,
                                              map_type & coupled_elements,
                                              BoundaryID secondary_boundary_id,
                                              SubdomainID secondary_subdomain_id,
                                              const AutomaticMortarGeneration & amg) const;

  void ghostHigherDNeighbors(const processor_id_type p,
                             const Elem * const query_elem,
                             map_type & coupled_elements,
                             BoundaryID secondary_boundary_id,
                             SubdomainID secondary_subdomain_id,
                             const AutomaticMortarGeneration & amg) const;

  const BoundaryName _primary_boundary_name;
  const BoundaryName _secondary_boundary_name;
  const SubdomainName _primary_subdomain_name;
  const SubdomainName _secondary_subdomain_name;

  const bool _is_coupling_functor;

  const bool _ghost_point_neighbors;

  const bool _ghost_higher_d_neighbors;

  const CouplingMatrix * const _null_mat = nullptr;
};