1 /*===---- avxintrin.h - AVX intrinsics -------------------------------------===
2 *
3 * Permission is hereby granted, free of charge, to any person obtaining a copy
4 * of this software and associated documentation files (the "Software"), to deal
5 * in the Software without restriction, including without limitation the rights
6 * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
7 * copies of the Software, and to permit persons to whom the Software is
8 * furnished to do so, subject to the following conditions:
9 *
10 * The above copyright notice and this permission notice shall be included in
11 * all copies or substantial portions of the Software.
12 *
13 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
14 * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
15 * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
16 * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
17 * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
18 * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
19 * THE SOFTWARE.
20 *
21 *===-----------------------------------------------------------------------===
22 */
23
24 #ifndef __IMMINTRIN_H
25 #error "Never use <avxintrin.h> directly; include <immintrin.h> instead."
26 #endif
27
28 #ifndef __AVXINTRIN_H
29 #define __AVXINTRIN_H
30
31 typedef double __v4df __attribute__ ((__vector_size__ (32)));
32 typedef float __v8sf __attribute__ ((__vector_size__ (32)));
33 typedef long long __v4di __attribute__ ((__vector_size__ (32)));
34 typedef int __v8si __attribute__ ((__vector_size__ (32)));
35 typedef short __v16hi __attribute__ ((__vector_size__ (32)));
36 typedef char __v32qi __attribute__ ((__vector_size__ (32)));
37
38 /* We need an explicitly signed variant for char. Note that this shouldn't
39 * appear in the interface though. */
40 typedef signed char __v32qs __attribute__((__vector_size__(32)));
41
42 typedef float __m256 __attribute__ ((__vector_size__ (32)));
43 typedef double __m256d __attribute__((__vector_size__(32)));
44 typedef long long __m256i __attribute__((__vector_size__(32)));
45
46 /* Define the default attributes for the functions in this file. */
47 #define __DEFAULT_FN_ATTRS __attribute__((__always_inline__, __nodebug__, __target__("avx")))
48
49 /* Arithmetic */
50 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_add_pd(__m256d __a,__m256d __b)51 _mm256_add_pd(__m256d __a, __m256d __b)
52 {
53 return __a+__b;
54 }
55
56 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_add_ps(__m256 __a,__m256 __b)57 _mm256_add_ps(__m256 __a, __m256 __b)
58 {
59 return __a+__b;
60 }
61
62 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_sub_pd(__m256d __a,__m256d __b)63 _mm256_sub_pd(__m256d __a, __m256d __b)
64 {
65 return __a-__b;
66 }
67
68 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_sub_ps(__m256 __a,__m256 __b)69 _mm256_sub_ps(__m256 __a, __m256 __b)
70 {
71 return __a-__b;
72 }
73
74 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_addsub_pd(__m256d __a,__m256d __b)75 _mm256_addsub_pd(__m256d __a, __m256d __b)
76 {
77 return (__m256d)__builtin_ia32_addsubpd256((__v4df)__a, (__v4df)__b);
78 }
79
80 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_addsub_ps(__m256 __a,__m256 __b)81 _mm256_addsub_ps(__m256 __a, __m256 __b)
82 {
83 return (__m256)__builtin_ia32_addsubps256((__v8sf)__a, (__v8sf)__b);
84 }
85
86 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_div_pd(__m256d __a,__m256d __b)87 _mm256_div_pd(__m256d __a, __m256d __b)
88 {
89 return __a / __b;
90 }
91
92 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_div_ps(__m256 __a,__m256 __b)93 _mm256_div_ps(__m256 __a, __m256 __b)
94 {
95 return __a / __b;
96 }
97
98 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_max_pd(__m256d __a,__m256d __b)99 _mm256_max_pd(__m256d __a, __m256d __b)
100 {
101 return (__m256d)__builtin_ia32_maxpd256((__v4df)__a, (__v4df)__b);
102 }
103
104 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_max_ps(__m256 __a,__m256 __b)105 _mm256_max_ps(__m256 __a, __m256 __b)
106 {
107 return (__m256)__builtin_ia32_maxps256((__v8sf)__a, (__v8sf)__b);
108 }
109
110 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_min_pd(__m256d __a,__m256d __b)111 _mm256_min_pd(__m256d __a, __m256d __b)
112 {
113 return (__m256d)__builtin_ia32_minpd256((__v4df)__a, (__v4df)__b);
114 }
115
116 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_min_ps(__m256 __a,__m256 __b)117 _mm256_min_ps(__m256 __a, __m256 __b)
118 {
119 return (__m256)__builtin_ia32_minps256((__v8sf)__a, (__v8sf)__b);
120 }
121
122 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_mul_pd(__m256d __a,__m256d __b)123 _mm256_mul_pd(__m256d __a, __m256d __b)
124 {
125 return __a * __b;
126 }
127
128 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_mul_ps(__m256 __a,__m256 __b)129 _mm256_mul_ps(__m256 __a, __m256 __b)
130 {
131 return __a * __b;
132 }
133
134 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_sqrt_pd(__m256d __a)135 _mm256_sqrt_pd(__m256d __a)
136 {
137 return (__m256d)__builtin_ia32_sqrtpd256((__v4df)__a);
138 }
139
140 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_sqrt_ps(__m256 __a)141 _mm256_sqrt_ps(__m256 __a)
142 {
143 return (__m256)__builtin_ia32_sqrtps256((__v8sf)__a);
144 }
145
146 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_rsqrt_ps(__m256 __a)147 _mm256_rsqrt_ps(__m256 __a)
148 {
149 return (__m256)__builtin_ia32_rsqrtps256((__v8sf)__a);
150 }
151
152 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_rcp_ps(__m256 __a)153 _mm256_rcp_ps(__m256 __a)
154 {
155 return (__m256)__builtin_ia32_rcpps256((__v8sf)__a);
156 }
157
158 #define _mm256_round_pd(V, M) __extension__ ({ \
159 (__m256d)__builtin_ia32_roundpd256((__v4df)(__m256d)(V), (M)); })
160
161 #define _mm256_round_ps(V, M) __extension__ ({ \
162 (__m256)__builtin_ia32_roundps256((__v8sf)(__m256)(V), (M)); })
163
164 #define _mm256_ceil_pd(V) _mm256_round_pd((V), _MM_FROUND_CEIL)
165 #define _mm256_floor_pd(V) _mm256_round_pd((V), _MM_FROUND_FLOOR)
166 #define _mm256_ceil_ps(V) _mm256_round_ps((V), _MM_FROUND_CEIL)
167 #define _mm256_floor_ps(V) _mm256_round_ps((V), _MM_FROUND_FLOOR)
168
169 /* Logical */
170 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_and_pd(__m256d __a,__m256d __b)171 _mm256_and_pd(__m256d __a, __m256d __b)
172 {
173 return (__m256d)((__v4di)__a & (__v4di)__b);
174 }
175
176 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_and_ps(__m256 __a,__m256 __b)177 _mm256_and_ps(__m256 __a, __m256 __b)
178 {
179 return (__m256)((__v8si)__a & (__v8si)__b);
180 }
181
182 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_andnot_pd(__m256d __a,__m256d __b)183 _mm256_andnot_pd(__m256d __a, __m256d __b)
184 {
185 return (__m256d)(~(__v4di)__a & (__v4di)__b);
186 }
187
188 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_andnot_ps(__m256 __a,__m256 __b)189 _mm256_andnot_ps(__m256 __a, __m256 __b)
190 {
191 return (__m256)(~(__v8si)__a & (__v8si)__b);
192 }
193
194 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_or_pd(__m256d __a,__m256d __b)195 _mm256_or_pd(__m256d __a, __m256d __b)
196 {
197 return (__m256d)((__v4di)__a | (__v4di)__b);
198 }
199
200 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_or_ps(__m256 __a,__m256 __b)201 _mm256_or_ps(__m256 __a, __m256 __b)
202 {
203 return (__m256)((__v8si)__a | (__v8si)__b);
204 }
205
206 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_xor_pd(__m256d __a,__m256d __b)207 _mm256_xor_pd(__m256d __a, __m256d __b)
208 {
209 return (__m256d)((__v4di)__a ^ (__v4di)__b);
210 }
211
212 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_xor_ps(__m256 __a,__m256 __b)213 _mm256_xor_ps(__m256 __a, __m256 __b)
214 {
215 return (__m256)((__v8si)__a ^ (__v8si)__b);
216 }
217
218 /* Horizontal arithmetic */
219 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_hadd_pd(__m256d __a,__m256d __b)220 _mm256_hadd_pd(__m256d __a, __m256d __b)
221 {
222 return (__m256d)__builtin_ia32_haddpd256((__v4df)__a, (__v4df)__b);
223 }
224
225 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_hadd_ps(__m256 __a,__m256 __b)226 _mm256_hadd_ps(__m256 __a, __m256 __b)
227 {
228 return (__m256)__builtin_ia32_haddps256((__v8sf)__a, (__v8sf)__b);
229 }
230
231 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_hsub_pd(__m256d __a,__m256d __b)232 _mm256_hsub_pd(__m256d __a, __m256d __b)
233 {
234 return (__m256d)__builtin_ia32_hsubpd256((__v4df)__a, (__v4df)__b);
235 }
236
237 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_hsub_ps(__m256 __a,__m256 __b)238 _mm256_hsub_ps(__m256 __a, __m256 __b)
239 {
240 return (__m256)__builtin_ia32_hsubps256((__v8sf)__a, (__v8sf)__b);
241 }
242
243 /* Vector permutations */
244 static __inline __m128d __DEFAULT_FN_ATTRS
_mm_permutevar_pd(__m128d __a,__m128i __c)245 _mm_permutevar_pd(__m128d __a, __m128i __c)
246 {
247 return (__m128d)__builtin_ia32_vpermilvarpd((__v2df)__a, (__v2di)__c);
248 }
249
250 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_permutevar_pd(__m256d __a,__m256i __c)251 _mm256_permutevar_pd(__m256d __a, __m256i __c)
252 {
253 return (__m256d)__builtin_ia32_vpermilvarpd256((__v4df)__a, (__v4di)__c);
254 }
255
256 static __inline __m128 __DEFAULT_FN_ATTRS
_mm_permutevar_ps(__m128 __a,__m128i __c)257 _mm_permutevar_ps(__m128 __a, __m128i __c)
258 {
259 return (__m128)__builtin_ia32_vpermilvarps((__v4sf)__a, (__v4si)__c);
260 }
261
262 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_permutevar_ps(__m256 __a,__m256i __c)263 _mm256_permutevar_ps(__m256 __a, __m256i __c)
264 {
265 return (__m256)__builtin_ia32_vpermilvarps256((__v8sf)__a, (__v8si)__c);
266 }
267
268 #define _mm_permute_pd(A, C) __extension__ ({ \
269 (__m128d)__builtin_shufflevector((__v2df)(__m128d)(A), \
270 (__v2df)_mm_setzero_pd(), \
271 (C) & 0x1, ((C) & 0x2) >> 1); })
272
273 #define _mm256_permute_pd(A, C) __extension__ ({ \
274 (__m256d)__builtin_shufflevector((__v4df)(__m256d)(A), \
275 (__v4df)_mm256_setzero_pd(), \
276 (C) & 0x1, ((C) & 0x2) >> 1, \
277 2 + (((C) & 0x4) >> 2), \
278 2 + (((C) & 0x8) >> 3)); })
279
280 #define _mm_permute_ps(A, C) __extension__ ({ \
281 (__m128)__builtin_shufflevector((__v4sf)(__m128)(A), \
282 (__v4sf)_mm_setzero_ps(), \
283 (C) & 0x3, ((C) & 0xc) >> 2, \
284 ((C) & 0x30) >> 4, ((C) & 0xc0) >> 6); })
285
286 #define _mm256_permute_ps(A, C) __extension__ ({ \
287 (__m256)__builtin_shufflevector((__v8sf)(__m256)(A), \
288 (__v8sf)_mm256_setzero_ps(), \
289 (C) & 0x3, ((C) & 0xc) >> 2, \
290 ((C) & 0x30) >> 4, ((C) & 0xc0) >> 6, \
291 4 + (((C) & 0x03) >> 0), \
292 4 + (((C) & 0x0c) >> 2), \
293 4 + (((C) & 0x30) >> 4), \
294 4 + (((C) & 0xc0) >> 6)); })
295
296 #define _mm256_permute2f128_pd(V1, V2, M) __extension__ ({ \
297 (__m256d)__builtin_ia32_vperm2f128_pd256((__v4df)(__m256d)(V1), \
298 (__v4df)(__m256d)(V2), (M)); })
299
300 #define _mm256_permute2f128_ps(V1, V2, M) __extension__ ({ \
301 (__m256)__builtin_ia32_vperm2f128_ps256((__v8sf)(__m256)(V1), \
302 (__v8sf)(__m256)(V2), (M)); })
303
304 #define _mm256_permute2f128_si256(V1, V2, M) __extension__ ({ \
305 (__m256i)__builtin_ia32_vperm2f128_si256((__v8si)(__m256i)(V1), \
306 (__v8si)(__m256i)(V2), (M)); })
307
308 /* Vector Blend */
309 #define _mm256_blend_pd(V1, V2, M) __extension__ ({ \
310 (__m256d)__builtin_shufflevector((__v4df)(__m256d)(V1), \
311 (__v4df)(__m256d)(V2), \
312 (((M) & 0x01) ? 4 : 0), \
313 (((M) & 0x02) ? 5 : 1), \
314 (((M) & 0x04) ? 6 : 2), \
315 (((M) & 0x08) ? 7 : 3)); })
316
317 #define _mm256_blend_ps(V1, V2, M) __extension__ ({ \
318 (__m256)__builtin_shufflevector((__v8sf)(__m256)(V1), \
319 (__v8sf)(__m256)(V2), \
320 (((M) & 0x01) ? 8 : 0), \
321 (((M) & 0x02) ? 9 : 1), \
322 (((M) & 0x04) ? 10 : 2), \
323 (((M) & 0x08) ? 11 : 3), \
324 (((M) & 0x10) ? 12 : 4), \
325 (((M) & 0x20) ? 13 : 5), \
326 (((M) & 0x40) ? 14 : 6), \
327 (((M) & 0x80) ? 15 : 7)); })
328
329 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_blendv_pd(__m256d __a,__m256d __b,__m256d __c)330 _mm256_blendv_pd(__m256d __a, __m256d __b, __m256d __c)
331 {
332 return (__m256d)__builtin_ia32_blendvpd256(
333 (__v4df)__a, (__v4df)__b, (__v4df)__c);
334 }
335
336 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_blendv_ps(__m256 __a,__m256 __b,__m256 __c)337 _mm256_blendv_ps(__m256 __a, __m256 __b, __m256 __c)
338 {
339 return (__m256)__builtin_ia32_blendvps256(
340 (__v8sf)__a, (__v8sf)__b, (__v8sf)__c);
341 }
342
343 /* Vector Dot Product */
344 #define _mm256_dp_ps(V1, V2, M) __extension__ ({ \
345 (__m256)__builtin_ia32_dpps256((__v8sf)(__m256)(V1), \
346 (__v8sf)(__m256)(V2), (M)); })
347
348 /* Vector shuffle */
349 #define _mm256_shuffle_ps(a, b, mask) __extension__ ({ \
350 (__m256)__builtin_shufflevector((__v8sf)(__m256)(a), \
351 (__v8sf)(__m256)(b), \
352 (mask) & 0x3, \
353 ((mask) & 0xc) >> 2, \
354 (((mask) & 0x30) >> 4) + 8, \
355 (((mask) & 0xc0) >> 6) + 8, \
356 ((mask) & 0x3) + 4, \
357 (((mask) & 0xc) >> 2) + 4, \
358 (((mask) & 0x30) >> 4) + 12, \
359 (((mask) & 0xc0) >> 6) + 12); })
360
361 #define _mm256_shuffle_pd(a, b, mask) __extension__ ({ \
362 (__m256d)__builtin_shufflevector((__v4df)(__m256d)(a), \
363 (__v4df)(__m256d)(b), \
364 (mask) & 0x1, \
365 (((mask) & 0x2) >> 1) + 4, \
366 (((mask) & 0x4) >> 2) + 2, \
367 (((mask) & 0x8) >> 3) + 6); })
368
369 /* Compare */
370 #define _CMP_EQ_OQ 0x00 /* Equal (ordered, non-signaling) */
371 #define _CMP_LT_OS 0x01 /* Less-than (ordered, signaling) */
372 #define _CMP_LE_OS 0x02 /* Less-than-or-equal (ordered, signaling) */
373 #define _CMP_UNORD_Q 0x03 /* Unordered (non-signaling) */
374 #define _CMP_NEQ_UQ 0x04 /* Not-equal (unordered, non-signaling) */
375 #define _CMP_NLT_US 0x05 /* Not-less-than (unordered, signaling) */
376 #define _CMP_NLE_US 0x06 /* Not-less-than-or-equal (unordered, signaling) */
377 #define _CMP_ORD_Q 0x07 /* Ordered (nonsignaling) */
378 #define _CMP_EQ_UQ 0x08 /* Equal (unordered, non-signaling) */
379 #define _CMP_NGE_US 0x09 /* Not-greater-than-or-equal (unord, signaling) */
380 #define _CMP_NGT_US 0x0a /* Not-greater-than (unordered, signaling) */
381 #define _CMP_FALSE_OQ 0x0b /* False (ordered, non-signaling) */
382 #define _CMP_NEQ_OQ 0x0c /* Not-equal (ordered, non-signaling) */
383 #define _CMP_GE_OS 0x0d /* Greater-than-or-equal (ordered, signaling) */
384 #define _CMP_GT_OS 0x0e /* Greater-than (ordered, signaling) */
385 #define _CMP_TRUE_UQ 0x0f /* True (unordered, non-signaling) */
386 #define _CMP_EQ_OS 0x10 /* Equal (ordered, signaling) */
387 #define _CMP_LT_OQ 0x11 /* Less-than (ordered, non-signaling) */
388 #define _CMP_LE_OQ 0x12 /* Less-than-or-equal (ordered, non-signaling) */
389 #define _CMP_UNORD_S 0x13 /* Unordered (signaling) */
390 #define _CMP_NEQ_US 0x14 /* Not-equal (unordered, signaling) */
391 #define _CMP_NLT_UQ 0x15 /* Not-less-than (unordered, non-signaling) */
392 #define _CMP_NLE_UQ 0x16 /* Not-less-than-or-equal (unord, non-signaling) */
393 #define _CMP_ORD_S 0x17 /* Ordered (signaling) */
394 #define _CMP_EQ_US 0x18 /* Equal (unordered, signaling) */
395 #define _CMP_NGE_UQ 0x19 /* Not-greater-than-or-equal (unord, non-sign) */
396 #define _CMP_NGT_UQ 0x1a /* Not-greater-than (unordered, non-signaling) */
397 #define _CMP_FALSE_OS 0x1b /* False (ordered, signaling) */
398 #define _CMP_NEQ_OS 0x1c /* Not-equal (ordered, signaling) */
399 #define _CMP_GE_OQ 0x1d /* Greater-than-or-equal (ordered, non-signaling) */
400 #define _CMP_GT_OQ 0x1e /* Greater-than (ordered, non-signaling) */
401 #define _CMP_TRUE_US 0x1f /* True (unordered, signaling) */
402
403 #define _mm_cmp_pd(a, b, c) __extension__ ({ \
404 (__m128d)__builtin_ia32_cmppd((__v2df)(__m128d)(a), \
405 (__v2df)(__m128d)(b), (c)); })
406
407 #define _mm_cmp_ps(a, b, c) __extension__ ({ \
408 (__m128)__builtin_ia32_cmpps((__v4sf)(__m128)(a), \
409 (__v4sf)(__m128)(b), (c)); })
410
411 #define _mm256_cmp_pd(a, b, c) __extension__ ({ \
412 (__m256d)__builtin_ia32_cmppd256((__v4df)(__m256d)(a), \
413 (__v4df)(__m256d)(b), (c)); })
414
415 #define _mm256_cmp_ps(a, b, c) __extension__ ({ \
416 (__m256)__builtin_ia32_cmpps256((__v8sf)(__m256)(a), \
417 (__v8sf)(__m256)(b), (c)); })
418
419 #define _mm_cmp_sd(a, b, c) __extension__ ({ \
420 (__m128d)__builtin_ia32_cmpsd((__v2df)(__m128d)(a), \
421 (__v2df)(__m128d)(b), (c)); })
422
423 #define _mm_cmp_ss(a, b, c) __extension__ ({ \
424 (__m128)__builtin_ia32_cmpss((__v4sf)(__m128)(a), \
425 (__v4sf)(__m128)(b), (c)); })
426
427 static __inline int __DEFAULT_FN_ATTRS
_mm256_extract_epi32(__m256i __a,const int __imm)428 _mm256_extract_epi32(__m256i __a, const int __imm)
429 {
430 __v8si __b = (__v8si)__a;
431 return __b[__imm & 7];
432 }
433
434 static __inline int __DEFAULT_FN_ATTRS
_mm256_extract_epi16(__m256i __a,const int __imm)435 _mm256_extract_epi16(__m256i __a, const int __imm)
436 {
437 __v16hi __b = (__v16hi)__a;
438 return __b[__imm & 15];
439 }
440
441 static __inline int __DEFAULT_FN_ATTRS
_mm256_extract_epi8(__m256i __a,const int __imm)442 _mm256_extract_epi8(__m256i __a, const int __imm)
443 {
444 __v32qi __b = (__v32qi)__a;
445 return __b[__imm & 31];
446 }
447
448 #ifdef __x86_64__
449 static __inline long long __DEFAULT_FN_ATTRS
_mm256_extract_epi64(__m256i __a,const int __imm)450 _mm256_extract_epi64(__m256i __a, const int __imm)
451 {
452 __v4di __b = (__v4di)__a;
453 return __b[__imm & 3];
454 }
455 #endif
456
457 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_insert_epi32(__m256i __a,int __b,int const __imm)458 _mm256_insert_epi32(__m256i __a, int __b, int const __imm)
459 {
460 __v8si __c = (__v8si)__a;
461 __c[__imm & 7] = __b;
462 return (__m256i)__c;
463 }
464
465 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_insert_epi16(__m256i __a,int __b,int const __imm)466 _mm256_insert_epi16(__m256i __a, int __b, int const __imm)
467 {
468 __v16hi __c = (__v16hi)__a;
469 __c[__imm & 15] = __b;
470 return (__m256i)__c;
471 }
472
473 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_insert_epi8(__m256i __a,int __b,int const __imm)474 _mm256_insert_epi8(__m256i __a, int __b, int const __imm)
475 {
476 __v32qi __c = (__v32qi)__a;
477 __c[__imm & 31] = __b;
478 return (__m256i)__c;
479 }
480
481 #ifdef __x86_64__
482 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_insert_epi64(__m256i __a,long long __b,int const __imm)483 _mm256_insert_epi64(__m256i __a, long long __b, int const __imm)
484 {
485 __v4di __c = (__v4di)__a;
486 __c[__imm & 3] = __b;
487 return (__m256i)__c;
488 }
489 #endif
490
491 /* Conversion */
492 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_cvtepi32_pd(__m128i __a)493 _mm256_cvtepi32_pd(__m128i __a)
494 {
495 return (__m256d)__builtin_ia32_cvtdq2pd256((__v4si) __a);
496 }
497
498 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_cvtepi32_ps(__m256i __a)499 _mm256_cvtepi32_ps(__m256i __a)
500 {
501 return (__m256)__builtin_ia32_cvtdq2ps256((__v8si) __a);
502 }
503
504 static __inline __m128 __DEFAULT_FN_ATTRS
_mm256_cvtpd_ps(__m256d __a)505 _mm256_cvtpd_ps(__m256d __a)
506 {
507 return (__m128)__builtin_ia32_cvtpd2ps256((__v4df) __a);
508 }
509
510 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_cvtps_epi32(__m256 __a)511 _mm256_cvtps_epi32(__m256 __a)
512 {
513 return (__m256i)__builtin_ia32_cvtps2dq256((__v8sf) __a);
514 }
515
516 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_cvtps_pd(__m128 __a)517 _mm256_cvtps_pd(__m128 __a)
518 {
519 return (__m256d)__builtin_ia32_cvtps2pd256((__v4sf) __a);
520 }
521
522 static __inline __m128i __DEFAULT_FN_ATTRS
_mm256_cvttpd_epi32(__m256d __a)523 _mm256_cvttpd_epi32(__m256d __a)
524 {
525 return (__m128i)__builtin_ia32_cvttpd2dq256((__v4df) __a);
526 }
527
528 static __inline __m128i __DEFAULT_FN_ATTRS
_mm256_cvtpd_epi32(__m256d __a)529 _mm256_cvtpd_epi32(__m256d __a)
530 {
531 return (__m128i)__builtin_ia32_cvtpd2dq256((__v4df) __a);
532 }
533
534 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_cvttps_epi32(__m256 __a)535 _mm256_cvttps_epi32(__m256 __a)
536 {
537 return (__m256i)__builtin_ia32_cvttps2dq256((__v8sf) __a);
538 }
539
540 /* Vector replicate */
541 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_movehdup_ps(__m256 __a)542 _mm256_movehdup_ps(__m256 __a)
543 {
544 return __builtin_shufflevector(__a, __a, 1, 1, 3, 3, 5, 5, 7, 7);
545 }
546
547 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_moveldup_ps(__m256 __a)548 _mm256_moveldup_ps(__m256 __a)
549 {
550 return __builtin_shufflevector(__a, __a, 0, 0, 2, 2, 4, 4, 6, 6);
551 }
552
553 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_movedup_pd(__m256d __a)554 _mm256_movedup_pd(__m256d __a)
555 {
556 return __builtin_shufflevector(__a, __a, 0, 0, 2, 2);
557 }
558
559 /* Unpack and Interleave */
560 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_unpackhi_pd(__m256d __a,__m256d __b)561 _mm256_unpackhi_pd(__m256d __a, __m256d __b)
562 {
563 return __builtin_shufflevector(__a, __b, 1, 5, 1+2, 5+2);
564 }
565
566 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_unpacklo_pd(__m256d __a,__m256d __b)567 _mm256_unpacklo_pd(__m256d __a, __m256d __b)
568 {
569 return __builtin_shufflevector(__a, __b, 0, 4, 0+2, 4+2);
570 }
571
572 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_unpackhi_ps(__m256 __a,__m256 __b)573 _mm256_unpackhi_ps(__m256 __a, __m256 __b)
574 {
575 return __builtin_shufflevector(__a, __b, 2, 10, 2+1, 10+1, 6, 14, 6+1, 14+1);
576 }
577
578 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_unpacklo_ps(__m256 __a,__m256 __b)579 _mm256_unpacklo_ps(__m256 __a, __m256 __b)
580 {
581 return __builtin_shufflevector(__a, __b, 0, 8, 0+1, 8+1, 4, 12, 4+1, 12+1);
582 }
583
584 /* Bit Test */
585 static __inline int __DEFAULT_FN_ATTRS
_mm_testz_pd(__m128d __a,__m128d __b)586 _mm_testz_pd(__m128d __a, __m128d __b)
587 {
588 return __builtin_ia32_vtestzpd((__v2df)__a, (__v2df)__b);
589 }
590
591 static __inline int __DEFAULT_FN_ATTRS
_mm_testc_pd(__m128d __a,__m128d __b)592 _mm_testc_pd(__m128d __a, __m128d __b)
593 {
594 return __builtin_ia32_vtestcpd((__v2df)__a, (__v2df)__b);
595 }
596
597 static __inline int __DEFAULT_FN_ATTRS
_mm_testnzc_pd(__m128d __a,__m128d __b)598 _mm_testnzc_pd(__m128d __a, __m128d __b)
599 {
600 return __builtin_ia32_vtestnzcpd((__v2df)__a, (__v2df)__b);
601 }
602
603 static __inline int __DEFAULT_FN_ATTRS
_mm_testz_ps(__m128 __a,__m128 __b)604 _mm_testz_ps(__m128 __a, __m128 __b)
605 {
606 return __builtin_ia32_vtestzps((__v4sf)__a, (__v4sf)__b);
607 }
608
609 static __inline int __DEFAULT_FN_ATTRS
_mm_testc_ps(__m128 __a,__m128 __b)610 _mm_testc_ps(__m128 __a, __m128 __b)
611 {
612 return __builtin_ia32_vtestcps((__v4sf)__a, (__v4sf)__b);
613 }
614
615 static __inline int __DEFAULT_FN_ATTRS
_mm_testnzc_ps(__m128 __a,__m128 __b)616 _mm_testnzc_ps(__m128 __a, __m128 __b)
617 {
618 return __builtin_ia32_vtestnzcps((__v4sf)__a, (__v4sf)__b);
619 }
620
621 static __inline int __DEFAULT_FN_ATTRS
_mm256_testz_pd(__m256d __a,__m256d __b)622 _mm256_testz_pd(__m256d __a, __m256d __b)
623 {
624 return __builtin_ia32_vtestzpd256((__v4df)__a, (__v4df)__b);
625 }
626
627 static __inline int __DEFAULT_FN_ATTRS
_mm256_testc_pd(__m256d __a,__m256d __b)628 _mm256_testc_pd(__m256d __a, __m256d __b)
629 {
630 return __builtin_ia32_vtestcpd256((__v4df)__a, (__v4df)__b);
631 }
632
633 static __inline int __DEFAULT_FN_ATTRS
_mm256_testnzc_pd(__m256d __a,__m256d __b)634 _mm256_testnzc_pd(__m256d __a, __m256d __b)
635 {
636 return __builtin_ia32_vtestnzcpd256((__v4df)__a, (__v4df)__b);
637 }
638
639 static __inline int __DEFAULT_FN_ATTRS
_mm256_testz_ps(__m256 __a,__m256 __b)640 _mm256_testz_ps(__m256 __a, __m256 __b)
641 {
642 return __builtin_ia32_vtestzps256((__v8sf)__a, (__v8sf)__b);
643 }
644
645 static __inline int __DEFAULT_FN_ATTRS
_mm256_testc_ps(__m256 __a,__m256 __b)646 _mm256_testc_ps(__m256 __a, __m256 __b)
647 {
648 return __builtin_ia32_vtestcps256((__v8sf)__a, (__v8sf)__b);
649 }
650
651 static __inline int __DEFAULT_FN_ATTRS
_mm256_testnzc_ps(__m256 __a,__m256 __b)652 _mm256_testnzc_ps(__m256 __a, __m256 __b)
653 {
654 return __builtin_ia32_vtestnzcps256((__v8sf)__a, (__v8sf)__b);
655 }
656
657 static __inline int __DEFAULT_FN_ATTRS
_mm256_testz_si256(__m256i __a,__m256i __b)658 _mm256_testz_si256(__m256i __a, __m256i __b)
659 {
660 return __builtin_ia32_ptestz256((__v4di)__a, (__v4di)__b);
661 }
662
663 static __inline int __DEFAULT_FN_ATTRS
_mm256_testc_si256(__m256i __a,__m256i __b)664 _mm256_testc_si256(__m256i __a, __m256i __b)
665 {
666 return __builtin_ia32_ptestc256((__v4di)__a, (__v4di)__b);
667 }
668
669 static __inline int __DEFAULT_FN_ATTRS
_mm256_testnzc_si256(__m256i __a,__m256i __b)670 _mm256_testnzc_si256(__m256i __a, __m256i __b)
671 {
672 return __builtin_ia32_ptestnzc256((__v4di)__a, (__v4di)__b);
673 }
674
675 /* Vector extract sign mask */
676 static __inline int __DEFAULT_FN_ATTRS
_mm256_movemask_pd(__m256d __a)677 _mm256_movemask_pd(__m256d __a)
678 {
679 return __builtin_ia32_movmskpd256((__v4df)__a);
680 }
681
682 static __inline int __DEFAULT_FN_ATTRS
_mm256_movemask_ps(__m256 __a)683 _mm256_movemask_ps(__m256 __a)
684 {
685 return __builtin_ia32_movmskps256((__v8sf)__a);
686 }
687
688 /* Vector __zero */
689 static __inline void __DEFAULT_FN_ATTRS
_mm256_zeroall(void)690 _mm256_zeroall(void)
691 {
692 __builtin_ia32_vzeroall();
693 }
694
695 static __inline void __DEFAULT_FN_ATTRS
_mm256_zeroupper(void)696 _mm256_zeroupper(void)
697 {
698 __builtin_ia32_vzeroupper();
699 }
700
701 /* Vector load with broadcast */
702 static __inline __m128 __DEFAULT_FN_ATTRS
_mm_broadcast_ss(float const * __a)703 _mm_broadcast_ss(float const *__a)
704 {
705 float __f = *__a;
706 return (__m128)(__v4sf){ __f, __f, __f, __f };
707 }
708
709 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_broadcast_sd(double const * __a)710 _mm256_broadcast_sd(double const *__a)
711 {
712 double __d = *__a;
713 return (__m256d)(__v4df){ __d, __d, __d, __d };
714 }
715
716 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_broadcast_ss(float const * __a)717 _mm256_broadcast_ss(float const *__a)
718 {
719 float __f = *__a;
720 return (__m256)(__v8sf){ __f, __f, __f, __f, __f, __f, __f, __f };
721 }
722
723 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_broadcast_pd(__m128d const * __a)724 _mm256_broadcast_pd(__m128d const *__a)
725 {
726 return (__m256d)__builtin_ia32_vbroadcastf128_pd256(__a);
727 }
728
729 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_broadcast_ps(__m128 const * __a)730 _mm256_broadcast_ps(__m128 const *__a)
731 {
732 return (__m256)__builtin_ia32_vbroadcastf128_ps256(__a);
733 }
734
735 /* SIMD load ops */
736 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_load_pd(double const * __p)737 _mm256_load_pd(double const *__p)
738 {
739 return *(__m256d *)__p;
740 }
741
742 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_load_ps(float const * __p)743 _mm256_load_ps(float const *__p)
744 {
745 return *(__m256 *)__p;
746 }
747
748 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_loadu_pd(double const * __p)749 _mm256_loadu_pd(double const *__p)
750 {
751 struct __loadu_pd {
752 __m256d __v;
753 } __attribute__((__packed__, __may_alias__));
754 return ((struct __loadu_pd*)__p)->__v;
755 }
756
757 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_loadu_ps(float const * __p)758 _mm256_loadu_ps(float const *__p)
759 {
760 struct __loadu_ps {
761 __m256 __v;
762 } __attribute__((__packed__, __may_alias__));
763 return ((struct __loadu_ps*)__p)->__v;
764 }
765
766 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_load_si256(__m256i const * __p)767 _mm256_load_si256(__m256i const *__p)
768 {
769 return *__p;
770 }
771
772 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_loadu_si256(__m256i const * __p)773 _mm256_loadu_si256(__m256i const *__p)
774 {
775 struct __loadu_si256 {
776 __m256i __v;
777 } __attribute__((__packed__, __may_alias__));
778 return ((struct __loadu_si256*)__p)->__v;
779 }
780
781 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_lddqu_si256(__m256i const * __p)782 _mm256_lddqu_si256(__m256i const *__p)
783 {
784 return (__m256i)__builtin_ia32_lddqu256((char const *)__p);
785 }
786
787 /* SIMD store ops */
788 static __inline void __DEFAULT_FN_ATTRS
_mm256_store_pd(double * __p,__m256d __a)789 _mm256_store_pd(double *__p, __m256d __a)
790 {
791 *(__m256d *)__p = __a;
792 }
793
794 static __inline void __DEFAULT_FN_ATTRS
_mm256_store_ps(float * __p,__m256 __a)795 _mm256_store_ps(float *__p, __m256 __a)
796 {
797 *(__m256 *)__p = __a;
798 }
799
800 static __inline void __DEFAULT_FN_ATTRS
_mm256_storeu_pd(double * __p,__m256d __a)801 _mm256_storeu_pd(double *__p, __m256d __a)
802 {
803 __builtin_ia32_storeupd256(__p, (__v4df)__a);
804 }
805
806 static __inline void __DEFAULT_FN_ATTRS
_mm256_storeu_ps(float * __p,__m256 __a)807 _mm256_storeu_ps(float *__p, __m256 __a)
808 {
809 __builtin_ia32_storeups256(__p, (__v8sf)__a);
810 }
811
812 static __inline void __DEFAULT_FN_ATTRS
_mm256_store_si256(__m256i * __p,__m256i __a)813 _mm256_store_si256(__m256i *__p, __m256i __a)
814 {
815 *__p = __a;
816 }
817
818 static __inline void __DEFAULT_FN_ATTRS
_mm256_storeu_si256(__m256i * __p,__m256i __a)819 _mm256_storeu_si256(__m256i *__p, __m256i __a)
820 {
821 __builtin_ia32_storedqu256((char *)__p, (__v32qi)__a);
822 }
823
824 /* Conditional load ops */
825 static __inline __m128d __DEFAULT_FN_ATTRS
_mm_maskload_pd(double const * __p,__m128i __m)826 _mm_maskload_pd(double const *__p, __m128i __m)
827 {
828 return (__m128d)__builtin_ia32_maskloadpd((const __v2df *)__p, (__v2di)__m);
829 }
830
831 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_maskload_pd(double const * __p,__m256i __m)832 _mm256_maskload_pd(double const *__p, __m256i __m)
833 {
834 return (__m256d)__builtin_ia32_maskloadpd256((const __v4df *)__p,
835 (__v4di)__m);
836 }
837
838 static __inline __m128 __DEFAULT_FN_ATTRS
_mm_maskload_ps(float const * __p,__m128i __m)839 _mm_maskload_ps(float const *__p, __m128i __m)
840 {
841 return (__m128)__builtin_ia32_maskloadps((const __v4sf *)__p, (__v4si)__m);
842 }
843
844 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_maskload_ps(float const * __p,__m256i __m)845 _mm256_maskload_ps(float const *__p, __m256i __m)
846 {
847 return (__m256)__builtin_ia32_maskloadps256((const __v8sf *)__p, (__v8si)__m);
848 }
849
850 /* Conditional store ops */
851 static __inline void __DEFAULT_FN_ATTRS
_mm256_maskstore_ps(float * __p,__m256i __m,__m256 __a)852 _mm256_maskstore_ps(float *__p, __m256i __m, __m256 __a)
853 {
854 __builtin_ia32_maskstoreps256((__v8sf *)__p, (__v8si)__m, (__v8sf)__a);
855 }
856
857 static __inline void __DEFAULT_FN_ATTRS
_mm_maskstore_pd(double * __p,__m128i __m,__m128d __a)858 _mm_maskstore_pd(double *__p, __m128i __m, __m128d __a)
859 {
860 __builtin_ia32_maskstorepd((__v2df *)__p, (__v2di)__m, (__v2df)__a);
861 }
862
863 static __inline void __DEFAULT_FN_ATTRS
_mm256_maskstore_pd(double * __p,__m256i __m,__m256d __a)864 _mm256_maskstore_pd(double *__p, __m256i __m, __m256d __a)
865 {
866 __builtin_ia32_maskstorepd256((__v4df *)__p, (__v4di)__m, (__v4df)__a);
867 }
868
869 static __inline void __DEFAULT_FN_ATTRS
_mm_maskstore_ps(float * __p,__m128i __m,__m128 __a)870 _mm_maskstore_ps(float *__p, __m128i __m, __m128 __a)
871 {
872 __builtin_ia32_maskstoreps((__v4sf *)__p, (__v4si)__m, (__v4sf)__a);
873 }
874
875 /* Cacheability support ops */
876 static __inline void __DEFAULT_FN_ATTRS
_mm256_stream_si256(__m256i * __a,__m256i __b)877 _mm256_stream_si256(__m256i *__a, __m256i __b)
878 {
879 __builtin_ia32_movntdq256((__v4di *)__a, (__v4di)__b);
880 }
881
882 static __inline void __DEFAULT_FN_ATTRS
_mm256_stream_pd(double * __a,__m256d __b)883 _mm256_stream_pd(double *__a, __m256d __b)
884 {
885 __builtin_ia32_movntpd256(__a, (__v4df)__b);
886 }
887
888 static __inline void __DEFAULT_FN_ATTRS
_mm256_stream_ps(float * __p,__m256 __a)889 _mm256_stream_ps(float *__p, __m256 __a)
890 {
891 __builtin_ia32_movntps256(__p, (__v8sf)__a);
892 }
893
894 /* Create vectors */
895 static __inline__ __m256d __DEFAULT_FN_ATTRS
_mm256_undefined_pd()896 _mm256_undefined_pd()
897 {
898 return (__m256d)__builtin_ia32_undef256();
899 }
900
901 static __inline__ __m256 __DEFAULT_FN_ATTRS
_mm256_undefined_ps()902 _mm256_undefined_ps()
903 {
904 return (__m256)__builtin_ia32_undef256();
905 }
906
907 static __inline__ __m256i __DEFAULT_FN_ATTRS
_mm256_undefined_si256()908 _mm256_undefined_si256()
909 {
910 return (__m256i)__builtin_ia32_undef256();
911 }
912
913 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_set_pd(double __a,double __b,double __c,double __d)914 _mm256_set_pd(double __a, double __b, double __c, double __d)
915 {
916 return (__m256d){ __d, __c, __b, __a };
917 }
918
919 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_set_ps(float __a,float __b,float __c,float __d,float __e,float __f,float __g,float __h)920 _mm256_set_ps(float __a, float __b, float __c, float __d,
921 float __e, float __f, float __g, float __h)
922 {
923 return (__m256){ __h, __g, __f, __e, __d, __c, __b, __a };
924 }
925
926 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_set_epi32(int __i0,int __i1,int __i2,int __i3,int __i4,int __i5,int __i6,int __i7)927 _mm256_set_epi32(int __i0, int __i1, int __i2, int __i3,
928 int __i4, int __i5, int __i6, int __i7)
929 {
930 return (__m256i)(__v8si){ __i7, __i6, __i5, __i4, __i3, __i2, __i1, __i0 };
931 }
932
933 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_set_epi16(short __w15,short __w14,short __w13,short __w12,short __w11,short __w10,short __w09,short __w08,short __w07,short __w06,short __w05,short __w04,short __w03,short __w02,short __w01,short __w00)934 _mm256_set_epi16(short __w15, short __w14, short __w13, short __w12,
935 short __w11, short __w10, short __w09, short __w08,
936 short __w07, short __w06, short __w05, short __w04,
937 short __w03, short __w02, short __w01, short __w00)
938 {
939 return (__m256i)(__v16hi){ __w00, __w01, __w02, __w03, __w04, __w05, __w06,
940 __w07, __w08, __w09, __w10, __w11, __w12, __w13, __w14, __w15 };
941 }
942
943 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_set_epi8(char __b31,char __b30,char __b29,char __b28,char __b27,char __b26,char __b25,char __b24,char __b23,char __b22,char __b21,char __b20,char __b19,char __b18,char __b17,char __b16,char __b15,char __b14,char __b13,char __b12,char __b11,char __b10,char __b09,char __b08,char __b07,char __b06,char __b05,char __b04,char __b03,char __b02,char __b01,char __b00)944 _mm256_set_epi8(char __b31, char __b30, char __b29, char __b28,
945 char __b27, char __b26, char __b25, char __b24,
946 char __b23, char __b22, char __b21, char __b20,
947 char __b19, char __b18, char __b17, char __b16,
948 char __b15, char __b14, char __b13, char __b12,
949 char __b11, char __b10, char __b09, char __b08,
950 char __b07, char __b06, char __b05, char __b04,
951 char __b03, char __b02, char __b01, char __b00)
952 {
953 return (__m256i)(__v32qi){
954 __b00, __b01, __b02, __b03, __b04, __b05, __b06, __b07,
955 __b08, __b09, __b10, __b11, __b12, __b13, __b14, __b15,
956 __b16, __b17, __b18, __b19, __b20, __b21, __b22, __b23,
957 __b24, __b25, __b26, __b27, __b28, __b29, __b30, __b31
958 };
959 }
960
961 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_set_epi64x(long long __a,long long __b,long long __c,long long __d)962 _mm256_set_epi64x(long long __a, long long __b, long long __c, long long __d)
963 {
964 return (__m256i)(__v4di){ __d, __c, __b, __a };
965 }
966
967 /* Create vectors with elements in reverse order */
968 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_setr_pd(double __a,double __b,double __c,double __d)969 _mm256_setr_pd(double __a, double __b, double __c, double __d)
970 {
971 return (__m256d){ __a, __b, __c, __d };
972 }
973
974 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_setr_ps(float __a,float __b,float __c,float __d,float __e,float __f,float __g,float __h)975 _mm256_setr_ps(float __a, float __b, float __c, float __d,
976 float __e, float __f, float __g, float __h)
977 {
978 return (__m256){ __a, __b, __c, __d, __e, __f, __g, __h };
979 }
980
981 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_setr_epi32(int __i0,int __i1,int __i2,int __i3,int __i4,int __i5,int __i6,int __i7)982 _mm256_setr_epi32(int __i0, int __i1, int __i2, int __i3,
983 int __i4, int __i5, int __i6, int __i7)
984 {
985 return (__m256i)(__v8si){ __i0, __i1, __i2, __i3, __i4, __i5, __i6, __i7 };
986 }
987
988 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_setr_epi16(short __w15,short __w14,short __w13,short __w12,short __w11,short __w10,short __w09,short __w08,short __w07,short __w06,short __w05,short __w04,short __w03,short __w02,short __w01,short __w00)989 _mm256_setr_epi16(short __w15, short __w14, short __w13, short __w12,
990 short __w11, short __w10, short __w09, short __w08,
991 short __w07, short __w06, short __w05, short __w04,
992 short __w03, short __w02, short __w01, short __w00)
993 {
994 return (__m256i)(__v16hi){ __w15, __w14, __w13, __w12, __w11, __w10, __w09,
995 __w08, __w07, __w06, __w05, __w04, __w03, __w02, __w01, __w00 };
996 }
997
998 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_setr_epi8(char __b31,char __b30,char __b29,char __b28,char __b27,char __b26,char __b25,char __b24,char __b23,char __b22,char __b21,char __b20,char __b19,char __b18,char __b17,char __b16,char __b15,char __b14,char __b13,char __b12,char __b11,char __b10,char __b09,char __b08,char __b07,char __b06,char __b05,char __b04,char __b03,char __b02,char __b01,char __b00)999 _mm256_setr_epi8(char __b31, char __b30, char __b29, char __b28,
1000 char __b27, char __b26, char __b25, char __b24,
1001 char __b23, char __b22, char __b21, char __b20,
1002 char __b19, char __b18, char __b17, char __b16,
1003 char __b15, char __b14, char __b13, char __b12,
1004 char __b11, char __b10, char __b09, char __b08,
1005 char __b07, char __b06, char __b05, char __b04,
1006 char __b03, char __b02, char __b01, char __b00)
1007 {
1008 return (__m256i)(__v32qi){
1009 __b31, __b30, __b29, __b28, __b27, __b26, __b25, __b24,
1010 __b23, __b22, __b21, __b20, __b19, __b18, __b17, __b16,
1011 __b15, __b14, __b13, __b12, __b11, __b10, __b09, __b08,
1012 __b07, __b06, __b05, __b04, __b03, __b02, __b01, __b00 };
1013 }
1014
1015 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_setr_epi64x(long long __a,long long __b,long long __c,long long __d)1016 _mm256_setr_epi64x(long long __a, long long __b, long long __c, long long __d)
1017 {
1018 return (__m256i)(__v4di){ __a, __b, __c, __d };
1019 }
1020
1021 /* Create vectors with repeated elements */
1022 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_set1_pd(double __w)1023 _mm256_set1_pd(double __w)
1024 {
1025 return (__m256d){ __w, __w, __w, __w };
1026 }
1027
1028 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_set1_ps(float __w)1029 _mm256_set1_ps(float __w)
1030 {
1031 return (__m256){ __w, __w, __w, __w, __w, __w, __w, __w };
1032 }
1033
1034 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_set1_epi32(int __i)1035 _mm256_set1_epi32(int __i)
1036 {
1037 return (__m256i)(__v8si){ __i, __i, __i, __i, __i, __i, __i, __i };
1038 }
1039
1040 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_set1_epi16(short __w)1041 _mm256_set1_epi16(short __w)
1042 {
1043 return (__m256i)(__v16hi){ __w, __w, __w, __w, __w, __w, __w, __w, __w, __w,
1044 __w, __w, __w, __w, __w, __w };
1045 }
1046
1047 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_set1_epi8(char __b)1048 _mm256_set1_epi8(char __b)
1049 {
1050 return (__m256i)(__v32qi){ __b, __b, __b, __b, __b, __b, __b, __b, __b, __b,
1051 __b, __b, __b, __b, __b, __b, __b, __b, __b, __b, __b, __b, __b, __b, __b,
1052 __b, __b, __b, __b, __b, __b, __b };
1053 }
1054
1055 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_set1_epi64x(long long __q)1056 _mm256_set1_epi64x(long long __q)
1057 {
1058 return (__m256i)(__v4di){ __q, __q, __q, __q };
1059 }
1060
1061 /* Create __zeroed vectors */
1062 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_setzero_pd(void)1063 _mm256_setzero_pd(void)
1064 {
1065 return (__m256d){ 0, 0, 0, 0 };
1066 }
1067
1068 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_setzero_ps(void)1069 _mm256_setzero_ps(void)
1070 {
1071 return (__m256){ 0, 0, 0, 0, 0, 0, 0, 0 };
1072 }
1073
1074 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_setzero_si256(void)1075 _mm256_setzero_si256(void)
1076 {
1077 return (__m256i){ 0LL, 0LL, 0LL, 0LL };
1078 }
1079
1080 /* Cast between vector types */
1081 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_castpd_ps(__m256d __a)1082 _mm256_castpd_ps(__m256d __a)
1083 {
1084 return (__m256)__a;
1085 }
1086
1087 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_castpd_si256(__m256d __a)1088 _mm256_castpd_si256(__m256d __a)
1089 {
1090 return (__m256i)__a;
1091 }
1092
1093 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_castps_pd(__m256 __a)1094 _mm256_castps_pd(__m256 __a)
1095 {
1096 return (__m256d)__a;
1097 }
1098
1099 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_castps_si256(__m256 __a)1100 _mm256_castps_si256(__m256 __a)
1101 {
1102 return (__m256i)__a;
1103 }
1104
1105 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_castsi256_ps(__m256i __a)1106 _mm256_castsi256_ps(__m256i __a)
1107 {
1108 return (__m256)__a;
1109 }
1110
1111 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_castsi256_pd(__m256i __a)1112 _mm256_castsi256_pd(__m256i __a)
1113 {
1114 return (__m256d)__a;
1115 }
1116
1117 static __inline __m128d __DEFAULT_FN_ATTRS
_mm256_castpd256_pd128(__m256d __a)1118 _mm256_castpd256_pd128(__m256d __a)
1119 {
1120 return __builtin_shufflevector(__a, __a, 0, 1);
1121 }
1122
1123 static __inline __m128 __DEFAULT_FN_ATTRS
_mm256_castps256_ps128(__m256 __a)1124 _mm256_castps256_ps128(__m256 __a)
1125 {
1126 return __builtin_shufflevector(__a, __a, 0, 1, 2, 3);
1127 }
1128
1129 static __inline __m128i __DEFAULT_FN_ATTRS
_mm256_castsi256_si128(__m256i __a)1130 _mm256_castsi256_si128(__m256i __a)
1131 {
1132 return __builtin_shufflevector(__a, __a, 0, 1);
1133 }
1134
1135 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_castpd128_pd256(__m128d __a)1136 _mm256_castpd128_pd256(__m128d __a)
1137 {
1138 return __builtin_shufflevector(__a, __a, 0, 1, -1, -1);
1139 }
1140
1141 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_castps128_ps256(__m128 __a)1142 _mm256_castps128_ps256(__m128 __a)
1143 {
1144 return __builtin_shufflevector(__a, __a, 0, 1, 2, 3, -1, -1, -1, -1);
1145 }
1146
1147 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_castsi128_si256(__m128i __a)1148 _mm256_castsi128_si256(__m128i __a)
1149 {
1150 return __builtin_shufflevector(__a, __a, 0, 1, -1, -1);
1151 }
1152
1153 /*
1154 Vector insert.
1155 We use macros rather than inlines because we only want to accept
1156 invocations where the immediate M is a constant expression.
1157 */
1158 #define _mm256_insertf128_ps(V1, V2, M) __extension__ ({ \
1159 (__m256)__builtin_shufflevector( \
1160 (__v8sf)(__m256)(V1), \
1161 (__v8sf)_mm256_castps128_ps256((__m128)(V2)), \
1162 (((M) & 1) ? 0 : 8), \
1163 (((M) & 1) ? 1 : 9), \
1164 (((M) & 1) ? 2 : 10), \
1165 (((M) & 1) ? 3 : 11), \
1166 (((M) & 1) ? 8 : 4), \
1167 (((M) & 1) ? 9 : 5), \
1168 (((M) & 1) ? 10 : 6), \
1169 (((M) & 1) ? 11 : 7) );})
1170
1171 #define _mm256_insertf128_pd(V1, V2, M) __extension__ ({ \
1172 (__m256d)__builtin_shufflevector( \
1173 (__v4df)(__m256d)(V1), \
1174 (__v4df)_mm256_castpd128_pd256((__m128d)(V2)), \
1175 (((M) & 1) ? 0 : 4), \
1176 (((M) & 1) ? 1 : 5), \
1177 (((M) & 1) ? 4 : 2), \
1178 (((M) & 1) ? 5 : 3) );})
1179
1180 #define _mm256_insertf128_si256(V1, V2, M) __extension__ ({ \
1181 (__m256i)__builtin_shufflevector( \
1182 (__v4di)(__m256i)(V1), \
1183 (__v4di)_mm256_castsi128_si256((__m128i)(V2)), \
1184 (((M) & 1) ? 0 : 4), \
1185 (((M) & 1) ? 1 : 5), \
1186 (((M) & 1) ? 4 : 2), \
1187 (((M) & 1) ? 5 : 3) );})
1188
1189 /*
1190 Vector extract.
1191 We use macros rather than inlines because we only want to accept
1192 invocations where the immediate M is a constant expression.
1193 */
1194 #define _mm256_extractf128_ps(V, M) __extension__ ({ \
1195 (__m128)__builtin_shufflevector( \
1196 (__v8sf)(__m256)(V), \
1197 (__v8sf)(_mm256_setzero_ps()), \
1198 (((M) & 1) ? 4 : 0), \
1199 (((M) & 1) ? 5 : 1), \
1200 (((M) & 1) ? 6 : 2), \
1201 (((M) & 1) ? 7 : 3) );})
1202
1203 #define _mm256_extractf128_pd(V, M) __extension__ ({ \
1204 (__m128d)__builtin_shufflevector( \
1205 (__v4df)(__m256d)(V), \
1206 (__v4df)(_mm256_setzero_pd()), \
1207 (((M) & 1) ? 2 : 0), \
1208 (((M) & 1) ? 3 : 1) );})
1209
1210 #define _mm256_extractf128_si256(V, M) __extension__ ({ \
1211 (__m128i)__builtin_shufflevector( \
1212 (__v4di)(__m256i)(V), \
1213 (__v4di)(_mm256_setzero_si256()), \
1214 (((M) & 1) ? 2 : 0), \
1215 (((M) & 1) ? 3 : 1) );})
1216
1217 /* SIMD load ops (unaligned) */
1218 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_loadu2_m128(float const * __addr_hi,float const * __addr_lo)1219 _mm256_loadu2_m128(float const *__addr_hi, float const *__addr_lo)
1220 {
1221 struct __loadu_ps {
1222 __m128 __v;
1223 } __attribute__((__packed__, __may_alias__));
1224
1225 __m256 __v256 = _mm256_castps128_ps256(((struct __loadu_ps*)__addr_lo)->__v);
1226 return _mm256_insertf128_ps(__v256, ((struct __loadu_ps*)__addr_hi)->__v, 1);
1227 }
1228
1229 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_loadu2_m128d(double const * __addr_hi,double const * __addr_lo)1230 _mm256_loadu2_m128d(double const *__addr_hi, double const *__addr_lo)
1231 {
1232 struct __loadu_pd {
1233 __m128d __v;
1234 } __attribute__((__packed__, __may_alias__));
1235
1236 __m256d __v256 = _mm256_castpd128_pd256(((struct __loadu_pd*)__addr_lo)->__v);
1237 return _mm256_insertf128_pd(__v256, ((struct __loadu_pd*)__addr_hi)->__v, 1);
1238 }
1239
1240 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_loadu2_m128i(__m128i const * __addr_hi,__m128i const * __addr_lo)1241 _mm256_loadu2_m128i(__m128i const *__addr_hi, __m128i const *__addr_lo)
1242 {
1243 struct __loadu_si128 {
1244 __m128i __v;
1245 } __attribute__((__packed__, __may_alias__));
1246 __m256i __v256 = _mm256_castsi128_si256(
1247 ((struct __loadu_si128*)__addr_lo)->__v);
1248 return _mm256_insertf128_si256(__v256,
1249 ((struct __loadu_si128*)__addr_hi)->__v, 1);
1250 }
1251
1252 /* SIMD store ops (unaligned) */
1253 static __inline void __DEFAULT_FN_ATTRS
_mm256_storeu2_m128(float * __addr_hi,float * __addr_lo,__m256 __a)1254 _mm256_storeu2_m128(float *__addr_hi, float *__addr_lo, __m256 __a)
1255 {
1256 __m128 __v128;
1257
1258 __v128 = _mm256_castps256_ps128(__a);
1259 __builtin_ia32_storeups(__addr_lo, __v128);
1260 __v128 = _mm256_extractf128_ps(__a, 1);
1261 __builtin_ia32_storeups(__addr_hi, __v128);
1262 }
1263
1264 static __inline void __DEFAULT_FN_ATTRS
_mm256_storeu2_m128d(double * __addr_hi,double * __addr_lo,__m256d __a)1265 _mm256_storeu2_m128d(double *__addr_hi, double *__addr_lo, __m256d __a)
1266 {
1267 __m128d __v128;
1268
1269 __v128 = _mm256_castpd256_pd128(__a);
1270 __builtin_ia32_storeupd(__addr_lo, __v128);
1271 __v128 = _mm256_extractf128_pd(__a, 1);
1272 __builtin_ia32_storeupd(__addr_hi, __v128);
1273 }
1274
1275 static __inline void __DEFAULT_FN_ATTRS
_mm256_storeu2_m128i(__m128i * __addr_hi,__m128i * __addr_lo,__m256i __a)1276 _mm256_storeu2_m128i(__m128i *__addr_hi, __m128i *__addr_lo, __m256i __a)
1277 {
1278 __m128i __v128;
1279
1280 __v128 = _mm256_castsi256_si128(__a);
1281 __builtin_ia32_storedqu((char *)__addr_lo, (__v16qi)__v128);
1282 __v128 = _mm256_extractf128_si256(__a, 1);
1283 __builtin_ia32_storedqu((char *)__addr_hi, (__v16qi)__v128);
1284 }
1285
1286 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_set_m128(__m128 __hi,__m128 __lo)1287 _mm256_set_m128 (__m128 __hi, __m128 __lo) {
1288 return (__m256) __builtin_shufflevector(__lo, __hi, 0, 1, 2, 3, 4, 5, 6, 7);
1289 }
1290
1291 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_set_m128d(__m128d __hi,__m128d __lo)1292 _mm256_set_m128d (__m128d __hi, __m128d __lo) {
1293 return (__m256d)_mm256_set_m128((__m128)__hi, (__m128)__lo);
1294 }
1295
1296 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_set_m128i(__m128i __hi,__m128i __lo)1297 _mm256_set_m128i (__m128i __hi, __m128i __lo) {
1298 return (__m256i)_mm256_set_m128((__m128)__hi, (__m128)__lo);
1299 }
1300
1301 static __inline __m256 __DEFAULT_FN_ATTRS
_mm256_setr_m128(__m128 __lo,__m128 __hi)1302 _mm256_setr_m128 (__m128 __lo, __m128 __hi) {
1303 return _mm256_set_m128(__hi, __lo);
1304 }
1305
1306 static __inline __m256d __DEFAULT_FN_ATTRS
_mm256_setr_m128d(__m128d __lo,__m128d __hi)1307 _mm256_setr_m128d (__m128d __lo, __m128d __hi) {
1308 return (__m256d)_mm256_set_m128((__m128)__hi, (__m128)__lo);
1309 }
1310
1311 static __inline __m256i __DEFAULT_FN_ATTRS
_mm256_setr_m128i(__m128i __lo,__m128i __hi)1312 _mm256_setr_m128i (__m128i __lo, __m128i __hi) {
1313 return (__m256i)_mm256_set_m128((__m128)__hi, (__m128)__lo);
1314 }
1315
1316 #undef __DEFAULT_FN_ATTRS
1317
1318 #endif /* __AVXINTRIN_H */
1319