1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
|
#pragma once
namespace rawaccel {
constexpr double minsd(double a, double b)
{
return (a < b) ? a : b;
}
constexpr double maxsd(double a, double b)
{
return (b < a) ? a : b;
}
constexpr double clampsd(double v, double lo, double hi)
{
return minsd(maxsd(v, lo), hi);
}
template <typename T>
constexpr const T& min(const T& a, const T& b)
{
return (b < a) ? b : a;
}
template <typename T>
constexpr const T& max(const T& a, const T& b)
{
return (b < a) ? a : b;
}
template <typename T>
constexpr const T& clamp(const T& v, const T& lo, const T& hi)
{
return (v < lo) ? lo : (hi < v) ? hi : v;
}
// returns the unbiased exponent of x if x is normal
inline int ilogb(double x)
{
union { double f; unsigned long long i; } u = { x };
return static_cast<int>((u.i >> 52) & 0x7ff) - 0x3ff;
}
// returns x * 2^n if n is in [-1022, 1023]
inline double scalbn(double x, int n)
{
union { double f; unsigned long long i; } u;
u.i = static_cast<unsigned long long>(0x3ff + n) << 52;
return x * u.f;
}
inline bool infnan(double x)
{
return ilogb(x) == 0x400;
}
struct noop {
template <typename... Ts>
constexpr void operator()(Ts&&...) const noexcept {}
};
template <typename T> struct remove_ref { using type = T; };
template <typename T> struct remove_ref<T&> { using type = T; };
template <typename T> struct remove_ref<T&&> { using type = T; };
template <typename T>
using remove_ref_t = typename remove_ref<T>::type;
template <typename T, typename U>
struct is_same { static constexpr bool value = false; };
template <typename T>
struct is_same<T, T> { static constexpr bool value = true; };
template <typename T, typename U>
inline constexpr bool is_same_v = is_same<T, U>::value;
template <class T> struct is_rvalue_ref { static constexpr bool value = false; };
template <class T> struct is_rvalue_ref<T&&> { static constexpr bool value = true; };
template <class T>
inline constexpr bool is_rvalue_ref_v = is_rvalue_ref<T>::value;
}
|