#include <CoarseMeshExtraElementIDGenerator.h>

Inheritance diagram for CoarseMeshExtraElementIDGenerator:

Public Types
typedef DataFileName	DataFileParameterType

Public Member Functions
	CoarseMeshExtraElementIDGenerator (const InputParameters &parameters)

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

std::unique_ptr< MeshBase >	generateInternal ()

const std::set< MeshGeneratorName > &	getRequestedMeshGenerators () const

const std::set< MeshGeneratorName > &	getRequestedMeshGeneratorsForSub () const

void	addParentMeshGenerator (const MeshGenerator &mg, const AddParentChildKey)

void	addChildMeshGenerator (const MeshGenerator &mg, const AddParentChildKey)

const std::set< const MeshGenerator *, Comparator > &	getParentMeshGenerators () const

const std::set< const MeshGenerator *, Comparator > &	getChildMeshGenerators () const

const std::set< const MeshGenerator *, Comparator > &	getSubMeshGenerators () const

bool	isParentMeshGenerator (const MeshGeneratorName &name, const bool direct=true) const

bool	isChildMeshGenerator (const MeshGeneratorName &name, const bool direct=true) const

bool	isNullMeshName (const MeshGeneratorName &name) const

bool	hasSaveMesh () const

bool	hasOutput () const

const std::string &	getSavedMeshName () const

bool	hasGenerateData () const

bool	isDataOnly () const

virtual bool	enabled () const

std::shared_ptr< MooseObject >	getSharedPtr ()

std::shared_ptr< const MooseObject >	getSharedPtr () const

MooseApp &	getMooseApp () const

const std::string &	type () const

virtual const std::string &	name () const

std::string	typeAndName () const

std::string	errorPrefix (const std::string &error_type) const

void	callMooseError (std::string msg, const bool with_prefix) const

MooseObjectParameterName	uniqueParameterName (const std::string &parameter_name) const

const InputParameters &	parameters () const

MooseObjectName	uniqueName () const

const T &	getParam (const std::string &name) const

std::vector< std::pair< T1, T2 > >	getParam (const std::string &param1, const std::string &param2) const

const T *	queryParam (const std::string &name) const

const T &	getRenamedParam (const std::string &old_name, const std::string &new_name) const

T	getCheckedPointerParam (const std::string &name, const std::string &error_string="") const

bool	isParamValid (const std::string &name) const

bool	isParamSetByUser (const std::string &nm) const

void	paramError (const std::string &param, Args... args) const

void	paramWarning (const std::string &param, Args... args) const

void	paramInfo (const std::string &param, Args... args) const

void	connectControllableParams (const std::string &parameter, const std::string &object_type, const std::string &object_name, const std::string &object_parameter) const

void	mooseError (Args &&... args) const

void	mooseErrorNonPrefixed (Args &&... args) const

void	mooseDocumentedError (const std::string &repo_name, const unsigned int issue_num, Args &&... args) const

void	mooseWarning (Args &&... args) const

void	mooseWarningNonPrefixed (Args &&... args) const

void	mooseDeprecated (Args &&... args) const

void	mooseInfo (Args &&... args) const

std::string	getDataFileName (const std::string &param) const

std::string	getDataFileNameByName (const std::string &relative_path) const

std::string	getDataFilePath (const std::string &relative_path) const

const Parallel::Communicator &	comm () const

processor_id_type	n_processors () const

processor_id_type	processor_id () const

Static Public Member Functions
static InputParameters	validParams ()

static bool	hasGenerateData (const InputParameters &params)

static void	setHasGenerateData (InputParameters &params)

Public Attributes
const ConsoleStream	_console

Static Public Attributes
static const std::string	data_only_param

static constexpr auto	SYSTEM

static constexpr auto	NAME

Protected Member Functions
virtual void	generateData ()

T &	copyMeshProperty (const std::string &target_data_name, const std::string &source_data_name, const std::string &source_mesh)

T &	copyMeshProperty (const std::string &source_data_name, const std::string &source_mesh)

std::unique_ptr< MeshBase > &	getMesh (const std::string &param_name, const bool allow_invalid=false)

std::vector< std::unique_ptr< MeshBase > *>	getMeshes (const std::string &param_name)

std::unique_ptr< MeshBase > &	getMeshByName (const MeshGeneratorName &mesh_generator_name)

std::vector< std::unique_ptr< MeshBase > *>	getMeshesByName (const std::vector< MeshGeneratorName > &mesh_generator_names)

void	declareMeshForSub (const std::string &param_name)

void	declareMeshesForSub (const std::string &param_name)

void	declareMeshForSubByName (const MeshGeneratorName &mesh_generator_name)

void	declareMeshesForSubByName (const std::vector< MeshGeneratorName > &mesh_generator_names)

std::unique_ptr< MeshBase >	buildMeshBaseObject (unsigned int dim=libMesh::invalid_uint)

std::unique_ptr< ReplicatedMesh >	buildReplicatedMesh (unsigned int dim=libMesh::invalid_uint)

std::unique_ptr< DistributedMesh >	buildDistributedMesh (unsigned int dim=libMesh::invalid_uint)

void	addMeshSubgenerator (const std::string &type, const std::string &name, Ts... extra_input_parameters)

void	addMeshSubgenerator (const std::string &type, const std::string &name, InputParameters params)

void	declareNullMeshName (const MeshGeneratorName &name)

const T &	getMeshProperty (const std::string &data_name, const std::string &prefix)

const T &	getMeshProperty (const std::string &data_name)

bool	hasMeshProperty (const std::string &data_name, const std::string &prefix) const

bool	hasMeshProperty (const std::string &data_name, const std::string &prefix) const

bool	hasMeshProperty (const std::string &data_name) const

bool	hasMeshProperty (const std::string &data_name) const

std::string	meshPropertyName (const std::string &data_name) const

T &	declareMeshProperty (const std::string &data_name, Args &&... args)

T &	declareMeshProperty (const std::string &data_name, const T &data_value)

T &	declareMeshProperty (const std::string &data_name, Args &&... args)

T &	declareMeshProperty (const std::string &data_name, const T &data_value)

T &	setMeshProperty (const std::string &data_name, Args &&... args)

T &	setMeshProperty (const std::string &data_name, const T &data_value)

T &	setMeshProperty (const std::string &data_name, Args &&... args)

T &	setMeshProperty (const std::string &data_name, const T &data_value)

Static Protected Member Functions
static std::string	meshPropertyName (const std::string &data_name, const std::string &prefix)

Protected Attributes
std::unique_ptr< MeshBase > &	_input
	input mesh for adding element IDs More...

std::unique_ptr< MeshBase > &	_coarse_mesh
	name of the coarse mesh file More...

const std::string &	_coarse_id_name
	coarse element ID name More...

const bool	_using_coarse_element_id
	whether or not using the coarse element ID for the extra element ID assignment More...

const bool &	_embedding_necessary
	whether the input mesh must be embedded in the coarse mesh More...

MooseMesh *const	_mesh

const bool &	_enabled

MooseApp &	_app

const std::string	_type

const std::string	_name

const InputParameters &	_pars

Factory &	_factory

ActionFactory &	_action_factory

const Parallel::Communicator &	_communicator

Detailed Description

Definition at line 14 of file CoarseMeshExtraElementIDGenerator.h.

Constructor & Destructor Documentation

◆ CoarseMeshExtraElementIDGenerator()

CoarseMeshExtraElementIDGenerator::CoarseMeshExtraElementIDGenerator ( const InputParameters & parameters )

Definition at line 45 of file CoarseMeshExtraElementIDGenerator.C.

   : MeshGenerator(params),
     _input(getMesh("input")),
     _coarse_mesh(getMesh("coarse_mesh")),
     _coarse_id_name(getParam<std::string>("extra_element_id_name")),
     _using_coarse_element_id(isParamValid("coarse_mesh_extra_element_id")),
     _embedding_necessary(getParam<bool>("enforce_mesh_embedding"))
 {
 }

Member Function Documentation

◆ generate()

std::unique_ptr< MeshBase > CoarseMeshExtraElementIDGenerator::generate ( )

overridevirtual

Implements MeshGenerator.

Definition at line 56 of file CoarseMeshExtraElementIDGenerator.C.

 {
   std::unique_ptr<MeshBase> mesh = std::move(_input);
 
   unsigned int coarse_id;
   const bool already_has_id = mesh->has_elem_integer(_coarse_id_name);
   if (!already_has_id)
     coarse_id = mesh->add_elem_integer(_coarse_id_name);
   else
     coarse_id = mesh->get_elem_integer_index(_coarse_id_name);
 
   // Get the requested subdomain IDs
   std::set<SubdomainID> included_subdomains;
   std::set<SubdomainName> bad_subdomains;
   for (const auto & snm : getParam<std::vector<SubdomainName>>("subdomains"))
   {
     auto sid = MooseMeshUtils::getSubdomainID(snm, *mesh);
     if (!MooseMeshUtils::hasSubdomainID(*mesh, sid))
       bad_subdomains.insert(snm);
     else
       included_subdomains.insert(sid);
   }
   if (!bad_subdomains.empty())
     paramError("subdomains",
                "The requested subdomains do not exist on the fine mesh: ",
                Moose::stringify(bad_subdomains));
 
   // If we are not going through the entire mesh and the ID already exists, we need to offset the
   // assigned ID so that we don't copy ones that are already there
   dof_id_type id_offset = 0;
   if (already_has_id && !included_subdomains.empty())
   {
     for (auto & elem : mesh->active_element_ptr_range())
     {
       dof_id_type elem_id = elem->get_extra_integer(coarse_id);
       if (elem_id != DofObject::invalid_id && elem_id > id_offset)
         id_offset = elem_id;
     }
     id_offset++;
   }
 
   std::unique_ptr<MeshBase> coarse_mesh = std::move(_coarse_mesh);
 
   bool using_subdomain_id = false;
   unsigned int id_for_assignment = 0;
   if (_using_coarse_element_id)
   {
     const auto & id_name = getParam<std::string>("coarse_mesh_extra_element_id");
     if (id_name == "subdomain_id")
       using_subdomain_id = true;
     else
     {
       using_subdomain_id = false;
       if (!coarse_mesh->has_elem_integer(id_name))
         paramError("coarse_mesh_extra_element_id",
                    "The extra element ID does not exist on the coarse mesh");
       else
         id_for_assignment = coarse_mesh->get_elem_integer_index(id_name);
     }
   }
 
   // a relative tolerance on checking if fine mesh is embedded in the coarse mesh
   Real aeps = 0.01;
 
   // The coarse mesh is serialized for distributed mesh to avoid boundary issues.
   // As it is coarse this should be cheap. This will be a null operation for a replicated mesh.
   MeshSerializer tm(*coarse_mesh);
 
   // build a point_locator on coarse mesh
   std::unique_ptr<PointLocatorBase> point_locator =
       PointLocatorBase::build(TREE_ELEMENTS, *coarse_mesh);
   point_locator->enable_out_of_mesh_mode();
 
   // loop through fine mesh elements and get element's centroid
   auto elem_range = included_subdomains.empty()
                         ? mesh->active_element_ptr_range()
                         : mesh->active_subdomain_set_elements_ptr_range(included_subdomains);
   for (auto & elem : elem_range)
   {
     // Get the centroid of the fine elem
     Point centroid = elem->true_centroid();
 
     // Find coarse elem
     const Elem * coarse_elem = (*point_locator)(centroid);
     if (!coarse_elem)
       mooseError("Could not find a coarse element containing a fine element with centroid ",
                  centroid);
 
     // get id from the coarse element
     dof_id_type elem_id;
     if (_using_coarse_element_id)
     {
       if (using_subdomain_id)
         elem_id = coarse_elem->subdomain_id();
       else
         elem_id = coarse_elem->get_extra_integer(id_for_assignment);
     }
     else
       elem_id = coarse_elem->id();
     elem_id += id_offset;
 
     // Check if the fine elem is nested in the coarse element
     for (unsigned int n = 0; n < elem->n_nodes(); n++)
     {
       // Get the node: we need to manually move it towards the centroid to
       // ensure that nothing weird happes due to round-off
       Point current_node = elem->point(n);
       current_node.add_scaled(current_node, -aeps);
       current_node.add_scaled(centroid, aeps);
 
       // Get the element this node is in and check if it is the same
       // as the coarse elem; however check if node_elem is valid as it
       // might not be in case the sub element it outside the coarse domain
       const Elem * node_elem = (*point_locator)(current_node);
       if (!node_elem)
         mooseError("Could not find a coarse element containing a node of fine element at ",
                    elem->point(n));
 
       // get another id from the coarse element with this node
       dof_id_type node_elem_id;
       if (_using_coarse_element_id)
       {
         if (using_subdomain_id)
           node_elem_id = node_elem->subdomain_id();
         else
           node_elem_id = node_elem->get_extra_integer(id_for_assignment);
       }
       else
         node_elem_id = node_elem->id();
 
       if (node_elem_id != elem_id)
         if (_embedding_necessary)
           mooseError(
               "Input mesh is not nested in the coarse mesh in CoarseMeshExtraElementIDGenerator.");
     }
 
     elem->set_extra_integer(coarse_id, elem_id);
   }
 
   return mesh;
 }

◆ validParams()

InputParameters CoarseMeshExtraElementIDGenerator::validParams ( )

static

Definition at line 22 of file CoarseMeshExtraElementIDGenerator.C.

 {
   InputParameters params = MeshGenerator::validParams();
   params.addRequiredParam<MeshGeneratorName>(
       "input", "Name of an existing mesh generator to which we assign coarse element IDs");
   params.addRequiredParam<MeshGeneratorName>(
       "coarse_mesh", "Name of an existing mesh generator as the coarse mesh");
   params.addRequiredParam<std::string>(
       "extra_element_id_name", "Name for the extra element ID that is added to the input mesh");
   params.addParam<std::string>(
       "coarse_mesh_extra_element_id",
       "Name for the extra element ID that is copied from the coarse mesh (default to element ID)");
   params.addParam<std::vector<SubdomainName>>(
       "subdomains", std::vector<SubdomainName>(), "Subdomains to apply extra element IDs to.");
   params.addParam<bool>("enforce_mesh_embedding",
                         false,
                         "True to error out when the input mesh is not embedded in the coarse mesh");
 
   params.addClassDescription("Assign coarse element IDs for elements on a "
                              "mesh based on a coarse mesh.");
   return params;
 }

Member Data Documentation

◆ _coarse_id_name

const std::string& CoarseMeshExtraElementIDGenerator::_coarse_id_name

protected

coarse element ID name

Definition at line 28 of file CoarseMeshExtraElementIDGenerator.h.

Referenced by generate().

◆ _coarse_mesh

std::unique_ptr<MeshBase>& CoarseMeshExtraElementIDGenerator::_coarse_mesh

protected

name of the coarse mesh file

Definition at line 26 of file CoarseMeshExtraElementIDGenerator.h.

Referenced by generate().

◆ _embedding_necessary

const bool& CoarseMeshExtraElementIDGenerator::_embedding_necessary

protected

whether the input mesh must be embedded in the coarse mesh

Definition at line 32 of file CoarseMeshExtraElementIDGenerator.h.

Referenced by generate().

◆ _input

std::unique_ptr<MeshBase>& CoarseMeshExtraElementIDGenerator::_input

protected

input mesh for adding element IDs

Definition at line 24 of file CoarseMeshExtraElementIDGenerator.h.

Referenced by generate().

◆ _using_coarse_element_id

const bool CoarseMeshExtraElementIDGenerator::_using_coarse_element_id

protected

whether or not using the coarse element ID for the extra element ID assignment

Definition at line 30 of file CoarseMeshExtraElementIDGenerator.h.

Referenced by generate().

The documentation for this class was generated from the following files:

reactor/include/meshgenerators/CoarseMeshExtraElementIDGenerator.h
reactor/src/meshgenerators/CoarseMeshExtraElementIDGenerator.C

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Static Public Attributes

Protected Member Functions

Static Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ CoarseMeshExtraElementIDGenerator()

Member Function Documentation

◆ generate()

◆ validParams()

Member Data Documentation

◆ _coarse_id_name

◆ _coarse_mesh

◆ _embedding_necessary

◆ _input

◆ _using_coarse_element_id