BreakMeshByBlockGenerator.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 "MeshGenerator.h"
#include <unordered_set>
 
/*
 * A mesh generator to split a mesh by a set of blocks
 */
class BreakMeshByBlockGenerator : public MeshGenerator
{
public:
  static InputParameters validParams();
 
  BreakMeshByBlockGenerator(const InputParameters & parameters);
 
  std::unique_ptr<MeshBase> generate() override;
 
private:
  void findBoundaryName(const MeshBase & mesh,
                        subdomain_id_type primaryBlockID,
                        subdomain_id_type secondaryBlockID,
                        BoundaryName & boundary_name,
                        boundary_id_type boundaryID,
                        BoundaryInfo & boundary_info);
 
  std::set<std::pair<std::string, BoundaryID>> _bName_bID_set;
 
  BoundaryName generateBoundaryName(const MeshBase & mesh,
                                    subdomain_id_type primaryBlockID,
                                    subdomain_id_type secondaryBlockID);
 
  void mapBoundaryIdAndBoundaryName(boundary_id_type boundaryID, const std::string & boundaryName);
 
  void addInterface(MeshBase & mesh);
 
  std::unordered_map<dof_id_type, std::set<subdomain_id_type>> syncConnectedBlocks(
      const std::unordered_map<dof_id_type, std::vector<dof_id_type>> & node_to_elem_map,
      const MeshBase & mesh);
 
  subdomain_id_type blockRestrictedElementSubdomainID(const Elem * elem);
 
  bool findBlockPairs(subdomain_id_type block_one, subdomain_id_type block_two);
 
  std::unique_ptr<MeshBase> & _input;
  std::unordered_set<std::pair<SubdomainID, SubdomainID>> _block_pairs;
  std::unordered_set<SubdomainID> _block_set;
  const bool _block_pairs_restricted;
  const bool _surrounding_blocks_restricted;
 
  using NodeConnectedBlocksPair = std::pair<dof_id_type, std::vector<subdomain_id_type>>;
  using SubdomainPair = std::pair<subdomain_id_type, subdomain_id_type>;
  using ElemSide = std::tuple<const Elem *, unsigned int>;
 
  std::set<SubdomainPair> _neighboring_block_list;
 
  std::unordered_map<SubdomainPair, std::set<ElemSide>> _subid_pairs_to_sides;
 
  std::unordered_map<std::pair<subdomain_id_type, subdomain_id_type>, boundary_id_type>
      _subid_pairs_to_boundary_id;
};