| /*===---- emmintrin.h - Implementation of SSE2 intrinsics on PowerPC -------=== |
| * |
| * Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions. |
| * See https://llvm.org/LICENSE.txt for license information. |
| * SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception |
| * |
| *===-----------------------------------------------------------------------=== |
| */ |
| |
| /* Implemented from the specification included in the Intel C++ Compiler |
| User Guide and Reference, version 9.0. */ |
| |
| #ifndef NO_WARN_X86_INTRINSICS |
| /* This header file is to help porting code using Intel intrinsics |
| explicitly from x86_64 to powerpc64/powerpc64le. |
| |
| Since X86 SSE2 intrinsics mainly handles __m128i and __m128d type, |
| PowerPC VMX/VSX ISA is a good match for vector float SIMD operations. |
| However scalar float operations in vector (XMM) registers require |
| the POWER8 VSX ISA (2.07) level. There are differences for data |
| format and placement of float scalars in the vector register, which |
| require extra steps to match SSE2 scalar float semantics on POWER. |
| |
| It should be noted that there's much difference between X86_64's |
| MXSCR and PowerISA's FPSCR/VSCR registers. It's recommended to use |
| portable <fenv.h> instead of access MXSCR directly. |
| |
| Most SSE2 scalar float intrinsic operations can be performed more |
| efficiently as C language float scalar operations or optimized to |
| use vector SIMD operations. We recommend this for new applications. |
| */ |
| #error "Please read comment above. Use -DNO_WARN_X86_INTRINSICS to disable this error." |
| #endif |
| |
| #ifndef EMMINTRIN_H_ |
| #define EMMINTRIN_H_ |
| |
| #if defined(__linux__) && defined(__ppc64__) |
| |
| #include <altivec.h> |
| |
| /* We need definitions from the SSE header files. */ |
| #include <xmmintrin.h> |
| |
| /* SSE2 */ |
| typedef __vector double __v2df; |
| typedef __vector long long __v2di; |
| typedef __vector unsigned long long __v2du; |
| typedef __vector int __v4si; |
| typedef __vector unsigned int __v4su; |
| typedef __vector short __v8hi; |
| typedef __vector unsigned short __v8hu; |
| typedef __vector signed char __v16qi; |
| typedef __vector unsigned char __v16qu; |
| |
| /* The Intel API is flexible enough that we must allow aliasing with other |
| vector types, and their scalar components. */ |
| typedef long long __m128i __attribute__ ((__vector_size__ (16), __may_alias__)); |
| typedef double __m128d __attribute__ ((__vector_size__ (16), __may_alias__)); |
| |
| /* Unaligned version of the same types. */ |
| typedef long long __m128i_u __attribute__ ((__vector_size__ (16), __may_alias__, __aligned__ (1))); |
| typedef double __m128d_u __attribute__ ((__vector_size__ (16), __may_alias__, __aligned__ (1))); |
| |
| /* Define two value permute mask. */ |
| #define _MM_SHUFFLE2(x,y) (((x) << 1) | (y)) |
| |
| /* Create a vector with element 0 as F and the rest zero. */ |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_set_sd (double __F) |
| { |
| return __extension__ (__m128d){ __F, 0.0 }; |
| } |
| |
| /* Create a vector with both elements equal to F. */ |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_set1_pd (double __F) |
| { |
| return __extension__ (__m128d){ __F, __F }; |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_set_pd1 (double __F) |
| { |
| return _mm_set1_pd (__F); |
| } |
| |
| /* Create a vector with the lower value X and upper value W. */ |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_set_pd (double __W, double __X) |
| { |
| return __extension__ (__m128d){ __X, __W }; |
| } |
| |
| /* Create a vector with the lower value W and upper value X. */ |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_setr_pd (double __W, double __X) |
| { |
| return __extension__ (__m128d){ __W, __X }; |
| } |
| |
| /* Create an undefined vector. */ |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_undefined_pd (void) |
| { |
| __m128d __Y = __Y; |
| return __Y; |
| } |
| |
| /* Create a vector of zeros. */ |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_setzero_pd (void) |
| { |
| return (__m128d) vec_splats (0); |
| } |
| |
| /* Sets the low DPFP value of A from the low value of B. */ |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_move_sd (__m128d __A, __m128d __B) |
| { |
| __v2df result = (__v2df) __A; |
| result [0] = ((__v2df) __B)[0]; |
| return (__m128d) result; |
| } |
| |
| /* Load two DPFP values from P. The address must be 16-byte aligned. */ |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_load_pd (double const *__P) |
| { |
| return ((__m128d)vec_ld(0, (__v16qu*)__P)); |
| } |
| |
| /* Load two DPFP values from P. The address need not be 16-byte aligned. */ |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_loadu_pd (double const *__P) |
| { |
| return (vec_vsx_ld(0, __P)); |
| } |
| |
| /* Create a vector with all two elements equal to *P. */ |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_load1_pd (double const *__P) |
| { |
| return (vec_splats (*__P)); |
| } |
| |
| /* Create a vector with element 0 as *P and the rest zero. */ |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_load_sd (double const *__P) |
| { |
| return _mm_set_sd (*__P); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_load_pd1 (double const *__P) |
| { |
| return _mm_load1_pd (__P); |
| } |
| |
| /* Load two DPFP values in reverse order. The address must be aligned. */ |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_loadr_pd (double const *__P) |
| { |
| __v2df __tmp = _mm_load_pd (__P); |
| return (__m128d)vec_xxpermdi (__tmp, __tmp, 2); |
| } |
| |
| /* Store two DPFP values. The address must be 16-byte aligned. */ |
| extern __inline void __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_store_pd (double *__P, __m128d __A) |
| { |
| vec_st((__v16qu)__A, 0, (__v16qu*)__P); |
| } |
| |
| /* Store two DPFP values. The address need not be 16-byte aligned. */ |
| extern __inline void __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_storeu_pd (double *__P, __m128d __A) |
| { |
| *(__m128d_u *)__P = __A; |
| } |
| |
| /* Stores the lower DPFP value. */ |
| extern __inline void __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_store_sd (double *__P, __m128d __A) |
| { |
| *__P = ((__v2df)__A)[0]; |
| } |
| |
| extern __inline double __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvtsd_f64 (__m128d __A) |
| { |
| return ((__v2df)__A)[0]; |
| } |
| |
| extern __inline void __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_storel_pd (double *__P, __m128d __A) |
| { |
| _mm_store_sd (__P, __A); |
| } |
| |
| /* Stores the upper DPFP value. */ |
| extern __inline void __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_storeh_pd (double *__P, __m128d __A) |
| { |
| *__P = ((__v2df)__A)[1]; |
| } |
| /* Store the lower DPFP value across two words. |
| The address must be 16-byte aligned. */ |
| extern __inline void __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_store1_pd (double *__P, __m128d __A) |
| { |
| _mm_store_pd (__P, vec_splat (__A, 0)); |
| } |
| |
| extern __inline void __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_store_pd1 (double *__P, __m128d __A) |
| { |
| _mm_store1_pd (__P, __A); |
| } |
| |
| /* Store two DPFP values in reverse order. The address must be aligned. */ |
| extern __inline void __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_storer_pd (double *__P, __m128d __A) |
| { |
| _mm_store_pd (__P, vec_xxpermdi (__A, __A, 2)); |
| } |
| |
| /* Intel intrinsic. */ |
| extern __inline long long __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvtsi128_si64 (__m128i __A) |
| { |
| return ((__v2di)__A)[0]; |
| } |
| |
| /* Microsoft intrinsic. */ |
| extern __inline long long __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvtsi128_si64x (__m128i __A) |
| { |
| return ((__v2di)__A)[0]; |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_add_pd (__m128d __A, __m128d __B) |
| { |
| return (__m128d) ((__v2df)__A + (__v2df)__B); |
| } |
| |
| /* Add the lower double-precision (64-bit) floating-point element in |
| a and b, store the result in the lower element of dst, and copy |
| the upper element from a to the upper element of dst. */ |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_add_sd (__m128d __A, __m128d __B) |
| { |
| __A[0] = __A[0] + __B[0]; |
| return (__A); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_sub_pd (__m128d __A, __m128d __B) |
| { |
| return (__m128d) ((__v2df)__A - (__v2df)__B); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_sub_sd (__m128d __A, __m128d __B) |
| { |
| __A[0] = __A[0] - __B[0]; |
| return (__A); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_mul_pd (__m128d __A, __m128d __B) |
| { |
| return (__m128d) ((__v2df)__A * (__v2df)__B); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_mul_sd (__m128d __A, __m128d __B) |
| { |
| __A[0] = __A[0] * __B[0]; |
| return (__A); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_div_pd (__m128d __A, __m128d __B) |
| { |
| return (__m128d) ((__v2df)__A / (__v2df)__B); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_div_sd (__m128d __A, __m128d __B) |
| { |
| __A[0] = __A[0] / __B[0]; |
| return (__A); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_sqrt_pd (__m128d __A) |
| { |
| return (vec_sqrt (__A)); |
| } |
| |
| /* Return pair {sqrt (B[0]), A[1]}. */ |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_sqrt_sd (__m128d __A, __m128d __B) |
| { |
| __v2df c; |
| c = vec_sqrt ((__v2df) _mm_set1_pd (__B[0])); |
| return (__m128d) _mm_setr_pd (c[0], __A[1]); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_min_pd (__m128d __A, __m128d __B) |
| { |
| return (vec_min (__A, __B)); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_min_sd (__m128d __A, __m128d __B) |
| { |
| __v2df a, b, c; |
| a = vec_splats (__A[0]); |
| b = vec_splats (__B[0]); |
| c = vec_min (a, b); |
| return (__m128d) _mm_setr_pd (c[0], __A[1]); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_max_pd (__m128d __A, __m128d __B) |
| { |
| return (vec_max (__A, __B)); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_max_sd (__m128d __A, __m128d __B) |
| { |
| __v2df a, b, c; |
| a = vec_splats (__A[0]); |
| b = vec_splats (__B[0]); |
| c = vec_max (a, b); |
| return (__m128d) _mm_setr_pd (c[0], __A[1]); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpeq_pd (__m128d __A, __m128d __B) |
| { |
| return ((__m128d)vec_cmpeq ((__v2df) __A, (__v2df) __B)); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmplt_pd (__m128d __A, __m128d __B) |
| { |
| return ((__m128d)vec_cmplt ((__v2df) __A, (__v2df) __B)); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmple_pd (__m128d __A, __m128d __B) |
| { |
| return ((__m128d)vec_cmple ((__v2df) __A, (__v2df) __B)); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpgt_pd (__m128d __A, __m128d __B) |
| { |
| return ((__m128d)vec_cmpgt ((__v2df) __A, (__v2df) __B)); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpge_pd (__m128d __A, __m128d __B) |
| { |
| return ((__m128d)vec_cmpge ((__v2df) __A,(__v2df) __B)); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpneq_pd (__m128d __A, __m128d __B) |
| { |
| __v2df temp = (__v2df) vec_cmpeq ((__v2df) __A, (__v2df)__B); |
| return ((__m128d)vec_nor (temp, temp)); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpnlt_pd (__m128d __A, __m128d __B) |
| { |
| return ((__m128d)vec_cmpge ((__v2df) __A, (__v2df) __B)); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpnle_pd (__m128d __A, __m128d __B) |
| { |
| return ((__m128d)vec_cmpgt ((__v2df) __A, (__v2df) __B)); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpngt_pd (__m128d __A, __m128d __B) |
| { |
| return ((__m128d)vec_cmple ((__v2df) __A, (__v2df) __B)); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpnge_pd (__m128d __A, __m128d __B) |
| { |
| return ((__m128d)vec_cmplt ((__v2df) __A, (__v2df) __B)); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpord_pd (__m128d __A, __m128d __B) |
| { |
| #if _ARCH_PWR8 |
| __v2du c, d; |
| /* Compare against self will return false (0's) if NAN. */ |
| c = (__v2du)vec_cmpeq (__A, __A); |
| d = (__v2du)vec_cmpeq (__B, __B); |
| #else |
| __v2du a, b; |
| __v2du c, d; |
| const __v2du double_exp_mask = {0x7ff0000000000000, 0x7ff0000000000000}; |
| a = (__v2du)vec_abs ((__v2df)__A); |
| b = (__v2du)vec_abs ((__v2df)__B); |
| c = (__v2du)vec_cmpgt (double_exp_mask, a); |
| d = (__v2du)vec_cmpgt (double_exp_mask, b); |
| #endif |
| /* A != NAN and B != NAN. */ |
| return ((__m128d)vec_and(c, d)); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpunord_pd (__m128d __A, __m128d __B) |
| { |
| #if _ARCH_PWR8 |
| __v2du c, d; |
| /* Compare against self will return false (0's) if NAN. */ |
| c = (__v2du)vec_cmpeq ((__v2df)__A, (__v2df)__A); |
| d = (__v2du)vec_cmpeq ((__v2df)__B, (__v2df)__B); |
| /* A == NAN OR B == NAN converts too: |
| NOT(A != NAN) OR NOT(B != NAN). */ |
| c = vec_nor (c, c); |
| return ((__m128d)vec_orc(c, d)); |
| #else |
| __v2du c, d; |
| /* Compare against self will return false (0's) if NAN. */ |
| c = (__v2du)vec_cmpeq ((__v2df)__A, (__v2df)__A); |
| d = (__v2du)vec_cmpeq ((__v2df)__B, (__v2df)__B); |
| /* Convert the true ('1's) is NAN. */ |
| c = vec_nor (c, c); |
| d = vec_nor (d, d); |
| return ((__m128d)vec_or(c, d)); |
| #endif |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpeq_sd(__m128d __A, __m128d __B) |
| { |
| __v2df a, b, c; |
| /* PowerISA VSX does not allow partial (for just lower double) |
| results. So to insure we don't generate spurious exceptions |
| (from the upper double values) we splat the lower double |
| before we do the operation. */ |
| a = vec_splats (__A[0]); |
| b = vec_splats (__B[0]); |
| c = (__v2df) vec_cmpeq(a, b); |
| /* Then we merge the lower double result with the original upper |
| double from __A. */ |
| return (__m128d) _mm_setr_pd (c[0], __A[1]); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmplt_sd (__m128d __A, __m128d __B) |
| { |
| __v2df a, b, c; |
| a = vec_splats (__A[0]); |
| b = vec_splats (__B[0]); |
| c = (__v2df) vec_cmplt(a, b); |
| return (__m128d) _mm_setr_pd (c[0], __A[1]); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmple_sd (__m128d __A, __m128d __B) |
| { |
| __v2df a, b, c; |
| a = vec_splats (__A[0]); |
| b = vec_splats (__B[0]); |
| c = (__v2df) vec_cmple(a, b); |
| return (__m128d) _mm_setr_pd (c[0], __A[1]); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpgt_sd (__m128d __A, __m128d __B) |
| { |
| __v2df a, b, c; |
| a = vec_splats (__A[0]); |
| b = vec_splats (__B[0]); |
| c = (__v2df) vec_cmpgt(a, b); |
| return (__m128d) _mm_setr_pd (c[0], __A[1]); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpge_sd (__m128d __A, __m128d __B) |
| { |
| __v2df a, b, c; |
| a = vec_splats (__A[0]); |
| b = vec_splats (__B[0]); |
| c = (__v2df) vec_cmpge(a, b); |
| return (__m128d) _mm_setr_pd (c[0], __A[1]); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpneq_sd (__m128d __A, __m128d __B) |
| { |
| __v2df a, b, c; |
| a = vec_splats (__A[0]); |
| b = vec_splats (__B[0]); |
| c = (__v2df) vec_cmpeq(a, b); |
| c = vec_nor (c, c); |
| return (__m128d) _mm_setr_pd (c[0], __A[1]); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpnlt_sd (__m128d __A, __m128d __B) |
| { |
| __v2df a, b, c; |
| a = vec_splats (__A[0]); |
| b = vec_splats (__B[0]); |
| /* Not less than is just greater than or equal. */ |
| c = (__v2df) vec_cmpge(a, b); |
| return (__m128d) _mm_setr_pd (c[0], __A[1]); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpnle_sd (__m128d __A, __m128d __B) |
| { |
| __v2df a, b, c; |
| a = vec_splats (__A[0]); |
| b = vec_splats (__B[0]); |
| /* Not less than or equal is just greater than. */ |
| c = (__v2df) vec_cmpge(a, b); |
| return (__m128d) _mm_setr_pd (c[0], __A[1]); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpngt_sd (__m128d __A, __m128d __B) |
| { |
| __v2df a, b, c; |
| a = vec_splats (__A[0]); |
| b = vec_splats (__B[0]); |
| /* Not greater than is just less than or equal. */ |
| c = (__v2df) vec_cmple(a, b); |
| return (__m128d) _mm_setr_pd (c[0], __A[1]); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpnge_sd (__m128d __A, __m128d __B) |
| { |
| __v2df a, b, c; |
| a = vec_splats (__A[0]); |
| b = vec_splats (__B[0]); |
| /* Not greater than or equal is just less than. */ |
| c = (__v2df) vec_cmplt(a, b); |
| return (__m128d) _mm_setr_pd (c[0], __A[1]); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpord_sd (__m128d __A, __m128d __B) |
| { |
| __v2df r; |
| r = (__v2df)_mm_cmpord_pd (vec_splats (__A[0]), vec_splats (__B[0])); |
| return (__m128d) _mm_setr_pd (r[0], ((__v2df)__A)[1]); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpunord_sd (__m128d __A, __m128d __B) |
| { |
| __v2df r; |
| r = _mm_cmpunord_pd (vec_splats (__A[0]), vec_splats (__B[0])); |
| return (__m128d) _mm_setr_pd (r[0], __A[1]); |
| } |
| |
| /* FIXME |
| The __mm_comi??_sd and __mm_ucomi??_sd implementations below are |
| exactly the same because GCC for PowerPC only generates unordered |
| compares (scalar and vector). |
| Technically __mm_comieq_sp et all should be using the ordered |
| compare and signal for QNaNs. The __mm_ucomieq_sd et all should |
| be OK. */ |
| extern __inline int __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_comieq_sd (__m128d __A, __m128d __B) |
| { |
| return (__A[0] == __B[0]); |
| } |
| |
| extern __inline int __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_comilt_sd (__m128d __A, __m128d __B) |
| { |
| return (__A[0] < __B[0]); |
| } |
| |
| extern __inline int __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_comile_sd (__m128d __A, __m128d __B) |
| { |
| return (__A[0] <= __B[0]); |
| } |
| |
| extern __inline int __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_comigt_sd (__m128d __A, __m128d __B) |
| { |
| return (__A[0] > __B[0]); |
| } |
| |
| extern __inline int __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_comige_sd (__m128d __A, __m128d __B) |
| { |
| return (__A[0] >= __B[0]); |
| } |
| |
| extern __inline int __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_comineq_sd (__m128d __A, __m128d __B) |
| { |
| return (__A[0] != __B[0]); |
| } |
| |
| extern __inline int __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_ucomieq_sd (__m128d __A, __m128d __B) |
| { |
| return (__A[0] == __B[0]); |
| } |
| |
| extern __inline int __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_ucomilt_sd (__m128d __A, __m128d __B) |
| { |
| return (__A[0] < __B[0]); |
| } |
| |
| extern __inline int __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_ucomile_sd (__m128d __A, __m128d __B) |
| { |
| return (__A[0] <= __B[0]); |
| } |
| |
| extern __inline int __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_ucomigt_sd (__m128d __A, __m128d __B) |
| { |
| return (__A[0] > __B[0]); |
| } |
| |
| extern __inline int __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_ucomige_sd (__m128d __A, __m128d __B) |
| { |
| return (__A[0] >= __B[0]); |
| } |
| |
| extern __inline int __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_ucomineq_sd (__m128d __A, __m128d __B) |
| { |
| return (__A[0] != __B[0]); |
| } |
| |
| /* Create a vector of Qi, where i is the element number. */ |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_set_epi64x (long long __q1, long long __q0) |
| { |
| return __extension__ (__m128i)(__v2di){ __q0, __q1 }; |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_set_epi64 (__m64 __q1, __m64 __q0) |
| { |
| return _mm_set_epi64x ((long long)__q1, (long long)__q0); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_set_epi32 (int __q3, int __q2, int __q1, int __q0) |
| { |
| return __extension__ (__m128i)(__v4si){ __q0, __q1, __q2, __q3 }; |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_set_epi16 (short __q7, short __q6, short __q5, short __q4, |
| short __q3, short __q2, short __q1, short __q0) |
| { |
| return __extension__ (__m128i)(__v8hi){ |
| __q0, __q1, __q2, __q3, __q4, __q5, __q6, __q7 }; |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_set_epi8 (char __q15, char __q14, char __q13, char __q12, |
| char __q11, char __q10, char __q09, char __q08, |
| char __q07, char __q06, char __q05, char __q04, |
| char __q03, char __q02, char __q01, char __q00) |
| { |
| return __extension__ (__m128i)(__v16qi){ |
| __q00, __q01, __q02, __q03, __q04, __q05, __q06, __q07, |
| __q08, __q09, __q10, __q11, __q12, __q13, __q14, __q15 |
| }; |
| } |
| |
| /* Set all of the elements of the vector to A. */ |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_set1_epi64x (long long __A) |
| { |
| return _mm_set_epi64x (__A, __A); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_set1_epi64 (__m64 __A) |
| { |
| return _mm_set_epi64 (__A, __A); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_set1_epi32 (int __A) |
| { |
| return _mm_set_epi32 (__A, __A, __A, __A); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_set1_epi16 (short __A) |
| { |
| return _mm_set_epi16 (__A, __A, __A, __A, __A, __A, __A, __A); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_set1_epi8 (char __A) |
| { |
| return _mm_set_epi8 (__A, __A, __A, __A, __A, __A, __A, __A, |
| __A, __A, __A, __A, __A, __A, __A, __A); |
| } |
| |
| /* Create a vector of Qi, where i is the element number. |
| The parameter order is reversed from the _mm_set_epi* functions. */ |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_setr_epi64 (__m64 __q0, __m64 __q1) |
| { |
| return _mm_set_epi64 (__q1, __q0); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_setr_epi32 (int __q0, int __q1, int __q2, int __q3) |
| { |
| return _mm_set_epi32 (__q3, __q2, __q1, __q0); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_setr_epi16 (short __q0, short __q1, short __q2, short __q3, |
| short __q4, short __q5, short __q6, short __q7) |
| { |
| return _mm_set_epi16 (__q7, __q6, __q5, __q4, __q3, __q2, __q1, __q0); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_setr_epi8 (char __q00, char __q01, char __q02, char __q03, |
| char __q04, char __q05, char __q06, char __q07, |
| char __q08, char __q09, char __q10, char __q11, |
| char __q12, char __q13, char __q14, char __q15) |
| { |
| return _mm_set_epi8 (__q15, __q14, __q13, __q12, __q11, __q10, __q09, __q08, |
| __q07, __q06, __q05, __q04, __q03, __q02, __q01, __q00); |
| } |
| |
| /* Create a vector with element 0 as *P and the rest zero. */ |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_load_si128 (__m128i const *__P) |
| { |
| return *__P; |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_loadu_si128 (__m128i_u const *__P) |
| { |
| return (__m128i) (vec_vsx_ld(0, (signed int const *)__P)); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_loadl_epi64 (__m128i_u const *__P) |
| { |
| return _mm_set_epi64 ((__m64)0LL, *(__m64 *)__P); |
| } |
| |
| extern __inline void __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_store_si128 (__m128i *__P, __m128i __B) |
| { |
| vec_st ((__v16qu) __B, 0, (__v16qu*)__P); |
| } |
| |
| extern __inline void __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_storeu_si128 (__m128i_u *__P, __m128i __B) |
| { |
| *__P = __B; |
| } |
| |
| extern __inline void __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_storel_epi64 (__m128i_u *__P, __m128i __B) |
| { |
| *(long long *)__P = ((__v2di)__B)[0]; |
| } |
| |
| extern __inline __m64 __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_movepi64_pi64 (__m128i_u __B) |
| { |
| return (__m64) ((__v2di)__B)[0]; |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_movpi64_epi64 (__m64 __A) |
| { |
| return _mm_set_epi64 ((__m64)0LL, __A); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_move_epi64 (__m128i __A) |
| { |
| return _mm_set_epi64 ((__m64)0LL, (__m64)__A[0]); |
| } |
| |
| /* Create an undefined vector. */ |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_undefined_si128 (void) |
| { |
| __m128i __Y = __Y; |
| return __Y; |
| } |
| |
| /* Create a vector of zeros. */ |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_setzero_si128 (void) |
| { |
| return __extension__ (__m128i)(__v4si){ 0, 0, 0, 0 }; |
| } |
| |
| #ifdef _ARCH_PWR8 |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvtepi32_pd (__m128i __A) |
| { |
| __v2di val; |
| /* For LE need to generate Vector Unpack Low Signed Word. |
| Which is generated from unpackh. */ |
| val = (__v2di)vec_unpackh ((__v4si)__A); |
| |
| return (__m128d)vec_ctf (val, 0); |
| } |
| #endif |
| |
| extern __inline __m128 __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvtepi32_ps (__m128i __A) |
| { |
| return ((__m128)vec_ctf((__v4si)__A, 0)); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvtpd_epi32 (__m128d __A) |
| { |
| __v2df rounded = vec_rint (__A); |
| __v4si result, temp; |
| const __v4si vzero = |
| { 0, 0, 0, 0 }; |
| |
| /* VSX Vector truncate Double-Precision to integer and Convert to |
| Signed Integer Word format with Saturate. */ |
| __asm__( |
| "xvcvdpsxws %x0,%x1" |
| : "=wa" (temp) |
| : "wa" (rounded) |
| : ); |
| |
| #ifdef _ARCH_PWR8 |
| temp = vec_mergeo (temp, temp); |
| result = (__v4si) vec_vpkudum ((__vector long long) temp, |
| (__vector long long) vzero); |
| #else |
| { |
| const __v16qu pkperm = {0x00, 0x01, 0x02, 0x03, 0x08, 0x09, 0x0a, 0x0b, |
| 0x14, 0x15, 0x16, 0x17, 0x1c, 0x1d, 0x1e, 0x1f }; |
| result = (__v4si) vec_perm ((__v16qu) temp, (__v16qu) vzero, pkperm); |
| } |
| #endif |
| return (__m128i) result; |
| } |
| |
| extern __inline __m64 __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvtpd_pi32 (__m128d __A) |
| { |
| __m128i result = _mm_cvtpd_epi32(__A); |
| |
| return (__m64) result[0]; |
| } |
| |
| extern __inline __m128 __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvtpd_ps (__m128d __A) |
| { |
| __v4sf result; |
| __v4si temp; |
| const __v4si vzero = { 0, 0, 0, 0 }; |
| |
| __asm__( |
| "xvcvdpsp %x0,%x1" |
| : "=wa" (temp) |
| : "wa" (__A) |
| : ); |
| |
| #ifdef _ARCH_PWR8 |
| temp = vec_mergeo (temp, temp); |
| result = (__v4sf) vec_vpkudum ((__vector long long) temp, |
| (__vector long long) vzero); |
| #else |
| { |
| const __v16qu pkperm = {0x00, 0x01, 0x02, 0x03, 0x08, 0x09, 0x0a, 0x0b, |
| 0x14, 0x15, 0x16, 0x17, 0x1c, 0x1d, 0x1e, 0x1f }; |
| result = (__v4sf) vec_perm ((__v16qu) temp, (__v16qu) vzero, pkperm); |
| } |
| #endif |
| return ((__m128)result); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvttpd_epi32 (__m128d __A) |
| { |
| __v4si result; |
| __v4si temp; |
| const __v4si vzero = { 0, 0, 0, 0 }; |
| |
| /* VSX Vector truncate Double-Precision to integer and Convert to |
| Signed Integer Word format with Saturate. */ |
| __asm__( |
| "xvcvdpsxws %x0,%x1" |
| : "=wa" (temp) |
| : "wa" (__A) |
| : ); |
| |
| #ifdef _ARCH_PWR8 |
| temp = vec_mergeo (temp, temp); |
| result = (__v4si) vec_vpkudum ((__vector long long) temp, |
| (__vector long long) vzero); |
| #else |
| { |
| const __v16qu pkperm = {0x00, 0x01, 0x02, 0x03, 0x08, 0x09, 0x0a, 0x0b, |
| 0x14, 0x15, 0x16, 0x17, 0x1c, 0x1d, 0x1e, 0x1f }; |
| result = (__v4si) vec_perm ((__v16qu) temp, (__v16qu) vzero, pkperm); |
| } |
| #endif |
| |
| return ((__m128i) result); |
| } |
| |
| extern __inline __m64 __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvttpd_pi32 (__m128d __A) |
| { |
| __m128i result = _mm_cvttpd_epi32 (__A); |
| |
| return (__m64) result[0]; |
| } |
| |
| extern __inline int __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvtsi128_si32 (__m128i __A) |
| { |
| return ((__v4si)__A)[0]; |
| } |
| |
| #ifdef _ARCH_PWR8 |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvtpi32_pd (__m64 __A) |
| { |
| __v4si temp; |
| __v2di tmp2; |
| __v2df result; |
| |
| temp = (__v4si)vec_splats (__A); |
| tmp2 = (__v2di)vec_unpackl (temp); |
| result = vec_ctf ((__vector signed long long) tmp2, 0); |
| return (__m128d)result; |
| } |
| #endif |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvtps_epi32 (__m128 __A) |
| { |
| __v4sf rounded; |
| __v4si result; |
| |
| rounded = vec_rint((__v4sf) __A); |
| result = vec_cts (rounded, 0); |
| return (__m128i) result; |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvttps_epi32 (__m128 __A) |
| { |
| __v4si result; |
| |
| result = vec_cts ((__v4sf) __A, 0); |
| return (__m128i) result; |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvtps_pd (__m128 __A) |
| { |
| /* Check if vec_doubleh is defined by <altivec.h>. If so use that. */ |
| #ifdef vec_doubleh |
| return (__m128d) vec_doubleh ((__v4sf)__A); |
| #else |
| /* Otherwise the compiler is not current and so need to generate the |
| equivalent code. */ |
| __v4sf a = (__v4sf)__A; |
| __v4sf temp; |
| __v2df result; |
| #ifdef __LITTLE_ENDIAN__ |
| /* The input float values are in elements {[0], [1]} but the convert |
| instruction needs them in elements {[1], [3]}, So we use two |
| shift left double vector word immediates to get the elements |
| lined up. */ |
| temp = __builtin_vsx_xxsldwi (a, a, 3); |
| temp = __builtin_vsx_xxsldwi (a, temp, 2); |
| #else |
| /* The input float values are in elements {[0], [1]} but the convert |
| instruction needs them in elements {[0], [2]}, So we use two |
| shift left double vector word immediates to get the elements |
| lined up. */ |
| temp = vec_vmrghw (a, a); |
| #endif |
| __asm__( |
| " xvcvspdp %x0,%x1" |
| : "=wa" (result) |
| : "wa" (temp) |
| : ); |
| return (__m128d) result; |
| #endif |
| } |
| |
| extern __inline int __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvtsd_si32 (__m128d __A) |
| { |
| __v2df rounded = vec_rint((__v2df) __A); |
| int result = ((__v2df)rounded)[0]; |
| |
| return result; |
| } |
| /* Intel intrinsic. */ |
| extern __inline long long __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvtsd_si64 (__m128d __A) |
| { |
| __v2df rounded = vec_rint ((__v2df) __A ); |
| long long result = ((__v2df) rounded)[0]; |
| |
| return result; |
| } |
| |
| /* Microsoft intrinsic. */ |
| extern __inline long long __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvtsd_si64x (__m128d __A) |
| { |
| return _mm_cvtsd_si64 ((__v2df)__A); |
| } |
| |
| extern __inline int __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvttsd_si32 (__m128d __A) |
| { |
| int result = ((__v2df)__A)[0]; |
| |
| return result; |
| } |
| |
| /* Intel intrinsic. */ |
| extern __inline long long __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvttsd_si64 (__m128d __A) |
| { |
| long long result = ((__v2df)__A)[0]; |
| |
| return result; |
| } |
| |
| /* Microsoft intrinsic. */ |
| extern __inline long long __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvttsd_si64x (__m128d __A) |
| { |
| return _mm_cvttsd_si64 (__A); |
| } |
| |
| extern __inline __m128 __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvtsd_ss (__m128 __A, __m128d __B) |
| { |
| __v4sf result = (__v4sf)__A; |
| |
| #ifdef __LITTLE_ENDIAN__ |
| __v4sf temp_s; |
| /* Copy double element[0] to element [1] for conversion. */ |
| __v2df temp_b = vec_splat((__v2df)__B, 0); |
| |
| /* Pre-rotate __A left 3 (logically right 1) elements. */ |
| result = __builtin_vsx_xxsldwi (result, result, 3); |
| /* Convert double to single float scalar in a vector. */ |
| __asm__( |
| "xscvdpsp %x0,%x1" |
| : "=wa" (temp_s) |
| : "wa" (temp_b) |
| : ); |
| /* Shift the resulting scalar into vector element [0]. */ |
| result = __builtin_vsx_xxsldwi (result, temp_s, 1); |
| #else |
| result [0] = ((__v2df)__B)[0]; |
| #endif |
| return (__m128) result; |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvtsi32_sd (__m128d __A, int __B) |
| { |
| __v2df result = (__v2df)__A; |
| double db = __B; |
| result [0] = db; |
| return (__m128d)result; |
| } |
| |
| /* Intel intrinsic. */ |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvtsi64_sd (__m128d __A, long long __B) |
| { |
| __v2df result = (__v2df)__A; |
| double db = __B; |
| result [0] = db; |
| return (__m128d)result; |
| } |
| |
| /* Microsoft intrinsic. */ |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvtsi64x_sd (__m128d __A, long long __B) |
| { |
| return _mm_cvtsi64_sd (__A, __B); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvtss_sd (__m128d __A, __m128 __B) |
| { |
| #ifdef __LITTLE_ENDIAN__ |
| /* Use splat to move element [0] into position for the convert. */ |
| __v4sf temp = vec_splat ((__v4sf)__B, 0); |
| __v2df res; |
| /* Convert single float scalar to double in a vector. */ |
| __asm__( |
| "xscvspdp %x0,%x1" |
| : "=wa" (res) |
| : "wa" (temp) |
| : ); |
| return (__m128d) vec_mergel (res, (__v2df)__A); |
| #else |
| __v2df res = (__v2df)__A; |
| res [0] = ((__v4sf)__B) [0]; |
| return (__m128d) res; |
| #endif |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_shuffle_pd(__m128d __A, __m128d __B, const int __mask) |
| { |
| __vector double result; |
| const int litmsk = __mask & 0x3; |
| |
| if (litmsk == 0) |
| result = vec_mergeh (__A, __B); |
| #if __GNUC__ < 6 |
| else if (litmsk == 1) |
| result = vec_xxpermdi (__B, __A, 2); |
| else if (litmsk == 2) |
| result = vec_xxpermdi (__B, __A, 1); |
| #else |
| else if (litmsk == 1) |
| result = vec_xxpermdi (__A, __B, 2); |
| else if (litmsk == 2) |
| result = vec_xxpermdi (__A, __B, 1); |
| #endif |
| else |
| result = vec_mergel (__A, __B); |
| |
| return result; |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_unpackhi_pd (__m128d __A, __m128d __B) |
| { |
| return (__m128d) vec_mergel ((__v2df)__A, (__v2df)__B); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_unpacklo_pd (__m128d __A, __m128d __B) |
| { |
| return (__m128d) vec_mergeh ((__v2df)__A, (__v2df)__B); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_loadh_pd (__m128d __A, double const *__B) |
| { |
| __v2df result = (__v2df)__A; |
| result [1] = *__B; |
| return (__m128d)result; |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_loadl_pd (__m128d __A, double const *__B) |
| { |
| __v2df result = (__v2df)__A; |
| result [0] = *__B; |
| return (__m128d)result; |
| } |
| |
| #ifdef _ARCH_PWR8 |
| /* Intrinsic functions that require PowerISA 2.07 minimum. */ |
| |
| /* Creates a 2-bit mask from the most significant bits of the DPFP values. */ |
| extern __inline int __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_movemask_pd (__m128d __A) |
| { |
| __vector unsigned long long result; |
| static const __vector unsigned int perm_mask = |
| { |
| #ifdef __LITTLE_ENDIAN__ |
| 0x80800040, 0x80808080, 0x80808080, 0x80808080 |
| #else |
| 0x80808080, 0x80808080, 0x80808080, 0x80804000 |
| #endif |
| }; |
| |
| result = ((__vector unsigned long long) |
| vec_vbpermq ((__vector unsigned char) __A, |
| (__vector unsigned char) perm_mask)); |
| |
| #ifdef __LITTLE_ENDIAN__ |
| return result[1]; |
| #else |
| return result[0]; |
| #endif |
| } |
| #endif /* _ARCH_PWR8 */ |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_packs_epi16 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_packs ((__v8hi) __A, (__v8hi)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_packs_epi32 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_packs ((__v4si)__A, (__v4si)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_packus_epi16 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_packsu ((__v8hi) __A, (__v8hi)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_unpackhi_epi8 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_mergel ((__v16qu)__A, (__v16qu)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_unpackhi_epi16 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_mergel ((__v8hu)__A, (__v8hu)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_unpackhi_epi32 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_mergel ((__v4su)__A, (__v4su)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_unpackhi_epi64 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_mergel ((__vector long long) __A, |
| (__vector long long) __B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_unpacklo_epi8 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_mergeh ((__v16qu)__A, (__v16qu)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_unpacklo_epi16 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_mergeh ((__v8hi)__A, (__v8hi)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_unpacklo_epi32 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_mergeh ((__v4si)__A, (__v4si)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_unpacklo_epi64 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_mergeh ((__vector long long) __A, |
| (__vector long long) __B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_add_epi8 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) ((__v16qu)__A + (__v16qu)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_add_epi16 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) ((__v8hu)__A + (__v8hu)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_add_epi32 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) ((__v4su)__A + (__v4su)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_add_epi64 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) ((__v2du)__A + (__v2du)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_adds_epi8 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_adds ((__v16qi)__A, (__v16qi)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_adds_epi16 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_adds ((__v8hi)__A, (__v8hi)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_adds_epu8 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_adds ((__v16qu)__A, (__v16qu)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_adds_epu16 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_adds ((__v8hu)__A, (__v8hu)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_sub_epi8 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) ((__v16qu)__A - (__v16qu)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_sub_epi16 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) ((__v8hu)__A - (__v8hu)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_sub_epi32 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) ((__v4su)__A - (__v4su)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_sub_epi64 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) ((__v2du)__A - (__v2du)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_subs_epi8 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_subs ((__v16qi)__A, (__v16qi)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_subs_epi16 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_subs ((__v8hi)__A, (__v8hi)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_subs_epu8 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_subs ((__v16qu)__A, (__v16qu)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_subs_epu16 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_subs ((__v8hu)__A, (__v8hu)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_madd_epi16 (__m128i __A, __m128i __B) |
| { |
| __vector signed int zero = {0, 0, 0, 0}; |
| |
| return (__m128i) vec_vmsumshm ((__v8hi)__A, (__v8hi)__B, zero); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_mulhi_epi16 (__m128i __A, __m128i __B) |
| { |
| __vector signed int w0, w1; |
| |
| __vector unsigned char xform1 = { |
| #ifdef __LITTLE_ENDIAN__ |
| 0x02, 0x03, 0x12, 0x13, 0x06, 0x07, 0x16, 0x17, |
| 0x0A, 0x0B, 0x1A, 0x1B, 0x0E, 0x0F, 0x1E, 0x1F |
| #else |
| 0x00, 0x01, 0x10, 0x11, 0x04, 0x05, 0x14, 0x15, |
| 0x08, 0x09, 0x18, 0x19, 0x0C, 0x0D, 0x1C, 0x1D |
| #endif |
| }; |
| |
| w0 = vec_vmulesh ((__v8hi)__A, (__v8hi)__B); |
| w1 = vec_vmulosh ((__v8hi)__A, (__v8hi)__B); |
| return (__m128i) vec_perm (w0, w1, xform1); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_mullo_epi16 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) ((__v8hi)__A * (__v8hi)__B); |
| } |
| |
| extern __inline __m64 __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_mul_su32 (__m64 __A, __m64 __B) |
| { |
| unsigned int a = __A; |
| unsigned int b = __B; |
| |
| return ((__m64)a * (__m64)b); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_mul_epu32 (__m128i __A, __m128i __B) |
| { |
| #if __GNUC__ < 8 |
| __v2du result; |
| |
| #ifdef __LITTLE_ENDIAN__ |
| /* VMX Vector Multiply Odd Unsigned Word. */ |
| __asm__( |
| "vmulouw %0,%1,%2" |
| : "=v" (result) |
| : "v" (__A), "v" (__B) |
| : ); |
| #else |
| /* VMX Vector Multiply Even Unsigned Word. */ |
| __asm__( |
| "vmuleuw %0,%1,%2" |
| : "=v" (result) |
| : "v" (__A), "v" (__B) |
| : ); |
| #endif |
| return (__m128i) result; |
| #else |
| return (__m128i) vec_mule ((__v4su)__A, (__v4su)__B); |
| #endif |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_slli_epi16 (__m128i __A, int __B) |
| { |
| __v8hu lshift; |
| __v8hi result = { 0, 0, 0, 0, 0, 0, 0, 0 }; |
| |
| if (__B >= 0 && __B < 16) |
| { |
| if (__builtin_constant_p(__B)) |
| lshift = (__v8hu) vec_splat_s16(__B); |
| else |
| lshift = vec_splats ((unsigned short) __B); |
| |
| result = vec_sl ((__v8hi) __A, lshift); |
| } |
| |
| return (__m128i) result; |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_slli_epi32 (__m128i __A, int __B) |
| { |
| __v4su lshift; |
| __v4si result = { 0, 0, 0, 0 }; |
| |
| if (__B >= 0 && __B < 32) |
| { |
| if (__builtin_constant_p(__B) && __B < 16) |
| lshift = (__v4su) vec_splat_s32(__B); |
| else |
| lshift = vec_splats ((unsigned int) __B); |
| |
| result = vec_sl ((__v4si) __A, lshift); |
| } |
| |
| return (__m128i) result; |
| } |
| |
| #ifdef _ARCH_PWR8 |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_slli_epi64 (__m128i __A, int __B) |
| { |
| __v2du lshift; |
| __v2di result = { 0, 0 }; |
| |
| if (__B >= 0 && __B < 64) |
| { |
| if (__builtin_constant_p(__B) && __B < 16) |
| lshift = (__v2du) vec_splat_s32(__B); |
| else |
| lshift = (__v2du) vec_splats ((unsigned int) __B); |
| |
| result = vec_sl ((__v2di) __A, lshift); |
| } |
| |
| return (__m128i) result; |
| } |
| #endif |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_srai_epi16 (__m128i __A, int __B) |
| { |
| __v8hu rshift = { 15, 15, 15, 15, 15, 15, 15, 15 }; |
| __v8hi result; |
| |
| if (__B < 16) |
| { |
| if (__builtin_constant_p(__B)) |
| rshift = (__v8hu) vec_splat_s16(__B); |
| else |
| rshift = vec_splats ((unsigned short) __B); |
| } |
| result = vec_sra ((__v8hi) __A, rshift); |
| |
| return (__m128i) result; |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_srai_epi32 (__m128i __A, int __B) |
| { |
| __v4su rshift = { 31, 31, 31, 31 }; |
| __v4si result; |
| |
| if (__B < 32) |
| { |
| if (__builtin_constant_p(__B)) |
| { |
| if (__B < 16) |
| rshift = (__v4su) vec_splat_s32(__B); |
| else |
| rshift = (__v4su) vec_splats((unsigned int)__B); |
| } |
| else |
| rshift = vec_splats ((unsigned int) __B); |
| } |
| result = vec_sra ((__v4si) __A, rshift); |
| |
| return (__m128i) result; |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_bslli_si128 (__m128i __A, const int __N) |
| { |
| __v16qu result; |
| const __v16qu zeros = { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 }; |
| |
| if (__N < 16) |
| result = vec_sld ((__v16qu) __A, zeros, __N); |
| else |
| result = zeros; |
| |
| return (__m128i) result; |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_bsrli_si128 (__m128i __A, const int __N) |
| { |
| __v16qu result; |
| const __v16qu zeros = { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 }; |
| |
| if (__N < 16) |
| #ifdef __LITTLE_ENDIAN__ |
| if (__builtin_constant_p(__N)) |
| /* Would like to use Vector Shift Left Double by Octet |
| Immediate here to use the immediate form and avoid |
| load of __N * 8 value into a separate VR. */ |
| result = vec_sld (zeros, (__v16qu) __A, (16 - __N)); |
| else |
| #endif |
| { |
| __v16qu shift = vec_splats((unsigned char)(__N*8)); |
| #ifdef __LITTLE_ENDIAN__ |
| result = vec_sro ((__v16qu)__A, shift); |
| #else |
| result = vec_slo ((__v16qu)__A, shift); |
| #endif |
| } |
| else |
| result = zeros; |
| |
| return (__m128i) result; |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_srli_si128 (__m128i __A, const int __N) |
| { |
| return _mm_bsrli_si128 (__A, __N); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_slli_si128 (__m128i __A, const int _imm5) |
| { |
| __v16qu result; |
| const __v16qu zeros = { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 }; |
| |
| if (_imm5 < 16) |
| #ifdef __LITTLE_ENDIAN__ |
| result = vec_sld ((__v16qu) __A, zeros, _imm5); |
| #else |
| result = vec_sld (zeros, (__v16qu) __A, (16 - _imm5)); |
| #endif |
| else |
| result = zeros; |
| |
| return (__m128i) result; |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| |
| _mm_srli_epi16 (__m128i __A, int __B) |
| { |
| __v8hu rshift; |
| __v8hi result = { 0, 0, 0, 0, 0, 0, 0, 0 }; |
| |
| if (__B < 16) |
| { |
| if (__builtin_constant_p(__B)) |
| rshift = (__v8hu) vec_splat_s16(__B); |
| else |
| rshift = vec_splats ((unsigned short) __B); |
| |
| result = vec_sr ((__v8hi) __A, rshift); |
| } |
| |
| return (__m128i) result; |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_srli_epi32 (__m128i __A, int __B) |
| { |
| __v4su rshift; |
| __v4si result = { 0, 0, 0, 0 }; |
| |
| if (__B < 32) |
| { |
| if (__builtin_constant_p(__B)) |
| { |
| if (__B < 16) |
| rshift = (__v4su) vec_splat_s32(__B); |
| else |
| rshift = (__v4su) vec_splats((unsigned int)__B); |
| } |
| else |
| rshift = vec_splats ((unsigned int) __B); |
| |
| result = vec_sr ((__v4si) __A, rshift); |
| } |
| |
| return (__m128i) result; |
| } |
| |
| #ifdef _ARCH_PWR8 |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_srli_epi64 (__m128i __A, int __B) |
| { |
| __v2du rshift; |
| __v2di result = { 0, 0 }; |
| |
| if (__B < 64) |
| { |
| if (__builtin_constant_p(__B)) |
| { |
| if (__B < 16) |
| rshift = (__v2du) vec_splat_s32(__B); |
| else |
| rshift = (__v2du) vec_splats((unsigned long long)__B); |
| } |
| else |
| rshift = (__v2du) vec_splats ((unsigned int) __B); |
| |
| result = vec_sr ((__v2di) __A, rshift); |
| } |
| |
| return (__m128i) result; |
| } |
| #endif |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_sll_epi16 (__m128i __A, __m128i __B) |
| { |
| __v8hu lshift; |
| __vector __bool short shmask; |
| const __v8hu shmax = { 15, 15, 15, 15, 15, 15, 15, 15 }; |
| __v8hu result; |
| |
| #ifdef __LITTLE_ENDIAN__ |
| lshift = vec_splat ((__v8hu) __B, 0); |
| #else |
| lshift = vec_splat ((__v8hu) __B, 3); |
| #endif |
| shmask = vec_cmple (lshift, shmax); |
| result = vec_sl ((__v8hu) __A, lshift); |
| result = vec_sel ((__v8hu) shmask, result, shmask); |
| |
| return (__m128i) result; |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_sll_epi32 (__m128i __A, __m128i __B) |
| { |
| __v4su lshift; |
| __vector __bool int shmask; |
| const __v4su shmax = { 32, 32, 32, 32 }; |
| __v4su result; |
| #ifdef __LITTLE_ENDIAN__ |
| lshift = vec_splat ((__v4su) __B, 0); |
| #else |
| lshift = vec_splat ((__v4su) __B, 1); |
| #endif |
| shmask = vec_cmplt (lshift, shmax); |
| result = vec_sl ((__v4su) __A, lshift); |
| result = vec_sel ((__v4su) shmask, result, shmask); |
| |
| return (__m128i) result; |
| } |
| |
| #ifdef _ARCH_PWR8 |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_sll_epi64 (__m128i __A, __m128i __B) |
| { |
| __v2du lshift; |
| __vector __bool long long shmask; |
| const __v2du shmax = { 64, 64 }; |
| __v2du result; |
| |
| lshift = vec_splat ((__v2du) __B, 0); |
| shmask = vec_cmplt (lshift, shmax); |
| result = vec_sl ((__v2du) __A, lshift); |
| result = vec_sel ((__v2du) shmask, result, shmask); |
| |
| return (__m128i) result; |
| } |
| #endif |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_sra_epi16 (__m128i __A, __m128i __B) |
| { |
| const __v8hu rshmax = { 15, 15, 15, 15, 15, 15, 15, 15 }; |
| __v8hu rshift; |
| __v8hi result; |
| |
| #ifdef __LITTLE_ENDIAN__ |
| rshift = vec_splat ((__v8hu)__B, 0); |
| #else |
| rshift = vec_splat ((__v8hu)__B, 3); |
| #endif |
| rshift = vec_min (rshift, rshmax); |
| result = vec_sra ((__v8hi) __A, rshift); |
| |
| return (__m128i) result; |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_sra_epi32 (__m128i __A, __m128i __B) |
| { |
| const __v4su rshmax = { 31, 31, 31, 31 }; |
| __v4su rshift; |
| __v4si result; |
| |
| #ifdef __LITTLE_ENDIAN__ |
| rshift = vec_splat ((__v4su)__B, 0); |
| #else |
| rshift = vec_splat ((__v4su)__B, 1); |
| #endif |
| rshift = vec_min (rshift, rshmax); |
| result = vec_sra ((__v4si) __A, rshift); |
| |
| return (__m128i) result; |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_srl_epi16 (__m128i __A, __m128i __B) |
| { |
| __v8hu rshift; |
| __vector __bool short shmask; |
| const __v8hu shmax = { 15, 15, 15, 15, 15, 15, 15, 15 }; |
| __v8hu result; |
| |
| #ifdef __LITTLE_ENDIAN__ |
| rshift = vec_splat ((__v8hu) __B, 0); |
| #else |
| rshift = vec_splat ((__v8hu) __B, 3); |
| #endif |
| shmask = vec_cmple (rshift, shmax); |
| result = vec_sr ((__v8hu) __A, rshift); |
| result = vec_sel ((__v8hu) shmask, result, shmask); |
| |
| return (__m128i) result; |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_srl_epi32 (__m128i __A, __m128i __B) |
| { |
| __v4su rshift; |
| __vector __bool int shmask; |
| const __v4su shmax = { 32, 32, 32, 32 }; |
| __v4su result; |
| |
| #ifdef __LITTLE_ENDIAN__ |
| rshift = vec_splat ((__v4su) __B, 0); |
| #else |
| rshift = vec_splat ((__v4su) __B, 1); |
| #endif |
| shmask = vec_cmplt (rshift, shmax); |
| result = vec_sr ((__v4su) __A, rshift); |
| result = vec_sel ((__v4su) shmask, result, shmask); |
| |
| return (__m128i) result; |
| } |
| |
| #ifdef _ARCH_PWR8 |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_srl_epi64 (__m128i __A, __m128i __B) |
| { |
| __v2du rshift; |
| __vector __bool long long shmask; |
| const __v2du shmax = { 64, 64 }; |
| __v2du result; |
| |
| rshift = vec_splat ((__v2du) __B, 0); |
| shmask = vec_cmplt (rshift, shmax); |
| result = vec_sr ((__v2du) __A, rshift); |
| result = vec_sel ((__v2du) shmask, result, shmask); |
| |
| return (__m128i) result; |
| } |
| #endif |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_and_pd (__m128d __A, __m128d __B) |
| { |
| return (vec_and ((__v2df) __A, (__v2df) __B)); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_andnot_pd (__m128d __A, __m128d __B) |
| { |
| return (vec_andc ((__v2df) __B, (__v2df) __A)); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_or_pd (__m128d __A, __m128d __B) |
| { |
| return (vec_or ((__v2df) __A, (__v2df) __B)); |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_xor_pd (__m128d __A, __m128d __B) |
| { |
| return (vec_xor ((__v2df) __A, (__v2df) __B)); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_and_si128 (__m128i __A, __m128i __B) |
| { |
| return (__m128i)vec_and ((__v2di) __A, (__v2di) __B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_andnot_si128 (__m128i __A, __m128i __B) |
| { |
| return (__m128i)vec_andc ((__v2di) __B, (__v2di) __A); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_or_si128 (__m128i __A, __m128i __B) |
| { |
| return (__m128i)vec_or ((__v2di) __A, (__v2di) __B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_xor_si128 (__m128i __A, __m128i __B) |
| { |
| return (__m128i)vec_xor ((__v2di) __A, (__v2di) __B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpeq_epi8 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_cmpeq ((__v16qi) __A, (__v16qi)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpeq_epi16 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_cmpeq ((__v8hi) __A, (__v8hi)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpeq_epi32 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_cmpeq ((__v4si) __A, (__v4si)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmplt_epi8 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_cmplt ((__v16qi) __A, (__v16qi)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmplt_epi16 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_cmplt ((__v8hi) __A, (__v8hi)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmplt_epi32 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_cmplt ((__v4si) __A, (__v4si)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpgt_epi8 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_cmpgt ((__v16qi) __A, (__v16qi)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpgt_epi16 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_cmpgt ((__v8hi) __A, (__v8hi)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cmpgt_epi32 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_cmpgt ((__v4si) __A, (__v4si)__B); |
| } |
| |
| extern __inline int __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_extract_epi16 (__m128i const __A, int const __N) |
| { |
| return (unsigned short) ((__v8hi)__A)[__N & 7]; |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_insert_epi16 (__m128i const __A, int const __D, int const __N) |
| { |
| __v8hi result = (__v8hi)__A; |
| |
| result [(__N & 7)] = __D; |
| |
| return (__m128i) result; |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_max_epi16 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_max ((__v8hi)__A, (__v8hi)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_max_epu8 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_max ((__v16qu) __A, (__v16qu)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_min_epi16 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_min ((__v8hi) __A, (__v8hi)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_min_epu8 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_min ((__v16qu) __A, (__v16qu)__B); |
| } |
| |
| |
| #ifdef _ARCH_PWR8 |
| /* Intrinsic functions that require PowerISA 2.07 minimum. */ |
| |
| /* Creates a 4-bit mask from the most significant bits of the SPFP values. */ |
| extern __inline int __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_movemask_epi8 (__m128i __A) |
| { |
| __vector unsigned long long result; |
| static const __vector unsigned char perm_mask = |
| { |
| 0x78, 0x70, 0x68, 0x60, 0x58, 0x50, 0x48, 0x40, |
| 0x38, 0x30, 0x28, 0x20, 0x18, 0x10, 0x08, 0x00 |
| }; |
| |
| result = ((__vector unsigned long long) |
| vec_vbpermq ((__vector unsigned char) __A, |
| (__vector unsigned char) perm_mask)); |
| |
| #ifdef __LITTLE_ENDIAN__ |
| return result[1]; |
| #else |
| return result[0]; |
| #endif |
| } |
| #endif /* _ARCH_PWR8 */ |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_mulhi_epu16 (__m128i __A, __m128i __B) |
| { |
| __v4su w0, w1; |
| __v16qu xform1 = { |
| #ifdef __LITTLE_ENDIAN__ |
| 0x02, 0x03, 0x12, 0x13, 0x06, 0x07, 0x16, 0x17, |
| 0x0A, 0x0B, 0x1A, 0x1B, 0x0E, 0x0F, 0x1E, 0x1F |
| #else |
| 0x00, 0x01, 0x10, 0x11, 0x04, 0x05, 0x14, 0x15, |
| 0x08, 0x09, 0x18, 0x19, 0x0C, 0x0D, 0x1C, 0x1D |
| #endif |
| }; |
| |
| w0 = vec_vmuleuh ((__v8hu)__A, (__v8hu)__B); |
| w1 = vec_vmulouh ((__v8hu)__A, (__v8hu)__B); |
| return (__m128i) vec_perm (w0, w1, xform1); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_shufflehi_epi16 (__m128i __A, const int __mask) |
| { |
| unsigned long element_selector_98 = __mask & 0x03; |
| unsigned long element_selector_BA = (__mask >> 2) & 0x03; |
| unsigned long element_selector_DC = (__mask >> 4) & 0x03; |
| unsigned long element_selector_FE = (__mask >> 6) & 0x03; |
| static const unsigned short permute_selectors[4] = |
| { |
| #ifdef __LITTLE_ENDIAN__ |
| 0x0908, 0x0B0A, 0x0D0C, 0x0F0E |
| #else |
| 0x0809, 0x0A0B, 0x0C0D, 0x0E0F |
| #endif |
| }; |
| __v2du pmask = |
| #ifdef __LITTLE_ENDIAN__ |
| { 0x1716151413121110UL, 0UL}; |
| #else |
| { 0x1011121314151617UL, 0UL}; |
| #endif |
| __m64_union t; |
| __v2du a, r; |
| |
| t.as_short[0] = permute_selectors[element_selector_98]; |
| t.as_short[1] = permute_selectors[element_selector_BA]; |
| t.as_short[2] = permute_selectors[element_selector_DC]; |
| t.as_short[3] = permute_selectors[element_selector_FE]; |
| pmask[1] = t.as_m64; |
| a = (__v2du)__A; |
| r = vec_perm (a, a, (__vector unsigned char)pmask); |
| return (__m128i) r; |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_shufflelo_epi16 (__m128i __A, const int __mask) |
| { |
| unsigned long element_selector_10 = __mask & 0x03; |
| unsigned long element_selector_32 = (__mask >> 2) & 0x03; |
| unsigned long element_selector_54 = (__mask >> 4) & 0x03; |
| unsigned long element_selector_76 = (__mask >> 6) & 0x03; |
| static const unsigned short permute_selectors[4] = |
| { |
| #ifdef __LITTLE_ENDIAN__ |
| 0x0100, 0x0302, 0x0504, 0x0706 |
| #else |
| 0x0001, 0x0203, 0x0405, 0x0607 |
| #endif |
| }; |
| __v2du pmask = |
| #ifdef __LITTLE_ENDIAN__ |
| { 0UL, 0x1f1e1d1c1b1a1918UL}; |
| #else |
| { 0UL, 0x18191a1b1c1d1e1fUL}; |
| #endif |
| __m64_union t; |
| __v2du a, r; |
| t.as_short[0] = permute_selectors[element_selector_10]; |
| t.as_short[1] = permute_selectors[element_selector_32]; |
| t.as_short[2] = permute_selectors[element_selector_54]; |
| t.as_short[3] = permute_selectors[element_selector_76]; |
| pmask[0] = t.as_m64; |
| a = (__v2du)__A; |
| r = vec_perm (a, a, (__vector unsigned char)pmask); |
| return (__m128i) r; |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_shuffle_epi32 (__m128i __A, const int __mask) |
| { |
| unsigned long element_selector_10 = __mask & 0x03; |
| unsigned long element_selector_32 = (__mask >> 2) & 0x03; |
| unsigned long element_selector_54 = (__mask >> 4) & 0x03; |
| unsigned long element_selector_76 = (__mask >> 6) & 0x03; |
| static const unsigned int permute_selectors[4] = |
| { |
| #ifdef __LITTLE_ENDIAN__ |
| 0x03020100, 0x07060504, 0x0B0A0908, 0x0F0E0D0C |
| #else |
| 0x00010203, 0x04050607, 0x08090A0B, 0x0C0D0E0F |
| #endif |
| }; |
| __v4su t; |
| |
| t[0] = permute_selectors[element_selector_10]; |
| t[1] = permute_selectors[element_selector_32]; |
| t[2] = permute_selectors[element_selector_54] + 0x10101010; |
| t[3] = permute_selectors[element_selector_76] + 0x10101010; |
| return (__m128i)vec_perm ((__v4si) __A, (__v4si)__A, (__vector unsigned char)t); |
| } |
| |
| extern __inline void __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_maskmoveu_si128 (__m128i __A, __m128i __B, char *__C) |
| { |
| __v2du hibit = { 0x7f7f7f7f7f7f7f7fUL, 0x7f7f7f7f7f7f7f7fUL}; |
| __v16qu mask, tmp; |
| __m128i_u *p = (__m128i_u*)__C; |
| |
| tmp = (__v16qu)_mm_loadu_si128(p); |
| mask = (__v16qu)vec_cmpgt ((__v16qu)__B, (__v16qu)hibit); |
| tmp = vec_sel (tmp, (__v16qu)__A, mask); |
| _mm_storeu_si128 (p, (__m128i)tmp); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_avg_epu8 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_avg ((__v16qu)__A, (__v16qu)__B); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_avg_epu16 (__m128i __A, __m128i __B) |
| { |
| return (__m128i) vec_avg ((__v8hu)__A, (__v8hu)__B); |
| } |
| |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_sad_epu8 (__m128i __A, __m128i __B) |
| { |
| __v16qu a, b; |
| __v16qu vmin, vmax, vabsdiff; |
| __v4si vsum; |
| const __v4su zero = { 0, 0, 0, 0 }; |
| __v4si result; |
| |
| a = (__v16qu) __A; |
| b = (__v16qu) __B; |
| vmin = vec_min (a, b); |
| vmax = vec_max (a, b); |
| vabsdiff = vec_sub (vmax, vmin); |
| /* Sum four groups of bytes into integers. */ |
| vsum = (__vector signed int) vec_sum4s (vabsdiff, zero); |
| /* Sum across four integers with two integer results. */ |
| result = vec_sum2s (vsum, (__vector signed int) zero); |
| /* Rotate the sums into the correct position. */ |
| #ifdef __LITTLE_ENDIAN__ |
| result = vec_sld (result, result, 4); |
| #else |
| result = vec_sld (result, result, 6); |
| #endif |
| /* Rotate the sums into the correct position. */ |
| return (__m128i) result; |
| } |
| |
| extern __inline void __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_stream_si32 (int *__A, int __B) |
| { |
| /* Use the data cache block touch for store transient. */ |
| __asm__ ( |
| "dcbtstt 0,%0" |
| : |
| : "b" (__A) |
| : "memory" |
| ); |
| *__A = __B; |
| } |
| |
| extern __inline void __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_stream_si64 (long long int *__A, long long int __B) |
| { |
| /* Use the data cache block touch for store transient. */ |
| __asm__ ( |
| " dcbtstt 0,%0" |
| : |
| : "b" (__A) |
| : "memory" |
| ); |
| *__A = __B; |
| } |
| |
| extern __inline void __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_stream_si128 (__m128i *__A, __m128i __B) |
| { |
| /* Use the data cache block touch for store transient. */ |
| __asm__ ( |
| "dcbtstt 0,%0" |
| : |
| : "b" (__A) |
| : "memory" |
| ); |
| *__A = __B; |
| } |
| |
| extern __inline void __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_stream_pd (double *__A, __m128d __B) |
| { |
| /* Use the data cache block touch for store transient. */ |
| __asm__ ( |
| "dcbtstt 0,%0" |
| : |
| : "b" (__A) |
| : "memory" |
| ); |
| *(__m128d*)__A = __B; |
| } |
| |
| extern __inline void __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_clflush (void const *__A) |
| { |
| /* Use the data cache block flush. */ |
| __asm__ ( |
| "dcbf 0,%0" |
| : |
| : "b" (__A) |
| : "memory" |
| ); |
| } |
| |
| extern __inline void __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_lfence (void) |
| { |
| /* Use light weight sync for load to load ordering. */ |
| __atomic_thread_fence (__ATOMIC_RELEASE); |
| } |
| |
| extern __inline void __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_mfence (void) |
| { |
| /* Use heavy weight sync for any to any ordering. */ |
| __atomic_thread_fence (__ATOMIC_SEQ_CST); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvtsi32_si128 (int __A) |
| { |
| return _mm_set_epi32 (0, 0, 0, __A); |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvtsi64_si128 (long long __A) |
| { |
| return __extension__ (__m128i)(__v2di){ __A, 0LL }; |
| } |
| |
| /* Microsoft intrinsic. */ |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_cvtsi64x_si128 (long long __A) |
| { |
| return __extension__ (__m128i)(__v2di){ __A, 0LL }; |
| } |
| |
| /* Casts between various SP, DP, INT vector types. Note that these do no |
| conversion of values, they just change the type. */ |
| extern __inline __m128 __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_castpd_ps(__m128d __A) |
| { |
| return (__m128) __A; |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_castpd_si128(__m128d __A) |
| { |
| return (__m128i) __A; |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_castps_pd(__m128 __A) |
| { |
| return (__m128d) __A; |
| } |
| |
| extern __inline __m128i __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_castps_si128(__m128 __A) |
| { |
| return (__m128i) __A; |
| } |
| |
| extern __inline __m128 __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_castsi128_ps(__m128i __A) |
| { |
| return (__m128) __A; |
| } |
| |
| extern __inline __m128d __attribute__((__gnu_inline__, __always_inline__, __artificial__)) |
| _mm_castsi128_pd(__m128i __A) |
| { |
| return (__m128d) __A; |
| } |
| |
| #else |
| #include_next <emmintrin.h> |
| #endif /* defined(__linux__) && defined(__ppc64__) */ |
| |
| #endif /* EMMINTRIN_H_ */ |
