Select Git revision
InterpolateItemValue.hpp
MeshNodeBoundary.hpp 12.91 KiB
#ifndef MESH_NODE_BOUNDARY_HPP
#define MESH_NODE_BOUNDARY_HPP
#include <Array.hpp>
#include <ItemValue.hpp>
#include <Kokkos_Vector.hpp>
#include <TinyVector.hpp>
#include <RefItemList.hpp>
#include <ConnectivityMatrix.hpp>
#include <IConnectivity.hpp>
#include <iostream>
#include <Messenger.hpp>
template <size_t Dimension>
class MeshNodeBoundary
{
protected:
Array<const NodeId> m_node_list;
public:
MeshNodeBoundary& operator=(const MeshNodeBoundary&) = default;
MeshNodeBoundary& operator=(MeshNodeBoundary&&) = default;
const Array<const NodeId>&
nodeList() const
{
return m_node_list;
}
template <typename MeshType>
MeshNodeBoundary(const MeshType& mesh, const RefFaceList& ref_face_list)
{
static_assert(Dimension == MeshType::Dimension);
const auto& face_to_cell_matrix = mesh.connectivity().faceToCellMatrix();
const Array<const FaceId>& face_list = ref_face_list.list();
parallel_for(face_list.size(), PUGS_LAMBDA(const int& l) {
const auto& face_cells = face_to_cell_matrix[face_list[l]];
if (face_cells.size() > 1) {
perr() << "internal faces cannot be used to define mesh boundaries\n";
std::exit(1);
}
});
Kokkos::vector<unsigned int> node_ids;
// not enough but should reduce significantly the number of resizing
node_ids.reserve(Dimension * face_list.size());
const auto& face_to_node_matrix = mesh.connectivity().faceToNodeMatrix();
for (size_t l = 0; l < face_list.size(); ++l) {
const FaceId face_number = face_list[l];
const auto& face_nodes = face_to_node_matrix[face_number];
for (size_t r = 0; r < face_nodes.size(); ++r) {
node_ids.push_back(face_nodes[r]);
}
}
std::sort(node_ids.begin(), node_ids.end());
auto last = std::unique(node_ids.begin(), node_ids.end());
node_ids.resize(std::distance(node_ids.begin(), last));
Array<NodeId> node_list(node_ids.size());
parallel_for(node_ids.size(),
PUGS_LAMBDA(const int& r) { node_list[r] = node_ids[r]; });
m_node_list = node_list;
}