/* * Copyright (c) 2021 Loongson Technology Corporation Limited * Contributed by Jin Bo * * 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); } } }