ExtraNodesetGenerator.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 "ExtraNodesetGenerator.h"
#include "MooseMesh.h"
#include "Conversion.h"
#include "MooseMeshUtils.h"
#include "CastUniquePointer.h"
 
#include "libmesh/elem.h"
 
registerMooseObject("MooseApp", ExtraNodesetGenerator);
 
InputParameters
ExtraNodesetGenerator::validParams()
{
  InputParameters params = MeshGenerator::validParams();
 
  params.addRequiredParam<MeshGeneratorName>("input", "The mesh we want to modify");
  params.addRequiredParam<std::vector<BoundaryName>>("new_boundary",
                                                     "The names of the boundaries to create");
 
  params.addParam<std::vector<unsigned int>>("nodes",
                                             "The nodes you want to be in the nodeset "
                                             "(Either this parameter or \"coord\" must be "
                                             "supplied).");
  params.addParam<std::vector<std::vector<Real>>>(
      "coord",
      {},
      "The nodes with coordinates you want to be in the "
      "nodeset. Separate multple coords with ';' (Either this parameter or \"nodes\" must be "
      "supplied).");
  params.addParam<Real>(
      "tolerance", TOLERANCE, "The tolerance in which two nodes are considered identical");
  params.addParam<bool>("use_closest_node", false, "Use the node closest to the coordinate.");
  params.addClassDescription(
      "Creates a new node set and a new boundary made with the nodes the user provides.");
 
  return params;
}
 
ExtraNodesetGenerator::ExtraNodesetGenerator(const InputParameters & parameters)
  : MeshGenerator(parameters), _input(getMesh("input"))
{
}
 
std::unique_ptr<MeshBase>
ExtraNodesetGenerator::generate()
{
  std::unique_ptr<MeshBase> mesh = std::move(_input);
 
  // Get the BoundaryIDs from the mesh
  std::vector<BoundaryName> boundary_names = getParam<std::vector<BoundaryName>>("new_boundary");
  std::vector<boundary_id_type> boundary_ids =
      MooseMeshUtils::getBoundaryIDs(*mesh, boundary_names, true);
 
  // Get a reference to our BoundaryInfo object
  BoundaryInfo & boundary_info = mesh->get_boundary_info();
 
  // add nodes with their ids
  if (isParamValid("nodes"))
    for (const auto & node_id : getParam<std::vector<unsigned int>>("nodes"))
    {
      // Our mesh may be distributed and this node may not exist on this process
      if (!mesh->query_node_ptr(node_id))
        continue;
 
      for (const auto & boundary_id : boundary_ids)
        boundary_info.add_node(node_id, boundary_id);
    }
 
  // add nodes with their coordinates
  const auto dim = mesh->mesh_dimension();
 
  std::unique_ptr<libMesh::PointLocatorBase> locator = mesh->sub_point_locator();
  locator->enable_out_of_mesh_mode();
 
  const auto tolerance = getParam<Real>("tolerance");
  const bool use_closest_node = getParam<bool>("use_closest_node");
  const auto coords = getParam<std::vector<std::vector<Real>>>("coord");
  if (use_closest_node && coords.empty())
    paramError("coord", "A coordinate should be specified to use 'use_closest_node'");
  for (const auto & c : coords)
  {
    Point p;
    if (c.size() < dim)
      paramError("coord",
                 "Coordinate ",
                 Moose::stringify(c),
                 " does not have enough components for a ",
                 dim,
                 "D mesh.");
 
    if (c.size() > 3)
      paramError("coord",
                 "Coordinate ",
                 Moose::stringify(c),
                 " has too many components. Did you maybe forget to separate multiple coordinates "
                 "with a ';'?");
 
    for (unsigned int j = 0; j < c.size(); ++j)
      p(j) = c[j];
 
    // locate candidate element
    bool on_node = false;
    bool found_elem = false;
    const Elem * elem = (*locator)(p);
    if (elem)
    {
      found_elem = true;
      for (unsigned int j = 0; j < elem->n_nodes(); ++j)
      {
        const Node * node = elem->node_ptr(j);
        if (p.absolute_fuzzy_equals(*node, tolerance))
        {
          for (const auto & boundary_id : boundary_ids)
            boundary_info.add_node(node, boundary_id);
 
          on_node = true;
          break;
        }
      }
    }
 
    // If we are on a distributed mesh, then any particular processor
    // may be unable to find any particular node, but *some* processor
    // should have found it.
    if (!mesh->is_replicated())
    {
      this->comm().max(found_elem);
      this->comm().max(on_node);
    }
 
    // only search for closest node if it is not found on a node
    if (use_closest_node && !on_node)
    {
      // these are always true when using closest node
      found_elem = true;
      on_node = true;
      Real dmin(std::numeric_limits<Real>::max());
      dof_id_type dmin_id(std::numeric_limits<dof_id_type>::max());
      // find this processors closest node and save distance and the node id
      for (const auto & node_iter : as_range(mesh->local_nodes_begin(), mesh->local_nodes_end()))
      {
        Real dist = (*node_iter - p).norm_sq();
        if (dist < dmin)
        {
          dmin = dist;
          dmin_id = node_iter->id();
        }
      }
 
      // get proc id with the global closest node and then communicate that procs closest node_id
      processor_id_type dmin_proc_id;
      this->comm().minloc(dmin, dmin_proc_id);
      this->comm().broadcast(dmin_id, dmin_proc_id);
 
      const Node * node = mesh->query_node_ptr(dmin_id);
      if (node)
        for (const auto & boundary_id : boundary_ids)
          boundary_info.add_node(node, boundary_id);
    }
 
    if (!found_elem)
      mooseError("Unable to locate the following point within the domain, please check its "
                 "coordinates:\n",
                 p);
 
    if (!on_node)
      mooseError("No node found at point:\n",
                 p,
                 "\nSet use_closest_node=true if you want to find the closest node.");
  }
 
  for (unsigned int i = 0; i < boundary_ids.size(); ++i)
    boundary_info.nodeset_name(boundary_ids[i]) = boundary_names[i];
 
  mesh->unset_is_prepared();
  return dynamic_pointer_cast<MeshBase>(mesh);
}