blob: 67fb4d32bd7a91767edc38c1d5ff4dad8ede632c [file] [log] [blame]
/*
* Copyright (c) 2018, Alliance for Open Media. All rights reserved
*
* This source code is subject to the terms of the BSD 2 Clause License and
* the Alliance for Open Media Patent License 1.0. If the BSD 2 Clause License
* was not distributed with this source code in the LICENSE file, you can
* obtain it at www.aomedia.org/license/software. If the Alliance for Open
* Media Patent License 1.0 was not distributed with this source code in the
* PATENTS file, you can obtain it at www.aomedia.org/license/patent.
*/
#include <smmintrin.h> // SSE4.1
#include <immintrin.h> // AVX2
#include <assert.h>
#include "aom/aom_integer.h"
#include "aom_ports/mem.h"
#include "aom_dsp/aom_dsp_common.h"
#include "aom_dsp/x86/synonyms.h"
#include "aom_dsp/x86/synonyms_avx2.h"
#include "aom_dsp/x86/blend_sse4.h"
#include "aom_dsp/x86/blend_mask_sse4.h"
#include "config/aom_dsp_rtcd.h"
static INLINE void blend_a64_d16_mask_w16_avx2(
uint8_t *dst, const CONV_BUF_TYPE *src0, const CONV_BUF_TYPE *src1,
const __m256i *m0, const __m256i *v_round_offset, const __m256i *v_maxval,
int shift) {
const __m256i max_minus_m0 = _mm256_sub_epi16(*v_maxval, *m0);
const __m256i s0_0 = yy_loadu_256(src0);
const __m256i s1_0 = yy_loadu_256(src1);
__m256i res0_lo = _mm256_madd_epi16(_mm256_unpacklo_epi16(s0_0, s1_0),
_mm256_unpacklo_epi16(*m0, max_minus_m0));
__m256i res0_hi = _mm256_madd_epi16(_mm256_unpackhi_epi16(s0_0, s1_0),
_mm256_unpackhi_epi16(*m0, max_minus_m0));
res0_lo =
_mm256_srai_epi32(_mm256_sub_epi32(res0_lo, *v_round_offset), shift);
res0_hi =
_mm256_srai_epi32(_mm256_sub_epi32(res0_hi, *v_round_offset), shift);
const __m256i res0 = _mm256_packs_epi32(res0_lo, res0_hi);
__m256i res = _mm256_packus_epi16(res0, res0);
res = _mm256_permute4x64_epi64(res, 0xd8);
_mm_storeu_si128((__m128i *)(dst), _mm256_castsi256_si128(res));
}
static INLINE void blend_a64_d16_mask_w32_avx2(
uint8_t *dst, const CONV_BUF_TYPE *src0, const CONV_BUF_TYPE *src1,
const __m256i *m0, const __m256i *m1, const __m256i *v_round_offset,
const __m256i *v_maxval, int shift) {
const __m256i max_minus_m0 = _mm256_sub_epi16(*v_maxval, *m0);
const __m256i max_minus_m1 = _mm256_sub_epi16(*v_maxval, *m1);
const __m256i s0_0 = yy_loadu_256(src0);
const __m256i s0_1 = yy_loadu_256(src0 + 16);
const __m256i s1_0 = yy_loadu_256(src1);
const __m256i s1_1 = yy_loadu_256(src1 + 16);
__m256i res0_lo = _mm256_madd_epi16(_mm256_unpacklo_epi16(s0_0, s1_0),
_mm256_unpacklo_epi16(*m0, max_minus_m0));
__m256i res0_hi = _mm256_madd_epi16(_mm256_unpackhi_epi16(s0_0, s1_0),
_mm256_unpackhi_epi16(*m0, max_minus_m0));
__m256i res1_lo = _mm256_madd_epi16(_mm256_unpacklo_epi16(s0_1, s1_1),
_mm256_unpacklo_epi16(*m1, max_minus_m1));
__m256i res1_hi = _mm256_madd_epi16(_mm256_unpackhi_epi16(s0_1, s1_1),
_mm256_unpackhi_epi16(*m1, max_minus_m1));
res0_lo =
_mm256_srai_epi32(_mm256_sub_epi32(res0_lo, *v_round_offset), shift);
res0_hi =
_mm256_srai_epi32(_mm256_sub_epi32(res0_hi, *v_round_offset), shift);
res1_lo =
_mm256_srai_epi32(_mm256_sub_epi32(res1_lo, *v_round_offset), shift);
res1_hi =
_mm256_srai_epi32(_mm256_sub_epi32(res1_hi, *v_round_offset), shift);
const __m256i res0 = _mm256_packs_epi32(res0_lo, res0_hi);
const __m256i res1 = _mm256_packs_epi32(res1_lo, res1_hi);
__m256i res = _mm256_packus_epi16(res0, res1);
res = _mm256_permute4x64_epi64(res, 0xd8);
_mm256_storeu_si256((__m256i *)(dst), res);
}
static INLINE void lowbd_blend_a64_d16_mask_subw0_subh0_w16_avx2(
uint8_t *dst, uint32_t dst_stride, const CONV_BUF_TYPE *src0,
uint32_t src0_stride, const CONV_BUF_TYPE *src1, uint32_t src1_stride,
const uint8_t *mask, uint32_t mask_stride, int h,
const __m256i *round_offset, int shift) {
const __m256i v_maxval = _mm256_set1_epi16(AOM_BLEND_A64_MAX_ALPHA);
for (int i = 0; i < h; ++i) {
const __m128i m = xx_loadu_128(mask);
const __m256i m0 = _mm256_cvtepu8_epi16(m);
blend_a64_d16_mask_w16_avx2(dst, src0, src1, &m0, round_offset, &v_maxval,
shift);
mask += mask_stride;
dst += dst_stride;
src0 += src0_stride;
src1 += src1_stride;
}
}
static INLINE void lowbd_blend_a64_d16_mask_subw0_subh0_w32_avx2(
uint8_t *dst, uint32_t dst_stride, const CONV_BUF_TYPE *src0,
uint32_t src0_stride, const CONV_BUF_TYPE *src1, uint32_t src1_stride,
const uint8_t *mask, uint32_t mask_stride, int h, int w,
const __m256i *round_offset, int shift) {
const __m256i v_maxval = _mm256_set1_epi16(AOM_BLEND_A64_MAX_ALPHA);
for (int i = 0; i < h; ++i) {
for (int j = 0; j < w; j += 32) {
const __m256i m = yy_loadu_256(mask + j);
const __m256i m0 = _mm256_cvtepu8_epi16(_mm256_castsi256_si128(m));
const __m256i m1 = _mm256_cvtepu8_epi16(_mm256_extracti128_si256(m, 1));
blend_a64_d16_mask_w32_avx2(dst + j, src0 + j, src1 + j, &m0, &m1,
round_offset, &v_maxval, shift);
}
mask += mask_stride;
dst += dst_stride;
src0 += src0_stride;
src1 += src1_stride;
}
}
static INLINE void lowbd_blend_a64_d16_mask_subw1_subh1_w16_avx2(
uint8_t *dst, uint32_t dst_stride, const CONV_BUF_TYPE *src0,
uint32_t src0_stride, const CONV_BUF_TYPE *src1, uint32_t src1_stride,
const uint8_t *mask, uint32_t mask_stride, int h,
const __m256i *round_offset, int shift) {
const __m256i v_maxval = _mm256_set1_epi16(AOM_BLEND_A64_MAX_ALPHA);
const __m256i one_b = _mm256_set1_epi8(1);
const __m256i two_w = _mm256_set1_epi16(2);
for (int i = 0; i < h; ++i) {
const __m256i m_i00 = yy_loadu_256(mask);
const __m256i m_i10 = yy_loadu_256(mask + mask_stride);
const __m256i m0_ac = _mm256_adds_epu8(m_i00, m_i10);
const __m256i m0_acbd = _mm256_maddubs_epi16(m0_ac, one_b);
const __m256i m0 = _mm256_srli_epi16(_mm256_add_epi16(m0_acbd, two_w), 2);
blend_a64_d16_mask_w16_avx2(dst, src0, src1, &m0, round_offset, &v_maxval,
shift);
mask += mask_stride << 1;
dst += dst_stride;
src0 += src0_stride;
src1 += src1_stride;
}
}
static INLINE void lowbd_blend_a64_d16_mask_subw1_subh1_w32_avx2(
uint8_t *dst, uint32_t dst_stride, const CONV_BUF_TYPE *src0,
uint32_t src0_stride, const CONV_BUF_TYPE *src1, uint32_t src1_stride,
const uint8_t *mask, uint32_t mask_stride, int h, int w,
const __m256i *round_offset, int shift) {
const __m256i v_maxval = _mm256_set1_epi16(AOM_BLEND_A64_MAX_ALPHA);
const __m256i one_b = _mm256_set1_epi8(1);
const __m256i two_w = _mm256_set1_epi16(2);
for (int i = 0; i < h; ++i) {
for (int j = 0; j < w; j += 32) {
const __m256i m_i00 = yy_loadu_256(mask + 2 * j);
const __m256i m_i01 = yy_loadu_256(mask + 2 * j + 32);
const __m256i m_i10 = yy_loadu_256(mask + mask_stride + 2 * j);
const __m256i m_i11 = yy_loadu_256(mask + mask_stride + 2 * j + 32);
const __m256i m0_ac = _mm256_adds_epu8(m_i00, m_i10);
const __m256i m1_ac = _mm256_adds_epu8(m_i01, m_i11);
const __m256i m0_acbd = _mm256_maddubs_epi16(m0_ac, one_b);
const __m256i m1_acbd = _mm256_maddubs_epi16(m1_ac, one_b);
const __m256i m0 = _mm256_srli_epi16(_mm256_add_epi16(m0_acbd, two_w), 2);
const __m256i m1 = _mm256_srli_epi16(_mm256_add_epi16(m1_acbd, two_w), 2);
blend_a64_d16_mask_w32_avx2(dst + j, src0 + j, src1 + j, &m0, &m1,
round_offset, &v_maxval, shift);
}
mask += mask_stride << 1;
dst += dst_stride;
src0 += src0_stride;
src1 += src1_stride;
}
}
static INLINE void lowbd_blend_a64_d16_mask_subw1_subh0_w16_avx2(
uint8_t *dst, uint32_t dst_stride, const CONV_BUF_TYPE *src0,
uint32_t src0_stride, const CONV_BUF_TYPE *src1, uint32_t src1_stride,
const uint8_t *mask, uint32_t mask_stride, int h, int w,
const __m256i *round_offset, int shift) {
const __m256i v_maxval = _mm256_set1_epi16(AOM_BLEND_A64_MAX_ALPHA);
const __m256i one_b = _mm256_set1_epi8(1);
const __m256i zeros = _mm256_setzero_si256();
for (int i = 0; i < h; ++i) {
for (int j = 0; j < w; j += 16) {
const __m256i m_i00 = yy_loadu_256(mask + 2 * j);
const __m256i m0_ac = _mm256_maddubs_epi16(m_i00, one_b);
const __m256i m0 = _mm256_avg_epu16(m0_ac, zeros);
blend_a64_d16_mask_w16_avx2(dst + j, src0 + j, src1 + j, &m0,
round_offset, &v_maxval, shift);
}
mask += mask_stride;
dst += dst_stride;
src0 += src0_stride;
src1 += src1_stride;
}
}
static INLINE void lowbd_blend_a64_d16_mask_subw1_subh0_w32_avx2(
uint8_t *dst, uint32_t dst_stride, const CONV_BUF_TYPE *src0,
uint32_t src0_stride, const CONV_BUF_TYPE *src1, uint32_t src1_stride,
const uint8_t *mask, uint32_t mask_stride, int h, int w,
const __m256i *round_offset, int shift) {
const __m256i v_maxval = _mm256_set1_epi16(AOM_BLEND_A64_MAX_ALPHA);
const __m256i one_b = _mm256_set1_epi8(1);
const __m256i zeros = _mm256_setzero_si256();
for (int i = 0; i < h; ++i) {
for (int j = 0; j < w; j += 32) {
const __m256i m_i00 = yy_loadu_256(mask + 2 * j);
const __m256i m_i01 = yy_loadu_256(mask + 2 * j + 32);
const __m256i m0_ac = _mm256_maddubs_epi16(m_i00, one_b);
const __m256i m1_ac = _mm256_maddubs_epi16(m_i01, one_b);
const __m256i m0 = _mm256_avg_epu16(m0_ac, zeros);
const __m256i m1 = _mm256_avg_epu16(m1_ac, zeros);
blend_a64_d16_mask_w32_avx2(dst + j, src0 + j, src1 + j, &m0, &m1,
round_offset, &v_maxval, shift);
}
mask += mask_stride;
dst += dst_stride;
src0 += src0_stride;
src1 += src1_stride;
}
}
static INLINE void lowbd_blend_a64_d16_mask_subw0_subh1_w16_avx2(
uint8_t *dst, uint32_t dst_stride, const CONV_BUF_TYPE *src0,
uint32_t src0_stride, const CONV_BUF_TYPE *src1, uint32_t src1_stride,
const uint8_t *mask, uint32_t mask_stride, int h, int w,
const __m256i *round_offset, int shift) {
const __m256i v_maxval = _mm256_set1_epi16(AOM_BLEND_A64_MAX_ALPHA);
const __m128i zeros = _mm_setzero_si128();
for (int i = 0; i < h; ++i) {
for (int j = 0; j < w; j += 16) {
const __m128i m_i00 = xx_loadu_128(mask + j);
const __m128i m_i10 = xx_loadu_128(mask + mask_stride + j);
const __m128i m_ac = _mm_avg_epu8(_mm_adds_epu8(m_i00, m_i10), zeros);
const __m256i m0 = _mm256_cvtepu8_epi16(m_ac);
blend_a64_d16_mask_w16_avx2(dst + j, src0 + j, src1 + j, &m0,
round_offset, &v_maxval, shift);
}
mask += mask_stride << 1;
dst += dst_stride;
src0 += src0_stride;
src1 += src1_stride;
}
}
static INLINE void lowbd_blend_a64_d16_mask_subw0_subh1_w32_avx2(
uint8_t *dst, uint32_t dst_stride, const CONV_BUF_TYPE *src0,
uint32_t src0_stride, const CONV_BUF_TYPE *src1, uint32_t src1_stride,
const uint8_t *mask, uint32_t mask_stride, int h, int w,
const __m256i *round_offset, int shift) {
const __m256i v_maxval = _mm256_set1_epi16(AOM_BLEND_A64_MAX_ALPHA);
const __m256i zeros = _mm256_setzero_si256();
for (int i = 0; i < h; ++i) {
for (int j = 0; j < w; j += 32) {
const __m256i m_i00 = yy_loadu_256(mask + j);
const __m256i m_i10 = yy_loadu_256(mask + mask_stride + j);
const __m256i m_ac =
_mm256_avg_epu8(_mm256_adds_epu8(m_i00, m_i10), zeros);
const __m256i m0 = _mm256_cvtepu8_epi16(_mm256_castsi256_si128(m_ac));
const __m256i m1 =
_mm256_cvtepu8_epi16(_mm256_extracti128_si256(m_ac, 1));
blend_a64_d16_mask_w32_avx2(dst + j, src0 + j, src1 + j, &m0, &m1,
round_offset, &v_maxval, shift);
}
mask += mask_stride << 1;
dst += dst_stride;
src0 += src0_stride;
src1 += src1_stride;
}
}
void aom_lowbd_blend_a64_d16_mask_avx2(
uint8_t *dst, uint32_t dst_stride, const CONV_BUF_TYPE *src0,
uint32_t src0_stride, const CONV_BUF_TYPE *src1, uint32_t src1_stride,
const uint8_t *mask, uint32_t mask_stride, int w, int h, int subw, int subh,
ConvolveParams *conv_params) {
const int bd = 8;
const int round_bits =
2 * FILTER_BITS - conv_params->round_0 - conv_params->round_1;
const int round_offset =
((1 << (round_bits + bd)) + (1 << (round_bits + bd - 1)) -
(1 << (round_bits - 1)))
<< AOM_BLEND_A64_ROUND_BITS;
const int shift = round_bits + AOM_BLEND_A64_ROUND_BITS;
assert(IMPLIES((void *)src0 == dst, src0_stride == dst_stride));
assert(IMPLIES((void *)src1 == dst, src1_stride == dst_stride));
assert(h >= 4);
assert(w >= 4);
assert(IS_POWER_OF_TWO(h));
assert(IS_POWER_OF_TWO(w));
const __m128i v_round_offset = _mm_set1_epi32(round_offset);
const __m256i y_round_offset = _mm256_set1_epi32(round_offset);
if (subw == 0 && subh == 0) {
switch (w) {
case 4:
aom_lowbd_blend_a64_d16_mask_subw0_subh0_w4_sse4_1(
dst, dst_stride, src0, src0_stride, src1, src1_stride, mask,
mask_stride, h, &v_round_offset, shift);
break;
case 8:
aom_lowbd_blend_a64_d16_mask_subw0_subh0_w8_sse4_1(
dst, dst_stride, src0, src0_stride, src1, src1_stride, mask,
mask_stride, h, &v_round_offset, shift);
break;
case 16:
lowbd_blend_a64_d16_mask_subw0_subh0_w16_avx2(
dst, dst_stride, src0, src0_stride, src1, src1_stride, mask,
mask_stride, h, &y_round_offset, shift);
break;
default:
lowbd_blend_a64_d16_mask_subw0_subh0_w32_avx2(
dst, dst_stride, src0, src0_stride, src1, src1_stride, mask,
mask_stride, h, w, &y_round_offset, shift);
break;
}
} else if (subw == 1 && subh == 1) {
switch (w) {
case 4:
aom_lowbd_blend_a64_d16_mask_subw1_subh1_w4_sse4_1(
dst, dst_stride, src0, src0_stride, src1, src1_stride, mask,
mask_stride, h, &v_round_offset, shift);
break;
case 8:
aom_lowbd_blend_a64_d16_mask_subw1_subh1_w8_sse4_1(
dst, dst_stride, src0, src0_stride, src1, src1_stride, mask,
mask_stride, h, &v_round_offset, shift);
break;
case 16:
lowbd_blend_a64_d16_mask_subw1_subh1_w16_avx2(
dst, dst_stride, src0, src0_stride, src1, src1_stride, mask,
mask_stride, h, &y_round_offset, shift);
break;
default:
lowbd_blend_a64_d16_mask_subw1_subh1_w32_avx2(
dst, dst_stride, src0, src0_stride, src1, src1_stride, mask,
mask_stride, h, w, &y_round_offset, shift);
break;
}
} else if (subw == 1 && subh == 0) {
switch (w) {
case 4:
aom_lowbd_blend_a64_d16_mask_subw1_subh0_w4_sse4_1(
dst, dst_stride, src0, src0_stride, src1, src1_stride, mask,
mask_stride, h, &v_round_offset, shift);
break;
case 8:
aom_lowbd_blend_a64_d16_mask_subw1_subh0_w8_sse4_1(
dst, dst_stride, src0, src0_stride, src1, src1_stride, mask,
mask_stride, h, &v_round_offset, shift);
break;
case 16:
lowbd_blend_a64_d16_mask_subw1_subh0_w16_avx2(
dst, dst_stride, src0, src0_stride, src1, src1_stride, mask,
mask_stride, h, w, &y_round_offset, shift);
break;
default:
lowbd_blend_a64_d16_mask_subw1_subh0_w32_avx2(
dst, dst_stride, src0, src0_stride, src1, src1_stride, mask,
mask_stride, h, w, &y_round_offset, shift);
break;
}
} else {
switch (w) {
case 4:
aom_lowbd_blend_a64_d16_mask_subw0_subh1_w4_sse4_1(
dst, dst_stride, src0, src0_stride, src1, src1_stride, mask,
mask_stride, h, &v_round_offset, shift);
break;
case 8:
aom_lowbd_blend_a64_d16_mask_subw0_subh1_w8_sse4_1(
dst, dst_stride, src0, src0_stride, src1, src1_stride, mask,
mask_stride, h, &v_round_offset, shift);
break;
case 16:
lowbd_blend_a64_d16_mask_subw0_subh1_w16_avx2(
dst, dst_stride, src0, src0_stride, src1, src1_stride, mask,
mask_stride, h, w, &y_round_offset, shift);
break;
default:
lowbd_blend_a64_d16_mask_subw0_subh1_w32_avx2(
dst, dst_stride, src0, src0_stride, src1, src1_stride, mask,
mask_stride, h, w, &y_round_offset, shift);
break;
}
}
}
static INLINE __m256i blend_16_u8_avx2(const uint8_t *src0, const uint8_t *src1,
const __m256i *v_m0_b,
const __m256i *v_m1_b,
const int32_t bits) {
const __m256i v_s0_b = _mm256_castsi128_si256(xx_loadu_128(src0));
const __m256i v_s1_b = _mm256_castsi128_si256(xx_loadu_128(src1));
const __m256i v_s0_s_b = _mm256_permute4x64_epi64(v_s0_b, 0xd8);
const __m256i v_s1_s_b = _mm256_permute4x64_epi64(v_s1_b, 0xd8);
const __m256i v_p0_w =
_mm256_maddubs_epi16(_mm256_unpacklo_epi8(v_s0_s_b, v_s1_s_b),
_mm256_unpacklo_epi8(*v_m0_b, *v_m1_b));
const __m256i v_res0_w = yy_roundn_epu16(v_p0_w, bits);
const __m256i v_res_b = _mm256_packus_epi16(v_res0_w, v_res0_w);
const __m256i v_res = _mm256_permute4x64_epi64(v_res_b, 0xd8);
return v_res;
}
static INLINE __m256i blend_32_u8_avx2(const uint8_t *src0, const uint8_t *src1,
const __m256i *v_m0_b,
const __m256i *v_m1_b,
const int32_t bits) {
const __m256i v_s0_b = yy_loadu_256(src0);
const __m256i v_s1_b = yy_loadu_256(src1);
const __m256i v_p0_w =
_mm256_maddubs_epi16(_mm256_unpacklo_epi8(v_s0_b, v_s1_b),
_mm256_unpacklo_epi8(*v_m0_b, *v_m1_b));
const __m256i v_p1_w =
_mm256_maddubs_epi16(_mm256_unpackhi_epi8(v_s0_b, v_s1_b),
_mm256_unpackhi_epi8(*v_m0_b, *v_m1_b));
const __m256i v_res0_w = yy_roundn_epu16(v_p0_w, bits);
const __m256i v_res1_w = yy_roundn_epu16(v_p1_w, bits);
const __m256i v_res = _mm256_packus_epi16(v_res0_w, v_res1_w);
return v_res;
}
static INLINE void blend_a64_mask_sx_sy_w16_avx2(
uint8_t *dst, uint32_t dst_stride, const uint8_t *src0,
uint32_t src0_stride, const uint8_t *src1, uint32_t src1_stride,
const uint8_t *mask, uint32_t mask_stride, int h) {
const __m256i v_zmask_b = _mm256_set1_epi16(0xFF);
const __m256i v_maxval_b = _mm256_set1_epi8(AOM_BLEND_A64_MAX_ALPHA);
do {
const __m256i v_ral_b = yy_loadu_256(mask);
const __m256i v_rbl_b = yy_loadu_256(mask + mask_stride);
const __m256i v_rvsl_b = _mm256_add_epi8(v_ral_b, v_rbl_b);
const __m256i v_rvsal_w = _mm256_and_si256(v_rvsl_b, v_zmask_b);
const __m256i v_rvsbl_w =
_mm256_and_si256(_mm256_srli_si256(v_rvsl_b, 1), v_zmask_b);
const __m256i v_rsl_w = _mm256_add_epi16(v_rvsal_w, v_rvsbl_w);
const __m256i v_m0_w = yy_roundn_epu16(v_rsl_w, 2);
const __m256i v_m0_b = _mm256_packus_epi16(v_m0_w, v_m0_w);
const __m256i v_m1_b = _mm256_sub_epi8(v_maxval_b, v_m0_b);
const __m256i y_res_b = blend_16_u8_avx2(src0, src1, &v_m0_b, &v_m1_b,
AOM_BLEND_A64_ROUND_BITS);
xx_storeu_128(dst, _mm256_castsi256_si128(y_res_b));
dst += dst_stride;
src0 += src0_stride;
src1 += src1_stride;
mask += 2 * mask_stride;
} while (--h);
}
static INLINE void blend_a64_mask_sx_sy_w32n_avx2(
uint8_t *dst, uint32_t dst_stride, const uint8_t *src0,
uint32_t src0_stride, const uint8_t *src1, uint32_t src1_stride,
const uint8_t *mask, uint32_t mask_stride, int w, int h) {
const __m256i v_maxval_b = _mm256_set1_epi8(AOM_BLEND_A64_MAX_ALPHA);
const __m256i v_zmask_b = _mm256_set1_epi16(0xFF);
do {
int c;
for (c = 0; c < w; c += 32) {
const __m256i v_ral_b = yy_loadu_256(mask + 2 * c);
const __m256i v_rah_b = yy_loadu_256(mask + 2 * c + 32);
const __m256i v_rbl_b = yy_loadu_256(mask + mask_stride + 2 * c);
const __m256i v_rbh_b = yy_loadu_256(mask + mask_stride + 2 * c + 32);
const __m256i v_rvsl_b = _mm256_add_epi8(v_ral_b, v_rbl_b);
const __m256i v_rvsh_b = _mm256_add_epi8(v_rah_b, v_rbh_b);
const __m256i v_rvsal_w = _mm256_and_si256(v_rvsl_b, v_zmask_b);
const __m256i v_rvsah_w = _mm256_and_si256(v_rvsh_b, v_zmask_b);
const __m256i v_rvsbl_w =
_mm256_and_si256(_mm256_srli_si256(v_rvsl_b, 1), v_zmask_b);
const __m256i v_rvsbh_w =
_mm256_and_si256(_mm256_srli_si256(v_rvsh_b, 1), v_zmask_b);
const __m256i v_rsl_w = _mm256_add_epi16(v_rvsal_w, v_rvsbl_w);
const __m256i v_rsh_w = _mm256_add_epi16(v_rvsah_w, v_rvsbh_w);
const __m256i v_m0l_w = yy_roundn_epu16(v_rsl_w, 2);
const __m256i v_m0h_w = yy_roundn_epu16(v_rsh_w, 2);
const __m256i v_m0_b =
_mm256_permute4x64_epi64(_mm256_packus_epi16(v_m0l_w, v_m0h_w), 0xd8);
const __m256i v_m1_b = _mm256_sub_epi8(v_maxval_b, v_m0_b);
const __m256i v_res_b = blend_32_u8_avx2(
src0 + c, src1 + c, &v_m0_b, &v_m1_b, AOM_BLEND_A64_ROUND_BITS);
yy_storeu_256(dst + c, v_res_b);
}
dst += dst_stride;
src0 += src0_stride;
src1 += src1_stride;
mask += 2 * mask_stride;
} while (--h);
}
static INLINE void blend_a64_mask_sx_sy_avx2(
uint8_t *dst, uint32_t dst_stride, const uint8_t *src0,
uint32_t src0_stride, const uint8_t *src1, uint32_t src1_stride,
const uint8_t *mask, uint32_t mask_stride, int w, int h) {
const __m128i v_shuffle_b = xx_loadu_128(g_blend_a64_mask_shuffle);
const __m128i v_maxval_b = _mm_set1_epi8(AOM_BLEND_A64_MAX_ALPHA);
const __m128i _r = _mm_set1_epi16(1 << (15 - AOM_BLEND_A64_ROUND_BITS));
switch (w) {
case 4:
do {
const __m128i v_ra_b = xx_loadl_64(mask);
const __m128i v_rb_b = xx_loadl_64(mask + mask_stride);
const __m128i v_rvs_b = _mm_add_epi8(v_ra_b, v_rb_b);
const __m128i v_r_s_b = _mm_shuffle_epi8(v_rvs_b, v_shuffle_b);
const __m128i v_r0_s_w = _mm_cvtepu8_epi16(v_r_s_b);
const __m128i v_r1_s_w = _mm_cvtepu8_epi16(_mm_srli_si128(v_r_s_b, 8));
const __m128i v_rs_w = _mm_add_epi16(v_r0_s_w, v_r1_s_w);
const __m128i v_m0_w = xx_roundn_epu16(v_rs_w, 2);
const __m128i v_m0_b = _mm_packus_epi16(v_m0_w, v_m0_w);
const __m128i v_m1_b = _mm_sub_epi8(v_maxval_b, v_m0_b);
const __m128i v_res_b = blend_4_u8(src0, src1, &v_m0_b, &v_m1_b, &_r);
xx_storel_32(dst, v_res_b);
dst += dst_stride;
src0 += src0_stride;
src1 += src1_stride;
mask += 2 * mask_stride;
} while (--h);
break;
case 8:
do {
const __m128i v_ra_b = xx_loadu_128(mask);
const __m128i v_rb_b = xx_loadu_128(mask + mask_stride);
const __m128i v_rvs_b = _mm_add_epi8(v_ra_b, v_rb_b);
const __m128i v_r_s_b = _mm_shuffle_epi8(v_rvs_b, v_shuffle_b);
const __m128i v_r0_s_w = _mm_cvtepu8_epi16(v_r_s_b);
const __m128i v_r1_s_w = _mm_cvtepu8_epi16(_mm_srli_si128(v_r_s_b, 8));
const __m128i v_rs_w = _mm_add_epi16(v_r0_s_w, v_r1_s_w);
const __m128i v_m0_w = xx_roundn_epu16(v_rs_w, 2);
const __m128i v_m0_b = _mm_packus_epi16(v_m0_w, v_m0_w);
const __m128i v_m1_b = _mm_sub_epi8(v_maxval_b, v_m0_b);
const __m128i v_res_b = blend_8_u8(src0, src1, &v_m0_b, &v_m1_b, &_r);
xx_storel_64(dst, v_res_b);
dst += dst_stride;
src0 += src0_stride;
src1 += src1_stride;
mask += 2 * mask_stride;
} while (--h);
break;
case 16:
blend_a64_mask_sx_sy_w16_avx2(dst, dst_stride, src0, src0_stride, src1,
src1_stride, mask, mask_stride, h);
break;
default:
blend_a64_mask_sx_sy_w32n_avx2(dst, dst_stride, src0, src0_stride, src1,
src1_stride, mask, mask_stride, w, h);
break;
}
}
static INLINE void blend_a64_mask_sx_w16_avx2(
uint8_t *dst, uint32_t dst_stride, const uint8_t *src0,
uint32_t src0_stride, const uint8_t *src1, uint32_t src1_stride,
const uint8_t *mask, uint32_t mask_stride, int h) {
const __m256i v_maxval_b = _mm256_set1_epi8(AOM_BLEND_A64_MAX_ALPHA);
const __m256i v_zmask_b = _mm256_set1_epi16(0xff);
do {
const __m256i v_rl_b = yy_loadu_256(mask);
const __m256i v_al_b =
_mm256_avg_epu8(v_rl_b, _mm256_srli_si256(v_rl_b, 1));
const __m256i v_m0_w = _mm256_and_si256(v_al_b, v_zmask_b);
const __m256i v_m0_b = _mm256_packus_epi16(v_m0_w, _mm256_setzero_si256());
const __m256i v_m1_b = _mm256_sub_epi8(v_maxval_b, v_m0_b);
const __m256i v_res_b = blend_16_u8_avx2(src0, src1, &v_m0_b, &v_m1_b,
AOM_BLEND_A64_ROUND_BITS);
xx_storeu_128(dst, _mm256_castsi256_si128(v_res_b));
dst += dst_stride;
src0 += src0_stride;
src1 += src1_stride;
mask += mask_stride;
} while (--h);
}
static INLINE void blend_a64_mask_sx_w32n_avx2(
uint8_t *dst, uint32_t dst_stride, const uint8_t *src0,
uint32_t src0_stride, const uint8_t *src1, uint32_t src1_stride,
const uint8_t *mask, uint32_t mask_stride, int w, int h) {
const __m256i v_shuffle_b = yy_loadu_256(g_blend_a64_mask_shuffle);
const __m256i v_maxval_b = _mm256_set1_epi8(AOM_BLEND_A64_MAX_ALPHA);
do {
int c;
for (c = 0; c < w; c += 32) {
const __m256i v_r0_b = yy_loadu_256(mask + 2 * c);
const __m256i v_r1_b = yy_loadu_256(mask + 2 * c + 32);
const __m256i v_r0_s_b = _mm256_shuffle_epi8(v_r0_b, v_shuffle_b);
const __m256i v_r1_s_b = _mm256_shuffle_epi8(v_r1_b, v_shuffle_b);
const __m256i v_al_b =
_mm256_avg_epu8(v_r0_s_b, _mm256_srli_si256(v_r0_s_b, 8));
const __m256i v_ah_b =
_mm256_avg_epu8(v_r1_s_b, _mm256_srli_si256(v_r1_s_b, 8));
const __m256i v_m0_b =
_mm256_permute4x64_epi64(_mm256_unpacklo_epi64(v_al_b, v_ah_b), 0xd8);
const __m256i v_m1_b = _mm256_sub_epi8(v_maxval_b, v_m0_b);
const __m256i v_res_b = blend_32_u8_avx2(
src0 + c, src1 + c, &v_m0_b, &v_m1_b, AOM_BLEND_A64_ROUND_BITS);
yy_storeu_256(dst + c, v_res_b);
}
dst += dst_stride;
src0 += src0_stride;
src1 += src1_stride;
mask += mask_stride;
} while (--h);
}
static INLINE void blend_a64_mask_sx_avx2(
uint8_t *dst, uint32_t dst_stride, const uint8_t *src0,
uint32_t src0_stride, const uint8_t *src1, uint32_t src1_stride,
const uint8_t *mask, uint32_t mask_stride, int w, int h) {
const __m128i v_shuffle_b = xx_loadu_128(g_blend_a64_mask_shuffle);
const __m128i v_maxval_b = _mm_set1_epi8(AOM_BLEND_A64_MAX_ALPHA);
const __m128i _r = _mm_set1_epi16(1 << (15 - AOM_BLEND_A64_ROUND_BITS));
switch (w) {
case 4:
do {
const __m128i v_r_b = xx_loadl_64(mask);
const __m128i v_r0_s_b = _mm_shuffle_epi8(v_r_b, v_shuffle_b);
const __m128i v_r_lo_b = _mm_unpacklo_epi64(v_r0_s_b, v_r0_s_b);
const __m128i v_r_hi_b = _mm_unpackhi_epi64(v_r0_s_b, v_r0_s_b);
const __m128i v_m0_b = _mm_avg_epu8(v_r_lo_b, v_r_hi_b);
const __m128i v_m1_b = _mm_sub_epi8(v_maxval_b, v_m0_b);
const __m128i v_res_b = blend_4_u8(src0, src1, &v_m0_b, &v_m1_b, &_r);
xx_storel_32(dst, v_res_b);
dst += dst_stride;
src0 += src0_stride;
src1 += src1_stride;
mask += mask_stride;
} while (--h);
break;
case 8:
do {
const __m128i v_r_b = xx_loadu_128(mask);
const __m128i v_r0_s_b = _mm_shuffle_epi8(v_r_b, v_shuffle_b);
const __m128i v_r_lo_b = _mm_unpacklo_epi64(v_r0_s_b, v_r0_s_b);
const __m128i v_r_hi_b = _mm_unpackhi_epi64(v_r0_s_b, v_r0_s_b);
const __m128i v_m0_b = _mm_avg_epu8(v_r_lo_b, v_r_hi_b);
const __m128i v_m1_b = _mm_sub_epi8(v_maxval_b, v_m0_b);
const __m128i v_res_b = blend_8_u8(src0, src1, &v_m0_b, &v_m1_b, &_r);
xx_storel_64(dst, v_res_b);
dst += dst_stride;
src0 += src0_stride;
src1 += src1_stride;
mask += mask_stride;
} while (--h);
break;
case 16:
blend_a64_mask_sx_w16_avx2(dst, dst_stride, src0, src0_stride, src1,
src1_stride, mask, mask_stride, h);
break;
default:
blend_a64_mask_sx_w32n_avx2(dst, dst_stride, src0, src0_stride, src1,
src1_stride, mask, mask_stride, w, h);
break;
}
}
static INLINE void blend_a64_mask_sy_w16_avx2(
uint8_t *dst, uint32_t dst_stride, const uint8_t *src0,
uint32_t src0_stride, const uint8_t *src1, uint32_t src1_stride,
const uint8_t *mask, uint32_t mask_stride, int h) {
const __m128i _r = _mm_set1_epi16(1 << (15 - AOM_BLEND_A64_ROUND_BITS));
const __m128i v_maxval_b = _mm_set1_epi8(AOM_BLEND_A64_MAX_ALPHA);
do {
const __m128i v_ra_b = xx_loadu_128(mask);
const __m128i v_rb_b = xx_loadu_128(mask + mask_stride);
const __m128i v_m0_b = _mm_avg_epu8(v_ra_b, v_rb_b);
const __m128i v_m1_b = _mm_sub_epi16(v_maxval_b, v_m0_b);
const __m128i v_res_b = blend_16_u8(src0, src1, &v_m0_b, &v_m1_b, &_r);
xx_storeu_128(dst, v_res_b);
dst += dst_stride;
src0 += src0_stride;
src1 += src1_stride;
mask += 2 * mask_stride;
} while (--h);
}
static INLINE void blend_a64_mask_sy_w32n_avx2(
uint8_t *dst, uint32_t dst_stride, const uint8_t *src0,
uint32_t src0_stride, const uint8_t *src1, uint32_t src1_stride,
const uint8_t *mask, uint32_t mask_stride, int w, int h) {
const __m256i v_maxval_b = _mm256_set1_epi8(AOM_BLEND_A64_MAX_ALPHA);
do {
int c;
for (c = 0; c < w; c += 32) {
const __m256i v_ra_b = yy_loadu_256(mask + c);
const __m256i v_rb_b = yy_loadu_256(mask + c + mask_stride);
const __m256i v_m0_b = _mm256_avg_epu8(v_ra_b, v_rb_b);
const __m256i v_m1_b = _mm256_sub_epi8(v_maxval_b, v_m0_b);
const __m256i v_res_b = blend_32_u8_avx2(
src0 + c, src1 + c, &v_m0_b, &v_m1_b, AOM_BLEND_A64_ROUND_BITS);
yy_storeu_256(dst + c, v_res_b);
}
dst += dst_stride;
src0 += src0_stride;
src1 += src1_stride;
mask += 2 * mask_stride;
} while (--h);
}
static INLINE void blend_a64_mask_sy_avx2(
uint8_t *dst, uint32_t dst_stride, const uint8_t *src0,
uint32_t src0_stride, const uint8_t *src1, uint32_t src1_stride,
const uint8_t *mask, uint32_t mask_stride, int w, int h) {
const __m128i _r = _mm_set1_epi16(1 << (15 - AOM_BLEND_A64_ROUND_BITS));
const __m128i v_maxval_b = _mm_set1_epi8(AOM_BLEND_A64_MAX_ALPHA);
switch (w) {
case 4:
do {
const __m128i v_ra_b = xx_loadl_32(mask);
const __m128i v_rb_b = xx_loadl_32(mask + mask_stride);
const __m128i v_m0_b = _mm_avg_epu8(v_ra_b, v_rb_b);
const __m128i v_m1_b = _mm_sub_epi8(v_maxval_b, v_m0_b);
const __m128i v_res_b = blend_4_u8(src0, src1, &v_m0_b, &v_m1_b, &_r);
xx_storel_32(dst, v_res_b);
dst += dst_stride;
src0 += src0_stride;
src1 += src1_stride;
mask += 2 * mask_stride;
} while (--h);
break;
case 8:
do {
const __m128i v_ra_b = xx_loadl_64(mask);
const __m128i v_rb_b = xx_loadl_64(mask + mask_stride);
const __m128i v_m0_b = _mm_avg_epu8(v_ra_b, v_rb_b);
const __m128i v_m1_b = _mm_sub_epi8(v_maxval_b, v_m0_b);
const __m128i v_res_b = blend_8_u8(src0, src1, &v_m0_b, &v_m1_b, &_r);
xx_storel_64(dst, v_res_b);
dst += dst_stride;
src0 += src0_stride;
src1 += src1_stride;
mask += 2 * mask_stride;
} while (--h);
break;
case 16:
blend_a64_mask_sy_w16_avx2(dst, dst_stride, src0, src0_stride, src1,
src1_stride, mask, mask_stride, h);
break;
default:
blend_a64_mask_sy_w32n_avx2(dst, dst_stride, src0, src0_stride, src1,
src1_stride, mask, mask_stride, w, h);
}
}
static INLINE void blend_a64_mask_w32n_avx2(
uint8_t *dst, uint32_t dst_stride, const uint8_t *src0,
uint32_t src0_stride, const uint8_t *src1, uint32_t src1_stride,
const uint8_t *mask, uint32_t mask_stride, int w, int h) {
const __m256i v_maxval_b = _mm256_set1_epi8(AOM_BLEND_A64_MAX_ALPHA);
do {
int c;
for (c = 0; c < w; c += 32) {
const __m256i v_m0_b = yy_loadu_256(mask + c);
const __m256i v_m1_b = _mm256_sub_epi8(v_maxval_b, v_m0_b);
const __m256i v_res_b = blend_32_u8_avx2(
src0 + c, src1 + c, &v_m0_b, &v_m1_b, AOM_BLEND_A64_ROUND_BITS);
yy_storeu_256(dst + c, v_res_b);
}
dst += dst_stride;
src0 += src0_stride;
src1 += src1_stride;
mask += mask_stride;
} while (--h);
}
static INLINE void blend_a64_mask_avx2(
uint8_t *dst, uint32_t dst_stride, const uint8_t *src0,
uint32_t src0_stride, const uint8_t *src1, uint32_t src1_stride,
const uint8_t *mask, uint32_t mask_stride, int w, int h) {
const __m128i v_maxval_b = _mm_set1_epi8(AOM_BLEND_A64_MAX_ALPHA);
const __m128i _r = _mm_set1_epi16(1 << (15 - AOM_BLEND_A64_ROUND_BITS));
switch (w) {
case 4:
do {
const __m128i v_m0_b = xx_loadl_32(mask);
const __m128i v_m1_b = _mm_sub_epi8(v_maxval_b, v_m0_b);
const __m128i v_res_b = blend_4_u8(src0, src1, &v_m0_b, &v_m1_b, &_r);
xx_storel_32(dst, v_res_b);
dst += dst_stride;
src0 += src0_stride;
src1 += src1_stride;
mask += mask_stride;
} while (--h);
break;
case 8:
do {
const __m128i v_m0_b = xx_loadl_64(mask);
const __m128i v_m1_b = _mm_sub_epi8(v_maxval_b, v_m0_b);
const __m128i v_res_b = blend_8_u8(src0, src1, &v_m0_b, &v_m1_b, &_r);
xx_storel_64(dst, v_res_b);
dst += dst_stride;
src0 += src0_stride;
src1 += src1_stride;
mask += mask_stride;
} while (--h);
break;
case 16:
do {
const __m128i v_m0_b = xx_loadu_128(mask);
const __m128i v_m1_b = _mm_sub_epi8(v_maxval_b, v_m0_b);
const __m128i v_res_b = blend_16_u8(src0, src1, &v_m0_b, &v_m1_b, &_r);
xx_storeu_128(dst, v_res_b);
dst += dst_stride;
src0 += src0_stride;
src1 += src1_stride;
mask += mask_stride;
} while (--h);
break;
default:
blend_a64_mask_w32n_avx2(dst, dst_stride, src0, src0_stride, src1,
src1_stride, mask, mask_stride, w, h);
}
}
void aom_blend_a64_mask_avx2(uint8_t *dst, uint32_t dst_stride,
const uint8_t *src0, uint32_t src0_stride,
const uint8_t *src1, uint32_t src1_stride,
const uint8_t *mask, uint32_t mask_stride, int w,
int h, int subx, int suby) {
assert(IMPLIES(src0 == dst, src0_stride == dst_stride));
assert(IMPLIES(src1 == dst, src1_stride == dst_stride));
assert(h >= 1);
assert(w >= 1);
assert(IS_POWER_OF_TWO(h));
assert(IS_POWER_OF_TWO(w));
if (UNLIKELY((h | w) & 3)) { // if (w <= 2 || h <= 2)
aom_blend_a64_mask_c(dst, dst_stride, src0, src0_stride, src1, src1_stride,
mask, mask_stride, w, h, subx, suby);
} else {
if (subx & suby) {
blend_a64_mask_sx_sy_avx2(dst, dst_stride, src0, src0_stride, src1,
src1_stride, mask, mask_stride, w, h);
} else if (subx) {
blend_a64_mask_sx_avx2(dst, dst_stride, src0, src0_stride, src1,
src1_stride, mask, mask_stride, w, h);
} else if (suby) {
blend_a64_mask_sy_avx2(dst, dst_stride, src0, src0_stride, src1,
src1_stride, mask, mask_stride, w, h);
} else {
blend_a64_mask_avx2(dst, dst_stride, src0, src0_stride, src1, src1_stride,
mask, mask_stride, w, h);
}
}
}