inBox.h

Go to the documentation of this file.

 
#pragma once
 
#include <array>
#include <type_traits>
 
namespace autopas::utils {
template <typename T>
bool inBox(const std::array<T, 3> &position, const std::array<T, 3> &low, const std::array<T, 3> &high) {
  static_assert(std::is_floating_point<T>::value, "inBox assumes floating point types");
 
  bool inBox = true;
  for (int d = 0; d < 3; ++d) {
    const bool isLargerThanLower = position[d] >= low[d];
    const bool isSmallerThanHigher = position[d] < high[d];
    inBox = inBox and isLargerThanLower and isSmallerThanHigher;
  }
  return inBox;
}
 
template <typename T>
bool notInBox(const std::array<T, 3> &position, const std::array<T, 3> &low, const std::array<T, 3> &high) {
  return not(inBox(position, low, high));
}
 
template <typename T>
bool boxesOverlap(const std::array<T, 3> &boxALow, const std::array<T, 3> &boxAHigh, const std::array<T, 3> &boxBLow,
                  const std::array<T, 3> &boxBHigh) {
  static_assert(std::is_floating_point_v<T>, "boxesOverlap assumes floating point types");
 
  auto overlap1D = [&](size_t dim) { return boxAHigh[dim] > boxBLow[dim] and boxBHigh[dim] > boxALow[dim]; };
 
  return overlap1D(0) and overlap1D(1) and overlap1D(2);
}
}  // namespace autopas::utils