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