BoundaryDeletionGenerator.C

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
 
#include "BoundaryDeletionGenerator.h"
 
#include "MooseMeshUtils.h"
#include "CastUniquePointer.h"
 
#include "libmesh/boundary_info.h"
 
registerMooseObject("MooseApp", BoundaryDeletionGenerator);
 
InputParameters
BoundaryDeletionGenerator::validParams()
{
  InputParameters params = MeshGenerator::validParams();
 
  params.addClassDescription("Mesh generator which removes side sets");
  params.addRequiredParam<MeshGeneratorName>("input", "The mesh we want to modify");
  params.addRequiredParam<std::vector<BoundaryName>>("boundary_names",
                                                     "The boundaries to be deleted / kept");
  MooseEnum opt("keep remove", "remove");
  params.addParam<MooseEnum>("operation", opt, "Whether to remove or keep the listed boundaries");
 
  return params;
}
 
BoundaryDeletionGenerator::BoundaryDeletionGenerator(const InputParameters & parameters)
  : MeshGenerator(parameters),
    _input(getMesh("input")),
    _boundary_names(getParam<std::vector<BoundaryName>>("boundary_names"))
{
}
 
std::unique_ptr<MeshBase>
BoundaryDeletionGenerator::generate()
{
  std::unique_ptr<MeshBase> mesh = std::move(_input);
 
  // We need prepared boundary id sets here
  if (!mesh->preparation().has_boundary_id_sets)
    mesh->get_boundary_info().regenerate_id_sets();
 
  // Gather the boundary ids of interest from the names
  std::set<boundary_id_type> ids;
  for (const auto & name : _boundary_names)
  {
    auto bid = MooseMeshUtils::getBoundaryID(name, *mesh);
    if (bid == BoundaryInfo::invalid_id)
      paramError("boundary_names", "The boundary '", name, "' was not found in the mesh");
 
    ids.insert(bid);
  }
 
  std::set<boundary_id_type> ids_to_remove;
  if (getParam<MooseEnum>("operation") == "keep")
  {
    ids_to_remove = mesh->get_boundary_info().get_boundary_ids();
    for (const auto & id : ids)
      ids_to_remove.erase(id);
  }
  else
    ids_to_remove = ids;
 
  for (const auto & id : ids_to_remove)
    mesh->get_boundary_info().remove_id(id);
 
  mesh->unset_is_prepared();
  return dynamic_pointer_cast<MeshBase>(mesh);
}