EBGeometry_SphereAoSoA.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_SPHEREAOSOA_HPP
#define EBGEOMETRY_SPHEREAOSOA_HPP
 
// Std includes
#include <cstddef>
#include <cstdint>
#include <type_traits>
 
// Our includes
#include "EBGeometry_Array.hpp"
#include "EBGeometry_PointAoSoA.hpp"
#include "EBGeometry_SphereCloud.hpp"
#include "EBGeometry_Vec.hpp"
 
namespace EBGeometry {
 
template <class T, size_t W = PointSoA::DefaultWidth<T>()>
struct SphereAoSoA
{
  static_assert(W > 0, "W must be positive");
  static_assert(std::is_floating_point_v<T>, "SphereAoSoA requires a floating-point type T");
 
public:
  EBGEOMETRY_HOST
  void
  pack(const SphereCloud::Sphere<T>* a_spheres, const uint32_t* a_sphereIds, uint32_t a_count) noexcept;
 
  EBGEOMETRY_HOST_DEVICE
  void
  packUnchecked(const SphereCloud::Sphere<T>* a_spheres, const uint32_t* a_sphereIds, uint32_t a_count) noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  Array<T, W>
  getSignedDistances(const Vec3T<T>& a_point) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  uint32_t
  getSphereId(size_t a_lane) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  uint32_t
  numValid() const noexcept
  {
    return m_centers.numValid();
  }
 
  template <class BV>
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  BV
  computeBoundingVolume() const noexcept;
 
protected:
  PointAoSoA<T, W> m_centers;
 
  Array<T, W> m_radii;
};
 
static_assert(std::is_trivially_copyable_v<SphereAoSoA<float>>, "SphereAoSoA<float> must be trivially copyable");
static_assert(std::is_trivially_copyable_v<SphereAoSoA<double>>, "SphereAoSoA<double> must be trivially copyable");
 
} // namespace EBGeometry
 
#include "EBGeometry_SphereAoSoAImplem.hpp"
 
#endif