12 #include "libmesh/replicated_mesh.h" 13 #include "libmesh/boundary_info.h" 62 const std::vector<BoundaryName> & boundary_name,
63 bool generate_unknown,
64 const std::set<BoundaryID> & mesh_boundary_ids);
76 const std::vector<BoundaryName> & boundary_name,
77 bool generate_unknown);
90 const std::vector<BoundaryName> & boundary_name,
91 bool generate_unknown);
117 const std::vector<SubdomainName> & subdomain_name);
119 const std::set<SubdomainName> & subdomain_name);
138 template <
typename P,
typename C>
153 "Must pass in a valid rz radial coordinate");
154 factor = 2 * M_PI * point(rz_radial_coord);
158 factor = 4 * M_PI * point(0) * point(0);
172 template <
typename P,
typename C>
176 return (point - origin).cross(direction).norm();
187 template <
typename P,
typename C>
190 const std::pair<Point, RealVectorValue> & axis,
193 factor = 2 * M_PI * computeDistanceToAxis<P, C>(point, axis.first, axis.second);
218 std::unordered_map<dof_id_type, dof_id_type>
220 const std::set<SubdomainID> & block_ids,
221 std::vector<ExtraElementIDName> extra_ids);
231 bool isCoPlanar(
const std::vector<Point> vec_pts,
const Point plane_nvec,
const Point fixed_pt);
239 bool isCoPlanar(
const std::vector<Point> vec_pts,
const Point plane_nvec);
246 bool isCoPlanar(
const std::vector<Point> vec_pts);
274 bool hasSubdomainName(
const MeshBase & input_mesh,
const SubdomainName & name);
288 bool hasBoundaryName(
const MeshBase & input_mesh,
const BoundaryName & name);
301 std::vector<dof_id_type> & elem_id_list,
302 std::vector<dof_id_type> & midpoint_node_list,
303 std::vector<dof_id_type> & ordered_node_list,
304 std::vector<dof_id_type> & ordered_elem_id_list);
315 std::vector<dof_id_type> & elem_id_list,
316 std::vector<dof_id_type> & ordered_node_list,
317 std::vector<dof_id_type> & ordered_elem_id_list);
326 template <
typename T,
typename Q>
332 "'name' ",
name,
" should only contain digits that can be converted to a numerical type.");
333 long long id = std::stoll(
name);
339 " is not within the numeric limits of the expected ID type ",
340 MooseUtils::prettyCppType<Q>(&id_Q),
352 void swapNodesInElem(Elem & elem,
const unsigned int nd1,
const unsigned int nd2);
363 template <
typename T>
366 const std::string & id_name,
367 const std::vector<std::vector<T>> & id_swaps,
368 std::vector<std::unordered_map<T, T>> & id_swap_pairs,
369 const unsigned int row_index_shift = 0)
371 id_swap_pairs.resize(id_swaps.size());
374 const auto & swaps = id_swaps[i];
375 auto & swap_pairs = id_swap_pairs[i];
377 if (swaps.size() % 2)
379 row_index_shift + i + 1,
384 " does not contain an even number of entries! Num entries: ",
387 swap_pairs.reserve(swaps.size() / 2);
388 for (
unsigned int j = 0; j < swaps.size(); j += 2)
389 swap_pairs[swaps[j]] = swaps[j + 1];
402 const std::string & class_name,
403 const unsigned int num_sections,
404 const unsigned int num_integers,
405 const std::vector<std::vector<std::vector<dof_id_type>>> & elem_integers_swaps,
406 std::vector<std::unordered_map<dof_id_type, dof_id_type>> & elem_integers_swap_pairs);
417 std::unique_ptr<ReplicatedMesh>
buildBoundaryMesh(
const ReplicatedMesh & input_mesh,
void swapNodesInElem(Elem &elem, const unsigned int nd1, const unsigned int nd2)
Swap two nodes within an element.
std::string name(const ElemQuality q)
void makeOrderedNodeList(std::vector< std::pair< dof_id_type, dof_id_type >> &node_assm, std::vector< dof_id_type > &elem_id_list, std::vector< dof_id_type > &midpoint_node_list, std::vector< dof_id_type > &ordered_node_list, std::vector< dof_id_type > &ordered_elem_id_list)
Convert a list of sides in the form of a vector of pairs of node ids into a list of ordered nodes bas...
std::unordered_map< dof_id_type, dof_id_type > getExtraIDUniqueCombinationMap(const MeshBase &mesh, const std::set< SubdomainID > &block_ids, std::vector< ExtraElementIDName > extra_ids)
Crate a new set of element-wise IDs by finding unique combinations of existing extra ID values...
std::set< BoundaryID > getBoundaryIDSet(const libMesh::MeshBase &mesh, const std::vector< BoundaryName > &boundary_name, bool generate_unknown)
Gets the boundary IDs into a set with their names.
std::unique_ptr< ReplicatedMesh > buildBoundaryMesh(const ReplicatedMesh &input_mesh, const boundary_id_type boundary_id)
Build a lower-dimensional mesh from a boundary of an input mesh Note: The lower-dimensional mesh will...
const unsigned int invalid_uint
bool hasBoundaryName(const MeshBase &input_mesh, const BoundaryName &name)
Whether a particular boundary name exists in the mesh.
void mooseError(Args &&... args)
Emit an error message with the given stringified, concatenated args and terminate the application...
void coordTransformFactor(const P &point, C &factor, const Moose::CoordinateSystemType coord_type, const unsigned int rz_radial_coord=libMesh::invalid_uint)
compute a coordinate transformation factor
const Point & faceCentroid() const
Returns the coordinates of the face centroid.
void coordTransformFactorRZGeneral(const P &point, const std::pair< Point, RealVectorValue > &axis, C &factor)
Computes a coordinate transformation factor for a general axisymmetric axis.
bool hasBoundaryID(const MeshBase &input_mesh, const BoundaryID id)
Whether a particular boundary ID exists in the mesh.
Real & faceCoord()
Sets/gets the coordinate transformation factor (for e.g.
std::vector< subdomain_id_type > getSubdomainIDs(const libMesh::MeshBase &mesh, const std::vector< SubdomainName > &subdomain_name)
Get the associated subdomainIDs for the subdomain names that are passed in.
SubdomainID getSubdomainID(const SubdomainName &subdomain_name, const MeshBase &mesh)
Gets the subdomain ID associated with the given SubdomainName.
auto max(const L &left, const R &right)
bool hasSubdomainID(const MeshBase &input_mesh, const SubdomainID &id)
Whether a particular subdomain ID exists in the mesh.
C computeDistanceToAxis(const P &point, const Point &origin, const RealVectorValue &direction)
Computes the distance to a general axis.
BoundaryID getBoundaryID(const BoundaryName &boundary_name, const MeshBase &mesh)
Gets the boundary ID associated with the given BoundaryName.
void computeFiniteVolumeCoords(FaceInfo &fi, const Moose::CoordinateSystemType coord_type, const unsigned int rz_radial_coord=libMesh::invalid_uint)
This data structure is used to store geometric and variable related metadata about each cell face in ...
boundary_id_type BoundaryID
void idSwapParametersProcessor(const std::string &class_name, const std::string &id_name, const std::vector< std::vector< T >> &id_swaps, std::vector< std::unordered_map< T, T >> &id_swap_pairs, const unsigned int row_index_shift=0)
Reprocess the swap related input parameters to make pairs out of them to ease further processing...
std::vector< BoundaryID > getBoundaryIDs(const libMesh::MeshBase &mesh, const std::vector< BoundaryName > &boundary_name, bool generate_unknown, const std::set< BoundaryID > &mesh_boundary_ids)
Gets the boundary IDs with their names.
bool isCoPlanar(const std::vector< Point > vec_pts, const Point plane_nvec, const Point fixed_pt)
Decides whether all the Points of a vector of Points are in a plane that is defined by a normal vecto...
bool hasSubdomainName(const MeshBase &input_mesh, const SubdomainName &name)
Whether a particular subdomain name exists in the mesh.
void changeSubdomainId(MeshBase &mesh, const subdomain_id_type old_id, const subdomain_id_type new_id)
Changes the old subdomain ID to a new ID in the mesh.
Provides a way for users to bail out of the current solve.
void extraElemIntegerSwapParametersProcessor(const std::string &class_name, const unsigned int num_sections, const unsigned int num_integers, const std::vector< std::vector< std::vector< dof_id_type >>> &elem_integers_swaps, std::vector< std::unordered_map< dof_id_type, dof_id_type >> &elem_integers_swap_pairs)
Reprocess the elem_integers_swaps into maps so they are easier to use.
Point meshCentroidCalculator(const MeshBase &mesh)
Calculates the centroid of a MeshBase.
void changeBoundaryId(MeshBase &mesh, const libMesh::boundary_id_type old_id, const libMesh::boundary_id_type new_id, bool delete_prev)
Changes the old boundary ID to a new ID in the mesh.
void mergeBoundaryIDsWithSameName(MeshBase &mesh)
Merges the boundary IDs of boundaries that have the same names but different IDs. ...
Q getIDFromName(const T &name)
Converts a given name (BoundaryName or SubdomainName) that is known to only contain digits into a cor...
SubdomainID getNextFreeSubdomainID(MeshBase &input_mesh)
Checks input mesh and returns max(block ID) + 1, which represents a block ID that is not currently in...
bool isDigits(const std::string &str)
Courtesy https://stackoverflow.com/a/8889045 and https://en.cppreference.com/w/cpp/string/byte/isdigi...
BoundaryID getNextFreeBoundaryID(MeshBase &input_mesh)
Checks input mesh and returns the largest boundary ID in the mesh plus one, which is a boundary ID in...
auto min(const L &left, const R &right)
auto index_range(const T &sizable)