EBGeometry_DCEL_Mesh.hpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: 2022 Robert Marskar <robert.marskar@sintef.no>
//
// SPDX-License-Identifier: GPL-3.0-or-later
 
#ifndef EBGEOMETRY_DCEL_MESH_HPP
#define EBGEOMETRY_DCEL_MESH_HPP
 
// Std includes
#include <cstddef>
#include <cstdint>
#include <map>
#include <memory>
#include <string>
#include <type_traits>
#include <vector>
 
// Our includes
#include "EBGeometry_DCEL.hpp"
#include "EBGeometry_DCEL_Edge.hpp"
#include "EBGeometry_DCEL_Face.hpp"
#include "EBGeometry_DCEL_Vertex.hpp"
#include "EBGeometry_GPU.hpp"
#include "EBGeometry_PODVector.hpp"
#include "EBGeometry_Pool.hpp"
 
namespace EBGeometry {
 
namespace DCEL {
 
template <class T>
class MeshT : public PoolHandle
{
  static_assert(std::is_floating_point_v<T>, "MeshT requires a floating-point T");
 
public:
  enum class SearchAlgorithm
  {
    Direct, 
    Direct2 
  };
 
  using Vec3 = Vec3T<T>;
 
  using Vertex = VertexT<T>;
 
  using Edge = EdgeT<T>;
 
  using Face = FaceT<T>;
 
  using Mesh = MeshT<T>;
 
  MeshT() noexcept = default;
 
  struct Layout
  {
    SearchAlgorithm m_algorithm = SearchAlgorithm::Direct2;
 
    PODVector<Vertex> m_vertices;
 
    PODVector<Edge> m_edges;
 
    PODVector<Face> m_faces;
  };
 
  EBGEOMETRY_HOST_DEVICE
  MeshT(const Layout& a_layout, const PoolLocation& a_location) noexcept
    : PoolHandle(a_location),
      m_algorithm(a_layout.m_algorithm),
      m_vertices(a_layout.m_vertices),
      m_edges(a_layout.m_edges),
      m_faces(a_layout.m_faces)
  {}
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  Layout
  layout() const noexcept
  {
    return Layout{m_algorithm, m_vertices, m_edges, m_faces};
  }
 
  MeshT(const Mesh& a_otherMesh) noexcept = default;
 
  MeshT(Mesh&& a_otherMesh) noexcept = default;
 
  ~MeshT() noexcept = default;
 
  Mesh&
  operator=(const Mesh& a_otherMesh) noexcept = default;
 
  Mesh&
  operator=(Mesh&& a_otherMesh) noexcept = default;
 
  [[nodiscard]] EBGEOMETRY_HOST
  inline Mesh
  deepCopy(Pool& a_dstPool) const;
 
  EBGEOMETRY_HOST
  inline void
  sanityCheck(const std::string a_id) const;
 
  EBGEOMETRY_HOST_DEVICE
  inline void
  setSearchAlgorithm(const SearchAlgorithm a_algorithm) noexcept;
 
  EBGEOMETRY_HOST_DEVICE
  inline void
  setInsideOutsideAlgorithm(InsideOutsideAlgorithm a_algorithm) noexcept;
 
  EBGEOMETRY_HOST
  inline void
  reconcile();
 
  EBGEOMETRY_HOST_DEVICE
  inline void
  flip() noexcept;
 
  EBGEOMETRY_HOST
  inline void
  reserveVertices(Pool& a_pool, uint32_t a_capacity);
 
  EBGEOMETRY_HOST
  inline void
  reserveEdges(Pool& a_pool, uint32_t a_capacity);
 
  EBGEOMETRY_HOST
  inline void
  reserveFaces(Pool& a_pool, uint32_t a_capacity);
 
  EBGEOMETRY_HOST
  inline uint32_t
  addVertex(Pool& a_pool, const Vertex& a_vertex);
 
  EBGEOMETRY_HOST
  inline uint32_t
  addEdge(Pool& a_pool, const Edge& a_edge);
 
  EBGEOMETRY_HOST
  inline uint32_t
  addFace(Pool& a_pool, const Face& a_face);
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline Vertex&
  getVertex(uint32_t a_index) noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline const Vertex&
  getVertex(uint32_t a_index) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline Edge&
  getEdge(uint32_t a_index) noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline const Edge&
  getEdge(uint32_t a_index) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline Face&
  getFace(uint32_t a_index) noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline const Face&
  getFace(uint32_t a_index) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline uint32_t
  numVertices() const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline uint32_t
  numEdges() const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline uint32_t
  numFaces() const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST
  inline std::vector<Vec3T<T>>
  getAllVertexCoordinates() const;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline T
  signedDistance(const Vec3& a_x0) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline T
  signedDistance(const Vec3& a_x0, SearchAlgorithm a_algorithm) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline T
  unsignedDistance2(const Vec3& a_x0) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST
  inline Mesh
  rebasedView(const Pool& a_pool) const noexcept;
 
protected:
  SearchAlgorithm m_algorithm = SearchAlgorithm::Direct2;
 
  PODVector<Vertex> m_vertices;
 
  PODVector<Edge> m_edges;
 
  PODVector<Face> m_faces;
 
  EBGEOMETRY_HOST
  inline void
  attachTo(const Pool& a_pool) noexcept;
 
  EBGEOMETRY_HOST
  inline void
  reconcileFaces() noexcept;
 
  EBGEOMETRY_HOST_DEVICE
  inline void
  reconcileEdges() noexcept;
 
  EBGEOMETRY_HOST
  inline void
  reconcileVertices();
 
  EBGEOMETRY_HOST_DEVICE
  inline void
  flipFaceNormals() noexcept;
 
  EBGEOMETRY_HOST_DEVICE
  inline void
  flipEdgeNormals() noexcept;
 
  EBGEOMETRY_HOST_DEVICE
  inline void
  flipVertexNormals() noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline T
  DirectSignedDistance(const Vec3& a_point) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline T
  DirectSignedDistance2(const Vec3& a_point) const noexcept;
 
  EBGEOMETRY_HOST
  inline void
  incrementWarning(std::map<std::string, size_t>& a_warnings, const std::string& a_warn) const;
 
  EBGEOMETRY_HOST
  inline void
  printWarnings(const std::map<std::string, size_t>& a_warnings, const std::string& a_id) const;
 
  EBGEOMETRY_HOST
  inline bool
  faceLoopCloses(uint32_t a_faceIndex, uint32_t& a_lastEdge, uint32_t& a_steps) const noexcept;
};
} // namespace DCEL
 
} // namespace EBGeometry
 
#include "EBGeometry_DCEL_MeshImplem.hpp"
 
#endif