FFmpeg/libavcodec/loongarch/vp9_lpf_lsx.c

/*
 * Copyright (c) 2021 Loongson Technology Corporation Limited
 * Contributed by Jin Bo <jinbo@loongson.cn>
 *
 * This file is part of FFmpeg.
 *
 * FFmpeg is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public
 * License as published by the Free Software Foundation; either
 * version 2.1 of the License, or (at your option) any later version.
 *
 * FFmpeg is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public
 * License along with FFmpeg; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
 */

#include "libavcodec/vp9dsp.h"
#include "libavutil/loongarch/loongson_intrinsics.h"
#include "libavutil/common.h"
#include "vp9dsp_loongarch.h"

#define LSX_LD_8(_src, _stride, _stride2, _stride3, _stride4, _in0, _in1, _in2, \
                 _in3, _in4, _in5, _in6, _in7)                                  \
{                                                                               \
    _in0 = __lsx_vld(_src, 0);                                                  \
    _in1 = __lsx_vldx(_src, _stride);                                           \
    _in2 = __lsx_vldx(_src, _stride2);                                          \
    _in3 = __lsx_vldx(_src, _stride3);                                          \
    _src += _stride4;                                                           \
    _in4 = __lsx_vld(_src, 0);                                                  \
    _in5 = __lsx_vldx(_src, _stride);                                           \
    _in6 = __lsx_vldx(_src, _stride2);                                          \
    _in7 = __lsx_vldx(_src, _stride3);                                          \
}

#define LSX_ST_8(_dst0, _dst1, _dst2, _dst3, _dst4, _dst5, _dst6, _dst7,        \
                 _dst, _stride, _stride2, _stride3, _stride4)                   \
{                                                                               \
    __lsx_vst(_dst0, _dst, 0);                                                  \
    __lsx_vstx(_dst1, _dst, _stride);                                           \
    __lsx_vstx(_dst2, _dst, _stride2);                                          \
    __lsx_vstx(_dst3, _dst, _stride3);                                          \
    _dst += _stride4;                                                           \
    __lsx_vst(_dst4, _dst, 0);                                                  \
    __lsx_vstx(_dst5, _dst, _stride);                                           \
    __lsx_vstx(_dst6, _dst, _stride2);                                          \
    __lsx_vstx(_dst7, _dst, _stride3);                                          \
}

#define VP9_LPF_FILTER4_4W(p1_src, p0_src, q0_src, q1_src, mask_src, hev_src, \
                           p1_dst, p0_dst, q0_dst, q1_dst)                    \
{                                                                             \
    __m128i p1_tmp, p0_tmp, q0_tmp, q1_tmp, q0_sub_p0, filt, filt1, filt2;    \
    const __m128i cnst3b = __lsx_vldi(3);                                     \
    const __m128i cnst4b = __lsx_vldi(4);                                     \
                                                                              \
    p1_tmp = __lsx_vxori_b(p1_src, 0x80);                                     \
    p0_tmp = __lsx_vxori_b(p0_src, 0x80);                                     \
    q0_tmp = __lsx_vxori_b(q0_src, 0x80);                                     \
    q1_tmp = __lsx_vxori_b(q1_src, 0x80);                                     \
                                                                              \
    filt = __lsx_vssub_b(p1_tmp, q1_tmp);                                     \
                                                                              \
    filt = filt & hev_src;                                                    \
                                                                              \
    q0_sub_p0 = __lsx_vssub_b(q0_tmp, p0_tmp);                                \
    filt = __lsx_vsadd_b(filt, q0_sub_p0);                                    \
    filt = __lsx_vsadd_b(filt, q0_sub_p0);                                    \
    filt = __lsx_vsadd_b(filt, q0_sub_p0);                                    \
    filt = filt & mask_src;                                                   \
                                                                              \
    filt1 = __lsx_vsadd_b(filt, cnst4b);                                      \
    filt1 = __lsx_vsrai_b(filt1, 3);                                          \
                                                                              \
    filt2 = __lsx_vsadd_b(filt, cnst3b);                                      \
    filt2 = __lsx_vsrai_b(filt2, 3);                                          \
                                                                              \
    q0_tmp = __lsx_vssub_b(q0_tmp, filt1);                                    \
    q0_dst = __lsx_vxori_b(q0_tmp, 0x80);                                     \
    p0_tmp = __lsx_vsadd_b(p0_tmp, filt2);                                    \
    p0_dst = __lsx_vxori_b(p0_tmp, 0x80);                                     \
                                                                              \
    filt = __lsx_vsrari_b(filt1, 1);                                          \
    hev_src = __lsx_vxori_b(hev_src, 0xff);                                   \
    filt = filt & hev_src;                                                    \
                                                                              \
    q1_tmp = __lsx_vssub_b(q1_tmp, filt);                                     \
    q1_dst = __lsx_vxori_b(q1_tmp, 0x80);                                     \
    p1_tmp = __lsx_vsadd_b(p1_tmp, filt);                                     \
    p1_dst = __lsx_vxori_b(p1_tmp, 0x80);                                     \
}

#define VP9_FLAT4(p3_src, p2_src, p0_src, q0_src, q2_src, q3_src, flat_dst)  \
{                                                                            \
    __m128i f_tmp = __lsx_vldi(1);                                           \
    __m128i p2_a_sub_p0, q2_a_sub_q0, p3_a_sub_p0, q3_a_sub_q0;              \
                                                                             \
    p2_a_sub_p0 = __lsx_vabsd_bu(p2_src, p0_src);                            \
    q2_a_sub_q0 = __lsx_vabsd_bu(q2_src, q0_src);                            \
    p3_a_sub_p0 = __lsx_vabsd_bu(p3_src, p0_src);                            \
    q3_a_sub_q0 = __lsx_vabsd_bu(q3_src, q0_src);                            \
                                                                             \
    p2_a_sub_p0 = __lsx_vmax_bu(p2_a_sub_p0, q2_a_sub_q0);                   \
    flat_dst = __lsx_vmax_bu(p2_a_sub_p0, flat_dst);                         \
    p3_a_sub_p0 = __lsx_vmax_bu(p3_a_sub_p0, q3_a_sub_q0);                   \
    flat_dst = __lsx_vmax_bu(p3_a_sub_p0, flat_dst);                         \
                                                                             \
    flat_dst = __lsx_vslt_bu(f_tmp, flat_dst);                               \
    flat_dst = __lsx_vxori_b(flat_dst, 0xff);                                \
    flat_dst = flat_dst & mask;                                              \
}

#define VP9_FLAT5(p7_src, p6_src, p5_src, p4_src, p0_src, q0_src, q4_src, \
                  q5_src, q6_src, q7_src, flat_src, flat2_dst)            \
{                                                                         \
    __m128i f_tmp = __lsx_vldi(1);                                        \
    __m128i p4_a_sub_p0, q4_a_sub_q0, p5_a_sub_p0, q5_a_sub_q0;           \
    __m128i p6_a_sub_p0, q6_a_sub_q0, p7_a_sub_p0, q7_a_sub_q0;           \
                                                                          \
    p4_a_sub_p0 = __lsx_vabsd_bu(p4_src, p0_src);                         \
    q4_a_sub_q0 = __lsx_vabsd_bu(q4_src, q0_src);                         \
    p5_a_sub_p0 = __lsx_vabsd_bu(p5_src, p0_src);                         \
    q5_a_sub_q0 = __lsx_vabsd_bu(q5_src, q0_src);                         \
    p6_a_sub_p0 = __lsx_vabsd_bu(p6_src, p0_src);                         \
    q6_a_sub_q0 = __lsx_vabsd_bu(q6_src, q0_src);                         \
    p7_a_sub_p0 = __lsx_vabsd_bu(p7_src, p0_src);                         \
    q7_a_sub_q0 = __lsx_vabsd_bu(q7_src, q0_src);                         \
                                                                          \
    p4_a_sub_p0 = __lsx_vmax_bu(p4_a_sub_p0, q4_a_sub_q0);                \
    flat2_dst = __lsx_vmax_bu(p5_a_sub_p0, q5_a_sub_q0);                  \
    flat2_dst = __lsx_vmax_bu(p4_a_sub_p0, flat2_dst);                    \
    p6_a_sub_p0 = __lsx_vmax_bu(p6_a_sub_p0, q6_a_sub_q0);                \
    flat2_dst = __lsx_vmax_bu(p6_a_sub_p0, flat2_dst);                    \
    p7_a_sub_p0 = __lsx_vmax_bu(p7_a_sub_p0, q7_a_sub_q0);                \
    flat2_dst = __lsx_vmax_bu(p7_a_sub_p0, flat2_dst);                    \
                                                                          \
    flat2_dst = __lsx_vslt_bu(f_tmp, flat2_dst);                          \
    flat2_dst = __lsx_vxori_b(flat2_dst, 0xff);                           \
    flat2_dst = flat2_dst & flat_src;                                     \
}

#define VP9_FILTER8(p3_src, p2_src, p1_src, p0_src,            \
                    q0_src, q1_src, q2_src, q3_src,            \
                    p2_filt8_dst, p1_filt8_dst, p0_filt8_dst,  \
                    q0_filt8_dst, q1_filt8_dst, q2_filt8_dst)  \
{                                                              \
    __m128i tmp0, tmp1, tmp2;                                  \
                                                               \
    tmp2 = __lsx_vadd_h(p2_src, p1_src);                       \
    tmp2 = __lsx_vadd_h(tmp2, p0_src);                         \
    tmp0 = __lsx_vslli_h(p3_src, 1);                           \
                                                               \
    tmp0 = __lsx_vadd_h(tmp0, tmp2);                           \
    tmp0 = __lsx_vadd_h(tmp0, q0_src);                         \
    tmp1 = __lsx_vadd_h(tmp0, p3_src);                         \
    tmp1 = __lsx_vadd_h(tmp1, p2_src);                         \
    p2_filt8_dst = __lsx_vsrari_h(tmp1, 3);                    \
                                                               \
    tmp1 = __lsx_vadd_h(tmp0, p1_src);                         \
    tmp1 = __lsx_vadd_h(tmp1, q1_src);                         \
    p1_filt8_dst = __lsx_vsrari_h(tmp1, 3);                    \
                                                               \
    tmp1 = __lsx_vadd_h(q2_src, q1_src);                       \
    tmp1 = __lsx_vadd_h(tmp1, q0_src);                         \
    tmp2 = __lsx_vadd_h(tmp2, tmp1);                           \
    tmp0 = __lsx_vadd_h(tmp2, p0_src);                         \
    tmp0 = __lsx_vadd_h(tmp0, p3_src);                         \
    p0_filt8_dst = __lsx_vsrari_h(tmp0, 3);                    \
                                                               \
    tmp0 = __lsx_vadd_h(q2_src, q3_src);                       \
    tmp0 = __lsx_vadd_h(tmp0, p0_src);                         \
    tmp0 = __lsx_vadd_h(tmp0, tmp1);                           \
    tmp1 = __lsx_vadd_h(q3_src, q3_src);                       \
    tmp1 = __lsx_vadd_h(tmp1, tmp0);                           \
    q2_filt8_dst = __lsx_vsrari_h(tmp1, 3);                    \
                                                               \
    tmp0 = __lsx_vadd_h(tmp2, q3_src);                         \
    tmp1 = __lsx_vadd_h(tmp0, q0_src);                         \
    q0_filt8_dst = __lsx_vsrari_h(tmp1, 3);                    \
                                                               \
    tmp1 = __lsx_vsub_h(tmp0, p2_src);                         \
    tmp0 = __lsx_vadd_h(q1_src, q3_src);                       \
    tmp1 = __lsx_vadd_h(tmp0, tmp1);                           \
    q1_filt8_dst = __lsx_vsrari_h(tmp1, 3);                    \
}

#define LPF_MASK_HEV(p3_src, p2_src, p1_src, p0_src, q0_src, q1_src,        \
                     q2_src, q3_src, limit_src, b_limit_src, thresh_src,    \
                     hev_dst, mask_dst, flat_dst)                           \
{                                                                           \
    __m128i p3_asub_p2_tmp, p2_asub_p1_tmp, p1_asub_p0_tmp, q1_asub_q0_tmp; \
    __m128i p1_asub_q1_tmp, p0_asub_q0_tmp, q3_asub_q2_tmp, q2_asub_q1_tmp; \
                                                                            \
    /* absolute subtraction of pixel values */                              \
    p3_asub_p2_tmp = __lsx_vabsd_bu(p3_src, p2_src);                        \
    p2_asub_p1_tmp = __lsx_vabsd_bu(p2_src, p1_src);                        \
    p1_asub_p0_tmp = __lsx_vabsd_bu(p1_src, p0_src);                        \
    q1_asub_q0_tmp = __lsx_vabsd_bu(q1_src, q0_src);                        \
    q2_asub_q1_tmp = __lsx_vabsd_bu(q2_src, q1_src);                        \
    q3_asub_q2_tmp = __lsx_vabsd_bu(q3_src, q2_src);                        \
    p0_asub_q0_tmp = __lsx_vabsd_bu(p0_src, q0_src);                        \
    p1_asub_q1_tmp = __lsx_vabsd_bu(p1_src, q1_src);                        \
                                                                            \
    /* calculation of hev */                                                \
    flat_dst = __lsx_vmax_bu(p1_asub_p0_tmp, q1_asub_q0_tmp);               \
    hev_dst = __lsx_vslt_bu(thresh_src, flat_dst);                          \
                                                                            \
    /* calculation of mask */                                               \
    p0_asub_q0_tmp = __lsx_vsadd_bu(p0_asub_q0_tmp, p0_asub_q0_tmp);        \
    p1_asub_q1_tmp = __lsx_vsrli_b(p1_asub_q1_tmp, 1);                      \
    p0_asub_q0_tmp = __lsx_vsadd_bu(p0_asub_q0_tmp, p1_asub_q1_tmp);        \
                                                                            \
    mask_dst = __lsx_vslt_bu(b_limit_src, p0_asub_q0_tmp);                  \
    mask_dst = __lsx_vmax_bu(flat_dst, mask_dst);                           \
    p3_asub_p2_tmp = __lsx_vmax_bu(p3_asub_p2_tmp, p2_asub_p1_tmp);         \
    mask_dst = __lsx_vmax_bu(p3_asub_p2_tmp, mask_dst);                     \
    q2_asub_q1_tmp = __lsx_vmax_bu(q2_asub_q1_tmp, q3_asub_q2_tmp);         \
    mask_dst = __lsx_vmax_bu(q2_asub_q1_tmp, mask_dst);                     \
                                                                            \
    mask_dst = __lsx_vslt_bu(limit_src, mask_dst);                          \
    mask_dst = __lsx_vxori_b(mask_dst, 0xff);                               \
}

void ff_loop_filter_v_4_8_lsx(uint8_t *dst, ptrdiff_t stride,
                              int32_t b_limit_ptr,
                              int32_t limit_ptr,
                              int32_t thresh_ptr)
{
    ptrdiff_t stride2 = stride << 1;
    ptrdiff_t stride3 = stride2 + stride;
    ptrdiff_t stride4 = stride2 << 1;
    __m128i mask, hev, flat, thresh, b_limit, limit;
    __m128i p3, p2, p1, p0, q3, q2, q1, q0, p1_out, p0_out, q0_out, q1_out;

    DUP4_ARG2(__lsx_vldx, dst, -stride4, dst, -stride3, dst, -stride2,
              dst, -stride, p3, p2, p1, p0);
    q0 = __lsx_vld(dst, 0);
    DUP2_ARG2(__lsx_vldx, dst, stride, dst, stride2, q1, q2);
    q3 = __lsx_vldx(dst, stride3);

    thresh  = __lsx_vreplgr2vr_b(thresh_ptr);
    b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
    limit   = __lsx_vreplgr2vr_b(limit_ptr);

    LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
                 hev, mask, flat);

    VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
                       q1_out);

    __lsx_vstelm_d(p1_out, dst - stride2, 0, 0);
    __lsx_vstelm_d(p0_out, dst -  stride, 0, 0);
    __lsx_vstelm_d(q0_out, dst          , 0, 0);
    __lsx_vstelm_d(q1_out, dst +  stride, 0, 0);
}

void ff_loop_filter_v_44_16_lsx(uint8_t *dst, ptrdiff_t stride,
                                int32_t b_limit_ptr,
                                int32_t limit_ptr,
                                int32_t thresh_ptr)
{
    ptrdiff_t stride2 = stride << 1;
    ptrdiff_t stride3 = stride2 + stride;
    ptrdiff_t stride4 = stride2 << 1;
    __m128i mask, hev, flat, thresh0, b_limit0;
    __m128i limit0, thresh1, b_limit1, limit1;
    __m128i p3, p2, p1, p0, q3, q2, q1, q0;

    DUP4_ARG2(__lsx_vldx, dst, -stride4, dst, -stride3, dst, -stride2,
              dst, -stride, p3, p2, p1, p0);
    q0 = __lsx_vld(dst, 0);
    DUP2_ARG2(__lsx_vldx, dst, stride, dst, stride2, q1, q2);
    q3 = __lsx_vldx(dst, stride3);

    thresh0 = __lsx_vreplgr2vr_b(thresh_ptr);
    thresh1 = __lsx_vreplgr2vr_b(thresh_ptr >> 8);
    thresh0 = __lsx_vilvl_d(thresh1, thresh0);

    b_limit0 = __lsx_vreplgr2vr_b(b_limit_ptr);
    b_limit1 = __lsx_vreplgr2vr_b(b_limit_ptr >> 8);
    b_limit0 = __lsx_vilvl_d(b_limit1, b_limit0);

    limit0 = __lsx_vreplgr2vr_b(limit_ptr);
    limit1 = __lsx_vreplgr2vr_b(limit_ptr >> 8);
    limit0 = __lsx_vilvl_d(limit1, limit0);

    LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit0, b_limit0, thresh0,
                 hev, mask, flat);
    VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1, p0, q0, q1);

    __lsx_vst(p1, dst - stride2, 0);
    __lsx_vst(p0, dst -  stride, 0);
    __lsx_vst(q0, dst          , 0);
    __lsx_vst(q1, dst +  stride, 0);
}

void ff_loop_filter_v_8_8_lsx(uint8_t *dst, ptrdiff_t stride,
                              int32_t b_limit_ptr,
                              int32_t limit_ptr,
                              int32_t thresh_ptr)
{
    ptrdiff_t stride2 = stride << 1;
    ptrdiff_t stride3 = stride2 + stride;
    ptrdiff_t stride4 = stride2 << 1;
    __m128i mask, hev, flat, thresh, b_limit, limit;
    __m128i p3, p2, p1, p0, q3, q2, q1, q0;
    __m128i p2_out, p1_out, p0_out, q0_out, q1_out, q2_out;
    __m128i p2_filter8, p1_filter8, p0_filter8;
    __m128i q0_filter8, q1_filter8, q2_filter8;
    __m128i p3_l, p2_l, p1_l, p0_l, q3_l, q2_l, q1_l, q0_l;
    __m128i zero = __lsx_vldi(0);

    DUP4_ARG2(__lsx_vldx, dst, -stride4, dst, -stride3, dst, -stride2,
              dst, -stride, p3, p2, p1, p0);
    q0 = __lsx_vld(dst, 0);
    DUP2_ARG2(__lsx_vldx, dst, stride, dst, stride2, q1, q2);
    q3 = __lsx_vldx(dst, stride3);

    thresh  = __lsx_vreplgr2vr_b(thresh_ptr);
    b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
    limit   = __lsx_vreplgr2vr_b(limit_ptr);

    LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
                 hev, mask, flat);
    VP9_FLAT4(p3, p2, p0, q0, q2, q3, flat);
    VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
                       q1_out);

    flat = __lsx_vilvl_d(zero, flat);

    /* if flat is zero for all pixels, then no need to calculate other filter */
    if (__lsx_bz_v(flat)) {
        __lsx_vstelm_d(p1_out, dst - stride2, 0, 0);
        __lsx_vstelm_d(p0_out, dst -  stride, 0, 0);
        __lsx_vstelm_d(q0_out, dst          , 0, 0);
        __lsx_vstelm_d(q1_out, dst +  stride, 0, 0);
    } else {
        DUP4_ARG2(__lsx_vilvl_b, zero, p3, zero, p2, zero, p1, zero, p0,
                  p3_l, p2_l, p1_l, p0_l);
        DUP4_ARG2(__lsx_vilvl_b, zero, q0, zero, q1, zero, q2, zero, q3,
                  q0_l, q1_l, q2_l, q3_l);
        VP9_FILTER8(p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l, p2_filter8,
                    p1_filter8, p0_filter8, q0_filter8, q1_filter8, q2_filter8);

        /* convert 16 bit output data into 8 bit */
        DUP4_ARG2(__lsx_vpickev_b, zero, p2_filter8, zero, p1_filter8,
                  zero, p0_filter8, zero, q0_filter8, p2_filter8,
                  p1_filter8, p0_filter8, q0_filter8);
        DUP2_ARG2(__lsx_vpickev_b, zero, q1_filter8, zero, q2_filter8,
                  q1_filter8, q2_filter8);

        /* store pixel values */
        p2_out = __lsx_vbitsel_v(p2, p2_filter8, flat);
        p1_out = __lsx_vbitsel_v(p1_out, p1_filter8, flat);
        p0_out = __lsx_vbitsel_v(p0_out, p0_filter8, flat);
        q0_out = __lsx_vbitsel_v(q0_out, q0_filter8, flat);
        q1_out = __lsx_vbitsel_v(q1_out, q1_filter8, flat);
        q2_out = __lsx_vbitsel_v(q2, q2_filter8, flat);

        __lsx_vstelm_d(p2_out, dst - stride3, 0, 0);
        __lsx_vstelm_d(p1_out, dst - stride2, 0, 0);
        __lsx_vstelm_d(p0_out, dst - stride, 0, 0);
        __lsx_vstelm_d(q0_out, dst, 0, 0);
        __lsx_vstelm_d(q1_out, dst + stride, 0, 0);
        __lsx_vstelm_d(q2_out, dst + stride2, 0, 0);
    }
}

void ff_loop_filter_v_88_16_lsx(uint8_t *dst, ptrdiff_t stride,
                                int32_t b_limit_ptr,
                                int32_t limit_ptr,
                                int32_t thresh_ptr)
{
    ptrdiff_t stride2 = stride << 1;
    ptrdiff_t stride3 = stride2 + stride;
    ptrdiff_t stride4 = stride2 << 1;
    __m128i p3, p2, p1, p0, q3, q2, q1, q0;
    __m128i p2_out, p1_out, p0_out, q0_out, q1_out, q2_out;
    __m128i flat, mask, hev, tmp, thresh, b_limit, limit;
    __m128i p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l;
    __m128i p3_h, p2_h, p1_h, p0_h, q0_h, q1_h, q2_h, q3_h;
    __m128i p2_filt8_l, p1_filt8_l, p0_filt8_l;
    __m128i q0_filt8_l, q1_filt8_l, q2_filt8_l;
    __m128i p2_filt8_h, p1_filt8_h, p0_filt8_h;
    __m128i q0_filt8_h, q1_filt8_h, q2_filt8_h;
    __m128i zero = __lsx_vldi(0);

    /* load vector elements */
    DUP4_ARG2(__lsx_vldx, dst, -stride4, dst, -stride3, dst, -stride2,
              dst, -stride, p3, p2, p1, p0);
    q0 = __lsx_vld(dst, 0);
    DUP2_ARG2(__lsx_vldx, dst, stride, dst, stride2, q1, q2);
    q3 = __lsx_vldx(dst, stride3);

    thresh = __lsx_vreplgr2vr_b(thresh_ptr);
    tmp    = __lsx_vreplgr2vr_b(thresh_ptr >> 8);
    thresh = __lsx_vilvl_d(tmp, thresh);

    b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
    tmp     = __lsx_vreplgr2vr_b(b_limit_ptr >> 8);
    b_limit = __lsx_vilvl_d(tmp, b_limit);

    limit = __lsx_vreplgr2vr_b(limit_ptr);
    tmp   = __lsx_vreplgr2vr_b(limit_ptr >> 8);
    limit = __lsx_vilvl_d(tmp, limit);

    /* mask and hev */
    LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
                 hev, mask, flat);
    VP9_FLAT4(p3, p2, p0, q0, q2, q3, flat);
    VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
                       q1_out);

    /* if flat is zero for all pixels, then no need to calculate other filter */
    if (__lsx_bz_v(flat)) {
        __lsx_vst(p1_out, dst - stride2, 0);
        __lsx_vst(p0_out, dst - stride, 0);
        __lsx_vst(q0_out, dst, 0);
        __lsx_vst(q1_out, dst + stride, 0);
    } else {
        DUP4_ARG2(__lsx_vilvl_b, zero, p3, zero, p2, zero, p1, zero, p0,
                  p3_l, p2_l, p1_l, p0_l);
        DUP4_ARG2(__lsx_vilvl_b, zero, q0, zero, q1, zero, q2, zero, q3,
                  q0_l, q1_l, q2_l, q3_l);
        VP9_FILTER8(p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l, p2_filt8_l,
                    p1_filt8_l, p0_filt8_l, q0_filt8_l, q1_filt8_l, q2_filt8_l);

        DUP4_ARG2(__lsx_vilvh_b, zero, p3, zero, p2, zero, p1, zero, p0,
                  p3_h, p2_h, p1_h, p0_h);
        DUP4_ARG2(__lsx_vilvh_b, zero, q0, zero, q1, zero, q2, zero, q3,
                  q0_h, q1_h, q2_h, q3_h);
        VP9_FILTER8(p3_h, p2_h, p1_h, p0_h, q0_h, q1_h, q2_h, q3_h, p2_filt8_h,
                    p1_filt8_h, p0_filt8_h, q0_filt8_h, q1_filt8_h, q2_filt8_h);

        /* convert 16 bit output data into 8 bit */
        DUP4_ARG2(__lsx_vpickev_b, p2_filt8_h, p2_filt8_l, p1_filt8_h,
                  p1_filt8_l, p0_filt8_h, p0_filt8_l, q0_filt8_h, q0_filt8_l,
                  p2_filt8_l, p1_filt8_l, p0_filt8_l, q0_filt8_l);
        DUP2_ARG2(__lsx_vpickev_b, q1_filt8_h, q1_filt8_l, q2_filt8_h,
                  q2_filt8_l, q1_filt8_l, q2_filt8_l);

        /* store pixel values */
        p2_out = __lsx_vbitsel_v(p2, p2_filt8_l, flat);
        p1_out = __lsx_vbitsel_v(p1_out, p1_filt8_l, flat);
        p0_out = __lsx_vbitsel_v(p0_out, p0_filt8_l, flat);
        q0_out = __lsx_vbitsel_v(q0_out, q0_filt8_l, flat);
        q1_out = __lsx_vbitsel_v(q1_out, q1_filt8_l, flat);
        q2_out = __lsx_vbitsel_v(q2, q2_filt8_l, flat);


        __lsx_vstx(p2_out, dst, -stride3);
        __lsx_vstx(p1_out, dst, -stride2);
        __lsx_vstx(p0_out, dst, -stride);
        __lsx_vst(q0_out, dst, 0);
        __lsx_vstx(q1_out, dst, stride);
        __lsx_vstx(q2_out, dst, stride2);
    }
}

void ff_loop_filter_v_84_16_lsx(uint8_t *dst, ptrdiff_t stride,
                                int32_t b_limit_ptr,
                                int32_t limit_ptr,
                                int32_t thresh_ptr)
{
    ptrdiff_t stride2 = stride << 1;
    ptrdiff_t stride3 = stride2 + stride;
    ptrdiff_t stride4 = stride2 << 1;
    __m128i p3, p2, p1, p0, q3, q2, q1, q0;
    __m128i p2_out, p1_out, p0_out, q0_out, q1_out, q2_out;
    __m128i flat, mask, hev, tmp, thresh, b_limit, limit;
    __m128i p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l;
    __m128i p2_filt8_l, p1_filt8_l, p0_filt8_l;
    __m128i q0_filt8_l, q1_filt8_l, q2_filt8_l;
    __m128i zero = __lsx_vldi(0);

    /* load vector elements */
    DUP4_ARG2(__lsx_vldx, dst, -stride4, dst, -stride3, dst, -stride2,
              dst, -stride, p3, p2, p1, p0);
    q0 = __lsx_vld(dst, 0);
    DUP2_ARG2(__lsx_vldx, dst, stride, dst, stride2, q1, q2);
    q3 = __lsx_vldx(dst, stride3);

    thresh = __lsx_vreplgr2vr_b(thresh_ptr);
    tmp    = __lsx_vreplgr2vr_b(thresh_ptr >> 8);
    thresh = __lsx_vilvl_d(tmp, thresh);

    b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
    tmp     = __lsx_vreplgr2vr_b(b_limit_ptr >> 8);
    b_limit = __lsx_vilvl_d(tmp, b_limit);

    limit = __lsx_vreplgr2vr_b(limit_ptr);
    tmp   = __lsx_vreplgr2vr_b(limit_ptr >> 8);
    limit = __lsx_vilvl_d(tmp, limit);

    /* mask and hev */
    LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
                 hev, mask, flat);
    VP9_FLAT4(p3, p2, p0, q0, q2, q3, flat);
    VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
                       q1_out);

    flat = __lsx_vilvl_d(zero, flat);

    /* if flat is zero for all pixels, then no need to calculate other filter */
    if (__lsx_bz_v(flat)) {
        __lsx_vstx(p1_out, dst, -stride2);
        __lsx_vstx(p0_out, dst, -stride);
        __lsx_vst(q0_out, dst, 0);
        __lsx_vstx(q1_out, dst, stride);
    } else {
        DUP4_ARG2(__lsx_vilvl_b, zero, p3, zero, p2, zero, p1, zero, p0,
                  p3_l, p2_l, p1_l, p0_l);
        DUP4_ARG2(__lsx_vilvl_b, zero, q0, zero, q1, zero, q2, zero, q3,
                  q0_l, q1_l, q2_l, q3_l);
        VP9_FILTER8(p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l, p2_filt8_l,
                    p1_filt8_l, p0_filt8_l, q0_filt8_l, q1_filt8_l, q2_filt8_l);

        /* convert 16 bit output data into 8 bit */
        DUP4_ARG2(__lsx_vpickev_b, p2_filt8_l, p2_filt8_l, p1_filt8_l,
                  p1_filt8_l, p0_filt8_l, p0_filt8_l, q0_filt8_l, q0_filt8_l,
                  p2_filt8_l, p1_filt8_l, p0_filt8_l, q0_filt8_l);
        DUP2_ARG2(__lsx_vpickev_b, q1_filt8_l, q1_filt8_l, q2_filt8_l,
                  q2_filt8_l, q1_filt8_l, q2_filt8_l);

        /* store pixel values */
        p2_out = __lsx_vbitsel_v(p2, p2_filt8_l, flat);
        p1_out = __lsx_vbitsel_v(p1_out, p1_filt8_l, flat);
        p0_out = __lsx_vbitsel_v(p0_out, p0_filt8_l, flat);
        q0_out = __lsx_vbitsel_v(q0_out, q0_filt8_l, flat);
        q1_out = __lsx_vbitsel_v(q1_out, q1_filt8_l, flat);
        q2_out = __lsx_vbitsel_v(q2, q2_filt8_l, flat);

        __lsx_vstx(p2_out, dst, -stride3);
        __lsx_vstx(p1_out, dst, -stride2);
        __lsx_vstx(p0_out, dst, -stride);
        __lsx_vst(q0_out, dst, 0);
        __lsx_vstx(q1_out, dst, stride);
        __lsx_vstx(q2_out, dst, stride2);
    }
}

void ff_loop_filter_v_48_16_lsx(uint8_t *dst, ptrdiff_t stride,
                                int32_t b_limit_ptr,
                                int32_t limit_ptr,
                                int32_t thresh_ptr)
{
    ptrdiff_t stride2 = stride << 1;
    ptrdiff_t stride3 = stride2 + stride;
    ptrdiff_t stride4 = stride2 << 1;
    __m128i p3, p2, p1, p0, q3, q2, q1, q0;
    __m128i p2_out, p1_out, p0_out, q0_out, q1_out, q2_out;
    __m128i flat, mask, hev, tmp, thresh, b_limit, limit;
    __m128i p3_h, p2_h, p1_h, p0_h, q0_h, q1_h, q2_h, q3_h;
    __m128i p2_filt8_h, p1_filt8_h, p0_filt8_h;
    __m128i q0_filt8_h, q1_filt8_h, q2_filt8_h;
    __m128i zero = { 0 };

    /* load vector elements */
    DUP4_ARG2(__lsx_vldx, dst, -stride4, dst, -stride3, dst, -stride2,
              dst, -stride, p3, p2, p1, p0);
    q0 = __lsx_vld(dst, 0);
    DUP2_ARG2(__lsx_vldx, dst, stride, dst, stride2, q1, q2);
    q3 = __lsx_vldx(dst, stride3);

    thresh = __lsx_vreplgr2vr_b(thresh_ptr);
    tmp    = __lsx_vreplgr2vr_b(thresh_ptr >> 8);
    thresh = __lsx_vilvl_d(tmp, thresh);

    b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
    tmp     = __lsx_vreplgr2vr_b(b_limit_ptr >> 8);
    b_limit = __lsx_vilvl_d(tmp, b_limit);

    limit = __lsx_vreplgr2vr_b(limit_ptr);
    tmp   = __lsx_vreplgr2vr_b(limit_ptr >> 8);
    limit = __lsx_vilvl_d(tmp, limit);

    /* mask and hev */
    LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
                 hev, mask, flat);
    VP9_FLAT4(p3, p2, p0, q0, q2, q3, flat);
    VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
                       q1_out);

    flat = __lsx_vilvh_d(flat, zero);

    /* if flat is zero for all pixels, then no need to calculate other filter */
    if (__lsx_bz_v(flat)) {
        __lsx_vstx(p1_out, dst, -stride2);
        __lsx_vstx(p0_out, dst, -stride);
        __lsx_vst(q0_out, dst, 0);
        __lsx_vstx(q1_out, dst, stride);
    } else {
        DUP4_ARG2(__lsx_vilvh_b, zero, p3, zero, p2, zero, p1, zero, p0,
                  p3_h, p2_h, p1_h, p0_h);
        DUP4_ARG2(__lsx_vilvh_b, zero, q0, zero, q1, zero, q2, zero, q3,
                  q0_h, q1_h, q2_h, q3_h);
        VP9_FILTER8(p3_h, p2_h, p1_h, p0_h, q0_h, q1_h, q2_h, q3_h, p2_filt8_h,
                    p1_filt8_h, p0_filt8_h, q0_filt8_h, q1_filt8_h, q2_filt8_h);

        /* convert 16 bit output data into 8 bit */
        DUP4_ARG2(__lsx_vpickev_b, p2_filt8_h, p2_filt8_h, p1_filt8_h,
                  p1_filt8_h, p0_filt8_h, p0_filt8_h, q0_filt8_h, q0_filt8_h,
                  p2_filt8_h, p1_filt8_h, p0_filt8_h, q0_filt8_h);
        DUP2_ARG2(__lsx_vpickev_b, q1_filt8_h, q1_filt8_h, q2_filt8_h,
                  q2_filt8_h, q1_filt8_h, q2_filt8_h);

        /* store pixel values */
        p2_out = __lsx_vbitsel_v(p2, p2_filt8_h, flat);
        p1_out = __lsx_vbitsel_v(p1_out, p1_filt8_h, flat);
        p0_out = __lsx_vbitsel_v(p0_out, p0_filt8_h, flat);
        q0_out = __lsx_vbitsel_v(q0_out, q0_filt8_h, flat);
        q1_out = __lsx_vbitsel_v(q1_out, q1_filt8_h, flat);
        q2_out = __lsx_vbitsel_v(q2, q2_filt8_h, flat);

        __lsx_vstx(p2_out, dst, -stride3);
        __lsx_vstx(p1_out, dst, -stride2);
        __lsx_vstx(p0_out, dst, -stride);
        __lsx_vst(q0_out, dst, 0);
        __lsx_vstx(q1_out, dst, stride);
        __lsx_vstx(q2_out, dst, stride2);
    }
}

static int32_t vp9_hz_lpf_t4_and_t8_16w(uint8_t *dst, ptrdiff_t stride,
                                        uint8_t *filter48,
                                        int32_t b_limit_ptr,
                                        int32_t limit_ptr,
                                        int32_t thresh_ptr)
{
    ptrdiff_t stride2 = stride << 1;
    ptrdiff_t stride3 = stride2 + stride;
    ptrdiff_t stride4 = stride2 << 1;
    __m128i p3, p2, p1, p0, q3, q2, q1, q0;
    __m128i p2_out, p1_out, p0_out, q0_out, q1_out, q2_out;
    __m128i flat, mask, hev, thresh, b_limit, limit;
    __m128i p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l;
    __m128i p3_h, p2_h, p1_h, p0_h, q0_h, q1_h, q2_h, q3_h;
    __m128i p2_filt8_l, p1_filt8_l, p0_filt8_l;
    __m128i q0_filt8_l, q1_filt8_l, q2_filt8_l;
    __m128i p2_filt8_h, p1_filt8_h, p0_filt8_h;
    __m128i q0_filt8_h, q1_filt8_h, q2_filt8_h;
    __m128i zero = __lsx_vldi(0);

    /* load vector elements */
    DUP4_ARG2(__lsx_vldx, dst, -stride4, dst, -stride3, dst, -stride2,
              dst, -stride, p3, p2, p1, p0);
    q0 = __lsx_vld(dst, 0);
    DUP2_ARG2(__lsx_vldx, dst, stride, dst, stride2, q1, q2);
    q3 = __lsx_vldx(dst, stride3);

    thresh  = __lsx_vreplgr2vr_b(thresh_ptr);
    b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
    limit   = __lsx_vreplgr2vr_b(limit_ptr);

    /* mask and hev */
    LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
                 hev, mask, flat);
    VP9_FLAT4(p3, p2, p0, q0, q2, q3, flat);
    VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
                       q1_out);

    /* if flat is zero for all pixels, then no need to calculate other filter */
    if (__lsx_bz_v(flat)) {
        __lsx_vstx(p1_out, dst, -stride2);
        __lsx_vstx(p0_out, dst, -stride);
        __lsx_vst(q0_out, dst, 0);
        __lsx_vstx(q1_out, dst, stride);
        return 1;
    } else {
        DUP4_ARG2(__lsx_vilvl_b, zero, p3, zero, p2, zero, p1, zero, p0,
                  p3_l, p2_l, p1_l, p0_l);
        DUP4_ARG2(__lsx_vilvl_b, zero, q0, zero, q1, zero, q2, zero, q3,
                  q0_l, q1_l, q2_l, q3_l);
        VP9_FILTER8(p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l, p2_filt8_l,
                    p1_filt8_l, p0_filt8_l, q0_filt8_l, q1_filt8_l, q2_filt8_l);

        DUP4_ARG2(__lsx_vilvh_b, zero, p3, zero, p2, zero, p1, zero, p0,
                  p3_h, p2_h, p1_h, p0_h);
        DUP4_ARG2(__lsx_vilvh_b, zero, q0, zero, q1, zero, q2, zero, q3,
                  q0_h, q1_h, q2_h, q3_h);
        VP9_FILTER8(p3_h, p2_h, p1_h, p0_h, q0_h, q1_h, q2_h, q3_h, p2_filt8_h,
                    p1_filt8_h, p0_filt8_h, q0_filt8_h, q1_filt8_h, q2_filt8_h);

        /* convert 16 bit output data into 8 bit */
        DUP4_ARG2(__lsx_vpickev_b, p2_filt8_h, p2_filt8_l, p1_filt8_h,
                  p1_filt8_l, p0_filt8_h, p0_filt8_l, q0_filt8_h, q0_filt8_l,
                  p2_filt8_l, p1_filt8_l, p0_filt8_l, q0_filt8_l);
        DUP2_ARG2(__lsx_vpickev_b, q1_filt8_h, q1_filt8_l, q2_filt8_h,
                  q2_filt8_l, q1_filt8_l, q2_filt8_l);

        /* store pixel values */
        p2_out = __lsx_vbitsel_v(p2, p2_filt8_l, flat);
        p1_out = __lsx_vbitsel_v(p1_out, p1_filt8_l, flat);
        p0_out = __lsx_vbitsel_v(p0_out, p0_filt8_l, flat);
        q0_out = __lsx_vbitsel_v(q0_out, q0_filt8_l, flat);
        q1_out = __lsx_vbitsel_v(q1_out, q1_filt8_l, flat);
        q2_out = __lsx_vbitsel_v(q2, q2_filt8_l, flat);

        __lsx_vst(p2_out, filter48, 0);
        __lsx_vst(p1_out, filter48, 16);
        __lsx_vst(p0_out, filter48, 32);
        __lsx_vst(q0_out, filter48, 48);
        __lsx_vst(q1_out, filter48, 64);
        __lsx_vst(q2_out, filter48, 80);
        __lsx_vst(flat, filter48, 96);

        return 0;
    }
}

static void vp9_hz_lpf_t16_16w(uint8_t *dst, ptrdiff_t stride,
                               uint8_t *filter48)
{
    ptrdiff_t stride2 = stride << 1;
    ptrdiff_t stride3 = stride2 + stride;
    ptrdiff_t stride4 = stride2 << 1;
    uint8_t *dst_tmp = dst - stride4;
    uint8_t *dst_tmp1 = dst + stride4;
    __m128i p7, p6, p5, p4, p3, p2, p1, p0, q0, q1, q2, q3, q4, q5, q6, q7;
    __m128i flat, flat2, filter8;
    __m128i zero = __lsx_vldi(0);
    __m128i out_h, out_l;
    v8u16 p7_l_in, p6_l_in, p5_l_in, p4_l_in;
    v8u16 p3_l_in, p2_l_in, p1_l_in, p0_l_in;
    v8u16 q7_l_in, q6_l_in, q5_l_in, q4_l_in;
    v8u16 q3_l_in, q2_l_in, q1_l_in, q0_l_in;
    v8u16 p7_h_in, p6_h_in, p5_h_in, p4_h_in;
    v8u16 p3_h_in, p2_h_in, p1_h_in, p0_h_in;
    v8u16 q7_h_in, q6_h_in, q5_h_in, q4_h_in;
    v8u16 q3_h_in, q2_h_in, q1_h_in, q0_h_in;
    v8u16 tmp0_l, tmp1_l, tmp0_h, tmp1_h;

    flat = __lsx_vld(filter48, 96);

    DUP4_ARG2(__lsx_vldx, dst_tmp, -stride4, dst_tmp, -stride3, dst_tmp,
              -stride2, dst_tmp, -stride, p7, p6, p5, p4);
    p3 = __lsx_vld(dst_tmp, 0);
    DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, p2, p1);
    p0 = __lsx_vldx(dst_tmp, stride3);

    q0 = __lsx_vld(dst, 0);
    DUP2_ARG2(__lsx_vldx, dst, stride, dst, stride2, q1, q2);
    q3 = __lsx_vldx(dst, stride3);

    q4 = __lsx_vld(dst_tmp1, 0);
    DUP2_ARG2(__lsx_vldx, dst_tmp1, stride, dst_tmp1, stride2, q5, q6);
    q7 = __lsx_vldx(dst_tmp1, stride3);
    VP9_FLAT5(p7, p6, p5, p4, p0, q0, q4, q5, q6, q7, flat, flat2);

    /* if flat2 is zero for all pixels, then no need to calculate other filter */
    if (__lsx_bz_v(flat2)) {
        DUP4_ARG2(__lsx_vld, filter48, 0, filter48, 16, filter48, 32, filter48,
                  48, p2, p1, p0, q0);
        DUP2_ARG2(__lsx_vld, filter48, 64, filter48, 80, q1, q2);

        __lsx_vstx(p2, dst, -stride3);
        __lsx_vstx(p1, dst, -stride2);
        __lsx_vstx(p0, dst, -stride);
        __lsx_vst(q0, dst, 0);
        __lsx_vstx(q1, dst, stride);
        __lsx_vstx(q2, dst, stride2);
    } else {
        dst = dst_tmp - stride3;

        p7_l_in = (v8u16)__lsx_vilvl_b(zero, p7);
        p6_l_in = (v8u16)__lsx_vilvl_b(zero, p6);
        p5_l_in = (v8u16)__lsx_vilvl_b(zero, p5);
        p4_l_in = (v8u16)__lsx_vilvl_b(zero, p4);
        p3_l_in = (v8u16)__lsx_vilvl_b(zero, p3);
        p2_l_in = (v8u16)__lsx_vilvl_b(zero, p2);
        p1_l_in = (v8u16)__lsx_vilvl_b(zero, p1);
        p0_l_in = (v8u16)__lsx_vilvl_b(zero, p0);

        q0_l_in = (v8u16)__lsx_vilvl_b(zero, q0);

        tmp0_l = p7_l_in << 3;
        tmp0_l -= p7_l_in;
        tmp0_l += p6_l_in;
        tmp0_l += q0_l_in;
        tmp1_l = p6_l_in + p5_l_in;
        tmp1_l += p4_l_in;
        tmp1_l += p3_l_in;
        tmp1_l += p2_l_in;
        tmp1_l += p1_l_in;
        tmp1_l += p0_l_in;
        tmp1_l += tmp0_l;

        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);

        p7_h_in = (v8u16)__lsx_vilvh_b(zero, p7);
        p6_h_in = (v8u16)__lsx_vilvh_b(zero, p6);
        p5_h_in = (v8u16)__lsx_vilvh_b(zero, p5);
        p4_h_in = (v8u16)__lsx_vilvh_b(zero, p4);

        p3_h_in = (v8u16)__lsx_vilvh_b(zero, p3);
        p2_h_in = (v8u16)__lsx_vilvh_b(zero, p2);
        p1_h_in = (v8u16)__lsx_vilvh_b(zero, p1);
        p0_h_in = (v8u16)__lsx_vilvh_b(zero, p0);
        q0_h_in = (v8u16)__lsx_vilvh_b(zero, q0);

        tmp0_h = p7_h_in << 3;
        tmp0_h -= p7_h_in;
        tmp0_h += p6_h_in;
        tmp0_h += q0_h_in;
        tmp1_h = p6_h_in + p5_h_in;
        tmp1_h += p4_h_in;
        tmp1_h += p3_h_in;
        tmp1_h += p2_h_in;
        tmp1_h += p1_h_in;
        tmp1_h += p0_h_in;
        tmp1_h += tmp0_h;

        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);

        out_l = __lsx_vpickev_b(out_h, out_l);
        p6 = __lsx_vbitsel_v(p6, out_l, flat2);
        __lsx_vst(p6, dst, 0);
        dst += stride;

        /* p5 */
        q1_l_in = (v8u16)__lsx_vilvl_b(zero, q1);
        tmp0_l = p5_l_in - p6_l_in;
        tmp0_l += q1_l_in;
        tmp0_l -= p7_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);

        q1_h_in = (v8u16)__lsx_vilvh_b(zero, q1);
        tmp0_h = p5_h_in - p6_h_in;
        tmp0_h += q1_h_in;
        tmp0_h -= p7_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);

        out_l = __lsx_vpickev_b(out_h, out_l);
        p5 = __lsx_vbitsel_v(p5, out_l, flat2);
        __lsx_vst(p5, dst, 0);
        dst += stride;

        /* p4 */
        q2_l_in = (v8u16)__lsx_vilvl_b(zero, q2);
        tmp0_l = p4_l_in - p5_l_in;
        tmp0_l += q2_l_in;
        tmp0_l -= p7_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);

        q2_h_in = (v8u16)__lsx_vilvh_b(zero, q2);
        tmp0_h = p4_h_in - p5_h_in;
        tmp0_h += q2_h_in;
        tmp0_h -= p7_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);

        out_l = __lsx_vpickev_b(out_h, out_l);
        p4 = __lsx_vbitsel_v(p4, out_l, flat2);
        __lsx_vst(p4, dst, 0);
        dst += stride;

        /* p3 */
        q3_l_in = (v8u16)__lsx_vilvl_b(zero, q3);
        tmp0_l = p3_l_in - p4_l_in;
        tmp0_l += q3_l_in;
        tmp0_l -= p7_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);

        q3_h_in = (v8u16)__lsx_vilvh_b(zero, q3);
        tmp0_h = p3_h_in - p4_h_in;
        tmp0_h += q3_h_in;
        tmp0_h -= p7_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);

        out_l = __lsx_vpickev_b(out_h, out_l);
        p3 = __lsx_vbitsel_v(p3, out_l, flat2);
        __lsx_vst(p3, dst, 0);
        dst += stride;

        /* p2 */
        q4_l_in = (v8u16)__lsx_vilvl_b(zero, q4);
        filter8 = __lsx_vld(filter48, 0);
        tmp0_l = p2_l_in - p3_l_in;
        tmp0_l += q4_l_in;
        tmp0_l -= p7_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);

        q4_h_in = (v8u16)__lsx_vilvh_b(zero, q4);
        tmp0_h = p2_h_in - p3_h_in;
        tmp0_h += q4_h_in;
        tmp0_h -= p7_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);

        out_l = __lsx_vpickev_b(out_h, out_l);
        filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
        __lsx_vst(filter8, dst, 0);
        dst += stride;

        /* p1 */
        q5_l_in = (v8u16)__lsx_vilvl_b(zero, q5);
        filter8 = __lsx_vld(filter48, 16);
        tmp0_l = p1_l_in - p2_l_in;
        tmp0_l += q5_l_in;
        tmp0_l -= p7_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);

        q5_h_in = (v8u16)__lsx_vilvh_b(zero, q5);
        tmp0_h = p1_h_in - p2_h_in;
        tmp0_h += q5_h_in;
        tmp0_h -= p7_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);

        out_l = __lsx_vpickev_b(out_h, out_l);
        filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
        __lsx_vst(filter8, dst, 0);
        dst += stride;

        /* p0 */
        q6_l_in = (v8u16)__lsx_vilvl_b(zero, q6);
        filter8 = __lsx_vld(filter48, 32);
        tmp0_l = p0_l_in - p1_l_in;
        tmp0_l += q6_l_in;
        tmp0_l -= p7_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);

        q6_h_in = (v8u16)__lsx_vilvh_b(zero, q6);
        tmp0_h = p0_h_in - p1_h_in;
        tmp0_h += q6_h_in;
        tmp0_h -= p7_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);

        out_l = __lsx_vpickev_b(out_h, out_l);
        filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
        __lsx_vst(filter8, dst, 0);
        dst += stride;

        /* q0 */
        q7_l_in = (v8u16)__lsx_vilvl_b(zero, q7);
        filter8 = __lsx_vld(filter48, 48);
        tmp0_l = q7_l_in - p0_l_in;
        tmp0_l += q0_l_in;
        tmp0_l -= p7_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);

        q7_h_in = (v8u16)__lsx_vilvh_b(zero, q7);
        tmp0_h = q7_h_in - p0_h_in;
        tmp0_h += q0_h_in;
        tmp0_h -= p7_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);

        out_l = __lsx_vpickev_b(out_h, out_l);
        filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
        __lsx_vst(filter8, dst, 0);
        dst += stride;

        /* q1 */
        filter8 = __lsx_vld(filter48, 64);
        tmp0_l = q7_l_in - q0_l_in;
        tmp0_l += q1_l_in;
        tmp0_l -= p6_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);

        tmp0_h = q7_h_in - q0_h_in;
        tmp0_h += q1_h_in;
        tmp0_h -= p6_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);

        out_l = __lsx_vpickev_b(out_h, out_l);
        filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
        __lsx_vst(filter8, dst, 0);
        dst += stride;

        /* q2 */
        filter8 = __lsx_vld(filter48, 80);
        tmp0_l = q7_l_in - q1_l_in;
        tmp0_l += q2_l_in;
        tmp0_l -= p5_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);

        tmp0_h = q7_h_in - q1_h_in;
        tmp0_h += q2_h_in;
        tmp0_h -= p5_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);

        out_l = __lsx_vpickev_b(out_h, out_l);
        filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
        __lsx_vst(filter8, dst, 0);
        dst += stride;

        /* q3 */
        tmp0_l = q7_l_in - q2_l_in;
        tmp0_l += q3_l_in;
        tmp0_l -= p4_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);

        tmp0_h = q7_h_in - q2_h_in;
        tmp0_h += q3_h_in;
        tmp0_h -= p4_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);

        out_l = __lsx_vpickev_b(out_h, out_l);
        q3 = __lsx_vbitsel_v(q3, out_l, flat2);
        __lsx_vst(q3, dst, 0);
        dst += stride;

        /* q4 */
        tmp0_l = q7_l_in - q3_l_in;
        tmp0_l += q4_l_in;
        tmp0_l -= p3_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);

        tmp0_h = q7_h_in - q3_h_in;
        tmp0_h += q4_h_in;
        tmp0_h -= p3_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);

        out_l = __lsx_vpickev_b(out_h, out_l);
        q4 = __lsx_vbitsel_v(q4, out_l, flat2);
        __lsx_vst(q4, dst, 0);
        dst += stride;

        /* q5 */
        tmp0_l = q7_l_in - q4_l_in;
        tmp0_l += q5_l_in;
        tmp0_l -= p2_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);

        tmp0_h = q7_h_in - q4_h_in;
        tmp0_h += q5_h_in;
        tmp0_h -= p2_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);

        out_l = __lsx_vpickev_b(out_h, out_l);
        q5 = __lsx_vbitsel_v(q5, out_l, flat2);
        __lsx_vst(q5, dst, 0);
        dst += stride;

        /* q6 */
        tmp0_l = q7_l_in - q5_l_in;
        tmp0_l += q6_l_in;
        tmp0_l -= p1_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);

        tmp0_h = q7_h_in - q5_h_in;
        tmp0_h += q6_h_in;
        tmp0_h -= p1_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);

        out_l = __lsx_vpickev_b(out_h, out_l);
        q6 = __lsx_vbitsel_v(q6, out_l, flat2);
        __lsx_vst(q6, dst, 0);
    }
}

void ff_loop_filter_v_16_16_lsx(uint8_t *dst, ptrdiff_t stride,
                                int32_t b_limit_ptr,
                                int32_t limit_ptr,
                                int32_t thresh_ptr)
{
    uint8_t filter48[16 * 8] __attribute__ ((aligned(16)));
    uint8_t early_exit = 0;

    early_exit = vp9_hz_lpf_t4_and_t8_16w(dst, stride, &filter48[0],
                                          b_limit_ptr, limit_ptr, thresh_ptr);

    if (0 == early_exit) {
        vp9_hz_lpf_t16_16w(dst, stride, filter48);
    }
}

void ff_loop_filter_v_16_8_lsx(uint8_t *dst, ptrdiff_t stride,
                               int32_t b_limit_ptr,
                               int32_t limit_ptr,
                               int32_t thresh_ptr)
{
    ptrdiff_t stride2 = stride << 1;
    ptrdiff_t stride3 = stride2 + stride;
    ptrdiff_t stride4 = stride2 << 1;
    uint8_t *dst_tmp = dst - stride4;
    uint8_t *dst_tmp1 = dst + stride4;
    __m128i zero = __lsx_vldi(0);
    __m128i flat2, mask, hev, flat, thresh, b_limit, limit;
    __m128i p3, p2, p1, p0, q3, q2, q1, q0, p7, p6, p5, p4, q4, q5, q6, q7;
    __m128i p2_out, p1_out, p0_out, q0_out, q1_out, q2_out;
    __m128i p0_filter16, p1_filter16;
    __m128i p2_filter8, p1_filter8, p0_filter8;
    __m128i q0_filter8, q1_filter8, q2_filter8;
    __m128i p7_l, p6_l, p5_l, p4_l, q7_l, q6_l, q5_l, q4_l;
    __m128i p3_l, p2_l, p1_l, p0_l, q3_l, q2_l, q1_l, q0_l;
    __m128i tmp0, tmp1, tmp2;

    /* load vector elements */
    DUP4_ARG2(__lsx_vldx, dst, -stride4, dst, -stride3, dst, -stride2,
              dst, -stride, p3, p2, p1, p0);
    q0 = __lsx_vld(dst, 0);
    DUP2_ARG2(__lsx_vldx, dst, stride, dst, stride2, q1, q2);
    q3 = __lsx_vldx(dst, stride3);

    thresh  = __lsx_vreplgr2vr_b(thresh_ptr);
    b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
    limit   = __lsx_vreplgr2vr_b(limit_ptr);

    LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
                 hev, mask, flat);
    VP9_FLAT4(p3, p2, p0, q0, q2, q3, flat);
    VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
                       q1_out);

    flat = __lsx_vilvl_d(zero, flat);

    /* if flat is zero for all pixels, then no need to calculate other filter */
    if (__lsx_bz_v(flat)) {
        __lsx_vstelm_d(p1_out, dst - stride2, 0, 0);
        __lsx_vstelm_d(p0_out, dst -   stride, 0, 0);
        __lsx_vstelm_d(q0_out, dst           , 0, 0);
        __lsx_vstelm_d(q1_out, dst +   stride, 0, 0);
    } else {
        /* convert 8 bit input data into 16 bit */
        DUP4_ARG2(__lsx_vilvl_b, zero, p3, zero, p2, zero, p1, zero, p0,
                  p3_l, p2_l, p1_l, p0_l);
        DUP4_ARG2(__lsx_vilvl_b, zero, q0, zero, q1, zero, q2, zero, q3,
                  q0_l, q1_l, q2_l, q3_l);
        VP9_FILTER8(p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l,
                    p2_filter8, p1_filter8, p0_filter8, q0_filter8,
                    q1_filter8, q2_filter8);

        /* convert 16 bit output data into 8 bit */
        DUP4_ARG2(__lsx_vpickev_b, zero, p2_filter8, zero, p1_filter8,
                  zero, p0_filter8, zero, q0_filter8, p2_filter8,
                  p1_filter8, p0_filter8, q0_filter8);
        DUP2_ARG2(__lsx_vpickev_b, zero, q1_filter8, zero, q2_filter8,
                  q1_filter8, q2_filter8);

        /* store pixel values */
        p2_out = __lsx_vbitsel_v(p2, p2_filter8, flat);
        p1_out = __lsx_vbitsel_v(p1_out, p1_filter8, flat);
        p0_out = __lsx_vbitsel_v(p0_out, p0_filter8, flat);
        q0_out = __lsx_vbitsel_v(q0_out, q0_filter8, flat);
        q1_out = __lsx_vbitsel_v(q1_out, q1_filter8, flat);
        q2_out = __lsx_vbitsel_v(q2, q2_filter8, flat);

        /* load 16 vector elements */
        DUP4_ARG2(__lsx_vld, dst_tmp - stride4, 0, dst_tmp - stride3, 0,
                  dst_tmp - stride2, 0, dst_tmp - stride, 0, p7, p6, p5, p4);
        DUP4_ARG2(__lsx_vld, dst_tmp1, 0, dst_tmp1 + stride, 0,
                dst_tmp1 + stride2, 0, dst_tmp1 + stride3, 0, q4, q5, q6, q7);

        VP9_FLAT5(p7, p6, p5, p4, p0, q0, q4, q5, q6, q7, flat, flat2);

        /* if flat2 is zero for all pixels, then no need to calculate other filter */
        if (__lsx_bz_v(flat2)) {
            dst -= stride3;
            __lsx_vstelm_d(p2_out, dst, 0, 0);
            dst += stride;
            __lsx_vstelm_d(p1_out, dst, 0, 0);
            dst += stride;
            __lsx_vstelm_d(p0_out, dst, 0, 0);
            dst += stride;
            __lsx_vstelm_d(q0_out, dst, 0, 0);
            dst += stride;
            __lsx_vstelm_d(q1_out, dst, 0, 0);
            dst += stride;
            __lsx_vstelm_d(q2_out, dst, 0, 0);
        } else {
            /* LSB(right) 8 pixel operation */
            DUP4_ARG2(__lsx_vilvl_b, zero, p7, zero, p6, zero, p5, zero, p4,
                      p7_l, p6_l, p5_l, p4_l);
            DUP4_ARG2(__lsx_vilvl_b, zero, q4, zero, q5, zero, q6, zero, q7,
                      q4_l, q5_l, q6_l, q7_l);

            tmp0 = __lsx_vslli_h(p7_l, 3);
            tmp0 = __lsx_vsub_h(tmp0, p7_l);
            tmp0 = __lsx_vadd_h(tmp0, p6_l);
            tmp0 = __lsx_vadd_h(tmp0, q0_l);

            dst = dst_tmp - stride3;

            /* calculation of p6 and p5 */
            tmp1 = __lsx_vadd_h(p6_l, p5_l);
            tmp1 = __lsx_vadd_h(tmp1, p4_l);
            tmp1 = __lsx_vadd_h(tmp1, p3_l);
            tmp1 = __lsx_vadd_h(tmp1, p2_l);
            tmp1 = __lsx_vadd_h(tmp1, p1_l);
            tmp1 = __lsx_vadd_h(tmp1, p0_l);
            tmp1 = __lsx_vadd_h(tmp1, tmp0);

            p0_filter16 = __lsx_vsrari_h(tmp1, 4);
            tmp0 = __lsx_vsub_h(p5_l, p6_l);
            tmp0 = __lsx_vadd_h(tmp0, q1_l);
            tmp0 = __lsx_vsub_h(tmp0, p7_l);
            tmp1 = __lsx_vadd_h(tmp1, tmp0);

            p1_filter16 = __lsx_vsrari_h(tmp1, 4);
            DUP2_ARG2(__lsx_vpickev_b, zero, p0_filter16, zero,
                      p1_filter16, p0_filter16, p1_filter16);
            p0_filter16 = __lsx_vbitsel_v(p6, p0_filter16, flat2);
            p1_filter16 = __lsx_vbitsel_v(p5, p1_filter16, flat2);
            __lsx_vstelm_d(p0_filter16, dst, 0, 0);
            dst += stride;
            __lsx_vstelm_d(p1_filter16, dst, 0, 0);
            dst += stride;

            /* calculation of p4 and p3 */
            tmp0 = __lsx_vsub_h(p4_l, p5_l);
            tmp0 = __lsx_vadd_h(tmp0, q2_l);
            tmp0 = __lsx_vsub_h(tmp0, p7_l);
            tmp2 = __lsx_vsub_h(p3_l, p4_l);
            tmp2 = __lsx_vadd_h(tmp2, q3_l);
            tmp2 = __lsx_vsub_h(tmp2, p7_l);
            tmp1 = __lsx_vadd_h(tmp1, tmp0);
            p0_filter16 = __lsx_vsrari_h(tmp1, 4);
            tmp1 = __lsx_vadd_h(tmp1, tmp2);
            p1_filter16 = __lsx_vsrari_h(tmp1, 4);
            DUP2_ARG2(__lsx_vpickev_b, zero, p0_filter16, zero,
                      p1_filter16, p0_filter16, p1_filter16);
            p0_filter16 = __lsx_vbitsel_v(p4, p0_filter16, flat2);
            p1_filter16 = __lsx_vbitsel_v(p3, p1_filter16, flat2);
            __lsx_vstelm_d(p0_filter16, dst, 0, 0);
            dst += stride;
            __lsx_vstelm_d(p1_filter16, dst, 0, 0);
            dst += stride;

            /* calculation of p2 and p1 */
            tmp0 = __lsx_vsub_h(p2_l, p3_l);
            tmp0 = __lsx_vadd_h(tmp0, q4_l);
            tmp0 = __lsx_vsub_h(tmp0, p7_l);
            tmp2 = __lsx_vsub_h(p1_l, p2_l);
            tmp2 = __lsx_vadd_h(tmp2, q5_l);
            tmp2 = __lsx_vsub_h(tmp2, p7_l);
            tmp1 = __lsx_vadd_h(tmp1, tmp0);
            p0_filter16 = __lsx_vsrari_h(tmp1, 4);
            tmp1 = __lsx_vadd_h(tmp1, tmp2);
            p1_filter16 = __lsx_vsrari_h(tmp1, 4);
            DUP2_ARG2(__lsx_vpickev_b, zero, p0_filter16, zero,
                      p1_filter16, p0_filter16, p1_filter16);
            p0_filter16 = __lsx_vbitsel_v(p2_out, p0_filter16, flat2);
            p1_filter16 = __lsx_vbitsel_v(p1_out, p1_filter16, flat2);
            __lsx_vstelm_d(p0_filter16, dst, 0, 0);
            dst += stride;
            __lsx_vstelm_d(p1_filter16, dst, 0, 0);
            dst += stride;

            /* calculation of p0 and q0 */
            tmp0 = __lsx_vsub_h(p0_l, p1_l);
            tmp0 = __lsx_vadd_h(tmp0, q6_l);
            tmp0 = __lsx_vsub_h(tmp0, p7_l);
            tmp2 = __lsx_vsub_h(q7_l, p0_l);
            tmp2 = __lsx_vadd_h(tmp2, q0_l);
            tmp2 = __lsx_vsub_h(tmp2, p7_l);
            tmp1 = __lsx_vadd_h(tmp1, tmp0);
            p0_filter16 = __lsx_vsrari_h((__m128i)tmp1, 4);
            tmp1 = __lsx_vadd_h(tmp1, tmp2);
            p1_filter16 = __lsx_vsrari_h((__m128i)tmp1, 4);
            DUP2_ARG2(__lsx_vpickev_b, zero, p0_filter16, zero,
                      p1_filter16, p0_filter16, p1_filter16);
            p0_filter16 = __lsx_vbitsel_v(p0_out, p0_filter16, flat2);
            p1_filter16 = __lsx_vbitsel_v(q0_out, p1_filter16, flat2);
            __lsx_vstelm_d(p0_filter16, dst, 0, 0);
            dst += stride;
            __lsx_vstelm_d(p1_filter16, dst, 0, 0);
            dst += stride;

            /* calculation of q1 and q2 */
            tmp0 = __lsx_vsub_h(q7_l, q0_l);
            tmp0 = __lsx_vadd_h(tmp0, q1_l);
            tmp0 = __lsx_vsub_h(tmp0, p6_l);
            tmp2 = __lsx_vsub_h(q7_l, q1_l);
            tmp2 = __lsx_vadd_h(tmp2, q2_l);
            tmp2 = __lsx_vsub_h(tmp2, p5_l);
            tmp1 = __lsx_vadd_h(tmp1, tmp0);
            p0_filter16 = __lsx_vsrari_h(tmp1, 4);
            tmp1 = __lsx_vadd_h(tmp1, tmp2);
            p1_filter16 = __lsx_vsrari_h(tmp1, 4);
            DUP2_ARG2(__lsx_vpickev_b, zero, p0_filter16, zero,
                      p1_filter16, p0_filter16, p1_filter16);
            p0_filter16 = __lsx_vbitsel_v(q1_out, p0_filter16, flat2);
            p1_filter16 = __lsx_vbitsel_v(q2_out, p1_filter16, flat2);
            __lsx_vstelm_d(p0_filter16, dst, 0, 0);
            dst += stride;
            __lsx_vstelm_d(p1_filter16, dst, 0, 0);
            dst += stride;

            /* calculation of q3 and q4 */
            tmp0 = __lsx_vsub_h(q7_l, q2_l);
            tmp0 = __lsx_vadd_h(tmp0, q3_l);
            tmp0 = __lsx_vsub_h(tmp0, p4_l);
            tmp2 = __lsx_vsub_h(q7_l, q3_l);
            tmp2 = __lsx_vadd_h(tmp2, q4_l);
            tmp2 = __lsx_vsub_h(tmp2, p3_l);
            tmp1 = __lsx_vadd_h(tmp1, tmp0);
            p0_filter16 = __lsx_vsrari_h(tmp1, 4);
            tmp1 = __lsx_vadd_h(tmp1, tmp2);
            p1_filter16 = __lsx_vsrari_h(tmp1, 4);
            DUP2_ARG2(__lsx_vpickev_b, zero, p0_filter16, zero,
                      p1_filter16, p0_filter16, p1_filter16);
            p0_filter16 = __lsx_vbitsel_v(q3, p0_filter16, flat2);
            p1_filter16 = __lsx_vbitsel_v(q4, p1_filter16, flat2);
            __lsx_vstelm_d(p0_filter16, dst, 0, 0);
            dst += stride;
            __lsx_vstelm_d(p1_filter16, dst, 0, 0);
            dst += stride;

            /* calculation of q5 and q6 */
            tmp0 = __lsx_vsub_h(q7_l, q4_l);
            tmp0 = __lsx_vadd_h(tmp0, q5_l);
            tmp0 = __lsx_vsub_h(tmp0, p2_l);
            tmp2 = __lsx_vsub_h(q7_l, q5_l);
            tmp2 = __lsx_vadd_h(tmp2, q6_l);
            tmp2 = __lsx_vsub_h(tmp2, p1_l);
            tmp1 = __lsx_vadd_h(tmp1, tmp0);
            p0_filter16 = __lsx_vsrari_h(tmp1, 4);
            tmp1 = __lsx_vadd_h(tmp1, tmp2);
            p1_filter16 = __lsx_vsrari_h(tmp1, 4);
            DUP2_ARG2(__lsx_vpickev_b, zero, p0_filter16, zero,
                      p1_filter16, p0_filter16, p1_filter16);
            p0_filter16 = __lsx_vbitsel_v(q5, p0_filter16, flat2);
            p1_filter16 = __lsx_vbitsel_v(q6, p1_filter16, flat2);
            __lsx_vstelm_d(p0_filter16, dst, 0, 0);
            dst += stride;
            __lsx_vstelm_d(p1_filter16, dst, 0, 0);
        }
    }
}

void ff_loop_filter_h_4_8_lsx(uint8_t *dst, ptrdiff_t stride,
                              int32_t b_limit_ptr,
                              int32_t limit_ptr,
                              int32_t thresh_ptr)
{
    ptrdiff_t stride2 = stride << 1;
    ptrdiff_t stride3 = stride2 + stride;
    ptrdiff_t stride4 = stride2 << 1;
    uint8_t *dst_tmp1 = dst - 4;
    uint8_t *dst_tmp2 = dst_tmp1 + stride4;
    __m128i mask, hev, flat, limit, thresh, b_limit;
    __m128i p3, p2, p1, p0, q3, q2, q1, q0;
    __m128i vec0, vec1, vec2, vec3;

    p3 = __lsx_vld(dst_tmp1, 0);
    DUP2_ARG2(__lsx_vldx, dst_tmp1, stride, dst_tmp1, stride2, p2, p1);
    p0 = __lsx_vldx(dst_tmp1, stride3);
    q0 = __lsx_vld(dst_tmp2, 0);
    DUP2_ARG2(__lsx_vldx, dst_tmp2, stride, dst_tmp2, stride2, q1, q2);
    q3 = __lsx_vldx(dst_tmp2, stride3);

    thresh  = __lsx_vreplgr2vr_b(thresh_ptr);
    b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
    limit   = __lsx_vreplgr2vr_b(limit_ptr);

    LSX_TRANSPOSE8x8_B(p3, p2, p1, p0, q0, q1, q2, q3,
                       p3, p2, p1, p0, q0, q1, q2, q3);
    LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
                 hev, mask, flat);
    VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1, p0, q0, q1);
    DUP2_ARG2(__lsx_vilvl_b, p0, p1, q1, q0, vec0, vec1);
    vec2 = __lsx_vilvl_h(vec1, vec0);
    vec3 = __lsx_vilvh_h(vec1, vec0);

    dst -= 2;
    __lsx_vstelm_w(vec2, dst, 0, 0);
    __lsx_vstelm_w(vec2, dst + stride, 0, 1);
    __lsx_vstelm_w(vec2, dst + stride2, 0, 2);
    __lsx_vstelm_w(vec2, dst + stride3, 0, 3);
    dst += stride4;
    __lsx_vstelm_w(vec3, dst, 0, 0);
    __lsx_vstelm_w(vec3, dst + stride, 0, 1);
    __lsx_vstelm_w(vec3, dst + stride2, 0, 2);
    __lsx_vstelm_w(vec3, dst + stride3, 0, 3);
}

void ff_loop_filter_h_44_16_lsx(uint8_t *dst, ptrdiff_t stride,
                                int32_t b_limit_ptr,
                                int32_t limit_ptr,
                                int32_t thresh_ptr)
{
    ptrdiff_t stride2 = stride << 1;
    ptrdiff_t stride3 = stride2 + stride;
    ptrdiff_t stride4 = stride2 << 1;
    uint8_t *dst_tmp = dst - 4;
    __m128i mask, hev, flat;
    __m128i thresh0, b_limit0, limit0, thresh1, b_limit1, limit1;
    __m128i p3, p2, p1, p0, q3, q2, q1, q0;
    __m128i row0, row1, row2, row3, row4, row5, row6, row7;
    __m128i row8, row9, row10, row11, row12, row13, row14, row15;
    __m128i tmp0, tmp1, tmp2, tmp3, tmp4, tmp5;

    row0 = __lsx_vld(dst_tmp, 0);
    DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, row1, row2);
    row3 = __lsx_vldx(dst_tmp, stride3);
    dst_tmp += stride4;
    row4 = __lsx_vld(dst_tmp, 0);
    DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, row5, row6);
    row7 = __lsx_vldx(dst_tmp, stride3);
    dst_tmp += stride4;
    row8 = __lsx_vld(dst_tmp, 0);
    DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, row9, row10);
    row11 = __lsx_vldx(dst_tmp, stride3);
    dst_tmp += stride4;
    row12 = __lsx_vld(dst_tmp, 0);
    DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, row13, row14);
    row15 = __lsx_vldx(dst_tmp, stride3);

    LSX_TRANSPOSE16x8_B(row0, row1, row2, row3, row4, row5, row6, row7,
                        row8, row9, row10, row11, row12, row13, row14, row15,
                        p3, p2, p1, p0, q0, q1, q2, q3);

    thresh0 = __lsx_vreplgr2vr_b(thresh_ptr);
    thresh1 = __lsx_vreplgr2vr_b(thresh_ptr >> 8);
    thresh0 = __lsx_vilvl_d(thresh1, thresh0);

    b_limit0 = __lsx_vreplgr2vr_b(b_limit_ptr);
    b_limit1 = __lsx_vreplgr2vr_b(b_limit_ptr >> 8);
    b_limit0 = __lsx_vilvl_d(b_limit1, b_limit0);

    limit0 = __lsx_vreplgr2vr_b(limit_ptr);
    limit1 = __lsx_vreplgr2vr_b(limit_ptr >> 8);
    limit0 = __lsx_vilvl_d(limit1, limit0);

    LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit0, b_limit0, thresh0,
                 hev, mask, flat);
    VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1, p0, q0, q1);
    DUP2_ARG2(__lsx_vilvl_b, p0, p1, q1, q0, tmp0, tmp1);
    tmp2 = __lsx_vilvl_h(tmp1, tmp0);
    tmp3 = __lsx_vilvh_h(tmp1, tmp0);
    DUP2_ARG2(__lsx_vilvh_b, p0, p1, q1, q0, tmp0, tmp1);
    tmp4 = __lsx_vilvl_h(tmp1, tmp0);
    tmp5 = __lsx_vilvh_h(tmp1, tmp0);

    dst -= 2;
    __lsx_vstelm_w(tmp2, dst, 0, 0);
    __lsx_vstelm_w(tmp2, dst + stride, 0, 1);
    __lsx_vstelm_w(tmp2, dst + stride2, 0, 2);
    __lsx_vstelm_w(tmp2, dst + stride3, 0, 3);
    dst += stride4;
    __lsx_vstelm_w(tmp3, dst, 0, 0);
    __lsx_vstelm_w(tmp3, dst + stride, 0, 1);
    __lsx_vstelm_w(tmp3, dst + stride2, 0, 2);
    __lsx_vstelm_w(tmp3, dst + stride3, 0, 3);
    dst += stride4;
    __lsx_vstelm_w(tmp4, dst, 0, 0);
    __lsx_vstelm_w(tmp4, dst + stride, 0, 1);
    __lsx_vstelm_w(tmp4, dst + stride2, 0, 2);
    __lsx_vstelm_w(tmp4, dst + stride3, 0, 3);
    dst += stride4;
    __lsx_vstelm_w(tmp5, dst, 0, 0);
    __lsx_vstelm_w(tmp5, dst + stride, 0, 1);
    __lsx_vstelm_w(tmp5, dst + stride2, 0, 2);
    __lsx_vstelm_w(tmp5, dst + stride3, 0, 3);
}

void ff_loop_filter_h_8_8_lsx(uint8_t *dst, ptrdiff_t stride,
                              int32_t b_limit_ptr,
                              int32_t limit_ptr,
                              int32_t thresh_ptr)
{
    ptrdiff_t stride2 = stride << 1;
    ptrdiff_t stride3 = stride2 + stride;
    ptrdiff_t stride4 = stride2 << 1;
    uint8_t *dst_tmp = dst - 4;
    __m128i p3, p2, p1, p0, q3, q2, q1, q0;
    __m128i p1_out, p0_out, q0_out, q1_out;
    __m128i flat, mask, hev, thresh, b_limit, limit;
    __m128i p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l;
    __m128i p2_filt8_l, p1_filt8_l, p0_filt8_l;
    __m128i q0_filt8_l, q1_filt8_l, q2_filt8_l;
    __m128i vec0, vec1, vec2, vec3, vec4;
    __m128i zero = __lsx_vldi(0);

    /* load vector elements */
    p3 = __lsx_vld(dst_tmp, 0);
    DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, p2, p1);
    p0 = __lsx_vldx(dst_tmp, stride3);
    dst_tmp += stride4;
    q0 = __lsx_vld(dst_tmp, 0);
    DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, q1, q2);
    q3 = __lsx_vldx(dst_tmp, stride3);

    LSX_TRANSPOSE8x8_B(p3, p2, p1, p0, q0, q1, q2, q3,
                       p3, p2, p1, p0, q0, q1, q2, q3);

    thresh  = __lsx_vreplgr2vr_b(thresh_ptr);
    b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
    limit   = __lsx_vreplgr2vr_b(limit_ptr);

    /* mask and hev */
    LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
                 hev, mask, flat);
    /* flat4 */
    VP9_FLAT4(p3, p2, p0, q0, q2, q3, flat);
    /* filter4 */
    VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
                       q1_out);

    flat = __lsx_vilvl_d(zero, flat);

    /* if flat is zero for all pixels, then no need to calculate other filter */
    if (__lsx_bz_v(flat)) {
        /* Store 4 pixels p1-_q1 */
        DUP2_ARG2(__lsx_vilvl_b, p0_out, p1_out, q1_out, q0_out, vec0, vec1);
        vec2 = __lsx_vilvl_h(vec1, vec0);
        vec3 = __lsx_vilvh_h(vec1, vec0);

        dst -= 2;
        __lsx_vstelm_w(vec2, dst, 0, 0);
        __lsx_vstelm_w(vec2, dst + stride, 0, 1);
        __lsx_vstelm_w(vec2, dst + stride2, 0, 2);
        __lsx_vstelm_w(vec2, dst + stride3, 0, 3);
        dst += stride4;
        __lsx_vstelm_w(vec3, dst, 0, 0);
        __lsx_vstelm_w(vec3, dst + stride, 0, 1);
        __lsx_vstelm_w(vec3, dst + stride2, 0, 2);
        __lsx_vstelm_w(vec3, dst + stride3, 0, 3);
    } else {
        DUP4_ARG2(__lsx_vilvl_b, zero, p3, zero, p2, zero, p1, zero, p0,
                  p3_l, p2_l, p1_l, p0_l);
        DUP4_ARG2(__lsx_vilvl_b, zero, q0, zero, q1, zero, q2, zero, q3,
                  q0_l, q1_l, q2_l, q3_l);
        VP9_FILTER8(p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l, p2_filt8_l,
                    p1_filt8_l, p0_filt8_l, q0_filt8_l, q1_filt8_l, q2_filt8_l);
        /* convert 16 bit output data into 8 bit */
        DUP4_ARG2(__lsx_vpickev_b, p2_filt8_l, p2_filt8_l, p1_filt8_l,
                  p1_filt8_l, p0_filt8_l, p0_filt8_l, q0_filt8_l,
                  q0_filt8_l, p2_filt8_l, p1_filt8_l, p0_filt8_l,
                  q0_filt8_l);
        DUP2_ARG2(__lsx_vpickev_b, q1_filt8_l, q1_filt8_l, q2_filt8_l,
                  q2_filt8_l, q1_filt8_l, q2_filt8_l);

        /* store pixel values */
        p2 = __lsx_vbitsel_v(p2, p2_filt8_l, flat);
        p1 = __lsx_vbitsel_v(p1_out, p1_filt8_l, flat);
        p0 = __lsx_vbitsel_v(p0_out, p0_filt8_l, flat);
        q0 = __lsx_vbitsel_v(q0_out, q0_filt8_l, flat);
        q1 = __lsx_vbitsel_v(q1_out, q1_filt8_l, flat);
        q2 = __lsx_vbitsel_v(q2, q2_filt8_l, flat);

        /* Store 6 pixels p2-_q2 */
        DUP2_ARG2(__lsx_vilvl_b, p1, p2, q0, p0, vec0, vec1);
        vec2 = __lsx_vilvl_h(vec1, vec0);
        vec3 = __lsx_vilvh_h(vec1, vec0);
        vec4 = __lsx_vilvl_b(q2, q1);

        dst -= 3;
        __lsx_vstelm_w(vec2, dst, 0, 0);
        __lsx_vstelm_h(vec4, dst, 4, 0);
        dst += stride;
        __lsx_vstelm_w(vec2, dst, 0, 1);
        __lsx_vstelm_h(vec4, dst, 4, 1);
        dst += stride;
        __lsx_vstelm_w(vec2, dst, 0, 2);
        __lsx_vstelm_h(vec4, dst, 4, 2);
        dst += stride;
        __lsx_vstelm_w(vec2, dst, 0, 3);
        __lsx_vstelm_h(vec4, dst, 4, 3);
        dst += stride;
        __lsx_vstelm_w(vec3, dst, 0, 0);
        __lsx_vstelm_h(vec4, dst, 4, 4);
        dst += stride;
        __lsx_vstelm_w(vec3, dst, 0, 1);
        __lsx_vstelm_h(vec4, dst, 4, 5);
        dst += stride;
        __lsx_vstelm_w(vec3, dst, 0, 2);
        __lsx_vstelm_h(vec4, dst, 4, 6);
        dst += stride;
        __lsx_vstelm_w(vec3, dst, 0, 3);
        __lsx_vstelm_h(vec4, dst, 4, 7);
    }
}

void ff_loop_filter_h_88_16_lsx(uint8_t *dst, ptrdiff_t stride,
                                int32_t b_limit_ptr,
                                int32_t limit_ptr,
                                int32_t thresh_ptr)
{
    ptrdiff_t stride2 = stride << 1;
    ptrdiff_t stride3 = stride2 + stride;
    ptrdiff_t stride4 = stride2 << 1;
    uint8_t *dst_tmp = dst - 4;
    __m128i p3, p2, p1, p0, q3, q2, q1, q0;
    __m128i p1_out, p0_out, q0_out, q1_out;
    __m128i flat, mask, hev, thresh, b_limit, limit;
    __m128i row4, row5, row6, row7, row12, row13, row14, row15;
    __m128i p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l;
    __m128i p3_h, p2_h, p1_h, p0_h, q0_h, q1_h, q2_h, q3_h;
    __m128i p2_filt8_l, p1_filt8_l, p0_filt8_l;
    __m128i q0_filt8_l, q1_filt8_l, q2_filt8_l;
    __m128i p2_filt8_h, p1_filt8_h, p0_filt8_h;
    __m128i q0_filt8_h, q1_filt8_h, q2_filt8_h;
    __m128i vec0, vec1, vec2, vec3, vec4, vec5, vec6, vec7;
    __m128i zero = __lsx_vldi(0);

    p0 = __lsx_vld(dst_tmp, 0);
    DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, p1, p2);
    p3 = __lsx_vldx(dst_tmp, stride3);
    dst_tmp += stride4;
    row4 = __lsx_vld(dst_tmp, 0);
    DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, row5, row6);
    row7 = __lsx_vldx(dst_tmp, stride3);
    dst_tmp += stride4;
    q3 = __lsx_vld(dst_tmp, 0);
    DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, q2, q1);
    q0 = __lsx_vldx(dst_tmp, stride3);
    dst_tmp += stride4;
    row12 = __lsx_vld(dst_tmp, 0);
    DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, row13, row14);
    row15 = __lsx_vldx(dst_tmp, stride3);

    /* transpose 16x8 matrix into 8x16 */
    LSX_TRANSPOSE16x8_B(p0, p1, p2, p3, row4, row5, row6, row7,
                        q3, q2, q1, q0, row12, row13, row14, row15,
                        p3, p2, p1, p0, q0, q1, q2, q3);

    thresh = __lsx_vreplgr2vr_b(thresh_ptr);
    vec0   = __lsx_vreplgr2vr_b(thresh_ptr >> 8);
    thresh = __lsx_vilvl_d(vec0, thresh);

    b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
    vec0    = __lsx_vreplgr2vr_b(b_limit_ptr >> 8);
    b_limit = __lsx_vilvl_d(vec0, b_limit);

    limit = __lsx_vreplgr2vr_b(limit_ptr);
    vec0  = __lsx_vreplgr2vr_b(limit_ptr >> 8);
    limit = __lsx_vilvl_d(vec0, limit);

    /* mask and hev */
    LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
                 hev, mask, flat);
    /* flat4 */
    VP9_FLAT4(p3, p2, p0, q0, q2, q3, flat);
    /* filter4 */
    VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
                       q1_out);

    /* if flat is zero for all pixels, then no need to calculate other filter */
    if (__lsx_bz_v(flat)) {
        DUP2_ARG2(__lsx_vilvl_b, p0_out, p1_out, q1_out, q0_out, vec0, vec1);
        vec2 = __lsx_vilvl_h(vec1, vec0);
        vec3 = __lsx_vilvh_h(vec1, vec0);
        DUP2_ARG2(__lsx_vilvh_b, p0_out, p1_out, q1_out, q0_out, vec0, vec1);
        vec4 = __lsx_vilvl_h(vec1, vec0);
        vec5 = __lsx_vilvh_h(vec1, vec0);

        dst -= 2;
        __lsx_vstelm_w(vec2, dst, 0, 0);
        __lsx_vstelm_w(vec2, dst + stride, 0, 1);
        __lsx_vstelm_w(vec2, dst + stride2, 0, 2);
        __lsx_vstelm_w(vec2, dst + stride3, 0, 3);
        dst += stride4;
        __lsx_vstelm_w(vec3, dst, 0, 0);
        __lsx_vstelm_w(vec3, dst + stride, 0, 1);
        __lsx_vstelm_w(vec3, dst + stride2, 0, 2);
        __lsx_vstelm_w(vec3, dst + stride3, 0, 3);
        dst += stride4;
        __lsx_vstelm_w(vec4, dst, 0, 0);
        __lsx_vstelm_w(vec4, dst + stride, 0, 1);
        __lsx_vstelm_w(vec4, dst + stride2, 0, 2);
        __lsx_vstelm_w(vec4, dst + stride3, 0, 3);
        dst += stride4;
        __lsx_vstelm_w(vec5, dst, 0, 0);
        __lsx_vstelm_w(vec5, dst + stride, 0, 1);
        __lsx_vstelm_w(vec5, dst + stride2, 0, 2);
        __lsx_vstelm_w(vec5, dst + stride3, 0, 3);
    } else {
        DUP4_ARG2(__lsx_vilvl_b, zero, p3, zero, p2, zero, p1, zero, p0,
                  p3_l, p2_l, p1_l, p0_l);
        DUP4_ARG2(__lsx_vilvl_b, zero, q0, zero, q1, zero, q2, zero, q3,
                  q0_l, q1_l, q2_l, q3_l);
        VP9_FILTER8(p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l, p2_filt8_l,
                    p1_filt8_l, p0_filt8_l, q0_filt8_l, q1_filt8_l, q2_filt8_l);

        DUP4_ARG2(__lsx_vilvh_b, zero, p3, zero, p2, zero, p1, zero, p0,
                  p3_h, p2_h, p1_h, p0_h);
        DUP4_ARG2(__lsx_vilvh_b, zero, q0, zero, q1, zero, q2, zero, q3,
                  q0_h, q1_h, q2_h, q3_h);

        /* filter8 */
        VP9_FILTER8(p3_h, p2_h, p1_h, p0_h, q0_h, q1_h, q2_h, q3_h, p2_filt8_h,
                    p1_filt8_h, p0_filt8_h, q0_filt8_h, q1_filt8_h, q2_filt8_h);

        /* convert 16 bit output data into 8 bit */
        DUP4_ARG2(__lsx_vpickev_b, p2_filt8_h, p2_filt8_l, p1_filt8_h,
                  p1_filt8_l, p0_filt8_h, p0_filt8_l, q0_filt8_h, q0_filt8_l,
                  p2_filt8_l, p1_filt8_l, p0_filt8_l, q0_filt8_l);
        DUP2_ARG2(__lsx_vpickev_b, q1_filt8_h, q1_filt8_l, q2_filt8_h,
                  q2_filt8_l, q1_filt8_l, q2_filt8_l);

        /* store pixel values */
        p2 = __lsx_vbitsel_v(p2, p2_filt8_l, flat);
        p1 = __lsx_vbitsel_v(p1_out, p1_filt8_l, flat);
        p0 = __lsx_vbitsel_v(p0_out, p0_filt8_l, flat);
        q0 = __lsx_vbitsel_v(q0_out, q0_filt8_l, flat);
        q1 = __lsx_vbitsel_v(q1_out, q1_filt8_l, flat);
        q2 = __lsx_vbitsel_v(q2, q2_filt8_l, flat);

        DUP2_ARG2(__lsx_vilvl_b, p1, p2, q0, p0, vec0, vec1);
        vec3 = __lsx_vilvl_h(vec1, vec0);
        vec4 = __lsx_vilvh_h(vec1, vec0);
        DUP2_ARG2(__lsx_vilvh_b, p1, p2, q0, p0, vec0, vec1);
        vec6 = __lsx_vilvl_h(vec1, vec0);
        vec7 = __lsx_vilvh_h(vec1, vec0);
        vec2 = __lsx_vilvl_b(q2, q1);
        vec5 = __lsx_vilvh_b(q2, q1);

        dst -= 3;
        __lsx_vstelm_w(vec3, dst, 0, 0);
        __lsx_vstelm_h(vec2, dst, 4, 0);
        dst += stride;
        __lsx_vstelm_w(vec3, dst, 0, 1);
        __lsx_vstelm_h(vec2, dst, 4, 1);
        dst += stride;
        __lsx_vstelm_w(vec3, dst, 0, 2);
        __lsx_vstelm_h(vec2, dst, 4, 2);
        dst += stride;
        __lsx_vstelm_w(vec3, dst, 0, 3);
        __lsx_vstelm_h(vec2, dst, 4, 3);
        dst += stride;
        __lsx_vstelm_w(vec4, dst, 0, 0);
        __lsx_vstelm_h(vec2, dst, 4, 4);
        dst += stride;
        __lsx_vstelm_w(vec4, dst, 0, 1);
        __lsx_vstelm_h(vec2, dst, 4, 5);
        dst += stride;
        __lsx_vstelm_w(vec4, dst, 0, 2);
        __lsx_vstelm_h(vec2, dst, 4, 6);
        dst += stride;
        __lsx_vstelm_w(vec4, dst, 0, 3);
        __lsx_vstelm_h(vec2, dst, 4, 7);
        dst += stride;
        __lsx_vstelm_w(vec6, dst, 0, 0);
        __lsx_vstelm_h(vec5, dst, 4, 0);
        dst += stride;
        __lsx_vstelm_w(vec6, dst, 0, 1);
        __lsx_vstelm_h(vec5, dst, 4, 1);
        dst += stride;
        __lsx_vstelm_w(vec6, dst, 0, 2);
        __lsx_vstelm_h(vec5, dst, 4, 2);
        dst += stride;
        __lsx_vstelm_w(vec6, dst, 0, 3);
        __lsx_vstelm_h(vec5, dst, 4, 3);
        dst += stride;
        __lsx_vstelm_w(vec7, dst, 0, 0);
        __lsx_vstelm_h(vec5, dst, 4, 4);
        dst += stride;
        __lsx_vstelm_w(vec7, dst, 0, 1);
        __lsx_vstelm_h(vec5, dst, 4, 5);
        dst += stride;
        __lsx_vstelm_w(vec7, dst, 0, 2);
        __lsx_vstelm_h(vec5, dst, 4, 6);
        dst += stride;
        __lsx_vstelm_w(vec7, dst, 0, 3);
        __lsx_vstelm_h(vec5, dst, 4, 7);
    }
}

void ff_loop_filter_h_84_16_lsx(uint8_t *dst, ptrdiff_t stride,
                                int32_t b_limit_ptr,
                                int32_t limit_ptr,
                                int32_t thresh_ptr)
{
    ptrdiff_t stride2 = stride << 1;
    ptrdiff_t stride3 = stride2 + stride;
    ptrdiff_t stride4 = stride2 << 1;
    uint8_t *dst_tmp = dst - 4;
    __m128i p3, p2, p1, p0, q3, q2, q1, q0;
    __m128i p1_out, p0_out, q0_out, q1_out;
    __m128i flat, mask, hev, thresh, b_limit, limit;
    __m128i row4, row5, row6, row7, row12, row13, row14, row15;
    __m128i p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l;
    __m128i p2_filt8_l, p1_filt8_l, p0_filt8_l;
    __m128i q0_filt8_l, q1_filt8_l, q2_filt8_l;
    __m128i vec0, vec1, vec2, vec3, vec4, vec5, vec6, vec7;
    __m128i zero = __lsx_vldi(0);

    p0 = __lsx_vld(dst_tmp, 0);
    DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, p1, p2);
    p3 = __lsx_vldx(dst_tmp, stride3);
    dst_tmp += stride4;
    row4 = __lsx_vld(dst_tmp, 0);
    DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, row5, row6);
    row7 = __lsx_vldx(dst_tmp, stride3);
    dst_tmp += stride4;
    q3 = __lsx_vld(dst_tmp, 0);
    DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, q2, q1);
    q0 = __lsx_vldx(dst_tmp, stride3);
    dst_tmp += stride4;
    row12 = __lsx_vld(dst_tmp, 0);
    DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, row13, row14);
    row15 = __lsx_vldx(dst_tmp, stride3);

    /* transpose 16x8 matrix into 8x16 */
    LSX_TRANSPOSE16x8_B(p0, p1, p2, p3, row4, row5, row6, row7,
                        q3, q2, q1, q0, row12, row13, row14, row15,
                        p3, p2, p1, p0, q0, q1, q2, q3);

    thresh = __lsx_vreplgr2vr_b(thresh_ptr);
    vec0   = __lsx_vreplgr2vr_b(thresh_ptr >> 8);
    thresh = __lsx_vilvl_d(vec0, thresh);

    b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
    vec0    = __lsx_vreplgr2vr_b(b_limit_ptr >> 8);
    b_limit = __lsx_vilvl_d(vec0, b_limit);

    limit = __lsx_vreplgr2vr_b(limit_ptr);
    vec0  = __lsx_vreplgr2vr_b(limit_ptr >> 8);
    limit = __lsx_vilvl_d(vec0, limit);

    /* mask and hev */
    LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
                 hev, mask, flat);
    /* flat4 */
    VP9_FLAT4(p3, p2, p0, q0, q2, q3, flat);
    /* filter4 */
    VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
                       q1_out);

    flat = __lsx_vilvl_d(zero, flat);

    /* if flat is zero for all pixels, then no need to calculate other filter */
    if (__lsx_bz_v(flat)) {
        DUP2_ARG2(__lsx_vilvl_b, p0_out, p1_out, q1_out, q0_out, vec0, vec1);
        vec2 = __lsx_vilvl_h(vec1, vec0);
        vec3 = __lsx_vilvh_h(vec1, vec0);
        DUP2_ARG2(__lsx_vilvh_b, p0_out, p1_out, q1_out, q0_out, vec0, vec1);
        vec4 = __lsx_vilvl_h(vec1, vec0);
        vec5 = __lsx_vilvh_h(vec1, vec0);

        dst -= 2;
        __lsx_vstelm_w(vec2, dst, 0, 0);
        __lsx_vstelm_w(vec2, dst + stride, 0, 1);
        __lsx_vstelm_w(vec2, dst + stride2, 0, 2);
        __lsx_vstelm_w(vec2, dst + stride3, 0, 3);
        dst += stride4;
        __lsx_vstelm_w(vec3, dst, 0, 0);
        __lsx_vstelm_w(vec3, dst + stride, 0, 1);
        __lsx_vstelm_w(vec3, dst + stride2, 0, 2);
        __lsx_vstelm_w(vec3, dst + stride3, 0, 3);
        dst += stride4;
        __lsx_vstelm_w(vec4, dst, 0, 0);
        __lsx_vstelm_w(vec4, dst + stride, 0, 1);
        __lsx_vstelm_w(vec4, dst + stride2, 0, 2);
        __lsx_vstelm_w(vec4, dst + stride3, 0, 3);
        dst += stride4;
        __lsx_vstelm_w(vec5, dst, 0, 0);
        __lsx_vstelm_w(vec5, dst + stride, 0, 1);
        __lsx_vstelm_w(vec5, dst + stride2, 0, 2);
        __lsx_vstelm_w(vec5, dst + stride3, 0, 3);
    } else {
        DUP4_ARG2(__lsx_vilvl_b, zero, p3, zero, p2, zero, p1, zero, p0,
                  p3_l, p2_l, p1_l, p0_l);
        DUP4_ARG2(__lsx_vilvl_b, zero, q0, zero, q1, zero, q2, zero, q3,
                  q0_l, q1_l, q2_l, q3_l);
        VP9_FILTER8(p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l, p2_filt8_l,
                    p1_filt8_l, p0_filt8_l, q0_filt8_l, q1_filt8_l, q2_filt8_l);

        /* convert 16 bit output data into 8 bit */
        DUP4_ARG2(__lsx_vpickev_b, p2_filt8_l, p2_filt8_l, p1_filt8_l, p1_filt8_l,
                  p0_filt8_l, p0_filt8_l, q0_filt8_l, q0_filt8_l, p2_filt8_l,
                  p1_filt8_l, p0_filt8_l, q0_filt8_l);
        DUP2_ARG2(__lsx_vpickev_b, q1_filt8_l, q1_filt8_l, q2_filt8_l, q2_filt8_l,
                  q1_filt8_l, q2_filt8_l);

        /* store pixel values */
        p2 = __lsx_vbitsel_v(p2, p2_filt8_l, flat);
        p1 = __lsx_vbitsel_v(p1_out, p1_filt8_l, flat);
        p0 = __lsx_vbitsel_v(p0_out, p0_filt8_l, flat);
        q0 = __lsx_vbitsel_v(q0_out, q0_filt8_l, flat);
        q1 = __lsx_vbitsel_v(q1_out, q1_filt8_l, flat);
        q2 = __lsx_vbitsel_v(q2, q2_filt8_l, flat);

        DUP2_ARG2(__lsx_vilvl_b, p1, p2, q0, p0, vec0, vec1);
        vec3 = __lsx_vilvl_h(vec1, vec0);
        vec4 = __lsx_vilvh_h(vec1, vec0);
        DUP2_ARG2(__lsx_vilvh_b, p1, p2, q0, p0, vec0, vec1);
        vec6 = __lsx_vilvl_h(vec1, vec0);
        vec7 = __lsx_vilvh_h(vec1, vec0);
        vec2 = __lsx_vilvl_b(q2, q1);
        vec5 = __lsx_vilvh_b(q2, q1);

        dst -= 3;
        __lsx_vstelm_w(vec3, dst, 0, 0);
        __lsx_vstelm_h(vec2, dst, 4, 0);
        dst += stride;
        __lsx_vstelm_w(vec3, dst, 0, 1);
        __lsx_vstelm_h(vec2, dst, 4, 1);
        dst += stride;
        __lsx_vstelm_w(vec3, dst, 0, 2);
        __lsx_vstelm_h(vec2, dst, 4, 2);
        dst += stride;
        __lsx_vstelm_w(vec3, dst, 0, 3);
        __lsx_vstelm_h(vec2, dst, 4, 3);
        dst += stride;
        __lsx_vstelm_w(vec4, dst, 0, 0);
        __lsx_vstelm_h(vec2, dst, 4, 4);
        dst += stride;
        __lsx_vstelm_w(vec4, dst, 0, 1);
        __lsx_vstelm_h(vec2, dst, 4, 5);
        dst += stride;
        __lsx_vstelm_w(vec4, dst, 0, 2);
        __lsx_vstelm_h(vec2, dst, 4, 6);
        dst += stride;
        __lsx_vstelm_w(vec4, dst, 0, 3);
        __lsx_vstelm_h(vec2, dst, 4, 7);
        dst += stride;
        __lsx_vstelm_w(vec6, dst, 0, 0);
        __lsx_vstelm_h(vec5, dst, 4, 0);
        dst += stride;
        __lsx_vstelm_w(vec6, dst, 0, 1);
        __lsx_vstelm_h(vec5, dst, 4, 1);
        dst += stride;
        __lsx_vstelm_w(vec6, dst, 0, 2);
        __lsx_vstelm_h(vec5, dst, 4, 2);
        dst += stride;
        __lsx_vstelm_w(vec6, dst, 0, 3);
        __lsx_vstelm_h(vec5, dst, 4, 3);
        dst += stride;
        __lsx_vstelm_w(vec7, dst, 0, 0);
        __lsx_vstelm_h(vec5, dst, 4, 4);
        dst += stride;
        __lsx_vstelm_w(vec7, dst, 0, 1);
        __lsx_vstelm_h(vec5, dst, 4, 5);
        dst += stride;
        __lsx_vstelm_w(vec7, dst, 0, 2);
        __lsx_vstelm_h(vec5, dst, 4, 6);
        dst += stride;
        __lsx_vstelm_w(vec7, dst, 0, 3);
        __lsx_vstelm_h(vec5, dst, 4, 7);
    }
}

void ff_loop_filter_h_48_16_lsx(uint8_t *dst, ptrdiff_t stride,
                                int32_t b_limit_ptr,
                                int32_t limit_ptr,
                                int32_t thresh_ptr)
{
    ptrdiff_t stride2 = stride << 1;
    ptrdiff_t stride3 = stride2 + stride;
    ptrdiff_t stride4 = stride2 << 1;
    uint8_t *dst_tmp = dst - 4;
    __m128i p3, p2, p1, p0, q3, q2, q1, q0;
    __m128i p1_out, p0_out, q0_out, q1_out;
    __m128i flat, mask, hev, thresh, b_limit, limit;
    __m128i row4, row5, row6, row7, row12, row13, row14, row15;
    __m128i p3_h, p2_h, p1_h, p0_h, q0_h, q1_h, q2_h, q3_h;
    __m128i p2_filt8_h, p1_filt8_h, p0_filt8_h;
    __m128i q0_filt8_h, q1_filt8_h, q2_filt8_h;
    __m128i vec0, vec1, vec2, vec3, vec4, vec5, vec6, vec7;
    __m128i zero = __lsx_vldi(0);

    p0 = __lsx_vld(dst_tmp, 0);
    DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, p1, p2);
    p3 = __lsx_vldx(dst_tmp, stride3);
    dst_tmp += stride4;
    row4 = __lsx_vld(dst_tmp, 0);
    DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, row5, row6);
    row7 = __lsx_vldx(dst_tmp, stride3);
    dst_tmp += stride4;
    q3 = __lsx_vld(dst_tmp, 0);
    DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, q2, q1);
    q0 = __lsx_vldx(dst_tmp, stride3);
    dst_tmp += stride4;
    row12 = __lsx_vld(dst_tmp, 0);
    DUP2_ARG2(__lsx_vldx, dst_tmp, stride, dst_tmp, stride2, row13, row14);
    row15 = __lsx_vldx(dst_tmp, stride3);

    /* transpose 16x8 matrix into 8x16 */
    LSX_TRANSPOSE16x8_B(p0, p1, p2, p3, row4, row5, row6, row7,
                        q3, q2, q1, q0, row12, row13, row14, row15,
                        p3, p2, p1, p0, q0, q1, q2, q3);

    thresh = __lsx_vreplgr2vr_b(thresh_ptr);
    vec0   = __lsx_vreplgr2vr_b(thresh_ptr >> 8);
    thresh = __lsx_vilvl_d(vec0, thresh);

    b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
    vec0    = __lsx_vreplgr2vr_b(b_limit_ptr >> 8);
    b_limit = __lsx_vilvl_d(vec0, b_limit);

    limit = __lsx_vreplgr2vr_b(limit_ptr);
    vec0  = __lsx_vreplgr2vr_b(limit_ptr >> 8);
    limit = __lsx_vilvl_d(vec0, limit);

    /* mask and hev */
    LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
                 hev, mask, flat);
    /* flat4 */
    VP9_FLAT4(p3, p2, p0, q0, q2, q3, flat);
    /* filter4 */
    VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
                       q1_out);

    flat = __lsx_vilvh_d(flat, zero);

    /* if flat is zero for all pixels, then no need to calculate other filter */
    if (__lsx_bz_v(flat)) {
        DUP2_ARG2(__lsx_vilvl_b, p0_out, p1_out, q1_out, q0_out, vec0, vec1);
        vec2 = __lsx_vilvl_h(vec1, vec0);
        vec3 = __lsx_vilvh_h(vec1, vec0);
        DUP2_ARG2(__lsx_vilvh_b, p0_out, p1_out, q1_out, q0_out, vec0, vec1);
        vec4 = __lsx_vilvl_h(vec1, vec0);
        vec5 = __lsx_vilvh_h(vec1, vec0);

        dst -= 2;
        __lsx_vstelm_w(vec2, dst, 0, 0);
        __lsx_vstelm_w(vec2, dst + stride, 0, 1);
        __lsx_vstelm_w(vec2, dst + stride2, 0, 2);
        __lsx_vstelm_w(vec2, dst + stride3, 0, 3);
        dst += stride4;
        __lsx_vstelm_w(vec3, dst, 0, 0);
        __lsx_vstelm_w(vec3, dst + stride, 0, 1);
        __lsx_vstelm_w(vec3, dst + stride2, 0, 2);
        __lsx_vstelm_w(vec3, dst + stride3, 0, 3);
        dst += stride4;
        __lsx_vstelm_w(vec4, dst, 0, 0);
        __lsx_vstelm_w(vec4, dst + stride, 0, 1);
        __lsx_vstelm_w(vec4, dst + stride2, 0, 2);
        __lsx_vstelm_w(vec4, dst + stride3, 0, 3);
        dst += stride4;
        __lsx_vstelm_w(vec5, dst, 0, 0);
        __lsx_vstelm_w(vec5, dst + stride, 0, 1);
        __lsx_vstelm_w(vec5, dst + stride2, 0, 2);
        __lsx_vstelm_w(vec5, dst + stride3, 0, 3);
    } else {
        DUP4_ARG2(__lsx_vilvh_b, zero, p3, zero, p2, zero, p1, zero, p0,
                  p3_h, p2_h, p1_h, p0_h);
        DUP4_ARG2(__lsx_vilvh_b, zero, q0, zero, q1, zero, q2, zero, q3,
                  q0_h, q1_h, q2_h, q3_h);

        VP9_FILTER8(p3_h, p2_h, p1_h, p0_h, q0_h, q1_h, q2_h, q3_h, p2_filt8_h,
                    p1_filt8_h, p0_filt8_h, q0_filt8_h, q1_filt8_h, q2_filt8_h);

        /* convert 16 bit output data into 8 bit */
        DUP4_ARG2(__lsx_vpickev_b, p2_filt8_h, p2_filt8_h, p1_filt8_h,
                  p1_filt8_h, p0_filt8_h, p0_filt8_h, q0_filt8_h, q0_filt8_h,
                  p2_filt8_h, p1_filt8_h, p0_filt8_h, q0_filt8_h);
        DUP2_ARG2(__lsx_vpickev_b, q1_filt8_h, q1_filt8_h, q2_filt8_h,
                  q2_filt8_h, q1_filt8_h, q2_filt8_h);

        /* store pixel values */
        p2 = __lsx_vbitsel_v(p2, p2_filt8_h, flat);
        p1 = __lsx_vbitsel_v(p1_out, p1_filt8_h, flat);
        p0 = __lsx_vbitsel_v(p0_out, p0_filt8_h, flat);
        q0 = __lsx_vbitsel_v(q0_out, q0_filt8_h, flat);
        q1 = __lsx_vbitsel_v(q1_out, q1_filt8_h, flat);
        q2 = __lsx_vbitsel_v(q2, q2_filt8_h, flat);

        DUP2_ARG2(__lsx_vilvl_b, p1, p2, q0, p0, vec0, vec1);
        vec3 = __lsx_vilvl_h(vec1, vec0);
        vec4 = __lsx_vilvh_h(vec1, vec0);
        DUP2_ARG2(__lsx_vilvh_b, p1, p2, q0, p0, vec0, vec1);
        vec6 = __lsx_vilvl_h(vec1, vec0);
        vec7 = __lsx_vilvh_h(vec1, vec0);
        vec2 = __lsx_vilvl_b(q2, q1);
        vec5 = __lsx_vilvh_b(q2, q1);

        dst -= 3;
        __lsx_vstelm_w(vec3, dst, 0, 0);
        __lsx_vstelm_h(vec2, dst, 4, 0);
        dst += stride;
        __lsx_vstelm_w(vec3, dst, 0, 1);
        __lsx_vstelm_h(vec2, dst, 4, 1);
        dst += stride;
        __lsx_vstelm_w(vec3, dst, 0, 2);
        __lsx_vstelm_h(vec2, dst, 4, 2);
        dst += stride;
        __lsx_vstelm_w(vec3, dst, 0, 3);
        __lsx_vstelm_h(vec2, dst, 4, 3);
        dst += stride;
        __lsx_vstelm_w(vec4, dst, 0, 0);
        __lsx_vstelm_h(vec2, dst, 4, 4);
        dst += stride;
        __lsx_vstelm_w(vec4, dst, 0, 1);
        __lsx_vstelm_h(vec2, dst, 4, 5);
        dst += stride;
        __lsx_vstelm_w(vec4, dst, 0, 2);
        __lsx_vstelm_h(vec2, dst, 4, 6);
        dst += stride;
        __lsx_vstelm_w(vec4, dst, 0, 3);
        __lsx_vstelm_h(vec2, dst, 4, 7);
        dst += stride;
        __lsx_vstelm_w(vec6, dst, 0, 0);
        __lsx_vstelm_h(vec5, dst, 4, 0);
        dst += stride;
        __lsx_vstelm_w(vec6, dst, 0, 1);
        __lsx_vstelm_h(vec5, dst, 4, 1);
        dst += stride;
        __lsx_vstelm_w(vec6, dst, 0, 2);
        __lsx_vstelm_h(vec5, dst, 4, 2);
        dst += stride;
        __lsx_vstelm_w(vec6, dst, 0, 3);
        __lsx_vstelm_h(vec5, dst, 4, 3);
        dst += stride;
        __lsx_vstelm_w(vec7, dst, 0, 0);
        __lsx_vstelm_h(vec5, dst, 4, 4);
        dst += stride;
        __lsx_vstelm_w(vec7, dst, 0, 1);
        __lsx_vstelm_h(vec5, dst, 4, 5);
        dst += stride;
        __lsx_vstelm_w(vec7, dst, 0, 2);
        __lsx_vstelm_h(vec5, dst, 4, 6);
        dst += stride;
        __lsx_vstelm_w(vec7, dst, 0, 3);
        __lsx_vstelm_h(vec5, dst, 4, 7);
    }
}

static void vp9_transpose_16x8_to_8x16(uint8_t *input, ptrdiff_t in_pitch,
                                       uint8_t *output)
{
    __m128i p7_org, p6_org, p5_org, p4_org, p3_org, p2_org, p1_org, p0_org;
    __m128i tmp0, tmp1, tmp2, tmp3, tmp4, tmp5, tmp6, tmp7;
    __m128i p7, p6, p5, p4, p3, p2, p1, p0, q0, q1, q2, q3, q4, q5, q6, q7;
    ptrdiff_t in_pitch2 = in_pitch << 1;
    ptrdiff_t in_pitch3 = in_pitch2 + in_pitch;
    ptrdiff_t in_pitch4 = in_pitch2 << 1;

    LSX_LD_8(input, in_pitch, in_pitch2, in_pitch3, in_pitch4,
             p7_org, p6_org, p5_org, p4_org, p3_org, p2_org, p1_org, p0_org);
    /* 8x8 transpose */
    LSX_TRANSPOSE8x8_B(p7_org, p6_org, p5_org, p4_org, p3_org, p2_org, p1_org,
                       p0_org, p7, p6, p5, p4, p3, p2, p1, p0);
    /* 8x8 transpose */
    DUP4_ARG2(__lsx_vilvh_b, p5_org, p7_org, p4_org, p6_org, p1_org,
              p3_org, p0_org, p2_org, tmp0, tmp1, tmp2, tmp3);
    DUP2_ARG2(__lsx_vilvl_b, tmp1, tmp0, tmp3, tmp2, tmp4, tmp6);
    DUP2_ARG2(__lsx_vilvh_b, tmp1, tmp0, tmp3, tmp2, tmp5, tmp7);
    DUP2_ARG2(__lsx_vilvl_w, tmp6, tmp4, tmp7, tmp5, q0, q4);
    DUP2_ARG2(__lsx_vilvh_w, tmp6, tmp4, tmp7, tmp5, q2, q6);
    DUP4_ARG2(__lsx_vbsrl_v, q0, 8, q2, 8, q4, 8, q6, 8, q1, q3, q5, q7);

    __lsx_vst(p7, output, 0);
    __lsx_vst(p6, output, 16);
    __lsx_vst(p5, output, 32);
    __lsx_vst(p4, output, 48);
    __lsx_vst(p3, output, 64);
    __lsx_vst(p2, output, 80);
    __lsx_vst(p1, output, 96);
    __lsx_vst(p0, output, 112);
    __lsx_vst(q0, output, 128);
    __lsx_vst(q1, output, 144);
    __lsx_vst(q2, output, 160);
    __lsx_vst(q3, output, 176);
    __lsx_vst(q4, output, 192);
    __lsx_vst(q5, output, 208);
    __lsx_vst(q6, output, 224);
    __lsx_vst(q7, output, 240);
}

static void vp9_transpose_8x16_to_16x8(uint8_t *input, uint8_t *output,
                                       ptrdiff_t out_pitch)
{
    __m128i p7_o, p6_o, p5_o, p4_o, p3_o, p2_o, p1_o, p0_o;
    __m128i p7, p6, p5, p4, p3, p2, p1, p0, q0, q1, q2, q3, q4, q5, q6, q7;
    ptrdiff_t out_pitch2 = out_pitch << 1;
    ptrdiff_t out_pitch3 = out_pitch2 + out_pitch;
    ptrdiff_t out_pitch4 = out_pitch2 << 1;

    DUP4_ARG2(__lsx_vld, input, 0, input, 16, input, 32, input, 48,
              p7, p6, p5, p4);
    DUP4_ARG2(__lsx_vld, input, 64, input, 80, input, 96, input, 112,
              p3, p2, p1, p0);
    DUP4_ARG2(__lsx_vld, input, 128, input, 144, input, 160, input, 176,
              q0, q1, q2, q3);
    DUP4_ARG2(__lsx_vld, input, 192, input, 208, input, 224, input, 240,
              q4, q5, q6, q7);
    LSX_TRANSPOSE16x8_B(p7, p6, p5, p4, p3, p2, p1, p0, q0, q1, q2, q3, q4, q5,
                        q6, q7, p7_o, p6_o, p5_o, p4_o, p3_o, p2_o, p1_o, p0_o);
    LSX_ST_8(p7_o, p6_o, p5_o, p4_o, p3_o, p2_o, p1_o, p0_o,
             output, out_pitch, out_pitch2, out_pitch3, out_pitch4);
}

static void vp9_transpose_16x16(uint8_t *input, int32_t in_stride,
                                uint8_t *output, int32_t out_stride)
{
    __m128i row0, row1, row2, row3, row4, row5, row6, row7;
    __m128i row8, row9, row10, row11, row12, row13, row14, row15;
    __m128i tmp0, tmp1, tmp4, tmp5, tmp6, tmp7;
    __m128i tmp2, tmp3;
    __m128i p7, p6, p5, p4, p3, p2, p1, p0, q0, q1, q2, q3, q4, q5, q6, q7;
    int32_t in_stride2 = in_stride << 1;
    int32_t in_stride3 = in_stride2 + in_stride;
    int32_t in_stride4 = in_stride2 << 1;
    int32_t out_stride2 = out_stride << 1;
    int32_t out_stride3 = out_stride2 + out_stride;
    int32_t out_stride4 = out_stride2 << 1;

    LSX_LD_8(input, in_stride, in_stride2, in_stride3, in_stride4,
             row0, row1, row2, row3, row4, row5, row6, row7);
    input += in_stride4;
    LSX_LD_8(input, in_stride, in_stride2, in_stride3, in_stride4,
             row8, row9, row10, row11, row12, row13, row14, row15);

    LSX_TRANSPOSE16x8_B(row0, row1, row2, row3, row4, row5, row6, row7,
                        row8, row9, row10, row11, row12, row13, row14, row15,
                        p7, p6, p5, p4, p3, p2, p1, p0);

    /* transpose 16x8 matrix into 8x16 */
    /* total 8 intermediate register and 32 instructions */
    q7 = __lsx_vpackod_d(row8, row0);
    q6 = __lsx_vpackod_d(row9, row1);
    q5 = __lsx_vpackod_d(row10, row2);
    q4 = __lsx_vpackod_d(row11, row3);
    q3 = __lsx_vpackod_d(row12, row4);
    q2 = __lsx_vpackod_d(row13, row5);
    q1 = __lsx_vpackod_d(row14, row6);
    q0 = __lsx_vpackod_d(row15, row7);

    DUP2_ARG2(__lsx_vpackev_b, q6, q7, q4, q5, tmp0, tmp1);
    DUP2_ARG2(__lsx_vpackod_b, q6, q7, q4, q5, tmp4, tmp5);

    DUP2_ARG2(__lsx_vpackev_b, q2, q3, q0, q1, q5, q7);
    DUP2_ARG2(__lsx_vpackod_b, q2, q3, q0, q1, tmp6, tmp7);

    DUP2_ARG2(__lsx_vpackev_h, tmp1, tmp0, q7, q5, tmp2, tmp3);
    q0 = __lsx_vpackev_w(tmp3, tmp2);
    q4 = __lsx_vpackod_w(tmp3, tmp2);

    tmp2 = __lsx_vpackod_h(tmp1, tmp0);
    tmp3 = __lsx_vpackod_h(q7, q5);
    q2 = __lsx_vpackev_w(tmp3, tmp2);
    q6 = __lsx_vpackod_w(tmp3, tmp2);

    DUP2_ARG2(__lsx_vpackev_h, tmp5, tmp4, tmp7, tmp6, tmp2, tmp3);
    q1 = __lsx_vpackev_w(tmp3, tmp2);
    q5 = __lsx_vpackod_w(tmp3, tmp2);

    tmp2 = __lsx_vpackod_h(tmp5, tmp4);
    tmp3 = __lsx_vpackod_h(tmp7, tmp6);
    q3 = __lsx_vpackev_w(tmp3, tmp2);
    q7 = __lsx_vpackod_w(tmp3, tmp2);

    LSX_ST_8(p7, p6, p5, p4, p3, p2, p1, p0, output, out_stride,
             out_stride2, out_stride3, out_stride4);
    output += out_stride4;
    LSX_ST_8(q0, q1, q2, q3, q4, q5, q6, q7, output, out_stride,
             out_stride2, out_stride3, out_stride4);
}

static int32_t vp9_vt_lpf_t4_and_t8_8w(uint8_t *src, uint8_t *filter48,
                                       uint8_t *src_org, int32_t pitch_org,
                                       int32_t b_limit_ptr,
                                       int32_t limit_ptr,
                                       int32_t thresh_ptr)
{
    __m128i p3, p2, p1, p0, q3, q2, q1, q0;
    __m128i p2_out, p1_out, p0_out, q0_out, q1_out, q2_out;
    __m128i flat, mask, hev, thresh, b_limit, limit;
    __m128i p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l;
    __m128i p2_filt8_l, p1_filt8_l, p0_filt8_l;
    __m128i q0_filt8_l, q1_filt8_l, q2_filt8_l;
    __m128i vec0, vec1, vec2, vec3;
    __m128i zero = __lsx_vldi(0);

    /* load vector elements */
    DUP4_ARG2(__lsx_vld, src, -64, src, -48, src, -32, src, -16,
              p3, p2, p1, p0);
    DUP4_ARG2(__lsx_vld, src, 0, src, 16, src, 32, src, 48, q0, q1, q2, q3);

    thresh  = __lsx_vreplgr2vr_b(thresh_ptr);
    b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
    limit   = __lsx_vreplgr2vr_b(limit_ptr);

    /* mask and hev */
    LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
                 hev, mask, flat);
    /* flat4 */
    VP9_FLAT4(p3, p2, p0, q0, q2, q3, flat);
    /* filter4 */
    VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
                       q1_out);

    flat = __lsx_vilvl_d(zero, flat);

    /* if flat is zero for all pixels, then no need to calculate other filter */
    if (__lsx_bz_v(flat)) {
        DUP2_ARG2(__lsx_vilvl_b, p0_out, p1_out, q1_out, q0_out, vec0, vec1);
        vec2 = __lsx_vilvl_h(vec1, vec0);
        vec3 = __lsx_vilvh_h(vec1, vec0);

        src_org -= 2;
        __lsx_vstelm_w(vec2, src_org, 0, 0);
        src_org += pitch_org;
        __lsx_vstelm_w(vec2, src_org, 0, 1);
        src_org += pitch_org;
        __lsx_vstelm_w(vec2, src_org, 0, 2);
        src_org += pitch_org;
        __lsx_vstelm_w(vec2, src_org, 0, 3);
        src_org += pitch_org;
        __lsx_vstelm_w(vec3, src_org, 0, 0);
        src_org += pitch_org;
        __lsx_vstelm_w(vec3, src_org, 0, 1);
        src_org += pitch_org;
        __lsx_vstelm_w(vec3, src_org, 0, 2);
        src_org += pitch_org;
        __lsx_vstelm_w(vec3, src_org, 0, 3);
        return 1;
    } else {
        DUP4_ARG2(__lsx_vilvl_b, zero, p3, zero, p2, zero, p1, zero, p0,
                  p3_l, p2_l, p1_l, p0_l);
        DUP4_ARG2(__lsx_vilvl_b, zero, q0, zero, q1, zero, q2, zero, q3,
                  q0_l, q1_l, q2_l, q3_l);
        VP9_FILTER8(p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l, p2_filt8_l,
                    p1_filt8_l, p0_filt8_l, q0_filt8_l, q1_filt8_l, q2_filt8_l);

        /* convert 16 bit output data into 8 bit */
        p2_l = __lsx_vpickev_b(p2_filt8_l, p2_filt8_l);
        p1_l = __lsx_vpickev_b(p1_filt8_l, p1_filt8_l);
        p0_l = __lsx_vpickev_b(p0_filt8_l, p0_filt8_l);
        q0_l = __lsx_vpickev_b(q0_filt8_l, q0_filt8_l);
        q1_l = __lsx_vpickev_b(q1_filt8_l, q1_filt8_l);
        q2_l = __lsx_vpickev_b(q2_filt8_l, q2_filt8_l);

        /* store pixel values */
        p2_out = __lsx_vbitsel_v(p2, p2_l, flat);
        p1_out = __lsx_vbitsel_v(p1_out, p1_l, flat);
        p0_out = __lsx_vbitsel_v(p0_out, p0_l, flat);
        q0_out = __lsx_vbitsel_v(q0_out, q0_l, flat);
        q1_out = __lsx_vbitsel_v(q1_out, q1_l, flat);
        q2_out = __lsx_vbitsel_v(q2, q2_l, flat);

        __lsx_vst(p2_out, filter48, 0);
        __lsx_vst(p1_out, filter48, 16);
        __lsx_vst(p0_out, filter48, 32);
        __lsx_vst(q0_out, filter48, 48);
        __lsx_vst(q1_out, filter48, 64);
        __lsx_vst(q2_out, filter48, 80);
        __lsx_vst(flat, filter48, 96);

        return 0;
    }
}

static int32_t vp9_vt_lpf_t16_8w(uint8_t *dst, uint8_t *dst_org,
                                 ptrdiff_t stride,
                                 uint8_t *filter48)
{
    __m128i zero = __lsx_vldi(0);
    __m128i filter8, flat, flat2;
    __m128i p7, p6, p5, p4, p3, p2, p1, p0, q0, q1, q2, q3, q4, q5, q6, q7;
    v8u16 p7_l_in, p6_l_in, p5_l_in, p4_l_in;
    v8u16 p3_l_in, p2_l_in, p1_l_in, p0_l_in;
    v8u16 q7_l_in, q6_l_in, q5_l_in, q4_l_in;
    v8u16 q3_l_in, q2_l_in, q1_l_in, q0_l_in;
    v8u16 tmp0_l, tmp1_l;
    __m128i out_l;
    uint8_t *dst_tmp = dst - 128;

    /* load vector elements */
    DUP4_ARG2(__lsx_vld, dst_tmp, 0, dst_tmp, 16, dst_tmp, 32,
              dst_tmp, 48, p7, p6, p5, p4);
    DUP4_ARG2(__lsx_vld, dst_tmp, 64, dst_tmp, 80, dst_tmp, 96,
              dst_tmp, 112, p3, p2, p1, p0);
    DUP4_ARG2(__lsx_vld, dst, 0, dst, 16, dst, 32, dst, 48, q0, q1, q2, q3);
    DUP4_ARG2(__lsx_vld, dst, 64, dst, 80, dst, 96, dst, 112, q4, q5, q6, q7);

    flat = __lsx_vld(filter48, 96);


    VP9_FLAT5(p7, p6, p5, p4, p0, q0, q4, q5, q6, q7, flat, flat2);

    /* if flat2 is zero for all pixels, then no need to calculate other filter */
    if (__lsx_bz_v(flat2)) {
        __m128i vec0, vec1, vec2, vec3, vec4;

        DUP4_ARG2(__lsx_vld, filter48, 0, filter48, 16, filter48, 32,
                  filter48, 48, p2, p1, p0, q0);
        DUP2_ARG2(__lsx_vld, filter48, 64, filter48, 80, q1, q2);

        DUP2_ARG2(__lsx_vilvl_b, p1, p2, q0, p0, vec0, vec1);
        vec3 = __lsx_vilvl_h(vec1, vec0);
        vec4 = __lsx_vilvh_h(vec1, vec0);
        vec2 = __lsx_vilvl_b(q2, q1);

        dst_org -= 3;
        __lsx_vstelm_w(vec3, dst_org, 0, 0);
        __lsx_vstelm_h(vec2, dst_org, 4, 0);
        dst_org += stride;
        __lsx_vstelm_w(vec3, dst_org, 0, 1);
        __lsx_vstelm_h(vec2, dst_org, 4, 1);
        dst_org += stride;
        __lsx_vstelm_w(vec3, dst_org, 0, 2);
        __lsx_vstelm_h(vec2, dst_org, 4, 2);
        dst_org += stride;
        __lsx_vstelm_w(vec3, dst_org, 0, 3);
        __lsx_vstelm_h(vec2, dst_org, 4, 3);
        dst_org += stride;
        __lsx_vstelm_w(vec4, dst_org, 0, 0);
        __lsx_vstelm_h(vec2, dst_org, 4, 4);
        dst_org += stride;
        __lsx_vstelm_w(vec4, dst_org, 0, 1);
        __lsx_vstelm_h(vec2, dst_org, 4, 5);
        dst_org += stride;
        __lsx_vstelm_w(vec4, dst_org, 0, 2);
        __lsx_vstelm_h(vec2, dst_org, 4, 6);
        dst_org += stride;
        __lsx_vstelm_w(vec4, dst_org, 0, 3);
        __lsx_vstelm_h(vec2, dst_org, 4, 7);
        return 1;
    } else {
        dst -= 7 * 16;

        p7_l_in = (v8u16)__lsx_vilvl_b(zero, p7);
        p6_l_in = (v8u16)__lsx_vilvl_b(zero, p6);
        p5_l_in = (v8u16)__lsx_vilvl_b(zero, p5);
        p4_l_in = (v8u16)__lsx_vilvl_b(zero, p4);
        p3_l_in = (v8u16)__lsx_vilvl_b(zero, p3);
        p2_l_in = (v8u16)__lsx_vilvl_b(zero, p2);
        p1_l_in = (v8u16)__lsx_vilvl_b(zero, p1);
        p0_l_in = (v8u16)__lsx_vilvl_b(zero, p0);
        q0_l_in = (v8u16)__lsx_vilvl_b(zero, q0);

        tmp0_l = p7_l_in << 3;
        tmp0_l -= p7_l_in;
        tmp0_l += p6_l_in;
        tmp0_l += q0_l_in;
        tmp1_l = p6_l_in + p5_l_in;
        tmp1_l += p4_l_in;
        tmp1_l += p3_l_in;
        tmp1_l += p2_l_in;
        tmp1_l += p1_l_in;
        tmp1_l += p0_l_in;
        tmp1_l += tmp0_l;

        out_l =__lsx_vsrari_h((__m128i)tmp1_l, 4);
        out_l =__lsx_vpickev_b(out_l, out_l);
        p6 = __lsx_vbitsel_v(p6, out_l, flat2);
        __lsx_vstelm_d(p6, dst, 0, 0);
        dst += 16;

        /* p5 */
        q1_l_in = (v8u16)__lsx_vilvl_b(zero, q1);
        tmp0_l = p5_l_in - p6_l_in;
        tmp0_l += q1_l_in;
        tmp0_l -= p7_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
        out_l = __lsx_vpickev_b(out_l, out_l);
        p5 = __lsx_vbitsel_v(p5, out_l, flat2);
        __lsx_vstelm_d(p5, dst, 0, 0);
        dst += 16;

        /* p4 */
        q2_l_in = (v8u16)__lsx_vilvl_b(zero, q2);
        tmp0_l = p4_l_in - p5_l_in;
        tmp0_l += q2_l_in;
        tmp0_l -= p7_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
        out_l = __lsx_vpickev_b(out_l, out_l);
        p4 = __lsx_vbitsel_v(p4, out_l, flat2);
        __lsx_vstelm_d(p4, dst, 0, 0);
        dst += 16;

        /* p3 */
        q3_l_in = (v8u16)__lsx_vilvl_b(zero, q3);
        tmp0_l = p3_l_in - p4_l_in;
        tmp0_l += q3_l_in;
        tmp0_l -= p7_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
        out_l = __lsx_vpickev_b(out_l, out_l);
        p3 = __lsx_vbitsel_v(p3, out_l, flat2);
        __lsx_vstelm_d(p3, dst, 0, 0);
        dst += 16;

        /* p2 */
        q4_l_in = (v8u16)__lsx_vilvl_b(zero, q4);
        filter8 = __lsx_vld(filter48, 0);
        tmp0_l = p2_l_in - p3_l_in;
        tmp0_l += q4_l_in;
        tmp0_l -= p7_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
        out_l = __lsx_vpickev_b(out_l, out_l);
        filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
        __lsx_vstelm_d(filter8, dst, 0, 0);
        dst += 16;

        /* p1 */
        q5_l_in = (v8u16)__lsx_vilvl_b(zero, q5);
        filter8 = __lsx_vld(filter48, 16);
        tmp0_l = p1_l_in - p2_l_in;
        tmp0_l += q5_l_in;
        tmp0_l -= p7_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
        out_l = __lsx_vpickev_b(out_l, out_l);
        filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
        __lsx_vstelm_d(filter8, dst, 0, 0);
        dst += 16;

        /* p0 */
        q6_l_in = (v8u16)__lsx_vilvl_b(zero, q6);
        filter8 = __lsx_vld(filter48, 32);
        tmp0_l = p0_l_in - p1_l_in;
        tmp0_l += q6_l_in;
        tmp0_l -= p7_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
        out_l = __lsx_vpickev_b(out_l, out_l);
        filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
        __lsx_vstelm_d(filter8, dst, 0, 0);
        dst += 16;

        /* q0 */
        q7_l_in = (v8u16)__lsx_vilvl_b(zero, q7);
        filter8 = __lsx_vld(filter48, 48);
        tmp0_l = q7_l_in - p0_l_in;
        tmp0_l += q0_l_in;
        tmp0_l -= p7_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((v8i16) tmp1_l, 4);
        out_l = __lsx_vpickev_b(out_l, out_l);
        filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
        __lsx_vstelm_d(filter8, dst, 0, 0);
        dst += 16;

        /* q1 */
        filter8 = __lsx_vld(filter48, 64);
        tmp0_l = q7_l_in - q0_l_in;
        tmp0_l += q1_l_in;
        tmp0_l -= p6_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
        out_l = __lsx_vpickev_b(out_l, out_l);
        filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
        __lsx_vstelm_d(filter8, dst, 0, 0);
        dst += 16;

        /* q2 */
        filter8 = __lsx_vld(filter48, 80);
        tmp0_l = q7_l_in - q1_l_in;
        tmp0_l += q2_l_in;
        tmp0_l -= p5_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
        out_l = __lsx_vpickev_b(out_l, out_l);
        filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
        __lsx_vstelm_d(filter8, dst, 0, 0);
        dst += 16;

        /* q3 */
        tmp0_l = q7_l_in - q2_l_in;
        tmp0_l += q3_l_in;
        tmp0_l -= p4_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
        out_l = __lsx_vpickev_b(out_l, out_l);
        q3 = __lsx_vbitsel_v(q3, out_l, flat2);
        __lsx_vstelm_d(q3, dst, 0, 0);
        dst += 16;

        /* q4 */
        tmp0_l = q7_l_in - q3_l_in;
        tmp0_l += q4_l_in;
        tmp0_l -= p3_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
        out_l = __lsx_vpickev_b(out_l, out_l);
        q4 = __lsx_vbitsel_v(q4, out_l, flat2);
        __lsx_vstelm_d(q4, dst, 0, 0);
        dst += 16;

        /* q5 */
        tmp0_l = q7_l_in - q4_l_in;
        tmp0_l += q5_l_in;
        tmp0_l -= p2_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
        out_l = __lsx_vpickev_b(out_l, out_l);
        q5 = __lsx_vbitsel_v(q5, out_l, flat2);
        __lsx_vstelm_d(q5, dst, 0, 0);
        dst += 16;

        /* q6 */
        tmp0_l = q7_l_in - q5_l_in;
        tmp0_l += q6_l_in;
        tmp0_l -= p1_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
        out_l = __lsx_vpickev_b(out_l, out_l);
        q6 = __lsx_vbitsel_v(q6, out_l, flat2);
        __lsx_vstelm_d(q6, dst, 0, 0);

        return 0;
    }
}

void ff_loop_filter_h_16_8_lsx(uint8_t *dst, ptrdiff_t stride,
                               int32_t b_limit_ptr,
                               int32_t limit_ptr,
                               int32_t thresh_ptr)
{
    uint8_t early_exit = 0;
    uint8_t transposed_input[16 * 24] __attribute__ ((aligned(16)));
    uint8_t *filter48 = &transposed_input[16 * 16];

    vp9_transpose_16x8_to_8x16(dst - 8, stride, transposed_input);

    early_exit = vp9_vt_lpf_t4_and_t8_8w((transposed_input + 16 * 8),
                                         &filter48[0], dst, stride,
                                         b_limit_ptr, limit_ptr, thresh_ptr);

    if (0 == early_exit) {
        early_exit = vp9_vt_lpf_t16_8w((transposed_input + 16 * 8), dst, stride,
                                       &filter48[0]);

        if (0 == early_exit) {
            vp9_transpose_8x16_to_16x8(transposed_input, dst - 8, stride);
        }
    }
}

static int32_t vp9_vt_lpf_t4_and_t8_16w(uint8_t *dst, uint8_t *filter48,
                                        uint8_t *dst_org, ptrdiff_t stride,
                                        int32_t b_limit_ptr,
                                        int32_t limit_ptr,
                                        int32_t thresh_ptr)
{
    ptrdiff_t stride2 = stride << 1;
    ptrdiff_t stride3 = stride2 + stride;
    ptrdiff_t stride4 = stride2 << 1;
    __m128i p3, p2, p1, p0, q3, q2, q1, q0;
    __m128i p2_out, p1_out, p0_out, q0_out, q1_out, q2_out;
    __m128i flat, mask, hev, thresh, b_limit, limit;
    __m128i p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l;
    __m128i p3_h, p2_h, p1_h, p0_h, q0_h, q1_h, q2_h, q3_h;
    __m128i p2_filt8_l, p1_filt8_l, p0_filt8_l;
    __m128i q0_filt8_l, q1_filt8_l, q2_filt8_l;
    __m128i p2_filt8_h, p1_filt8_h, p0_filt8_h;
    __m128i q0_filt8_h, q1_filt8_h, q2_filt8_h;
    __m128i vec0, vec1, vec2, vec3, vec4, vec5;
    __m128i zero = __lsx_vldi(0);

    /* load vector elements */
    DUP4_ARG2(__lsx_vld, dst, -64, dst, -48, dst, -32, dst, -16,
              p3, p2, p1, p0);
    DUP4_ARG2(__lsx_vld, dst, 0, dst, 16, dst, 32, dst, 48, q0, q1, q2, q3);

    thresh  = __lsx_vreplgr2vr_b(thresh_ptr);
    b_limit = __lsx_vreplgr2vr_b(b_limit_ptr);
    limit   = __lsx_vreplgr2vr_b(limit_ptr);

    /* mask and hev */
    LPF_MASK_HEV(p3, p2, p1, p0, q0, q1, q2, q3, limit, b_limit, thresh,
                 hev, mask, flat);
    /* flat4 */
    VP9_FLAT4(p3, p2, p0, q0, q2, q3, flat);
    /* filter4 */
    VP9_LPF_FILTER4_4W(p1, p0, q0, q1, mask, hev, p1_out, p0_out, q0_out,
                       q1_out);

    /* if flat is zero for all pixels, then no need to calculate other filter */
    if (__lsx_bz_v(flat)) {
        DUP2_ARG2(__lsx_vilvl_b, p0_out, p1_out, q1_out, q0_out, vec0, vec1);
        vec2 = __lsx_vilvl_h(vec1, vec0);
        vec3 = __lsx_vilvh_h(vec1, vec0);
        DUP2_ARG2(__lsx_vilvh_b, p0_out, p1_out, q1_out, q0_out, vec0, vec1);
        vec4 = __lsx_vilvl_h(vec1, vec0);
        vec5 = __lsx_vilvh_h(vec1, vec0);

        dst_org -= 2;
        __lsx_vstelm_w(vec2, dst_org, 0, 0);
        __lsx_vstelm_w(vec2, dst_org + stride, 0, 1);
        __lsx_vstelm_w(vec2, dst_org + stride2, 0, 2);
        __lsx_vstelm_w(vec2, dst_org + stride3, 0, 3);
        dst_org += stride4;
        __lsx_vstelm_w(vec3, dst_org, 0, 0);
        __lsx_vstelm_w(vec3, dst_org + stride, 0, 1);
        __lsx_vstelm_w(vec3, dst_org + stride2, 0, 2);
        __lsx_vstelm_w(vec3, dst_org + stride3, 0, 3);
        dst_org += stride4;
        __lsx_vstelm_w(vec4, dst_org, 0, 0);
        __lsx_vstelm_w(vec4, dst_org + stride, 0, 1);
        __lsx_vstelm_w(vec4, dst_org + stride2, 0, 2);
        __lsx_vstelm_w(vec4, dst_org + stride3, 0, 3);
        dst_org += stride4;
        __lsx_vstelm_w(vec5, dst_org, 0, 0);
        __lsx_vstelm_w(vec5, dst_org + stride, 0, 1);
        __lsx_vstelm_w(vec5, dst_org + stride2, 0, 2);
        __lsx_vstelm_w(vec5, dst_org + stride3, 0, 3);

        return 1;
    } else {
        DUP4_ARG2(__lsx_vilvl_b, zero, p3, zero, p2, zero, p1, zero, p0,
                  p3_l, p2_l, p1_l, p0_l);
        DUP4_ARG2(__lsx_vilvl_b, zero, q0, zero, q1, zero, q2, zero, q3,
                  q0_l, q1_l, q2_l, q3_l);
        VP9_FILTER8(p3_l, p2_l, p1_l, p0_l, q0_l, q1_l, q2_l, q3_l, p2_filt8_l,
                    p1_filt8_l, p0_filt8_l, q0_filt8_l, q1_filt8_l, q2_filt8_l);
        DUP4_ARG2(__lsx_vilvh_b, zero, p3, zero, p2, zero, p1, zero, p0,
                  p3_h, p2_h, p1_h, p0_h);
        DUP4_ARG2(__lsx_vilvh_b, zero, q0, zero, q1, zero, q2, zero, q3,
                  q0_h, q1_h, q2_h, q3_h);
        VP9_FILTER8(p3_h, p2_h, p1_h, p0_h, q0_h, q1_h, q2_h, q3_h, p2_filt8_h,
                    p1_filt8_h, p0_filt8_h, q0_filt8_h, q1_filt8_h, q2_filt8_h);

        /* convert 16 bit output data into 8 bit */
        DUP4_ARG2(__lsx_vpickev_b, p2_filt8_h, p2_filt8_l, p1_filt8_h,
                      p1_filt8_l, p0_filt8_h, p0_filt8_l, q0_filt8_h,
                      q0_filt8_l, p2_filt8_l, p1_filt8_l, p0_filt8_l,
                      q0_filt8_l);
        DUP2_ARG2(__lsx_vpickev_b, q1_filt8_h, q1_filt8_l, q2_filt8_h,
                  q2_filt8_l, q1_filt8_l, q2_filt8_l);

        /* store pixel values */
        p2_out = __lsx_vbitsel_v(p2, p2_filt8_l, flat);
        p1_out = __lsx_vbitsel_v(p1_out, p1_filt8_l, flat);
        p0_out = __lsx_vbitsel_v(p0_out, p0_filt8_l, flat);
        q0_out = __lsx_vbitsel_v(q0_out, q0_filt8_l, flat);
        q1_out = __lsx_vbitsel_v(q1_out, q1_filt8_l, flat);
        q2_out = __lsx_vbitsel_v(q2, q2_filt8_l, flat);

        __lsx_vst(p2_out, filter48, 0);
        __lsx_vst(p1_out, filter48, 16);
        __lsx_vst(p0_out, filter48, 32);
        __lsx_vst(q0_out, filter48, 48);
        __lsx_vst(q1_out, filter48, 64);
        __lsx_vst(q2_out, filter48, 80);
        __lsx_vst(flat, filter48, 96);

        return 0;
    }
}

static int32_t vp9_vt_lpf_t16_16w(uint8_t *dst, uint8_t *dst_org,
                                  ptrdiff_t stride,
                                  uint8_t *filter48)
{
    __m128i zero = __lsx_vldi(0);
    __m128i flat, flat2, filter8;
    __m128i p7, p6, p5, p4, p3, p2, p1, p0, q0, q1, q2, q3, q4, q5, q6, q7;
    v8u16 p7_l_in, p6_l_in, p5_l_in, p4_l_in;
    v8u16 p3_l_in, p2_l_in, p1_l_in, p0_l_in;
    v8u16 q7_l_in, q6_l_in, q5_l_in, q4_l_in;
    v8u16 q3_l_in, q2_l_in, q1_l_in, q0_l_in;
    v8u16 p7_h_in, p6_h_in, p5_h_in, p4_h_in;
    v8u16 p3_h_in, p2_h_in, p1_h_in, p0_h_in;
    v8u16 q7_h_in, q6_h_in, q5_h_in, q4_h_in;
    v8u16 q3_h_in, q2_h_in, q1_h_in, q0_h_in;
    v8u16 tmp0_l, tmp1_l, tmp0_h, tmp1_h;
    __m128i out_l, out_h;
    uint8_t *dst_tmp = dst - 128;

    flat = __lsx_vld(filter48, 96);

    DUP4_ARG2(__lsx_vld, dst_tmp, 0, dst_tmp, 16, dst_tmp, 32,
              dst_tmp, 48, p7, p6, p5, p4);
    DUP4_ARG2(__lsx_vld, dst_tmp, 64, dst_tmp, 80, dst_tmp, 96,
              dst_tmp, 112, p3, p2, p1, p0);
    DUP4_ARG2(__lsx_vld, dst, 0, dst, 16, dst, 32, dst, 48, q0, q1, q2, q3);
    DUP4_ARG2(__lsx_vld, dst, 64, dst, 80, dst, 96, dst, 112, q4, q5, q6, q7);

    VP9_FLAT5(p7, p6, p5, p4, p0, q0, q4, q5, q6, q7, flat, flat2);

    /* if flat2 is zero for all pixels, then no need to calculate other filter */
    if (__lsx_bz_v(flat2)) {
        __m128i vec0, vec1, vec2, vec3, vec4, vec5, vec6, vec7;

        DUP4_ARG2(__lsx_vld, filter48, 0, filter48, 16, filter48, 32,
                  filter48, 48, p2, p1, p0, q0);
        DUP2_ARG2(__lsx_vld, filter48, 64, filter48, 80, q1, q2);

        DUP2_ARG2(__lsx_vilvl_b, p1, p2, q0, p0, vec0, vec1);
        vec3 = __lsx_vilvl_h(vec1, vec0);
        vec4 = __lsx_vilvh_h(vec1, vec0);
        DUP2_ARG2(__lsx_vilvh_b, p1, p2, q0, p0, vec0, vec1);
        vec6 = __lsx_vilvl_h(vec1, vec0);
        vec7 = __lsx_vilvh_h(vec1, vec0);
        vec2 = __lsx_vilvl_b(q2, q1);
        vec5 = __lsx_vilvh_b(q2, q1);

        dst_org -= 3;
        __lsx_vstelm_w(vec3, dst_org, 0, 0);
        __lsx_vstelm_h(vec2, dst_org, 4, 0);
        dst_org += stride;
        __lsx_vstelm_w(vec3, dst_org, 0, 1);
        __lsx_vstelm_h(vec2, dst_org, 4, 1);
        dst_org += stride;
        __lsx_vstelm_w(vec3, dst_org, 0, 2);
        __lsx_vstelm_h(vec2, dst_org, 4, 2);
        dst_org += stride;
        __lsx_vstelm_w(vec3, dst_org, 0, 3);
        __lsx_vstelm_h(vec2, dst_org, 4, 3);
        dst_org += stride;
        __lsx_vstelm_w(vec4, dst_org, 0, 0);
        __lsx_vstelm_h(vec2, dst_org, 4, 4);
        dst_org += stride;
        __lsx_vstelm_w(vec4, dst_org, 0, 1);
        __lsx_vstelm_h(vec2, dst_org, 4, 5);
        dst_org += stride;
        __lsx_vstelm_w(vec4, dst_org, 0, 2);
        __lsx_vstelm_h(vec2, dst_org, 4, 6);
        dst_org += stride;
        __lsx_vstelm_w(vec4, dst_org, 0, 3);
        __lsx_vstelm_h(vec2, dst_org, 4, 7);
        dst_org += stride;
        __lsx_vstelm_w(vec6, dst_org, 0, 0);
        __lsx_vstelm_h(vec5, dst_org, 4, 0);
        dst_org += stride;
        __lsx_vstelm_w(vec6, dst_org, 0, 1);
        __lsx_vstelm_h(vec5, dst_org, 4, 1);
        dst_org += stride;
        __lsx_vstelm_w(vec6, dst_org, 0, 2);
        __lsx_vstelm_h(vec5, dst_org, 4, 2);
        dst_org += stride;
        __lsx_vstelm_w(vec6, dst_org, 0, 3);
        __lsx_vstelm_h(vec5, dst_org, 4, 3);
        dst_org += stride;
        __lsx_vstelm_w(vec7, dst_org, 0, 0);
        __lsx_vstelm_h(vec5, dst_org, 4, 4);
        dst_org += stride;
        __lsx_vstelm_w(vec7, dst_org, 0, 1);
        __lsx_vstelm_h(vec5, dst_org, 4, 5);
        dst_org += stride;
        __lsx_vstelm_w(vec7, dst_org, 0, 2);
        __lsx_vstelm_h(vec5, dst_org, 4, 6);
        dst_org += stride;
        __lsx_vstelm_w(vec7, dst_org, 0, 3);
        __lsx_vstelm_h(vec5, dst_org, 4, 7);

        return 1;
    } else {
        dst -= 7 * 16;

        p7_l_in = (v8u16)__lsx_vilvl_b(zero, p7);
        p6_l_in = (v8u16)__lsx_vilvl_b(zero, p6);
        p5_l_in = (v8u16)__lsx_vilvl_b(zero, p5);
        p4_l_in = (v8u16)__lsx_vilvl_b(zero, p4);
        p3_l_in = (v8u16)__lsx_vilvl_b(zero, p3);
        p2_l_in = (v8u16)__lsx_vilvl_b(zero, p2);
        p1_l_in = (v8u16)__lsx_vilvl_b(zero, p1);
        p0_l_in = (v8u16)__lsx_vilvl_b(zero, p0);
        q0_l_in = (v8u16)__lsx_vilvl_b(zero, q0);

        tmp0_l = p7_l_in << 3;
        tmp0_l -= p7_l_in;
        tmp0_l += p6_l_in;
        tmp0_l += q0_l_in;
        tmp1_l = p6_l_in + p5_l_in;
        tmp1_l += p4_l_in;
        tmp1_l += p3_l_in;
        tmp1_l += p2_l_in;
        tmp1_l += p1_l_in;
        tmp1_l += p0_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);

        p7_h_in = (v8u16)__lsx_vilvh_b(zero, p7);
        p6_h_in = (v8u16)__lsx_vilvh_b(zero, p6);
        p5_h_in = (v8u16)__lsx_vilvh_b(zero, p5);
        p4_h_in = (v8u16)__lsx_vilvh_b(zero, p4);
        p3_h_in = (v8u16)__lsx_vilvh_b(zero, p3);
        p2_h_in = (v8u16)__lsx_vilvh_b(zero, p2);
        p1_h_in = (v8u16)__lsx_vilvh_b(zero, p1);
        p0_h_in = (v8u16)__lsx_vilvh_b(zero, p0);
        q0_h_in = (v8u16)__lsx_vilvh_b(zero, q0);

        tmp0_h = p7_h_in << 3;
        tmp0_h -= p7_h_in;
        tmp0_h += p6_h_in;
        tmp0_h += q0_h_in;
        tmp1_h = p6_h_in + p5_h_in;
        tmp1_h += p4_h_in;
        tmp1_h += p3_h_in;
        tmp1_h += p2_h_in;
        tmp1_h += p1_h_in;
        tmp1_h += p0_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);

        out_l = __lsx_vpickev_b(out_h, out_l);
        p6 = __lsx_vbitsel_v(p6, out_l, flat2);
        __lsx_vst(p6, dst, 0);

        /* p5 */
        q1_l_in = (v8u16)__lsx_vilvl_b(zero, q1);
        tmp0_l = p5_l_in - p6_l_in;
        tmp0_l += q1_l_in;
        tmp0_l -= p7_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
        q1_h_in = (v8u16)__lsx_vilvh_b(zero, q1);
        tmp0_h = p5_h_in - p6_h_in;
        tmp0_h += q1_h_in;
        tmp0_h -= p7_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
        out_l = __lsx_vpickev_b(out_h, out_l);
        p5 = __lsx_vbitsel_v(p5, out_l, flat2);
        __lsx_vst(p5, dst, 16);

        /* p4 */
        q2_l_in = (v8u16)__lsx_vilvl_b(zero, q2);
        tmp0_l = p4_l_in - p5_l_in;
        tmp0_l += q2_l_in;
        tmp0_l -= p7_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
        q2_h_in = (v8u16)__lsx_vilvh_b(zero, q2);
        tmp0_h = p4_h_in - p5_h_in;
        tmp0_h += q2_h_in;
        tmp0_h -= p7_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
        out_l = __lsx_vpickev_b(out_h, out_l);
        p4 = __lsx_vbitsel_v(p4, out_l, flat2);
        __lsx_vst(p4, dst, 16*2);

        /* p3 */
        q3_l_in = (v8u16)__lsx_vilvl_b(zero, q3);
        tmp0_l = p3_l_in - p4_l_in;
        tmp0_l += q3_l_in;
        tmp0_l -= p7_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
        q3_h_in = (v8u16)__lsx_vilvh_b(zero, q3);
        tmp0_h = p3_h_in - p4_h_in;
        tmp0_h += q3_h_in;
        tmp0_h -= p7_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
        out_l = __lsx_vpickev_b(out_h, out_l);
        p3 = __lsx_vbitsel_v(p3, out_l, flat2);
        __lsx_vst(p3, dst, 16*3);

        /* p2 */
        q4_l_in = (v8u16)__lsx_vilvl_b(zero, q4);
        filter8 = __lsx_vld(filter48, 0);
        tmp0_l = p2_l_in - p3_l_in;
        tmp0_l += q4_l_in;
        tmp0_l -= p7_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
        q4_h_in = (v8u16)__lsx_vilvh_b(zero, q4);
        tmp0_h = p2_h_in - p3_h_in;
        tmp0_h += q4_h_in;
        tmp0_h -= p7_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
        out_l = __lsx_vpickev_b(out_h, out_l);
        filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
        __lsx_vst(filter8, dst, 16*4);

        /* p1 */
        q5_l_in = (v8u16)__lsx_vilvl_b(zero, q5);
        filter8 = __lsx_vld(filter48, 16);
        tmp0_l = p1_l_in - p2_l_in;
        tmp0_l += q5_l_in;
        tmp0_l -= p7_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
        q5_h_in = (v8u16)__lsx_vilvh_b(zero, q5);
        tmp0_h = p1_h_in - p2_h_in;
        tmp0_h += q5_h_in;
        tmp0_h -= p7_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)(tmp1_h), 4);
        out_l = __lsx_vpickev_b(out_h, out_l);
        filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
        __lsx_vst(filter8, dst, 16*5);

        /* p0 */
        q6_l_in = (v8u16)__lsx_vilvl_b(zero, q6);
        filter8 = __lsx_vld(filter48, 32);
        tmp0_l = p0_l_in - p1_l_in;
        tmp0_l += q6_l_in;
        tmp0_l -= p7_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
        q6_h_in = (v8u16)__lsx_vilvh_b(zero, q6);
        tmp0_h = p0_h_in - p1_h_in;
        tmp0_h += q6_h_in;
        tmp0_h -= p7_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
        out_l = __lsx_vpickev_b(out_h, out_l);
        filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
        __lsx_vst(filter8, dst, 16*6);

        /* q0 */
        q7_l_in = (v8u16)__lsx_vilvl_b(zero, q7);
        filter8 = __lsx_vld(filter48, 48);
        tmp0_l = q7_l_in - p0_l_in;
        tmp0_l += q0_l_in;
        tmp0_l -= p7_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
        q7_h_in = (v8u16)__lsx_vilvh_b(zero, q7);
        tmp0_h = q7_h_in - p0_h_in;
        tmp0_h += q0_h_in;
        tmp0_h -= p7_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
        out_l = __lsx_vpickev_b(out_h, out_l);
        filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
        __lsx_vst(filter8, dst, 16*7);

        /* q1 */
        filter8 = __lsx_vld(filter48, 64);
        tmp0_l = q7_l_in - q0_l_in;
        tmp0_l += q1_l_in;
        tmp0_l -= p6_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
        tmp0_h = q7_h_in - q0_h_in;
        tmp0_h += q1_h_in;
        tmp0_h -= p6_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
        out_l = __lsx_vpickev_b(out_h, out_l);
        filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
        __lsx_vst(filter8, dst, 16*8);

        /* q2 */
        filter8 = __lsx_vld(filter48, 80);
        tmp0_l = q7_l_in - q1_l_in;
        tmp0_l += q2_l_in;
        tmp0_l -= p5_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
        tmp0_h = q7_h_in - q1_h_in;
        tmp0_h += q2_h_in;
        tmp0_h -= p5_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
        out_l = __lsx_vpickev_b(out_h, out_l);
        filter8 = __lsx_vbitsel_v(filter8, out_l, flat2);
        __lsx_vst(filter8, dst, 16*9);

        /* q3 */
        tmp0_l = q7_l_in - q2_l_in;
        tmp0_l += q3_l_in;
        tmp0_l -= p4_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
        tmp0_h = q7_h_in - q2_h_in;
        tmp0_h += q3_h_in;
        tmp0_h -= p4_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
        out_l = __lsx_vpickev_b(out_h, out_l);
        q3 = __lsx_vbitsel_v(q3, out_l, flat2);
        __lsx_vst(q3, dst, 16*10);

        /* q4 */
        tmp0_l = q7_l_in - q3_l_in;
        tmp0_l += q4_l_in;
        tmp0_l -= p3_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
        tmp0_h = q7_h_in - q3_h_in;
        tmp0_h += q4_h_in;
        tmp0_h -= p3_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
        out_l = __lsx_vpickev_b(out_h, out_l);
        q4 = __lsx_vbitsel_v(q4, out_l, flat2);
        __lsx_vst(q4, dst, 16*11);

        /* q5 */
        tmp0_l = q7_l_in - q4_l_in;
        tmp0_l += q5_l_in;
        tmp0_l -= p2_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
        tmp0_h = q7_h_in - q4_h_in;
        tmp0_h += q5_h_in;
        tmp0_h -= p2_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
        out_l = __lsx_vpickev_b(out_h, out_l);
        q5 = __lsx_vbitsel_v(q5, out_l, flat2);
        __lsx_vst(q5, dst, 16*12);

        /* q6 */
        tmp0_l = q7_l_in - q5_l_in;
        tmp0_l += q6_l_in;
        tmp0_l -= p1_l_in;
        tmp1_l += tmp0_l;
        out_l = __lsx_vsrari_h((__m128i)tmp1_l, 4);
        tmp0_h = q7_h_in - q5_h_in;
        tmp0_h += q6_h_in;
        tmp0_h -= p1_h_in;
        tmp1_h += tmp0_h;
        out_h = __lsx_vsrari_h((__m128i)tmp1_h, 4);
        out_l = __lsx_vpickev_b(out_h, out_l);
        q6 = __lsx_vbitsel_v(q6, out_l, flat2);
        __lsx_vst(q6, dst, 16*13);

        return 0;
    }
}

void ff_loop_filter_h_16_16_lsx(uint8_t *dst, ptrdiff_t stride,
                                int32_t b_limit_ptr,
                                int32_t limit_ptr,
                                int32_t thresh_ptr)
{
    uint8_t early_exit = 0;
    uint8_t transposed_input[16 * 24] __attribute__ ((aligned(16)));
    uint8_t *filter48 = &transposed_input[16 * 16];

    vp9_transpose_16x16((dst - 8), stride, &transposed_input[0], 16);

    early_exit = vp9_vt_lpf_t4_and_t8_16w((transposed_input + 16 * 8),
                                          &filter48[0], dst, stride,
                                          b_limit_ptr, limit_ptr, thresh_ptr);

    if (0 == early_exit) {
        early_exit = vp9_vt_lpf_t16_16w((transposed_input + 16 * 8), dst,
                                         stride, &filter48[0]);

        if (0 == early_exit) {
            vp9_transpose_16x16(transposed_input, 16, (dst - 8), stride);
        }
    }
}