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
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
|
// Author: Martin C. Frith 2019
// SPDX-License-Identifier: MPL-2.0
#ifndef MCF_SIMD_HH
#define MCF_SIMD_HH
#define SIMDE_ENABLE_NATIVE_ALIASES
#include "simde/x86/avx2.h"
#include <stddef.h> // size_t
namespace mcf {
typedef __m256i SimdInt;
typedef __m256i SimdUint1;
typedef __m256d SimdDbl;
const int simdBytes = 32;
static inline SimdInt simdZero() {
return _mm256_setzero_si256();
}
static inline SimdInt simdZero1() {
return _mm256_setzero_si256();
}
static inline SimdDbl simdZeroDbl() {
return _mm256_setzero_pd();
}
static inline SimdInt simdOnes1() {
return _mm256_set1_epi32(-1);
}
static inline SimdInt simdLoad(const void *p) {
return _mm256_loadu_si256((const SimdInt *)p);
}
static inline SimdInt simdLoad1(const void *p) {
return _mm256_loadu_si256((const SimdInt *)p);
}
static inline SimdDbl simdLoadDbl(const double *p) {
return _mm256_loadu_pd(p);
}
static inline void simdStore(void *p, SimdInt x) {
_mm256_storeu_si256((SimdInt *)p, x);
}
static inline void simdStore1(void *p, SimdInt x) {
_mm256_storeu_si256((SimdInt *)p, x);
}
static inline void simdStoreDbl(double *p, SimdDbl x) {
_mm256_storeu_pd(p, x);
}
static inline SimdInt simdOr1(SimdInt x, SimdInt y) {
return _mm256_or_si256(x, y);
}
static inline SimdInt simdBlend(SimdInt x, SimdInt y, SimdInt mask) {
return _mm256_blendv_epi8(x, y, mask);
}
const int simdLen = 8;
const int simdDblLen = 4;
static inline SimdInt simdSet(int i7, int i6, int i5, int i4,
int i3, int i2, int i1, int i0) {
return _mm256_set_epi32(i7, i6, i5, i4, i3, i2, i1, i0);
}
static inline SimdInt simdSet1(char jF, char jE, char jD, char jC,
char jB, char jA, char j9, char j8,
char j7, char j6, char j5, char j4,
char j3, char j2, char j1, char j0,
char iF, char iE, char iD, char iC,
char iB, char iA, char i9, char i8,
char i7, char i6, char i5, char i4,
char i3, char i2, char i1, char i0) {
return _mm256_set_epi8(jF, jE, jD, jC, jB, jA, j9, j8,
j7, j6, j5, j4, j3, j2, j1, j0,
iF, iE, iD, iC, iB, iA, i9, i8,
i7, i6, i5, i4, i3, i2, i1, i0);
}
static inline SimdDbl simdSetDbl(double i3, double i2, double i1, double i0) {
return _mm256_set_pd(i3, i2, i1, i0);
}
static inline SimdInt simdFill(int x) {
return _mm256_set1_epi32(x);
}
static inline SimdInt simdFill1(char x) {
return _mm256_set1_epi8(x);
}
static inline SimdDbl simdFillDbl(double x) {
return _mm256_set1_pd(x);
}
static inline SimdInt simdGt(SimdInt x, SimdInt y) {
return _mm256_cmpgt_epi32(x, y);
}
static inline SimdInt simdGe1(SimdInt x, SimdInt y) {
return _mm256_cmpeq_epi8(_mm256_min_epu8(x, y), y);
}
static inline SimdInt simdAdd(SimdInt x, SimdInt y) {
return _mm256_add_epi32(x, y);
}
static inline SimdInt simdAdd1(SimdInt x, SimdInt y) {
return _mm256_add_epi8(x, y);
}
static inline SimdInt simdAdds1(SimdInt x, SimdInt y) {
return _mm256_adds_epu8(x, y);
}
static inline SimdDbl simdAddDbl(SimdDbl x, SimdDbl y) {
return _mm256_add_pd(x, y);
}
static inline SimdInt simdSub(SimdInt x, SimdInt y) {
return _mm256_sub_epi32(x, y);
}
static inline SimdInt simdSub1(SimdInt x, SimdInt y) {
return _mm256_sub_epi8(x, y);
}
static inline SimdDbl simdMulDbl(SimdDbl x, SimdDbl y) {
return _mm256_mul_pd(x, y);
}
static inline SimdInt simdQuadruple1(SimdInt x) {
return _mm256_slli_epi32(x, 2);
}
static inline SimdInt simdMax(SimdInt x, SimdInt y) {
return _mm256_max_epi32(x, y);
}
static inline SimdInt simdMin1(SimdInt x, SimdInt y) {
return _mm256_min_epu8(x, y);
}
static inline int simdHorizontalMax(SimdInt x) {
__m128i z = _mm256_castsi256_si128(x);
z = _mm_max_epi32(z, _mm256_extracti128_si256(x, 1));
z = _mm_max_epi32(z, _mm_shuffle_epi32(z, 0x4E));
z = _mm_max_epi32(z, _mm_shuffle_epi32(z, 0xB1));
return _mm_cvtsi128_si32(z);
}
static inline int simdHorizontalMin1(SimdInt x) {
__m128i z = _mm256_castsi256_si128(x);
z = _mm_min_epu8(z, _mm256_extracti128_si256(x, 1));
z = _mm_min_epu8(z, _mm_srli_epi16(z, 8));
z = _mm_minpos_epu16(z);
return _mm_extract_epi16(z, 0);
}
static inline double simdHorizontalAddDbl(SimdDbl x) {
__m128d z = _mm256_castpd256_pd128(x);
z = _mm_add_pd(z, _mm256_extractf128_pd(x, 1));
return _mm_cvtsd_f64(_mm_hadd_pd(z, z));
}
static inline SimdInt simdChoose1(SimdInt items, SimdInt choices) {
return _mm256_shuffle_epi8(items, choices);
}
}
#endif
|
