EBGeometry_SFC.hpp

Go to the documentation of this file.

/* EBGeometry
 * Copyright © 2024 Robert Marskar
 * Please refer to Copyright.txt and LICENSE in the EBGeometry root directory.
 */
 
#ifndef EBGeometry_SFC
#define EBGeometry_SFC
 
// Std includes
#include <cstdint>
 
// Our includes
#include "EBGeometry_NamespaceHeader.hpp"
 
namespace SFC {
 
  using Code = uint64_t;
 
  using Index = std::array<unsigned int, 3>;
 
  static constexpr unsigned int ValidBits = 21;
 
  static constexpr Code ValidSpan = (static_cast<uint64_t>(1) << ValidBits) - 1;
 
#if __cplusplus >= 202002L
  template <typename S>
  concept Encodable = requires(const Index& point, const SFC::Code code) {
    { S::encode(point) } -> std::same_as<SFC::Code>;
    { S::decode(code) } -> std::same_as<Index>;
  };
#endif
 
  struct Morton
  {
    inline static uint64_t
    encode(const Index& a_point) noexcept;
 
    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
  {
    inline static uint64_t
    encode(const Index& a_point) noexcept;
 
    inline static Index
    decode(const uint64_t& a_code) noexcept;
  };
} // namespace SFC
 
#include "EBGeometry_NamespaceFooter.hpp"
 
#if __cplusplus >= 202002L
static_assert(EBGeometry::SFC::Encodable<EBGeometry::SFC::Morton>);
static_assert(EBGeometry::SFC::Encodable<EBGeometry::SFC::Nested>);
#endif
 
#include "EBGeometry_SFCImplem.hpp"
 
#endif