traits.hpp

#ifndef __COMMON_TRAITS_HPP__
#define __COMMON_TRAITS_HPP__
#include <cmath>
#include <type_traits>
 
#if defined __GNUC__
#  define LIKELY(EXPR) __builtin_expect(!!(EXPR), 1)
#else
#  define LIKELY(EXPR) (!!(EXPR))
#endif
 
#if defined NDEBUG
#  define X_ASSERT(CHECK) void(0)
#else
#  define X_ASSERT(CHECK)                                                      \
    (LIKELY(CHECK) ? void(0) : [] { assert(!(#CHECK)); }())
#endif
 
template <typename T>
concept FloatingPointType = std::is_floating_point_v<std::remove_cvref_t<T>>;
 
static_assert(FloatingPointType<double>, "double ok");
static_assert(FloatingPointType<float>, "float ok");
static_assert(!FloatingPointType<int>, "int ok");
 
constexpr double tolerance_equality_float = 1e-15;
 
template <typename T>
concept IntegerType = requires(T n) {
  requires std::is_integral_v<std::remove_cvref_t<T>>;
  requires !std::is_same_v<std::remove_cvref_t<T>, bool>;
  requires std::is_arithmetic_v<decltype(n + 1)>;
  requires !std::is_pointer_v<std::remove_cvref_t<T>>;
};
 
template <typename T>
concept NumberType = IntegerType<T> || FloatingPointType<T>;
 
// General case: by value
template <NumberType T>
inline bool almost_equal(T val, T val2, T tolerance = tolerance_equality_float)
{
  using CommonT = std::common_type_t<T, T>;
  return std::abs(static_cast<CommonT>(val) - static_cast<CommonT>(val2)) <
         static_cast<CommonT>(tolerance);
}
 
// Overload for references
template <NumberType T>
inline bool almost_equal(const T& val,
                         const T& val2,
                         T tolerance = tolerance_equality_float)
{
  using CommonT = std::common_type_t<T, T>;
  return std::abs(static_cast<CommonT>(val) - static_cast<CommonT>(val2)) <
         static_cast<CommonT>(tolerance);
}
 
// Overload for pointers
template <NumberType T>
inline bool almost_equal(const T* val,
                         const T* val2,
                         T tolerance = tolerance_equality_float)
{
  if (!val || !val2)
  {
    return false; // Null pointer check
  }
  using CommonT = std::common_type_t<T, T>;
  return std::abs(static_cast<CommonT>(*val) - static_cast<CommonT>(*val2)) <
         static_cast<CommonT>(tolerance);
}
 
// Overload for rvalue references
template <NumberType T>
inline bool
almost_equal(T&& val, T&& val2, T tolerance = tolerance_equality_float)
{
  using CommonT = std::common_type_t<T, T>;
  return std::abs(static_cast<CommonT>(val) - static_cast<CommonT>(val2)) <
         static_cast<CommonT>(tolerance);
}
 
#endif