LCC08CellHandlerUtility.h

Go to the documentation of this file.

 
#pragma once
 
#include <array>
#include <tuple>
#include <type_traits>
#include <utility>
#include <vector>
 
namespace autopas::LCC08CellHandlerUtility {
 
using OffsetPairSorting = std::tuple<unsigned long, unsigned long, std::array<double, 3>>;
 
using OffsetTripletSorting = std::tuple<unsigned long, unsigned long, unsigned long, std::array<double, 3>>;
 
using OffsetPair = std::pair<unsigned long, unsigned long>;
 
using OffsetTriplet = std::tuple<unsigned long, unsigned long, unsigned long>;
 
using OffsetPairVector = std::vector<OffsetPair>;
 
using OffsetTripletVector = std::vector<OffsetTriplet>;
 
enum class C08OffsetMode {
  noSorting = 0,
  sorting = 1,
  c04NoSorting = 2,
};
 
template <C08OffsetMode Mode>
using OffsetPairType = std::vector<
    std::conditional_t<Mode == C08OffsetMode::sorting, OffsetPairSorting,
                       std::conditional_t<Mode == C08OffsetMode::c04NoSorting, OffsetPairVector, OffsetPair>>>;
 
template <C08OffsetMode Mode>
using OffsetTripletType = std::vector<
    std::conditional_t<Mode == C08OffsetMode::sorting, OffsetTripletSorting,
                       std::conditional_t<Mode == C08OffsetMode::c04NoSorting, OffsetTripletVector, OffsetTriplet>>>;
 
namespace internal {
 
enum class C08CellDirection : int {
  frontLeft,
  backLeft,
  frontRight,
  backRight,
};
 
constexpr inline std::array<C08CellDirection, 4> ALL_DIRECTIONS{
    {C08CellDirection::frontLeft, C08CellDirection::backLeft, C08CellDirection::frontRight,
     C08CellDirection::backRight}};
 
constexpr inline char ENUM_EXTENSION_EXCEPTION[]{
    "Enum C08CellDirection was extended, but its associated switch-case statements was not!"};
 
constexpr std::pair<int, int> toMaskXY(const C08CellDirection &direction);
 
constexpr bool includeCellPair(const C08CellDirection &direction, const std::array<int, 3> &overlap, int x, int y,
                               int z);
 
std::array<double, 3> computeSortingDirection(const std::array<double, 3> &offset1Vector,
                                              const std::array<double, 3> &offset2Vector,
                                              const std::array<double, 3> &cellLength);
 
}  // namespace internal
 
template <C08OffsetMode Mode>
OffsetPairType<Mode> computePairwiseCellOffsetsC08(const std::array<unsigned long, 3> &cellsPerDimension,
                                                   const std::array<double, 3> &cellLength, double interactionLength);
 
template <C08OffsetMode Mode>
OffsetTripletType<Mode> computeTriwiseCellOffsetsC08(const std::array<unsigned long, 3> &cellsPerDimension,
                                                     const std::array<double, 3> &cellLength, double interactionLength);
 
}  // namespace autopas::LCC08CellHandlerUtility