MeshAlignment.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

#include "MeshAlignmentBase.h"

class Assembly;

class MeshAlignment : public MeshAlignmentBase
{
public:
  MeshAlignment(const MooseMesh & mesh);

  void initialize(
      const std::vector<dof_id_type> & primary_elem_ids,
      const std::vector<std::tuple<dof_id_type, unsigned short int>> & secondary_boundary_info);

  void initialize(
      const std::vector<std::tuple<dof_id_type, unsigned short int>> & primary_boundary_info,
      const std::vector<std::tuple<dof_id_type, unsigned short int>> & secondary_boundary_info);

  void buildCoupledElemQpIndexMap(Assembly & assembly);

  bool meshesAreCoincident() const { return _meshes_are_coincident; }

  bool hasCoupledElemID(const dof_id_type & elem_id) const;

  const dof_id_type & getCoupledElemID(const dof_id_type & elem_id) const;

  bool hasCoupledNodeID(const dof_id_type & node_id) const;

  const dof_id_type & getCoupledNodeID(const dof_id_type & node_id) const;

  unsigned int getCoupledElemQpIndex(const dof_id_type & elem_id, const unsigned int & qp) const;

protected:
  void buildMapping();

  std::map<dof_id_type, std::vector<Point>>
  getLocalQuadraturePointMap(Assembly & assembly,
                             const std::vector<dof_id_type> & elem_ids,
                             const std::vector<unsigned short int> & side_ids) const;

  std::map<dof_id_type, dof_id_type> _coupled_elem_ids;
  std::map<dof_id_type, dof_id_type> _coupled_node_ids;
  std::map<dof_id_type, std::vector<unsigned int>> _coupled_elem_qp_indices;

  bool _meshes_are_coincident;
};