avxintrin.h revision 296417
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
38/* We need an explicitly signed variant for char. Note that this shouldn't
39 * appear in the interface though. */
40typedef signed char __v32qs __attribute__((__vector_size__(32)));
41
42typedef float __m256 __attribute__ ((__vector_size__ (32)));
43typedef double __m256d __attribute__((__vector_size__(32)));
44typedef 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 */
50static __inline __m256d __DEFAULT_FN_ATTRS
51_mm256_add_pd(__m256d __a, __m256d __b)
52{
53  return __a+__b;
54}
55
56static __inline __m256 __DEFAULT_FN_ATTRS
57_mm256_add_ps(__m256 __a, __m256 __b)
58{
59  return __a+__b;
60}
61
62static __inline __m256d __DEFAULT_FN_ATTRS
63_mm256_sub_pd(__m256d __a, __m256d __b)
64{
65  return __a-__b;
66}
67
68static __inline __m256 __DEFAULT_FN_ATTRS
69_mm256_sub_ps(__m256 __a, __m256 __b)
70{
71  return __a-__b;
72}
73
74static __inline __m256d __DEFAULT_FN_ATTRS
75_mm256_addsub_pd(__m256d __a, __m256d __b)
76{
77  return (__m256d)__builtin_ia32_addsubpd256((__v4df)__a, (__v4df)__b);
78}
79
80static __inline __m256 __DEFAULT_FN_ATTRS
81_mm256_addsub_ps(__m256 __a, __m256 __b)
82{
83  return (__m256)__builtin_ia32_addsubps256((__v8sf)__a, (__v8sf)__b);
84}
85
86static __inline __m256d __DEFAULT_FN_ATTRS
87_mm256_div_pd(__m256d __a, __m256d __b)
88{
89  return __a / __b;
90}
91
92static __inline __m256 __DEFAULT_FN_ATTRS
93_mm256_div_ps(__m256 __a, __m256 __b)
94{
95  return __a / __b;
96}
97
98static __inline __m256d __DEFAULT_FN_ATTRS
99_mm256_max_pd(__m256d __a, __m256d __b)
100{
101  return (__m256d)__builtin_ia32_maxpd256((__v4df)__a, (__v4df)__b);
102}
103
104static __inline __m256 __DEFAULT_FN_ATTRS
105_mm256_max_ps(__m256 __a, __m256 __b)
106{
107  return (__m256)__builtin_ia32_maxps256((__v8sf)__a, (__v8sf)__b);
108}
109
110static __inline __m256d __DEFAULT_FN_ATTRS
111_mm256_min_pd(__m256d __a, __m256d __b)
112{
113  return (__m256d)__builtin_ia32_minpd256((__v4df)__a, (__v4df)__b);
114}
115
116static __inline __m256 __DEFAULT_FN_ATTRS
117_mm256_min_ps(__m256 __a, __m256 __b)
118{
119  return (__m256)__builtin_ia32_minps256((__v8sf)__a, (__v8sf)__b);
120}
121
122static __inline __m256d __DEFAULT_FN_ATTRS
123_mm256_mul_pd(__m256d __a, __m256d __b)
124{
125  return __a * __b;
126}
127
128static __inline __m256 __DEFAULT_FN_ATTRS
129_mm256_mul_ps(__m256 __a, __m256 __b)
130{
131  return __a * __b;
132}
133
134static __inline __m256d __DEFAULT_FN_ATTRS
135_mm256_sqrt_pd(__m256d __a)
136{
137  return (__m256d)__builtin_ia32_sqrtpd256((__v4df)__a);
138}
139
140static __inline __m256 __DEFAULT_FN_ATTRS
141_mm256_sqrt_ps(__m256 __a)
142{
143  return (__m256)__builtin_ia32_sqrtps256((__v8sf)__a);
144}
145
146static __inline __m256 __DEFAULT_FN_ATTRS
147_mm256_rsqrt_ps(__m256 __a)
148{
149  return (__m256)__builtin_ia32_rsqrtps256((__v8sf)__a);
150}
151
152static __inline __m256 __DEFAULT_FN_ATTRS
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 */
170static __inline __m256d __DEFAULT_FN_ATTRS
171_mm256_and_pd(__m256d __a, __m256d __b)
172{
173  return (__m256d)((__v4di)__a & (__v4di)__b);
174}
175
176static __inline __m256 __DEFAULT_FN_ATTRS
177_mm256_and_ps(__m256 __a, __m256 __b)
178{
179  return (__m256)((__v8si)__a & (__v8si)__b);
180}
181
182static __inline __m256d __DEFAULT_FN_ATTRS
183_mm256_andnot_pd(__m256d __a, __m256d __b)
184{
185  return (__m256d)(~(__v4di)__a & (__v4di)__b);
186}
187
188static __inline __m256 __DEFAULT_FN_ATTRS
189_mm256_andnot_ps(__m256 __a, __m256 __b)
190{
191  return (__m256)(~(__v8si)__a & (__v8si)__b);
192}
193
194static __inline __m256d __DEFAULT_FN_ATTRS
195_mm256_or_pd(__m256d __a, __m256d __b)
196{
197  return (__m256d)((__v4di)__a | (__v4di)__b);
198}
199
200static __inline __m256 __DEFAULT_FN_ATTRS
201_mm256_or_ps(__m256 __a, __m256 __b)
202{
203  return (__m256)((__v8si)__a | (__v8si)__b);
204}
205
206static __inline __m256d __DEFAULT_FN_ATTRS
207_mm256_xor_pd(__m256d __a, __m256d __b)
208{
209  return (__m256d)((__v4di)__a ^ (__v4di)__b);
210}
211
212static __inline __m256 __DEFAULT_FN_ATTRS
213_mm256_xor_ps(__m256 __a, __m256 __b)
214{
215  return (__m256)((__v8si)__a ^ (__v8si)__b);
216}
217
218/* Horizontal arithmetic */
219static __inline __m256d __DEFAULT_FN_ATTRS
220_mm256_hadd_pd(__m256d __a, __m256d __b)
221{
222  return (__m256d)__builtin_ia32_haddpd256((__v4df)__a, (__v4df)__b);
223}
224
225static __inline __m256 __DEFAULT_FN_ATTRS
226_mm256_hadd_ps(__m256 __a, __m256 __b)
227{
228  return (__m256)__builtin_ia32_haddps256((__v8sf)__a, (__v8sf)__b);
229}
230
231static __inline __m256d __DEFAULT_FN_ATTRS
232_mm256_hsub_pd(__m256d __a, __m256d __b)
233{
234  return (__m256d)__builtin_ia32_hsubpd256((__v4df)__a, (__v4df)__b);
235}
236
237static __inline __m256 __DEFAULT_FN_ATTRS
238_mm256_hsub_ps(__m256 __a, __m256 __b)
239{
240  return (__m256)__builtin_ia32_hsubps256((__v8sf)__a, (__v8sf)__b);
241}
242
243/* Vector permutations */
244static __inline __m128d __DEFAULT_FN_ATTRS
245_mm_permutevar_pd(__m128d __a, __m128i __c)
246{
247  return (__m128d)__builtin_ia32_vpermilvarpd((__v2df)__a, (__v2di)__c);
248}
249
250static __inline __m256d __DEFAULT_FN_ATTRS
251_mm256_permutevar_pd(__m256d __a, __m256i __c)
252{
253  return (__m256d)__builtin_ia32_vpermilvarpd256((__v4df)__a, (__v4di)__c);
254}
255
256static __inline __m128 __DEFAULT_FN_ATTRS
257_mm_permutevar_ps(__m128 __a, __m128i __c)
258{
259  return (__m128)__builtin_ia32_vpermilvarps((__v4sf)__a, (__v4si)__c);
260}
261
262static __inline __m256 __DEFAULT_FN_ATTRS
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
329static __inline __m256d __DEFAULT_FN_ATTRS
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
336static __inline __m256 __DEFAULT_FN_ATTRS
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
427static __inline int __DEFAULT_FN_ATTRS
428_mm256_extract_epi32(__m256i __a, const int __imm)
429{
430  __v8si __b = (__v8si)__a;
431  return __b[__imm & 7];
432}
433
434static __inline int __DEFAULT_FN_ATTRS
435_mm256_extract_epi16(__m256i __a, const int __imm)
436{
437  __v16hi __b = (__v16hi)__a;
438  return __b[__imm & 15];
439}
440
441static __inline int __DEFAULT_FN_ATTRS
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__
449static __inline long long  __DEFAULT_FN_ATTRS
450_mm256_extract_epi64(__m256i __a, const int __imm)
451{
452  __v4di __b = (__v4di)__a;
453  return __b[__imm & 3];
454}
455#endif
456
457static __inline __m256i __DEFAULT_FN_ATTRS
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
465static __inline __m256i __DEFAULT_FN_ATTRS
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
473static __inline __m256i __DEFAULT_FN_ATTRS
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__
482static __inline __m256i __DEFAULT_FN_ATTRS
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 */
492static __inline __m256d __DEFAULT_FN_ATTRS
493_mm256_cvtepi32_pd(__m128i __a)
494{
495  return (__m256d)__builtin_ia32_cvtdq2pd256((__v4si) __a);
496}
497
498static __inline __m256 __DEFAULT_FN_ATTRS
499_mm256_cvtepi32_ps(__m256i __a)
500{
501  return (__m256)__builtin_ia32_cvtdq2ps256((__v8si) __a);
502}
503
504static __inline __m128 __DEFAULT_FN_ATTRS
505_mm256_cvtpd_ps(__m256d __a)
506{
507  return (__m128)__builtin_ia32_cvtpd2ps256((__v4df) __a);
508}
509
510static __inline __m256i __DEFAULT_FN_ATTRS
511_mm256_cvtps_epi32(__m256 __a)
512{
513  return (__m256i)__builtin_ia32_cvtps2dq256((__v8sf) __a);
514}
515
516static __inline __m256d __DEFAULT_FN_ATTRS
517_mm256_cvtps_pd(__m128 __a)
518{
519  return (__m256d)__builtin_ia32_cvtps2pd256((__v4sf) __a);
520}
521
522static __inline __m128i __DEFAULT_FN_ATTRS
523_mm256_cvttpd_epi32(__m256d __a)
524{
525  return (__m128i)__builtin_ia32_cvttpd2dq256((__v4df) __a);
526}
527
528static __inline __m128i __DEFAULT_FN_ATTRS
529_mm256_cvtpd_epi32(__m256d __a)
530{
531  return (__m128i)__builtin_ia32_cvtpd2dq256((__v4df) __a);
532}
533
534static __inline __m256i __DEFAULT_FN_ATTRS
535_mm256_cvttps_epi32(__m256 __a)
536{
537  return (__m256i)__builtin_ia32_cvttps2dq256((__v8sf) __a);
538}
539
540/* Vector replicate */
541static __inline __m256 __DEFAULT_FN_ATTRS
542_mm256_movehdup_ps(__m256 __a)
543{
544  return __builtin_shufflevector(__a, __a, 1, 1, 3, 3, 5, 5, 7, 7);
545}
546
547static __inline __m256 __DEFAULT_FN_ATTRS
548_mm256_moveldup_ps(__m256 __a)
549{
550  return __builtin_shufflevector(__a, __a, 0, 0, 2, 2, 4, 4, 6, 6);
551}
552
553static __inline __m256d __DEFAULT_FN_ATTRS
554_mm256_movedup_pd(__m256d __a)
555{
556  return __builtin_shufflevector(__a, __a, 0, 0, 2, 2);
557}
558
559/* Unpack and Interleave */
560static __inline __m256d __DEFAULT_FN_ATTRS
561_mm256_unpackhi_pd(__m256d __a, __m256d __b)
562{
563  return __builtin_shufflevector(__a, __b, 1, 5, 1+2, 5+2);
564}
565
566static __inline __m256d __DEFAULT_FN_ATTRS
567_mm256_unpacklo_pd(__m256d __a, __m256d __b)
568{
569  return __builtin_shufflevector(__a, __b, 0, 4, 0+2, 4+2);
570}
571
572static __inline __m256 __DEFAULT_FN_ATTRS
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
578static __inline __m256 __DEFAULT_FN_ATTRS
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 */
585static __inline int __DEFAULT_FN_ATTRS
586_mm_testz_pd(__m128d __a, __m128d __b)
587{
588  return __builtin_ia32_vtestzpd((__v2df)__a, (__v2df)__b);
589}
590
591static __inline int __DEFAULT_FN_ATTRS
592_mm_testc_pd(__m128d __a, __m128d __b)
593{
594  return __builtin_ia32_vtestcpd((__v2df)__a, (__v2df)__b);
595}
596
597static __inline int __DEFAULT_FN_ATTRS
598_mm_testnzc_pd(__m128d __a, __m128d __b)
599{
600  return __builtin_ia32_vtestnzcpd((__v2df)__a, (__v2df)__b);
601}
602
603static __inline int __DEFAULT_FN_ATTRS
604_mm_testz_ps(__m128 __a, __m128 __b)
605{
606  return __builtin_ia32_vtestzps((__v4sf)__a, (__v4sf)__b);
607}
608
609static __inline int __DEFAULT_FN_ATTRS
610_mm_testc_ps(__m128 __a, __m128 __b)
611{
612  return __builtin_ia32_vtestcps((__v4sf)__a, (__v4sf)__b);
613}
614
615static __inline int __DEFAULT_FN_ATTRS
616_mm_testnzc_ps(__m128 __a, __m128 __b)
617{
618  return __builtin_ia32_vtestnzcps((__v4sf)__a, (__v4sf)__b);
619}
620
621static __inline int __DEFAULT_FN_ATTRS
622_mm256_testz_pd(__m256d __a, __m256d __b)
623{
624  return __builtin_ia32_vtestzpd256((__v4df)__a, (__v4df)__b);
625}
626
627static __inline int __DEFAULT_FN_ATTRS
628_mm256_testc_pd(__m256d __a, __m256d __b)
629{
630  return __builtin_ia32_vtestcpd256((__v4df)__a, (__v4df)__b);
631}
632
633static __inline int __DEFAULT_FN_ATTRS
634_mm256_testnzc_pd(__m256d __a, __m256d __b)
635{
636  return __builtin_ia32_vtestnzcpd256((__v4df)__a, (__v4df)__b);
637}
638
639static __inline int __DEFAULT_FN_ATTRS
640_mm256_testz_ps(__m256 __a, __m256 __b)
641{
642  return __builtin_ia32_vtestzps256((__v8sf)__a, (__v8sf)__b);
643}
644
645static __inline int __DEFAULT_FN_ATTRS
646_mm256_testc_ps(__m256 __a, __m256 __b)
647{
648  return __builtin_ia32_vtestcps256((__v8sf)__a, (__v8sf)__b);
649}
650
651static __inline int __DEFAULT_FN_ATTRS
652_mm256_testnzc_ps(__m256 __a, __m256 __b)
653{
654  return __builtin_ia32_vtestnzcps256((__v8sf)__a, (__v8sf)__b);
655}
656
657static __inline int __DEFAULT_FN_ATTRS
658_mm256_testz_si256(__m256i __a, __m256i __b)
659{
660  return __builtin_ia32_ptestz256((__v4di)__a, (__v4di)__b);
661}
662
663static __inline int __DEFAULT_FN_ATTRS
664_mm256_testc_si256(__m256i __a, __m256i __b)
665{
666  return __builtin_ia32_ptestc256((__v4di)__a, (__v4di)__b);
667}
668
669static __inline int __DEFAULT_FN_ATTRS
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 */
676static __inline int __DEFAULT_FN_ATTRS
677_mm256_movemask_pd(__m256d __a)
678{
679  return __builtin_ia32_movmskpd256((__v4df)__a);
680}
681
682static __inline int __DEFAULT_FN_ATTRS
683_mm256_movemask_ps(__m256 __a)
684{
685  return __builtin_ia32_movmskps256((__v8sf)__a);
686}
687
688/* Vector __zero */
689static __inline void __DEFAULT_FN_ATTRS
690_mm256_zeroall(void)
691{
692  __builtin_ia32_vzeroall();
693}
694
695static __inline void __DEFAULT_FN_ATTRS
696_mm256_zeroupper(void)
697{
698  __builtin_ia32_vzeroupper();
699}
700
701/* Vector load with broadcast */
702static __inline __m128 __DEFAULT_FN_ATTRS
703_mm_broadcast_ss(float const *__a)
704{
705  float __f = *__a;
706  return (__m128)(__v4sf){ __f, __f, __f, __f };
707}
708
709static __inline __m256d __DEFAULT_FN_ATTRS
710_mm256_broadcast_sd(double const *__a)
711{
712  double __d = *__a;
713  return (__m256d)(__v4df){ __d, __d, __d, __d };
714}
715
716static __inline __m256 __DEFAULT_FN_ATTRS
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
723static __inline __m256d __DEFAULT_FN_ATTRS
724_mm256_broadcast_pd(__m128d const *__a)
725{
726  return (__m256d)__builtin_ia32_vbroadcastf128_pd256(__a);
727}
728
729static __inline __m256 __DEFAULT_FN_ATTRS
730_mm256_broadcast_ps(__m128 const *__a)
731{
732  return (__m256)__builtin_ia32_vbroadcastf128_ps256(__a);
733}
734
735/* SIMD load ops */
736static __inline __m256d __DEFAULT_FN_ATTRS
737_mm256_load_pd(double const *__p)
738{
739  return *(__m256d *)__p;
740}
741
742static __inline __m256 __DEFAULT_FN_ATTRS
743_mm256_load_ps(float const *__p)
744{
745  return *(__m256 *)__p;
746}
747
748static __inline __m256d __DEFAULT_FN_ATTRS
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
757static __inline __m256 __DEFAULT_FN_ATTRS
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
766static __inline __m256i __DEFAULT_FN_ATTRS
767_mm256_load_si256(__m256i const *__p)
768{
769  return *__p;
770}
771
772static __inline __m256i __DEFAULT_FN_ATTRS
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
781static __inline __m256i __DEFAULT_FN_ATTRS
782_mm256_lddqu_si256(__m256i const *__p)
783{
784  return (__m256i)__builtin_ia32_lddqu256((char const *)__p);
785}
786
787/* SIMD store ops */
788static __inline void __DEFAULT_FN_ATTRS
789_mm256_store_pd(double *__p, __m256d __a)
790{
791  *(__m256d *)__p = __a;
792}
793
794static __inline void __DEFAULT_FN_ATTRS
795_mm256_store_ps(float *__p, __m256 __a)
796{
797  *(__m256 *)__p = __a;
798}
799
800static __inline void __DEFAULT_FN_ATTRS
801_mm256_storeu_pd(double *__p, __m256d __a)
802{
803  __builtin_ia32_storeupd256(__p, (__v4df)__a);
804}
805
806static __inline void __DEFAULT_FN_ATTRS
807_mm256_storeu_ps(float *__p, __m256 __a)
808{
809  __builtin_ia32_storeups256(__p, (__v8sf)__a);
810}
811
812static __inline void __DEFAULT_FN_ATTRS
813_mm256_store_si256(__m256i *__p, __m256i __a)
814{
815  *__p = __a;
816}
817
818static __inline void __DEFAULT_FN_ATTRS
819_mm256_storeu_si256(__m256i *__p, __m256i __a)
820{
821  __builtin_ia32_storedqu256((char *)__p, (__v32qi)__a);
822}
823
824/* Conditional load ops */
825static __inline __m128d __DEFAULT_FN_ATTRS
826_mm_maskload_pd(double const *__p, __m128i __m)
827{
828  return (__m128d)__builtin_ia32_maskloadpd((const __v2df *)__p, (__v2di)__m);
829}
830
831static __inline __m256d __DEFAULT_FN_ATTRS
832_mm256_maskload_pd(double const *__p, __m256i __m)
833{
834  return (__m256d)__builtin_ia32_maskloadpd256((const __v4df *)__p,
835                                               (__v4di)__m);
836}
837
838static __inline __m128 __DEFAULT_FN_ATTRS
839_mm_maskload_ps(float const *__p, __m128i __m)
840{
841  return (__m128)__builtin_ia32_maskloadps((const __v4sf *)__p, (__v4si)__m);
842}
843
844static __inline __m256 __DEFAULT_FN_ATTRS
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 */
851static __inline void __DEFAULT_FN_ATTRS
852_mm256_maskstore_ps(float *__p, __m256i __m, __m256 __a)
853{
854  __builtin_ia32_maskstoreps256((__v8sf *)__p, (__v8si)__m, (__v8sf)__a);
855}
856
857static __inline void __DEFAULT_FN_ATTRS
858_mm_maskstore_pd(double *__p, __m128i __m, __m128d __a)
859{
860  __builtin_ia32_maskstorepd((__v2df *)__p, (__v2di)__m, (__v2df)__a);
861}
862
863static __inline void __DEFAULT_FN_ATTRS
864_mm256_maskstore_pd(double *__p, __m256i __m, __m256d __a)
865{
866  __builtin_ia32_maskstorepd256((__v4df *)__p, (__v4di)__m, (__v4df)__a);
867}
868
869static __inline void __DEFAULT_FN_ATTRS
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 */
876static __inline void __DEFAULT_FN_ATTRS
877_mm256_stream_si256(__m256i *__a, __m256i __b)
878{
879  __builtin_ia32_movntdq256((__v4di *)__a, (__v4di)__b);
880}
881
882static __inline void __DEFAULT_FN_ATTRS
883_mm256_stream_pd(double *__a, __m256d __b)
884{
885  __builtin_ia32_movntpd256(__a, (__v4df)__b);
886}
887
888static __inline void __DEFAULT_FN_ATTRS
889_mm256_stream_ps(float *__p, __m256 __a)
890{
891  __builtin_ia32_movntps256(__p, (__v8sf)__a);
892}
893
894/* Create vectors */
895static __inline__ __m256d __DEFAULT_FN_ATTRS
896_mm256_undefined_pd()
897{
898  return (__m256d)__builtin_ia32_undef256();
899}
900
901static __inline__ __m256 __DEFAULT_FN_ATTRS
902_mm256_undefined_ps()
903{
904  return (__m256)__builtin_ia32_undef256();
905}
906
907static __inline__ __m256i __DEFAULT_FN_ATTRS
908_mm256_undefined_si256()
909{
910  return (__m256i)__builtin_ia32_undef256();
911}
912
913static __inline __m256d __DEFAULT_FN_ATTRS
914_mm256_set_pd(double __a, double __b, double __c, double __d)
915{
916  return (__m256d){ __d, __c, __b, __a };
917}
918
919static __inline __m256 __DEFAULT_FN_ATTRS
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
926static __inline __m256i __DEFAULT_FN_ATTRS
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
933static __inline __m256i __DEFAULT_FN_ATTRS
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
943static __inline __m256i __DEFAULT_FN_ATTRS
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
961static __inline __m256i __DEFAULT_FN_ATTRS
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 */
968static __inline __m256d __DEFAULT_FN_ATTRS
969_mm256_setr_pd(double __a, double __b, double __c, double __d)
970{
971  return (__m256d){ __a, __b, __c, __d };
972}
973
974static __inline __m256 __DEFAULT_FN_ATTRS
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
981static __inline __m256i __DEFAULT_FN_ATTRS
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
988static __inline __m256i __DEFAULT_FN_ATTRS
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
998static __inline __m256i __DEFAULT_FN_ATTRS
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
1015static __inline __m256i __DEFAULT_FN_ATTRS
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 */
1022static __inline __m256d __DEFAULT_FN_ATTRS
1023_mm256_set1_pd(double __w)
1024{
1025  return (__m256d){ __w, __w, __w, __w };
1026}
1027
1028static __inline __m256 __DEFAULT_FN_ATTRS
1029_mm256_set1_ps(float __w)
1030{
1031  return (__m256){ __w, __w, __w, __w, __w, __w, __w, __w };
1032}
1033
1034static __inline __m256i __DEFAULT_FN_ATTRS
1035_mm256_set1_epi32(int __i)
1036{
1037  return (__m256i)(__v8si){ __i, __i, __i, __i, __i, __i, __i, __i };
1038}
1039
1040static __inline __m256i __DEFAULT_FN_ATTRS
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
1047static __inline __m256i __DEFAULT_FN_ATTRS
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
1055static __inline __m256i __DEFAULT_FN_ATTRS
1056_mm256_set1_epi64x(long long __q)
1057{
1058  return (__m256i)(__v4di){ __q, __q, __q, __q };
1059}
1060
1061/* Create __zeroed vectors */
1062static __inline __m256d __DEFAULT_FN_ATTRS
1063_mm256_setzero_pd(void)
1064{
1065  return (__m256d){ 0, 0, 0, 0 };
1066}
1067
1068static __inline __m256 __DEFAULT_FN_ATTRS
1069_mm256_setzero_ps(void)
1070{
1071  return (__m256){ 0, 0, 0, 0, 0, 0, 0, 0 };
1072}
1073
1074static __inline __m256i __DEFAULT_FN_ATTRS
1075_mm256_setzero_si256(void)
1076{
1077  return (__m256i){ 0LL, 0LL, 0LL, 0LL };
1078}
1079
1080/* Cast between vector types */
1081static __inline __m256 __DEFAULT_FN_ATTRS
1082_mm256_castpd_ps(__m256d __a)
1083{
1084  return (__m256)__a;
1085}
1086
1087static __inline __m256i __DEFAULT_FN_ATTRS
1088_mm256_castpd_si256(__m256d __a)
1089{
1090  return (__m256i)__a;
1091}
1092
1093static __inline __m256d __DEFAULT_FN_ATTRS
1094_mm256_castps_pd(__m256 __a)
1095{
1096  return (__m256d)__a;
1097}
1098
1099static __inline __m256i __DEFAULT_FN_ATTRS
1100_mm256_castps_si256(__m256 __a)
1101{
1102  return (__m256i)__a;
1103}
1104
1105static __inline __m256 __DEFAULT_FN_ATTRS
1106_mm256_castsi256_ps(__m256i __a)
1107{
1108  return (__m256)__a;
1109}
1110
1111static __inline __m256d __DEFAULT_FN_ATTRS
1112_mm256_castsi256_pd(__m256i __a)
1113{
1114  return (__m256d)__a;
1115}
1116
1117static __inline __m128d __DEFAULT_FN_ATTRS
1118_mm256_castpd256_pd128(__m256d __a)
1119{
1120  return __builtin_shufflevector(__a, __a, 0, 1);
1121}
1122
1123static __inline __m128 __DEFAULT_FN_ATTRS
1124_mm256_castps256_ps128(__m256 __a)
1125{
1126  return __builtin_shufflevector(__a, __a, 0, 1, 2, 3);
1127}
1128
1129static __inline __m128i __DEFAULT_FN_ATTRS
1130_mm256_castsi256_si128(__m256i __a)
1131{
1132  return __builtin_shufflevector(__a, __a, 0, 1);
1133}
1134
1135static __inline __m256d __DEFAULT_FN_ATTRS
1136_mm256_castpd128_pd256(__m128d __a)
1137{
1138  return __builtin_shufflevector(__a, __a, 0, 1, -1, -1);
1139}
1140
1141static __inline __m256 __DEFAULT_FN_ATTRS
1142_mm256_castps128_ps256(__m128 __a)
1143{
1144  return __builtin_shufflevector(__a, __a, 0, 1, 2, 3, -1, -1, -1, -1);
1145}
1146
1147static __inline __m256i __DEFAULT_FN_ATTRS
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) */
1218static __inline __m256 __DEFAULT_FN_ATTRS
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
1229static __inline __m256d __DEFAULT_FN_ATTRS
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
1240static __inline __m256i __DEFAULT_FN_ATTRS
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) */
1253static __inline void __DEFAULT_FN_ATTRS
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
1264static __inline void __DEFAULT_FN_ATTRS
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
1275static __inline void __DEFAULT_FN_ATTRS
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
1286static __inline __m256 __DEFAULT_FN_ATTRS
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
1291static __inline __m256d __DEFAULT_FN_ATTRS
1292_mm256_set_m128d (__m128d __hi, __m128d __lo) {
1293  return (__m256d)_mm256_set_m128((__m128)__hi, (__m128)__lo);
1294}
1295
1296static __inline __m256i __DEFAULT_FN_ATTRS
1297_mm256_set_m128i (__m128i __hi, __m128i __lo) {
1298  return (__m256i)_mm256_set_m128((__m128)__hi, (__m128)__lo);
1299}
1300
1301static __inline __m256 __DEFAULT_FN_ATTRS
1302_mm256_setr_m128 (__m128 __lo, __m128 __hi) {
1303  return _mm256_set_m128(__hi, __lo);
1304}
1305
1306static __inline __m256d __DEFAULT_FN_ATTRS
1307_mm256_setr_m128d (__m128d __lo, __m128d __hi) {
1308  return (__m256d)_mm256_set_m128((__m128)__hi, (__m128)__lo);
1309}
1310
1311static __inline __m256i __DEFAULT_FN_ATTRS
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