GeometricSearchData.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* 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
 
// MOOSE includes
#include "MooseTypes.h"
 
// libMesh includes
#include "libmesh/enum_order.h"
 
// C++ includes
#include <map>
 
// Forward Declarations
class MooseMesh;
class SubProblem;
class PenetrationLocator;
class NearestNodeLocator;
class ElementPairLocator;
 
class GeometricSearchData
{
public:
  enum GeometricSearchType
  {
    ALL,
    NEAREST_NODE,
    PENETRATION,
    QUADRATURE,
    ELEMENTPAIR
  };
 
  GeometricSearchData(SubProblem & subproblem, MooseMesh & mesh);
  virtual ~GeometricSearchData();
 
  PenetrationLocator & getPenetrationLocator(const BoundaryName & master,
                                             const BoundaryName & slave,
                                             Order order = FIRST);
  PenetrationLocator & getQuadraturePenetrationLocator(const BoundaryName & master,
                                                       const BoundaryName & slave,
                                                       Order order = FIRST);
 
  NearestNodeLocator & getNearestNodeLocator(const BoundaryName & master,
                                             const BoundaryName & slave);
  NearestNodeLocator & getNearestNodeLocator(const unsigned int master_id,
                                             const unsigned int slave_id);
 
  NearestNodeLocator & getQuadratureNearestNodeLocator(const BoundaryName & master,
                                                       const BoundaryName & slave);
  NearestNodeLocator & getQuadratureNearestNodeLocator(const unsigned int master_id,
                                                       const unsigned int slave_id);
 
  const std::map<std::pair<unsigned int, unsigned int>, PenetrationLocator *> &
  getPenetrationLocators() const
  {
    return _penetration_locators;
  }
 
  void addElementPairLocator(const unsigned int & interface_id,
                             std::shared_ptr<ElementPairLocator> epl);
 
  void update(GeometricSearchType type = ALL);
 
  void reinit();
 
  void clearNearestNodeLocators();
 
  Real maxPatchPercentage();
 
  void updateGhostedElems();
 
  // protected:
  SubProblem & _subproblem;
  MooseMesh & _mesh;
  std::map<std::pair<unsigned int, unsigned int>, PenetrationLocator *> _penetration_locators;
  std::map<std::pair<unsigned int, unsigned int>, NearestNodeLocator *> _nearest_node_locators;
  std::map<unsigned int, std::shared_ptr<ElementPairLocator>> _element_pair_locators;
 
protected:
  std::set<unsigned int> _quadrature_boundaries;
 
  std::map<unsigned int, unsigned int> _slave_to_qslave;
 
private:
  void
  generateQuadratureNodes(unsigned int slave_id, unsigned int qslave_id, bool reiniting = false);
 
  void updateQuadratureNodes(unsigned int slave_id);
 
  void reinitQuadratureNodes(unsigned int slave_id);
 
  bool _first;
};