Skip to content
Snippets Groups Projects
Select Git revision
  • abae271d222231b934d7e5a0191b8e58b10a0fe7
  • develop default protected
  • feature/variational-hydro
  • origin/stage/bouguettaia
  • feature/gmsh-reader
  • feature/reconstruction
  • save_clemence
  • feature/kinetic-schemes
  • feature/local-dt-fsi
  • feature/composite-scheme-sources
  • feature/composite-scheme-other-fluxes
  • feature/serraille
  • feature/composite-scheme
  • hyperplastic
  • feature/polynomials
  • feature/gks
  • feature/implicit-solver-o2
  • feature/coupling_module
  • feature/implicit-solver
  • feature/merge-local-dt-fsi
  • master protected
  • v0.5.0 protected
  • v0.4.1 protected
  • v0.4.0 protected
  • v0.3.0 protected
  • v0.2.0 protected
  • v0.1.0 protected
  • Kidder
  • v0.0.4 protected
  • v0.0.3 protected
  • v0.0.2 protected
  • v0 protected
  • v0.0.1 protected
33 results

InterpolateItemValue.hpp

Blame
  • 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;
      }