hashirama/uvg266
1
0
Fork 0
mirror of https://github.com/ultravideo/uvg266.git synced 2024-11-27 19:24:06 +00:00
Code Issues Actions Projects Releases Wiki Activity

[SIMD] Convert planar and DC mode PDPC loops to AVX2

Browse source
This commit is contained in:
Ari Lemmetti 2021-09-07 00:45:44 +03:00
parent ad35d4a4c8
commit 171b9c60b3
1 changed files with 61 additions and 20 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

81
src/strategies/avx2/intra-avx2.c
View file

@ -973,16 +973,35 @@ static void kvz_pdpc_planar_dc_avx2(
{
assert(mode == 0 || mode == 1); // planar or DC
__m256i shuf_mask_byte = _mm256_setr_epi8(
0, -1, 0, -1, 0, -1, 0, -1,
1, -1, 1, -1, 1, -1, 1, -1,
2, -1, 2, -1, 2, -1, 2, -1,
3, -1, 3, -1, 3, -1, 3, -1
);
__m256i shuf_mask_word = _mm256_setr_epi8(
0, 1, 0, 1, 0, 1, 0, 1,
2, 3, 2, 3, 2, 3, 2, 3,
4, 5, 4, 5, 4, 5, 4, 5,
6, 7, 6, 7, 6, 7, 6, 7
);
// TODO: replace latter log2_width with log2_height
const int scale = ((log2_width - 2 + log2_width - 2 + 2) >> 2);
int16_t w[LCU_WIDTH];
int16_t left[4][4];
int16_t top [4][4];
// Same weights regardless of axis, compute once
for (int i = 0; i < width; ++i) {
w[i] = 32 >> MIN(31, ((i << 1) >> scale));
int16_t w[LCU_WIDTH];
for (int i = 0; i < width; i += 4) {
__m128i base = _mm_set1_epi32(i);
__m128i offs = _mm_setr_epi32(0, 1, 2, 3);
__m128i idxs = _mm_add_epi32(base, offs);
__m128i unclipped = _mm_slli_epi32(idxs, 1);
unclipped = _mm_srli_epi32(unclipped, scale);
__m128i clipped = _mm_min_epi32( _mm_set1_epi32(31), unclipped);
__m128i weights = _mm_srlv_epi32(_mm_set1_epi32(32), clipped);
weights = _mm_packus_epi32(weights, weights);
_mm_storel_epi64((__m128i*)&w[i], weights);
}
// Process in 4x4 blocks
@ -990,21 +1009,43 @@ static void kvz_pdpc_planar_dc_avx2(
for (int y = 0; y < width; y += 4) {
for (int x = 0; x < width; x += 4) {
for (int yy = 0; yy < 4; ++yy) {
for (int xx = 0; xx < 4; ++xx) {
left[yy][xx] = used_ref->left[(y + yy) + 1];
top [yy][xx] = used_ref->top [(x + xx) + 1];
}
}
uint32_t dw_left;
uint32_t dw_top;
memcpy(&dw_left, &used_ref->left[y + 1], sizeof(dw_left));
memcpy(&dw_top , &used_ref->top [x + 1], sizeof(dw_top));
__m256i vleft = _mm256_set1_epi32(dw_left);
__m256i vtop = _mm256_set1_epi32(dw_top);
vleft = _mm256_shuffle_epi8(vleft, shuf_mask_byte);
vtop = _mm256_cvtepu8_epi16(_mm256_castsi256_si128(vtop));
for (int yy = 0; yy < 4; ++yy) {
for (int xx = 0; xx < 4; ++xx) {
dst[(x + xx) + (y + yy) * width] += ((
w[(x + xx)] * (left[yy][xx] - dst[(x + xx) + (y + yy) * width]) +
w[(y + yy)] * (top [yy][xx] - dst[(x + xx) + (y + yy) * width]) +
32) >> 6);
}
}
__m128i vseq = _mm_setr_epi32(0, 1, 2, 3);
__m128i vidx = _mm_slli_epi32(vseq, log2_width);
__m128i vdst = _mm_i32gather_epi32((uint32_t*)(dst + y * width + x), vidx, 1);
__m256i vdst16 = _mm256_cvtepu8_epi16(vdst);
uint64_t quad_wL;
uint64_t quad_wT;
memcpy(&quad_wL, &w[x], sizeof(quad_wL));
memcpy(&quad_wT, &w[y], sizeof(quad_wT));
__m256i vwL = _mm256_set1_epi64x(quad_wL);
__m256i vwT = _mm256_set1_epi64x(quad_wT);
vwT = _mm256_shuffle_epi8(vwT, shuf_mask_word);
__m256i diff_left = _mm256_sub_epi16(vleft, vdst16);
__m256i diff_top = _mm256_sub_epi16(vtop , vdst16);
__m256i prod_left = _mm256_mullo_epi16(vwL, diff_left);
__m256i prod_top = _mm256_mullo_epi16(vwT, diff_top);
__m256i accu = _mm256_add_epi16(prod_left, prod_top);
accu = _mm256_add_epi16(accu, _mm256_set1_epi16(32));
accu = _mm256_srai_epi16(accu, 6);
accu = _mm256_add_epi16(vdst16, accu);
__m128i lo = _mm256_castsi256_si128(accu);
__m128i hi = _mm256_extracti128_si256(accu, 1);
vdst = _mm_packus_epi16(lo, hi);
*(uint32_t*)(dst + (y + 0) * width + x) = _mm_extract_epi32(vdst, 0);
*(uint32_t*)(dst + (y + 1) * width + x) = _mm_extract_epi32(vdst, 1);
*(uint32_t*)(dst + (y + 2) * width + x) = _mm_extract_epi32(vdst, 2);
*(uint32_t*)(dst + (y + 3) * width + x) = _mm_extract_epi32(vdst, 3);
}
}
}