EBGeometry_Math.hpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: 2026 Robert Marskar <robert.marskar@sintef.no>
//
// SPDX-License-Identifier: GPL-3.0-or-later
 
#ifndef EBGEOMETRY_MATH_HPP
#define EBGEOMETRY_MATH_HPP
 
#include <cfloat>
#include <climits>
#include <cmath>
 
#include "EBGeometry_GPU.hpp"
 
namespace EBGeometry {
 
namespace Math {
 
template <class T>
[[nodiscard]] EBGEOMETRY_HOST_DEVICE
constexpr T
min(const T& a_a, const T& a_b) noexcept
{
  return (a_b < a_a) ? a_b : a_a;
}
 
template <class T>
[[nodiscard]] EBGEOMETRY_HOST_DEVICE
constexpr T
max(const T& a_a, const T& a_b) noexcept
{
  return (a_a < a_b) ? a_b : a_a;
}
 
template <class T>
[[nodiscard]] EBGEOMETRY_HOST_DEVICE
constexpr T
clamp(const T& a_value, const T& a_lo, const T& a_hi) noexcept
{
  return (a_value < a_lo) ? a_lo : ((a_hi < a_value) ? a_hi : a_value);
}
 
template <class T>
struct Limits;
 
template <>
struct Limits<float>
{
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  static constexpr float
  max() noexcept
  {
    return FLT_MAX;
  }
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  static constexpr float
  lowest() noexcept
  {
    return -FLT_MAX;
  }
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  static constexpr float
  min() noexcept
  {
    return FLT_MIN;
  }
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  static constexpr float
  epsilon() noexcept
  {
    return FLT_EPSILON;
  }
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  static constexpr float
  infinity() noexcept
  {
    return HUGE_VALF;
  }
};
 
template <>
struct Limits<double>
{
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  static constexpr double
  max() noexcept
  {
    return DBL_MAX;
  }
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  static constexpr double
  lowest() noexcept
  {
    return -DBL_MAX;
  }
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  static constexpr double
  min() noexcept
  {
    return DBL_MIN;
  }
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  static constexpr double
  epsilon() noexcept
  {
    return DBL_EPSILON;
  }
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  static constexpr double
  infinity() noexcept
  {
    return HUGE_VAL;
  }
};
 
template <class T, T Lowest, T Max>
struct IntegerLimits
{
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  static constexpr T
  max() noexcept
  {
    return Max;
  }
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  static constexpr T
  lowest() noexcept
  {
    return Lowest;
  }
 
  [[nodiscard]] EBGEOMETRY_HOST_DEVICE
  static constexpr T
  min() noexcept
  {
    return Lowest;
  }
};
 
template <>
struct Limits<int> : IntegerLimits<int, INT_MIN, INT_MAX>
{
};
 
template <>
struct Limits<unsigned int> : IntegerLimits<unsigned int, 0U, UINT_MAX>
{
};
 
template <>
struct Limits<long> : IntegerLimits<long, LONG_MIN, LONG_MAX>
{
};
 
template <>
struct Limits<unsigned long> : IntegerLimits<unsigned long, 0UL, ULONG_MAX>
{
};
 
template <>
struct Limits<long long> : IntegerLimits<long long, LLONG_MIN, LLONG_MAX>
{
};
 
template <>
struct Limits<unsigned long long> : IntegerLimits<unsigned long long, 0ULL, ULLONG_MAX>
{
};
 
} // namespace Math
 
} // namespace EBGeometry
 
#endif