EBGeometry_TriangleAoSoA.hpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: 2026 Robert Marskar <robert.marskar@sintef.no>
//
// SPDX-License-Identifier: GPL-3.0-or-later
 
#ifndef EBGEOMETRY_TRIANGLEAOSOA_HPP
#define EBGEOMETRY_TRIANGLEAOSOA_HPP
 
// Std includes
#include <cstddef>
#include <cstdint>
#include <type_traits>
 
// Our includes
#include "EBGeometry_Array.hpp"
#include "EBGeometry_Triangle.hpp"
#include "EBGeometry_TriangleSoA.hpp"
#include "EBGeometry_Vec.hpp"
 
namespace EBGeometry {
 
template <class T, size_t W = TriangleSoA::DefaultWidth<T>()>
struct TriangleAoSoA
{
  static_assert(W > 0, "W must be positive");
  static_assert(std::is_floating_point_v<T>, "TriangleAoSoA requires a floating-point type T");
 
public:
  EBGEOMETRY_HOST
  void
  pack(const Triangle<T>* a_triangles, uint32_t a_count) noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  T
  signedDistance(const Vec3T<T>& a_point) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  T
  signedDistance(const Vec3T<T>& a_point, uint32_t& a_closestFaceId) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  uint32_t
  nearestLane(const Vec3T<T>& a_point, T& a_d2) const noexcept
  {
    return m_triangles.nearestLane(a_point, a_d2);
  }
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  T
  signedDistanceAt(uint32_t a_lane, const Vec3T<T>& a_point) const noexcept
  {
    return m_triangles.signedDistanceAt(a_lane, a_point);
  }
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  T
  rayCrossing(const Vec3T<T>& a_origin,
              const Vec3T<T>& a_direction,
              T               a_tMin,
              T               a_tMax,
              uint32_t&       a_faceId,
              Vec3T<T>&       a_normal) const noexcept;
 
  EBGEOMETRY_HOST_DEVICE
  void
  rayCrossings(
    const Vec3T<T>& a_origin, const Vec3T<T>& a_direction, T a_tMin, T a_tMax, T* a_crossings) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  Vec3T<T>
  getNormal(size_t a_lane) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  uint32_t
  getFaceId(size_t a_lane) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline Vec3T<T>
  getVertex(size_t a_lane, size_t a_vertex) const noexcept
  {
    return m_triangles.getVertex(a_lane, a_vertex);
  }
 
  template <class BV>
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  BV
  computeBoundingVolume() const noexcept;
 
protected:
  TriangleSoAT<T, W> m_triangles;
 
  Array<uint32_t, W> m_faceIds;
 
  uint32_t m_validCount = 0;
};
 
static_assert(std::is_trivially_copyable_v<TriangleAoSoA<float>>, "TriangleAoSoA<float> must be trivially copyable");
static_assert(std::is_trivially_copyable_v<TriangleAoSoA<double>>, "TriangleAoSoA<double> must be trivially copyable");
 
} // namespace EBGeometry
 
#include "EBGeometry_TriangleAoSoAImplem.hpp"
 
#endif