MooseMeshElementConversionUtils.h

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

#pragma once

#include "libmesh/replicated_mesh.h"

#include "MooseTypes.h"

namespace MooseMeshElementConversionUtils
{
// Define for comparing the key value of BCTuple
struct BCTupleKeyComp
{
  bool operator()(const libMesh::BoundaryInfo::BCTuple & s, dof_id_type i) const
  {
    return std::get<0>(s) < i;
  }
  bool operator()(dof_id_type i, const libMesh::BoundaryInfo::BCTuple & s) const
  {
    return i < std::get<0>(s);
  }
};

void hexElemSplitter(ReplicatedMesh & mesh,
                     const std::vector<libMesh::BoundaryInfo::BCTuple> & bdry_side_list,
                     const dof_id_type elem_id,
                     std::vector<dof_id_type> & converted_elems_ids);

void pyramidElemSplitter(ReplicatedMesh & mesh,
                         const std::vector<libMesh::BoundaryInfo::BCTuple> & bdry_side_list,
                         const dof_id_type elem_id,
                         std::vector<dof_id_type> & converted_elems_ids);

void prismElemSplitter(ReplicatedMesh & mesh,
                       const std::vector<libMesh::BoundaryInfo::BCTuple> & bdry_side_list,
                       const dof_id_type elem_id,
                       std::vector<dof_id_type> & converted_elems_ids);

void nodeRotationHEX8(const unsigned int min_id_index,
                      const unsigned int sec_min_pos,
                      std::vector<unsigned int> & face_rotation,
                      std::vector<unsigned int> & node_rotation);

std::vector<unsigned int> neighborNodeIndicesHEX8(unsigned int min_id_index);

void hexNodesToTetNodesDeterminer(std::vector<const Node *> & hex_nodes,
                                  std::vector<std::vector<unsigned int>> & rotated_tet_face_indices,
                                  std::vector<std::vector<const Node *>> & tet_nodes_list);

std::vector<bool> quadFaceDiagonalDirectionsHex(const std::vector<const Node *> & hex_nodes);

bool quadFaceDiagonalDirection(const std::vector<const Node *> & quad_nodes);

std::vector<std::vector<unsigned int>>
tetNodesForHex(const std::vector<bool> diagonal_directions,
               std::vector<std::vector<unsigned int>> & tet_face_indices);

void nodeRotationPRISM6(unsigned int min_id_index,
                        std::vector<unsigned int> & face_rotation,
                        std::vector<unsigned int> & node_rotation);

void
prismNodesToTetNodesDeterminer(std::vector<const Node *> & prism_nodes,
                               std::vector<std::vector<unsigned int>> & rotated_tet_face_indices,
                               std::vector<std::vector<const Node *>> & tet_nodes_list);

std::vector<std::vector<unsigned int>>
tetNodesForPrism(const bool diagonal_direction,
                 std::vector<std::vector<unsigned int>> & tet_face_indices);

void nodeRotationPYRAMID5(unsigned int min_id_index,
                          std::vector<unsigned int> & face_rotation,
                          std::vector<unsigned int> & node_rotation);

void
pyramidNodesToTetNodesDeterminer(std::vector<const Node *> & pyramid_nodes,
                                 std::vector<std::vector<unsigned int>> & rotated_tet_face_indices,
                                 std::vector<std::vector<const Node *>> & tet_nodes_list);

void convert3DMeshToAllTet4(ReplicatedMesh & mesh,
                            const std::vector<std::pair<dof_id_type, bool>> & elems_to_process,
                            std::vector<dof_id_type> & converted_elems_ids_to_track,
                            const subdomain_id_type block_id_to_remove,
                            const bool delete_block_to_remove);

void convert3DMeshToAllTet4(ReplicatedMesh & mesh);

void
elementBoundaryInfoCollector(const std::vector<libMesh::BoundaryInfo::BCTuple> & bdry_side_list,
                             const dof_id_type elem_id,
                             const unsigned short n_elem_sides,
                             std::vector<std::vector<boundary_id_type>> & elem_side_list);
}