EBGeometry_PointAoSoA.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_POINTAOSOA_HPP
#define EBGEOMETRY_POINTAOSOA_HPP
 
// Std includes
#include <cstddef>
#include <cstdint>
#include <type_traits>
 
// Our includes
#include "EBGeometry_Array.hpp"
#include "EBGeometry_PointSoA.hpp"
#include "EBGeometry_Vec.hpp"
 
namespace EBGeometry {
 
template <class T, size_t W = PointSoA::DefaultWidth<T>()>
struct PointAoSoA
{
  static_assert(W > 0, "W must be positive");
  static_assert(std::is_floating_point_v<T>, "PointAoSoA requires a floating-point type T");
 
public:
  EBGEOMETRY_HOST
  void
  pack(const Vec3T<T>* a_positions, const uint32_t* a_pointIds, uint32_t a_count) noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  Array<T, W>
  getDistances2(const Vec3T<T>& a_point) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  Array<T, W>
  getDistances(const Vec3T<T>& a_point) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  T
  getMinimumDistance2(const Vec3T<T>& a_point) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  T
  getMinimumDistance(const Vec3T<T>& a_point) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  T
  getMaximumDistance2(const Vec3T<T>& a_point) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  T
  getMaximumDistance(const Vec3T<T>& a_point) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  uint32_t
  getPointId(size_t a_lane) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  uint32_t
  numValid() const noexcept
  {
    return m_positions.numValid();
  }
 
  template <class BV>
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  BV
  computeBoundingVolume() const noexcept;
 
protected:
  PointSoAT<T, W> m_positions;
 
  Array<uint32_t, W> m_pointIds;
};
 
static_assert(std::is_trivially_copyable_v<PointAoSoA<float>>, "PointAoSoA<float> must be trivially copyable");
static_assert(std::is_trivially_copyable_v<PointAoSoA<double>>, "PointAoSoA<double> must be trivially copyable");
 
} // namespace EBGeometry
 
#include "EBGeometry_PointAoSoAImplem.hpp"
 
#endif