EBGeometry_SFC.hpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: 2024 Robert Marskar <robert.marskar@sintef.no>
//
// SPDX-License-Identifier: GPL-3.0-or-later
 
#ifndef EBGEOMETRY_SFC_HPP
#define EBGEOMETRY_SFC_HPP
 
// Std includes
#include <cstdint>
#include <vector>
 
// Our includes
#include "EBGeometry_Array.hpp"
#include "EBGeometry_GPU.hpp"
#include "EBGeometry_Vec.hpp"
 
namespace EBGeometry {
 
namespace SFC {
 
using Code = uint64_t;
 
using Index = Array<unsigned int, 3>;
 
static constexpr unsigned int ValidBits = 21;
 
static constexpr Code ValidSpan = (static_cast<uint64_t>(1) << ValidBits) - 1;
 
struct Morton
{
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  inline static uint64_t
  encode(const Index& a_point) noexcept;
 
  [[nodiscard]] inline static Index
  decode(const uint64_t& a_code) noexcept;
 
protected:
  static constexpr uint_fast64_t Mask_64[6]{
    0x1fffff, 0x1f00000000ffff, 0x1f0000ff0000ff, 0x100f00f00f00f00f, 0x10c30c30c30c30c3, 0x1249249249249249};
};
 
struct Nested
{
  [[nodiscard]] inline static uint64_t
  encode(const Index& a_point) noexcept;
 
  [[nodiscard]] inline static Index
  decode(const uint64_t& a_code) noexcept;
};
 
struct Hilbert
{
  [[nodiscard]] inline static uint64_t
  encode(const Index& a_point) noexcept;
 
  [[nodiscard]] inline static Index
  decode(const uint64_t& a_code) noexcept;
};
 
template <class T>
[[nodiscard]] inline std::vector<Index>
computeBins(const std::vector<Vec3T<T>>& a_points);
 
template <class Curve = Morton, class T>
[[nodiscard]] inline std::vector<uint32_t>
order(const std::vector<Vec3T<T>>& a_points);
 
} // namespace SFC
 
} // namespace EBGeometry
 
#include "EBGeometry_SFCImplem.hpp"
 
#endif