Line data    Source code

       1             : //* This file is part of the MOOSE framework
       2             : //* https://mooseframework.inl.gov
       3             : //*
       4             : //* All rights reserved, see COPYRIGHT for full restrictions
       5             : //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
       6             : //*
       7             : //* Licensed under LGPL 2.1, please see LICENSE for details
       8             : //* https://www.gnu.org/licenses/lgpl-2.1.html
       9             : 
      10             : #include "LowerDBlockFromSidesetGenerator.h"
      11             : #include "InputParameters.h"
      12             : #include "MooseTypes.h"
      13             : #include "CastUniquePointer.h"
      14             : #include "MooseMeshUtils.h"
      15             : 
      16             : #include "libmesh/distributed_mesh.h"
      17             : #include "libmesh/elem.h"
      18             : #include "libmesh/parallel_elem.h"
      19             : #include "libmesh/parallel_node.h"
      20             : #include "libmesh/compare_elems_by_level.h"
      21             : #include "libmesh/mesh_communication.h"
      22             : 
      23             : #include "timpi/parallel_sync.h"
      24             : 
      25             : #include <set>
      26             : #include <typeinfo>
      27             : 
      28             : registerMooseObject("MooseApp", LowerDBlockFromSidesetGenerator);
      29             : 
      30             : InputParameters
      31       19385 : LowerDBlockFromSidesetGenerator::validParams()
      32             : {
      33       19385 :   InputParameters params = MeshGenerator::validParams();
      34             : 
      35       77540 :   params.addRequiredParam<MeshGeneratorName>("input", "The mesh we want to modify");
      36       77540 :   params.addParam<SubdomainID>("new_block_id", "The lower dimensional block id to create");
      37       77540 :   params.addParam<SubdomainName>("new_block_name",
      38             :                                  "The lower dimensional block name to create (optional)");
      39       77540 :   params.addRequiredParam<std::vector<BoundaryName>>(
      40             :       "sidesets", "The sidesets from which to create the new block");
      41             : 
      42       19385 :   params.addClassDescription("Adds lower dimensional elements on the specified sidesets.");
      43             : 
      44       19385 :   return params;
      45           0 : }
      46             : 
      47        2320 : LowerDBlockFromSidesetGenerator::LowerDBlockFromSidesetGenerator(const InputParameters & parameters)
      48             :   : MeshGenerator(parameters),
      49        2320 :     _input(getMesh("input")),
      50        6960 :     _sideset_names(getParam<std::vector<BoundaryName>>("sidesets"))
      51             : {
      52        2320 : }
      53             : 
      54             : std::unique_ptr<MeshBase>
      55        2183 : LowerDBlockFromSidesetGenerator::generate()
      56             : {
      57        2183 :   std::unique_ptr<MeshBase> mesh = std::move(_input);
      58             : 
      59             :   // Generate a new block id if one isn't supplied.
      60        4366 :   SubdomainID new_block_id = isParamValid("new_block_id")
      61        8342 :                                  ? getParam<SubdomainID>("new_block_id")
      62        2183 :                                  : MooseMeshUtils::getNextFreeSubdomainID(*mesh);
      63             :   try
      64             :   {
      65        4378 :     MooseMeshUtils::createSubdomainFromSidesets(mesh,
      66        2183 :                                                 _sideset_names,
      67             :                                                 new_block_id,
      68        4374 :                                                 isParamValid("new_block_name")
      69        9523 :                                                     ? getParam<SubdomainName>("new_block_name")
      70             :                                                     : SubdomainName(),
      71        2183 :                                                 type());
      72             :   }
      73           4 :   catch (MooseException & e)
      74             :   {
      75           8 :     if (((std::string)e.what()).compare(0, 12, "The sideset ") == 0)
      76           8 :       paramError("sidesets", e.what());
      77             :     else
      78           0 :       mooseError(e.what());
      79           0 :   }
      80             : 
      81        2175 :   return mesh;
      82           0 : }