EBGeometry_BoundingVolumes.hpp

Go to the documentation of this file.

/* EBGeometry
 * Copyright © 2022 Robert Marskar
 * Please refer to Copyright.txt and LICENSE in the EBGeometry root directory.
 */
 
#ifndef EBGeometry_BoundingVolumes_H
#define EBGeometry_BoundingVolumes_H
 
// Std includes
#include <vector>
#include <iostream>
 
// Our includes
#include "EBGeometry_Vec.hpp"
#include "EBGeometry_NamespaceHeader.hpp"
 
namespace BoundingVolumes {
 
  template <class T>
  class BoundingSphereT
  {
  public:
    friend std::ostream&
    operator<<(std::ostream& os, const BoundingSphereT<T>& sphere)
    {
      os << '(' << sphere.getCentroid() << ", " << sphere.getRadius() << ')';
 
      return os;
    }
 
    enum class BoundingVolumeAlgorithm
    {
      Ritter,
    };
 
    using Vec3 = Vec3T<T>;
 
    inline BoundingSphereT() noexcept;
 
    inline BoundingSphereT(const Vec3T<T>& a_center, const T& a_radius) noexcept;
 
    inline BoundingSphereT(const std::vector<BoundingSphereT<T>>& a_otherSpheres) noexcept;
 
    inline BoundingSphereT(const BoundingSphereT& a_other) noexcept;
 
    template <class P>
    inline BoundingSphereT(const std::vector<Vec3T<P>>&   a_points,
                           const BoundingVolumeAlgorithm& a_alg = BoundingVolumeAlgorithm::Ritter) noexcept;
 
    virtual ~BoundingSphereT() noexcept;
 
    BoundingSphereT&
    operator=(const BoundingSphereT& a_other) = default;
 
    template <class P>
    inline void
    define(const std::vector<Vec3T<P>>& a_points, const BoundingVolumeAlgorithm& a_alg) noexcept;
 
    inline bool
    intersects(const BoundingSphereT& a_other) const noexcept;
 
    inline T&
    getRadius() noexcept;
 
    inline const T&
    getRadius() const noexcept;
 
    inline Vec3&
    getCentroid() noexcept;
 
    inline const Vec3&
    getCentroid() const noexcept;
 
    inline T
    getOverlappingVolume(const BoundingSphereT<T>& a_other) const noexcept;
 
    inline T
    getDistance(const Vec3& a_x0) const noexcept;
 
    inline T
    getVolume() const noexcept;
 
    inline T
    getArea() const noexcept;
 
  protected:
    T m_radius;
 
    Vec3 m_center;
 
    template <class P>
    inline void
    buildRitter(const std::vector<Vec3T<P>>& a_points) noexcept;
  };
 
  template <class T>
  class AABBT
  {
  public:
    friend std::ostream&
    operator<<(std::ostream& os, const AABBT<T>& aabb)
    {
      os << '(' << aabb.getLowCorner() << ", " << aabb.getHighCorner() << ')';
 
      return os;
    }
 
    using Vec3 = Vec3T<T>;
 
    inline AABBT() noexcept;
 
    inline AABBT(const Vec3T<T>& a_lo, const Vec3T<T>& a_hi) noexcept;
 
    inline AABBT(const AABBT& a_other) noexcept;
 
    inline AABBT(const std::vector<AABBT<T>>& a_others) noexcept;
 
    template <class P>
    inline AABBT(const std::vector<Vec3T<P>>& a_points) noexcept;
 
    virtual ~AABBT() noexcept;
 
    AABBT&
    operator=(const AABBT<T>& a_other) = default;
 
    template <class P>
    inline void
    define(const std::vector<Vec3T<P>>& a_points) noexcept;
 
    inline bool
    intersects(const AABBT& a_other) const noexcept;
 
    inline Vec3T<T>&
    getLowCorner() noexcept;
 
    inline const Vec3T<T>&
    getLowCorner() const noexcept;
 
    inline Vec3T<T>&
    getHighCorner() noexcept;
 
    inline const Vec3T<T>&
    getHighCorner() const noexcept;
 
    inline Vec3
    getCentroid() const noexcept;
 
    inline T
    getOverlappingVolume(const AABBT<T>& a_other) const noexcept;
 
    inline T
    getDistance(const Vec3& a_x0) const noexcept;
 
    inline T
    getVolume() const noexcept;
 
    inline T
    getArea() const noexcept;
 
  protected:
    Vec3 m_loCorner;
 
    Vec3 m_hiCorner;
  };
 
  template <class T>
  bool
  intersects(const BoundingSphereT<T>& a_u, const BoundingSphereT<T>& a_v) noexcept;
 
  template <class T>
  bool
  intersects(const AABBT<T>& a_u, const AABBT<T>& a_v) noexcept;
 
  template <class T>
  T
  getOverlappingVolume(const BoundingSphereT<T>& a_u, const BoundingSphereT<T>& a_v) noexcept;
 
  template <class T>
  T
  getOverlappingVolume(const AABBT<T>& a_u, const AABBT<T>& a_v) noexcept;
} // namespace BoundingVolumes
 
#include "EBGeometry_NamespaceFooter.hpp"
 
#include "EBGeometry_BoundingVolumesImplem.hpp"
 
#endif