EBGeometry_BoundingVolumes.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_BOUNDINGVOLUMES_HPP
#define EBGEOMETRY_BOUNDINGVOLUMES_HPP
 
// Std includes
#include <cstddef>
#include <type_traits>
#include <vector>
 
// Our includes
#include "EBGeometry_GPU.hpp"
#include "EBGeometry_Macros.hpp"
#include "EBGeometry_Vec.hpp"
 
namespace EBGeometry {
 
namespace BoundingVolumes {
 
template <class T>
class AABBT
{
  static_assert(std::is_floating_point_v<T>, "AABBT<T>: T must be a floating-point type.");
 
public:
  using Vec3 = Vec3T<T>;
 
  AABBT() noexcept = default;
 
  EBGEOMETRY_HOST_DEVICE
  inline AABBT(const Vec3T<T>& a_lo, const Vec3T<T>& a_hi) noexcept;
 
  AABBT(const AABBT& a_other) noexcept = default;
 
  EBGEOMETRY_HOST
  inline AABBT(const std::vector<AABBT<T>>& a_others) noexcept;
 
  template <class P>
  EBGEOMETRY_HOST
  inline AABBT(const std::vector<Vec3T<P>>& a_points) noexcept;
 
  template <class P>
  EBGEOMETRY_HOST_DEVICE
  inline AABBT(const Vec3T<P>* a_points, size_t a_numPoints) noexcept;
 
  ~AABBT() noexcept = default;
 
  AABBT&
  operator=(const AABBT<T>& a_other) = default;
 
  AABBT(AABBT&& a_other) noexcept = default;
 
  AABBT&
  operator=(AABBT&& a_other) noexcept = default;
 
  template <class P>
  EBGEOMETRY_HOST
  inline void
  define(const std::vector<Vec3T<P>>& a_points) noexcept;
 
  template <class P>
  EBGEOMETRY_HOST_DEVICE
  inline void
  define(const Vec3T<P>* a_points, size_t a_numPoints) noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline bool
  intersects(const AABBT& a_other) const noexcept;
 
  EBGEOMETRY_HOST_DEVICE
  inline Vec3T<T>&
  getLowCorner() noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline const Vec3T<T>&
  getLowCorner() const noexcept;
 
  EBGEOMETRY_HOST_DEVICE
  inline Vec3T<T>&
  getHighCorner() noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline const Vec3T<T>&
  getHighCorner() const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline Vec3
  getCentroid() const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline T
  getOverlappingVolume(const AABBT<T>& a_other) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline T
  getDistance(const Vec3& a_x0) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline T
  getDistance2(const Vec3& a_x0) const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline T
  getVolume() const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline T
  getArea() const noexcept;
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline AABBT<T>
  merged(const AABBT<T>& a_other) const noexcept;
 
protected:
  Vec3 m_loCorner = Vec3::infinity();
 
  Vec3 m_hiCorner = -Vec3::infinity();
};
 
template <class T>
[[nodiscard]] EBGEOMETRY_HOST_DEVICE
bool
intersects(const AABBT<T>& a_u, const AABBT<T>& a_v) noexcept;
 
template <class T>
[[nodiscard]] EBGEOMETRY_HOST_DEVICE
T
getOverlappingVolume(const AABBT<T>& a_u, const AABBT<T>& a_v) noexcept;
 
static_assert(std::is_trivially_copyable_v<AABBT<float>>, "AABBT<float> must be trivially copyable");
static_assert(std::is_trivially_copyable_v<AABBT<double>>, "AABBT<double> must be trivially copyable");
} // namespace BoundingVolumes
 
} // namespace EBGeometry
 
#include "EBGeometry_BoundingVolumesImplem.hpp"
 
#endif