1
0
mirror of https://github.com/FFmpeg/FFmpeg.git synced 2024-11-26 19:01:44 +02:00
FFmpeg/libavcodec/loongarch/vp9_lpf_lsx.c
Jin Bo fea299f876 avcodec: [loongarch] Optimize vp9_lpf/idct with LSX.
ffmpeg -i ../10_vp9_1080p_30fps_3Mbps.webm -f rawvideo -y /dev/null -an
before:294fps
after :567fps

Reviewed-by: Shiyou Yin <yinshiyou-hf@loongson.cn>
Signed-off-by: Michael Niedermayer <michael@niedermayer.cc>
2021-12-23 12:28:54 +01:00

3142 lines
125 KiB
C

/*
* Copyright (c) 2021 Loongson Technology Corporation Limited
* Contributed by Jin Bo <jinbo@loongson.cn>
*
* 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 "libavcodec/vp9dsp.h"
#include "libavutil/loongarch/loongson_intrinsics.h"
#include "libavutil/common.h"
#include "vp9dsp_loongarch.h"
#define LSX_LD_8(_src, _stride, _stride2, _stride3, _stride4, _in0, _in1, _in2, \
_in3, _in4, _in5, _in6, _in7) \
{ \
_in0 = __lsx_vld(_src, 0); \
_in1 = __lsx_vldx(_src, _stride); \
_in2 = __lsx_vldx(_src, _stride2); \
_in3 = __lsx_vldx(_src, _stride3); \
_src += _stride4; \
_in4 = __lsx_vld(_src, 0); \
_in5 = __lsx_vldx(_src, _stride); \
_in6 = __lsx_vldx(_src, _stride2); \
_in7 = __lsx_vldx(_src, _stride3); \
}
#define LSX_ST_8(_dst0, _dst1, _dst2, _dst3, _dst4, _dst5, _dst6, _dst7, \
_dst, _stride, _stride2, _stride3, _stride4) \
{ \
__lsx_vst(_dst0, _dst, 0); \
__lsx_vstx(_dst1, _dst, _stride); \
__lsx_vstx(_dst2, _dst, _stride2); \
__lsx_vstx(_dst3, _dst, _stride3); \
_dst += _stride4; \
__lsx_vst(_dst4, _dst, 0); \
__lsx_vstx(_dst5, _dst, _stride); \
__lsx_vstx(_dst6, _dst, _stride2); \
__lsx_vstx(_dst7, _dst, _stride3); \
}
#define VP9_LPF_FILTER4_4W(p1_src, p0_src, q0_src, q1_src, mask_src, hev_src, \
p1_dst, p0_dst, q0_dst, q1_dst) \
{ \
__m128i p1_tmp, p0_tmp, q0_tmp, q1_tmp, q0_sub_p0, filt, filt1, filt2; \
const __m128i cnst3b = __lsx_vldi(3); \
const __m128i cnst4b = __lsx_vldi(4); \
\
p1_tmp = __lsx_vxori_b(p1_src, 0x80); \
p0_tmp = __lsx_vxori_b(p0_src, 0x80); \
q0_tmp = __lsx_vxori_b(q0_src, 0x80); \
q1_tmp = __lsx_vxori_b(q1_src, 0x80); \
\
filt = __lsx_vssub_b(p1_tmp, q1_tmp); \
\
filt = filt & hev_src; \
\
q0_sub_p0 = __lsx_vssub_b(q0_tmp, p0_tmp); \
filt = __lsx_vsadd_b(filt, q0_sub_p0); \
filt = __lsx_vsadd_b(filt, q0_sub_p0); \
filt = __lsx_vsadd_b(filt, q0_sub_p0); \
filt = filt & mask_src; \
\
filt1 = __lsx_vsadd_b(filt, cnst4b); \
filt1 = __lsx_vsrai_b(filt1, 3); \
\
filt2 = __lsx_vsadd_b(filt, cnst3b); \
filt2 = __lsx_vsrai_b(filt2, 3); \
\
q0_tmp = __lsx_vssub_b(q0_tmp, filt1); \
q0_dst = __lsx_vxori_b(q0_tmp, 0x80); \
p0_tmp = __lsx_vsadd_b(p0_tmp, filt2); \
p0_dst = __lsx_vxori_b(p0_tmp, 0x80); \
\
filt = __lsx_vsrari_b(filt1, 1); \
hev_src = __lsx_vxori_b(hev_src, 0xff); \
filt = filt & hev_src; \
\
q1_tmp = __lsx_vssub_b(q1_tmp, filt); \
q1_dst = __lsx_vxori_b(q1_tmp, 0x80); \
p1_tmp = __lsx_vsadd_b(p1_tmp, filt); \
p1_dst = __lsx_vxori_b(p1_tmp, 0x80); \
}
#define VP9_FLAT4(p3_src, p2_src, p0_src, q0_src, q2_src, q3_src, flat_dst) \
{ \
__m128i f_tmp = __lsx_vldi(1); \
__m128i p2_a_sub_p0, q2_a_sub_q0, p3_a_sub_p0, q3_a_sub_q0; \
\
p2_a_sub_p0 = __lsx_vabsd_bu(p2_src, p0_src); \
q2_a_sub_q0 = __lsx_vabsd_bu(q2_src, q0_src); \
p3_a_sub_p0 = __lsx_vabsd_bu(p3_src, p0_src); \
q3_a_sub_q0 = __lsx_vabsd_bu(q3_src, q0_src); \
\
p2_a_sub_p0 = __lsx_vmax_bu(p2_a_sub_p0, q2_a_sub_q0); \
flat_dst = __lsx_vmax_bu(p2_a_sub_p0, flat_dst); \
p3_a_sub_p0 = __lsx_vmax_bu(p3_a_sub_p0, q3_a_sub_q0); \
flat_dst = __lsx_vmax_bu(p3_a_sub_p0, flat_dst); \
\
flat_dst = __lsx_vslt_bu(f_tmp, flat_dst); \
flat_dst = __lsx_vxori_b(flat_dst, 0xff); \
flat_dst = flat_dst & mask; \
}
#define VP9_FLAT5(p7_src, p6_src, p5_src, p4_src, p0_src, q0_src, q4_src, \
q5_src, q6_src, q7_src, flat_src, flat2_dst) \
{ \
__m128i f_tmp = __lsx_vldi(1); \
__m128i p4_a_sub_p0, q4_a_sub_q0, p5_a_sub_p0, q5_a_sub_q0; \
__m128i p6_a_sub_p0, q6_a_sub_q0, p7_a_sub_p0, q7_a_sub_q0; \
\
p4_a_sub_p0 = __lsx_vabsd_bu(p4_src, p0_src); \
q4_a_sub_q0 = __lsx_vabsd_bu(q4_src, q0_src); \
p5_a_sub_p0 = __lsx_vabsd_bu(p5_src, p0_src); \
q5_a_sub_q0 = __lsx_vabsd_bu(q5_src, q0_src); \
p6_a_sub_p0 = __lsx_vabsd_bu(p6_src, p0_src); \
q6_a_sub_q0 = __lsx_vabsd_bu(q6_src, q0_src); \
p7_a_sub_p0 = __lsx_vabsd_bu(p7_src, p0_src); \
q7_a_sub_q0 = __lsx_vabsd_bu(q7_src, q0_src); \
\
p4_a_sub_p0 = __lsx_vmax_bu(p4_a_sub_p0, q4_a_sub_q0); \
flat2_dst = __lsx_vmax_bu(p5_a_sub_p0, q5_a_sub_q0); \
flat2_dst = __lsx_vmax_bu(p4_a_sub_p0, flat2_dst); \
p6_a_sub_p0 = __lsx_vmax_bu(p6_a_sub_p0, q6_a_sub_q0); \
flat2_dst = __lsx_vmax_bu(p6_a_sub_p0, flat2_dst); \
p7_a_sub_p0 = __lsx_vmax_bu(p7_a_sub_p0, q7_a_sub_q0); \
flat2_dst = __lsx_vmax_bu(p7_a_sub_p0, flat2_dst); \
\
flat2_dst = __lsx_vslt_bu(f_tmp, flat2_dst); \
flat2_dst = __lsx_vxori_b(flat2_dst, 0xff); \
flat2_dst = flat2_dst & flat_src; \
}
#define VP9_FILTER8(p3_src, p2_src, p1_src, p0_src, \
q0_src, q1_src, q2_src, q3_src, \
p2_filt8_dst, p1_filt8_dst, p0_filt8_dst, \
q0_filt8_dst, q1_filt8_dst, q2_filt8_dst) \
{ \
__m128i tmp0, tmp1, tmp2; \
\
tmp2 = __lsx_vadd_h(p2_src, p1_src); \
tmp2 = __lsx_vadd_h(tmp2, p0_src); \
tmp0 = __lsx_vslli_h(p3_src, 1); \
\
tmp0 = __lsx_vadd_h(tmp0, tmp2); \
tmp0 = __lsx_vadd_h(tmp0, q0_src); \
tmp1 = __lsx_vadd_h(tmp0, p3_src); \
tmp1 = __lsx_vadd_h(tmp1, p2_src); \
p2_filt8_dst = __lsx_vsrari_h(tmp1, 3); \
\
tmp1 = __lsx_vadd_h(tmp0, p1_src); \
tmp1 = __lsx_vadd_h(tmp1, q1_src); \
p1_filt8_dst = __lsx_vsrari_h(tmp1, 3); \
\
tmp1 = __lsx_vadd_h(q2_src, q1_src); \
tmp1 = __lsx_vadd_h(tmp1, q0_src); \
tmp2 = __lsx_vadd_h(tmp2, tmp1); \
tmp0 = __lsx_vadd_h(tmp2, p0_src); \
tmp0 = __lsx_vadd_h(tmp0, p3_src); \
p0_filt8_dst = __lsx_vsrari_h(tmp0, 3); \
\
tmp0 = __lsx_vadd_h(q2_src, q3_src); \
tmp0 = __lsx_vadd_h(tmp0, p0_src); \
tmp0 = __lsx_vadd_h(tmp0, tmp1); \
tmp1 = __lsx_vadd_h(q3_src, q3_src); \
tmp1 = __lsx_vadd_h(tmp1, tmp0); \
q2_filt8_dst = __lsx_vsrari_h(tmp1, 3); \
\
tmp0 = __lsx_vadd_h(tmp2, q3_src); \
tmp1 = __lsx_vadd_h(tmp0, q0_src); \
q0_filt8_dst = __lsx_vsrari_h(tmp1, 3); \
\
tmp1 = __lsx_vsub_h(tmp0, p2_src); \
tmp0 = __lsx_vadd_h(q1_src, q3_src); \
tmp1 = __lsx_vadd_h(tmp0, tmp1); \
q1_filt8_dst = __lsx_vsrari_h(tmp1, 3); \
}
#define LPF_MASK_HEV(p3_src, p2_src, p1_src, p0_src, q0_src, q1_src, \
q2_src, q3_src, limit_src, b_limit_src, thresh_src, \
hev_dst, mask_dst, flat_dst) \
{ \
__m128i p3_asub_p2_tmp, p2_asub_p1_tmp, p1_asub_p0_tmp, q1_asub_q0_tmp; \
__m128i p1_asub_q1_tmp, p0_asub_q0_tmp, q3_asub_q2_tmp, q2_asub_q1_tmp; \
\
/* absolute subtraction of pixel values */ \
p3_asub_p2_tmp = __lsx_vabsd_bu(p3_src, p2_src); \
p2_asub_p1_tmp = __lsx_vabsd_bu(p2_src, p1_src); \
p1_asub_p0_tmp = __lsx_vabsd_bu(p1_src, p0_src); \
q1_asub_q0_tmp = __lsx_vabsd_bu(q1_src, q0_src); \
q2_asub_q1_tmp = __lsx_vabsd_bu(q2_src, q1_src); \
q3_asub_q2_tmp = __lsx_vabsd_bu(q3_src, q2_src); \
p0_asub_q0_tmp = __lsx_vabsd_bu(p0_src, q0_src); \
p1_asub_q1_tmp = __lsx_vabsd_bu(p1_src, q1_src); \
\
/* calculation of hev */ \
flat_dst = __lsx_vmax_bu(p1_asub_p0_tmp, q1_asub_q0_tmp); \
hev_dst = __lsx_vslt_bu(thresh_src, flat_dst); \
\
/* calculation of mask */ \
p0_asub_q0_tmp = __lsx_vsadd_bu(p0_asub_q0_tmp, p0_asub_q0_tmp); \
p1_asub_q1_tmp = __lsx_vsrli_b(p1_asub_q1_tmp, 1); \
p0_asub_q0_tmp = __lsx_vsadd_bu(p0_asub_q0_tmp, p1_asub_q1_tmp); \
\
mask_dst = __lsx_vslt_bu(b_limit_src, p0_asub_q0_tmp); \
mask_dst = __lsx_vmax_bu(flat_dst, mask_dst); \
p3_asub_p2_tmp = __lsx_vmax_bu(p3_asub_p2_tmp, p2_asub_p1_tmp); \
mask_dst = __lsx_vmax_bu(p3_asub_p2_tmp, mask_dst); \
q2_asub_q1_tmp = __lsx_vmax_bu(q2_asub_q1_tmp, q3_asub_q2_tmp); \
mask_dst = __lsx_vmax_bu(q2_asub_q1_tmp, mask_dst); \
\
mask_dst = __lsx_vslt_bu(limit_src, mask_dst); \
mask_dst = __lsx_vxori_b(mask_dst, 0xff); \
}
void ff_loop_filter_v_4_8_lsx(uint8_t *dst, ptrdiff_t stride,
int32_t b_limit_ptr,
int32_t limit_ptr,
int32_t thresh_ptr)
{
ptrdiff_t stride2 = stride << 1;
ptrdiff_t stride3 = stride2 + stride;
ptrdiff_t stride4 = stride2 << 1;
__m128i mask, hev, flat, thresh, b_limit, limit;
__m128i p3, p2, p1, p0, q3, q2, q1, q0, p1_out, p0_out, q0_out, q1_out;
DUP4_ARG2(__lsx_vldx, dst, -stride4, dst, -stride3, dst, -stride2,
dst, -stride, p3, p2, p1, p0);
q0 = __lsx_vld(dst, 0);
DUP2_ARG2(__lsx_vldx, dst, stride, dst, stride2, q1, q2);
q3 = __lsx_vldx(dst, stride3);
thresh = __lsx_vreplgr2vr_b(thresh_ptr);
b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
limit = __lsx_vreplgr2vr_b(limit_ptr);
LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
hev, mask, flat);
VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
q1_out);
__lsx_vstelm_d(p1_out, dst - stride2, 0, 0);
__lsx_vstelm_d(p0_out, dst - stride, 0, 0);
__lsx_vstelm_d(q0_out, dst , 0, 0);
__lsx_vstelm_d(q1_out, dst + stride, 0, 0);
}
void ff_loop_filter_v_44_16_lsx(uint8_t *dst, ptrdiff_t stride,
int32_t b_limit_ptr,
int32_t limit_ptr,
int32_t thresh_ptr)
{
ptrdiff_t stride2 = stride << 1;
ptrdiff_t stride3 = stride2 + stride;
ptrdiff_t stride4 = stride2 << 1;
__m128i mask, hev, flat, thresh0, b_limit0;
__m128i limit0, thresh1, b_limit1, limit1;
__m128i p3, p2, p1, p0, q3, q2, q1, q0;
DUP4_ARG2(__lsx_vldx, dst, -stride4, dst, -stride3, dst, -stride2,
dst, -stride, p3, p2, p1, p0);
q0 = __lsx_vld(dst, 0);
DUP2_ARG2(__lsx_vldx, dst, stride, dst, stride2, q1, q2);
q3 = __lsx_vldx(dst, stride3);
thresh0 = __lsx_vreplgr2vr_b(thresh_ptr);
thresh1 = __lsx_vreplgr2vr_b(thresh_ptr >> 8);
thresh0 = __lsx_vilvl_d(thresh1, thresh0);
b_limit0 = __lsx_vreplgr2vr_b(b_limit_ptr);
b_limit1 = __lsx_vreplgr2vr_b(b_limit_ptr >> 8);
b_limit0 = __lsx_vilvl_d(b_limit1, b_limit0);
limit0 = __lsx_vreplgr2vr_b(limit_ptr);
limit1 = __lsx_vreplgr2vr_b(limit_ptr >> 8);
limit0 = __lsx_vilvl_d(limit1, limit0);
LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit0, b_limit0, thresh0,
hev, mask, flat);
VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1, p0, q0, q1);
__lsx_vst(p1, dst - stride2, 0);
__lsx_vst(p0, dst - stride, 0);
__lsx_vst(q0, dst , 0);
__lsx_vst(q1, dst + stride, 0);
}
void ff_loop_filter_v_8_8_lsx(uint8_t *dst, ptrdiff_t stride,
int32_t b_limit_ptr,
int32_t limit_ptr,
int32_t thresh_ptr)
{
ptrdiff_t stride2 = stride << 1;
ptrdiff_t stride3 = stride2 + stride;
ptrdiff_t stride4 = stride2 << 1;
__m128i mask, hev, flat, thresh, b_limit, limit;
__m128i p3, p2, p1, p0, q3, q2, q1, q0;
__m128i p2_out, p1_out, p0_out, q0_out, q1_out, q2_out;
__m128i p2_filter8, p1_filter8, p0_filter8;
__m128i q0_filter8, q1_filter8, q2_filter8;
__m128i p3_l, p2_l, p1_l, p0_l, q3_l, q2_l, q1_l, q0_l;
__m128i zero = __lsx_vldi(0);
DUP4_ARG2(__lsx_vldx, dst, -stride4, dst, -stride3, dst, -stride2,
dst, -stride, p3, p2, p1, p0);
q0 = __lsx_vld(dst, 0);
DUP2_ARG2(__lsx_vldx, dst, stride, dst, stride2, q1, q2);
q3 = __lsx_vldx(dst, stride3);
thresh = __lsx_vreplgr2vr_b(thresh_ptr);
b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
limit = __lsx_vreplgr2vr_b(limit_ptr);
LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
hev, mask, flat);
VP9_FLAT4(p3, p2, p0, q0, q2, q3, flat);
VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
q1_out);
flat = __lsx_vilvl_d(zero, flat);
/* if flat is zero for all pixels, then no need to calculate other filter */
if (__lsx_bz_v(flat)) {
__lsx_vstelm_d(p1_out, dst - stride2, 0, 0);
__lsx_vstelm_d(p0_out, dst - stride, 0, 0);
__lsx_vstelm_d(q0_out, dst , 0, 0);
__lsx_vstelm_d(q1_out, dst + stride, 0, 0);
} else {
DUP4_ARG2(__lsx_vilvl_b, zero, p3, zero, p2, zero, p1, zero, p0,
p3_l, p2_l, p1_l, p0_l);
DUP4_ARG2(__lsx_vilvl_b, zero, q0, zero, q1, zero, q2, zero, q3,
q0_l, q1_l, q2_l, q3_l);
VP9_FILTER8(p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l, p2_filter8,
p1_filter8, p0_filter8, q0_filter8, q1_filter8, q2_filter8);
/* convert 16 bit output data into 8 bit */
DUP4_ARG2(__lsx_vpickev_b, zero, p2_filter8, zero, p1_filter8,
zero, p0_filter8, zero, q0_filter8, p2_filter8,
p1_filter8, p0_filter8, q0_filter8);
DUP2_ARG2(__lsx_vpickev_b, zero, q1_filter8, zero, q2_filter8,
q1_filter8, q2_filter8);
/* store pixel values */
p2_out = __lsx_vbitsel_v(p2, p2_filter8, flat);
p1_out = __lsx_vbitsel_v(p1_out, p1_filter8, flat);
p0_out = __lsx_vbitsel_v(p0_out, p0_filter8, flat);
q0_out = __lsx_vbitsel_v(q0_out, q0_filter8, flat);
q1_out = __lsx_vbitsel_v(q1_out, q1_filter8, flat);
q2_out = __lsx_vbitsel_v(q2, q2_filter8, flat);
__lsx_vstelm_d(p2_out, dst - stride3, 0, 0);
__lsx_vstelm_d(p1_out, dst - stride2, 0, 0);
__lsx_vstelm_d(p0_out, dst - stride, 0, 0);
__lsx_vstelm_d(q0_out, dst, 0, 0);
__lsx_vstelm_d(q1_out, dst + stride, 0, 0);
__lsx_vstelm_d(q2_out, dst + stride2, 0, 0);
}
}
void ff_loop_filter_v_88_16_lsx(uint8_t *dst, ptrdiff_t stride,
int32_t b_limit_ptr,
int32_t limit_ptr,
int32_t thresh_ptr)
{
ptrdiff_t stride2 = stride << 1;
ptrdiff_t stride3 = stride2 + stride;
ptrdiff_t stride4 = stride2 << 1;
__m128i p3, p2, p1, p0, q3, q2, q1, q0;
__m128i p2_out, p1_out, p0_out, q0_out, q1_out, q2_out;
__m128i flat, mask, hev, tmp, thresh, b_limit, limit;
__m128i p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l;
__m128i p3_h, p2_h, p1_h, p0_h, q0_h, q1_h, q2_h, q3_h;
__m128i p2_filt8_l, p1_filt8_l, p0_filt8_l;
__m128i q0_filt8_l, q1_filt8_l, q2_filt8_l;
__m128i p2_filt8_h, p1_filt8_h, p0_filt8_h;
__m128i q0_filt8_h, q1_filt8_h, q2_filt8_h;
__m128i zero = __lsx_vldi(0);
/* load vector elements */
DUP4_ARG2(__lsx_vldx, dst, -stride4, dst, -stride3, dst, -stride2,
dst, -stride, p3, p2, p1, p0);
q0 = __lsx_vld(dst, 0);
DUP2_ARG2(__lsx_vldx, dst, stride, dst, stride2, q1, q2);
q3 = __lsx_vldx(dst, stride3);
thresh = __lsx_vreplgr2vr_b(thresh_ptr);
tmp = __lsx_vreplgr2vr_b(thresh_ptr >> 8);
thresh = __lsx_vilvl_d(tmp, thresh);
b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
tmp = __lsx_vreplgr2vr_b(b_limit_ptr >> 8);
b_limit = __lsx_vilvl_d(tmp, b_limit);
limit = __lsx_vreplgr2vr_b(limit_ptr);
tmp = __lsx_vreplgr2vr_b(limit_ptr >> 8);
limit = __lsx_vilvl_d(tmp, limit);
/* mask and hev */
LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
hev, mask, flat);
VP9_FLAT4(p3, p2, p0, q0, q2, q3, flat);
VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
q1_out);
/* if flat is zero for all pixels, then no need to calculate other filter */
if (__lsx_bz_v(flat)) {
__lsx_vst(p1_out, dst - stride2, 0);
__lsx_vst(p0_out, dst - stride, 0);
__lsx_vst(q0_out, dst, 0);
__lsx_vst(q1_out, dst + stride, 0);
} else {
DUP4_ARG2(__lsx_vilvl_b, zero, p3, zero, p2, zero, p1, zero, p0,
p3_l, p2_l, p1_l, p0_l);
DUP4_ARG2(__lsx_vilvl_b, zero, q0, zero, q1, zero, q2, zero, q3,
q0_l, q1_l, q2_l, q3_l);
VP9_FILTER8(p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l, p2_filt8_l,
p1_filt8_l, p0_filt8_l, q0_filt8_l, q1_filt8_l, q2_filt8_l);
DUP4_ARG2(__lsx_vilvh_b, zero, p3, zero, p2, zero, p1, zero, p0,
p3_h, p2_h, p1_h, p0_h);
DUP4_ARG2(__lsx_vilvh_b, zero, q0, zero, q1, zero, q2, zero, q3,
q0_h, q1_h, q2_h, q3_h);
VP9_FILTER8(p3_h, p2_h, p1_h, p0_h, q0_h, q1_h, q2_h, q3_h, p2_filt8_h,
p1_filt8_h, p0_filt8_h, q0_filt8_h, q1_filt8_h, q2_filt8_h);
/* convert 16 bit output data into 8 bit */
DUP4_ARG2(__lsx_vpickev_b, p2_filt8_h, p2_filt8_l, p1_filt8_h,
p1_filt8_l, p0_filt8_h, p0_filt8_l, q0_filt8_h, q0_filt8_l,
p2_filt8_l, p1_filt8_l, p0_filt8_l, q0_filt8_l);
DUP2_ARG2(__lsx_vpickev_b, q1_filt8_h, q1_filt8_l, q2_filt8_h,
q2_filt8_l, q1_filt8_l, q2_filt8_l);
/* store pixel values */
p2_out = __lsx_vbitsel_v(p2, p2_filt8_l, flat);
p1_out = __lsx_vbitsel_v(p1_out, p1_filt8_l, flat);
p0_out = __lsx_vbitsel_v(p0_out, p0_filt8_l, flat);
q0_out = __lsx_vbitsel_v(q0_out, q0_filt8_l, flat);
q1_out = __lsx_vbitsel_v(q1_out, q1_filt8_l, flat);
q2_out = __lsx_vbitsel_v(q2, q2_filt8_l, flat);
__lsx_vstx(p2_out, dst, -stride3);
__lsx_vstx(p1_out, dst, -stride2);
__lsx_vstx(p0_out, dst, -stride);
__lsx_vst(q0_out, dst, 0);
__lsx_vstx(q1_out, dst, stride);
__lsx_vstx(q2_out, dst, stride2);
}
}
void ff_loop_filter_v_84_16_lsx(uint8_t *dst, ptrdiff_t stride,
int32_t b_limit_ptr,
int32_t limit_ptr,
int32_t thresh_ptr)
{
ptrdiff_t stride2 = stride << 1;
ptrdiff_t stride3 = stride2 + stride;
ptrdiff_t stride4 = stride2 << 1;
__m128i p3, p2, p1, p0, q3, q2, q1, q0;
__m128i p2_out, p1_out, p0_out, q0_out, q1_out, q2_out;
__m128i flat, mask, hev, tmp, thresh, b_limit, limit;
__m128i p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l;
__m128i p2_filt8_l, p1_filt8_l, p0_filt8_l;
__m128i q0_filt8_l, q1_filt8_l, q2_filt8_l;
__m128i zero = __lsx_vldi(0);
/* load vector elements */
DUP4_ARG2(__lsx_vldx, dst, -stride4, dst, -stride3, dst, -stride2,
dst, -stride, p3, p2, p1, p0);
q0 = __lsx_vld(dst, 0);
DUP2_ARG2(__lsx_vldx, dst, stride, dst, stride2, q1, q2);
q3 = __lsx_vldx(dst, stride3);
thresh = __lsx_vreplgr2vr_b(thresh_ptr);
tmp = __lsx_vreplgr2vr_b(thresh_ptr >> 8);
thresh = __lsx_vilvl_d(tmp, thresh);
b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
tmp = __lsx_vreplgr2vr_b(b_limit_ptr >> 8);
b_limit = __lsx_vilvl_d(tmp, b_limit);
limit = __lsx_vreplgr2vr_b(limit_ptr);
tmp = __lsx_vreplgr2vr_b(limit_ptr >> 8);
limit = __lsx_vilvl_d(tmp, limit);
/* mask and hev */
LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
hev, mask, flat);
VP9_FLAT4(p3, p2, p0, q0, q2, q3, flat);
VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
q1_out);
flat = __lsx_vilvl_d(zero, flat);
/* if flat is zero for all pixels, then no need to calculate other filter */
if (__lsx_bz_v(flat)) {
__lsx_vstx(p1_out, dst, -stride2);
__lsx_vstx(p0_out, dst, -stride);
__lsx_vst(q0_out, dst, 0);
__lsx_vstx(q1_out, dst, stride);
} else {
DUP4_ARG2(__lsx_vilvl_b, zero, p3, zero, p2, zero, p1, zero, p0,
p3_l, p2_l, p1_l, p0_l);
DUP4_ARG2(__lsx_vilvl_b, zero, q0, zero, q1, zero, q2, zero, q3,
q0_l, q1_l, q2_l, q3_l);
VP9_FILTER8(p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l, p2_filt8_l,
p1_filt8_l, p0_filt8_l, q0_filt8_l, q1_filt8_l, q2_filt8_l);
/* convert 16 bit output data into 8 bit */
DUP4_ARG2(__lsx_vpickev_b, p2_filt8_l, p2_filt8_l, p1_filt8_l,
p1_filt8_l, p0_filt8_l, p0_filt8_l, q0_filt8_l, q0_filt8_l,
p2_filt8_l, p1_filt8_l, p0_filt8_l, q0_filt8_l);
DUP2_ARG2(__lsx_vpickev_b, q1_filt8_l, q1_filt8_l, q2_filt8_l,
q2_filt8_l, q1_filt8_l, q2_filt8_l);
/* store pixel values */
p2_out = __lsx_vbitsel_v(p2, p2_filt8_l, flat);
p1_out = __lsx_vbitsel_v(p1_out, p1_filt8_l, flat);
p0_out = __lsx_vbitsel_v(p0_out, p0_filt8_l, flat);
q0_out = __lsx_vbitsel_v(q0_out, q0_filt8_l, flat);
q1_out = __lsx_vbitsel_v(q1_out, q1_filt8_l, flat);
q2_out = __lsx_vbitsel_v(q2, q2_filt8_l, flat);
__lsx_vstx(p2_out, dst, -stride3);
__lsx_vstx(p1_out, dst, -stride2);
__lsx_vstx(p0_out, dst, -stride);
__lsx_vst(q0_out, dst, 0);
__lsx_vstx(q1_out, dst, stride);
__lsx_vstx(q2_out, dst, stride2);
}
}
void ff_loop_filter_v_48_16_lsx(uint8_t *dst, ptrdiff_t stride,
int32_t b_limit_ptr,
int32_t limit_ptr,
int32_t thresh_ptr)
{
ptrdiff_t stride2 = stride << 1;
ptrdiff_t stride3 = stride2 + stride;
ptrdiff_t stride4 = stride2 << 1;
__m128i p3, p2, p1, p0, q3, q2, q1, q0;
__m128i p2_out, p1_out, p0_out, q0_out, q1_out, q2_out;
__m128i flat, mask, hev, tmp, thresh, b_limit, limit;
__m128i p3_h, p2_h, p1_h, p0_h, q0_h, q1_h, q2_h, q3_h;
__m128i p2_filt8_h, p1_filt8_h, p0_filt8_h;
__m128i q0_filt8_h, q1_filt8_h, q2_filt8_h;
__m128i zero = { 0 };
/* load vector elements */
DUP4_ARG2(__lsx_vldx, dst, -stride4, dst, -stride3, dst, -stride2,
dst, -stride, p3, p2, p1, p0);
q0 = __lsx_vld(dst, 0);
DUP2_ARG2(__lsx_vldx, dst, stride, dst, stride2, q1, q2);
q3 = __lsx_vldx(dst, stride3);
thresh = __lsx_vreplgr2vr_b(thresh_ptr);
tmp = __lsx_vreplgr2vr_b(thresh_ptr >> 8);
thresh = __lsx_vilvl_d(tmp, thresh);
b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
tmp = __lsx_vreplgr2vr_b(b_limit_ptr >> 8);
b_limit = __lsx_vilvl_d(tmp, b_limit);
limit = __lsx_vreplgr2vr_b(limit_ptr);
tmp = __lsx_vreplgr2vr_b(limit_ptr >> 8);
limit = __lsx_vilvl_d(tmp, limit);
/* mask and hev */
LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
hev, mask, flat);
VP9_FLAT4(p3, p2, p0, q0, q2, q3, flat);
VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
q1_out);
flat = __lsx_vilvh_d(flat, zero);
/* if flat is zero for all pixels, then no need to calculate other filter */
if (__lsx_bz_v(flat)) {
__lsx_vstx(p1_out, dst, -stride2);
__lsx_vstx(p0_out, dst, -stride);
__lsx_vst(q0_out, dst, 0);
__lsx_vstx(q1_out, dst, stride);
} else {
DUP4_ARG2(__lsx_vilvh_b, zero, p3, zero, p2, zero, p1, zero, p0,
p3_h, p2_h, p1_h, p0_h);
DUP4_ARG2(__lsx_vilvh_b, zero, q0, zero, q1, zero, q2, zero, q3,
q0_h, q1_h, q2_h, q3_h);
VP9_FILTER8(p3_h, p2_h, p1_h, p0_h, q0_h, q1_h, q2_h, q3_h, p2_filt8_h,
p1_filt8_h, p0_filt8_h, q0_filt8_h, q1_filt8_h, q2_filt8_h);
/* convert 16 bit output data into 8 bit */
DUP4_ARG2(__lsx_vpickev_b, p2_filt8_h, p2_filt8_h, p1_filt8_h,
p1_filt8_h, p0_filt8_h, p0_filt8_h, q0_filt8_h, q0_filt8_h,
p2_filt8_h, p1_filt8_h, p0_filt8_h, q0_filt8_h);
DUP2_ARG2(__lsx_vpickev_b, q1_filt8_h, q1_filt8_h, q2_filt8_h,
q2_filt8_h, q1_filt8_h, q2_filt8_h);
/* store pixel values */
p2_out = __lsx_vbitsel_v(p2, p2_filt8_h, flat);
p1_out = __lsx_vbitsel_v(p1_out, p1_filt8_h, flat);
p0_out = __lsx_vbitsel_v(p0_out, p0_filt8_h, flat);
q0_out = __lsx_vbitsel_v(q0_out, q0_filt8_h, flat);
q1_out = __lsx_vbitsel_v(q1_out, q1_filt8_h, flat);
q2_out = __lsx_vbitsel_v(q2, q2_filt8_h, flat);
__lsx_vstx(p2_out, dst, -stride3);
__lsx_vstx(p1_out, dst, -stride2);
__lsx_vstx(p0_out, dst, -stride);
__lsx_vst(q0_out, dst, 0);
__lsx_vstx(q1_out, dst, stride);
__lsx_vstx(q2_out, dst, stride2);
}
}
static int32_t vp9_hz_lpf_t4_and_t8_16w(uint8_t *dst, ptrdiff_t stride,
uint8_t *filter48,
int32_t b_limit_ptr,
int32_t limit_ptr,
int32_t thresh_ptr)
{
ptrdiff_t stride2 = stride << 1;
ptrdiff_t stride3 = stride2 + stride;
ptrdiff_t stride4 = stride2 << 1;
__m128i p3, p2, p1, p0, q3, q2, q1, q0;
__m128i p2_out, p1_out, p0_out, q0_out, q1_out, q2_out;
__m128i flat, mask, hev, thresh, b_limit, limit;
__m128i p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l;
__m128i p3_h, p2_h, p1_h, p0_h, q0_h, q1_h, q2_h, q3_h;
__m128i p2_filt8_l, p1_filt8_l, p0_filt8_l;
__m128i q0_filt8_l, q1_filt8_l, q2_filt8_l;
__m128i p2_filt8_h, p1_filt8_h, p0_filt8_h;
__m128i q0_filt8_h, q1_filt8_h, q2_filt8_h;
__m128i zero = __lsx_vldi(0);
/* load vector elements */
DUP4_ARG2(__lsx_vldx, dst, -stride4, dst, -stride3, dst, -stride2,
dst, -stride, p3, p2, p1, p0);
q0 = __lsx_vld(dst, 0);
DUP2_ARG2(__lsx_vldx, dst, stride, dst, stride2, q1, q2);
q3 = __lsx_vldx(dst, stride3);
thresh = __lsx_vreplgr2vr_b(thresh_ptr);
b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
limit = __lsx_vreplgr2vr_b(limit_ptr);
/* mask and hev */
LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
hev, mask, flat);
VP9_FLAT4(p3, p2, p0, q0, q2, q3, flat);
VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
q1_out);
/* if flat is zero for all pixels, then no need to calculate other filter */
if (__lsx_bz_v(flat)) {
__lsx_vstx(p1_out, dst, -stride2);
__lsx_vstx(p0_out, dst, -stride);
__lsx_vst(q0_out, dst, 0);
__lsx_vstx(q1_out, dst, stride);
return 1;
} else {
DUP4_ARG2(__lsx_vilvl_b, zero, p3, zero, p2, zero, p1, zero, p0,
p3_l, p2_l, p1_l, p0_l);
DUP4_ARG2(__lsx_vilvl_b, zero, q0, zero, q1, zero, q2, zero, q3,
q0_l, q1_l, q2_l, q3_l);
VP9_FILTER8(p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l, p2_filt8_l,
p1_filt8_l, p0_filt8_l, q0_filt8_l, q1_filt8_l, q2_filt8_l);
DUP4_ARG2(__lsx_vilvh_b, zero, p3, zero, p2, zero, p1, zero, p0,
p3_h, p2_h, p1_h, p0_h);
DUP4_ARG2(__lsx_vilvh_b, zero, q0, zero, q1, zero, q2, zero, q3,
q0_h, q1_h, q2_h, q3_h);
VP9_FILTER8(p3_h, p2_h, p1_h, p0_h, q0_h, q1_h, q2_h, q3_h, p2_filt8_h,
p1_filt8_h, p0_filt8_h, q0_filt8_h, q1_filt8_h, q2_filt8_h);
/* convert 16 bit output data into 8 bit */
DUP4_ARG2(__lsx_vpickev_b, p2_filt8_h, p2_filt8_l, p1_filt8_h,
p1_filt8_l, p0_filt8_h, p0_filt8_l, q0_filt8_h, q0_filt8_l,
p2_filt8_l, p1_filt8_l, p0_filt8_l, q0_filt8_l);
DUP2_ARG2(__lsx_vpickev_b, q1_filt8_h, q1_filt8_l, q2_filt8_h,
q2_filt8_l, q1_filt8_l, q2_filt8_l);
/* store pixel values */
p2_out = __lsx_vbitsel_v(p2, p2_filt8_l, flat);
p1_out = __lsx_vbitsel_v(p1_out, p1_filt8_l, flat);
p0_out = __lsx_vbitsel_v(p0_out, p0_filt8_l, flat);
q0_out = __lsx_vbitsel_v(q0_out, q0_filt8_l, flat);
q1_out = __lsx_vbitsel_v(q1_out, q1_filt8_l, flat);
q2_out = __lsx_vbitsel_v(q2, q2_filt8_l, flat);
__lsx_vst(p2_out, filter48, 0);
__lsx_vst(p1_out, filter48, 16);
__lsx_vst(p0_out, filter48, 32);
__lsx_vst(q0_out, filter48, 48);
__lsx_vst(q1_out, filter48, 64);
__lsx_vst(q2_out, filter48, 80);
__lsx_vst(flat, filter48, 96);
return 0;
}
}
static void vp9_hz_lpf_t16_16w(uint8_t *dst, ptrdiff_t stride,
uint8_t *filter48)
{
ptrdiff_t stride2 = stride << 1;
ptrdiff_t stride3 = stride2 + stride;
ptrdiff_t stride4 = stride2 << 1;
uint8_t *dst_tmp = dst - stride4;
uint8_t *dst_tmp1 = dst + stride4;
__m128i p7, p6, p5, p4, p3, p2, p1, p0, q0, q1, q2, q3, q4, q5, q6, q7;
__m128i flat, flat2, filter8;
__m128i zero = __lsx_vldi(0);
__m128i out_h, out_l;
v8u16 p7_l_in, p6_l_in, p5_l_in, p4_l_in;
v8u16 p3_l_in, p2_l_in, p1_l_in, p0_l_in;
v8u16 q7_l_in, q6_l_in, q5_l_in, q4_l_in;
v8u16 q3_l_in, q2_l_in, q1_l_in, q0_l_in;
v8u16 p7_h_in, p6_h_in, p5_h_in, p4_h_in;
v8u16 p3_h_in, p2_h_in, p1_h_in, p0_h_in;
v8u16 q7_h_in, q6_h_in, q5_h_in, q4_h_in;
v8u16 q3_h_in, q2_h_in, q1_h_in, q0_h_in;
v8u16 tmp0_l, tmp1_l, tmp0_h, tmp1_h;
flat = __lsx_vld(filter48, 96);
DUP4_ARG2(__lsx_vldx, dst_tmp, -stride4, dst_tmp, -stride3, dst_tmp,
-stride2, dst_tmp, -stride, p7, p6, p5, p4);
p3 = __lsx_vld(dst_tmp, 0);
DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, p2, p1);
p0 = __lsx_vldx(dst_tmp, stride3);
q0 = __lsx_vld(dst, 0);
DUP2_ARG2(__lsx_vldx, dst, stride, dst, stride2, q1, q2);
q3 = __lsx_vldx(dst, stride3);
q4 = __lsx_vld(dst_tmp1, 0);
DUP2_ARG2(__lsx_vldx, dst_tmp1, stride, dst_tmp1, stride2, q5, q6);
q7 = __lsx_vldx(dst_tmp1, stride3);
VP9_FLAT5(p7, p6, p5, p4, p0, q0, q4, q5, q6, q7, flat, flat2);
/* if flat2 is zero for all pixels, then no need to calculate other filter */
if (__lsx_bz_v(flat2)) {
DUP4_ARG2(__lsx_vld, filter48, 0, filter48, 16, filter48, 32, filter48,
48, p2, p1, p0, q0);
DUP2_ARG2(__lsx_vld, filter48, 64, filter48, 80, q1, q2);
__lsx_vstx(p2, dst, -stride3);
__lsx_vstx(p1, dst, -stride2);
__lsx_vstx(p0, dst, -stride);
__lsx_vst(q0, dst, 0);
__lsx_vstx(q1, dst, stride);
__lsx_vstx(q2, dst, stride2);
} else {
dst = dst_tmp - stride3;
p7_l_in = (v8u16)__lsx_vilvl_b(zero, p7);
p6_l_in = (v8u16)__lsx_vilvl_b(zero, p6);
p5_l_in = (v8u16)__lsx_vilvl_b(zero, p5);
p4_l_in = (v8u16)__lsx_vilvl_b(zero, p4);
p3_l_in = (v8u16)__lsx_vilvl_b(zero, p3);
p2_l_in = (v8u16)__lsx_vilvl_b(zero, p2);
p1_l_in = (v8u16)__lsx_vilvl_b(zero, p1);
p0_l_in = (v8u16)__lsx_vilvl_b(zero, p0);
q0_l_in = (v8u16)__lsx_vilvl_b(zero, q0);
tmp0_l = p7_l_in << 3;
tmp0_l -= p7_l_in;
tmp0_l += p6_l_in;
tmp0_l += q0_l_in;
tmp1_l = p6_l_in + p5_l_in;
tmp1_l += p4_l_in;
tmp1_l += p3_l_in;
tmp1_l += p2_l_in;
tmp1_l += p1_l_in;
tmp1_l += p0_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
p7_h_in = (v8u16)__lsx_vilvh_b(zero, p7);
p6_h_in = (v8u16)__lsx_vilvh_b(zero, p6);
p5_h_in = (v8u16)__lsx_vilvh_b(zero, p5);
p4_h_in = (v8u16)__lsx_vilvh_b(zero, p4);
p3_h_in = (v8u16)__lsx_vilvh_b(zero, p3);
p2_h_in = (v8u16)__lsx_vilvh_b(zero, p2);
p1_h_in = (v8u16)__lsx_vilvh_b(zero, p1);
p0_h_in = (v8u16)__lsx_vilvh_b(zero, p0);
q0_h_in = (v8u16)__lsx_vilvh_b(zero, q0);
tmp0_h = p7_h_in << 3;
tmp0_h -= p7_h_in;
tmp0_h += p6_h_in;
tmp0_h += q0_h_in;
tmp1_h = p6_h_in + p5_h_in;
tmp1_h += p4_h_in;
tmp1_h += p3_h_in;
tmp1_h += p2_h_in;
tmp1_h += p1_h_in;
tmp1_h += p0_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
p6 = __lsx_vbitsel_v(p6, out_l, flat2);
__lsx_vst(p6, dst, 0);
dst += stride;
/* p5 */
q1_l_in = (v8u16)__lsx_vilvl_b(zero, q1);
tmp0_l = p5_l_in - p6_l_in;
tmp0_l += q1_l_in;
tmp0_l -= p7_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
q1_h_in = (v8u16)__lsx_vilvh_b(zero, q1);
tmp0_h = p5_h_in - p6_h_in;
tmp0_h += q1_h_in;
tmp0_h -= p7_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
p5 = __lsx_vbitsel_v(p5, out_l, flat2);
__lsx_vst(p5, dst, 0);
dst += stride;
/* p4 */
q2_l_in = (v8u16)__lsx_vilvl_b(zero, q2);
tmp0_l = p4_l_in - p5_l_in;
tmp0_l += q2_l_in;
tmp0_l -= p7_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
q2_h_in = (v8u16)__lsx_vilvh_b(zero, q2);
tmp0_h = p4_h_in - p5_h_in;
tmp0_h += q2_h_in;
tmp0_h -= p7_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
p4 = __lsx_vbitsel_v(p4, out_l, flat2);
__lsx_vst(p4, dst, 0);
dst += stride;
/* p3 */
q3_l_in = (v8u16)__lsx_vilvl_b(zero, q3);
tmp0_l = p3_l_in - p4_l_in;
tmp0_l += q3_l_in;
tmp0_l -= p7_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
q3_h_in = (v8u16)__lsx_vilvh_b(zero, q3);
tmp0_h = p3_h_in - p4_h_in;
tmp0_h += q3_h_in;
tmp0_h -= p7_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
p3 = __lsx_vbitsel_v(p3, out_l, flat2);
__lsx_vst(p3, dst, 0);
dst += stride;
/* p2 */
q4_l_in = (v8u16)__lsx_vilvl_b(zero, q4);
filter8 = __lsx_vld(filter48, 0);
tmp0_l = p2_l_in - p3_l_in;
tmp0_l += q4_l_in;
tmp0_l -= p7_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
q4_h_in = (v8u16)__lsx_vilvh_b(zero, q4);
tmp0_h = p2_h_in - p3_h_in;
tmp0_h += q4_h_in;
tmp0_h -= p7_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
__lsx_vst(filter8, dst, 0);
dst += stride;
/* p1 */
q5_l_in = (v8u16)__lsx_vilvl_b(zero, q5);
filter8 = __lsx_vld(filter48, 16);
tmp0_l = p1_l_in - p2_l_in;
tmp0_l += q5_l_in;
tmp0_l -= p7_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
q5_h_in = (v8u16)__lsx_vilvh_b(zero, q5);
tmp0_h = p1_h_in - p2_h_in;
tmp0_h += q5_h_in;
tmp0_h -= p7_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
__lsx_vst(filter8, dst, 0);
dst += stride;
/* p0 */
q6_l_in = (v8u16)__lsx_vilvl_b(zero, q6);
filter8 = __lsx_vld(filter48, 32);
tmp0_l = p0_l_in - p1_l_in;
tmp0_l += q6_l_in;
tmp0_l -= p7_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
q6_h_in = (v8u16)__lsx_vilvh_b(zero, q6);
tmp0_h = p0_h_in - p1_h_in;
tmp0_h += q6_h_in;
tmp0_h -= p7_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
__lsx_vst(filter8, dst, 0);
dst += stride;
/* q0 */
q7_l_in = (v8u16)__lsx_vilvl_b(zero, q7);
filter8 = __lsx_vld(filter48, 48);
tmp0_l = q7_l_in - p0_l_in;
tmp0_l += q0_l_in;
tmp0_l -= p7_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
q7_h_in = (v8u16)__lsx_vilvh_b(zero, q7);
tmp0_h = q7_h_in - p0_h_in;
tmp0_h += q0_h_in;
tmp0_h -= p7_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
__lsx_vst(filter8, dst, 0);
dst += stride;
/* q1 */
filter8 = __lsx_vld(filter48, 64);
tmp0_l = q7_l_in - q0_l_in;
tmp0_l += q1_l_in;
tmp0_l -= p6_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
tmp0_h = q7_h_in - q0_h_in;
tmp0_h += q1_h_in;
tmp0_h -= p6_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
__lsx_vst(filter8, dst, 0);
dst += stride;
/* q2 */
filter8 = __lsx_vld(filter48, 80);
tmp0_l = q7_l_in - q1_l_in;
tmp0_l += q2_l_in;
tmp0_l -= p5_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
tmp0_h = q7_h_in - q1_h_in;
tmp0_h += q2_h_in;
tmp0_h -= p5_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
__lsx_vst(filter8, dst, 0);
dst += stride;
/* q3 */
tmp0_l = q7_l_in - q2_l_in;
tmp0_l += q3_l_in;
tmp0_l -= p4_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
tmp0_h = q7_h_in - q2_h_in;
tmp0_h += q3_h_in;
tmp0_h -= p4_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
q3 = __lsx_vbitsel_v(q3, out_l, flat2);
__lsx_vst(q3, dst, 0);
dst += stride;
/* q4 */
tmp0_l = q7_l_in - q3_l_in;
tmp0_l += q4_l_in;
tmp0_l -= p3_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
tmp0_h = q7_h_in - q3_h_in;
tmp0_h += q4_h_in;
tmp0_h -= p3_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
q4 = __lsx_vbitsel_v(q4, out_l, flat2);
__lsx_vst(q4, dst, 0);
dst += stride;
/* q5 */
tmp0_l = q7_l_in - q4_l_in;
tmp0_l += q5_l_in;
tmp0_l -= p2_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
tmp0_h = q7_h_in - q4_h_in;
tmp0_h += q5_h_in;
tmp0_h -= p2_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
q5 = __lsx_vbitsel_v(q5, out_l, flat2);
__lsx_vst(q5, dst, 0);
dst += stride;
/* q6 */
tmp0_l = q7_l_in - q5_l_in;
tmp0_l += q6_l_in;
tmp0_l -= p1_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
tmp0_h = q7_h_in - q5_h_in;
tmp0_h += q6_h_in;
tmp0_h -= p1_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
q6 = __lsx_vbitsel_v(q6, out_l, flat2);
__lsx_vst(q6, dst, 0);
}
}
void ff_loop_filter_v_16_16_lsx(uint8_t *dst, ptrdiff_t stride,
int32_t b_limit_ptr,
int32_t limit_ptr,
int32_t thresh_ptr)
{
uint8_t filter48[16 * 8] __attribute__ ((aligned(16)));
uint8_t early_exit = 0;
early_exit = vp9_hz_lpf_t4_and_t8_16w(dst, stride, &filter48[0],
b_limit_ptr, limit_ptr, thresh_ptr);
if (0 == early_exit) {
vp9_hz_lpf_t16_16w(dst, stride, filter48);
}
}
void ff_loop_filter_v_16_8_lsx(uint8_t *dst, ptrdiff_t stride,
int32_t b_limit_ptr,
int32_t limit_ptr,
int32_t thresh_ptr)
{
ptrdiff_t stride2 = stride << 1;
ptrdiff_t stride3 = stride2 + stride;
ptrdiff_t stride4 = stride2 << 1;
uint8_t *dst_tmp = dst - stride4;
uint8_t *dst_tmp1 = dst + stride4;
__m128i zero = __lsx_vldi(0);
__m128i flat2, mask, hev, flat, thresh, b_limit, limit;
__m128i p3, p2, p1, p0, q3, q2, q1, q0, p7, p6, p5, p4, q4, q5, q6, q7;
__m128i p2_out, p1_out, p0_out, q0_out, q1_out, q2_out;
__m128i p0_filter16, p1_filter16;
__m128i p2_filter8, p1_filter8, p0_filter8;
__m128i q0_filter8, q1_filter8, q2_filter8;
__m128i p7_l, p6_l, p5_l, p4_l, q7_l, q6_l, q5_l, q4_l;
__m128i p3_l, p2_l, p1_l, p0_l, q3_l, q2_l, q1_l, q0_l;
__m128i tmp0, tmp1, tmp2;
/* load vector elements */
DUP4_ARG2(__lsx_vldx, dst, -stride4, dst, -stride3, dst, -stride2,
dst, -stride, p3, p2, p1, p0);
q0 = __lsx_vld(dst, 0);
DUP2_ARG2(__lsx_vldx, dst, stride, dst, stride2, q1, q2);
q3 = __lsx_vldx(dst, stride3);
thresh = __lsx_vreplgr2vr_b(thresh_ptr);
b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
limit = __lsx_vreplgr2vr_b(limit_ptr);
LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
hev, mask, flat);
VP9_FLAT4(p3, p2, p0, q0, q2, q3, flat);
VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
q1_out);
flat = __lsx_vilvl_d(zero, flat);
/* if flat is zero for all pixels, then no need to calculate other filter */
if (__lsx_bz_v(flat)) {
__lsx_vstelm_d(p1_out, dst - stride2, 0, 0);
__lsx_vstelm_d(p0_out, dst - stride, 0, 0);
__lsx_vstelm_d(q0_out, dst , 0, 0);
__lsx_vstelm_d(q1_out, dst + stride, 0, 0);
} else {
/* convert 8 bit input data into 16 bit */
DUP4_ARG2(__lsx_vilvl_b, zero, p3, zero, p2, zero, p1, zero, p0,
p3_l, p2_l, p1_l, p0_l);
DUP4_ARG2(__lsx_vilvl_b, zero, q0, zero, q1, zero, q2, zero, q3,
q0_l, q1_l, q2_l, q3_l);
VP9_FILTER8(p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l,
p2_filter8, p1_filter8, p0_filter8, q0_filter8,
q1_filter8, q2_filter8);
/* convert 16 bit output data into 8 bit */
DUP4_ARG2(__lsx_vpickev_b, zero, p2_filter8, zero, p1_filter8,
zero, p0_filter8, zero, q0_filter8, p2_filter8,
p1_filter8, p0_filter8, q0_filter8);
DUP2_ARG2(__lsx_vpickev_b, zero, q1_filter8, zero, q2_filter8,
q1_filter8, q2_filter8);
/* store pixel values */
p2_out = __lsx_vbitsel_v(p2, p2_filter8, flat);
p1_out = __lsx_vbitsel_v(p1_out, p1_filter8, flat);
p0_out = __lsx_vbitsel_v(p0_out, p0_filter8, flat);
q0_out = __lsx_vbitsel_v(q0_out, q0_filter8, flat);
q1_out = __lsx_vbitsel_v(q1_out, q1_filter8, flat);
q2_out = __lsx_vbitsel_v(q2, q2_filter8, flat);
/* load 16 vector elements */
DUP4_ARG2(__lsx_vld, dst_tmp - stride4, 0, dst_tmp - stride3, 0,
dst_tmp - stride2, 0, dst_tmp - stride, 0, p7, p6, p5, p4);
DUP4_ARG2(__lsx_vld, dst_tmp1, 0, dst_tmp1 + stride, 0,
dst_tmp1 + stride2, 0, dst_tmp1 + stride3, 0, q4, q5, q6, q7);
VP9_FLAT5(p7, p6, p5, p4, p0, q0, q4, q5, q6, q7, flat, flat2);
/* if flat2 is zero for all pixels, then no need to calculate other filter */
if (__lsx_bz_v(flat2)) {
dst -= stride3;
__lsx_vstelm_d(p2_out, dst, 0, 0);
dst += stride;
__lsx_vstelm_d(p1_out, dst, 0, 0);
dst += stride;
__lsx_vstelm_d(p0_out, dst, 0, 0);
dst += stride;
__lsx_vstelm_d(q0_out, dst, 0, 0);
dst += stride;
__lsx_vstelm_d(q1_out, dst, 0, 0);
dst += stride;
__lsx_vstelm_d(q2_out, dst, 0, 0);
} else {
/* LSB(right) 8 pixel operation */
DUP4_ARG2(__lsx_vilvl_b, zero, p7, zero, p6, zero, p5, zero, p4,
p7_l, p6_l, p5_l, p4_l);
DUP4_ARG2(__lsx_vilvl_b, zero, q4, zero, q5, zero, q6, zero, q7,
q4_l, q5_l, q6_l, q7_l);
tmp0 = __lsx_vslli_h(p7_l, 3);
tmp0 = __lsx_vsub_h(tmp0, p7_l);
tmp0 = __lsx_vadd_h(tmp0, p6_l);
tmp0 = __lsx_vadd_h(tmp0, q0_l);
dst = dst_tmp - stride3;
/* calculation of p6 and p5 */
tmp1 = __lsx_vadd_h(p6_l, p5_l);
tmp1 = __lsx_vadd_h(tmp1, p4_l);
tmp1 = __lsx_vadd_h(tmp1, p3_l);
tmp1 = __lsx_vadd_h(tmp1, p2_l);
tmp1 = __lsx_vadd_h(tmp1, p1_l);
tmp1 = __lsx_vadd_h(tmp1, p0_l);
tmp1 = __lsx_vadd_h(tmp1, tmp0);
p0_filter16 = __lsx_vsrari_h(tmp1, 4);
tmp0 = __lsx_vsub_h(p5_l, p6_l);
tmp0 = __lsx_vadd_h(tmp0, q1_l);
tmp0 = __lsx_vsub_h(tmp0, p7_l);
tmp1 = __lsx_vadd_h(tmp1, tmp0);
p1_filter16 = __lsx_vsrari_h(tmp1, 4);
DUP2_ARG2(__lsx_vpickev_b, zero, p0_filter16, zero,
p1_filter16, p0_filter16, p1_filter16);
p0_filter16 = __lsx_vbitsel_v(p6, p0_filter16, flat2);
p1_filter16 = __lsx_vbitsel_v(p5, p1_filter16, flat2);
__lsx_vstelm_d(p0_filter16, dst, 0, 0);
dst += stride;
__lsx_vstelm_d(p1_filter16, dst, 0, 0);
dst += stride;
/* calculation of p4 and p3 */
tmp0 = __lsx_vsub_h(p4_l, p5_l);
tmp0 = __lsx_vadd_h(tmp0, q2_l);
tmp0 = __lsx_vsub_h(tmp0, p7_l);
tmp2 = __lsx_vsub_h(p3_l, p4_l);
tmp2 = __lsx_vadd_h(tmp2, q3_l);
tmp2 = __lsx_vsub_h(tmp2, p7_l);
tmp1 = __lsx_vadd_h(tmp1, tmp0);
p0_filter16 = __lsx_vsrari_h(tmp1, 4);
tmp1 = __lsx_vadd_h(tmp1, tmp2);
p1_filter16 = __lsx_vsrari_h(tmp1, 4);
DUP2_ARG2(__lsx_vpickev_b, zero, p0_filter16, zero,
p1_filter16, p0_filter16, p1_filter16);
p0_filter16 = __lsx_vbitsel_v(p4, p0_filter16, flat2);
p1_filter16 = __lsx_vbitsel_v(p3, p1_filter16, flat2);
__lsx_vstelm_d(p0_filter16, dst, 0, 0);
dst += stride;
__lsx_vstelm_d(p1_filter16, dst, 0, 0);
dst += stride;
/* calculation of p2 and p1 */
tmp0 = __lsx_vsub_h(p2_l, p3_l);
tmp0 = __lsx_vadd_h(tmp0, q4_l);
tmp0 = __lsx_vsub_h(tmp0, p7_l);
tmp2 = __lsx_vsub_h(p1_l, p2_l);
tmp2 = __lsx_vadd_h(tmp2, q5_l);
tmp2 = __lsx_vsub_h(tmp2, p7_l);
tmp1 = __lsx_vadd_h(tmp1, tmp0);
p0_filter16 = __lsx_vsrari_h(tmp1, 4);
tmp1 = __lsx_vadd_h(tmp1, tmp2);
p1_filter16 = __lsx_vsrari_h(tmp1, 4);
DUP2_ARG2(__lsx_vpickev_b, zero, p0_filter16, zero,
p1_filter16, p0_filter16, p1_filter16);
p0_filter16 = __lsx_vbitsel_v(p2_out, p0_filter16, flat2);
p1_filter16 = __lsx_vbitsel_v(p1_out, p1_filter16, flat2);
__lsx_vstelm_d(p0_filter16, dst, 0, 0);
dst += stride;
__lsx_vstelm_d(p1_filter16, dst, 0, 0);
dst += stride;
/* calculation of p0 and q0 */
tmp0 = __lsx_vsub_h(p0_l, p1_l);
tmp0 = __lsx_vadd_h(tmp0, q6_l);
tmp0 = __lsx_vsub_h(tmp0, p7_l);
tmp2 = __lsx_vsub_h(q7_l, p0_l);
tmp2 = __lsx_vadd_h(tmp2, q0_l);
tmp2 = __lsx_vsub_h(tmp2, p7_l);
tmp1 = __lsx_vadd_h(tmp1, tmp0);
p0_filter16 = __lsx_vsrari_h((__m128i)tmp1, 4);
tmp1 = __lsx_vadd_h(tmp1, tmp2);
p1_filter16 = __lsx_vsrari_h((__m128i)tmp1, 4);
DUP2_ARG2(__lsx_vpickev_b, zero, p0_filter16, zero,
p1_filter16, p0_filter16, p1_filter16);
p0_filter16 = __lsx_vbitsel_v(p0_out, p0_filter16, flat2);
p1_filter16 = __lsx_vbitsel_v(q0_out, p1_filter16, flat2);
__lsx_vstelm_d(p0_filter16, dst, 0, 0);
dst += stride;
__lsx_vstelm_d(p1_filter16, dst, 0, 0);
dst += stride;
/* calculation of q1 and q2 */
tmp0 = __lsx_vsub_h(q7_l, q0_l);
tmp0 = __lsx_vadd_h(tmp0, q1_l);
tmp0 = __lsx_vsub_h(tmp0, p6_l);
tmp2 = __lsx_vsub_h(q7_l, q1_l);
tmp2 = __lsx_vadd_h(tmp2, q2_l);
tmp2 = __lsx_vsub_h(tmp2, p5_l);
tmp1 = __lsx_vadd_h(tmp1, tmp0);
p0_filter16 = __lsx_vsrari_h(tmp1, 4);
tmp1 = __lsx_vadd_h(tmp1, tmp2);
p1_filter16 = __lsx_vsrari_h(tmp1, 4);
DUP2_ARG2(__lsx_vpickev_b, zero, p0_filter16, zero,
p1_filter16, p0_filter16, p1_filter16);
p0_filter16 = __lsx_vbitsel_v(q1_out, p0_filter16, flat2);
p1_filter16 = __lsx_vbitsel_v(q2_out, p1_filter16, flat2);
__lsx_vstelm_d(p0_filter16, dst, 0, 0);
dst += stride;
__lsx_vstelm_d(p1_filter16, dst, 0, 0);
dst += stride;
/* calculation of q3 and q4 */
tmp0 = __lsx_vsub_h(q7_l, q2_l);
tmp0 = __lsx_vadd_h(tmp0, q3_l);
tmp0 = __lsx_vsub_h(tmp0, p4_l);
tmp2 = __lsx_vsub_h(q7_l, q3_l);
tmp2 = __lsx_vadd_h(tmp2, q4_l);
tmp2 = __lsx_vsub_h(tmp2, p3_l);
tmp1 = __lsx_vadd_h(tmp1, tmp0);
p0_filter16 = __lsx_vsrari_h(tmp1, 4);
tmp1 = __lsx_vadd_h(tmp1, tmp2);
p1_filter16 = __lsx_vsrari_h(tmp1, 4);
DUP2_ARG2(__lsx_vpickev_b, zero, p0_filter16, zero,
p1_filter16, p0_filter16, p1_filter16);
p0_filter16 = __lsx_vbitsel_v(q3, p0_filter16, flat2);
p1_filter16 = __lsx_vbitsel_v(q4, p1_filter16, flat2);
__lsx_vstelm_d(p0_filter16, dst, 0, 0);
dst += stride;
__lsx_vstelm_d(p1_filter16, dst, 0, 0);
dst += stride;
/* calculation of q5 and q6 */
tmp0 = __lsx_vsub_h(q7_l, q4_l);
tmp0 = __lsx_vadd_h(tmp0, q5_l);
tmp0 = __lsx_vsub_h(tmp0, p2_l);
tmp2 = __lsx_vsub_h(q7_l, q5_l);
tmp2 = __lsx_vadd_h(tmp2, q6_l);
tmp2 = __lsx_vsub_h(tmp2, p1_l);
tmp1 = __lsx_vadd_h(tmp1, tmp0);
p0_filter16 = __lsx_vsrari_h(tmp1, 4);
tmp1 = __lsx_vadd_h(tmp1, tmp2);
p1_filter16 = __lsx_vsrari_h(tmp1, 4);
DUP2_ARG2(__lsx_vpickev_b, zero, p0_filter16, zero,
p1_filter16, p0_filter16, p1_filter16);
p0_filter16 = __lsx_vbitsel_v(q5, p0_filter16, flat2);
p1_filter16 = __lsx_vbitsel_v(q6, p1_filter16, flat2);
__lsx_vstelm_d(p0_filter16, dst, 0, 0);
dst += stride;
__lsx_vstelm_d(p1_filter16, dst, 0, 0);
}
}
}
void ff_loop_filter_h_4_8_lsx(uint8_t *dst, ptrdiff_t stride,
int32_t b_limit_ptr,
int32_t limit_ptr,
int32_t thresh_ptr)
{
ptrdiff_t stride2 = stride << 1;
ptrdiff_t stride3 = stride2 + stride;
ptrdiff_t stride4 = stride2 << 1;
uint8_t *dst_tmp1 = dst - 4;
uint8_t *dst_tmp2 = dst_tmp1 + stride4;
__m128i mask, hev, flat, limit, thresh, b_limit;
__m128i p3, p2, p1, p0, q3, q2, q1, q0;
__m128i vec0, vec1, vec2, vec3;
p3 = __lsx_vld(dst_tmp1, 0);
DUP2_ARG2(__lsx_vldx, dst_tmp1, stride, dst_tmp1, stride2, p2, p1);
p0 = __lsx_vldx(dst_tmp1, stride3);
q0 = __lsx_vld(dst_tmp2, 0);
DUP2_ARG2(__lsx_vldx, dst_tmp2, stride, dst_tmp2, stride2, q1, q2);
q3 = __lsx_vldx(dst_tmp2, stride3);
thresh = __lsx_vreplgr2vr_b(thresh_ptr);
b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
limit = __lsx_vreplgr2vr_b(limit_ptr);
LSX_TRANSPOSE8x8_B(p3, p2, p1, p0, q0, q1, q2, q3,
p3, p2, p1, p0, q0, q1, q2, q3);
LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
hev, mask, flat);
VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1, p0, q0, q1);
DUP2_ARG2(__lsx_vilvl_b, p0, p1, q1, q0, vec0, vec1);
vec2 = __lsx_vilvl_h(vec1, vec0);
vec3 = __lsx_vilvh_h(vec1, vec0);
dst -= 2;
__lsx_vstelm_w(vec2, dst, 0, 0);
__lsx_vstelm_w(vec2, dst + stride, 0, 1);
__lsx_vstelm_w(vec2, dst + stride2, 0, 2);
__lsx_vstelm_w(vec2, dst + stride3, 0, 3);
dst += stride4;
__lsx_vstelm_w(vec3, dst, 0, 0);
__lsx_vstelm_w(vec3, dst + stride, 0, 1);
__lsx_vstelm_w(vec3, dst + stride2, 0, 2);
__lsx_vstelm_w(vec3, dst + stride3, 0, 3);
}
void ff_loop_filter_h_44_16_lsx(uint8_t *dst, ptrdiff_t stride,
int32_t b_limit_ptr,
int32_t limit_ptr,
int32_t thresh_ptr)
{
ptrdiff_t stride2 = stride << 1;
ptrdiff_t stride3 = stride2 + stride;
ptrdiff_t stride4 = stride2 << 1;
uint8_t *dst_tmp = dst - 4;
__m128i mask, hev, flat;
__m128i thresh0, b_limit0, limit0, thresh1, b_limit1, limit1;
__m128i p3, p2, p1, p0, q3, q2, q1, q0;
__m128i row0, row1, row2, row3, row4, row5, row6, row7;
__m128i row8, row9, row10, row11, row12, row13, row14, row15;
__m128i tmp0, tmp1, tmp2, tmp3, tmp4, tmp5;
row0 = __lsx_vld(dst_tmp, 0);
DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, row1, row2);
row3 = __lsx_vldx(dst_tmp, stride3);
dst_tmp += stride4;
row4 = __lsx_vld(dst_tmp, 0);
DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, row5, row6);
row7 = __lsx_vldx(dst_tmp, stride3);
dst_tmp += stride4;
row8 = __lsx_vld(dst_tmp, 0);
DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, row9, row10);
row11 = __lsx_vldx(dst_tmp, stride3);
dst_tmp += stride4;
row12 = __lsx_vld(dst_tmp, 0);
DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, row13, row14);
row15 = __lsx_vldx(dst_tmp, stride3);
LSX_TRANSPOSE16x8_B(row0, row1, row2, row3, row4, row5, row6, row7,
row8, row9, row10, row11, row12, row13, row14, row15,
p3, p2, p1, p0, q0, q1, q2, q3);
thresh0 = __lsx_vreplgr2vr_b(thresh_ptr);
thresh1 = __lsx_vreplgr2vr_b(thresh_ptr >> 8);
thresh0 = __lsx_vilvl_d(thresh1, thresh0);
b_limit0 = __lsx_vreplgr2vr_b(b_limit_ptr);
b_limit1 = __lsx_vreplgr2vr_b(b_limit_ptr >> 8);
b_limit0 = __lsx_vilvl_d(b_limit1, b_limit0);
limit0 = __lsx_vreplgr2vr_b(limit_ptr);
limit1 = __lsx_vreplgr2vr_b(limit_ptr >> 8);
limit0 = __lsx_vilvl_d(limit1, limit0);
LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit0, b_limit0, thresh0,
hev, mask, flat);
VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1, p0, q0, q1);
DUP2_ARG2(__lsx_vilvl_b, p0, p1, q1, q0, tmp0, tmp1);
tmp2 = __lsx_vilvl_h(tmp1, tmp0);
tmp3 = __lsx_vilvh_h(tmp1, tmp0);
DUP2_ARG2(__lsx_vilvh_b, p0, p1, q1, q0, tmp0, tmp1);
tmp4 = __lsx_vilvl_h(tmp1, tmp0);
tmp5 = __lsx_vilvh_h(tmp1, tmp0);
dst -= 2;
__lsx_vstelm_w(tmp2, dst, 0, 0);
__lsx_vstelm_w(tmp2, dst + stride, 0, 1);
__lsx_vstelm_w(tmp2, dst + stride2, 0, 2);
__lsx_vstelm_w(tmp2, dst + stride3, 0, 3);
dst += stride4;
__lsx_vstelm_w(tmp3, dst, 0, 0);
__lsx_vstelm_w(tmp3, dst + stride, 0, 1);
__lsx_vstelm_w(tmp3, dst + stride2, 0, 2);
__lsx_vstelm_w(tmp3, dst + stride3, 0, 3);
dst += stride4;
__lsx_vstelm_w(tmp4, dst, 0, 0);
__lsx_vstelm_w(tmp4, dst + stride, 0, 1);
__lsx_vstelm_w(tmp4, dst + stride2, 0, 2);
__lsx_vstelm_w(tmp4, dst + stride3, 0, 3);
dst += stride4;
__lsx_vstelm_w(tmp5, dst, 0, 0);
__lsx_vstelm_w(tmp5, dst + stride, 0, 1);
__lsx_vstelm_w(tmp5, dst + stride2, 0, 2);
__lsx_vstelm_w(tmp5, dst + stride3, 0, 3);
}
void ff_loop_filter_h_8_8_lsx(uint8_t *dst, ptrdiff_t stride,
int32_t b_limit_ptr,
int32_t limit_ptr,
int32_t thresh_ptr)
{
ptrdiff_t stride2 = stride << 1;
ptrdiff_t stride3 = stride2 + stride;
ptrdiff_t stride4 = stride2 << 1;
uint8_t *dst_tmp = dst - 4;
__m128i p3, p2, p1, p0, q3, q2, q1, q0;
__m128i p1_out, p0_out, q0_out, q1_out;
__m128i flat, mask, hev, thresh, b_limit, limit;
__m128i p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l;
__m128i p2_filt8_l, p1_filt8_l, p0_filt8_l;
__m128i q0_filt8_l, q1_filt8_l, q2_filt8_l;
__m128i vec0, vec1, vec2, vec3, vec4;
__m128i zero = __lsx_vldi(0);
/* load vector elements */
p3 = __lsx_vld(dst_tmp, 0);
DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, p2, p1);
p0 = __lsx_vldx(dst_tmp, stride3);
dst_tmp += stride4;
q0 = __lsx_vld(dst_tmp, 0);
DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, q1, q2);
q3 = __lsx_vldx(dst_tmp, stride3);
LSX_TRANSPOSE8x8_B(p3, p2, p1, p0, q0, q1, q2, q3,
p3, p2, p1, p0, q0, q1, q2, q3);
thresh = __lsx_vreplgr2vr_b(thresh_ptr);
b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
limit = __lsx_vreplgr2vr_b(limit_ptr);
/* mask and hev */
LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
hev, mask, flat);
/* flat4 */
VP9_FLAT4(p3, p2, p0, q0, q2, q3, flat);
/* filter4 */
VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
q1_out);
flat = __lsx_vilvl_d(zero, flat);
/* if flat is zero for all pixels, then no need to calculate other filter */
if (__lsx_bz_v(flat)) {
/* Store 4 pixels p1-_q1 */
DUP2_ARG2(__lsx_vilvl_b, p0_out, p1_out, q1_out, q0_out, vec0, vec1);
vec2 = __lsx_vilvl_h(vec1, vec0);
vec3 = __lsx_vilvh_h(vec1, vec0);
dst -= 2;
__lsx_vstelm_w(vec2, dst, 0, 0);
__lsx_vstelm_w(vec2, dst + stride, 0, 1);
__lsx_vstelm_w(vec2, dst + stride2, 0, 2);
__lsx_vstelm_w(vec2, dst + stride3, 0, 3);
dst += stride4;
__lsx_vstelm_w(vec3, dst, 0, 0);
__lsx_vstelm_w(vec3, dst + stride, 0, 1);
__lsx_vstelm_w(vec3, dst + stride2, 0, 2);
__lsx_vstelm_w(vec3, dst + stride3, 0, 3);
} else {
DUP4_ARG2(__lsx_vilvl_b, zero, p3, zero, p2, zero, p1, zero, p0,
p3_l, p2_l, p1_l, p0_l);
DUP4_ARG2(__lsx_vilvl_b, zero, q0, zero, q1, zero, q2, zero, q3,
q0_l, q1_l, q2_l, q3_l);
VP9_FILTER8(p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l, p2_filt8_l,
p1_filt8_l, p0_filt8_l, q0_filt8_l, q1_filt8_l, q2_filt8_l);
/* convert 16 bit output data into 8 bit */
DUP4_ARG2(__lsx_vpickev_b, p2_filt8_l, p2_filt8_l, p1_filt8_l,
p1_filt8_l, p0_filt8_l, p0_filt8_l, q0_filt8_l,
q0_filt8_l, p2_filt8_l, p1_filt8_l, p0_filt8_l,
q0_filt8_l);
DUP2_ARG2(__lsx_vpickev_b, q1_filt8_l, q1_filt8_l, q2_filt8_l,
q2_filt8_l, q1_filt8_l, q2_filt8_l);
/* store pixel values */
p2 = __lsx_vbitsel_v(p2, p2_filt8_l, flat);
p1 = __lsx_vbitsel_v(p1_out, p1_filt8_l, flat);
p0 = __lsx_vbitsel_v(p0_out, p0_filt8_l, flat);
q0 = __lsx_vbitsel_v(q0_out, q0_filt8_l, flat);
q1 = __lsx_vbitsel_v(q1_out, q1_filt8_l, flat);
q2 = __lsx_vbitsel_v(q2, q2_filt8_l, flat);
/* Store 6 pixels p2-_q2 */
DUP2_ARG2(__lsx_vilvl_b, p1, p2, q0, p0, vec0, vec1);
vec2 = __lsx_vilvl_h(vec1, vec0);
vec3 = __lsx_vilvh_h(vec1, vec0);
vec4 = __lsx_vilvl_b(q2, q1);
dst -= 3;
__lsx_vstelm_w(vec2, dst, 0, 0);
__lsx_vstelm_h(vec4, dst, 4, 0);
dst += stride;
__lsx_vstelm_w(vec2, dst, 0, 1);
__lsx_vstelm_h(vec4, dst, 4, 1);
dst += stride;
__lsx_vstelm_w(vec2, dst, 0, 2);
__lsx_vstelm_h(vec4, dst, 4, 2);
dst += stride;
__lsx_vstelm_w(vec2, dst, 0, 3);
__lsx_vstelm_h(vec4, dst, 4, 3);
dst += stride;
__lsx_vstelm_w(vec3, dst, 0, 0);
__lsx_vstelm_h(vec4, dst, 4, 4);
dst += stride;
__lsx_vstelm_w(vec3, dst, 0, 1);
__lsx_vstelm_h(vec4, dst, 4, 5);
dst += stride;
__lsx_vstelm_w(vec3, dst, 0, 2);
__lsx_vstelm_h(vec4, dst, 4, 6);
dst += stride;
__lsx_vstelm_w(vec3, dst, 0, 3);
__lsx_vstelm_h(vec4, dst, 4, 7);
}
}
void ff_loop_filter_h_88_16_lsx(uint8_t *dst, ptrdiff_t stride,
int32_t b_limit_ptr,
int32_t limit_ptr,
int32_t thresh_ptr)
{
ptrdiff_t stride2 = stride << 1;
ptrdiff_t stride3 = stride2 + stride;
ptrdiff_t stride4 = stride2 << 1;
uint8_t *dst_tmp = dst - 4;
__m128i p3, p2, p1, p0, q3, q2, q1, q0;
__m128i p1_out, p0_out, q0_out, q1_out;
__m128i flat, mask, hev, thresh, b_limit, limit;
__m128i row4, row5, row6, row7, row12, row13, row14, row15;
__m128i p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l;
__m128i p3_h, p2_h, p1_h, p0_h, q0_h, q1_h, q2_h, q3_h;
__m128i p2_filt8_l, p1_filt8_l, p0_filt8_l;
__m128i q0_filt8_l, q1_filt8_l, q2_filt8_l;
__m128i p2_filt8_h, p1_filt8_h, p0_filt8_h;
__m128i q0_filt8_h, q1_filt8_h, q2_filt8_h;
__m128i vec0, vec1, vec2, vec3, vec4, vec5, vec6, vec7;
__m128i zero = __lsx_vldi(0);
p0 = __lsx_vld(dst_tmp, 0);
DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, p1, p2);
p3 = __lsx_vldx(dst_tmp, stride3);
dst_tmp += stride4;
row4 = __lsx_vld(dst_tmp, 0);
DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, row5, row6);
row7 = __lsx_vldx(dst_tmp, stride3);
dst_tmp += stride4;
q3 = __lsx_vld(dst_tmp, 0);
DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, q2, q1);
q0 = __lsx_vldx(dst_tmp, stride3);
dst_tmp += stride4;
row12 = __lsx_vld(dst_tmp, 0);
DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, row13, row14);
row15 = __lsx_vldx(dst_tmp, stride3);
/* transpose 16x8 matrix into 8x16 */
LSX_TRANSPOSE16x8_B(p0, p1, p2, p3, row4, row5, row6, row7,
q3, q2, q1, q0, row12, row13, row14, row15,
p3, p2, p1, p0, q0, q1, q2, q3);
thresh = __lsx_vreplgr2vr_b(thresh_ptr);
vec0 = __lsx_vreplgr2vr_b(thresh_ptr >> 8);
thresh = __lsx_vilvl_d(vec0, thresh);
b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
vec0 = __lsx_vreplgr2vr_b(b_limit_ptr >> 8);
b_limit = __lsx_vilvl_d(vec0, b_limit);
limit = __lsx_vreplgr2vr_b(limit_ptr);
vec0 = __lsx_vreplgr2vr_b(limit_ptr >> 8);
limit = __lsx_vilvl_d(vec0, limit);
/* mask and hev */
LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
hev, mask, flat);
/* flat4 */
VP9_FLAT4(p3, p2, p0, q0, q2, q3, flat);
/* filter4 */
VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
q1_out);
/* if flat is zero for all pixels, then no need to calculate other filter */
if (__lsx_bz_v(flat)) {
DUP2_ARG2(__lsx_vilvl_b, p0_out, p1_out, q1_out, q0_out, vec0, vec1);
vec2 = __lsx_vilvl_h(vec1, vec0);
vec3 = __lsx_vilvh_h(vec1, vec0);
DUP2_ARG2(__lsx_vilvh_b, p0_out, p1_out, q1_out, q0_out, vec0, vec1);
vec4 = __lsx_vilvl_h(vec1, vec0);
vec5 = __lsx_vilvh_h(vec1, vec0);
dst -= 2;
__lsx_vstelm_w(vec2, dst, 0, 0);
__lsx_vstelm_w(vec2, dst + stride, 0, 1);
__lsx_vstelm_w(vec2, dst + stride2, 0, 2);
__lsx_vstelm_w(vec2, dst + stride3, 0, 3);
dst += stride4;
__lsx_vstelm_w(vec3, dst, 0, 0);
__lsx_vstelm_w(vec3, dst + stride, 0, 1);
__lsx_vstelm_w(vec3, dst + stride2, 0, 2);
__lsx_vstelm_w(vec3, dst + stride3, 0, 3);
dst += stride4;
__lsx_vstelm_w(vec4, dst, 0, 0);
__lsx_vstelm_w(vec4, dst + stride, 0, 1);
__lsx_vstelm_w(vec4, dst + stride2, 0, 2);
__lsx_vstelm_w(vec4, dst + stride3, 0, 3);
dst += stride4;
__lsx_vstelm_w(vec5, dst, 0, 0);
__lsx_vstelm_w(vec5, dst + stride, 0, 1);
__lsx_vstelm_w(vec5, dst + stride2, 0, 2);
__lsx_vstelm_w(vec5, dst + stride3, 0, 3);
} else {
DUP4_ARG2(__lsx_vilvl_b, zero, p3, zero, p2, zero, p1, zero, p0,
p3_l, p2_l, p1_l, p0_l);
DUP4_ARG2(__lsx_vilvl_b, zero, q0, zero, q1, zero, q2, zero, q3,
q0_l, q1_l, q2_l, q3_l);
VP9_FILTER8(p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l, p2_filt8_l,
p1_filt8_l, p0_filt8_l, q0_filt8_l, q1_filt8_l, q2_filt8_l);
DUP4_ARG2(__lsx_vilvh_b, zero, p3, zero, p2, zero, p1, zero, p0,
p3_h, p2_h, p1_h, p0_h);
DUP4_ARG2(__lsx_vilvh_b, zero, q0, zero, q1, zero, q2, zero, q3,
q0_h, q1_h, q2_h, q3_h);
/* filter8 */
VP9_FILTER8(p3_h, p2_h, p1_h, p0_h, q0_h, q1_h, q2_h, q3_h, p2_filt8_h,
p1_filt8_h, p0_filt8_h, q0_filt8_h, q1_filt8_h, q2_filt8_h);
/* convert 16 bit output data into 8 bit */
DUP4_ARG2(__lsx_vpickev_b, p2_filt8_h, p2_filt8_l, p1_filt8_h,
p1_filt8_l, p0_filt8_h, p0_filt8_l, q0_filt8_h, q0_filt8_l,
p2_filt8_l, p1_filt8_l, p0_filt8_l, q0_filt8_l);
DUP2_ARG2(__lsx_vpickev_b, q1_filt8_h, q1_filt8_l, q2_filt8_h,
q2_filt8_l, q1_filt8_l, q2_filt8_l);
/* store pixel values */
p2 = __lsx_vbitsel_v(p2, p2_filt8_l, flat);
p1 = __lsx_vbitsel_v(p1_out, p1_filt8_l, flat);
p0 = __lsx_vbitsel_v(p0_out, p0_filt8_l, flat);
q0 = __lsx_vbitsel_v(q0_out, q0_filt8_l, flat);
q1 = __lsx_vbitsel_v(q1_out, q1_filt8_l, flat);
q2 = __lsx_vbitsel_v(q2, q2_filt8_l, flat);
DUP2_ARG2(__lsx_vilvl_b, p1, p2, q0, p0, vec0, vec1);
vec3 = __lsx_vilvl_h(vec1, vec0);
vec4 = __lsx_vilvh_h(vec1, vec0);
DUP2_ARG2(__lsx_vilvh_b, p1, p2, q0, p0, vec0, vec1);
vec6 = __lsx_vilvl_h(vec1, vec0);
vec7 = __lsx_vilvh_h(vec1, vec0);
vec2 = __lsx_vilvl_b(q2, q1);
vec5 = __lsx_vilvh_b(q2, q1);
dst -= 3;
__lsx_vstelm_w(vec3, dst, 0, 0);
__lsx_vstelm_h(vec2, dst, 4, 0);
dst += stride;
__lsx_vstelm_w(vec3, dst, 0, 1);
__lsx_vstelm_h(vec2, dst, 4, 1);
dst += stride;
__lsx_vstelm_w(vec3, dst, 0, 2);
__lsx_vstelm_h(vec2, dst, 4, 2);
dst += stride;
__lsx_vstelm_w(vec3, dst, 0, 3);
__lsx_vstelm_h(vec2, dst, 4, 3);
dst += stride;
__lsx_vstelm_w(vec4, dst, 0, 0);
__lsx_vstelm_h(vec2, dst, 4, 4);
dst += stride;
__lsx_vstelm_w(vec4, dst, 0, 1);
__lsx_vstelm_h(vec2, dst, 4, 5);
dst += stride;
__lsx_vstelm_w(vec4, dst, 0, 2);
__lsx_vstelm_h(vec2, dst, 4, 6);
dst += stride;
__lsx_vstelm_w(vec4, dst, 0, 3);
__lsx_vstelm_h(vec2, dst, 4, 7);
dst += stride;
__lsx_vstelm_w(vec6, dst, 0, 0);
__lsx_vstelm_h(vec5, dst, 4, 0);
dst += stride;
__lsx_vstelm_w(vec6, dst, 0, 1);
__lsx_vstelm_h(vec5, dst, 4, 1);
dst += stride;
__lsx_vstelm_w(vec6, dst, 0, 2);
__lsx_vstelm_h(vec5, dst, 4, 2);
dst += stride;
__lsx_vstelm_w(vec6, dst, 0, 3);
__lsx_vstelm_h(vec5, dst, 4, 3);
dst += stride;
__lsx_vstelm_w(vec7, dst, 0, 0);
__lsx_vstelm_h(vec5, dst, 4, 4);
dst += stride;
__lsx_vstelm_w(vec7, dst, 0, 1);
__lsx_vstelm_h(vec5, dst, 4, 5);
dst += stride;
__lsx_vstelm_w(vec7, dst, 0, 2);
__lsx_vstelm_h(vec5, dst, 4, 6);
dst += stride;
__lsx_vstelm_w(vec7, dst, 0, 3);
__lsx_vstelm_h(vec5, dst, 4, 7);
}
}
void ff_loop_filter_h_84_16_lsx(uint8_t *dst, ptrdiff_t stride,
int32_t b_limit_ptr,
int32_t limit_ptr,
int32_t thresh_ptr)
{
ptrdiff_t stride2 = stride << 1;
ptrdiff_t stride3 = stride2 + stride;
ptrdiff_t stride4 = stride2 << 1;
uint8_t *dst_tmp = dst - 4;
__m128i p3, p2, p1, p0, q3, q2, q1, q0;
__m128i p1_out, p0_out, q0_out, q1_out;
__m128i flat, mask, hev, thresh, b_limit, limit;
__m128i row4, row5, row6, row7, row12, row13, row14, row15;
__m128i p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l;
__m128i p2_filt8_l, p1_filt8_l, p0_filt8_l;
__m128i q0_filt8_l, q1_filt8_l, q2_filt8_l;
__m128i vec0, vec1, vec2, vec3, vec4, vec5, vec6, vec7;
__m128i zero = __lsx_vldi(0);
p0 = __lsx_vld(dst_tmp, 0);
DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, p1, p2);
p3 = __lsx_vldx(dst_tmp, stride3);
dst_tmp += stride4;
row4 = __lsx_vld(dst_tmp, 0);
DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, row5, row6);
row7 = __lsx_vldx(dst_tmp, stride3);
dst_tmp += stride4;
q3 = __lsx_vld(dst_tmp, 0);
DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, q2, q1);
q0 = __lsx_vldx(dst_tmp, stride3);
dst_tmp += stride4;
row12 = __lsx_vld(dst_tmp, 0);
DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, row13, row14);
row15 = __lsx_vldx(dst_tmp, stride3);
/* transpose 16x8 matrix into 8x16 */
LSX_TRANSPOSE16x8_B(p0, p1, p2, p3, row4, row5, row6, row7,
q3, q2, q1, q0, row12, row13, row14, row15,
p3, p2, p1, p0, q0, q1, q2, q3);
thresh = __lsx_vreplgr2vr_b(thresh_ptr);
vec0 = __lsx_vreplgr2vr_b(thresh_ptr >> 8);
thresh = __lsx_vilvl_d(vec0, thresh);
b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
vec0 = __lsx_vreplgr2vr_b(b_limit_ptr >> 8);
b_limit = __lsx_vilvl_d(vec0, b_limit);
limit = __lsx_vreplgr2vr_b(limit_ptr);
vec0 = __lsx_vreplgr2vr_b(limit_ptr >> 8);
limit = __lsx_vilvl_d(vec0, limit);
/* mask and hev */
LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
hev, mask, flat);
/* flat4 */
VP9_FLAT4(p3, p2, p0, q0, q2, q3, flat);
/* filter4 */
VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
q1_out);
flat = __lsx_vilvl_d(zero, flat);
/* if flat is zero for all pixels, then no need to calculate other filter */
if (__lsx_bz_v(flat)) {
DUP2_ARG2(__lsx_vilvl_b, p0_out, p1_out, q1_out, q0_out, vec0, vec1);
vec2 = __lsx_vilvl_h(vec1, vec0);
vec3 = __lsx_vilvh_h(vec1, vec0);
DUP2_ARG2(__lsx_vilvh_b, p0_out, p1_out, q1_out, q0_out, vec0, vec1);
vec4 = __lsx_vilvl_h(vec1, vec0);
vec5 = __lsx_vilvh_h(vec1, vec0);
dst -= 2;
__lsx_vstelm_w(vec2, dst, 0, 0);
__lsx_vstelm_w(vec2, dst + stride, 0, 1);
__lsx_vstelm_w(vec2, dst + stride2, 0, 2);
__lsx_vstelm_w(vec2, dst + stride3, 0, 3);
dst += stride4;
__lsx_vstelm_w(vec3, dst, 0, 0);
__lsx_vstelm_w(vec3, dst + stride, 0, 1);
__lsx_vstelm_w(vec3, dst + stride2, 0, 2);
__lsx_vstelm_w(vec3, dst + stride3, 0, 3);
dst += stride4;
__lsx_vstelm_w(vec4, dst, 0, 0);
__lsx_vstelm_w(vec4, dst + stride, 0, 1);
__lsx_vstelm_w(vec4, dst + stride2, 0, 2);
__lsx_vstelm_w(vec4, dst + stride3, 0, 3);
dst += stride4;
__lsx_vstelm_w(vec5, dst, 0, 0);
__lsx_vstelm_w(vec5, dst + stride, 0, 1);
__lsx_vstelm_w(vec5, dst + stride2, 0, 2);
__lsx_vstelm_w(vec5, dst + stride3, 0, 3);
} else {
DUP4_ARG2(__lsx_vilvl_b, zero, p3, zero, p2, zero, p1, zero, p0,
p3_l, p2_l, p1_l, p0_l);
DUP4_ARG2(__lsx_vilvl_b, zero, q0, zero, q1, zero, q2, zero, q3,
q0_l, q1_l, q2_l, q3_l);
VP9_FILTER8(p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l, p2_filt8_l,
p1_filt8_l, p0_filt8_l, q0_filt8_l, q1_filt8_l, q2_filt8_l);
/* convert 16 bit output data into 8 bit */
DUP4_ARG2(__lsx_vpickev_b, p2_filt8_l, p2_filt8_l, p1_filt8_l, p1_filt8_l,
p0_filt8_l, p0_filt8_l, q0_filt8_l, q0_filt8_l, p2_filt8_l,
p1_filt8_l, p0_filt8_l, q0_filt8_l);
DUP2_ARG2(__lsx_vpickev_b, q1_filt8_l, q1_filt8_l, q2_filt8_l, q2_filt8_l,
q1_filt8_l, q2_filt8_l);
/* store pixel values */
p2 = __lsx_vbitsel_v(p2, p2_filt8_l, flat);
p1 = __lsx_vbitsel_v(p1_out, p1_filt8_l, flat);
p0 = __lsx_vbitsel_v(p0_out, p0_filt8_l, flat);
q0 = __lsx_vbitsel_v(q0_out, q0_filt8_l, flat);
q1 = __lsx_vbitsel_v(q1_out, q1_filt8_l, flat);
q2 = __lsx_vbitsel_v(q2, q2_filt8_l, flat);
DUP2_ARG2(__lsx_vilvl_b, p1, p2, q0, p0, vec0, vec1);
vec3 = __lsx_vilvl_h(vec1, vec0);
vec4 = __lsx_vilvh_h(vec1, vec0);
DUP2_ARG2(__lsx_vilvh_b, p1, p2, q0, p0, vec0, vec1);
vec6 = __lsx_vilvl_h(vec1, vec0);
vec7 = __lsx_vilvh_h(vec1, vec0);
vec2 = __lsx_vilvl_b(q2, q1);
vec5 = __lsx_vilvh_b(q2, q1);
dst -= 3;
__lsx_vstelm_w(vec3, dst, 0, 0);
__lsx_vstelm_h(vec2, dst, 4, 0);
dst += stride;
__lsx_vstelm_w(vec3, dst, 0, 1);
__lsx_vstelm_h(vec2, dst, 4, 1);
dst += stride;
__lsx_vstelm_w(vec3, dst, 0, 2);
__lsx_vstelm_h(vec2, dst, 4, 2);
dst += stride;
__lsx_vstelm_w(vec3, dst, 0, 3);
__lsx_vstelm_h(vec2, dst, 4, 3);
dst += stride;
__lsx_vstelm_w(vec4, dst, 0, 0);
__lsx_vstelm_h(vec2, dst, 4, 4);
dst += stride;
__lsx_vstelm_w(vec4, dst, 0, 1);
__lsx_vstelm_h(vec2, dst, 4, 5);
dst += stride;
__lsx_vstelm_w(vec4, dst, 0, 2);
__lsx_vstelm_h(vec2, dst, 4, 6);
dst += stride;
__lsx_vstelm_w(vec4, dst, 0, 3);
__lsx_vstelm_h(vec2, dst, 4, 7);
dst += stride;
__lsx_vstelm_w(vec6, dst, 0, 0);
__lsx_vstelm_h(vec5, dst, 4, 0);
dst += stride;
__lsx_vstelm_w(vec6, dst, 0, 1);
__lsx_vstelm_h(vec5, dst, 4, 1);
dst += stride;
__lsx_vstelm_w(vec6, dst, 0, 2);
__lsx_vstelm_h(vec5, dst, 4, 2);
dst += stride;
__lsx_vstelm_w(vec6, dst, 0, 3);
__lsx_vstelm_h(vec5, dst, 4, 3);
dst += stride;
__lsx_vstelm_w(vec7, dst, 0, 0);
__lsx_vstelm_h(vec5, dst, 4, 4);
dst += stride;
__lsx_vstelm_w(vec7, dst, 0, 1);
__lsx_vstelm_h(vec5, dst, 4, 5);
dst += stride;
__lsx_vstelm_w(vec7, dst, 0, 2);
__lsx_vstelm_h(vec5, dst, 4, 6);
dst += stride;
__lsx_vstelm_w(vec7, dst, 0, 3);
__lsx_vstelm_h(vec5, dst, 4, 7);
}
}
void ff_loop_filter_h_48_16_lsx(uint8_t *dst, ptrdiff_t stride,
int32_t b_limit_ptr,
int32_t limit_ptr,
int32_t thresh_ptr)
{
ptrdiff_t stride2 = stride << 1;
ptrdiff_t stride3 = stride2 + stride;
ptrdiff_t stride4 = stride2 << 1;
uint8_t *dst_tmp = dst - 4;
__m128i p3, p2, p1, p0, q3, q2, q1, q0;
__m128i p1_out, p0_out, q0_out, q1_out;
__m128i flat, mask, hev, thresh, b_limit, limit;
__m128i row4, row5, row6, row7, row12, row13, row14, row15;
__m128i p3_h, p2_h, p1_h, p0_h, q0_h, q1_h, q2_h, q3_h;
__m128i p2_filt8_h, p1_filt8_h, p0_filt8_h;
__m128i q0_filt8_h, q1_filt8_h, q2_filt8_h;
__m128i vec0, vec1, vec2, vec3, vec4, vec5, vec6, vec7;
__m128i zero = __lsx_vldi(0);
p0 = __lsx_vld(dst_tmp, 0);
DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, p1, p2);
p3 = __lsx_vldx(dst_tmp, stride3);
dst_tmp += stride4;
row4 = __lsx_vld(dst_tmp, 0);
DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, row5, row6);
row7 = __lsx_vldx(dst_tmp, stride3);
dst_tmp += stride4;
q3 = __lsx_vld(dst_tmp, 0);
DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, q2, q1);
q0 = __lsx_vldx(dst_tmp, stride3);
dst_tmp += stride4;
row12 = __lsx_vld(dst_tmp, 0);
DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, row13, row14);
row15 = __lsx_vldx(dst_tmp, stride3);
/* transpose 16x8 matrix into 8x16 */
LSX_TRANSPOSE16x8_B(p0, p1, p2, p3, row4, row5, row6, row7,
q3, q2, q1, q0, row12, row13, row14, row15,
p3, p2, p1, p0, q0, q1, q2, q3);
thresh = __lsx_vreplgr2vr_b(thresh_ptr);
vec0 = __lsx_vreplgr2vr_b(thresh_ptr >> 8);
thresh = __lsx_vilvl_d(vec0, thresh);
b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
vec0 = __lsx_vreplgr2vr_b(b_limit_ptr >> 8);
b_limit = __lsx_vilvl_d(vec0, b_limit);
limit = __lsx_vreplgr2vr_b(limit_ptr);
vec0 = __lsx_vreplgr2vr_b(limit_ptr >> 8);
limit = __lsx_vilvl_d(vec0, limit);
/* mask and hev */
LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
hev, mask, flat);
/* flat4 */
VP9_FLAT4(p3, p2, p0, q0, q2, q3, flat);
/* filter4 */
VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
q1_out);
flat = __lsx_vilvh_d(flat, zero);
/* if flat is zero for all pixels, then no need to calculate other filter */
if (__lsx_bz_v(flat)) {
DUP2_ARG2(__lsx_vilvl_b, p0_out, p1_out, q1_out, q0_out, vec0, vec1);
vec2 = __lsx_vilvl_h(vec1, vec0);
vec3 = __lsx_vilvh_h(vec1, vec0);
DUP2_ARG2(__lsx_vilvh_b, p0_out, p1_out, q1_out, q0_out, vec0, vec1);
vec4 = __lsx_vilvl_h(vec1, vec0);
vec5 = __lsx_vilvh_h(vec1, vec0);
dst -= 2;
__lsx_vstelm_w(vec2, dst, 0, 0);
__lsx_vstelm_w(vec2, dst + stride, 0, 1);
__lsx_vstelm_w(vec2, dst + stride2, 0, 2);
__lsx_vstelm_w(vec2, dst + stride3, 0, 3);
dst += stride4;
__lsx_vstelm_w(vec3, dst, 0, 0);
__lsx_vstelm_w(vec3, dst + stride, 0, 1);
__lsx_vstelm_w(vec3, dst + stride2, 0, 2);
__lsx_vstelm_w(vec3, dst + stride3, 0, 3);
dst += stride4;
__lsx_vstelm_w(vec4, dst, 0, 0);
__lsx_vstelm_w(vec4, dst + stride, 0, 1);
__lsx_vstelm_w(vec4, dst + stride2, 0, 2);
__lsx_vstelm_w(vec4, dst + stride3, 0, 3);
dst += stride4;
__lsx_vstelm_w(vec5, dst, 0, 0);
__lsx_vstelm_w(vec5, dst + stride, 0, 1);
__lsx_vstelm_w(vec5, dst + stride2, 0, 2);
__lsx_vstelm_w(vec5, dst + stride3, 0, 3);
} else {
DUP4_ARG2(__lsx_vilvh_b, zero, p3, zero, p2, zero, p1, zero, p0,
p3_h, p2_h, p1_h, p0_h);
DUP4_ARG2(__lsx_vilvh_b, zero, q0, zero, q1, zero, q2, zero, q3,
q0_h, q1_h, q2_h, q3_h);
VP9_FILTER8(p3_h, p2_h, p1_h, p0_h, q0_h, q1_h, q2_h, q3_h, p2_filt8_h,
p1_filt8_h, p0_filt8_h, q0_filt8_h, q1_filt8_h, q2_filt8_h);
/* convert 16 bit output data into 8 bit */
DUP4_ARG2(__lsx_vpickev_b, p2_filt8_h, p2_filt8_h, p1_filt8_h,
p1_filt8_h, p0_filt8_h, p0_filt8_h, q0_filt8_h, q0_filt8_h,
p2_filt8_h, p1_filt8_h, p0_filt8_h, q0_filt8_h);
DUP2_ARG2(__lsx_vpickev_b, q1_filt8_h, q1_filt8_h, q2_filt8_h,
q2_filt8_h, q1_filt8_h, q2_filt8_h);
/* store pixel values */
p2 = __lsx_vbitsel_v(p2, p2_filt8_h, flat);
p1 = __lsx_vbitsel_v(p1_out, p1_filt8_h, flat);
p0 = __lsx_vbitsel_v(p0_out, p0_filt8_h, flat);
q0 = __lsx_vbitsel_v(q0_out, q0_filt8_h, flat);
q1 = __lsx_vbitsel_v(q1_out, q1_filt8_h, flat);
q2 = __lsx_vbitsel_v(q2, q2_filt8_h, flat);
DUP2_ARG2(__lsx_vilvl_b, p1, p2, q0, p0, vec0, vec1);
vec3 = __lsx_vilvl_h(vec1, vec0);
vec4 = __lsx_vilvh_h(vec1, vec0);
DUP2_ARG2(__lsx_vilvh_b, p1, p2, q0, p0, vec0, vec1);
vec6 = __lsx_vilvl_h(vec1, vec0);
vec7 = __lsx_vilvh_h(vec1, vec0);
vec2 = __lsx_vilvl_b(q2, q1);
vec5 = __lsx_vilvh_b(q2, q1);
dst -= 3;
__lsx_vstelm_w(vec3, dst, 0, 0);
__lsx_vstelm_h(vec2, dst, 4, 0);
dst += stride;
__lsx_vstelm_w(vec3, dst, 0, 1);
__lsx_vstelm_h(vec2, dst, 4, 1);
dst += stride;
__lsx_vstelm_w(vec3, dst, 0, 2);
__lsx_vstelm_h(vec2, dst, 4, 2);
dst += stride;
__lsx_vstelm_w(vec3, dst, 0, 3);
__lsx_vstelm_h(vec2, dst, 4, 3);
dst += stride;
__lsx_vstelm_w(vec4, dst, 0, 0);
__lsx_vstelm_h(vec2, dst, 4, 4);
dst += stride;
__lsx_vstelm_w(vec4, dst, 0, 1);
__lsx_vstelm_h(vec2, dst, 4, 5);
dst += stride;
__lsx_vstelm_w(vec4, dst, 0, 2);
__lsx_vstelm_h(vec2, dst, 4, 6);
dst += stride;
__lsx_vstelm_w(vec4, dst, 0, 3);
__lsx_vstelm_h(vec2, dst, 4, 7);
dst += stride;
__lsx_vstelm_w(vec6, dst, 0, 0);
__lsx_vstelm_h(vec5, dst, 4, 0);
dst += stride;
__lsx_vstelm_w(vec6, dst, 0, 1);
__lsx_vstelm_h(vec5, dst, 4, 1);
dst += stride;
__lsx_vstelm_w(vec6, dst, 0, 2);
__lsx_vstelm_h(vec5, dst, 4, 2);
dst += stride;
__lsx_vstelm_w(vec6, dst, 0, 3);
__lsx_vstelm_h(vec5, dst, 4, 3);
dst += stride;
__lsx_vstelm_w(vec7, dst, 0, 0);
__lsx_vstelm_h(vec5, dst, 4, 4);
dst += stride;
__lsx_vstelm_w(vec7, dst, 0, 1);
__lsx_vstelm_h(vec5, dst, 4, 5);
dst += stride;
__lsx_vstelm_w(vec7, dst, 0, 2);
__lsx_vstelm_h(vec5, dst, 4, 6);
dst += stride;
__lsx_vstelm_w(vec7, dst, 0, 3);
__lsx_vstelm_h(vec5, dst, 4, 7);
}
}
static void vp9_transpose_16x8_to_8x16(uint8_t *input, ptrdiff_t in_pitch,
uint8_t *output)
{
__m128i p7_org, p6_org, p5_org, p4_org, p3_org, p2_org, p1_org, p0_org;
__m128i tmp0, tmp1, tmp2, tmp3, tmp4, tmp5, tmp6, tmp7;
__m128i p7, p6, p5, p4, p3, p2, p1, p0, q0, q1, q2, q3, q4, q5, q6, q7;
ptrdiff_t in_pitch2 = in_pitch << 1;
ptrdiff_t in_pitch3 = in_pitch2 + in_pitch;
ptrdiff_t in_pitch4 = in_pitch2 << 1;
LSX_LD_8(input, in_pitch, in_pitch2, in_pitch3, in_pitch4,
p7_org, p6_org, p5_org, p4_org, p3_org, p2_org, p1_org, p0_org);
/* 8x8 transpose */
LSX_TRANSPOSE8x8_B(p7_org, p6_org, p5_org, p4_org, p3_org, p2_org, p1_org,
p0_org, p7, p6, p5, p4, p3, p2, p1, p0);
/* 8x8 transpose */
DUP4_ARG2(__lsx_vilvh_b, p5_org, p7_org, p4_org, p6_org, p1_org,
p3_org, p0_org, p2_org, tmp0, tmp1, tmp2, tmp3);
DUP2_ARG2(__lsx_vilvl_b, tmp1, tmp0, tmp3, tmp2, tmp4, tmp6);
DUP2_ARG2(__lsx_vilvh_b, tmp1, tmp0, tmp3, tmp2, tmp5, tmp7);
DUP2_ARG2(__lsx_vilvl_w, tmp6, tmp4, tmp7, tmp5, q0, q4);
DUP2_ARG2(__lsx_vilvh_w, tmp6, tmp4, tmp7, tmp5, q2, q6);
DUP4_ARG2(__lsx_vbsrl_v, q0, 8, q2, 8, q4, 8, q6, 8, q1, q3, q5, q7);
__lsx_vst(p7, output, 0);
__lsx_vst(p6, output, 16);
__lsx_vst(p5, output, 32);
__lsx_vst(p4, output, 48);
__lsx_vst(p3, output, 64);
__lsx_vst(p2, output, 80);
__lsx_vst(p1, output, 96);
__lsx_vst(p0, output, 112);
__lsx_vst(q0, output, 128);
__lsx_vst(q1, output, 144);
__lsx_vst(q2, output, 160);
__lsx_vst(q3, output, 176);
__lsx_vst(q4, output, 192);
__lsx_vst(q5, output, 208);
__lsx_vst(q6, output, 224);
__lsx_vst(q7, output, 240);
}
static void vp9_transpose_8x16_to_16x8(uint8_t *input, uint8_t *output,
ptrdiff_t out_pitch)
{
__m128i p7_o, p6_o, p5_o, p4_o, p3_o, p2_o, p1_o, p0_o;
__m128i p7, p6, p5, p4, p3, p2, p1, p0, q0, q1, q2, q3, q4, q5, q6, q7;
ptrdiff_t out_pitch2 = out_pitch << 1;
ptrdiff_t out_pitch3 = out_pitch2 + out_pitch;
ptrdiff_t out_pitch4 = out_pitch2 << 1;
DUP4_ARG2(__lsx_vld, input, 0, input, 16, input, 32, input, 48,
p7, p6, p5, p4);
DUP4_ARG2(__lsx_vld, input, 64, input, 80, input, 96, input, 112,
p3, p2, p1, p0);
DUP4_ARG2(__lsx_vld, input, 128, input, 144, input, 160, input, 176,
q0, q1, q2, q3);
DUP4_ARG2(__lsx_vld, input, 192, input, 208, input, 224, input, 240,
q4, q5, q6, q7);
LSX_TRANSPOSE16x8_B(p7, p6, p5, p4, p3, p2, p1, p0, q0, q1, q2, q3, q4, q5,
q6, q7, p7_o, p6_o, p5_o, p4_o, p3_o, p2_o, p1_o, p0_o);
LSX_ST_8(p7_o, p6_o, p5_o, p4_o, p3_o, p2_o, p1_o, p0_o,
output, out_pitch, out_pitch2, out_pitch3, out_pitch4);
}
static void vp9_transpose_16x16(uint8_t *input, int32_t in_stride,
uint8_t *output, int32_t out_stride)
{
__m128i row0, row1, row2, row3, row4, row5, row6, row7;
__m128i row8, row9, row10, row11, row12, row13, row14, row15;
__m128i tmp0, tmp1, tmp4, tmp5, tmp6, tmp7;
__m128i tmp2, tmp3;
__m128i p7, p6, p5, p4, p3, p2, p1, p0, q0, q1, q2, q3, q4, q5, q6, q7;
int32_t in_stride2 = in_stride << 1;
int32_t in_stride3 = in_stride2 + in_stride;
int32_t in_stride4 = in_stride2 << 1;
int32_t out_stride2 = out_stride << 1;
int32_t out_stride3 = out_stride2 + out_stride;
int32_t out_stride4 = out_stride2 << 1;
LSX_LD_8(input, in_stride, in_stride2, in_stride3, in_stride4,
row0, row1, row2, row3, row4, row5, row6, row7);
input += in_stride4;
LSX_LD_8(input, in_stride, in_stride2, in_stride3, in_stride4,
row8, row9, row10, row11, row12, row13, row14, row15);
LSX_TRANSPOSE16x8_B(row0, row1, row2, row3, row4, row5, row6, row7,
row8, row9, row10, row11, row12, row13, row14, row15,
p7, p6, p5, p4, p3, p2, p1, p0);
/* transpose 16x8 matrix into 8x16 */
/* total 8 intermediate register and 32 instructions */
q7 = __lsx_vpackod_d(row8, row0);
q6 = __lsx_vpackod_d(row9, row1);
q5 = __lsx_vpackod_d(row10, row2);
q4 = __lsx_vpackod_d(row11, row3);
q3 = __lsx_vpackod_d(row12, row4);
q2 = __lsx_vpackod_d(row13, row5);
q1 = __lsx_vpackod_d(row14, row6);
q0 = __lsx_vpackod_d(row15, row7);
DUP2_ARG2(__lsx_vpackev_b, q6, q7, q4, q5, tmp0, tmp1);
DUP2_ARG2(__lsx_vpackod_b, q6, q7, q4, q5, tmp4, tmp5);
DUP2_ARG2(__lsx_vpackev_b, q2, q3, q0, q1, q5, q7);
DUP2_ARG2(__lsx_vpackod_b, q2, q3, q0, q1, tmp6, tmp7);
DUP2_ARG2(__lsx_vpackev_h, tmp1, tmp0, q7, q5, tmp2, tmp3);
q0 = __lsx_vpackev_w(tmp3, tmp2);
q4 = __lsx_vpackod_w(tmp3, tmp2);
tmp2 = __lsx_vpackod_h(tmp1, tmp0);
tmp3 = __lsx_vpackod_h(q7, q5);
q2 = __lsx_vpackev_w(tmp3, tmp2);
q6 = __lsx_vpackod_w(tmp3, tmp2);
DUP2_ARG2(__lsx_vpackev_h, tmp5, tmp4, tmp7, tmp6, tmp2, tmp3);
q1 = __lsx_vpackev_w(tmp3, tmp2);
q5 = __lsx_vpackod_w(tmp3, tmp2);
tmp2 = __lsx_vpackod_h(tmp5, tmp4);
tmp3 = __lsx_vpackod_h(tmp7, tmp6);
q3 = __lsx_vpackev_w(tmp3, tmp2);
q7 = __lsx_vpackod_w(tmp3, tmp2);
LSX_ST_8(p7, p6, p5, p4, p3, p2, p1, p0, output, out_stride,
out_stride2, out_stride3, out_stride4);
output += out_stride4;
LSX_ST_8(q0, q1, q2, q3, q4, q5, q6, q7, output, out_stride,
out_stride2, out_stride3, out_stride4);
}
static int32_t vp9_vt_lpf_t4_and_t8_8w(uint8_t *src, uint8_t *filter48,
uint8_t *src_org, int32_t pitch_org,
int32_t b_limit_ptr,
int32_t limit_ptr,
int32_t thresh_ptr)
{
__m128i p3, p2, p1, p0, q3, q2, q1, q0;
__m128i p2_out, p1_out, p0_out, q0_out, q1_out, q2_out;
__m128i flat, mask, hev, thresh, b_limit, limit;
__m128i p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l;
__m128i p2_filt8_l, p1_filt8_l, p0_filt8_l;
__m128i q0_filt8_l, q1_filt8_l, q2_filt8_l;
__m128i vec0, vec1, vec2, vec3;
__m128i zero = __lsx_vldi(0);
/* load vector elements */
DUP4_ARG2(__lsx_vld, src, -64, src, -48, src, -32, src, -16,
p3, p2, p1, p0);
DUP4_ARG2(__lsx_vld, src, 0, src, 16, src, 32, src, 48, q0, q1, q2, q3);
thresh = __lsx_vreplgr2vr_b(thresh_ptr);
b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
limit = __lsx_vreplgr2vr_b(limit_ptr);
/* mask and hev */
LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
hev, mask, flat);
/* flat4 */
VP9_FLAT4(p3, p2, p0, q0, q2, q3, flat);
/* filter4 */
VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
q1_out);
flat = __lsx_vilvl_d(zero, flat);
/* if flat is zero for all pixels, then no need to calculate other filter */
if (__lsx_bz_v(flat)) {
DUP2_ARG2(__lsx_vilvl_b, p0_out, p1_out, q1_out, q0_out, vec0, vec1);
vec2 = __lsx_vilvl_h(vec1, vec0);
vec3 = __lsx_vilvh_h(vec1, vec0);
src_org -= 2;
__lsx_vstelm_w(vec2, src_org, 0, 0);
src_org += pitch_org;
__lsx_vstelm_w(vec2, src_org, 0, 1);
src_org += pitch_org;
__lsx_vstelm_w(vec2, src_org, 0, 2);
src_org += pitch_org;
__lsx_vstelm_w(vec2, src_org, 0, 3);
src_org += pitch_org;
__lsx_vstelm_w(vec3, src_org, 0, 0);
src_org += pitch_org;
__lsx_vstelm_w(vec3, src_org, 0, 1);
src_org += pitch_org;
__lsx_vstelm_w(vec3, src_org, 0, 2);
src_org += pitch_org;
__lsx_vstelm_w(vec3, src_org, 0, 3);
return 1;
} else {
DUP4_ARG2(__lsx_vilvl_b, zero, p3, zero, p2, zero, p1, zero, p0,
p3_l, p2_l, p1_l, p0_l);
DUP4_ARG2(__lsx_vilvl_b, zero, q0, zero, q1, zero, q2, zero, q3,
q0_l, q1_l, q2_l, q3_l);
VP9_FILTER8(p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l, p2_filt8_l,
p1_filt8_l, p0_filt8_l, q0_filt8_l, q1_filt8_l, q2_filt8_l);
/* convert 16 bit output data into 8 bit */
p2_l = __lsx_vpickev_b(p2_filt8_l, p2_filt8_l);
p1_l = __lsx_vpickev_b(p1_filt8_l, p1_filt8_l);
p0_l = __lsx_vpickev_b(p0_filt8_l, p0_filt8_l);
q0_l = __lsx_vpickev_b(q0_filt8_l, q0_filt8_l);
q1_l = __lsx_vpickev_b(q1_filt8_l, q1_filt8_l);
q2_l = __lsx_vpickev_b(q2_filt8_l, q2_filt8_l);
/* store pixel values */
p2_out = __lsx_vbitsel_v(p2, p2_l, flat);
p1_out = __lsx_vbitsel_v(p1_out, p1_l, flat);
p0_out = __lsx_vbitsel_v(p0_out, p0_l, flat);
q0_out = __lsx_vbitsel_v(q0_out, q0_l, flat);
q1_out = __lsx_vbitsel_v(q1_out, q1_l, flat);
q2_out = __lsx_vbitsel_v(q2, q2_l, flat);
__lsx_vst(p2_out, filter48, 0);
__lsx_vst(p1_out, filter48, 16);
__lsx_vst(p0_out, filter48, 32);
__lsx_vst(q0_out, filter48, 48);
__lsx_vst(q1_out, filter48, 64);
__lsx_vst(q2_out, filter48, 80);
__lsx_vst(flat, filter48, 96);
return 0;
}
}
static int32_t vp9_vt_lpf_t16_8w(uint8_t *dst, uint8_t *dst_org,
ptrdiff_t stride,
uint8_t *filter48)
{
__m128i zero = __lsx_vldi(0);
__m128i filter8, flat, flat2;
__m128i p7, p6, p5, p4, p3, p2, p1, p0, q0, q1, q2, q3, q4, q5, q6, q7;
v8u16 p7_l_in, p6_l_in, p5_l_in, p4_l_in;
v8u16 p3_l_in, p2_l_in, p1_l_in, p0_l_in;
v8u16 q7_l_in, q6_l_in, q5_l_in, q4_l_in;
v8u16 q3_l_in, q2_l_in, q1_l_in, q0_l_in;
v8u16 tmp0_l, tmp1_l;
__m128i out_l;
uint8_t *dst_tmp = dst - 128;
/* load vector elements */
DUP4_ARG2(__lsx_vld, dst_tmp, 0, dst_tmp, 16, dst_tmp, 32,
dst_tmp, 48, p7, p6, p5, p4);
DUP4_ARG2(__lsx_vld, dst_tmp, 64, dst_tmp, 80, dst_tmp, 96,
dst_tmp, 112, p3, p2, p1, p0);
DUP4_ARG2(__lsx_vld, dst, 0, dst, 16, dst, 32, dst, 48, q0, q1, q2, q3);
DUP4_ARG2(__lsx_vld, dst, 64, dst, 80, dst, 96, dst, 112, q4, q5, q6, q7);
flat = __lsx_vld(filter48, 96);
VP9_FLAT5(p7, p6, p5, p4, p0, q0, q4, q5, q6, q7, flat, flat2);
/* if flat2 is zero for all pixels, then no need to calculate other filter */
if (__lsx_bz_v(flat2)) {
__m128i vec0, vec1, vec2, vec3, vec4;
DUP4_ARG2(__lsx_vld, filter48, 0, filter48, 16, filter48, 32,
filter48, 48, p2, p1, p0, q0);
DUP2_ARG2(__lsx_vld, filter48, 64, filter48, 80, q1, q2);
DUP2_ARG2(__lsx_vilvl_b, p1, p2, q0, p0, vec0, vec1);
vec3 = __lsx_vilvl_h(vec1, vec0);
vec4 = __lsx_vilvh_h(vec1, vec0);
vec2 = __lsx_vilvl_b(q2, q1);
dst_org -= 3;
__lsx_vstelm_w(vec3, dst_org, 0, 0);
__lsx_vstelm_h(vec2, dst_org, 4, 0);
dst_org += stride;
__lsx_vstelm_w(vec3, dst_org, 0, 1);
__lsx_vstelm_h(vec2, dst_org, 4, 1);
dst_org += stride;
__lsx_vstelm_w(vec3, dst_org, 0, 2);
__lsx_vstelm_h(vec2, dst_org, 4, 2);
dst_org += stride;
__lsx_vstelm_w(vec3, dst_org, 0, 3);
__lsx_vstelm_h(vec2, dst_org, 4, 3);
dst_org += stride;
__lsx_vstelm_w(vec4, dst_org, 0, 0);
__lsx_vstelm_h(vec2, dst_org, 4, 4);
dst_org += stride;
__lsx_vstelm_w(vec4, dst_org, 0, 1);
__lsx_vstelm_h(vec2, dst_org, 4, 5);
dst_org += stride;
__lsx_vstelm_w(vec4, dst_org, 0, 2);
__lsx_vstelm_h(vec2, dst_org, 4, 6);
dst_org += stride;
__lsx_vstelm_w(vec4, dst_org, 0, 3);
__lsx_vstelm_h(vec2, dst_org, 4, 7);
return 1;
} else {
dst -= 7 * 16;
p7_l_in = (v8u16)__lsx_vilvl_b(zero, p7);
p6_l_in = (v8u16)__lsx_vilvl_b(zero, p6);
p5_l_in = (v8u16)__lsx_vilvl_b(zero, p5);
p4_l_in = (v8u16)__lsx_vilvl_b(zero, p4);
p3_l_in = (v8u16)__lsx_vilvl_b(zero, p3);
p2_l_in = (v8u16)__lsx_vilvl_b(zero, p2);
p1_l_in = (v8u16)__lsx_vilvl_b(zero, p1);
p0_l_in = (v8u16)__lsx_vilvl_b(zero, p0);
q0_l_in = (v8u16)__lsx_vilvl_b(zero, q0);
tmp0_l = p7_l_in << 3;
tmp0_l -= p7_l_in;
tmp0_l += p6_l_in;
tmp0_l += q0_l_in;
tmp1_l = p6_l_in + p5_l_in;
tmp1_l += p4_l_in;
tmp1_l += p3_l_in;
tmp1_l += p2_l_in;
tmp1_l += p1_l_in;
tmp1_l += p0_l_in;
tmp1_l += tmp0_l;
out_l =__lsx_vsrari_h((__m128i)tmp1_l, 4);
out_l =__lsx_vpickev_b(out_l, out_l);
p6 = __lsx_vbitsel_v(p6, out_l, flat2);
__lsx_vstelm_d(p6, dst, 0, 0);
dst += 16;
/* p5 */
q1_l_in = (v8u16)__lsx_vilvl_b(zero, q1);
tmp0_l = p5_l_in - p6_l_in;
tmp0_l += q1_l_in;
tmp0_l -= p7_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
out_l = __lsx_vpickev_b(out_l, out_l);
p5 = __lsx_vbitsel_v(p5, out_l, flat2);
__lsx_vstelm_d(p5, dst, 0, 0);
dst += 16;
/* p4 */
q2_l_in = (v8u16)__lsx_vilvl_b(zero, q2);
tmp0_l = p4_l_in - p5_l_in;
tmp0_l += q2_l_in;
tmp0_l -= p7_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
out_l = __lsx_vpickev_b(out_l, out_l);
p4 = __lsx_vbitsel_v(p4, out_l, flat2);
__lsx_vstelm_d(p4, dst, 0, 0);
dst += 16;
/* p3 */
q3_l_in = (v8u16)__lsx_vilvl_b(zero, q3);
tmp0_l = p3_l_in - p4_l_in;
tmp0_l += q3_l_in;
tmp0_l -= p7_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
out_l = __lsx_vpickev_b(out_l, out_l);
p3 = __lsx_vbitsel_v(p3, out_l, flat2);
__lsx_vstelm_d(p3, dst, 0, 0);
dst += 16;
/* p2 */
q4_l_in = (v8u16)__lsx_vilvl_b(zero, q4);
filter8 = __lsx_vld(filter48, 0);
tmp0_l = p2_l_in - p3_l_in;
tmp0_l += q4_l_in;
tmp0_l -= p7_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
out_l = __lsx_vpickev_b(out_l, out_l);
filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
__lsx_vstelm_d(filter8, dst, 0, 0);
dst += 16;
/* p1 */
q5_l_in = (v8u16)__lsx_vilvl_b(zero, q5);
filter8 = __lsx_vld(filter48, 16);
tmp0_l = p1_l_in - p2_l_in;
tmp0_l += q5_l_in;
tmp0_l -= p7_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
out_l = __lsx_vpickev_b(out_l, out_l);
filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
__lsx_vstelm_d(filter8, dst, 0, 0);
dst += 16;
/* p0 */
q6_l_in = (v8u16)__lsx_vilvl_b(zero, q6);
filter8 = __lsx_vld(filter48, 32);
tmp0_l = p0_l_in - p1_l_in;
tmp0_l += q6_l_in;
tmp0_l -= p7_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
out_l = __lsx_vpickev_b(out_l, out_l);
filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
__lsx_vstelm_d(filter8, dst, 0, 0);
dst += 16;
/* q0 */
q7_l_in = (v8u16)__lsx_vilvl_b(zero, q7);
filter8 = __lsx_vld(filter48, 48);
tmp0_l = q7_l_in - p0_l_in;
tmp0_l += q0_l_in;
tmp0_l -= p7_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((v8i16) tmp1_l, 4);
out_l = __lsx_vpickev_b(out_l, out_l);
filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
__lsx_vstelm_d(filter8, dst, 0, 0);
dst += 16;
/* q1 */
filter8 = __lsx_vld(filter48, 64);
tmp0_l = q7_l_in - q0_l_in;
tmp0_l += q1_l_in;
tmp0_l -= p6_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
out_l = __lsx_vpickev_b(out_l, out_l);
filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
__lsx_vstelm_d(filter8, dst, 0, 0);
dst += 16;
/* q2 */
filter8 = __lsx_vld(filter48, 80);
tmp0_l = q7_l_in - q1_l_in;
tmp0_l += q2_l_in;
tmp0_l -= p5_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
out_l = __lsx_vpickev_b(out_l, out_l);
filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
__lsx_vstelm_d(filter8, dst, 0, 0);
dst += 16;
/* q3 */
tmp0_l = q7_l_in - q2_l_in;
tmp0_l += q3_l_in;
tmp0_l -= p4_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
out_l = __lsx_vpickev_b(out_l, out_l);
q3 = __lsx_vbitsel_v(q3, out_l, flat2);
__lsx_vstelm_d(q3, dst, 0, 0);
dst += 16;
/* q4 */
tmp0_l = q7_l_in - q3_l_in;
tmp0_l += q4_l_in;
tmp0_l -= p3_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
out_l = __lsx_vpickev_b(out_l, out_l);
q4 = __lsx_vbitsel_v(q4, out_l, flat2);
__lsx_vstelm_d(q4, dst, 0, 0);
dst += 16;
/* q5 */
tmp0_l = q7_l_in - q4_l_in;
tmp0_l += q5_l_in;
tmp0_l -= p2_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
out_l = __lsx_vpickev_b(out_l, out_l);
q5 = __lsx_vbitsel_v(q5, out_l, flat2);
__lsx_vstelm_d(q5, dst, 0, 0);
dst += 16;
/* q6 */
tmp0_l = q7_l_in - q5_l_in;
tmp0_l += q6_l_in;
tmp0_l -= p1_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
out_l = __lsx_vpickev_b(out_l, out_l);
q6 = __lsx_vbitsel_v(q6, out_l, flat2);
__lsx_vstelm_d(q6, dst, 0, 0);
return 0;
}
}
void ff_loop_filter_h_16_8_lsx(uint8_t *dst, ptrdiff_t stride,
int32_t b_limit_ptr,
int32_t limit_ptr,
int32_t thresh_ptr)
{
uint8_t early_exit = 0;
uint8_t transposed_input[16 * 24] __attribute__ ((aligned(16)));
uint8_t *filter48 = &transposed_input[16 * 16];
vp9_transpose_16x8_to_8x16(dst - 8, stride, transposed_input);
early_exit = vp9_vt_lpf_t4_and_t8_8w((transposed_input + 16 * 8),
&filter48[0], dst, stride,
b_limit_ptr, limit_ptr, thresh_ptr);
if (0 == early_exit) {
early_exit = vp9_vt_lpf_t16_8w((transposed_input + 16 * 8), dst, stride,
&filter48[0]);
if (0 == early_exit) {
vp9_transpose_8x16_to_16x8(transposed_input, dst - 8, stride);
}
}
}
static int32_t vp9_vt_lpf_t4_and_t8_16w(uint8_t *dst, uint8_t *filter48,
uint8_t *dst_org, ptrdiff_t stride,
int32_t b_limit_ptr,
int32_t limit_ptr,
int32_t thresh_ptr)
{
ptrdiff_t stride2 = stride << 1;
ptrdiff_t stride3 = stride2 + stride;
ptrdiff_t stride4 = stride2 << 1;
__m128i p3, p2, p1, p0, q3, q2, q1, q0;
__m128i p2_out, p1_out, p0_out, q0_out, q1_out, q2_out;
__m128i flat, mask, hev, thresh, b_limit, limit;
__m128i p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l;
__m128i p3_h, p2_h, p1_h, p0_h, q0_h, q1_h, q2_h, q3_h;
__m128i p2_filt8_l, p1_filt8_l, p0_filt8_l;
__m128i q0_filt8_l, q1_filt8_l, q2_filt8_l;
__m128i p2_filt8_h, p1_filt8_h, p0_filt8_h;
__m128i q0_filt8_h, q1_filt8_h, q2_filt8_h;
__m128i vec0, vec1, vec2, vec3, vec4, vec5;
__m128i zero = __lsx_vldi(0);
/* load vector elements */
DUP4_ARG2(__lsx_vld, dst, -64, dst, -48, dst, -32, dst, -16,
p3, p2, p1, p0);
DUP4_ARG2(__lsx_vld, dst, 0, dst, 16, dst, 32, dst, 48, q0, q1, q2, q3);
thresh = __lsx_vreplgr2vr_b(thresh_ptr);
b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
limit = __lsx_vreplgr2vr_b(limit_ptr);
/* mask and hev */
LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
hev, mask, flat);
/* flat4 */
VP9_FLAT4(p3, p2, p0, q0, q2, q3, flat);
/* filter4 */
VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
q1_out);
/* if flat is zero for all pixels, then no need to calculate other filter */
if (__lsx_bz_v(flat)) {
DUP2_ARG2(__lsx_vilvl_b, p0_out, p1_out, q1_out, q0_out, vec0, vec1);
vec2 = __lsx_vilvl_h(vec1, vec0);
vec3 = __lsx_vilvh_h(vec1, vec0);
DUP2_ARG2(__lsx_vilvh_b, p0_out, p1_out, q1_out, q0_out, vec0, vec1);
vec4 = __lsx_vilvl_h(vec1, vec0);
vec5 = __lsx_vilvh_h(vec1, vec0);
dst_org -= 2;
__lsx_vstelm_w(vec2, dst_org, 0, 0);
__lsx_vstelm_w(vec2, dst_org + stride, 0, 1);
__lsx_vstelm_w(vec2, dst_org + stride2, 0, 2);
__lsx_vstelm_w(vec2, dst_org + stride3, 0, 3);
dst_org += stride4;
__lsx_vstelm_w(vec3, dst_org, 0, 0);
__lsx_vstelm_w(vec3, dst_org + stride, 0, 1);
__lsx_vstelm_w(vec3, dst_org + stride2, 0, 2);
__lsx_vstelm_w(vec3, dst_org + stride3, 0, 3);
dst_org += stride4;
__lsx_vstelm_w(vec4, dst_org, 0, 0);
__lsx_vstelm_w(vec4, dst_org + stride, 0, 1);
__lsx_vstelm_w(vec4, dst_org + stride2, 0, 2);
__lsx_vstelm_w(vec4, dst_org + stride3, 0, 3);
dst_org += stride4;
__lsx_vstelm_w(vec5, dst_org, 0, 0);
__lsx_vstelm_w(vec5, dst_org + stride, 0, 1);
__lsx_vstelm_w(vec5, dst_org + stride2, 0, 2);
__lsx_vstelm_w(vec5, dst_org + stride3, 0, 3);
return 1;
} else {
DUP4_ARG2(__lsx_vilvl_b, zero, p3, zero, p2, zero, p1, zero, p0,
p3_l, p2_l, p1_l, p0_l);
DUP4_ARG2(__lsx_vilvl_b, zero, q0, zero, q1, zero, q2, zero, q3,
q0_l, q1_l, q2_l, q3_l);
VP9_FILTER8(p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l, p2_filt8_l,
p1_filt8_l, p0_filt8_l, q0_filt8_l, q1_filt8_l, q2_filt8_l);
DUP4_ARG2(__lsx_vilvh_b, zero, p3, zero, p2, zero, p1, zero, p0,
p3_h, p2_h, p1_h, p0_h);
DUP4_ARG2(__lsx_vilvh_b, zero, q0, zero, q1, zero, q2, zero, q3,
q0_h, q1_h, q2_h, q3_h);
VP9_FILTER8(p3_h, p2_h, p1_h, p0_h, q0_h, q1_h, q2_h, q3_h, p2_filt8_h,
p1_filt8_h, p0_filt8_h, q0_filt8_h, q1_filt8_h, q2_filt8_h);
/* convert 16 bit output data into 8 bit */
DUP4_ARG2(__lsx_vpickev_b, p2_filt8_h, p2_filt8_l, p1_filt8_h,
p1_filt8_l, p0_filt8_h, p0_filt8_l, q0_filt8_h,
q0_filt8_l, p2_filt8_l, p1_filt8_l, p0_filt8_l,
q0_filt8_l);
DUP2_ARG2(__lsx_vpickev_b, q1_filt8_h, q1_filt8_l, q2_filt8_h,
q2_filt8_l, q1_filt8_l, q2_filt8_l);
/* store pixel values */
p2_out = __lsx_vbitsel_v(p2, p2_filt8_l, flat);
p1_out = __lsx_vbitsel_v(p1_out, p1_filt8_l, flat);
p0_out = __lsx_vbitsel_v(p0_out, p0_filt8_l, flat);
q0_out = __lsx_vbitsel_v(q0_out, q0_filt8_l, flat);
q1_out = __lsx_vbitsel_v(q1_out, q1_filt8_l, flat);
q2_out = __lsx_vbitsel_v(q2, q2_filt8_l, flat);
__lsx_vst(p2_out, filter48, 0);
__lsx_vst(p1_out, filter48, 16);
__lsx_vst(p0_out, filter48, 32);
__lsx_vst(q0_out, filter48, 48);
__lsx_vst(q1_out, filter48, 64);
__lsx_vst(q2_out, filter48, 80);
__lsx_vst(flat, filter48, 96);
return 0;
}
}
static int32_t vp9_vt_lpf_t16_16w(uint8_t *dst, uint8_t *dst_org,
ptrdiff_t stride,
uint8_t *filter48)
{
__m128i zero = __lsx_vldi(0);
__m128i flat, flat2, filter8;
__m128i p7, p6, p5, p4, p3, p2, p1, p0, q0, q1, q2, q3, q4, q5, q6, q7;
v8u16 p7_l_in, p6_l_in, p5_l_in, p4_l_in;
v8u16 p3_l_in, p2_l_in, p1_l_in, p0_l_in;
v8u16 q7_l_in, q6_l_in, q5_l_in, q4_l_in;
v8u16 q3_l_in, q2_l_in, q1_l_in, q0_l_in;
v8u16 p7_h_in, p6_h_in, p5_h_in, p4_h_in;
v8u16 p3_h_in, p2_h_in, p1_h_in, p0_h_in;
v8u16 q7_h_in, q6_h_in, q5_h_in, q4_h_in;
v8u16 q3_h_in, q2_h_in, q1_h_in, q0_h_in;
v8u16 tmp0_l, tmp1_l, tmp0_h, tmp1_h;
__m128i out_l, out_h;
uint8_t *dst_tmp = dst - 128;
flat = __lsx_vld(filter48, 96);
DUP4_ARG2(__lsx_vld, dst_tmp, 0, dst_tmp, 16, dst_tmp, 32,
dst_tmp, 48, p7, p6, p5, p4);
DUP4_ARG2(__lsx_vld, dst_tmp, 64, dst_tmp, 80, dst_tmp, 96,
dst_tmp, 112, p3, p2, p1, p0);
DUP4_ARG2(__lsx_vld, dst, 0, dst, 16, dst, 32, dst, 48, q0, q1, q2, q3);
DUP4_ARG2(__lsx_vld, dst, 64, dst, 80, dst, 96, dst, 112, q4, q5, q6, q7);
VP9_FLAT5(p7, p6, p5, p4, p0, q0, q4, q5, q6, q7, flat, flat2);
/* if flat2 is zero for all pixels, then no need to calculate other filter */
if (__lsx_bz_v(flat2)) {
__m128i vec0, vec1, vec2, vec3, vec4, vec5, vec6, vec7;
DUP4_ARG2(__lsx_vld, filter48, 0, filter48, 16, filter48, 32,
filter48, 48, p2, p1, p0, q0);
DUP2_ARG2(__lsx_vld, filter48, 64, filter48, 80, q1, q2);
DUP2_ARG2(__lsx_vilvl_b, p1, p2, q0, p0, vec0, vec1);
vec3 = __lsx_vilvl_h(vec1, vec0);
vec4 = __lsx_vilvh_h(vec1, vec0);
DUP2_ARG2(__lsx_vilvh_b, p1, p2, q0, p0, vec0, vec1);
vec6 = __lsx_vilvl_h(vec1, vec0);
vec7 = __lsx_vilvh_h(vec1, vec0);
vec2 = __lsx_vilvl_b(q2, q1);
vec5 = __lsx_vilvh_b(q2, q1);
dst_org -= 3;
__lsx_vstelm_w(vec3, dst_org, 0, 0);
__lsx_vstelm_h(vec2, dst_org, 4, 0);
dst_org += stride;
__lsx_vstelm_w(vec3, dst_org, 0, 1);
__lsx_vstelm_h(vec2, dst_org, 4, 1);
dst_org += stride;
__lsx_vstelm_w(vec3, dst_org, 0, 2);
__lsx_vstelm_h(vec2, dst_org, 4, 2);
dst_org += stride;
__lsx_vstelm_w(vec3, dst_org, 0, 3);
__lsx_vstelm_h(vec2, dst_org, 4, 3);
dst_org += stride;
__lsx_vstelm_w(vec4, dst_org, 0, 0);
__lsx_vstelm_h(vec2, dst_org, 4, 4);
dst_org += stride;
__lsx_vstelm_w(vec4, dst_org, 0, 1);
__lsx_vstelm_h(vec2, dst_org, 4, 5);
dst_org += stride;
__lsx_vstelm_w(vec4, dst_org, 0, 2);
__lsx_vstelm_h(vec2, dst_org, 4, 6);
dst_org += stride;
__lsx_vstelm_w(vec4, dst_org, 0, 3);
__lsx_vstelm_h(vec2, dst_org, 4, 7);
dst_org += stride;
__lsx_vstelm_w(vec6, dst_org, 0, 0);
__lsx_vstelm_h(vec5, dst_org, 4, 0);
dst_org += stride;
__lsx_vstelm_w(vec6, dst_org, 0, 1);
__lsx_vstelm_h(vec5, dst_org, 4, 1);
dst_org += stride;
__lsx_vstelm_w(vec6, dst_org, 0, 2);
__lsx_vstelm_h(vec5, dst_org, 4, 2);
dst_org += stride;
__lsx_vstelm_w(vec6, dst_org, 0, 3);
__lsx_vstelm_h(vec5, dst_org, 4, 3);
dst_org += stride;
__lsx_vstelm_w(vec7, dst_org, 0, 0);
__lsx_vstelm_h(vec5, dst_org, 4, 4);
dst_org += stride;
__lsx_vstelm_w(vec7, dst_org, 0, 1);
__lsx_vstelm_h(vec5, dst_org, 4, 5);
dst_org += stride;
__lsx_vstelm_w(vec7, dst_org, 0, 2);
__lsx_vstelm_h(vec5, dst_org, 4, 6);
dst_org += stride;
__lsx_vstelm_w(vec7, dst_org, 0, 3);
__lsx_vstelm_h(vec5, dst_org, 4, 7);
return 1;
} else {
dst -= 7 * 16;
p7_l_in = (v8u16)__lsx_vilvl_b(zero, p7);
p6_l_in = (v8u16)__lsx_vilvl_b(zero, p6);
p5_l_in = (v8u16)__lsx_vilvl_b(zero, p5);
p4_l_in = (v8u16)__lsx_vilvl_b(zero, p4);
p3_l_in = (v8u16)__lsx_vilvl_b(zero, p3);
p2_l_in = (v8u16)__lsx_vilvl_b(zero, p2);
p1_l_in = (v8u16)__lsx_vilvl_b(zero, p1);
p0_l_in = (v8u16)__lsx_vilvl_b(zero, p0);
q0_l_in = (v8u16)__lsx_vilvl_b(zero, q0);
tmp0_l = p7_l_in << 3;
tmp0_l -= p7_l_in;
tmp0_l += p6_l_in;
tmp0_l += q0_l_in;
tmp1_l = p6_l_in + p5_l_in;
tmp1_l += p4_l_in;
tmp1_l += p3_l_in;
tmp1_l += p2_l_in;
tmp1_l += p1_l_in;
tmp1_l += p0_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
p7_h_in = (v8u16)__lsx_vilvh_b(zero, p7);
p6_h_in = (v8u16)__lsx_vilvh_b(zero, p6);
p5_h_in = (v8u16)__lsx_vilvh_b(zero, p5);
p4_h_in = (v8u16)__lsx_vilvh_b(zero, p4);
p3_h_in = (v8u16)__lsx_vilvh_b(zero, p3);
p2_h_in = (v8u16)__lsx_vilvh_b(zero, p2);
p1_h_in = (v8u16)__lsx_vilvh_b(zero, p1);
p0_h_in = (v8u16)__lsx_vilvh_b(zero, p0);
q0_h_in = (v8u16)__lsx_vilvh_b(zero, q0);
tmp0_h = p7_h_in << 3;
tmp0_h -= p7_h_in;
tmp0_h += p6_h_in;
tmp0_h += q0_h_in;
tmp1_h = p6_h_in + p5_h_in;
tmp1_h += p4_h_in;
tmp1_h += p3_h_in;
tmp1_h += p2_h_in;
tmp1_h += p1_h_in;
tmp1_h += p0_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
p6 = __lsx_vbitsel_v(p6, out_l, flat2);
__lsx_vst(p6, dst, 0);
/* p5 */
q1_l_in = (v8u16)__lsx_vilvl_b(zero, q1);
tmp0_l = p5_l_in - p6_l_in;
tmp0_l += q1_l_in;
tmp0_l -= p7_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
q1_h_in = (v8u16)__lsx_vilvh_b(zero, q1);
tmp0_h = p5_h_in - p6_h_in;
tmp0_h += q1_h_in;
tmp0_h -= p7_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
p5 = __lsx_vbitsel_v(p5, out_l, flat2);
__lsx_vst(p5, dst, 16);
/* p4 */
q2_l_in = (v8u16)__lsx_vilvl_b(zero, q2);
tmp0_l = p4_l_in - p5_l_in;
tmp0_l += q2_l_in;
tmp0_l -= p7_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
q2_h_in = (v8u16)__lsx_vilvh_b(zero, q2);
tmp0_h = p4_h_in - p5_h_in;
tmp0_h += q2_h_in;
tmp0_h -= p7_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
p4 = __lsx_vbitsel_v(p4, out_l, flat2);
__lsx_vst(p4, dst, 16*2);
/* p3 */
q3_l_in = (v8u16)__lsx_vilvl_b(zero, q3);
tmp0_l = p3_l_in - p4_l_in;
tmp0_l += q3_l_in;
tmp0_l -= p7_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
q3_h_in = (v8u16)__lsx_vilvh_b(zero, q3);
tmp0_h = p3_h_in - p4_h_in;
tmp0_h += q3_h_in;
tmp0_h -= p7_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
p3 = __lsx_vbitsel_v(p3, out_l, flat2);
__lsx_vst(p3, dst, 16*3);
/* p2 */
q4_l_in = (v8u16)__lsx_vilvl_b(zero, q4);
filter8 = __lsx_vld(filter48, 0);
tmp0_l = p2_l_in - p3_l_in;
tmp0_l += q4_l_in;
tmp0_l -= p7_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
q4_h_in = (v8u16)__lsx_vilvh_b(zero, q4);
tmp0_h = p2_h_in - p3_h_in;
tmp0_h += q4_h_in;
tmp0_h -= p7_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
__lsx_vst(filter8, dst, 16*4);
/* p1 */
q5_l_in = (v8u16)__lsx_vilvl_b(zero, q5);
filter8 = __lsx_vld(filter48, 16);
tmp0_l = p1_l_in - p2_l_in;
tmp0_l += q5_l_in;
tmp0_l -= p7_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
q5_h_in = (v8u16)__lsx_vilvh_b(zero, q5);
tmp0_h = p1_h_in - p2_h_in;
tmp0_h += q5_h_in;
tmp0_h -= p7_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)(tmp1_h), 4);
out_l = __lsx_vpickev_b(out_h, out_l);
filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
__lsx_vst(filter8, dst, 16*5);
/* p0 */
q6_l_in = (v8u16)__lsx_vilvl_b(zero, q6);
filter8 = __lsx_vld(filter48, 32);
tmp0_l = p0_l_in - p1_l_in;
tmp0_l += q6_l_in;
tmp0_l -= p7_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
q6_h_in = (v8u16)__lsx_vilvh_b(zero, q6);
tmp0_h = p0_h_in - p1_h_in;
tmp0_h += q6_h_in;
tmp0_h -= p7_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
__lsx_vst(filter8, dst, 16*6);
/* q0 */
q7_l_in = (v8u16)__lsx_vilvl_b(zero, q7);
filter8 = __lsx_vld(filter48, 48);
tmp0_l = q7_l_in - p0_l_in;
tmp0_l += q0_l_in;
tmp0_l -= p7_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
q7_h_in = (v8u16)__lsx_vilvh_b(zero, q7);
tmp0_h = q7_h_in - p0_h_in;
tmp0_h += q0_h_in;
tmp0_h -= p7_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
__lsx_vst(filter8, dst, 16*7);
/* q1 */
filter8 = __lsx_vld(filter48, 64);
tmp0_l = q7_l_in - q0_l_in;
tmp0_l += q1_l_in;
tmp0_l -= p6_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
tmp0_h = q7_h_in - q0_h_in;
tmp0_h += q1_h_in;
tmp0_h -= p6_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
__lsx_vst(filter8, dst, 16*8);
/* q2 */
filter8 = __lsx_vld(filter48, 80);
tmp0_l = q7_l_in - q1_l_in;
tmp0_l += q2_l_in;
tmp0_l -= p5_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
tmp0_h = q7_h_in - q1_h_in;
tmp0_h += q2_h_in;
tmp0_h -= p5_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
__lsx_vst(filter8, dst, 16*9);
/* q3 */
tmp0_l = q7_l_in - q2_l_in;
tmp0_l += q3_l_in;
tmp0_l -= p4_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
tmp0_h = q7_h_in - q2_h_in;
tmp0_h += q3_h_in;
tmp0_h -= p4_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
q3 = __lsx_vbitsel_v(q3, out_l, flat2);
__lsx_vst(q3, dst, 16*10);
/* q4 */
tmp0_l = q7_l_in - q3_l_in;
tmp0_l += q4_l_in;
tmp0_l -= p3_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
tmp0_h = q7_h_in - q3_h_in;
tmp0_h += q4_h_in;
tmp0_h -= p3_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
q4 = __lsx_vbitsel_v(q4, out_l, flat2);
__lsx_vst(q4, dst, 16*11);
/* q5 */
tmp0_l = q7_l_in - q4_l_in;
tmp0_l += q5_l_in;
tmp0_l -= p2_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
tmp0_h = q7_h_in - q4_h_in;
tmp0_h += q5_h_in;
tmp0_h -= p2_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
q5 = __lsx_vbitsel_v(q5, out_l, flat2);
__lsx_vst(q5, dst, 16*12);
/* q6 */
tmp0_l = q7_l_in - q5_l_in;
tmp0_l += q6_l_in;
tmp0_l -= p1_l_in;
tmp1_l += tmp0_l;
out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
tmp0_h = q7_h_in - q5_h_in;
tmp0_h += q6_h_in;
tmp0_h -= p1_h_in;
tmp1_h += tmp0_h;
out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
out_l = __lsx_vpickev_b(out_h, out_l);
q6 = __lsx_vbitsel_v(q6, out_l, flat2);
__lsx_vst(q6, dst, 16*13);
return 0;
}
}
void ff_loop_filter_h_16_16_lsx(uint8_t *dst, ptrdiff_t stride,
int32_t b_limit_ptr,
int32_t limit_ptr,
int32_t thresh_ptr)
{
uint8_t early_exit = 0;
uint8_t transposed_input[16 * 24] __attribute__ ((aligned(16)));
uint8_t *filter48 = &transposed_input[16 * 16];
vp9_transpose_16x16((dst - 8), stride, &transposed_input[0], 16);
early_exit = vp9_vt_lpf_t4_and_t8_16w((transposed_input + 16 * 8),
&filter48[0], dst, stride,
b_limit_ptr, limit_ptr, thresh_ptr);
if (0 == early_exit) {
early_exit = vp9_vt_lpf_t16_16w((transposed_input + 16 * 8), dst,
stride, &filter48[0]);
if (0 == early_exit) {
vp9_transpose_16x16(transposed_input, 16, (dst - 8), stride);
}
}
}