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
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
|
// -*-C++-*-
#ifndef MATHFUNCS_FABS_H
#define MATHFUNCS_FABS_H
#include "mathfuncs_base.h"
#include <cassert>
#include <cmath>
namespace vecmathlib {
template<typename realvec_t>
realvec_t mathfuncs<realvec_t>::vml_copysign(realvec_t x, realvec_t y)
{
intvec_t value = as_int(x) & IV(~FP::signbit_mask);
intvec_t sign = as_int(y) & IV(FP::signbit_mask);
return as_float(sign | value);
}
template<typename realvec_t>
realvec_t mathfuncs<realvec_t>::vml_fabs(realvec_t x)
{
return as_float(as_int(x) & IV(~FP::signbit_mask));
}
template<typename realvec_t>
realvec_t mathfuncs<realvec_t>::vml_fdim(realvec_t x, realvec_t y)
{
// return ifthen(x > y, x - y, RV(0.0));
return x - fmin(x, y);
}
template<typename realvec_t>
realvec_t mathfuncs<realvec_t>::vml_fma(realvec_t x, realvec_t y, realvec_t z)
{
return x * y + z;
}
template<typename realvec_t>
realvec_t mathfuncs<realvec_t>::vml_fmax(realvec_t x, realvec_t y)
{
return ifthen(x < y, y, x);
}
template<typename realvec_t>
realvec_t mathfuncs<realvec_t>::vml_fmin(realvec_t x, realvec_t y)
{
return ifthen(y < x, y, x);
}
template<typename realvec_t>
auto mathfuncs<realvec_t>::vml_ilogb(realvec_t x) -> intvec_t
{
return
lsr(as_int(x) & IV(FP::exponent_mask), FP::mantissa_bits) -
IV(FP::exponent_offset);
}
template<typename realvec_t>
auto mathfuncs<realvec_t>::vml_isfinite(realvec_t x) -> boolvec_t
{
#if defined VML_HAVE_INF || defined VML_HAVE_NAN
return (as_int(x) & IV(FP::exponent_mask)) != IV(FP::exponent_mask);
#else
return BV(true);
#endif
}
template<typename realvec_t>
auto mathfuncs<realvec_t>::vml_isinf(realvec_t x) -> boolvec_t
{
#if defined VML_HAVE_INF
return (as_int(x) & IV(~FP::signbit_mask)) == IV(FP::exponent_mask);
#else
return BV(false);
#endif
}
template<typename realvec_t>
auto mathfuncs<realvec_t>::vml_isnan(realvec_t x) -> boolvec_t
{
#if defined VML_HAVE_NAN
return
(as_int(x) & IV(FP::exponent_mask)) == IV(FP::exponent_mask) &&
(as_int(x) & IV(FP::mantissa_mask)) != IV(I(0));
#else
return BV(false);
#endif
}
template<typename realvec_t>
auto mathfuncs<realvec_t>::vml_isnormal(realvec_t x) -> boolvec_t
{
#if defined VML_HAVE_DENORMALS || defined VML_HAVE_INF || defined VML_HAVE_NAN
return
(as_int(x) & IV(FP::exponent_mask)) != IV(FP::exponent_mask) &&
(as_int(x) & IV(FP::exponent_mask)) != IV(I(0));
#else
return BV(true);
#endif
}
template<typename realvec_t>
realvec_t mathfuncs<realvec_t>::vml_scalbn(realvec_t x, intvec_t n)
{
return as_float(as_int(x) + (n << I(FP::mantissa_bits)));
// return x * as_float((n + exponent_offset) << mantissa_bits);
}
template<typename realvec_t>
auto mathfuncs<realvec_t>::vml_signbit(realvec_t x) -> boolvec_t
{
return convert_bool(as_int(x) & IV(FP::signbit_mask));
}
}; // namespace vecmathlib
#endif // #ifndef MATHFUNCS_FABS_H
|