ref: b85e391ac8926467fb305eeddf2ada5608398b59
dir: /vpx_dsp/x86/highbd_intrapred_intrin_ssse3.c/
/*
* Copyright (c) 2017 The WebM project authors. All Rights Reserved.
*
* Use of this source code is governed by a BSD-style license
* that can be found in the LICENSE file in the root of the source
* tree. An additional intellectual property rights grant can be found
* in the file PATENTS. All contributing project authors may
* be found in the AUTHORS file in the root of the source tree.
*/
#include <tmmintrin.h>
#include "./vpx_config.h"
#include "./vpx_dsp_rtcd.h"
#include "vpx/vpx_integer.h"
// -----------------------------------------------------------------------------
/*
; ------------------------------------------
; input: x, y, z, result
;
; trick from pascal
; (x+2y+z+2)>>2 can be calculated as:
; result = avg(x,z)
; result -= xor(x,z) & 1
; result = avg(result,y)
; ------------------------------------------
*/
static INLINE __m128i avg3_epu16(const __m128i *x, const __m128i *y,
const __m128i *z) {
const __m128i one = _mm_set1_epi16(1);
const __m128i a = _mm_avg_epu16(*x, *z);
const __m128i b =
_mm_subs_epu16(a, _mm_and_si128(_mm_xor_si128(*x, *z), one));
return _mm_avg_epu16(b, *y);
}
static INLINE void d207_store_4x8(uint16_t **dst, const ptrdiff_t stride,
const __m128i *a, const __m128i *b) {
_mm_store_si128((__m128i *)*dst, *a);
*dst += stride;
_mm_store_si128((__m128i *)*dst, _mm_alignr_epi8(*b, *a, 4));
*dst += stride;
_mm_store_si128((__m128i *)*dst, _mm_alignr_epi8(*b, *a, 8));
*dst += stride;
_mm_store_si128((__m128i *)*dst, _mm_alignr_epi8(*b, *a, 12));
*dst += stride;
}
void vpx_highbd_d207_predictor_8x8_ssse3(uint16_t *dst, ptrdiff_t stride,
const uint16_t *above,
const uint16_t *left, int bd) {
const __m128i ABCDEFGH = _mm_load_si128((const __m128i *)left);
const __m128i ABCDHHHH = _mm_shufflehi_epi16(ABCDEFGH, 0xff);
const __m128i HHHHHHHH = _mm_unpackhi_epi64(ABCDHHHH, ABCDHHHH);
const __m128i BCDEFGHH = _mm_alignr_epi8(HHHHHHHH, ABCDEFGH, 2);
const __m128i CDEFGHHH = _mm_alignr_epi8(HHHHHHHH, ABCDEFGH, 4);
const __m128i avg3 = avg3_epu16(&ABCDEFGH, &BCDEFGHH, &CDEFGHHH);
const __m128i avg2 = _mm_avg_epu16(ABCDEFGH, BCDEFGHH);
const __m128i out_a = _mm_unpacklo_epi16(avg2, avg3);
const __m128i out_b = _mm_unpackhi_epi16(avg2, avg3);
(void)above;
(void)bd;
d207_store_4x8(&dst, stride, &out_a, &out_b);
d207_store_4x8(&dst, stride, &out_b, &HHHHHHHH);
}
static INLINE void d207_store_4x16(uint16_t **dst, const ptrdiff_t stride,
const __m128i *a, const __m128i *b,
const __m128i *c) {
_mm_store_si128((__m128i *)*dst, *a);
_mm_store_si128((__m128i *)(*dst + 8), *b);
*dst += stride;
_mm_store_si128((__m128i *)*dst, _mm_alignr_epi8(*b, *a, 4));
_mm_store_si128((__m128i *)(*dst + 8), _mm_alignr_epi8(*c, *b, 4));
*dst += stride;
_mm_store_si128((__m128i *)*dst, _mm_alignr_epi8(*b, *a, 8));
_mm_store_si128((__m128i *)(*dst + 8), _mm_alignr_epi8(*c, *b, 8));
*dst += stride;
_mm_store_si128((__m128i *)*dst, _mm_alignr_epi8(*b, *a, 12));
_mm_store_si128((__m128i *)(*dst + 8), _mm_alignr_epi8(*c, *b, 12));
*dst += stride;
}
void vpx_highbd_d207_predictor_16x16_ssse3(uint16_t *dst, ptrdiff_t stride,
const uint16_t *above,
const uint16_t *left, int bd) {
const __m128i A0 = _mm_load_si128((const __m128i *)left);
const __m128i A1 = _mm_load_si128((const __m128i *)(left + 8));
const __m128i LR0 = _mm_shufflehi_epi16(A1, 0xff);
const __m128i LR = _mm_unpackhi_epi64(LR0, LR0);
const __m128i B0 = _mm_alignr_epi8(A1, A0, 2);
const __m128i B1 = _mm_alignr_epi8(LR, A1, 2);
const __m128i C0 = _mm_alignr_epi8(A1, A0, 4);
const __m128i C1 = _mm_alignr_epi8(LR, A1, 4);
const __m128i avg3_0 = avg3_epu16(&A0, &B0, &C0);
const __m128i avg3_1 = avg3_epu16(&A1, &B1, &C1);
const __m128i avg2_0 = _mm_avg_epu16(A0, B0);
const __m128i avg2_1 = _mm_avg_epu16(A1, B1);
const __m128i out_a = _mm_unpacklo_epi16(avg2_0, avg3_0);
const __m128i out_b = _mm_unpackhi_epi16(avg2_0, avg3_0);
const __m128i out_c = _mm_unpacklo_epi16(avg2_1, avg3_1);
const __m128i out_d = _mm_unpackhi_epi16(avg2_1, avg3_1);
(void)above;
(void)bd;
d207_store_4x16(&dst, stride, &out_a, &out_b, &out_c);
d207_store_4x16(&dst, stride, &out_b, &out_c, &out_d);
d207_store_4x16(&dst, stride, &out_c, &out_d, &LR);
d207_store_4x16(&dst, stride, &out_d, &LR, &LR);
}
static INLINE void d207_store_4x32(uint16_t **dst, const ptrdiff_t stride,
const __m128i *a, const __m128i *b,
const __m128i *c, const __m128i *d,
const __m128i *e) {
_mm_store_si128((__m128i *)*dst, *a);
_mm_store_si128((__m128i *)(*dst + 8), *b);
_mm_store_si128((__m128i *)(*dst + 16), *c);
_mm_store_si128((__m128i *)(*dst + 24), *d);
*dst += stride;
_mm_store_si128((__m128i *)*dst, _mm_alignr_epi8(*b, *a, 4));
_mm_store_si128((__m128i *)(*dst + 8), _mm_alignr_epi8(*c, *b, 4));
_mm_store_si128((__m128i *)(*dst + 16), _mm_alignr_epi8(*d, *c, 4));
_mm_store_si128((__m128i *)(*dst + 24), _mm_alignr_epi8(*e, *d, 4));
*dst += stride;
_mm_store_si128((__m128i *)*dst, _mm_alignr_epi8(*b, *a, 8));
_mm_store_si128((__m128i *)(*dst + 8), _mm_alignr_epi8(*c, *b, 8));
_mm_store_si128((__m128i *)(*dst + 16), _mm_alignr_epi8(*d, *c, 8));
_mm_store_si128((__m128i *)(*dst + 24), _mm_alignr_epi8(*e, *d, 8));
*dst += stride;
_mm_store_si128((__m128i *)*dst, _mm_alignr_epi8(*b, *a, 12));
_mm_store_si128((__m128i *)(*dst + 8), _mm_alignr_epi8(*c, *b, 12));
_mm_store_si128((__m128i *)(*dst + 16), _mm_alignr_epi8(*d, *c, 12));
_mm_store_si128((__m128i *)(*dst + 24), _mm_alignr_epi8(*e, *d, 12));
*dst += stride;
}
void vpx_highbd_d207_predictor_32x32_ssse3(uint16_t *dst, ptrdiff_t stride,
const uint16_t *above,
const uint16_t *left, int bd) {
const __m128i A0 = _mm_load_si128((const __m128i *)left);
const __m128i A1 = _mm_load_si128((const __m128i *)(left + 8));
const __m128i A2 = _mm_load_si128((const __m128i *)(left + 16));
const __m128i A3 = _mm_load_si128((const __m128i *)(left + 24));
const __m128i LR0 = _mm_shufflehi_epi16(A3, 0xff);
const __m128i LR = _mm_unpackhi_epi64(LR0, LR0);
const __m128i B0 = _mm_alignr_epi8(A1, A0, 2);
const __m128i B1 = _mm_alignr_epi8(A2, A1, 2);
const __m128i B2 = _mm_alignr_epi8(A3, A2, 2);
const __m128i B3 = _mm_alignr_epi8(LR, A3, 2);
const __m128i C0 = _mm_alignr_epi8(A1, A0, 4);
const __m128i C1 = _mm_alignr_epi8(A2, A1, 4);
const __m128i C2 = _mm_alignr_epi8(A3, A2, 4);
const __m128i C3 = _mm_alignr_epi8(LR, A3, 4);
const __m128i avg3_0 = avg3_epu16(&A0, &B0, &C0);
const __m128i avg3_1 = avg3_epu16(&A1, &B1, &C1);
const __m128i avg3_2 = avg3_epu16(&A2, &B2, &C2);
const __m128i avg3_3 = avg3_epu16(&A3, &B3, &C3);
const __m128i avg2_0 = _mm_avg_epu16(A0, B0);
const __m128i avg2_1 = _mm_avg_epu16(A1, B1);
const __m128i avg2_2 = _mm_avg_epu16(A2, B2);
const __m128i avg2_3 = _mm_avg_epu16(A3, B3);
const __m128i out_a = _mm_unpacklo_epi16(avg2_0, avg3_0);
const __m128i out_b = _mm_unpackhi_epi16(avg2_0, avg3_0);
const __m128i out_c = _mm_unpacklo_epi16(avg2_1, avg3_1);
const __m128i out_d = _mm_unpackhi_epi16(avg2_1, avg3_1);
const __m128i out_e = _mm_unpacklo_epi16(avg2_2, avg3_2);
const __m128i out_f = _mm_unpackhi_epi16(avg2_2, avg3_2);
const __m128i out_g = _mm_unpacklo_epi16(avg2_3, avg3_3);
const __m128i out_h = _mm_unpackhi_epi16(avg2_3, avg3_3);
(void)above;
(void)bd;
d207_store_4x32(&dst, stride, &out_a, &out_b, &out_c, &out_d, &out_e);
d207_store_4x32(&dst, stride, &out_b, &out_c, &out_d, &out_e, &out_f);
d207_store_4x32(&dst, stride, &out_c, &out_d, &out_e, &out_f, &out_g);
d207_store_4x32(&dst, stride, &out_d, &out_e, &out_f, &out_g, &out_h);
d207_store_4x32(&dst, stride, &out_e, &out_f, &out_g, &out_h, &LR);
d207_store_4x32(&dst, stride, &out_f, &out_g, &out_h, &LR, &LR);
d207_store_4x32(&dst, stride, &out_g, &out_h, &LR, &LR, &LR);
d207_store_4x32(&dst, stride, &out_h, &LR, &LR, &LR, &LR);
}
static INLINE void d63_store_4x8(uint16_t **dst, const ptrdiff_t stride,
__m128i *a, __m128i *b, const __m128i *ar) {
_mm_store_si128((__m128i *)*dst, *a);
*dst += stride;
_mm_store_si128((__m128i *)*dst, *b);
*dst += stride;
*a = _mm_alignr_epi8(*ar, *a, 2);
*b = _mm_alignr_epi8(*ar, *b, 2);
_mm_store_si128((__m128i *)*dst, *a);
*dst += stride;
_mm_store_si128((__m128i *)*dst, *b);
*dst += stride;
*a = _mm_alignr_epi8(*ar, *a, 2);
*b = _mm_alignr_epi8(*ar, *b, 2);
}
void vpx_highbd_d63_predictor_8x8_ssse3(uint16_t *dst, ptrdiff_t stride,
const uint16_t *above,
const uint16_t *left, int bd) {
const __m128i ABCDEFGH = _mm_load_si128((const __m128i *)above);
const __m128i ABCDHHHH = _mm_shufflehi_epi16(ABCDEFGH, 0xff);
const __m128i HHHHHHHH = _mm_unpackhi_epi64(ABCDHHHH, ABCDHHHH);
const __m128i BCDEFGHH = _mm_alignr_epi8(HHHHHHHH, ABCDEFGH, 2);
const __m128i CDEFGHHH = _mm_alignr_epi8(HHHHHHHH, ABCDEFGH, 4);
__m128i avg3 = avg3_epu16(&ABCDEFGH, &BCDEFGHH, &CDEFGHHH);
__m128i avg2 = _mm_avg_epu16(ABCDEFGH, BCDEFGHH);
(void)left;
(void)bd;
d63_store_4x8(&dst, stride, &avg2, &avg3, &HHHHHHHH);
d63_store_4x8(&dst, stride, &avg2, &avg3, &HHHHHHHH);
}
void vpx_highbd_d63_predictor_16x16_ssse3(uint16_t *dst, ptrdiff_t stride,
const uint16_t *above,
const uint16_t *left, int bd) {
const __m128i A0 = _mm_load_si128((const __m128i *)above);
const __m128i A1 = _mm_load_si128((const __m128i *)(above + 8));
const __m128i AR0 = _mm_shufflehi_epi16(A1, 0xff);
const __m128i AR = _mm_unpackhi_epi64(AR0, AR0);
const __m128i B0 = _mm_alignr_epi8(A1, A0, 2);
const __m128i B1 = _mm_alignr_epi8(AR, A1, 2);
const __m128i C0 = _mm_alignr_epi8(A1, A0, 4);
const __m128i C1 = _mm_alignr_epi8(AR, A1, 4);
__m128i avg3_0 = avg3_epu16(&A0, &B0, &C0);
__m128i avg3_1 = avg3_epu16(&A1, &B1, &C1);
__m128i avg2_0 = _mm_avg_epu16(A0, B0);
__m128i avg2_1 = _mm_avg_epu16(A1, B1);
int i;
(void)left;
(void)bd;
for (i = 0; i < 14; i += 2) {
_mm_store_si128((__m128i *)dst, avg2_0);
_mm_store_si128((__m128i *)(dst + 8), avg2_1);
dst += stride;
_mm_store_si128((__m128i *)dst, avg3_0);
_mm_store_si128((__m128i *)(dst + 8), avg3_1);
dst += stride;
avg2_0 = _mm_alignr_epi8(avg2_1, avg2_0, 2);
avg2_1 = _mm_alignr_epi8(AR, avg2_1, 2);
avg3_0 = _mm_alignr_epi8(avg3_1, avg3_0, 2);
avg3_1 = _mm_alignr_epi8(AR, avg3_1, 2);
}
_mm_store_si128((__m128i *)dst, avg2_0);
_mm_store_si128((__m128i *)(dst + 8), avg2_1);
dst += stride;
_mm_store_si128((__m128i *)dst, avg3_0);
_mm_store_si128((__m128i *)(dst + 8), avg3_1);
}
void vpx_highbd_d63_predictor_32x32_ssse3(uint16_t *dst, ptrdiff_t stride,
const uint16_t *above,
const uint16_t *left, int bd) {
const __m128i A0 = _mm_load_si128((const __m128i *)above);
const __m128i A1 = _mm_load_si128((const __m128i *)(above + 8));
const __m128i A2 = _mm_load_si128((const __m128i *)(above + 16));
const __m128i A3 = _mm_load_si128((const __m128i *)(above + 24));
const __m128i AR0 = _mm_shufflehi_epi16(A3, 0xff);
const __m128i AR = _mm_unpackhi_epi64(AR0, AR0);
const __m128i B0 = _mm_alignr_epi8(A1, A0, 2);
const __m128i B1 = _mm_alignr_epi8(A2, A1, 2);
const __m128i B2 = _mm_alignr_epi8(A3, A2, 2);
const __m128i B3 = _mm_alignr_epi8(AR, A3, 2);
const __m128i C0 = _mm_alignr_epi8(A1, A0, 4);
const __m128i C1 = _mm_alignr_epi8(A2, A1, 4);
const __m128i C2 = _mm_alignr_epi8(A3, A2, 4);
const __m128i C3 = _mm_alignr_epi8(AR, A3, 4);
__m128i avg3_0 = avg3_epu16(&A0, &B0, &C0);
__m128i avg3_1 = avg3_epu16(&A1, &B1, &C1);
__m128i avg3_2 = avg3_epu16(&A2, &B2, &C2);
__m128i avg3_3 = avg3_epu16(&A3, &B3, &C3);
__m128i avg2_0 = _mm_avg_epu16(A0, B0);
__m128i avg2_1 = _mm_avg_epu16(A1, B1);
__m128i avg2_2 = _mm_avg_epu16(A2, B2);
__m128i avg2_3 = _mm_avg_epu16(A3, B3);
int i;
(void)left;
(void)bd;
for (i = 0; i < 30; i += 2) {
_mm_store_si128((__m128i *)dst, avg2_0);
_mm_store_si128((__m128i *)(dst + 8), avg2_1);
_mm_store_si128((__m128i *)(dst + 16), avg2_2);
_mm_store_si128((__m128i *)(dst + 24), avg2_3);
dst += stride;
_mm_store_si128((__m128i *)dst, avg3_0);
_mm_store_si128((__m128i *)(dst + 8), avg3_1);
_mm_store_si128((__m128i *)(dst + 16), avg3_2);
_mm_store_si128((__m128i *)(dst + 24), avg3_3);
dst += stride;
avg2_0 = _mm_alignr_epi8(avg2_1, avg2_0, 2);
avg2_1 = _mm_alignr_epi8(avg2_2, avg2_1, 2);
avg2_2 = _mm_alignr_epi8(avg2_3, avg2_2, 2);
avg2_3 = _mm_alignr_epi8(AR, avg2_3, 2);
avg3_0 = _mm_alignr_epi8(avg3_1, avg3_0, 2);
avg3_1 = _mm_alignr_epi8(avg3_2, avg3_1, 2);
avg3_2 = _mm_alignr_epi8(avg3_3, avg3_2, 2);
avg3_3 = _mm_alignr_epi8(AR, avg3_3, 2);
}
_mm_store_si128((__m128i *)dst, avg2_0);
_mm_store_si128((__m128i *)(dst + 8), avg2_1);
_mm_store_si128((__m128i *)(dst + 16), avg2_2);
_mm_store_si128((__m128i *)(dst + 24), avg2_3);
dst += stride;
_mm_store_si128((__m128i *)dst, avg3_0);
_mm_store_si128((__m128i *)(dst + 8), avg3_1);
_mm_store_si128((__m128i *)(dst + 16), avg3_2);
_mm_store_si128((__m128i *)(dst + 24), avg3_3);
}