BreakMeshByElementGenerator.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 "libmesh/elem.h"
 
/*
 * A mesh generator to split a mesh by breaking all element-element interfaces in the
 * specified subdomains
 */
class BreakMeshByElementGenerator : public MeshGenerator
{
public:
  static InputParameters validParams();
 
  BreakMeshByElementGenerator(const InputParameters & parameters);
 
  std::unique_ptr<MeshBase> generate() override;
 
protected:
  // It is important to maintain the exact same ordering across processors. Otherwise nodes with the
  // same ID might have different processor IDs. std::map and std::set of dof_id_type are used here
  // to ensure the same ordering.
  typedef std::map<const dof_id_type, std::set<dof_id_type>> NodeToElemMapType;
 
  NodeToElemMapType
  buildSubdomainRestrictedNodeToElemMap(std::unique_ptr<MeshBase> & mesh,
                                        const std::vector<SubdomainID> & subdomains) const;
 
  void duplicateNodes(std::unique_ptr<MeshBase> & mesh,
                      const NodeToElemMapType & node_to_elem_map) const;
 
  void duplicateNode(std::unique_ptr<MeshBase> & mesh, Elem * elem, const Node * node) const;
 
  void createInterface(MeshBase & mesh, const NodeToElemMapType & node_to_elem_map) const;
 
  std::unique_ptr<MeshBase> & _input;
 
  // The subdomains to explode
  const std::vector<SubdomainID> & _subdomains;
 
  // The name of the new boundary
  const BoundaryName _interface_name;
 
  // Whether to put that boundary on 0, 1, or 2 sides of each new
  // element break
  const unsigned int _interface_sides;
};