ffmpeg/libavcodec/mips/hevcdsp_msa.c
Shivraj Patil 97f074f134 avcodec/mips: MSA (MIPS-SIMD-Arch) optimizations for HEVC copy and hv mc functions
Incorporated review comment.
Removed "__" from volatile.

Signed-off-by: Shivraj Patil <shivraj.patil@imgtec.com>
Reviewed-by: Nedeljko Babic <Nedeljko.Babic@imgtec.com>
Signed-off-by: Michael Niedermayer <michaelni@gmx.at>
2015-04-24 17:12:04 +02:00

2358 lines
85 KiB
C

/*
* Copyright (c) 2015 Manojkumar Bhosale (Manojkumar.Bhosale@imgtec.com)
*
* This file is part of FFmpeg.
*
* FFmpeg is free software; you can redistribute it and/or
* modify it under the terms of the GNU Lesser General Public
* License as published by the Free Software Foundation; either
* version 2.1 of the License, or (at your option) any later version.
*
* FFmpeg is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
* Lesser General Public License for more details.
*
* You should have received a copy of the GNU Lesser General Public
* License along with FFmpeg; if not, write to the Free Software
* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
*/
#include "libavutil/mips/generic_macros_msa.h"
#include "libavcodec/mips/hevcdsp_mips.h"
#define HEVC_FILT_8TAP_DPADD_W(vec0, vec1, vec2, vec3, \
filt0, filt1, filt2, filt3) \
( { \
v4i32 out; \
\
out = __msa_dotp_s_w((v8i16) (vec0), (v8i16) (filt0)); \
out = __msa_dpadd_s_w(out, (v8i16) (vec1), (v8i16) (filt1)); \
out = __msa_dpadd_s_w(out, (v8i16) (vec2), (v8i16) (filt2)); \
out = __msa_dpadd_s_w(out, (v8i16) (vec3), (v8i16) (filt3)); \
out; \
} )
#define HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3, \
filt0, filt1, filt2, filt3, \
var_in) \
( { \
v8i16 out; \
\
out = __msa_dpadd_s_h((v8i16) (var_in), (v16i8) (vec0), (v16i8) (filt0)); \
out = __msa_dpadd_s_h(out, (v16i8) (vec1), (v16i8) (filt1)); \
out = __msa_dpadd_s_h(out, (v16i8) (vec2), (v16i8) (filt2)); \
out = __msa_dpadd_s_h(out, (v16i8) (vec3), (v16i8) (filt3)); \
out; \
} )
static void hevc_copy_4w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
int32_t height)
{
v16i8 zero = { 0 };
if (2 == height) {
uint64_t out0, out1;
v16i8 src0, src1;
v8i16 input0;
LOAD_2VECS_SB(src, src_stride, src0, src1);
src0 = (v16i8) __msa_ilvr_w((v4i32) src1, (v4i32) src0);
input0 = (v8i16) __msa_ilvr_b(zero, src0);
input0 <<= 6;
out0 = __msa_copy_u_d((v2i64) input0, 0);
out1 = __msa_copy_u_d((v2i64) input0, 1);
STORE_DWORD(dst, out0);
dst += dst_stride;
STORE_DWORD(dst, out1);
} else if (4 == height) {
uint64_t out0, out1, out2, out3;
v16i8 src0, src1, src2, src3;
v8i16 input0, input1;
LOAD_4VECS_SB(src, src_stride, src0, src1, src2, src3);
src0 = (v16i8) __msa_ilvr_w((v4i32) src1, (v4i32) src0);
src1 = (v16i8) __msa_ilvr_w((v4i32) src3, (v4i32) src2);
input0 = (v8i16) __msa_ilvr_b(zero, src0);
input1 = (v8i16) __msa_ilvr_b(zero, src1);
input0 <<= 6;
input1 <<= 6;
out0 = __msa_copy_u_d((v2i64) input0, 0);
out1 = __msa_copy_u_d((v2i64) input0, 1);
out2 = __msa_copy_u_d((v2i64) input1, 0);
out3 = __msa_copy_u_d((v2i64) input1, 1);
STORE_DWORD(dst, out0);
dst += dst_stride;
STORE_DWORD(dst, out1);
dst += dst_stride;
STORE_DWORD(dst, out2);
dst += dst_stride;
STORE_DWORD(dst, out3);
} else if (0 == height % 8) {
v16i8 src0, src1, src2, src3, src4, src5, src6, src7;
uint64_t out0, out1, out2, out3, out4, out5, out6, out7;
v8i16 input0, input1, input2, input3;
uint32_t loop_cnt;
for (loop_cnt = (height >> 3); loop_cnt--;) {
LOAD_8VECS_SB(src, src_stride,
src0, src1, src2, src3, src4, src5, src6, src7);
src += (8 * src_stride);
src0 = (v16i8) __msa_ilvr_w((v4i32) src1, (v4i32) src0);
src1 = (v16i8) __msa_ilvr_w((v4i32) src3, (v4i32) src2);
src2 = (v16i8) __msa_ilvr_w((v4i32) src5, (v4i32) src4);
src3 = (v16i8) __msa_ilvr_w((v4i32) src7, (v4i32) src6);
input0 = (v8i16) __msa_ilvr_b(zero, src0);
input1 = (v8i16) __msa_ilvr_b(zero, src1);
input2 = (v8i16) __msa_ilvr_b(zero, src2);
input3 = (v8i16) __msa_ilvr_b(zero, src3);
input0 <<= 6;
input1 <<= 6;
input2 <<= 6;
input3 <<= 6;
out0 = __msa_copy_u_d((v2i64) input0, 0);
out1 = __msa_copy_u_d((v2i64) input0, 1);
out2 = __msa_copy_u_d((v2i64) input1, 0);
out3 = __msa_copy_u_d((v2i64) input1, 1);
out4 = __msa_copy_u_d((v2i64) input2, 0);
out5 = __msa_copy_u_d((v2i64) input2, 1);
out6 = __msa_copy_u_d((v2i64) input3, 0);
out7 = __msa_copy_u_d((v2i64) input3, 1);
STORE_DWORD(dst, out0);
dst += dst_stride;
STORE_DWORD(dst, out1);
dst += dst_stride;
STORE_DWORD(dst, out2);
dst += dst_stride;
STORE_DWORD(dst, out3);
dst += dst_stride;
STORE_DWORD(dst, out4);
dst += dst_stride;
STORE_DWORD(dst, out5);
dst += dst_stride;
STORE_DWORD(dst, out6);
dst += dst_stride;
STORE_DWORD(dst, out7);
dst += dst_stride;
}
}
}
static void hevc_copy_6w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
int32_t height)
{
uint32_t loop_cnt;
uint64_t out0, out1, out2, out3, out4, out5, out6, out7;
uint32_t out8, out9, out10, out11, out12, out13, out14, out15;
v16i8 zero = { 0 };
v16i8 src0, src1, src2, src3, src4, src5, src6, src7;
v8i16 input0, input1, input2, input3, input4, input5, input6, input7;
for (loop_cnt = (height >> 3); loop_cnt--;) {
LOAD_8VECS_SB(src, src_stride,
src0, src1, src2, src3, src4, src5, src6, src7);
src += (8 * src_stride);
input0 = (v8i16) __msa_ilvr_b(zero, src0);
input1 = (v8i16) __msa_ilvr_b(zero, src1);
input2 = (v8i16) __msa_ilvr_b(zero, src2);
input3 = (v8i16) __msa_ilvr_b(zero, src3);
input4 = (v8i16) __msa_ilvr_b(zero, src4);
input5 = (v8i16) __msa_ilvr_b(zero, src5);
input6 = (v8i16) __msa_ilvr_b(zero, src6);
input7 = (v8i16) __msa_ilvr_b(zero, src7);
input0 <<= 6;
input1 <<= 6;
input2 <<= 6;
input3 <<= 6;
input4 <<= 6;
input5 <<= 6;
input6 <<= 6;
input7 <<= 6;
out0 = __msa_copy_u_d((v2i64) input0, 0);
out1 = __msa_copy_u_d((v2i64) input1, 0);
out2 = __msa_copy_u_d((v2i64) input2, 0);
out3 = __msa_copy_u_d((v2i64) input3, 0);
out4 = __msa_copy_u_d((v2i64) input4, 0);
out5 = __msa_copy_u_d((v2i64) input5, 0);
out6 = __msa_copy_u_d((v2i64) input6, 0);
out7 = __msa_copy_u_d((v2i64) input7, 0);
out8 = __msa_copy_u_w((v4i32) input0, 2);
out9 = __msa_copy_u_w((v4i32) input1, 2);
out10 = __msa_copy_u_w((v4i32) input2, 2);
out11 = __msa_copy_u_w((v4i32) input3, 2);
out12 = __msa_copy_u_w((v4i32) input4, 2);
out13 = __msa_copy_u_w((v4i32) input5, 2);
out14 = __msa_copy_u_w((v4i32) input6, 2);
out15 = __msa_copy_u_w((v4i32) input7, 2);
STORE_DWORD(dst, out0);
STORE_WORD(dst + 4, out8);
dst += dst_stride;
STORE_DWORD(dst, out1);
STORE_WORD(dst + 4, out9);
dst += dst_stride;
STORE_DWORD(dst, out2);
STORE_WORD(dst + 4, out10);
dst += dst_stride;
STORE_DWORD(dst, out3);
STORE_WORD(dst + 4, out11);
dst += dst_stride;
STORE_DWORD(dst, out4);
STORE_WORD(dst + 4, out12);
dst += dst_stride;
STORE_DWORD(dst, out5);
STORE_WORD(dst + 4, out13);
dst += dst_stride;
STORE_DWORD(dst, out6);
STORE_WORD(dst + 4, out14);
dst += dst_stride;
STORE_DWORD(dst, out7);
STORE_WORD(dst + 4, out15);
dst += dst_stride;
}
}
static void hevc_copy_8w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
int32_t height)
{
v16i8 zero = { 0 };
if (2 == height) {
v16i8 src0, src1;
v8i16 input0, input1;
LOAD_2VECS_SB(src, src_stride, src0, src1);
input0 = (v8i16) __msa_ilvr_b(zero, src0);
input1 = (v8i16) __msa_ilvr_b(zero, src1);
input0 <<= 6;
input1 <<= 6;
STORE_2VECS_SH(dst, dst_stride, input0, input1);
} else if (4 == height) {
v16i8 src0, src1, src2, src3;
v8i16 input0, input1, input2, input3;
LOAD_4VECS_SB(src, src_stride, src0, src1, src2, src3);
input0 = (v8i16) __msa_ilvr_b(zero, src0);
input1 = (v8i16) __msa_ilvr_b(zero, src1);
input2 = (v8i16) __msa_ilvr_b(zero, src2);
input3 = (v8i16) __msa_ilvr_b(zero, src3);
input0 <<= 6;
input1 <<= 6;
input2 <<= 6;
input3 <<= 6;
STORE_4VECS_SH(dst, dst_stride, input0, input1, input2, input3);
} else if (6 == height) {
v16i8 src0, src1, src2, src3, src4, src5;
v8i16 input0, input1, input2, input3, input4, input5;
LOAD_6VECS_SB(src, src_stride, src0, src1, src2, src3, src4, src5);
input0 = (v8i16) __msa_ilvr_b(zero, src0);
input1 = (v8i16) __msa_ilvr_b(zero, src1);
input2 = (v8i16) __msa_ilvr_b(zero, src2);
input3 = (v8i16) __msa_ilvr_b(zero, src3);
input4 = (v8i16) __msa_ilvr_b(zero, src4);
input5 = (v8i16) __msa_ilvr_b(zero, src5);
input0 <<= 6;
input1 <<= 6;
input2 <<= 6;
input3 <<= 6;
input4 <<= 6;
input5 <<= 6;
STORE_6VECS_SH(dst, dst_stride,
input0, input1, input2, input3, input4, input5);
} else if (0 == height % 8) {
uint32_t loop_cnt;
v16i8 src0, src1, src2, src3, src4, src5, src6, src7;
v8i16 input0, input1, input2, input3;
v8i16 input4, input5, input6, input7;
for (loop_cnt = (height >> 3); loop_cnt--;) {
LOAD_8VECS_SB(src, src_stride,
src0, src1, src2, src3, src4, src5, src6, src7);
src += (8 * src_stride);
input0 = (v8i16) __msa_ilvr_b(zero, src0);
input1 = (v8i16) __msa_ilvr_b(zero, src1);
input2 = (v8i16) __msa_ilvr_b(zero, src2);
input3 = (v8i16) __msa_ilvr_b(zero, src3);
input4 = (v8i16) __msa_ilvr_b(zero, src4);
input5 = (v8i16) __msa_ilvr_b(zero, src5);
input6 = (v8i16) __msa_ilvr_b(zero, src6);
input7 = (v8i16) __msa_ilvr_b(zero, src7);
input0 <<= 6;
input1 <<= 6;
input2 <<= 6;
input3 <<= 6;
input4 <<= 6;
input5 <<= 6;
input6 <<= 6;
input7 <<= 6;
STORE_8VECS_SH(dst, dst_stride,
input0, input1, input2, input3,
input4, input5, input6, input7);
dst += (8 * dst_stride);
}
}
}
static void hevc_copy_12w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
int32_t height)
{
uint32_t loop_cnt;
uint64_t dst_val0, dst_val1, dst_val2, dst_val3;
v16i8 zero = { 0 };
v16i8 src0, src1, src2, src3, src4, src5, src6, src7;
v8i16 input0, input1;
v8i16 input0_r, input1_r, input2_r, input3_r;
for (loop_cnt = (height >> 3); loop_cnt--;) {
LOAD_8VECS_SB(src, src_stride,
src0, src1, src2, src3, src4, src5, src6, src7);
src += (8 * src_stride);
input0_r = (v8i16) __msa_ilvr_b(zero, src0);
input1_r = (v8i16) __msa_ilvr_b(zero, src1);
input2_r = (v8i16) __msa_ilvr_b(zero, src2);
input3_r = (v8i16) __msa_ilvr_b(zero, src3);
input0_r <<= 6;
input1_r <<= 6;
input2_r <<= 6;
input3_r <<= 6;
src0 = (v16i8) __msa_ilvl_w((v4i32) src1, (v4i32) src0);
src1 = (v16i8) __msa_ilvl_w((v4i32) src3, (v4i32) src2);
input0 = (v8i16) __msa_ilvr_b(zero, src0);
input1 = (v8i16) __msa_ilvr_b(zero, src1);
input0 <<= 6;
input1 <<= 6;
dst_val0 = __msa_copy_u_d((v2i64) input0, 0);
dst_val1 = __msa_copy_u_d((v2i64) input0, 1);
dst_val2 = __msa_copy_u_d((v2i64) input1, 0);
dst_val3 = __msa_copy_u_d((v2i64) input1, 1);
STORE_4VECS_SH(dst, dst_stride, input0_r, input1_r, input2_r, input3_r);
STORE_DWORD(dst + 8, dst_val0);
dst += dst_stride;
STORE_DWORD(dst + 8, dst_val1);
dst += dst_stride;
STORE_DWORD(dst + 8, dst_val2);
dst += dst_stride;
STORE_DWORD(dst + 8, dst_val3);
dst += dst_stride;
input0_r = (v8i16) __msa_ilvr_b(zero, src4);
input1_r = (v8i16) __msa_ilvr_b(zero, src5);
input2_r = (v8i16) __msa_ilvr_b(zero, src6);
input3_r = (v8i16) __msa_ilvr_b(zero, src7);
input0_r <<= 6;
input1_r <<= 6;
input2_r <<= 6;
input3_r <<= 6;
src0 = (v16i8) __msa_ilvl_w((v4i32) src5, (v4i32) src4);
src1 = (v16i8) __msa_ilvl_w((v4i32) src7, (v4i32) src6);
input0 = (v8i16) __msa_ilvr_b(zero, src0);
input1 = (v8i16) __msa_ilvr_b(zero, src1);
input0 <<= 6;
input1 <<= 6;
dst_val0 = __msa_copy_u_d((v2i64) input0, 0);
dst_val1 = __msa_copy_u_d((v2i64) input0, 1);
dst_val2 = __msa_copy_u_d((v2i64) input1, 0);
dst_val3 = __msa_copy_u_d((v2i64) input1, 1);
STORE_4VECS_SH(dst, dst_stride, input0_r, input1_r, input2_r, input3_r);
STORE_DWORD(dst + 8, dst_val0);
dst += dst_stride;
STORE_DWORD(dst + 8, dst_val1);
dst += dst_stride;
STORE_DWORD(dst + 8, dst_val2);
dst += dst_stride;
STORE_DWORD(dst + 8, dst_val3);
dst += dst_stride;
}
}
static void hevc_copy_16multx8mult_msa(uint8_t * __restrict src,
int32_t src_stride,
int16_t * __restrict dst,
int32_t dst_stride,
int32_t height,
int32_t width)
{
uint8_t *src_tmp;
int16_t *dst_tmp;
uint32_t loop_cnt, cnt;
v16i8 zero = { 0 };
v16i8 src0, src1, src2, src3, src4, src5, src6, src7;
v8i16 input0_r, input1_r, input2_r, input3_r;
v8i16 input0_l, input1_l, input2_l, input3_l;
for (cnt = (width >> 4); cnt--;) {
src_tmp = src;
dst_tmp = dst;
for (loop_cnt = (height >> 3); loop_cnt--;) {
LOAD_8VECS_SB(src_tmp, src_stride,
src0, src1, src2, src3, src4, src5, src6, src7);
src_tmp += (8 * src_stride);
input0_r = (v8i16) __msa_ilvr_b(zero, src0);
input0_l = (v8i16) __msa_ilvl_b(zero, src0);
input1_r = (v8i16) __msa_ilvr_b(zero, src1);
input1_l = (v8i16) __msa_ilvl_b(zero, src1);
input2_r = (v8i16) __msa_ilvr_b(zero, src2);
input2_l = (v8i16) __msa_ilvl_b(zero, src2);
input3_r = (v8i16) __msa_ilvr_b(zero, src3);
input3_l = (v8i16) __msa_ilvl_b(zero, src3);
input0_r <<= 6;
input0_l <<= 6;
input1_r <<= 6;
input1_l <<= 6;
input2_r <<= 6;
input2_l <<= 6;
input3_r <<= 6;
input3_l <<= 6;
STORE_4VECS_SH(dst_tmp, dst_stride,
input0_r, input1_r, input2_r, input3_r);
STORE_4VECS_SH((dst_tmp + 8), dst_stride,
input0_l, input1_l, input2_l, input3_l);
dst_tmp += (4 * dst_stride);
input0_r = (v8i16) __msa_ilvr_b(zero, src4);
input0_l = (v8i16) __msa_ilvl_b(zero, src4);
input1_r = (v8i16) __msa_ilvr_b(zero, src5);
input1_l = (v8i16) __msa_ilvl_b(zero, src5);
input2_r = (v8i16) __msa_ilvr_b(zero, src6);
input2_l = (v8i16) __msa_ilvl_b(zero, src6);
input3_r = (v8i16) __msa_ilvr_b(zero, src7);
input3_l = (v8i16) __msa_ilvl_b(zero, src7);
input0_r <<= 6;
input0_l <<= 6;
input1_r <<= 6;
input1_l <<= 6;
input2_r <<= 6;
input2_l <<= 6;
input3_r <<= 6;
input3_l <<= 6;
STORE_4VECS_SH(dst_tmp, dst_stride,
input0_r, input1_r, input2_r, input3_r);
STORE_4VECS_SH((dst_tmp + 8), dst_stride,
input0_l, input1_l, input2_l, input3_l);
dst_tmp += (4 * dst_stride);
}
src += 16;
dst += 16;
}
}
static void hevc_copy_16w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
int32_t height)
{
v16i8 zero = { 0 };
if (4 == height) {
v16i8 src0, src1, src2, src3;
v8i16 input0_r, input1_r, input2_r, input3_r;
v8i16 input0_l, input1_l, input2_l, input3_l;
LOAD_4VECS_SB(src, src_stride, src0, src1, src2, src3);
input0_r = (v8i16) __msa_ilvr_b(zero, src0);
input0_l = (v8i16) __msa_ilvl_b(zero, src0);
input1_r = (v8i16) __msa_ilvr_b(zero, src1);
input1_l = (v8i16) __msa_ilvl_b(zero, src1);
input2_r = (v8i16) __msa_ilvr_b(zero, src2);
input2_l = (v8i16) __msa_ilvl_b(zero, src2);
input3_r = (v8i16) __msa_ilvr_b(zero, src3);
input3_l = (v8i16) __msa_ilvl_b(zero, src3);
input0_r <<= 6;
input0_l <<= 6;
input1_r <<= 6;
input1_l <<= 6;
input2_r <<= 6;
input2_l <<= 6;
input3_r <<= 6;
input3_l <<= 6;
STORE_4VECS_SH(dst, dst_stride, input0_r, input1_r, input2_r, input3_r);
STORE_4VECS_SH((dst + 8), dst_stride,
input0_l, input1_l, input2_l, input3_l);
} else if (12 == height) {
v16i8 src0, src1, src2, src3, src4, src5, src6, src7;
v16i8 src8, src9, src10, src11;
v8i16 input0_r, input1_r, input2_r, input3_r;
v8i16 input0_l, input1_l, input2_l, input3_l;
LOAD_8VECS_SB(src, src_stride,
src0, src1, src2, src3, src4, src5, src6, src7);
src += (8 * src_stride);
LOAD_4VECS_SB(src, src_stride, src8, src9, src10, src11);
input0_r = (v8i16) __msa_ilvr_b(zero, src0);
input0_l = (v8i16) __msa_ilvl_b(zero, src0);
input1_r = (v8i16) __msa_ilvr_b(zero, src1);
input1_l = (v8i16) __msa_ilvl_b(zero, src1);
input2_r = (v8i16) __msa_ilvr_b(zero, src2);
input2_l = (v8i16) __msa_ilvl_b(zero, src2);
input3_r = (v8i16) __msa_ilvr_b(zero, src3);
input3_l = (v8i16) __msa_ilvl_b(zero, src3);
input0_r <<= 6;
input0_l <<= 6;
input1_r <<= 6;
input1_l <<= 6;
input2_r <<= 6;
input2_l <<= 6;
input3_r <<= 6;
input3_l <<= 6;
STORE_4VECS_SH(dst, dst_stride, input0_r, input1_r, input2_r, input3_r);
STORE_4VECS_SH((dst + 8), dst_stride,
input0_l, input1_l, input2_l, input3_l);
dst += (4 * dst_stride);
input0_r = (v8i16) __msa_ilvr_b(zero, src4);
input0_l = (v8i16) __msa_ilvl_b(zero, src4);
input1_r = (v8i16) __msa_ilvr_b(zero, src5);
input1_l = (v8i16) __msa_ilvl_b(zero, src5);
input2_r = (v8i16) __msa_ilvr_b(zero, src6);
input2_l = (v8i16) __msa_ilvl_b(zero, src6);
input3_r = (v8i16) __msa_ilvr_b(zero, src7);
input3_l = (v8i16) __msa_ilvl_b(zero, src7);
input0_r <<= 6;
input0_l <<= 6;
input1_r <<= 6;
input1_l <<= 6;
input2_r <<= 6;
input2_l <<= 6;
input3_r <<= 6;
input3_l <<= 6;
STORE_4VECS_SH(dst, dst_stride, input0_r, input1_r, input2_r, input3_r);
STORE_4VECS_SH((dst + 8), dst_stride,
input0_l, input1_l, input2_l, input3_l);
dst += (4 * dst_stride);
input0_r = (v8i16) __msa_ilvr_b(zero, src8);
input0_l = (v8i16) __msa_ilvl_b(zero, src8);
input1_r = (v8i16) __msa_ilvr_b(zero, src9);
input1_l = (v8i16) __msa_ilvl_b(zero, src9);
input2_r = (v8i16) __msa_ilvr_b(zero, src10);
input2_l = (v8i16) __msa_ilvl_b(zero, src10);
input3_r = (v8i16) __msa_ilvr_b(zero, src11);
input3_l = (v8i16) __msa_ilvl_b(zero, src11);
input0_r <<= 6;
input0_l <<= 6;
input1_r <<= 6;
input1_l <<= 6;
input2_r <<= 6;
input2_l <<= 6;
input3_r <<= 6;
input3_l <<= 6;
STORE_4VECS_SH(dst, dst_stride, input0_r, input1_r, input2_r, input3_r);
STORE_4VECS_SH((dst + 8), dst_stride,
input0_l, input1_l, input2_l, input3_l);
} else if (0 == (height % 8)) {
hevc_copy_16multx8mult_msa(src, src_stride, dst, dst_stride,
height, 16);
}
}
static void hevc_copy_24w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
int32_t height)
{
hevc_copy_16multx8mult_msa(src, src_stride, dst, dst_stride, height, 16);
hevc_copy_8w_msa(src + 16, src_stride, dst + 16, dst_stride, height);
}
static void hevc_copy_32w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
int32_t height)
{
hevc_copy_16multx8mult_msa(src, src_stride, dst, dst_stride, height, 32);
}
static void hevc_copy_48w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
int32_t height)
{
hevc_copy_16multx8mult_msa(src, src_stride, dst, dst_stride, height, 48);
}
static void hevc_copy_64w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
int32_t height)
{
hevc_copy_16multx8mult_msa(src, src_stride, dst, dst_stride, height, 64);
}
static void hevc_hz_8t_4w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
const int8_t * __restrict filter, int32_t height)
{
uint32_t loop_cnt;
uint64_t out0, out1, out2, out3, out4, out5, out6, out7;
v16i8 src0, src1, src2, src3, src4, src5, src6, src7;
v8i16 filt0, filt1, filt2, filt3;
v16i8 mask1, mask2, mask3;
v8u16 const_vec;
v16i8 vec0, vec1, vec2, vec3;
v8i16 dst0, dst1, dst2, dst3;
v8i16 filter_vec;
v16i8 mask0 = { 0, 1, 1, 2, 2, 3, 3, 4, 16, 17, 17, 18, 18, 19, 19, 20 };
src -= 3;
const_vec = (v8u16) __msa_ldi_h(128);
const_vec <<= 6;
filter_vec = LOAD_SH(filter);
filt0 = __msa_splati_h(filter_vec, 0);
filt1 = __msa_splati_h(filter_vec, 1);
filt2 = __msa_splati_h(filter_vec, 2);
filt3 = __msa_splati_h(filter_vec, 3);
mask1 = mask0 + 2;
mask2 = mask0 + 4;
mask3 = mask0 + 6;
for (loop_cnt = (height >> 3); loop_cnt--;) {
LOAD_8VECS_SB(src, src_stride,
src0, src1, src2, src3, src4, src5, src6, src7);
src += (8 * src_stride);
XORI_B_8VECS_SB(src0, src1, src2, src3, src4, src5, src6, src7,
src0, src1, src2, src3, src4, src5, src6, src7, 128);
vec0 = __msa_vshf_b(mask0, src1, src0);
vec1 = __msa_vshf_b(mask1, src1, src0);
vec2 = __msa_vshf_b(mask2, src1, src0);
vec3 = __msa_vshf_b(mask3, src1, src0);
dst0 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask0, src3, src2);
vec1 = __msa_vshf_b(mask1, src3, src2);
vec2 = __msa_vshf_b(mask2, src3, src2);
vec3 = __msa_vshf_b(mask3, src3, src2);
dst1 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask0, src5, src4);
vec1 = __msa_vshf_b(mask1, src5, src4);
vec2 = __msa_vshf_b(mask2, src5, src4);
vec3 = __msa_vshf_b(mask3, src5, src4);
dst2 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask0, src7, src6);
vec1 = __msa_vshf_b(mask1, src7, src6);
vec2 = __msa_vshf_b(mask2, src7, src6);
vec3 = __msa_vshf_b(mask3, src7, src6);
dst3 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
out0 = __msa_copy_u_d((v2i64) dst0, 0);
out1 = __msa_copy_u_d((v2i64) dst0, 1);
out2 = __msa_copy_u_d((v2i64) dst1, 0);
out3 = __msa_copy_u_d((v2i64) dst1, 1);
out4 = __msa_copy_u_d((v2i64) dst2, 0);
out5 = __msa_copy_u_d((v2i64) dst2, 1);
out6 = __msa_copy_u_d((v2i64) dst3, 0);
out7 = __msa_copy_u_d((v2i64) dst3, 1);
STORE_DWORD(dst, out0);
dst += dst_stride;
STORE_DWORD(dst, out1);
dst += dst_stride;
STORE_DWORD(dst, out2);
dst += dst_stride;
STORE_DWORD(dst, out3);
dst += dst_stride;
STORE_DWORD(dst, out4);
dst += dst_stride;
STORE_DWORD(dst, out5);
dst += dst_stride;
STORE_DWORD(dst, out6);
dst += dst_stride;
STORE_DWORD(dst, out7);
dst += dst_stride;
}
}
static void hevc_hz_8t_8w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
const int8_t * __restrict filter, int32_t height)
{
uint32_t loop_cnt;
v16i8 src0, src1, src2, src3;
v8i16 filt0, filt1, filt2, filt3;
v16i8 mask1, mask2, mask3;
v8u16 const_vec;
v16i8 vec0, vec1, vec2, vec3;
v8i16 dst0, dst1, dst2, dst3;
v8i16 filter_vec;
v16i8 mask0 = { 0, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 7, 8 };
src -= 3;
const_vec = (v8u16) __msa_ldi_h(128);
const_vec <<= 6;
filter_vec = LOAD_SH(filter);
filt0 = __msa_splati_h(filter_vec, 0);
filt1 = __msa_splati_h(filter_vec, 1);
filt2 = __msa_splati_h(filter_vec, 2);
filt3 = __msa_splati_h(filter_vec, 3);
mask1 = mask0 + 2;
mask2 = mask0 + 4;
mask3 = mask0 + 6;
for (loop_cnt = (height >> 2); loop_cnt--;) {
LOAD_4VECS_SB(src, src_stride, src0, src1, src2, src3);
src += (4 * src_stride);
XORI_B_4VECS_SB(src0, src1, src2, src3, src0, src1, src2, src3, 128);
vec0 = __msa_vshf_b(mask0, src0, src0);
vec1 = __msa_vshf_b(mask1, src0, src0);
vec2 = __msa_vshf_b(mask2, src0, src0);
vec3 = __msa_vshf_b(mask3, src0, src0);
dst0 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask0, src1, src1);
vec1 = __msa_vshf_b(mask1, src1, src1);
vec2 = __msa_vshf_b(mask2, src1, src1);
vec3 = __msa_vshf_b(mask3, src1, src1);
dst1 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask0, src2, src2);
vec1 = __msa_vshf_b(mask1, src2, src2);
vec2 = __msa_vshf_b(mask2, src2, src2);
vec3 = __msa_vshf_b(mask3, src2, src2);
dst2 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask0, src3, src3);
vec1 = __msa_vshf_b(mask1, src3, src3);
vec2 = __msa_vshf_b(mask2, src3, src3);
vec3 = __msa_vshf_b(mask3, src3, src3);
dst3 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
STORE_SH(dst0, dst);
dst += dst_stride;
STORE_SH(dst1, dst);
dst += dst_stride;
STORE_SH(dst2, dst);
dst += dst_stride;
STORE_SH(dst3, dst);
dst += dst_stride;
}
}
static void hevc_hz_8t_12w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
const int8_t * __restrict filter, int32_t height)
{
hevc_hz_8t_8w_msa(src, src_stride, dst, dst_stride, filter, height);
hevc_hz_8t_4w_msa(src + 8, src_stride, dst + 8, dst_stride, filter, height);
}
static void hevc_hz_8t_16w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
const int8_t * __restrict filter, int32_t height)
{
uint32_t loop_cnt;
v16i8 src0, src1, src2, src3, src4, src5, src6, src7;
v8i16 filt0, filt1, filt2, filt3;
v16i8 mask1, mask2, mask3;
v8u16 const_vec;
v16i8 vec0, vec1, vec2, vec3;
v8i16 dst0, dst1, dst2, dst3, dst4, dst5, dst6, dst7;
v8i16 filter_vec;
v16i8 mask0 = { 0, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 7, 8 };
src -= 3;
const_vec = (v8u16) __msa_ldi_h(128);
const_vec <<= 6;
filter_vec = LOAD_SH(filter);
filt0 = __msa_splati_h(filter_vec, 0);
filt1 = __msa_splati_h(filter_vec, 1);
filt2 = __msa_splati_h(filter_vec, 2);
filt3 = __msa_splati_h(filter_vec, 3);
mask1 = mask0 + 2;
mask2 = mask0 + 4;
mask3 = mask0 + 6;
for (loop_cnt = (height >> 2); loop_cnt--;) {
LOAD_4VECS_SB(src, src_stride, src0, src2, src4, src6);
LOAD_4VECS_SB(src + 8, src_stride, src1, src3, src5, src7);
src += (4 * src_stride);
XORI_B_8VECS_SB(src0, src1, src2, src3, src4, src5, src6, src7,
src0, src1, src2, src3, src4, src5, src6, src7, 128);
vec0 = __msa_vshf_b(mask0, src0, src0);
vec1 = __msa_vshf_b(mask1, src0, src0);
vec2 = __msa_vshf_b(mask2, src0, src0);
vec3 = __msa_vshf_b(mask3, src0, src0);
dst0 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask0, src1, src1);
vec1 = __msa_vshf_b(mask1, src1, src1);
vec2 = __msa_vshf_b(mask2, src1, src1);
vec3 = __msa_vshf_b(mask3, src1, src1);
dst1 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask0, src2, src2);
vec1 = __msa_vshf_b(mask1, src2, src2);
vec2 = __msa_vshf_b(mask2, src2, src2);
vec3 = __msa_vshf_b(mask3, src2, src2);
dst2 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask0, src3, src3);
vec1 = __msa_vshf_b(mask1, src3, src3);
vec2 = __msa_vshf_b(mask2, src3, src3);
vec3 = __msa_vshf_b(mask3, src3, src3);
dst3 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask0, src4, src4);
vec1 = __msa_vshf_b(mask1, src4, src4);
vec2 = __msa_vshf_b(mask2, src4, src4);
vec3 = __msa_vshf_b(mask3, src4, src4);
dst4 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask0, src5, src5);
vec1 = __msa_vshf_b(mask1, src5, src5);
vec2 = __msa_vshf_b(mask2, src5, src5);
vec3 = __msa_vshf_b(mask3, src5, src5);
dst5 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask0, src6, src6);
vec1 = __msa_vshf_b(mask1, src6, src6);
vec2 = __msa_vshf_b(mask2, src6, src6);
vec3 = __msa_vshf_b(mask3, src6, src6);
dst6 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask0, src7, src7);
vec1 = __msa_vshf_b(mask1, src7, src7);
vec2 = __msa_vshf_b(mask2, src7, src7);
vec3 = __msa_vshf_b(mask3, src7, src7);
dst7 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
STORE_SH(dst0, dst);
STORE_SH(dst1, dst + 8);
dst += dst_stride;
STORE_SH(dst2, dst);
STORE_SH(dst3, dst + 8);
dst += dst_stride;
STORE_SH(dst4, dst);
STORE_SH(dst5, dst + 8);
dst += dst_stride;
STORE_SH(dst6, dst);
STORE_SH(dst7, dst + 8);
dst += dst_stride;
}
}
static void hevc_hz_8t_24w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
const int8_t * __restrict filter, int32_t height)
{
uint32_t loop_cnt;
v16i8 src0, src1, src2, src3;
v8i16 filt0, filt1, filt2, filt3;
v16i8 mask1, mask2, mask3, mask4, mask5, mask6, mask7;
v16i8 vec0, vec1, vec2, vec3;
v8i16 dst0, dst1, dst2, dst3, dst4, dst5;
v8i16 filter_vec;
v8u16 const_vec;
v16i8 mask0 = { 0, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 7, 8 };
src -= 3;
filter_vec = LOAD_SH(filter);
filt0 = __msa_splati_h(filter_vec, 0);
filt1 = __msa_splati_h(filter_vec, 1);
filt2 = __msa_splati_h(filter_vec, 2);
filt3 = __msa_splati_h(filter_vec, 3);
mask1 = mask0 + 2;
mask2 = mask0 + 4;
mask3 = mask0 + 6;
mask4 = mask0 + 8;
mask5 = mask0 + 10;
mask6 = mask0 + 12;
mask7 = mask0 + 14;
const_vec = (v8u16) __msa_ldi_h(128);
const_vec <<= 6;
for (loop_cnt = (height >> 1); loop_cnt--;) {
src0 = LOAD_SB(src);
src1 = LOAD_SB(src + 16);
src += src_stride;
src2 = LOAD_SB(src);
src3 = LOAD_SB(src + 16);
src += src_stride;
XORI_B_4VECS_SB(src0, src1, src2, src3, src0, src1, src2, src3, 128);
vec0 = __msa_vshf_b(mask0, src0, src0);
vec1 = __msa_vshf_b(mask1, src0, src0);
vec2 = __msa_vshf_b(mask2, src0, src0);
vec3 = __msa_vshf_b(mask3, src0, src0);
dst0 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask4, src1, src0);
vec1 = __msa_vshf_b(mask5, src1, src0);
vec2 = __msa_vshf_b(mask6, src1, src0);
vec3 = __msa_vshf_b(mask7, src1, src0);
dst1 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask0, src1, src1);
vec1 = __msa_vshf_b(mask1, src1, src1);
vec2 = __msa_vshf_b(mask2, src1, src1);
vec3 = __msa_vshf_b(mask3, src1, src1);
dst2 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask0, src2, src2);
vec1 = __msa_vshf_b(mask1, src2, src2);
vec2 = __msa_vshf_b(mask2, src2, src2);
vec3 = __msa_vshf_b(mask3, src2, src2);
dst3 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask4, src3, src2);
vec1 = __msa_vshf_b(mask5, src3, src2);
vec2 = __msa_vshf_b(mask6, src3, src2);
vec3 = __msa_vshf_b(mask7, src3, src2);
dst4 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask0, src3, src3);
vec1 = __msa_vshf_b(mask1, src3, src3);
vec2 = __msa_vshf_b(mask2, src3, src3);
vec3 = __msa_vshf_b(mask3, src3, src3);
dst5 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
STORE_SH(dst0, dst);
STORE_SH(dst1, dst + 8);
STORE_SH(dst2, dst + 16);
dst += dst_stride;
STORE_SH(dst3, dst);
STORE_SH(dst4, dst + 8);
STORE_SH(dst5, dst + 16);
dst += dst_stride;
}
}
static void hevc_hz_8t_32w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
const int8_t * __restrict filter, int32_t height)
{
uint32_t loop_cnt;
v16i8 src0, src1, src2;
v8i16 filt0, filt1, filt2, filt3;
v16i8 mask1, mask2, mask3, mask4, mask5, mask6, mask7;
v16i8 vec0, vec1, vec2, vec3;
v8i16 dst0, dst1, dst2, dst3;
v8i16 filter_vec;
v8u16 const_vec;
v16i8 mask0 = { 0, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 7, 8 };
src -= 3;
filter_vec = LOAD_SH(filter);
filt0 = __msa_splati_h(filter_vec, 0);
filt1 = __msa_splati_h(filter_vec, 1);
filt2 = __msa_splati_h(filter_vec, 2);
filt3 = __msa_splati_h(filter_vec, 3);
mask1 = mask0 + 2;
mask2 = mask0 + 4;
mask3 = mask0 + 6;
mask4 = mask0 + 8;
mask5 = mask0 + 10;
mask6 = mask0 + 12;
mask7 = mask0 + 14;
const_vec = (v8u16) __msa_ldi_h(128);
const_vec <<= 6;
for (loop_cnt = height; loop_cnt--;) {
src0 = LOAD_SB(src);
src1 = LOAD_SB(src + 16);
src2 = LOAD_SB(src + 24);
src += src_stride;
XORI_B_3VECS_SB(src0, src1, src2, src0, src1, src2, 128);
vec0 = __msa_vshf_b(mask0, src0, src0);
vec1 = __msa_vshf_b(mask1, src0, src0);
vec2 = __msa_vshf_b(mask2, src0, src0);
vec3 = __msa_vshf_b(mask3, src0, src0);
dst0 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask4, src1, src0);
vec1 = __msa_vshf_b(mask5, src1, src0);
vec2 = __msa_vshf_b(mask6, src1, src0);
vec3 = __msa_vshf_b(mask7, src1, src0);
dst1 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask0, src1, src1);
vec1 = __msa_vshf_b(mask1, src1, src1);
vec2 = __msa_vshf_b(mask2, src1, src1);
vec3 = __msa_vshf_b(mask3, src1, src1);
dst2 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask0, src2, src2);
vec1 = __msa_vshf_b(mask1, src2, src2);
vec2 = __msa_vshf_b(mask2, src2, src2);
vec3 = __msa_vshf_b(mask3, src2, src2);
dst3 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
STORE_SH(dst0, dst);
STORE_SH(dst1, dst + 8);
STORE_SH(dst2, dst + 16);
STORE_SH(dst3, dst + 24);
dst += dst_stride;
}
}
static void hevc_hz_8t_48w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
const int8_t * __restrict filter, int32_t height)
{
uint32_t loop_cnt;
v16i8 src0, src1, src2, src3;
v8i16 filt0, filt1, filt2, filt3;
v16i8 mask1, mask2, mask3, mask4, mask5, mask6, mask7;
v16i8 vec0, vec1, vec2, vec3;
v8i16 dst0, dst1, dst2, dst3, dst4, dst5;
v8i16 filter_vec;
v8u16 const_vec;
v16i8 mask0 = { 0, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 7, 8 };
src -= 3;
filter_vec = LOAD_SH(filter);
filt0 = __msa_splati_h(filter_vec, 0);
filt1 = __msa_splati_h(filter_vec, 1);
filt2 = __msa_splati_h(filter_vec, 2);
filt3 = __msa_splati_h(filter_vec, 3);
mask1 = mask0 + 2;
mask2 = mask0 + 4;
mask3 = mask0 + 6;
mask4 = mask0 + 8;
mask5 = mask0 + 10;
mask6 = mask0 + 12;
mask7 = mask0 + 14;
const_vec = (v8u16) __msa_ldi_h(128);
const_vec <<= 6;
for (loop_cnt = height; loop_cnt--;) {
src0 = LOAD_SB(src);
src1 = LOAD_SB(src + 16);
src2 = LOAD_SB(src + 32);
src3 = LOAD_SB(src + 40);
src += src_stride;
XORI_B_4VECS_SB(src0, src1, src2, src3, src0, src1, src2, src3, 128);
vec0 = __msa_vshf_b(mask0, src0, src0);
vec1 = __msa_vshf_b(mask1, src0, src0);
vec2 = __msa_vshf_b(mask2, src0, src0);
vec3 = __msa_vshf_b(mask3, src0, src0);
dst0 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask4, src1, src0);
vec1 = __msa_vshf_b(mask5, src1, src0);
vec2 = __msa_vshf_b(mask6, src1, src0);
vec3 = __msa_vshf_b(mask7, src1, src0);
dst1 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask0, src1, src1);
vec1 = __msa_vshf_b(mask1, src1, src1);
vec2 = __msa_vshf_b(mask2, src1, src1);
vec3 = __msa_vshf_b(mask3, src1, src1);
dst2 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask4, src2, src1);
vec1 = __msa_vshf_b(mask5, src2, src1);
vec2 = __msa_vshf_b(mask6, src2, src1);
vec3 = __msa_vshf_b(mask7, src2, src1);
dst3 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask0, src2, src2);
vec1 = __msa_vshf_b(mask1, src2, src2);
vec2 = __msa_vshf_b(mask2, src2, src2);
vec3 = __msa_vshf_b(mask3, src2, src2);
dst4 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
vec0 = __msa_vshf_b(mask0, src3, src3);
vec1 = __msa_vshf_b(mask1, src3, src3);
vec2 = __msa_vshf_b(mask2, src3, src3);
vec3 = __msa_vshf_b(mask3, src3, src3);
dst5 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
STORE_SH(dst0, dst);
STORE_SH(dst1, dst + 8);
STORE_SH(dst2, dst + 16);
STORE_SH(dst3, dst + 24);
STORE_SH(dst4, dst + 32);
STORE_SH(dst5, dst + 40);
dst += dst_stride;
}
}
static void hevc_hz_8t_64w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
const int8_t * __restrict filter, int32_t height)
{
uint32_t loop_cnt;
v16i8 src0, src1, src2, src3, src4;
v8i16 filt0, filt1, filt2, filt3;
v16i8 mask1, mask2, mask3, mask4, mask5, mask6, mask7;
v16i8 vec0, vec1, vec2, vec3;
v8i16 dst0, dst1, dst2, dst3, dst4, dst5, dst6, dst7;
v8i16 filter_vec;
v8u16 const_vec;
v16i8 mask0 = { 0, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 7, 8 };
src -= 3;
filter_vec = LOAD_SH(filter);
filt0 = __msa_splati_h(filter_vec, 0);
filt1 = __msa_splati_h(filter_vec, 1);
filt2 = __msa_splati_h(filter_vec, 2);
filt3 = __msa_splati_h(filter_vec, 3);
mask1 = mask0 + 2;
mask2 = mask0 + 4;
mask3 = mask0 + 6;
mask4 = mask0 + 8;
mask5 = mask0 + 10;
mask6 = mask0 + 12;
mask7 = mask0 + 14;
const_vec = (v8u16) __msa_ldi_h(128);
const_vec <<= 6;
for (loop_cnt = height; loop_cnt--;) {
src0 = LOAD_SB(src);
src1 = LOAD_SB(src + 16);
src2 = LOAD_SB(src + 32);
src3 = LOAD_SB(src + 48);
src4 = LOAD_SB(src + 56);
src += src_stride;
XORI_B_5VECS_SB(src0, src1, src2, src3, src4,
src0, src1, src2, src3, src4, 128);
vec0 = __msa_vshf_b(mask0, src0, src0);
vec1 = __msa_vshf_b(mask1, src0, src0);
vec2 = __msa_vshf_b(mask2, src0, src0);
vec3 = __msa_vshf_b(mask3, src0, src0);
dst0 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
STORE_SH(dst0, dst);
vec0 = __msa_vshf_b(mask4, src1, src0);
vec1 = __msa_vshf_b(mask5, src1, src0);
vec2 = __msa_vshf_b(mask6, src1, src0);
vec3 = __msa_vshf_b(mask7, src1, src0);
dst1 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
STORE_SH(dst1, dst + 8);
vec0 = __msa_vshf_b(mask0, src1, src1);
vec1 = __msa_vshf_b(mask1, src1, src1);
vec2 = __msa_vshf_b(mask2, src1, src1);
vec3 = __msa_vshf_b(mask3, src1, src1);
dst2 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
STORE_SH(dst2, dst + 16);
vec0 = __msa_vshf_b(mask4, src2, src1);
vec1 = __msa_vshf_b(mask5, src2, src1);
vec2 = __msa_vshf_b(mask6, src2, src1);
vec3 = __msa_vshf_b(mask7, src2, src1);
dst3 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
STORE_SH(dst3, dst + 24);
vec0 = __msa_vshf_b(mask0, src2, src2);
vec1 = __msa_vshf_b(mask1, src2, src2);
vec2 = __msa_vshf_b(mask2, src2, src2);
vec3 = __msa_vshf_b(mask3, src2, src2);
dst4 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
STORE_SH(dst4, dst + 32);
vec0 = __msa_vshf_b(mask4, src3, src2);
vec1 = __msa_vshf_b(mask5, src3, src2);
vec2 = __msa_vshf_b(mask6, src3, src2);
vec3 = __msa_vshf_b(mask7, src3, src2);
dst5 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
STORE_SH(dst5, dst + 40);
vec0 = __msa_vshf_b(mask0, src3, src3);
vec1 = __msa_vshf_b(mask1, src3, src3);
vec2 = __msa_vshf_b(mask2, src3, src3);
vec3 = __msa_vshf_b(mask3, src3, src3);
dst6 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
STORE_SH(dst6, dst + 48);
vec0 = __msa_vshf_b(mask0, src4, src4);
vec1 = __msa_vshf_b(mask1, src4, src4);
vec2 = __msa_vshf_b(mask2, src4, src4);
vec3 = __msa_vshf_b(mask3, src4, src4);
dst7 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
STORE_SH(dst7, dst + 56);
dst += dst_stride;
}
}
static void hevc_vt_8t_4w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
const int8_t * __restrict filter, int32_t height)
{
int32_t loop_cnt;
uint64_t out0, out1, out2, out3;
v16i8 src0, src1, src2, src3, src4, src5, src6, src7, src8;
v16i8 src9, src10, src11, src12, src13, src14;
v16i8 src10_r, src32_r, src54_r, src76_r, src98_r;
v16i8 src21_r, src43_r, src65_r, src87_r, src109_r;
v16i8 src1110_r, src1211_r, src1312_r, src1413_r;
v16i8 src2110, src4332, src6554, src8776, src10998;
v16i8 src12111110, src14131312;
v8i16 dst10, dst32, dst54, dst76;
v8i16 filter_vec;
v8i16 filt0, filt1, filt2, filt3;
v8u16 const_vec;
src -= (3 * src_stride);
const_vec = (v8u16) __msa_ldi_h(128);
const_vec <<= 6;
filter_vec = LOAD_SH(filter);
filt0 = __msa_splati_h(filter_vec, 0);
filt1 = __msa_splati_h(filter_vec, 1);
filt2 = __msa_splati_h(filter_vec, 2);
filt3 = __msa_splati_h(filter_vec, 3);
LOAD_7VECS_SB(src, src_stride,
src0, src1, src2, src3, src4, src5, src6);
src += (7 * src_stride);
ILVR_B_6VECS_SB(src0, src2, src4, src1, src3, src5,
src1, src3, src5, src2, src4, src6,
src10_r, src32_r, src54_r, src21_r, src43_r, src65_r);
ILVR_D_3VECS_SB(src2110, src21_r, src10_r, src4332, src43_r, src32_r,
src6554, src65_r, src54_r);
XORI_B_3VECS_SB(src2110, src4332, src6554, src2110, src4332, src6554, 128);
for (loop_cnt = (height >> 3); loop_cnt--;) {
LOAD_8VECS_SB(src, src_stride,
src7, src8, src9, src10, src11, src12, src13, src14);
src += (8 * src_stride);
ILVR_B_8VECS_SB(src6, src7, src8, src9, src10, src11, src12, src13,
src7, src8, src9, src10, src11, src12, src13, src14,
src76_r, src87_r, src98_r, src109_r,
src1110_r, src1211_r, src1312_r, src1413_r);
ILVR_D_4VECS_SB(src8776, src87_r, src76_r, src10998, src109_r, src98_r,
src12111110, src1211_r, src1110_r,
src14131312, src1413_r, src1312_r);
XORI_B_4VECS_SB(src8776, src10998, src12111110, src14131312,
src8776, src10998, src12111110, src14131312, 128);
dst10 = HEVC_FILT_8TAP_DPADD_H(src2110, src4332, src6554, src8776,
filt0, filt1, filt2, filt3, const_vec);
dst32 = HEVC_FILT_8TAP_DPADD_H(src4332, src6554, src8776, src10998,
filt0, filt1, filt2, filt3, const_vec);
dst54 = HEVC_FILT_8TAP_DPADD_H(src6554, src8776, src10998, src12111110,
filt0, filt1, filt2, filt3, const_vec);
dst76 = HEVC_FILT_8TAP_DPADD_H(src8776, src10998,
src12111110, src14131312,
filt0, filt1, filt2, filt3, const_vec);
out0 = __msa_copy_u_d((v2i64) dst10, 0);
out1 = __msa_copy_u_d((v2i64) dst10, 1);
out2 = __msa_copy_u_d((v2i64) dst32, 0);
out3 = __msa_copy_u_d((v2i64) dst32, 1);
STORE_DWORD(dst, out0);
dst += dst_stride;
STORE_DWORD(dst, out1);
dst += dst_stride;
STORE_DWORD(dst, out2);
dst += dst_stride;
STORE_DWORD(dst, out3);
dst += dst_stride;
out0 = __msa_copy_u_d((v2i64) dst54, 0);
out1 = __msa_copy_u_d((v2i64) dst54, 1);
out2 = __msa_copy_u_d((v2i64) dst76, 0);
out3 = __msa_copy_u_d((v2i64) dst76, 1);
STORE_DWORD(dst, out0);
dst += dst_stride;
STORE_DWORD(dst, out1);
dst += dst_stride;
STORE_DWORD(dst, out2);
dst += dst_stride;
STORE_DWORD(dst, out3);
dst += dst_stride;
src2110 = src10998;
src4332 = src12111110;
src6554 = src14131312;
src6 = src14;
}
}
static void hevc_vt_8t_8w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
const int8_t * __restrict filter, int32_t height)
{
int32_t loop_cnt;
v16i8 src0, src1, src2, src3, src4, src5, src6, src7, src8, src9, src10;
v16i8 src10_r, src32_r, src54_r, src76_r, src98_r;
v16i8 src21_r, src43_r, src65_r, src87_r, src109_r;
v8i16 dst0_r, dst1_r, dst2_r, dst3_r;
v8i16 filter_vec;
v8i16 filt0, filt1, filt2, filt3;
v8u16 const_vec;
src -= (3 * src_stride);
const_vec = (v8u16) __msa_ldi_h(128);
const_vec <<= 6;
filter_vec = LOAD_SH(filter);
filt0 = __msa_splati_h(filter_vec, 0);
filt1 = __msa_splati_h(filter_vec, 1);
filt2 = __msa_splati_h(filter_vec, 2);
filt3 = __msa_splati_h(filter_vec, 3);
LOAD_7VECS_SB(src, src_stride,
src0, src1, src2, src3, src4, src5, src6);
src += (7 * src_stride);
XORI_B_7VECS_SB(src0, src1, src2, src3, src4, src5, src6,
src0, src1, src2, src3, src4, src5, src6, 128);
ILVR_B_6VECS_SB(src0, src2, src4, src1, src3, src5,
src1, src3, src5, src2, src4, src6,
src10_r, src32_r, src54_r, src21_r, src43_r, src65_r);
for (loop_cnt = (height >> 2); loop_cnt--;) {
LOAD_4VECS_SB(src, src_stride, src7, src8, src9, src10);
src += (4 * src_stride);
XORI_B_4VECS_SB(src7, src8, src9, src10, src7, src8, src9, src10, 128);
ILVR_B_4VECS_SB(src6, src7, src8, src9, src7, src8, src9, src10,
src76_r, src87_r, src98_r, src109_r);
dst0_r = HEVC_FILT_8TAP_DPADD_H(src10_r, src32_r, src54_r, src76_r,
filt0, filt1, filt2, filt3, const_vec);
dst1_r = HEVC_FILT_8TAP_DPADD_H(src21_r, src43_r, src65_r, src87_r,
filt0, filt1, filt2, filt3, const_vec);
dst2_r = HEVC_FILT_8TAP_DPADD_H(src32_r, src54_r, src76_r, src98_r,
filt0, filt1, filt2, filt3, const_vec);
dst3_r = HEVC_FILT_8TAP_DPADD_H(src43_r, src65_r, src87_r, src109_r,
filt0, filt1, filt2, filt3, const_vec);
STORE_SH(dst0_r, dst);
dst += dst_stride;
STORE_SH(dst1_r, dst);
dst += dst_stride;
STORE_SH(dst2_r, dst);
dst += dst_stride;
STORE_SH(dst3_r, dst);
dst += dst_stride;
src10_r = src54_r;
src32_r = src76_r;
src54_r = src98_r;
src21_r = src65_r;
src43_r = src87_r;
src65_r = src109_r;
src6 = src10;
}
}
static void hevc_vt_8t_12w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
const int8_t * __restrict filter, int32_t height)
{
int32_t loop_cnt;
uint64_t out0, out1, out2, out3;
v16i8 src0, src1, src2, src3, src4, src5, src6, src7, src8, src9, src10;
v16i8 src10_r, src32_r, src54_r, src76_r, src98_r;
v16i8 src21_r, src43_r, src65_r, src87_r, src109_r;
v8i16 dst0_r, dst1_r, dst2_r, dst3_r;
v16i8 src10_l, src32_l, src54_l, src76_l, src98_l;
v16i8 src21_l, src43_l, src65_l, src87_l, src109_l;
v16i8 src2110, src4332, src6554, src8776, src10998;
v8i16 dst0_l, dst1_l;
v8i16 filter_vec;
v8i16 filt0, filt1, filt2, filt3;
v8u16 const_vec;
src -= (3 * src_stride);
const_vec = (v8u16) __msa_ldi_h(128);
const_vec <<= 6;
filter_vec = LOAD_SH(filter);
filt0 = __msa_splati_h(filter_vec, 0);
filt1 = __msa_splati_h(filter_vec, 1);
filt2 = __msa_splati_h(filter_vec, 2);
filt3 = __msa_splati_h(filter_vec, 3);
LOAD_7VECS_SB(src, src_stride,
src0, src1, src2, src3, src4, src5, src6);
src += (7 * src_stride);
XORI_B_7VECS_SB(src0, src1, src2, src3, src4, src5, src6,
src0, src1, src2, src3, src4, src5, src6, 128);
ILVR_B_6VECS_SB(src0, src2, src4, src1, src3, src5,
src1, src3, src5, src2, src4, src6,
src10_r, src32_r, src54_r, src21_r, src43_r, src65_r);
ILVL_B_6VECS_SB(src0, src2, src4, src1, src3, src5,
src1, src3, src5, src2, src4, src6,
src10_l, src32_l, src54_l, src21_l, src43_l, src65_l);
ILVR_D_3VECS_SB(src2110, src21_l, src10_l, src4332, src43_l, src32_l,
src6554, src65_l, src54_l);
for (loop_cnt = (height >> 2); loop_cnt--;) {
LOAD_4VECS_SB(src, src_stride, src7, src8, src9, src10);
src += (4 * src_stride);
XORI_B_4VECS_SB(src7, src8, src9, src10, src7, src8, src9, src10, 128);
ILVR_B_4VECS_SB(src6, src7, src8, src9, src7, src8, src9, src10,
src76_r, src87_r, src98_r, src109_r);
ILVL_B_4VECS_SB(src6, src7, src8, src9, src7, src8, src9, src10,
src76_l, src87_l, src98_l, src109_l);
ILVR_D_2VECS_SB(src8776, src87_l, src76_l, src10998, src109_l, src98_l);
dst0_r = HEVC_FILT_8TAP_DPADD_H(src10_r, src32_r, src54_r, src76_r,
filt0, filt1, filt2, filt3, const_vec);
dst1_r = HEVC_FILT_8TAP_DPADD_H(src21_r, src43_r, src65_r, src87_r,
filt0, filt1, filt2, filt3, const_vec);
dst2_r = HEVC_FILT_8TAP_DPADD_H(src32_r, src54_r, src76_r, src98_r,
filt0, filt1, filt2, filt3, const_vec);
dst3_r = HEVC_FILT_8TAP_DPADD_H(src43_r, src65_r, src87_r, src109_r,
filt0, filt1, filt2, filt3, const_vec);
dst0_l = HEVC_FILT_8TAP_DPADD_H(src2110, src4332, src6554, src8776,
filt0, filt1, filt2, filt3, const_vec);
dst1_l = HEVC_FILT_8TAP_DPADD_H(src4332, src6554, src8776, src10998,
filt0, filt1, filt2, filt3, const_vec);
out0 = __msa_copy_u_d((v2i64) dst0_l, 0);
out1 = __msa_copy_u_d((v2i64) dst0_l, 1);
out2 = __msa_copy_u_d((v2i64) dst1_l, 0);
out3 = __msa_copy_u_d((v2i64) dst1_l, 1);
STORE_SH(dst0_r, dst);
STORE_DWORD(dst + 8, out0);
dst += dst_stride;
STORE_SH(dst1_r, dst);
STORE_DWORD(dst + 8, out1);
dst += dst_stride;
STORE_SH(dst2_r, dst);
STORE_DWORD(dst + 8, out2);
dst += dst_stride;
STORE_SH(dst3_r, dst);
STORE_DWORD(dst + 8, out3);
dst += dst_stride;
src10_r = src54_r;
src32_r = src76_r;
src54_r = src98_r;
src21_r = src65_r;
src43_r = src87_r;
src65_r = src109_r;
src2110 = src6554;
src4332 = src8776;
src6554 = src10998;
src6 = src10;
}
}
static void hevc_vt_8t_16multx4mult_msa(uint8_t * __restrict src,
int32_t src_stride,
int16_t * __restrict dst,
int32_t dst_stride,
const int8_t * __restrict filter,
int32_t height,
int32_t width)
{
uint8_t *src_tmp;
int16_t *dst_tmp;
int32_t loop_cnt, cnt;
v16i8 src0, src1, src2, src3, src4, src5, src6, src7, src8, src9, src10;
v16i8 src10_r, src32_r, src54_r, src76_r, src98_r;
v16i8 src21_r, src43_r, src65_r, src87_r, src109_r;
v8i16 dst0_r, dst1_r, dst2_r, dst3_r;
v16i8 src10_l, src32_l, src54_l, src76_l, src98_l;
v16i8 src21_l, src43_l, src65_l, src87_l, src109_l;
v8i16 dst0_l, dst1_l, dst2_l, dst3_l;
v8i16 filter_vec;
v8i16 filt0, filt1, filt2, filt3;
v8u16 const_vec;
src -= (3 * src_stride);
const_vec = (v8u16) __msa_ldi_h(128);
const_vec <<= 6;
filter_vec = LOAD_SH(filter);
filt0 = __msa_splati_h(filter_vec, 0);
filt1 = __msa_splati_h(filter_vec, 1);
filt2 = __msa_splati_h(filter_vec, 2);
filt3 = __msa_splati_h(filter_vec, 3);
for (cnt = width >> 4; cnt--;) {
src_tmp = src;
dst_tmp = dst;
LOAD_7VECS_SB(src_tmp, src_stride,
src0, src1, src2, src3, src4, src5, src6);
src_tmp += (7 * src_stride);
XORI_B_7VECS_SB(src0, src1, src2, src3, src4, src5, src6,
src0, src1, src2, src3, src4, src5, src6, 128);
ILVR_B_6VECS_SB(src0, src2, src4, src1, src3, src5,
src1, src3, src5, src2, src4, src6,
src10_r, src32_r, src54_r, src21_r, src43_r, src65_r);
ILVL_B_6VECS_SB(src0, src2, src4, src1, src3, src5,
src1, src3, src5, src2, src4, src6,
src10_l, src32_l, src54_l, src21_l, src43_l, src65_l);
for (loop_cnt = (height >> 2); loop_cnt--;) {
LOAD_4VECS_SB(src_tmp, src_stride, src7, src8, src9, src10);
src_tmp += (4 * src_stride);
XORI_B_4VECS_SB(src7, src8, src9, src10,
src7, src8, src9, src10, 128);
ILVR_B_4VECS_SB(src6, src7, src8, src9, src7, src8, src9, src10,
src76_r, src87_r, src98_r, src109_r);
ILVL_B_4VECS_SB(src6, src7, src8, src9, src7, src8, src9, src10,
src76_l, src87_l, src98_l, src109_l);
dst0_r = HEVC_FILT_8TAP_DPADD_H(src10_r, src32_r, src54_r, src76_r,
filt0, filt1, filt2, filt3,
const_vec);
dst1_r = HEVC_FILT_8TAP_DPADD_H(src21_r, src43_r, src65_r, src87_r,
filt0, filt1, filt2, filt3,
const_vec);
dst2_r = HEVC_FILT_8TAP_DPADD_H(src32_r, src54_r, src76_r, src98_r,
filt0, filt1, filt2, filt3,
const_vec);
dst3_r = HEVC_FILT_8TAP_DPADD_H(src43_r, src65_r, src87_r, src109_r,
filt0, filt1, filt2, filt3,
const_vec);
dst0_l = HEVC_FILT_8TAP_DPADD_H(src10_l, src32_l, src54_l, src76_l,
filt0, filt1, filt2, filt3,
const_vec);
dst1_l = HEVC_FILT_8TAP_DPADD_H(src21_l, src43_l, src65_l, src87_l,
filt0, filt1, filt2, filt3,
const_vec);
dst2_l = HEVC_FILT_8TAP_DPADD_H(src32_l, src54_l, src76_l, src98_l,
filt0, filt1, filt2, filt3,
const_vec);
dst3_l = HEVC_FILT_8TAP_DPADD_H(src43_l, src65_l, src87_l, src109_l,
filt0, filt1, filt2, filt3,
const_vec);
STORE_SH(dst0_r, dst_tmp);
STORE_SH(dst0_l, dst_tmp + 8);
dst_tmp += dst_stride;
STORE_SH(dst1_r, dst_tmp);
STORE_SH(dst1_l, dst_tmp + 8);
dst_tmp += dst_stride;
STORE_SH(dst2_r, dst_tmp);
STORE_SH(dst2_l, dst_tmp + 8);
dst_tmp += dst_stride;
STORE_SH(dst3_r, dst_tmp);
STORE_SH(dst3_l, dst_tmp + 8);
dst_tmp += dst_stride;
src10_r = src54_r;
src32_r = src76_r;
src54_r = src98_r;
src21_r = src65_r;
src43_r = src87_r;
src65_r = src109_r;
src10_l = src54_l;
src32_l = src76_l;
src54_l = src98_l;
src21_l = src65_l;
src43_l = src87_l;
src65_l = src109_l;
src6 = src10;
}
src += 16;
dst += 16;
}
}
static void hevc_vt_8t_16w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
const int8_t * __restrict filter, int32_t height)
{
hevc_vt_8t_16multx4mult_msa(src, src_stride, dst, dst_stride,
filter, height, 16);
}
static void hevc_vt_8t_24w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
const int8_t * __restrict filter, int32_t height)
{
hevc_vt_8t_16multx4mult_msa(src, src_stride, dst, dst_stride,
filter, height, 16);
hevc_vt_8t_8w_msa(src + 16, src_stride, dst + 16, dst_stride,
filter, height);
}
static void hevc_vt_8t_32w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
const int8_t * __restrict filter, int32_t height)
{
hevc_vt_8t_16multx4mult_msa(src, src_stride, dst, dst_stride,
filter, height, 32);
}
static void hevc_vt_8t_48w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
const int8_t * __restrict filter, int32_t height)
{
hevc_vt_8t_16multx4mult_msa(src, src_stride, dst, dst_stride,
filter, height, 48);
}
static void hevc_vt_8t_64w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
const int8_t * __restrict filter, int32_t height)
{
hevc_vt_8t_16multx4mult_msa(src, src_stride, dst, dst_stride,
filter, height, 64);
}
static void hevc_hv_8t_4w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
const int8_t * __restrict filter_x,
const int8_t * __restrict filter_y,
int32_t height)
{
uint32_t loop_cnt;
uint64_t out0, out1;
v16i8 src0, src1, src2, src3, src4, src5, src6, src7, src8;
v8i16 filt0, filt1, filt2, filt3, filter_vec;
v4i32 filt_h0, filt_h1, filt_h2, filt_h3;
v16i8 mask1, mask2, mask3;
v8u16 const_vec;
v16i8 vec0, vec1, vec2, vec3, vec4, vec5, vec6, vec7;
v16i8 vec8, vec9, vec10, vec11, vec12, vec13, vec14, vec15;
v8i16 dst30, dst41, dst52, dst63, dst66, dst87;
v4i32 dst0_r, dst1_r;
v8i16 dst10_r, dst32_r, dst54_r, dst76_r;
v8i16 dst21_r, dst43_r, dst65_r, dst87_r;
v16i8 tmp;
v16i8 mask0 = {
0, 1, 1, 2, 2, 3, 3, 4, 16, 17, 17, 18, 18, 19, 19, 20
};
v8u16 mask4 = { 0, 4, 1, 5, 2, 6, 3, 7 };
src -= ((3 * src_stride) + 3);
filter_vec = LOAD_SH(filter_x);
filt0 = __msa_splati_h(filter_vec, 0);
filt1 = __msa_splati_h(filter_vec, 1);
filt2 = __msa_splati_h(filter_vec, 2);
filt3 = __msa_splati_h(filter_vec, 3);
filter_vec = LOAD_SH(filter_y);
tmp = __msa_clti_s_b((v16i8) filter_vec, 0);
filter_vec = (v8i16) __msa_ilvr_b(tmp, (v16i8) filter_vec);
filt_h0 = __msa_splati_w((v4i32) filter_vec, 0);
filt_h1 = __msa_splati_w((v4i32) filter_vec, 1);
filt_h2 = __msa_splati_w((v4i32) filter_vec, 2);
filt_h3 = __msa_splati_w((v4i32) filter_vec, 3);
mask1 = mask0 + 2;
mask2 = mask0 + 4;
mask3 = mask0 + 6;
const_vec = (v8u16) __msa_ldi_h(128);
const_vec <<= 6;
LOAD_7VECS_SB(src, src_stride,
src0, src1, src2, src3, src4, src5, src6);
src += (7 * src_stride);
XORI_B_7VECS_SB(src0, src1, src2, src3, src4, src5, src6,
src0, src1, src2, src3, src4, src5, src6, 128);
/* Row 0 Row 1 Row 2 Row 3 */
vec0 = __msa_vshf_b(mask0, src3, src0);
vec1 = __msa_vshf_b(mask1, src3, src0);
vec2 = __msa_vshf_b(mask2, src3, src0);
vec3 = __msa_vshf_b(mask3, src3, src0);
vec4 = __msa_vshf_b(mask0, src4, src1);
vec5 = __msa_vshf_b(mask1, src4, src1);
vec6 = __msa_vshf_b(mask2, src4, src1);
vec7 = __msa_vshf_b(mask3, src4, src1);
vec8 = __msa_vshf_b(mask0, src5, src2);
vec9 = __msa_vshf_b(mask1, src5, src2);
vec10 = __msa_vshf_b(mask2, src5, src2);
vec11 = __msa_vshf_b(mask3, src5, src2);
vec12 = __msa_vshf_b(mask0, src6, src3);
vec13 = __msa_vshf_b(mask1, src6, src3);
vec14 = __msa_vshf_b(mask2, src6, src3);
vec15 = __msa_vshf_b(mask3, src6, src3);
dst30 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
dst41 = HEVC_FILT_8TAP_DPADD_H(vec4, vec5, vec6, vec7,
filt0, filt1, filt2, filt3, const_vec);
dst52 = HEVC_FILT_8TAP_DPADD_H(vec8, vec9, vec10, vec11,
filt0, filt1, filt2, filt3, const_vec);
dst63 = HEVC_FILT_8TAP_DPADD_H(vec12, vec13, vec14, vec15,
filt0, filt1, filt2, filt3, const_vec);
dst10_r = __msa_ilvr_h(dst41, dst30);
dst21_r = __msa_ilvr_h(dst52, dst41);
dst32_r = __msa_ilvr_h(dst63, dst52);
dst43_r = __msa_ilvl_h(dst41, dst30);
dst54_r = __msa_ilvl_h(dst52, dst41);
dst65_r = __msa_ilvl_h(dst63, dst52);
dst66 = (v8i16) __msa_splati_d((v2i64) dst63, 1);
for (loop_cnt = height >> 1; loop_cnt--;) {
LOAD_2VECS_SB(src, src_stride, src7, src8);
src += (2 * src_stride);
src7 = (v16i8) __msa_xori_b((v16u8) src7, 128);
src8 = (v16i8) __msa_xori_b((v16u8) src8, 128);
vec0 = __msa_vshf_b(mask0, src8, src7);
vec1 = __msa_vshf_b(mask1, src8, src7);
vec2 = __msa_vshf_b(mask2, src8, src7);
vec3 = __msa_vshf_b(mask3, src8, src7);
dst87 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
dst76_r = __msa_ilvr_h(dst87, dst66);
dst0_r = HEVC_FILT_8TAP_DPADD_W(dst10_r, dst32_r, dst54_r, dst76_r,
filt_h0, filt_h1, filt_h2, filt_h3);
dst87_r = __msa_vshf_h((v8i16) mask4, dst87, dst87);
dst1_r = HEVC_FILT_8TAP_DPADD_W(dst21_r, dst43_r, dst65_r, dst87_r,
filt_h0, filt_h1, filt_h2, filt_h3);
dst0_r >>= 6;
dst1_r >>= 6;
dst0_r = (v4i32) __msa_pckev_h((v8i16) dst1_r, (v8i16) dst0_r);
out0 = __msa_copy_u_d((v2i64) dst0_r, 0);
out1 = __msa_copy_u_d((v2i64) dst0_r, 1);
STORE_DWORD(dst, out0);
dst += dst_stride;
STORE_DWORD(dst, out1);
dst += dst_stride;
dst10_r = dst32_r;
dst32_r = dst54_r;
dst54_r = dst76_r;
dst21_r = dst43_r;
dst43_r = dst65_r;
dst65_r = dst87_r;
dst66 = (v8i16) __msa_splati_d((v2i64) dst87, 1);
}
}
static void hevc_hv_8t_8multx2mult_msa(uint8_t * __restrict src,
int32_t src_stride,
int16_t * __restrict dst,
int32_t dst_stride,
const int8_t * __restrict filter_x,
const int8_t * __restrict filter_y,
int32_t height, int32_t width)
{
uint32_t loop_cnt, cnt;
uint8_t *src_tmp;
int16_t *dst_tmp;
v16i8 src0, src1, src2, src3, src4, src5, src6, src7, src8;
v8i16 filt0, filt1, filt2, filt3, filter_vec;
v4i32 filt_h0, filt_h1, filt_h2, filt_h3;
v16i8 mask1, mask2, mask3;
v8u16 const_vec;
v16i8 vec0, vec1, vec2, vec3, vec4, vec5, vec6, vec7;
v16i8 vec8, vec9, vec10, vec11, vec12, vec13, vec14, vec15;
v8i16 dst0, dst1, dst2, dst3, dst4, dst5, dst6, dst7, dst8;
v4i32 dst0_r, dst0_l;
v8i16 dst10_r, dst32_r, dst54_r, dst76_r;
v8i16 dst10_l, dst32_l, dst54_l, dst76_l;
v8i16 dst21_r, dst43_r, dst65_r, dst87_r;
v8i16 dst21_l, dst43_l, dst65_l, dst87_l;
v16i8 tmp;
v16i8 mask0 = { 0, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 7, 8 };
src -= ((3 * src_stride) + 3);
filter_vec = LOAD_SH(filter_x);
filt0 = __msa_splati_h(filter_vec, 0);
filt1 = __msa_splati_h(filter_vec, 1);
filt2 = __msa_splati_h(filter_vec, 2);
filt3 = __msa_splati_h(filter_vec, 3);
filter_vec = LOAD_SH(filter_y);
tmp = __msa_clti_s_b((v16i8) filter_vec, 0);
filter_vec = (v8i16) __msa_ilvr_b(tmp, (v16i8) filter_vec);
filt_h0 = __msa_splati_w((v4i32) filter_vec, 0);
filt_h1 = __msa_splati_w((v4i32) filter_vec, 1);
filt_h2 = __msa_splati_w((v4i32) filter_vec, 2);
filt_h3 = __msa_splati_w((v4i32) filter_vec, 3);
mask1 = mask0 + 2;
mask2 = mask0 + 4;
mask3 = mask0 + 6;
const_vec = (v8u16) __msa_ldi_h(128);
const_vec <<= 6;
for (cnt = width >> 3; cnt--;) {
src_tmp = src;
dst_tmp = dst;
LOAD_7VECS_SB(src_tmp, src_stride,
src0, src1, src2, src3, src4, src5, src6);
src_tmp += (7 * src_stride);
XORI_B_7VECS_SB(src0, src1, src2, src3, src4, src5, src6,
src0, src1, src2, src3, src4, src5, src6, 128);
/* Row 0 Row 1 Row 2 Row 3 */
vec0 = __msa_vshf_b(mask0, src0, src0);
vec1 = __msa_vshf_b(mask1, src0, src0);
vec2 = __msa_vshf_b(mask2, src0, src0);
vec3 = __msa_vshf_b(mask3, src0, src0);
vec4 = __msa_vshf_b(mask0, src1, src1);
vec5 = __msa_vshf_b(mask1, src1, src1);
vec6 = __msa_vshf_b(mask2, src1, src1);
vec7 = __msa_vshf_b(mask3, src1, src1);
vec8 = __msa_vshf_b(mask0, src2, src2);
vec9 = __msa_vshf_b(mask1, src2, src2);
vec10 = __msa_vshf_b(mask2, src2, src2);
vec11 = __msa_vshf_b(mask3, src2, src2);
vec12 = __msa_vshf_b(mask0, src3, src3);
vec13 = __msa_vshf_b(mask1, src3, src3);
vec14 = __msa_vshf_b(mask2, src3, src3);
vec15 = __msa_vshf_b(mask3, src3, src3);
dst0 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
dst1 = HEVC_FILT_8TAP_DPADD_H(vec4, vec5, vec6, vec7,
filt0, filt1, filt2, filt3, const_vec);
dst2 = HEVC_FILT_8TAP_DPADD_H(vec8, vec9, vec10, vec11,
filt0, filt1, filt2, filt3, const_vec);
dst3 = HEVC_FILT_8TAP_DPADD_H(vec12, vec13, vec14, vec15,
filt0, filt1, filt2, filt3, const_vec);
/* Row 4 Row 5 Row 6 */
vec0 = __msa_vshf_b(mask0, src4, src4);
vec1 = __msa_vshf_b(mask1, src4, src4);
vec2 = __msa_vshf_b(mask2, src4, src4);
vec3 = __msa_vshf_b(mask3, src4, src4);
vec4 = __msa_vshf_b(mask0, src5, src5);
vec5 = __msa_vshf_b(mask1, src5, src5);
vec6 = __msa_vshf_b(mask2, src5, src5);
vec7 = __msa_vshf_b(mask3, src5, src5);
vec8 = __msa_vshf_b(mask0, src6, src6);
vec9 = __msa_vshf_b(mask1, src6, src6);
vec10 = __msa_vshf_b(mask2, src6, src6);
vec11 = __msa_vshf_b(mask3, src6, src6);
dst4 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3, const_vec);
dst5 = HEVC_FILT_8TAP_DPADD_H(vec4, vec5, vec6, vec7,
filt0, filt1, filt2, filt3, const_vec);
dst6 = HEVC_FILT_8TAP_DPADD_H(vec8, vec9, vec10, vec11,
filt0, filt1, filt2, filt3, const_vec);
ILVR_H_6VECS_SH(dst0, dst2, dst4, dst1, dst3, dst5,
dst1, dst3, dst5, dst2, dst4, dst6,
dst10_r, dst32_r, dst54_r, dst21_r, dst43_r, dst65_r);
ILVL_H_6VECS_SH(dst0, dst2, dst4, dst1, dst3, dst5,
dst1, dst3, dst5, dst2, dst4, dst6,
dst10_l, dst32_l, dst54_l, dst21_l, dst43_l, dst65_l);
for (loop_cnt = height >> 1; loop_cnt--;) {
src7 = LOAD_SB(src_tmp);
src_tmp += src_stride;
src7 = (v16i8) __msa_xori_b((v16u8) src7, 128);
vec0 = __msa_vshf_b(mask0, src7, src7);
vec1 = __msa_vshf_b(mask1, src7, src7);
vec2 = __msa_vshf_b(mask2, src7, src7);
vec3 = __msa_vshf_b(mask3, src7, src7);
dst7 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3,
const_vec);
dst76_r = __msa_ilvr_h(dst7, dst6);
dst76_l = __msa_ilvl_h(dst7, dst6);
dst0_r = HEVC_FILT_8TAP_DPADD_W(dst10_r, dst32_r, dst54_r, dst76_r,
filt_h0, filt_h1, filt_h2, filt_h3);
dst0_l = HEVC_FILT_8TAP_DPADD_W(dst10_l, dst32_l, dst54_l, dst76_l,
filt_h0, filt_h1, filt_h2, filt_h3);
dst0_r >>= 6;
dst0_l >>= 6;
dst0_r = (v4i32) __msa_pckev_h((v8i16) dst0_l, (v8i16) dst0_r);
STORE_SW(dst0_r, dst_tmp);
dst_tmp += dst_stride;
/* Next row */
src8 = LOAD_SB(src_tmp);
src_tmp += src_stride;
src8 = (v16i8) __msa_xori_b((v16u8) src8, 128);
vec0 = __msa_vshf_b(mask0, src8, src8);
vec1 = __msa_vshf_b(mask1, src8, src8);
vec2 = __msa_vshf_b(mask2, src8, src8);
vec3 = __msa_vshf_b(mask3, src8, src8);
dst8 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
filt0, filt1, filt2, filt3,
const_vec);
dst87_r = __msa_ilvr_h(dst8, dst7);
dst87_l = __msa_ilvl_h(dst8, dst7);
dst6 = dst8;
dst0_r = HEVC_FILT_8TAP_DPADD_W(dst21_r, dst43_r, dst65_r, dst87_r,
filt_h0, filt_h1, filt_h2, filt_h3);
dst0_l = HEVC_FILT_8TAP_DPADD_W(dst21_l, dst43_l, dst65_l, dst87_l,
filt_h0, filt_h1, filt_h2, filt_h3);
dst0_r >>= 6;
dst0_l >>= 6;
dst0_r = (v4i32) __msa_pckev_h((v8i16) dst0_l, (v8i16) dst0_r);
STORE_SW(dst0_r, dst_tmp);
dst_tmp += dst_stride;
dst10_r = dst32_r;
dst32_r = dst54_r;
dst54_r = dst76_r;
dst10_l = dst32_l;
dst32_l = dst54_l;
dst54_l = dst76_l;
dst21_r = dst43_r;
dst43_r = dst65_r;
dst65_r = dst87_r;
dst21_l = dst43_l;
dst43_l = dst65_l;
dst65_l = dst87_l;
}
src += 8;
dst += 8;
}
}
static void hevc_hv_8t_8w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
const int8_t * __restrict filter_x,
const int8_t * __restrict filter_y,
int32_t height)
{
hevc_hv_8t_8multx2mult_msa(src, src_stride, dst, dst_stride,
filter_x, filter_y, height, 8);
}
static void hevc_hv_8t_12w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
const int8_t * __restrict filter_x,
const int8_t * __restrict filter_y,
int32_t height)
{
hevc_hv_8t_8multx2mult_msa(src, src_stride, dst, dst_stride,
filter_x, filter_y, height, 8);
hevc_hv_8t_4w_msa(src + 8, src_stride, dst + 8, dst_stride,
filter_x, filter_y, height);
}
static void hevc_hv_8t_16w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
const int8_t * __restrict filter_x,
const int8_t * __restrict filter_y,
int32_t height)
{
hevc_hv_8t_8multx2mult_msa(src, src_stride, dst, dst_stride,
filter_x, filter_y, height, 16);
}
static void hevc_hv_8t_24w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
const int8_t * __restrict filter_x,
const int8_t * __restrict filter_y,
int32_t height)
{
hevc_hv_8t_8multx2mult_msa(src, src_stride, dst, dst_stride,
filter_x, filter_y, height, 24);
}
static void hevc_hv_8t_32w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
const int8_t * __restrict filter_x,
const int8_t * __restrict filter_y,
int32_t height)
{
hevc_hv_8t_8multx2mult_msa(src, src_stride, dst, dst_stride,
filter_x, filter_y, height, 32);
}
static void hevc_hv_8t_48w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
const int8_t * __restrict filter_x,
const int8_t * __restrict filter_y,
int32_t height)
{
hevc_hv_8t_8multx2mult_msa(src, src_stride, dst, dst_stride,
filter_x, filter_y, height, 48);
}
static void hevc_hv_8t_64w_msa(uint8_t * __restrict src, int32_t src_stride,
int16_t * __restrict dst, int32_t dst_stride,
const int8_t * __restrict filter_x,
const int8_t * __restrict filter_y,
int32_t height)
{
hevc_hv_8t_8multx2mult_msa(src, src_stride, dst, dst_stride,
filter_x, filter_y, height, 64);
}
#define MC_COPY(WIDTH) \
void ff_hevc_put_hevc_pel_pixels##WIDTH##_8_msa(int16_t *dst, \
uint8_t *src, \
ptrdiff_t src_stride, \
int height, \
intptr_t mx, \
intptr_t my, \
int width) \
{ \
hevc_copy_##WIDTH##w_msa(src, src_stride, dst, MAX_PB_SIZE, height); \
}
MC_COPY(4);
MC_COPY(6);
MC_COPY(8);
MC_COPY(12);
MC_COPY(16);
MC_COPY(24);
MC_COPY(32);
MC_COPY(48);
MC_COPY(64);
#undef MC_COPY
#define MC(PEL, DIR, WIDTH, TAP, DIR1, FILT_DIR) \
void ff_hevc_put_hevc_##PEL##_##DIR####WIDTH##_8_msa(int16_t *dst, \
uint8_t *src, \
ptrdiff_t src_stride, \
int height, \
intptr_t mx, \
intptr_t my, \
int width) \
{ \
const int8_t *filter = ff_hevc_##PEL##_filters[FILT_DIR - 1]; \
\
hevc_##DIR1##_##TAP##t_##WIDTH##w_msa(src, src_stride, dst, \
MAX_PB_SIZE, filter, height); \
}
MC(qpel, h, 4, 8, hz, mx);
MC(qpel, h, 8, 8, hz, mx);
MC(qpel, h, 12, 8, hz, mx);
MC(qpel, h, 16, 8, hz, mx);
MC(qpel, h, 24, 8, hz, mx);
MC(qpel, h, 32, 8, hz, mx);
MC(qpel, h, 48, 8, hz, mx);
MC(qpel, h, 64, 8, hz, mx);
MC(qpel, v, 4, 8, vt, my);
MC(qpel, v, 8, 8, vt, my);
MC(qpel, v, 12, 8, vt, my);
MC(qpel, v, 16, 8, vt, my);
MC(qpel, v, 24, 8, vt, my);
MC(qpel, v, 32, 8, vt, my);
MC(qpel, v, 48, 8, vt, my);
MC(qpel, v, 64, 8, vt, my);
#undef MC
#define MC_HV(PEL, DIR, WIDTH, TAP, DIR1) \
void ff_hevc_put_hevc_##PEL##_##DIR####WIDTH##_8_msa(int16_t *dst, \
uint8_t *src, \
ptrdiff_t src_stride, \
int height, \
intptr_t mx, \
intptr_t my, \
int width) \
{ \
const int8_t *filter_x = ff_hevc_##PEL##_filters[mx - 1]; \
const int8_t *filter_y = ff_hevc_##PEL##_filters[my - 1]; \
\
hevc_##DIR1##_##TAP##t_##WIDTH##w_msa(src, src_stride, dst, MAX_PB_SIZE, \
filter_x, filter_y, height); \
}
MC_HV(qpel, hv, 4, 8, hv);
MC_HV(qpel, hv, 8, 8, hv);
MC_HV(qpel, hv, 12, 8, hv);
MC_HV(qpel, hv, 16, 8, hv);
MC_HV(qpel, hv, 24, 8, hv);
MC_HV(qpel, hv, 32, 8, hv);
MC_HV(qpel, hv, 48, 8, hv);
MC_HV(qpel, hv, 64, 8, hv);
#undef MC_HV