EBGeometry_DCEL_Edge.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_EDGE_HPP
#define EBGEOMETRY_DCEL_EDGE_HPP
 
// Std includes
#include <cstdint>
#include <type_traits>
 
// Our includes
#include "EBGeometry_DCEL.hpp"
#include "EBGeometry_DCEL_Face.hpp"
#include "EBGeometry_GPU.hpp"
#include "EBGeometry_Vec.hpp"
 
namespace EBGeometry {
 
namespace DCEL {
 
template <class T>
class EdgeT
{
  static_assert(std::is_floating_point_v<T>, "EdgeT<T>: T must be a floating-point type");
 
public:
  using Vec3 = Vec3T<T>;
 
  using Vertex = VertexT<T>;
 
  using Edge = EdgeT<T>;
 
  using Face = FaceT<T>;
 
  using Mesh = MeshT<T>;
 
  EdgeT() noexcept = default;
 
  EdgeT(const Edge& a_otherEdge) = default;
 
  EdgeT(Edge&& a_otherEdge) = default;
 
  EBGEOMETRY_HOST_DEVICE
  EdgeT(const uint32_t a_vertexIndex) noexcept;
 
  ~EdgeT() = default;
 
  Edge&
  operator=(const Edge& a_otherEdge) = default;
 
  Edge&
  operator=(Edge&& a_otherEdge) = default;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline size_t
  size() const noexcept;
 
  EBGEOMETRY_HOST_DEVICE
  inline void
  define(const uint32_t a_vertexIndex, const uint32_t a_pairEdgeIndex, const uint32_t a_nextEdgeIndex) noexcept;
 
  EBGEOMETRY_HOST_DEVICE
  inline void
  reconcile(const Mesh& a_mesh) noexcept;
 
  EBGEOMETRY_HOST_DEVICE
  inline void
  flipNormal() noexcept;
 
  EBGEOMETRY_HOST_DEVICE
  inline void
  setVertex(const uint32_t a_vertexIndex) noexcept;
 
  EBGEOMETRY_HOST_DEVICE
  inline void
  setPairEdge(const uint32_t a_pairEdgeIndex) noexcept;
 
  EBGEOMETRY_HOST_DEVICE
  inline void
  setNextEdge(const uint32_t a_nextEdgeIndex) noexcept;
 
  EBGEOMETRY_HOST_DEVICE
  inline void
  setFace(const uint32_t a_faceIndex) noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline uint32_t
  getVertexIndex() const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline uint32_t
  getPairEdgeIndex() const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline uint32_t
  getNextEdgeIndex() const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline uint32_t
  getFaceIndex() const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline Vertex&
  getVertex(Mesh& a_mesh) noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline const Vertex&
  getVertex(const Mesh& a_mesh) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline Vertex&
  getOtherVertex(Mesh& a_mesh) noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline const Vertex&
  getOtherVertex(const Mesh& a_mesh) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline Edge&
  getPairEdge(Mesh& a_mesh) noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline const Edge&
  getPairEdge(const Mesh& a_mesh) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline Edge&
  getNextEdge(Mesh& a_mesh) noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline const Edge&
  getNextEdge(const Mesh& a_mesh) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline Vec3T<T>
  computeNormal(const Mesh& a_mesh) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline Vec3T<T>&
  getNormal() noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline const Vec3T<T>&
  getNormal() const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline Face&
  getFace(Mesh& a_mesh) noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline const Face&
  getFace(const Mesh& a_mesh) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline T
  signedDistance(const Vec3& a_x0, const Mesh& a_mesh) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline T
  unsignedDistance2(const Vec3& a_x0, const Mesh& a_mesh) const noexcept;
 
protected:
  Vec3 m_normal = Vec3::zeros();
 
  uint32_t m_vertex = UINT32_MAX;
 
  uint32_t m_pairEdge = UINT32_MAX;
 
  uint32_t m_nextEdge = UINT32_MAX;
 
  uint32_t m_face = UINT32_MAX;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline T
  projectPointToEdge(const Vec3& a_x0, const Mesh& a_mesh) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline Vec3T<T>
  getX2X1(const Mesh& a_mesh) const noexcept;
};
 
static_assert(std::is_trivially_copyable_v<EdgeT<float>>, "EdgeT<float> must be trivially copyable");
static_assert(std::is_trivially_copyable_v<EdgeT<double>>, "EdgeT<double> must be trivially copyable");
 
} // namespace DCEL
 
} // namespace EBGeometry
 
#include "EBGeometry_DCEL_EdgeImplem.hpp"
 
#endif