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 "BreakMeshByBlockGeneratorBase.h"
#include <unordered_set>

/*
 * A mesh generator to split a mesh by a set of blocks
 */
class BreakMeshByBlockGenerator : public BreakMeshByBlockGeneratorBase
{
public:
  static InputParameters validParams();

  BreakMeshByBlockGenerator(const InputParameters & parameters);

  std::unique_ptr<MeshBase> generate() override;

protected:
  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;
  const bool _add_transition_interface;
  const bool _split_transition_interface;
  const BoundaryName _interface_transition_name;
  const bool _add_interface_on_two_sides;

private:
  void addInterfaceBoundary(MeshBase & mesh);

  std::set<std::pair<subdomain_id_type, subdomain_id_type>> _neighboring_block_list;
  std::map<std::pair<subdomain_id_type, subdomain_id_type>,
           std::set<std::pair<dof_id_type, unsigned int>>>
      _new_boundary_sides_map;
};