EBGeometry_MeshDistanceFunctions.hpp

Go to the documentation of this file.

/* EBGeometry
 * Copyright © 2023 Robert Marskar
 * Please refer to Copyright.txt and LICENSE in the EBGeometry root directory.
 */
 
#ifndef EBGeometry_MeshDistanceFunctions
#define EBGeometry_MeshDistanceFunctions
 
// Our includes
#include "EBGeometry_BoundingVolumes.hpp"
#include "EBGeometry_SignedDistanceFunction.hpp"
#include "EBGeometry_DCEL_Mesh.hpp"
#include "EBGeometry_BVH.hpp"
#include "EBGeometry_Triangle.hpp"
#include "EBGeometry_NamespaceHeader.hpp"
 
namespace DCEL {
 
  template <class T, class Meta, class BV, size_t K>
  std::shared_ptr<EBGeometry::BVH::NodeT<T, FaceT<T, Meta>, BV, K>>
  buildFullBVH(const std::shared_ptr<EBGeometry::DCEL::MeshT<T, Meta>>& a_dcelMesh,
               const BVH::Build                                         a_build = BVH::Build::TopDown) noexcept;
} // namespace DCEL
 
template <class T, class Meta = DCEL::DefaultMetaData>
class MeshSDF : public SignedDistanceFunction<T>
{
public:
  using Mesh = EBGeometry::DCEL::MeshT<T, Meta>;
 
  MeshSDF() = delete;
 
  MeshSDF(const std::shared_ptr<Mesh>& a_mesh) noexcept;
 
  virtual ~MeshSDF() = default;
 
  virtual T
  signedDistance(const Vec3T<T>& a_point) const noexcept override;
 
  const std::shared_ptr<Mesh>
  getMesh() const noexcept;
 
  template <class BV>
  BV
  computeBoundingVolume() const noexcept;
 
protected:
  std::shared_ptr<Mesh> m_mesh;
};
 
template <class T, class Meta, class BV, size_t K>
class FastMeshSDF : public SignedDistanceFunction<T>
{
public:
  using Face = typename EBGeometry::DCEL::FaceT<T, Meta>;
 
  using Mesh = typename EBGeometry::DCEL::MeshT<T, Meta>;
 
  using Node = EBGeometry::BVH::NodeT<T, Face, BV, K>;
 
  FastMeshSDF() = delete;
 
  FastMeshSDF(const std::shared_ptr<Mesh>& a_mesh, const BVH::Build a_build = BVH::Build::TopDown) noexcept;
 
  virtual ~FastMeshSDF() = default;
 
  virtual T
  signedDistance(const Vec3T<T>& a_point) const noexcept override;
 
  virtual std::vector<std::pair<std::shared_ptr<const Face>, T>>
  getClosestFaces(const Vec3T<T>& a_point, const bool a_sorted) const noexcept;
 
  virtual std::shared_ptr<Node>&
  getBVH() noexcept;
 
  virtual const std::shared_ptr<Node>&
  getBVH() const noexcept;
 
  BV
  computeBoundingVolume() const noexcept;
 
protected:
  std::shared_ptr<Node> m_bvh;
};
 
template <class T, class Meta, class BV, size_t K>
class FastCompactMeshSDF : public SignedDistanceFunction<T>
{
public:
  using Face = typename EBGeometry::DCEL::FaceT<T, Meta>;
 
  using Mesh = typename EBGeometry::DCEL::MeshT<T, Meta>;
 
  using Root = typename EBGeometry::BVH::LinearBVH<T, Face, BV, K>;
 
  using Node = typename Root::LinearNode;
 
  FastCompactMeshSDF() = delete;
 
  FastCompactMeshSDF(const std::shared_ptr<Mesh>& a_mesh, const BVH::Build a_build = BVH::Build::TopDown) noexcept;
 
  virtual ~FastCompactMeshSDF() = default;
 
  virtual T
  signedDistance(const Vec3T<T>& a_point) const noexcept override;
 
  virtual std::vector<std::pair<std::shared_ptr<const Face>, T>>
  getClosestFaces(const Vec3T<T>& a_point, const bool a_sorted) const noexcept;
 
  virtual std::shared_ptr<Root>&
  getRoot() noexcept;
 
  virtual const std::shared_ptr<Root>&
  getRoot() const noexcept;
 
  BV
  computeBoundingVolume() const noexcept;
 
protected:
  std::shared_ptr<Root> m_bvh;
};
 
template <class T, class Meta, class BV, size_t K>
class FastTriMeshSDF : public SignedDistanceFunction<T>
{
public:
  using Mesh = EBGeometry::DCEL::MeshT<T, Meta>;
 
  using Tri = typename EBGeometry::Triangle<T, Meta>;
 
  using Root = typename EBGeometry::BVH::LinearBVH<T, Tri, BV, K>;
 
  using Node = typename Root::LinearNode;
 
  FastTriMeshSDF() = delete;
 
  FastTriMeshSDF(const std::shared_ptr<Mesh>& a_mesh, const BVH::Build a_build = BVH::Build::TopDown) noexcept;
 
  FastTriMeshSDF(const std::vector<std::shared_ptr<Tri>>& a_triangles,
                 const BVH::Build                         a_build = BVH::Build::TopDown) noexcept;
 
  virtual ~FastTriMeshSDF() = default;
 
  virtual T
  signedDistance(const Vec3T<T>& a_point) const noexcept override;
 
  virtual std::vector<std::pair<std::shared_ptr<const Tri>, T>>
  getClosestTriangles(const Vec3T<T>& a_point, const bool a_sorted) const noexcept;
 
  virtual std::shared_ptr<Root>&
  getRoot() noexcept;
 
  virtual const std::shared_ptr<Root>&
  getRoot() const noexcept;
 
  BV
  computeBoundingVolume() const noexcept;
 
protected:
  std::shared_ptr<Root> m_bvh;
};
 
#include "EBGeometry_NamespaceFooter.hpp"
 
#include "EBGeometry_MeshDistanceFunctionsImplem.hpp"
 
#endif