blob: d7b63cfd7d4512f5c3c127e4f99237e3c05b912d (
plain) (
blame)
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
|
#pragma once
#ifdef _MANAGED
#include <math.h>
inline double sqrtsd(double val) { return sqrt(val); }
#else
#include <emmintrin.h>
inline double sqrtsd(double val) {
__m128d src = _mm_load_sd(&val);
__m128d dst = _mm_sqrt_sd(src, src);
_mm_store_sd(&val, dst);
return val;
}
#endif
inline constexpr double minsd(double a, double b) {
return (a < b) ? a : b;
}
inline constexpr double maxsd(double a, double b) {
return (b < a) ? a : b;
}
inline constexpr double clampsd(double v, double lo, double hi) {
return minsd(maxsd(v, lo), hi);
}
// 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;
}
inline bool infnan(double x)
{
return ilogb(x) == 0x400;
}
|