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

#include "libavutil/mips/generic_macros_msa.h"
#include "libavcodec/mips/hevcdsp_mips.h"

#define HEVC_FILT_8TAP_DPADD_W(vec0, vec1, vec2, vec3,            \
                               filt0, filt1, filt2, filt3)        \
( {                                                               \
    v4i32 out;                                                    \
                                                                  \
    out = __msa_dotp_s_w((v8i16) (vec0), (v8i16) (filt0));        \
    out = __msa_dpadd_s_w(out, (v8i16) (vec1), (v8i16) (filt1));  \
    out = __msa_dpadd_s_w(out, (v8i16) (vec2), (v8i16) (filt2));  \
    out = __msa_dpadd_s_w(out, (v8i16) (vec3), (v8i16) (filt3));  \
    out;                                                          \
} )

#define HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,                         \
                               filt0, filt1, filt2, filt3,                     \
                               var_in)                                         \
( {                                                                            \
    v8i16 out;                                                                 \
                                                                               \
    out = __msa_dpadd_s_h((v8i16) (var_in), (v16i8) (vec0), (v16i8) (filt0));  \
    out = __msa_dpadd_s_h(out, (v16i8) (vec1), (v16i8) (filt1));               \
    out = __msa_dpadd_s_h(out, (v16i8) (vec2), (v16i8) (filt2));               \
    out = __msa_dpadd_s_h(out, (v16i8) (vec3), (v16i8) (filt3));               \
    out;                                                                       \
} )

static void hevc_copy_4w_msa(uint8_t * __restrict src, int32_t src_stride,
                             int16_t * __restrict dst, int32_t dst_stride,
                             int32_t height)
{
    v16i8 zero = { 0 };

    if (2 == height) {
        uint64_t out0, out1;
        v16i8 src0, src1;
        v8i16 input0;

        LOAD_2VECS_SB(src, src_stride, src0, src1);

        src0 = (v16i8) __msa_ilvr_w((v4i32) src1, (v4i32) src0);

        input0 = (v8i16) __msa_ilvr_b(zero, src0);

        input0 <<= 6;

        out0 = __msa_copy_u_d((v2i64) input0, 0);
        out1 = __msa_copy_u_d((v2i64) input0, 1);

        STORE_DWORD(dst, out0);
        dst += dst_stride;
        STORE_DWORD(dst, out1);
    } else if (4 == height) {
        uint64_t out0, out1, out2, out3;
        v16i8 src0, src1, src2, src3;
        v8i16 input0, input1;

        LOAD_4VECS_SB(src, src_stride, src0, src1, src2, src3);

        src0 = (v16i8) __msa_ilvr_w((v4i32) src1, (v4i32) src0);
        src1 = (v16i8) __msa_ilvr_w((v4i32) src3, (v4i32) src2);

        input0 = (v8i16) __msa_ilvr_b(zero, src0);
        input1 = (v8i16) __msa_ilvr_b(zero, src1);

        input0 <<= 6;
        input1 <<= 6;

        out0 = __msa_copy_u_d((v2i64) input0, 0);
        out1 = __msa_copy_u_d((v2i64) input0, 1);
        out2 = __msa_copy_u_d((v2i64) input1, 0);
        out3 = __msa_copy_u_d((v2i64) input1, 1);

        STORE_DWORD(dst, out0);
        dst += dst_stride;
        STORE_DWORD(dst, out1);
        dst += dst_stride;
        STORE_DWORD(dst, out2);
        dst += dst_stride;
        STORE_DWORD(dst, out3);
    } else if (0 == height % 8) {
        v16i8 src0, src1, src2, src3, src4, src5, src6, src7;
        uint64_t out0, out1, out2, out3, out4, out5, out6, out7;
        v8i16 input0, input1, input2, input3;
        uint32_t loop_cnt;

        for (loop_cnt = (height >> 3); loop_cnt--;) {
            LOAD_8VECS_SB(src, src_stride,
                          src0, src1, src2, src3, src4, src5, src6, src7);
            src += (8 * src_stride);

            src0 = (v16i8) __msa_ilvr_w((v4i32) src1, (v4i32) src0);
            src1 = (v16i8) __msa_ilvr_w((v4i32) src3, (v4i32) src2);
            src2 = (v16i8) __msa_ilvr_w((v4i32) src5, (v4i32) src4);
            src3 = (v16i8) __msa_ilvr_w((v4i32) src7, (v4i32) src6);

            input0 = (v8i16) __msa_ilvr_b(zero, src0);
            input1 = (v8i16) __msa_ilvr_b(zero, src1);
            input2 = (v8i16) __msa_ilvr_b(zero, src2);
            input3 = (v8i16) __msa_ilvr_b(zero, src3);

            input0 <<= 6;
            input1 <<= 6;
            input2 <<= 6;
            input3 <<= 6;

            out0 = __msa_copy_u_d((v2i64) input0, 0);
            out1 = __msa_copy_u_d((v2i64) input0, 1);
            out2 = __msa_copy_u_d((v2i64) input1, 0);
            out3 = __msa_copy_u_d((v2i64) input1, 1);
            out4 = __msa_copy_u_d((v2i64) input2, 0);
            out5 = __msa_copy_u_d((v2i64) input2, 1);
            out6 = __msa_copy_u_d((v2i64) input3, 0);
            out7 = __msa_copy_u_d((v2i64) input3, 1);

            STORE_DWORD(dst, out0);
            dst += dst_stride;
            STORE_DWORD(dst, out1);
            dst += dst_stride;
            STORE_DWORD(dst, out2);
            dst += dst_stride;
            STORE_DWORD(dst, out3);
            dst += dst_stride;
            STORE_DWORD(dst, out4);
            dst += dst_stride;
            STORE_DWORD(dst, out5);
            dst += dst_stride;
            STORE_DWORD(dst, out6);
            dst += dst_stride;
            STORE_DWORD(dst, out7);
            dst += dst_stride;
        }
    }
}

static void hevc_copy_6w_msa(uint8_t * __restrict src, int32_t src_stride,
                             int16_t * __restrict dst, int32_t dst_stride,
                             int32_t height)
{
    uint32_t loop_cnt;
    uint64_t out0, out1, out2, out3, out4, out5, out6, out7;
    uint32_t out8, out9, out10, out11, out12, out13, out14, out15;
    v16i8 zero = { 0 };
    v16i8 src0, src1, src2, src3, src4, src5, src6, src7;
    v8i16 input0, input1, input2, input3, input4, input5, input6, input7;

    for (loop_cnt = (height >> 3); loop_cnt--;) {
        LOAD_8VECS_SB(src, src_stride,
                      src0, src1, src2, src3, src4, src5, src6, src7);
        src += (8 * src_stride);

        input0 = (v8i16) __msa_ilvr_b(zero, src0);
        input1 = (v8i16) __msa_ilvr_b(zero, src1);
        input2 = (v8i16) __msa_ilvr_b(zero, src2);
        input3 = (v8i16) __msa_ilvr_b(zero, src3);
        input4 = (v8i16) __msa_ilvr_b(zero, src4);
        input5 = (v8i16) __msa_ilvr_b(zero, src5);
        input6 = (v8i16) __msa_ilvr_b(zero, src6);
        input7 = (v8i16) __msa_ilvr_b(zero, src7);

        input0 <<= 6;
        input1 <<= 6;
        input2 <<= 6;
        input3 <<= 6;
        input4 <<= 6;
        input5 <<= 6;
        input6 <<= 6;
        input7 <<= 6;

        out0 = __msa_copy_u_d((v2i64) input0, 0);
        out1 = __msa_copy_u_d((v2i64) input1, 0);
        out2 = __msa_copy_u_d((v2i64) input2, 0);
        out3 = __msa_copy_u_d((v2i64) input3, 0);
        out4 = __msa_copy_u_d((v2i64) input4, 0);
        out5 = __msa_copy_u_d((v2i64) input5, 0);
        out6 = __msa_copy_u_d((v2i64) input6, 0);
        out7 = __msa_copy_u_d((v2i64) input7, 0);

        out8 =  __msa_copy_u_w((v4i32) input0, 2);
        out9 =  __msa_copy_u_w((v4i32) input1, 2);
        out10 = __msa_copy_u_w((v4i32) input2, 2);
        out11 = __msa_copy_u_w((v4i32) input3, 2);
        out12 = __msa_copy_u_w((v4i32) input4, 2);
        out13 = __msa_copy_u_w((v4i32) input5, 2);
        out14 = __msa_copy_u_w((v4i32) input6, 2);
        out15 = __msa_copy_u_w((v4i32) input7, 2);

        STORE_DWORD(dst, out0);
        STORE_WORD(dst + 4, out8);
        dst += dst_stride;
        STORE_DWORD(dst, out1);
        STORE_WORD(dst + 4, out9);
        dst += dst_stride;
        STORE_DWORD(dst, out2);
        STORE_WORD(dst + 4, out10);
        dst += dst_stride;
        STORE_DWORD(dst, out3);
        STORE_WORD(dst + 4, out11);
        dst += dst_stride;
        STORE_DWORD(dst, out4);
        STORE_WORD(dst + 4, out12);
        dst += dst_stride;
        STORE_DWORD(dst, out5);
        STORE_WORD(dst + 4, out13);
        dst += dst_stride;
        STORE_DWORD(dst, out6);
        STORE_WORD(dst + 4, out14);
        dst += dst_stride;
        STORE_DWORD(dst, out7);
        STORE_WORD(dst + 4, out15);
        dst += dst_stride;
    }
}

static void hevc_copy_8w_msa(uint8_t * __restrict src, int32_t src_stride,
                             int16_t * __restrict dst, int32_t dst_stride,
                             int32_t height)
{
    v16i8 zero = { 0 };

    if (2 == height) {
        v16i8 src0, src1;
        v8i16 input0, input1;

        LOAD_2VECS_SB(src, src_stride, src0, src1);

        input0 = (v8i16) __msa_ilvr_b(zero, src0);
        input1 = (v8i16) __msa_ilvr_b(zero, src1);

        input0 <<= 6;
        input1 <<= 6;

        STORE_2VECS_SH(dst, dst_stride, input0, input1);
    } else if (4 == height) {
        v16i8 src0, src1, src2, src3;
        v8i16 input0, input1, input2, input3;

        LOAD_4VECS_SB(src, src_stride, src0, src1, src2, src3);

        input0 = (v8i16) __msa_ilvr_b(zero, src0);
        input1 = (v8i16) __msa_ilvr_b(zero, src1);
        input2 = (v8i16) __msa_ilvr_b(zero, src2);
        input3 = (v8i16) __msa_ilvr_b(zero, src3);

        input0 <<= 6;
        input1 <<= 6;
        input2 <<= 6;
        input3 <<= 6;

        STORE_4VECS_SH(dst, dst_stride, input0, input1, input2, input3);
    } else if (6 == height) {
        v16i8 src0, src1, src2, src3, src4, src5;
        v8i16 input0, input1, input2, input3, input4, input5;

        LOAD_6VECS_SB(src, src_stride, src0, src1, src2, src3, src4, src5);

        input0 = (v8i16) __msa_ilvr_b(zero, src0);
        input1 = (v8i16) __msa_ilvr_b(zero, src1);
        input2 = (v8i16) __msa_ilvr_b(zero, src2);
        input3 = (v8i16) __msa_ilvr_b(zero, src3);
        input4 = (v8i16) __msa_ilvr_b(zero, src4);
        input5 = (v8i16) __msa_ilvr_b(zero, src5);

        input0 <<= 6;
        input1 <<= 6;
        input2 <<= 6;
        input3 <<= 6;
        input4 <<= 6;
        input5 <<= 6;

        STORE_6VECS_SH(dst, dst_stride,
                       input0, input1, input2, input3, input4, input5);
    } else if (0 == height % 8) {
        uint32_t loop_cnt;
        v16i8 src0, src1, src2, src3, src4, src5, src6, src7;
        v8i16 input0, input1, input2, input3;
        v8i16 input4, input5, input6, input7;

        for (loop_cnt = (height >> 3); loop_cnt--;) {
            LOAD_8VECS_SB(src, src_stride,
                          src0, src1, src2, src3, src4, src5, src6, src7);
            src += (8 * src_stride);

            input0 = (v8i16) __msa_ilvr_b(zero, src0);
            input1 = (v8i16) __msa_ilvr_b(zero, src1);
            input2 = (v8i16) __msa_ilvr_b(zero, src2);
            input3 = (v8i16) __msa_ilvr_b(zero, src3);
            input4 = (v8i16) __msa_ilvr_b(zero, src4);
            input5 = (v8i16) __msa_ilvr_b(zero, src5);
            input6 = (v8i16) __msa_ilvr_b(zero, src6);
            input7 = (v8i16) __msa_ilvr_b(zero, src7);

            input0 <<= 6;
            input1 <<= 6;
            input2 <<= 6;
            input3 <<= 6;
            input4 <<= 6;
            input5 <<= 6;
            input6 <<= 6;
            input7 <<= 6;

            STORE_8VECS_SH(dst, dst_stride,
                           input0, input1, input2, input3,
                           input4, input5, input6, input7);
            dst += (8 * dst_stride);
        }
    }
}

static void hevc_copy_12w_msa(uint8_t * __restrict src, int32_t src_stride,
                              int16_t * __restrict dst, int32_t dst_stride,
                              int32_t height)
{
    uint32_t loop_cnt;
    uint64_t dst_val0, dst_val1, dst_val2, dst_val3;
    v16i8 zero = { 0 };
    v16i8 src0, src1, src2, src3, src4, src5, src6, src7;
    v8i16 input0, input1;
    v8i16 input0_r, input1_r, input2_r, input3_r;

    for (loop_cnt = (height >> 3); loop_cnt--;) {
        LOAD_8VECS_SB(src, src_stride,
                      src0, src1, src2, src3, src4, src5, src6, src7);
        src += (8 * src_stride);

        input0_r = (v8i16) __msa_ilvr_b(zero, src0);
        input1_r = (v8i16) __msa_ilvr_b(zero, src1);
        input2_r = (v8i16) __msa_ilvr_b(zero, src2);
        input3_r = (v8i16) __msa_ilvr_b(zero, src3);

        input0_r <<= 6;
        input1_r <<= 6;
        input2_r <<= 6;
        input3_r <<= 6;

        src0 = (v16i8) __msa_ilvl_w((v4i32) src1, (v4i32) src0);
        src1 = (v16i8) __msa_ilvl_w((v4i32) src3, (v4i32) src2);

        input0 = (v8i16) __msa_ilvr_b(zero, src0);
        input1 = (v8i16) __msa_ilvr_b(zero, src1);

        input0 <<= 6;
        input1 <<= 6;

        dst_val0 = __msa_copy_u_d((v2i64) input0, 0);
        dst_val1 = __msa_copy_u_d((v2i64) input0, 1);
        dst_val2 = __msa_copy_u_d((v2i64) input1, 0);
        dst_val3 = __msa_copy_u_d((v2i64) input1, 1);

        STORE_4VECS_SH(dst, dst_stride, input0_r, input1_r, input2_r, input3_r);

        STORE_DWORD(dst + 8, dst_val0);
        dst += dst_stride;
        STORE_DWORD(dst + 8, dst_val1);
        dst += dst_stride;
        STORE_DWORD(dst + 8, dst_val2);
        dst += dst_stride;
        STORE_DWORD(dst + 8, dst_val3);
        dst += dst_stride;

        input0_r = (v8i16) __msa_ilvr_b(zero, src4);
        input1_r = (v8i16) __msa_ilvr_b(zero, src5);
        input2_r = (v8i16) __msa_ilvr_b(zero, src6);
        input3_r = (v8i16) __msa_ilvr_b(zero, src7);

        input0_r <<= 6;
        input1_r <<= 6;
        input2_r <<= 6;
        input3_r <<= 6;

        src0 = (v16i8) __msa_ilvl_w((v4i32) src5, (v4i32) src4);
        src1 = (v16i8) __msa_ilvl_w((v4i32) src7, (v4i32) src6);

        input0 = (v8i16) __msa_ilvr_b(zero, src0);
        input1 = (v8i16) __msa_ilvr_b(zero, src1);

        input0 <<= 6;
        input1 <<= 6;

        dst_val0 = __msa_copy_u_d((v2i64) input0, 0);
        dst_val1 = __msa_copy_u_d((v2i64) input0, 1);
        dst_val2 = __msa_copy_u_d((v2i64) input1, 0);
        dst_val3 = __msa_copy_u_d((v2i64) input1, 1);

        STORE_4VECS_SH(dst, dst_stride, input0_r, input1_r, input2_r, input3_r);

        STORE_DWORD(dst + 8, dst_val0);
        dst += dst_stride;
        STORE_DWORD(dst + 8, dst_val1);
        dst += dst_stride;
        STORE_DWORD(dst + 8, dst_val2);
        dst += dst_stride;
        STORE_DWORD(dst + 8, dst_val3);
        dst += dst_stride;
    }
}

static void hevc_copy_16multx8mult_msa(uint8_t * __restrict src,
                                       int32_t src_stride,
                                       int16_t * __restrict dst,
                                       int32_t dst_stride,
                                       int32_t height,
                                       int32_t width)
{
    uint8_t *src_tmp;
    int16_t *dst_tmp;
    uint32_t loop_cnt, cnt;
    v16i8 zero = { 0 };
    v16i8 src0, src1, src2, src3, src4, src5, src6, src7;
    v8i16 input0_r, input1_r, input2_r, input3_r;
    v8i16 input0_l, input1_l, input2_l, input3_l;

    for (cnt = (width >> 4); cnt--;) {
        src_tmp = src;
        dst_tmp = dst;

        for (loop_cnt = (height >> 3); loop_cnt--;) {
            LOAD_8VECS_SB(src_tmp, src_stride,
                          src0, src1, src2, src3, src4, src5, src6, src7);
            src_tmp += (8 * src_stride);

            input0_r = (v8i16) __msa_ilvr_b(zero, src0);
            input0_l = (v8i16) __msa_ilvl_b(zero, src0);
            input1_r = (v8i16) __msa_ilvr_b(zero, src1);
            input1_l = (v8i16) __msa_ilvl_b(zero, src1);
            input2_r = (v8i16) __msa_ilvr_b(zero, src2);
            input2_l = (v8i16) __msa_ilvl_b(zero, src2);
            input3_r = (v8i16) __msa_ilvr_b(zero, src3);
            input3_l = (v8i16) __msa_ilvl_b(zero, src3);

            input0_r <<= 6;
            input0_l <<= 6;
            input1_r <<= 6;
            input1_l <<= 6;
            input2_r <<= 6;
            input2_l <<= 6;
            input3_r <<= 6;
            input3_l <<= 6;

            STORE_4VECS_SH(dst_tmp, dst_stride,
                           input0_r, input1_r, input2_r, input3_r);
            STORE_4VECS_SH((dst_tmp + 8), dst_stride,
                           input0_l, input1_l, input2_l, input3_l);
            dst_tmp += (4 * dst_stride);

            input0_r = (v8i16) __msa_ilvr_b(zero, src4);
            input0_l = (v8i16) __msa_ilvl_b(zero, src4);
            input1_r = (v8i16) __msa_ilvr_b(zero, src5);
            input1_l = (v8i16) __msa_ilvl_b(zero, src5);
            input2_r = (v8i16) __msa_ilvr_b(zero, src6);
            input2_l = (v8i16) __msa_ilvl_b(zero, src6);
            input3_r = (v8i16) __msa_ilvr_b(zero, src7);
            input3_l = (v8i16) __msa_ilvl_b(zero, src7);

            input0_r <<= 6;
            input0_l <<= 6;
            input1_r <<= 6;
            input1_l <<= 6;
            input2_r <<= 6;
            input2_l <<= 6;
            input3_r <<= 6;
            input3_l <<= 6;

            STORE_4VECS_SH(dst_tmp, dst_stride,
                           input0_r, input1_r, input2_r, input3_r);
            STORE_4VECS_SH((dst_tmp + 8), dst_stride,
                           input0_l, input1_l, input2_l, input3_l);
            dst_tmp += (4 * dst_stride);
        }

        src += 16;
        dst += 16;
    }
}

static void hevc_copy_16w_msa(uint8_t * __restrict src, int32_t src_stride,
                              int16_t * __restrict dst, int32_t dst_stride,
                              int32_t height)
{
    v16i8 zero = { 0 };

    if (4 == height) {
        v16i8 src0, src1, src2, src3;
        v8i16 input0_r, input1_r, input2_r, input3_r;
        v8i16 input0_l, input1_l, input2_l, input3_l;

        LOAD_4VECS_SB(src, src_stride, src0, src1, src2, src3);

        input0_r = (v8i16) __msa_ilvr_b(zero, src0);
        input0_l = (v8i16) __msa_ilvl_b(zero, src0);
        input1_r = (v8i16) __msa_ilvr_b(zero, src1);
        input1_l = (v8i16) __msa_ilvl_b(zero, src1);
        input2_r = (v8i16) __msa_ilvr_b(zero, src2);
        input2_l = (v8i16) __msa_ilvl_b(zero, src2);
        input3_r = (v8i16) __msa_ilvr_b(zero, src3);
        input3_l = (v8i16) __msa_ilvl_b(zero, src3);

        input0_r <<= 6;
        input0_l <<= 6;
        input1_r <<= 6;
        input1_l <<= 6;
        input2_r <<= 6;
        input2_l <<= 6;
        input3_r <<= 6;
        input3_l <<= 6;

        STORE_4VECS_SH(dst, dst_stride, input0_r, input1_r, input2_r, input3_r);
        STORE_4VECS_SH((dst + 8), dst_stride,
                       input0_l, input1_l, input2_l, input3_l);
    } else if (12 == height) {
        v16i8 src0, src1, src2, src3, src4, src5, src6, src7;
        v16i8 src8, src9, src10, src11;
        v8i16 input0_r, input1_r, input2_r, input3_r;
        v8i16 input0_l, input1_l, input2_l, input3_l;

        LOAD_8VECS_SB(src, src_stride,
                      src0, src1, src2, src3, src4, src5, src6, src7);
        src += (8 * src_stride);

        LOAD_4VECS_SB(src, src_stride, src8, src9, src10, src11);

        input0_r = (v8i16) __msa_ilvr_b(zero, src0);
        input0_l = (v8i16) __msa_ilvl_b(zero, src0);
        input1_r = (v8i16) __msa_ilvr_b(zero, src1);
        input1_l = (v8i16) __msa_ilvl_b(zero, src1);
        input2_r = (v8i16) __msa_ilvr_b(zero, src2);
        input2_l = (v8i16) __msa_ilvl_b(zero, src2);
        input3_r = (v8i16) __msa_ilvr_b(zero, src3);
        input3_l = (v8i16) __msa_ilvl_b(zero, src3);

        input0_r <<= 6;
        input0_l <<= 6;
        input1_r <<= 6;
        input1_l <<= 6;
        input2_r <<= 6;
        input2_l <<= 6;
        input3_r <<= 6;
        input3_l <<= 6;

        STORE_4VECS_SH(dst, dst_stride, input0_r, input1_r, input2_r, input3_r);
        STORE_4VECS_SH((dst + 8), dst_stride,
                       input0_l, input1_l, input2_l, input3_l);
        dst += (4 * dst_stride);

        input0_r = (v8i16) __msa_ilvr_b(zero, src4);
        input0_l = (v8i16) __msa_ilvl_b(zero, src4);
        input1_r = (v8i16) __msa_ilvr_b(zero, src5);
        input1_l = (v8i16) __msa_ilvl_b(zero, src5);
        input2_r = (v8i16) __msa_ilvr_b(zero, src6);
        input2_l = (v8i16) __msa_ilvl_b(zero, src6);
        input3_r = (v8i16) __msa_ilvr_b(zero, src7);
        input3_l = (v8i16) __msa_ilvl_b(zero, src7);

        input0_r <<= 6;
        input0_l <<= 6;
        input1_r <<= 6;
        input1_l <<= 6;
        input2_r <<= 6;
        input2_l <<= 6;
        input3_r <<= 6;
        input3_l <<= 6;

        STORE_4VECS_SH(dst, dst_stride, input0_r, input1_r, input2_r, input3_r);
        STORE_4VECS_SH((dst + 8), dst_stride,
                       input0_l, input1_l, input2_l, input3_l);
        dst += (4 * dst_stride);

        input0_r = (v8i16) __msa_ilvr_b(zero, src8);
        input0_l = (v8i16) __msa_ilvl_b(zero, src8);
        input1_r = (v8i16) __msa_ilvr_b(zero, src9);
        input1_l = (v8i16) __msa_ilvl_b(zero, src9);
        input2_r = (v8i16) __msa_ilvr_b(zero, src10);
        input2_l = (v8i16) __msa_ilvl_b(zero, src10);
        input3_r = (v8i16) __msa_ilvr_b(zero, src11);
        input3_l = (v8i16) __msa_ilvl_b(zero, src11);

        input0_r <<= 6;
        input0_l <<= 6;
        input1_r <<= 6;
        input1_l <<= 6;
        input2_r <<= 6;
        input2_l <<= 6;
        input3_r <<= 6;
        input3_l <<= 6;

        STORE_4VECS_SH(dst, dst_stride, input0_r, input1_r, input2_r, input3_r);
        STORE_4VECS_SH((dst + 8), dst_stride,
                       input0_l, input1_l, input2_l, input3_l);
    } else if (0 == (height % 8)) {
        hevc_copy_16multx8mult_msa(src, src_stride, dst, dst_stride,
                                   height, 16);
    }
}

static void hevc_copy_24w_msa(uint8_t * __restrict src, int32_t src_stride,
                              int16_t * __restrict dst, int32_t dst_stride,
                              int32_t height)
{
    hevc_copy_16multx8mult_msa(src, src_stride, dst, dst_stride, height, 16);

    hevc_copy_8w_msa(src + 16, src_stride, dst + 16, dst_stride, height);
}

static void hevc_copy_32w_msa(uint8_t * __restrict src, int32_t src_stride,
                              int16_t * __restrict dst, int32_t dst_stride,
                              int32_t height)
{
    hevc_copy_16multx8mult_msa(src, src_stride, dst, dst_stride, height, 32);
}

static void hevc_copy_48w_msa(uint8_t * __restrict src, int32_t src_stride,
                              int16_t * __restrict dst, int32_t dst_stride,
                              int32_t height)
{
    hevc_copy_16multx8mult_msa(src, src_stride, dst, dst_stride, height, 48);
}

static void hevc_copy_64w_msa(uint8_t * __restrict src, int32_t src_stride,
                              int16_t * __restrict dst, int32_t dst_stride,
                              int32_t height)
{
    hevc_copy_16multx8mult_msa(src, src_stride, dst, dst_stride, height, 64);
}

static void hevc_hz_8t_4w_msa(uint8_t * __restrict src, int32_t src_stride,
                              int16_t * __restrict dst, int32_t dst_stride,
                              const int8_t * __restrict filter, int32_t height)
{
    uint32_t loop_cnt;
    uint64_t out0, out1, out2, out3, out4, out5, out6, out7;
    v16i8 src0, src1, src2, src3, src4, src5, src6, src7;
    v8i16 filt0, filt1, filt2, filt3;
    v16i8 mask1, mask2, mask3;
    v8u16 const_vec;
    v16i8 vec0, vec1, vec2, vec3;
    v8i16 dst0, dst1, dst2, dst3;
    v8i16 filter_vec;
    v16i8 mask0 = { 0, 1, 1, 2, 2, 3, 3, 4, 16, 17, 17, 18, 18, 19, 19, 20 };

    src -= 3;

    const_vec = (v8u16) __msa_ldi_h(128);
    const_vec <<= 6;

    filter_vec = LOAD_SH(filter);
    filt0 = __msa_splati_h(filter_vec, 0);
    filt1 = __msa_splati_h(filter_vec, 1);
    filt2 = __msa_splati_h(filter_vec, 2);
    filt3 = __msa_splati_h(filter_vec, 3);

    mask1 = mask0 + 2;
    mask2 = mask0 + 4;
    mask3 = mask0 + 6;

    for (loop_cnt = (height >> 3); loop_cnt--;) {
        LOAD_8VECS_SB(src, src_stride,
                      src0, src1, src2, src3, src4, src5, src6, src7);
        src += (8 * src_stride);

        XORI_B_8VECS_SB(src0, src1, src2, src3, src4, src5, src6, src7,
                        src0, src1, src2, src3, src4, src5, src6, src7, 128);

        vec0 = __msa_vshf_b(mask0, src1, src0);
        vec1 = __msa_vshf_b(mask1, src1, src0);
        vec2 = __msa_vshf_b(mask2, src1, src0);
        vec3 = __msa_vshf_b(mask3, src1, src0);

        dst0 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask0, src3, src2);
        vec1 = __msa_vshf_b(mask1, src3, src2);
        vec2 = __msa_vshf_b(mask2, src3, src2);
        vec3 = __msa_vshf_b(mask3, src3, src2);

        dst1 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask0, src5, src4);
        vec1 = __msa_vshf_b(mask1, src5, src4);
        vec2 = __msa_vshf_b(mask2, src5, src4);
        vec3 = __msa_vshf_b(mask3, src5, src4);

        dst2 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask0, src7, src6);
        vec1 = __msa_vshf_b(mask1, src7, src6);
        vec2 = __msa_vshf_b(mask2, src7, src6);
        vec3 = __msa_vshf_b(mask3, src7, src6);

        dst3 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        out0 = __msa_copy_u_d((v2i64) dst0, 0);
        out1 = __msa_copy_u_d((v2i64) dst0, 1);
        out2 = __msa_copy_u_d((v2i64) dst1, 0);
        out3 = __msa_copy_u_d((v2i64) dst1, 1);
        out4 = __msa_copy_u_d((v2i64) dst2, 0);
        out5 = __msa_copy_u_d((v2i64) dst2, 1);
        out6 = __msa_copy_u_d((v2i64) dst3, 0);
        out7 = __msa_copy_u_d((v2i64) dst3, 1);

        STORE_DWORD(dst, out0);
        dst += dst_stride;
        STORE_DWORD(dst, out1);
        dst += dst_stride;
        STORE_DWORD(dst, out2);
        dst += dst_stride;
        STORE_DWORD(dst, out3);
        dst += dst_stride;
        STORE_DWORD(dst, out4);
        dst += dst_stride;
        STORE_DWORD(dst, out5);
        dst += dst_stride;
        STORE_DWORD(dst, out6);
        dst += dst_stride;
        STORE_DWORD(dst, out7);
        dst += dst_stride;
    }
}

static void hevc_hz_8t_8w_msa(uint8_t * __restrict src, int32_t src_stride,
                              int16_t * __restrict dst, int32_t dst_stride,
                              const int8_t * __restrict filter, int32_t height)
{
    uint32_t loop_cnt;
    v16i8 src0, src1, src2, src3;
    v8i16 filt0, filt1, filt2, filt3;
    v16i8 mask1, mask2, mask3;
    v8u16 const_vec;
    v16i8 vec0, vec1, vec2, vec3;
    v8i16 dst0, dst1, dst2, dst3;
    v8i16 filter_vec;
    v16i8 mask0 = { 0, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 7, 8 };

    src -= 3;

    const_vec = (v8u16) __msa_ldi_h(128);
    const_vec <<= 6;

    filter_vec = LOAD_SH(filter);
    filt0 = __msa_splati_h(filter_vec, 0);
    filt1 = __msa_splati_h(filter_vec, 1);
    filt2 = __msa_splati_h(filter_vec, 2);
    filt3 = __msa_splati_h(filter_vec, 3);

    mask1 = mask0 + 2;
    mask2 = mask0 + 4;
    mask3 = mask0 + 6;

    for (loop_cnt = (height >> 2); loop_cnt--;) {
        LOAD_4VECS_SB(src, src_stride, src0, src1, src2, src3);
        src += (4 * src_stride);

        XORI_B_4VECS_SB(src0, src1, src2, src3, src0, src1, src2, src3, 128);

        vec0 = __msa_vshf_b(mask0, src0, src0);
        vec1 = __msa_vshf_b(mask1, src0, src0);
        vec2 = __msa_vshf_b(mask2, src0, src0);
        vec3 = __msa_vshf_b(mask3, src0, src0);

        dst0 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask0, src1, src1);
        vec1 = __msa_vshf_b(mask1, src1, src1);
        vec2 = __msa_vshf_b(mask2, src1, src1);
        vec3 = __msa_vshf_b(mask3, src1, src1);

        dst1 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask0, src2, src2);
        vec1 = __msa_vshf_b(mask1, src2, src2);
        vec2 = __msa_vshf_b(mask2, src2, src2);
        vec3 = __msa_vshf_b(mask3, src2, src2);

        dst2 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask0, src3, src3);
        vec1 = __msa_vshf_b(mask1, src3, src3);
        vec2 = __msa_vshf_b(mask2, src3, src3);
        vec3 = __msa_vshf_b(mask3, src3, src3);

        dst3 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        STORE_SH(dst0, dst);
        dst += dst_stride;
        STORE_SH(dst1, dst);
        dst += dst_stride;
        STORE_SH(dst2, dst);
        dst += dst_stride;
        STORE_SH(dst3, dst);
        dst += dst_stride;
    }
}

static void hevc_hz_8t_12w_msa(uint8_t * __restrict src, int32_t src_stride,
                               int16_t * __restrict dst, int32_t dst_stride,
                               const int8_t * __restrict filter, int32_t height)
{
    hevc_hz_8t_8w_msa(src, src_stride, dst, dst_stride, filter, height);

    hevc_hz_8t_4w_msa(src + 8, src_stride, dst + 8, dst_stride, filter, height);
}

static void hevc_hz_8t_16w_msa(uint8_t * __restrict src, int32_t src_stride,
                               int16_t * __restrict dst, int32_t dst_stride,
                               const int8_t * __restrict filter, int32_t height)
{
    uint32_t loop_cnt;
    v16i8 src0, src1, src2, src3, src4, src5, src6, src7;
    v8i16 filt0, filt1, filt2, filt3;
    v16i8 mask1, mask2, mask3;
    v8u16 const_vec;
    v16i8 vec0, vec1, vec2, vec3;
    v8i16 dst0, dst1, dst2, dst3, dst4, dst5, dst6, dst7;
    v8i16 filter_vec;
    v16i8 mask0 = { 0, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 7, 8 };

    src -= 3;

    const_vec = (v8u16) __msa_ldi_h(128);
    const_vec <<= 6;

    filter_vec = LOAD_SH(filter);
    filt0 = __msa_splati_h(filter_vec, 0);
    filt1 = __msa_splati_h(filter_vec, 1);
    filt2 = __msa_splati_h(filter_vec, 2);
    filt3 = __msa_splati_h(filter_vec, 3);

    mask1 = mask0 + 2;
    mask2 = mask0 + 4;
    mask3 = mask0 + 6;

    for (loop_cnt = (height >> 2); loop_cnt--;) {
        LOAD_4VECS_SB(src, src_stride, src0, src2, src4, src6);
        LOAD_4VECS_SB(src + 8, src_stride, src1, src3, src5, src7);
        src += (4 * src_stride);

        XORI_B_8VECS_SB(src0, src1, src2, src3, src4, src5, src6, src7,
                        src0, src1, src2, src3, src4, src5, src6, src7, 128);

        vec0 = __msa_vshf_b(mask0, src0, src0);
        vec1 = __msa_vshf_b(mask1, src0, src0);
        vec2 = __msa_vshf_b(mask2, src0, src0);
        vec3 = __msa_vshf_b(mask3, src0, src0);

        dst0 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask0, src1, src1);
        vec1 = __msa_vshf_b(mask1, src1, src1);
        vec2 = __msa_vshf_b(mask2, src1, src1);
        vec3 = __msa_vshf_b(mask3, src1, src1);

        dst1 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask0, src2, src2);
        vec1 = __msa_vshf_b(mask1, src2, src2);
        vec2 = __msa_vshf_b(mask2, src2, src2);
        vec3 = __msa_vshf_b(mask3, src2, src2);

        dst2 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask0, src3, src3);
        vec1 = __msa_vshf_b(mask1, src3, src3);
        vec2 = __msa_vshf_b(mask2, src3, src3);
        vec3 = __msa_vshf_b(mask3, src3, src3);

        dst3 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask0, src4, src4);
        vec1 = __msa_vshf_b(mask1, src4, src4);
        vec2 = __msa_vshf_b(mask2, src4, src4);
        vec3 = __msa_vshf_b(mask3, src4, src4);

        dst4 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask0, src5, src5);
        vec1 = __msa_vshf_b(mask1, src5, src5);
        vec2 = __msa_vshf_b(mask2, src5, src5);
        vec3 = __msa_vshf_b(mask3, src5, src5);

        dst5 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask0, src6, src6);
        vec1 = __msa_vshf_b(mask1, src6, src6);
        vec2 = __msa_vshf_b(mask2, src6, src6);
        vec3 = __msa_vshf_b(mask3, src6, src6);

        dst6 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask0, src7, src7);
        vec1 = __msa_vshf_b(mask1, src7, src7);
        vec2 = __msa_vshf_b(mask2, src7, src7);
        vec3 = __msa_vshf_b(mask3, src7, src7);

        dst7 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        STORE_SH(dst0, dst);
        STORE_SH(dst1, dst + 8);
        dst += dst_stride;
        STORE_SH(dst2, dst);
        STORE_SH(dst3, dst + 8);
        dst += dst_stride;
        STORE_SH(dst4, dst);
        STORE_SH(dst5, dst + 8);
        dst += dst_stride;
        STORE_SH(dst6, dst);
        STORE_SH(dst7, dst + 8);
        dst += dst_stride;
    }
}

static void hevc_hz_8t_24w_msa(uint8_t * __restrict src, int32_t src_stride,
                               int16_t * __restrict dst, int32_t dst_stride,
                               const int8_t * __restrict filter, int32_t height)
{
    uint32_t loop_cnt;
    v16i8 src0, src1, src2, src3;
    v8i16 filt0, filt1, filt2, filt3;
    v16i8 mask1, mask2, mask3, mask4, mask5, mask6, mask7;
    v16i8 vec0, vec1, vec2, vec3;
    v8i16 dst0, dst1, dst2, dst3, dst4, dst5;
    v8i16 filter_vec;
    v8u16 const_vec;
    v16i8 mask0 = { 0, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 7, 8 };

    src -= 3;

    filter_vec = LOAD_SH(filter);
    filt0 = __msa_splati_h(filter_vec, 0);
    filt1 = __msa_splati_h(filter_vec, 1);
    filt2 = __msa_splati_h(filter_vec, 2);
    filt3 = __msa_splati_h(filter_vec, 3);

    mask1 = mask0 + 2;
    mask2 = mask0 + 4;
    mask3 = mask0 + 6;
    mask4 = mask0 + 8;
    mask5 = mask0 + 10;
    mask6 = mask0 + 12;
    mask7 = mask0 + 14;

    const_vec = (v8u16) __msa_ldi_h(128);
    const_vec <<= 6;

    for (loop_cnt = (height >> 1); loop_cnt--;) {
        src0 = LOAD_SB(src);
        src1 = LOAD_SB(src + 16);
        src += src_stride;
        src2 = LOAD_SB(src);
        src3 = LOAD_SB(src + 16);
        src += src_stride;

        XORI_B_4VECS_SB(src0, src1, src2, src3, src0, src1, src2, src3, 128);

        vec0 = __msa_vshf_b(mask0, src0, src0);
        vec1 = __msa_vshf_b(mask1, src0, src0);
        vec2 = __msa_vshf_b(mask2, src0, src0);
        vec3 = __msa_vshf_b(mask3, src0, src0);

        dst0 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask4, src1, src0);
        vec1 = __msa_vshf_b(mask5, src1, src0);
        vec2 = __msa_vshf_b(mask6, src1, src0);
        vec3 = __msa_vshf_b(mask7, src1, src0);

        dst1 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask0, src1, src1);
        vec1 = __msa_vshf_b(mask1, src1, src1);
        vec2 = __msa_vshf_b(mask2, src1, src1);
        vec3 = __msa_vshf_b(mask3, src1, src1);

        dst2 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask0, src2, src2);
        vec1 = __msa_vshf_b(mask1, src2, src2);
        vec2 = __msa_vshf_b(mask2, src2, src2);
        vec3 = __msa_vshf_b(mask3, src2, src2);

        dst3 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask4, src3, src2);
        vec1 = __msa_vshf_b(mask5, src3, src2);
        vec2 = __msa_vshf_b(mask6, src3, src2);
        vec3 = __msa_vshf_b(mask7, src3, src2);

        dst4 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask0, src3, src3);
        vec1 = __msa_vshf_b(mask1, src3, src3);
        vec2 = __msa_vshf_b(mask2, src3, src3);
        vec3 = __msa_vshf_b(mask3, src3, src3);

        dst5 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        STORE_SH(dst0, dst);
        STORE_SH(dst1, dst + 8);
        STORE_SH(dst2, dst + 16);
        dst += dst_stride;
        STORE_SH(dst3, dst);
        STORE_SH(dst4, dst + 8);
        STORE_SH(dst5, dst + 16);
        dst += dst_stride;
    }
}

static void hevc_hz_8t_32w_msa(uint8_t * __restrict src, int32_t src_stride,
                               int16_t * __restrict dst, int32_t dst_stride,
                               const int8_t * __restrict filter, int32_t height)
{
    uint32_t loop_cnt;
    v16i8 src0, src1, src2;
    v8i16 filt0, filt1, filt2, filt3;
    v16i8 mask1, mask2, mask3, mask4, mask5, mask6, mask7;
    v16i8 vec0, vec1, vec2, vec3;
    v8i16 dst0, dst1, dst2, dst3;
    v8i16 filter_vec;
    v8u16 const_vec;
    v16i8 mask0 = { 0, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 7, 8 };

    src -= 3;

    filter_vec = LOAD_SH(filter);
    filt0 = __msa_splati_h(filter_vec, 0);
    filt1 = __msa_splati_h(filter_vec, 1);
    filt2 = __msa_splati_h(filter_vec, 2);
    filt3 = __msa_splati_h(filter_vec, 3);

    mask1 = mask0 + 2;
    mask2 = mask0 + 4;
    mask3 = mask0 + 6;
    mask4 = mask0 + 8;
    mask5 = mask0 + 10;
    mask6 = mask0 + 12;
    mask7 = mask0 + 14;

    const_vec = (v8u16) __msa_ldi_h(128);
    const_vec <<= 6;

    for (loop_cnt = height; loop_cnt--;) {
        src0 = LOAD_SB(src);
        src1 = LOAD_SB(src + 16);
        src2 = LOAD_SB(src + 24);
        src += src_stride;

        XORI_B_3VECS_SB(src0, src1, src2, src0, src1, src2, 128);

        vec0 = __msa_vshf_b(mask0, src0, src0);
        vec1 = __msa_vshf_b(mask1, src0, src0);
        vec2 = __msa_vshf_b(mask2, src0, src0);
        vec3 = __msa_vshf_b(mask3, src0, src0);

        dst0 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask4, src1, src0);
        vec1 = __msa_vshf_b(mask5, src1, src0);
        vec2 = __msa_vshf_b(mask6, src1, src0);
        vec3 = __msa_vshf_b(mask7, src1, src0);

        dst1 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask0, src1, src1);
        vec1 = __msa_vshf_b(mask1, src1, src1);
        vec2 = __msa_vshf_b(mask2, src1, src1);
        vec3 = __msa_vshf_b(mask3, src1, src1);

        dst2 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask0, src2, src2);
        vec1 = __msa_vshf_b(mask1, src2, src2);
        vec2 = __msa_vshf_b(mask2, src2, src2);
        vec3 = __msa_vshf_b(mask3, src2, src2);

        dst3 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        STORE_SH(dst0, dst);
        STORE_SH(dst1, dst + 8);
        STORE_SH(dst2, dst + 16);
        STORE_SH(dst3, dst + 24);
        dst += dst_stride;
    }
}

static void hevc_hz_8t_48w_msa(uint8_t * __restrict src, int32_t src_stride,
                               int16_t * __restrict dst, int32_t dst_stride,
                               const int8_t * __restrict filter, int32_t height)
{
    uint32_t loop_cnt;
    v16i8 src0, src1, src2, src3;
    v8i16 filt0, filt1, filt2, filt3;
    v16i8 mask1, mask2, mask3, mask4, mask5, mask6, mask7;
    v16i8 vec0, vec1, vec2, vec3;
    v8i16 dst0, dst1, dst2, dst3, dst4, dst5;
    v8i16 filter_vec;
    v8u16 const_vec;
    v16i8 mask0 = { 0, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 7, 8 };

    src -= 3;

    filter_vec = LOAD_SH(filter);
    filt0 = __msa_splati_h(filter_vec, 0);
    filt1 = __msa_splati_h(filter_vec, 1);
    filt2 = __msa_splati_h(filter_vec, 2);
    filt3 = __msa_splati_h(filter_vec, 3);

    mask1 = mask0 + 2;
    mask2 = mask0 + 4;
    mask3 = mask0 + 6;
    mask4 = mask0 + 8;
    mask5 = mask0 + 10;
    mask6 = mask0 + 12;
    mask7 = mask0 + 14;

    const_vec = (v8u16) __msa_ldi_h(128);
    const_vec <<= 6;

    for (loop_cnt = height; loop_cnt--;) {
        src0 = LOAD_SB(src);
        src1 = LOAD_SB(src + 16);
        src2 = LOAD_SB(src + 32);
        src3 = LOAD_SB(src + 40);
        src += src_stride;

        XORI_B_4VECS_SB(src0, src1, src2, src3, src0, src1, src2, src3, 128);

        vec0 = __msa_vshf_b(mask0, src0, src0);
        vec1 = __msa_vshf_b(mask1, src0, src0);
        vec2 = __msa_vshf_b(mask2, src0, src0);
        vec3 = __msa_vshf_b(mask3, src0, src0);

        dst0 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask4, src1, src0);
        vec1 = __msa_vshf_b(mask5, src1, src0);
        vec2 = __msa_vshf_b(mask6, src1, src0);
        vec3 = __msa_vshf_b(mask7, src1, src0);

        dst1 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask0, src1, src1);
        vec1 = __msa_vshf_b(mask1, src1, src1);
        vec2 = __msa_vshf_b(mask2, src1, src1);
        vec3 = __msa_vshf_b(mask3, src1, src1);

        dst2 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask4, src2, src1);
        vec1 = __msa_vshf_b(mask5, src2, src1);
        vec2 = __msa_vshf_b(mask6, src2, src1);
        vec3 = __msa_vshf_b(mask7, src2, src1);

        dst3 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask0, src2, src2);
        vec1 = __msa_vshf_b(mask1, src2, src2);
        vec2 = __msa_vshf_b(mask2, src2, src2);
        vec3 = __msa_vshf_b(mask3, src2, src2);

        dst4 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        vec0 = __msa_vshf_b(mask0, src3, src3);
        vec1 = __msa_vshf_b(mask1, src3, src3);
        vec2 = __msa_vshf_b(mask2, src3, src3);
        vec3 = __msa_vshf_b(mask3, src3, src3);

        dst5 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        STORE_SH(dst0, dst);
        STORE_SH(dst1, dst + 8);
        STORE_SH(dst2, dst + 16);
        STORE_SH(dst3, dst + 24);
        STORE_SH(dst4, dst + 32);
        STORE_SH(dst5, dst + 40);
        dst += dst_stride;
    }
}

static void hevc_hz_8t_64w_msa(uint8_t * __restrict src, int32_t src_stride,
                               int16_t * __restrict dst, int32_t dst_stride,
                               const int8_t * __restrict filter, int32_t height)
{
    uint32_t loop_cnt;
    v16i8 src0, src1, src2, src3, src4;
    v8i16 filt0, filt1, filt2, filt3;
    v16i8 mask1, mask2, mask3, mask4, mask5, mask6, mask7;
    v16i8 vec0, vec1, vec2, vec3;
    v8i16 dst0, dst1, dst2, dst3, dst4, dst5, dst6, dst7;
    v8i16 filter_vec;
    v8u16 const_vec;
    v16i8 mask0 = { 0, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 7, 8 };

    src -= 3;

    filter_vec = LOAD_SH(filter);
    filt0 = __msa_splati_h(filter_vec, 0);
    filt1 = __msa_splati_h(filter_vec, 1);
    filt2 = __msa_splati_h(filter_vec, 2);
    filt3 = __msa_splati_h(filter_vec, 3);

    mask1 = mask0 + 2;
    mask2 = mask0 + 4;
    mask3 = mask0 + 6;
    mask4 = mask0 + 8;
    mask5 = mask0 + 10;
    mask6 = mask0 + 12;
    mask7 = mask0 + 14;

    const_vec = (v8u16) __msa_ldi_h(128);
    const_vec <<= 6;

    for (loop_cnt = height; loop_cnt--;) {
        src0 = LOAD_SB(src);
        src1 = LOAD_SB(src + 16);
        src2 = LOAD_SB(src + 32);
        src3 = LOAD_SB(src + 48);
        src4 = LOAD_SB(src + 56);
        src += src_stride;

        XORI_B_5VECS_SB(src0, src1, src2, src3, src4,
                        src0, src1, src2, src3, src4, 128);

        vec0 = __msa_vshf_b(mask0, src0, src0);
        vec1 = __msa_vshf_b(mask1, src0, src0);
        vec2 = __msa_vshf_b(mask2, src0, src0);
        vec3 = __msa_vshf_b(mask3, src0, src0);

        dst0 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        STORE_SH(dst0, dst);

        vec0 = __msa_vshf_b(mask4, src1, src0);
        vec1 = __msa_vshf_b(mask5, src1, src0);
        vec2 = __msa_vshf_b(mask6, src1, src0);
        vec3 = __msa_vshf_b(mask7, src1, src0);

        dst1 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        STORE_SH(dst1, dst + 8);

        vec0 = __msa_vshf_b(mask0, src1, src1);
        vec1 = __msa_vshf_b(mask1, src1, src1);
        vec2 = __msa_vshf_b(mask2, src1, src1);
        vec3 = __msa_vshf_b(mask3, src1, src1);

        dst2 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        STORE_SH(dst2, dst + 16);

        vec0 = __msa_vshf_b(mask4, src2, src1);
        vec1 = __msa_vshf_b(mask5, src2, src1);
        vec2 = __msa_vshf_b(mask6, src2, src1);
        vec3 = __msa_vshf_b(mask7, src2, src1);

        dst3 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        STORE_SH(dst3, dst + 24);

        vec0 = __msa_vshf_b(mask0, src2, src2);
        vec1 = __msa_vshf_b(mask1, src2, src2);
        vec2 = __msa_vshf_b(mask2, src2, src2);
        vec3 = __msa_vshf_b(mask3, src2, src2);

        dst4 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        STORE_SH(dst4, dst + 32);

        vec0 = __msa_vshf_b(mask4, src3, src2);
        vec1 = __msa_vshf_b(mask5, src3, src2);
        vec2 = __msa_vshf_b(mask6, src3, src2);
        vec3 = __msa_vshf_b(mask7, src3, src2);

        dst5 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        STORE_SH(dst5, dst + 40);

        vec0 = __msa_vshf_b(mask0, src3, src3);
        vec1 = __msa_vshf_b(mask1, src3, src3);
        vec2 = __msa_vshf_b(mask2, src3, src3);
        vec3 = __msa_vshf_b(mask3, src3, src3);

        dst6 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        STORE_SH(dst6, dst + 48);

        vec0 = __msa_vshf_b(mask0, src4, src4);
        vec1 = __msa_vshf_b(mask1, src4, src4);
        vec2 = __msa_vshf_b(mask2, src4, src4);
        vec3 = __msa_vshf_b(mask3, src4, src4);

        dst7 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        STORE_SH(dst7, dst + 56);

        dst += dst_stride;
    }
}

static void hevc_vt_8t_4w_msa(uint8_t * __restrict src, int32_t src_stride,
                              int16_t * __restrict dst, int32_t dst_stride,
                              const int8_t * __restrict filter, int32_t height)
{
    int32_t loop_cnt;
    uint64_t out0, out1, out2, out3;
    v16i8 src0, src1, src2, src3, src4, src5, src6, src7, src8;
    v16i8 src9, src10, src11, src12, src13, src14;
    v16i8 src10_r, src32_r, src54_r, src76_r, src98_r;
    v16i8 src21_r, src43_r, src65_r, src87_r, src109_r;
    v16i8 src1110_r, src1211_r, src1312_r, src1413_r;
    v16i8 src2110, src4332, src6554, src8776, src10998;
    v16i8 src12111110, src14131312;
    v8i16 dst10, dst32, dst54, dst76;
    v8i16 filter_vec;
    v8i16 filt0, filt1, filt2, filt3;
    v8u16 const_vec;

    src -= (3 * src_stride);

    const_vec = (v8u16) __msa_ldi_h(128);
    const_vec <<= 6;

    filter_vec = LOAD_SH(filter);
    filt0 = __msa_splati_h(filter_vec, 0);
    filt1 = __msa_splati_h(filter_vec, 1);
    filt2 = __msa_splati_h(filter_vec, 2);
    filt3 = __msa_splati_h(filter_vec, 3);

    LOAD_7VECS_SB(src, src_stride,
                  src0, src1, src2, src3, src4, src5, src6);
    src += (7 * src_stride);

    ILVR_B_6VECS_SB(src0, src2, src4, src1, src3, src5,
                    src1, src3, src5, src2, src4, src6,
                    src10_r, src32_r, src54_r, src21_r, src43_r, src65_r);

    ILVR_D_3VECS_SB(src2110, src21_r, src10_r, src4332, src43_r, src32_r,
                    src6554, src65_r, src54_r);

    XORI_B_3VECS_SB(src2110, src4332, src6554, src2110, src4332, src6554, 128);

    for (loop_cnt = (height >> 3); loop_cnt--;) {
        LOAD_8VECS_SB(src, src_stride,
                      src7, src8, src9, src10, src11, src12, src13, src14);
        src += (8 * src_stride);

        ILVR_B_8VECS_SB(src6, src7, src8, src9, src10, src11, src12, src13,
                        src7, src8, src9, src10, src11, src12, src13, src14,
                        src76_r, src87_r, src98_r, src109_r,
                        src1110_r, src1211_r, src1312_r, src1413_r);

        ILVR_D_4VECS_SB(src8776, src87_r, src76_r, src10998, src109_r, src98_r,
                        src12111110, src1211_r, src1110_r,
                        src14131312, src1413_r, src1312_r);

        XORI_B_4VECS_SB(src8776, src10998, src12111110, src14131312,
                        src8776, src10998, src12111110, src14131312, 128);

        dst10 = HEVC_FILT_8TAP_DPADD_H(src2110, src4332, src6554, src8776,
                                       filt0, filt1, filt2, filt3, const_vec);

        dst32 = HEVC_FILT_8TAP_DPADD_H(src4332, src6554, src8776, src10998,
                                       filt0, filt1, filt2, filt3, const_vec);

        dst54 = HEVC_FILT_8TAP_DPADD_H(src6554, src8776, src10998, src12111110,
                                       filt0, filt1, filt2, filt3, const_vec);

        dst76 = HEVC_FILT_8TAP_DPADD_H(src8776, src10998,
                                       src12111110, src14131312,
                                       filt0, filt1, filt2, filt3, const_vec);

        out0 = __msa_copy_u_d((v2i64) dst10, 0);
        out1 = __msa_copy_u_d((v2i64) dst10, 1);
        out2 = __msa_copy_u_d((v2i64) dst32, 0);
        out3 = __msa_copy_u_d((v2i64) dst32, 1);

        STORE_DWORD(dst, out0);
        dst += dst_stride;
        STORE_DWORD(dst, out1);
        dst += dst_stride;
        STORE_DWORD(dst, out2);
        dst += dst_stride;
        STORE_DWORD(dst, out3);
        dst += dst_stride;

        out0 = __msa_copy_u_d((v2i64) dst54, 0);
        out1 = __msa_copy_u_d((v2i64) dst54, 1);
        out2 = __msa_copy_u_d((v2i64) dst76, 0);
        out3 = __msa_copy_u_d((v2i64) dst76, 1);

        STORE_DWORD(dst, out0);
        dst += dst_stride;
        STORE_DWORD(dst, out1);
        dst += dst_stride;
        STORE_DWORD(dst, out2);
        dst += dst_stride;
        STORE_DWORD(dst, out3);
        dst += dst_stride;

        src2110 = src10998;
        src4332 = src12111110;
        src6554 = src14131312;

        src6 = src14;
    }
}

static void hevc_vt_8t_8w_msa(uint8_t * __restrict src, int32_t src_stride,
                              int16_t * __restrict dst, int32_t dst_stride,
                              const int8_t * __restrict filter, int32_t height)
{
    int32_t loop_cnt;
    v16i8 src0, src1, src2, src3, src4, src5, src6, src7, src8, src9, src10;
    v16i8 src10_r, src32_r, src54_r, src76_r, src98_r;
    v16i8 src21_r, src43_r, src65_r, src87_r, src109_r;
    v8i16 dst0_r, dst1_r, dst2_r, dst3_r;
    v8i16 filter_vec;
    v8i16 filt0, filt1, filt2, filt3;
    v8u16 const_vec;

    src -= (3 * src_stride);

    const_vec = (v8u16) __msa_ldi_h(128);
    const_vec <<= 6;

    filter_vec = LOAD_SH(filter);
    filt0 = __msa_splati_h(filter_vec, 0);
    filt1 = __msa_splati_h(filter_vec, 1);
    filt2 = __msa_splati_h(filter_vec, 2);
    filt3 = __msa_splati_h(filter_vec, 3);

    LOAD_7VECS_SB(src, src_stride,
                  src0, src1, src2, src3, src4, src5, src6);
    src += (7 * src_stride);

    XORI_B_7VECS_SB(src0, src1, src2, src3, src4, src5, src6,
                    src0, src1, src2, src3, src4, src5, src6, 128);

    ILVR_B_6VECS_SB(src0, src2, src4, src1, src3, src5,
                    src1, src3, src5, src2, src4, src6,
                    src10_r, src32_r, src54_r, src21_r, src43_r, src65_r);

    for (loop_cnt = (height >> 2); loop_cnt--;) {
        LOAD_4VECS_SB(src, src_stride, src7, src8, src9, src10);
        src += (4 * src_stride);

        XORI_B_4VECS_SB(src7, src8, src9, src10, src7, src8, src9, src10, 128);

        ILVR_B_4VECS_SB(src6, src7, src8, src9, src7, src8, src9, src10,
                        src76_r, src87_r, src98_r, src109_r);

        dst0_r = HEVC_FILT_8TAP_DPADD_H(src10_r, src32_r, src54_r, src76_r,
                                        filt0, filt1, filt2, filt3, const_vec);

        dst1_r = HEVC_FILT_8TAP_DPADD_H(src21_r, src43_r, src65_r, src87_r,
                                        filt0, filt1, filt2, filt3, const_vec);

        dst2_r = HEVC_FILT_8TAP_DPADD_H(src32_r, src54_r, src76_r, src98_r,
                                        filt0, filt1, filt2, filt3, const_vec);

        dst3_r = HEVC_FILT_8TAP_DPADD_H(src43_r, src65_r, src87_r, src109_r,
                                        filt0, filt1, filt2, filt3, const_vec);

        STORE_SH(dst0_r, dst);
        dst += dst_stride;
        STORE_SH(dst1_r, dst);
        dst += dst_stride;
        STORE_SH(dst2_r, dst);
        dst += dst_stride;
        STORE_SH(dst3_r, dst);
        dst += dst_stride;

        src10_r = src54_r;
        src32_r = src76_r;
        src54_r = src98_r;

        src21_r = src65_r;
        src43_r = src87_r;
        src65_r = src109_r;

        src6 = src10;
    }
}

static void hevc_vt_8t_12w_msa(uint8_t * __restrict src, int32_t src_stride,
                               int16_t * __restrict dst, int32_t dst_stride,
                               const int8_t * __restrict filter, int32_t height)
{
    int32_t loop_cnt;
    uint64_t out0, out1, out2, out3;
    v16i8 src0, src1, src2, src3, src4, src5, src6, src7, src8, src9, src10;
    v16i8 src10_r, src32_r, src54_r, src76_r, src98_r;
    v16i8 src21_r, src43_r, src65_r, src87_r, src109_r;
    v8i16 dst0_r, dst1_r, dst2_r, dst3_r;
    v16i8 src10_l, src32_l, src54_l, src76_l, src98_l;
    v16i8 src21_l, src43_l, src65_l, src87_l, src109_l;
    v16i8 src2110, src4332, src6554, src8776, src10998;
    v8i16 dst0_l, dst1_l;
    v8i16 filter_vec;
    v8i16 filt0, filt1, filt2, filt3;
    v8u16 const_vec;

    src -= (3 * src_stride);

    const_vec = (v8u16) __msa_ldi_h(128);
    const_vec <<= 6;

    filter_vec = LOAD_SH(filter);
    filt0 = __msa_splati_h(filter_vec, 0);
    filt1 = __msa_splati_h(filter_vec, 1);
    filt2 = __msa_splati_h(filter_vec, 2);
    filt3 = __msa_splati_h(filter_vec, 3);

    LOAD_7VECS_SB(src, src_stride,
                  src0, src1, src2, src3, src4, src5, src6);
    src += (7 * src_stride);

    XORI_B_7VECS_SB(src0, src1, src2, src3, src4, src5, src6,
                    src0, src1, src2, src3, src4, src5, src6, 128);

    ILVR_B_6VECS_SB(src0, src2, src4, src1, src3, src5,
                    src1, src3, src5, src2, src4, src6,
                    src10_r, src32_r, src54_r, src21_r, src43_r, src65_r);

    ILVL_B_6VECS_SB(src0, src2, src4, src1, src3, src5,
                    src1, src3, src5, src2, src4, src6,
                    src10_l, src32_l, src54_l, src21_l, src43_l, src65_l);

    ILVR_D_3VECS_SB(src2110, src21_l, src10_l, src4332, src43_l, src32_l,
                    src6554, src65_l, src54_l);

    for (loop_cnt = (height >> 2); loop_cnt--;) {
        LOAD_4VECS_SB(src, src_stride, src7, src8, src9, src10);
        src += (4 * src_stride);

        XORI_B_4VECS_SB(src7, src8, src9, src10, src7, src8, src9, src10, 128);

        ILVR_B_4VECS_SB(src6, src7, src8, src9, src7, src8, src9, src10,
                        src76_r, src87_r, src98_r, src109_r);

        ILVL_B_4VECS_SB(src6, src7, src8, src9, src7, src8, src9, src10,
                        src76_l, src87_l, src98_l, src109_l);

        ILVR_D_2VECS_SB(src8776, src87_l, src76_l, src10998, src109_l, src98_l);

        dst0_r = HEVC_FILT_8TAP_DPADD_H(src10_r, src32_r, src54_r, src76_r,
                                        filt0, filt1, filt2, filt3, const_vec);

        dst1_r = HEVC_FILT_8TAP_DPADD_H(src21_r, src43_r, src65_r, src87_r,
                                        filt0, filt1, filt2, filt3, const_vec);

        dst2_r = HEVC_FILT_8TAP_DPADD_H(src32_r, src54_r, src76_r, src98_r,
                                        filt0, filt1, filt2, filt3, const_vec);

        dst3_r = HEVC_FILT_8TAP_DPADD_H(src43_r, src65_r, src87_r, src109_r,
                                        filt0, filt1, filt2, filt3, const_vec);

        dst0_l = HEVC_FILT_8TAP_DPADD_H(src2110, src4332, src6554, src8776,
                                        filt0, filt1, filt2, filt3, const_vec);

        dst1_l = HEVC_FILT_8TAP_DPADD_H(src4332, src6554, src8776, src10998,
                                        filt0, filt1, filt2, filt3, const_vec);

        out0 = __msa_copy_u_d((v2i64) dst0_l, 0);
        out1 = __msa_copy_u_d((v2i64) dst0_l, 1);
        out2 = __msa_copy_u_d((v2i64) dst1_l, 0);
        out3 = __msa_copy_u_d((v2i64) dst1_l, 1);

        STORE_SH(dst0_r, dst);
        STORE_DWORD(dst + 8, out0);
        dst += dst_stride;
        STORE_SH(dst1_r, dst);
        STORE_DWORD(dst + 8, out1);
        dst += dst_stride;

        STORE_SH(dst2_r, dst);
        STORE_DWORD(dst + 8, out2);
        dst += dst_stride;
        STORE_SH(dst3_r, dst);
        STORE_DWORD(dst + 8, out3);
        dst += dst_stride;

        src10_r = src54_r;
        src32_r = src76_r;
        src54_r = src98_r;

        src21_r = src65_r;
        src43_r = src87_r;
        src65_r = src109_r;

        src2110 = src6554;
        src4332 = src8776;
        src6554 = src10998;

        src6 = src10;
    }
}

static void hevc_vt_8t_16multx4mult_msa(uint8_t * __restrict src,
                                        int32_t src_stride,
                                        int16_t * __restrict dst,
                                        int32_t dst_stride,
                                        const int8_t * __restrict filter,
                                        int32_t height,
                                        int32_t width)
{
    uint8_t *src_tmp;
    int16_t *dst_tmp;
    int32_t loop_cnt, cnt;
    v16i8 src0, src1, src2, src3, src4, src5, src6, src7, src8, src9, src10;
    v16i8 src10_r, src32_r, src54_r, src76_r, src98_r;
    v16i8 src21_r, src43_r, src65_r, src87_r, src109_r;
    v8i16 dst0_r, dst1_r, dst2_r, dst3_r;
    v16i8 src10_l, src32_l, src54_l, src76_l, src98_l;
    v16i8 src21_l, src43_l, src65_l, src87_l, src109_l;
    v8i16 dst0_l, dst1_l, dst2_l, dst3_l;
    v8i16 filter_vec;
    v8i16 filt0, filt1, filt2, filt3;
    v8u16 const_vec;

    src -= (3 * src_stride);

    const_vec = (v8u16) __msa_ldi_h(128);
    const_vec <<= 6;

    filter_vec = LOAD_SH(filter);
    filt0 = __msa_splati_h(filter_vec, 0);
    filt1 = __msa_splati_h(filter_vec, 1);
    filt2 = __msa_splati_h(filter_vec, 2);
    filt3 = __msa_splati_h(filter_vec, 3);

    for (cnt = width >> 4; cnt--;) {
        src_tmp = src;
        dst_tmp = dst;

        LOAD_7VECS_SB(src_tmp, src_stride,
                      src0, src1, src2, src3, src4, src5, src6);
        src_tmp += (7 * src_stride);

        XORI_B_7VECS_SB(src0, src1, src2, src3, src4, src5, src6,
                        src0, src1, src2, src3, src4, src5, src6, 128);

        ILVR_B_6VECS_SB(src0, src2, src4, src1, src3, src5,
                        src1, src3, src5, src2, src4, src6,
                        src10_r, src32_r, src54_r, src21_r, src43_r, src65_r);

        ILVL_B_6VECS_SB(src0, src2, src4, src1, src3, src5,
                        src1, src3, src5, src2, src4, src6,
                        src10_l, src32_l, src54_l, src21_l, src43_l, src65_l);

        for (loop_cnt = (height >> 2); loop_cnt--;) {
            LOAD_4VECS_SB(src_tmp, src_stride, src7, src8, src9, src10);
            src_tmp += (4 * src_stride);

            XORI_B_4VECS_SB(src7, src8, src9, src10,
                            src7, src8, src9, src10, 128);

            ILVR_B_4VECS_SB(src6, src7, src8, src9, src7, src8, src9, src10,
                            src76_r, src87_r, src98_r, src109_r);

            ILVL_B_4VECS_SB(src6, src7, src8, src9, src7, src8, src9, src10,
                            src76_l, src87_l, src98_l, src109_l);

            dst0_r = HEVC_FILT_8TAP_DPADD_H(src10_r, src32_r, src54_r, src76_r,
                                            filt0, filt1, filt2, filt3,
                                            const_vec);

            dst1_r = HEVC_FILT_8TAP_DPADD_H(src21_r, src43_r, src65_r, src87_r,
                                            filt0, filt1, filt2, filt3,
                                            const_vec);

            dst2_r = HEVC_FILT_8TAP_DPADD_H(src32_r, src54_r, src76_r, src98_r,
                                            filt0, filt1, filt2, filt3,
                                            const_vec);

            dst3_r = HEVC_FILT_8TAP_DPADD_H(src43_r, src65_r, src87_r, src109_r,
                                            filt0, filt1, filt2, filt3,
                                            const_vec);

            dst0_l = HEVC_FILT_8TAP_DPADD_H(src10_l, src32_l, src54_l, src76_l,
                                            filt0, filt1, filt2, filt3,
                                            const_vec);

            dst1_l = HEVC_FILT_8TAP_DPADD_H(src21_l, src43_l, src65_l, src87_l,
                                            filt0, filt1, filt2, filt3,
                                            const_vec);

            dst2_l = HEVC_FILT_8TAP_DPADD_H(src32_l, src54_l, src76_l, src98_l,
                                            filt0, filt1, filt2, filt3,
                                            const_vec);

            dst3_l = HEVC_FILT_8TAP_DPADD_H(src43_l, src65_l, src87_l, src109_l,
                                            filt0, filt1, filt2, filt3,
                                            const_vec);

            STORE_SH(dst0_r, dst_tmp);
            STORE_SH(dst0_l, dst_tmp + 8);
            dst_tmp += dst_stride;
            STORE_SH(dst1_r, dst_tmp);
            STORE_SH(dst1_l, dst_tmp + 8);
            dst_tmp += dst_stride;

            STORE_SH(dst2_r, dst_tmp);
            STORE_SH(dst2_l, dst_tmp + 8);
            dst_tmp += dst_stride;
            STORE_SH(dst3_r, dst_tmp);
            STORE_SH(dst3_l, dst_tmp + 8);
            dst_tmp += dst_stride;

            src10_r = src54_r;
            src32_r = src76_r;
            src54_r = src98_r;

            src21_r = src65_r;
            src43_r = src87_r;
            src65_r = src109_r;

            src10_l = src54_l;
            src32_l = src76_l;
            src54_l = src98_l;

            src21_l = src65_l;
            src43_l = src87_l;
            src65_l = src109_l;

            src6 = src10;
        }

        src += 16;
        dst += 16;
    }
}

static void hevc_vt_8t_16w_msa(uint8_t * __restrict src, int32_t src_stride,
                               int16_t * __restrict dst, int32_t dst_stride,
                               const int8_t * __restrict filter, int32_t height)
{
    hevc_vt_8t_16multx4mult_msa(src, src_stride, dst, dst_stride,
                                filter, height, 16);
}

static void hevc_vt_8t_24w_msa(uint8_t * __restrict src, int32_t src_stride,
                               int16_t * __restrict dst, int32_t dst_stride,
                               const int8_t * __restrict filter, int32_t height)
{
    hevc_vt_8t_16multx4mult_msa(src, src_stride, dst, dst_stride,
                                filter, height, 16);

    hevc_vt_8t_8w_msa(src + 16, src_stride, dst + 16, dst_stride,
                      filter, height);
}

static void hevc_vt_8t_32w_msa(uint8_t * __restrict src, int32_t src_stride,
                               int16_t * __restrict dst, int32_t dst_stride,
                               const int8_t * __restrict filter, int32_t height)
{
    hevc_vt_8t_16multx4mult_msa(src, src_stride, dst, dst_stride,
                                filter, height, 32);
}

static void hevc_vt_8t_48w_msa(uint8_t * __restrict src, int32_t src_stride,
                               int16_t * __restrict dst, int32_t dst_stride,
                               const int8_t * __restrict filter, int32_t height)
{
    hevc_vt_8t_16multx4mult_msa(src, src_stride, dst, dst_stride,
                                filter, height, 48);
}

static void hevc_vt_8t_64w_msa(uint8_t * __restrict src, int32_t src_stride,
                               int16_t * __restrict dst, int32_t dst_stride,
                               const int8_t * __restrict filter, int32_t height)
{
    hevc_vt_8t_16multx4mult_msa(src, src_stride, dst, dst_stride,
                                filter, height, 64);
}

static void hevc_hv_8t_4w_msa(uint8_t * __restrict src, int32_t src_stride,
                              int16_t * __restrict dst, int32_t dst_stride,
                              const int8_t * __restrict filter_x,
                              const int8_t * __restrict filter_y,
                              int32_t height)
{
    uint32_t loop_cnt;
    uint64_t out0, out1;
    v16i8 src0, src1, src2, src3, src4, src5, src6, src7, src8;
    v8i16 filt0, filt1, filt2, filt3, filter_vec;
    v4i32 filt_h0, filt_h1, filt_h2, filt_h3;
    v16i8 mask1, mask2, mask3;
    v8u16 const_vec;
    v16i8 vec0, vec1, vec2, vec3, vec4, vec5, vec6, vec7;
    v16i8 vec8, vec9, vec10, vec11, vec12, vec13, vec14, vec15;
    v8i16 dst30, dst41, dst52, dst63, dst66, dst87;
    v4i32 dst0_r, dst1_r;
    v8i16 dst10_r, dst32_r, dst54_r, dst76_r;
    v8i16 dst21_r, dst43_r, dst65_r, dst87_r;
    v16i8 tmp;
    v16i8 mask0 = {
        0, 1, 1, 2, 2, 3, 3, 4, 16, 17, 17, 18, 18, 19, 19, 20
    };
    v8u16 mask4 = { 0, 4, 1, 5, 2, 6, 3, 7 };

    src -= ((3 * src_stride) + 3);

    filter_vec = LOAD_SH(filter_x);
    filt0 = __msa_splati_h(filter_vec, 0);
    filt1 = __msa_splati_h(filter_vec, 1);
    filt2 = __msa_splati_h(filter_vec, 2);
    filt3 = __msa_splati_h(filter_vec, 3);

    filter_vec = LOAD_SH(filter_y);
    tmp = __msa_clti_s_b((v16i8) filter_vec, 0);
    filter_vec = (v8i16) __msa_ilvr_b(tmp, (v16i8) filter_vec);

    filt_h0 = __msa_splati_w((v4i32) filter_vec, 0);
    filt_h1 = __msa_splati_w((v4i32) filter_vec, 1);
    filt_h2 = __msa_splati_w((v4i32) filter_vec, 2);
    filt_h3 = __msa_splati_w((v4i32) filter_vec, 3);

    mask1 = mask0 + 2;
    mask2 = mask0 + 4;
    mask3 = mask0 + 6;

    const_vec = (v8u16) __msa_ldi_h(128);
    const_vec <<= 6;

    LOAD_7VECS_SB(src, src_stride,
                  src0, src1, src2, src3, src4, src5, src6);
    src += (7 * src_stride);

    XORI_B_7VECS_SB(src0, src1, src2, src3, src4, src5, src6,
                    src0, src1, src2, src3, src4, src5, src6, 128);

    /* Row 0 Row 1 Row 2 Row 3 */
    vec0 = __msa_vshf_b(mask0, src3, src0);
    vec1 = __msa_vshf_b(mask1, src3, src0);
    vec2 = __msa_vshf_b(mask2, src3, src0);
    vec3 = __msa_vshf_b(mask3, src3, src0);

    vec4 = __msa_vshf_b(mask0, src4, src1);
    vec5 = __msa_vshf_b(mask1, src4, src1);
    vec6 = __msa_vshf_b(mask2, src4, src1);
    vec7 = __msa_vshf_b(mask3, src4, src1);

    vec8 = __msa_vshf_b(mask0, src5, src2);
    vec9 = __msa_vshf_b(mask1, src5, src2);
    vec10 = __msa_vshf_b(mask2, src5, src2);
    vec11 = __msa_vshf_b(mask3, src5, src2);

    vec12 = __msa_vshf_b(mask0, src6, src3);
    vec13 = __msa_vshf_b(mask1, src6, src3);
    vec14 = __msa_vshf_b(mask2, src6, src3);
    vec15 = __msa_vshf_b(mask3, src6, src3);

    dst30 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                   filt0, filt1, filt2, filt3, const_vec);

    dst41 = HEVC_FILT_8TAP_DPADD_H(vec4, vec5, vec6, vec7,
                                   filt0, filt1, filt2, filt3, const_vec);

    dst52 = HEVC_FILT_8TAP_DPADD_H(vec8, vec9, vec10, vec11,
                                   filt0, filt1, filt2, filt3, const_vec);

    dst63 = HEVC_FILT_8TAP_DPADD_H(vec12, vec13, vec14, vec15,
                                   filt0, filt1, filt2, filt3, const_vec);

    dst10_r = __msa_ilvr_h(dst41, dst30);
    dst21_r = __msa_ilvr_h(dst52, dst41);
    dst32_r = __msa_ilvr_h(dst63, dst52);

    dst43_r = __msa_ilvl_h(dst41, dst30);
    dst54_r = __msa_ilvl_h(dst52, dst41);
    dst65_r = __msa_ilvl_h(dst63, dst52);

    dst66 = (v8i16) __msa_splati_d((v2i64) dst63, 1);

    for (loop_cnt = height >> 1; loop_cnt--;) {
        LOAD_2VECS_SB(src, src_stride, src7, src8);
        src += (2 * src_stride);

        src7 = (v16i8) __msa_xori_b((v16u8) src7, 128);
        src8 = (v16i8) __msa_xori_b((v16u8) src8, 128);

        vec0 = __msa_vshf_b(mask0, src8, src7);
        vec1 = __msa_vshf_b(mask1, src8, src7);
        vec2 = __msa_vshf_b(mask2, src8, src7);
        vec3 = __msa_vshf_b(mask3, src8, src7);

        dst87 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                       filt0, filt1, filt2, filt3, const_vec);

        dst76_r = __msa_ilvr_h(dst87, dst66);

        dst0_r = HEVC_FILT_8TAP_DPADD_W(dst10_r, dst32_r, dst54_r, dst76_r,
                                        filt_h0, filt_h1, filt_h2, filt_h3);

        dst87_r = __msa_vshf_h((v8i16) mask4, dst87, dst87);

        dst1_r = HEVC_FILT_8TAP_DPADD_W(dst21_r, dst43_r, dst65_r, dst87_r,
                                        filt_h0, filt_h1, filt_h2, filt_h3);

        dst0_r >>= 6;
        dst1_r >>= 6;

        dst0_r = (v4i32) __msa_pckev_h((v8i16) dst1_r, (v8i16) dst0_r);

        out0 = __msa_copy_u_d((v2i64) dst0_r, 0);
        out1 = __msa_copy_u_d((v2i64) dst0_r, 1);

        STORE_DWORD(dst, out0);
        dst += dst_stride;
        STORE_DWORD(dst, out1);
        dst += dst_stride;

        dst10_r = dst32_r;
        dst32_r = dst54_r;
        dst54_r = dst76_r;

        dst21_r = dst43_r;
        dst43_r = dst65_r;
        dst65_r = dst87_r;

        dst66 = (v8i16) __msa_splati_d((v2i64) dst87, 1);
    }
}

static void hevc_hv_8t_8multx2mult_msa(uint8_t * __restrict src,
                                       int32_t src_stride,
                                       int16_t * __restrict dst,
                                       int32_t dst_stride,
                                       const int8_t * __restrict filter_x,
                                       const int8_t * __restrict filter_y,
                                       int32_t height, int32_t width)
{
    uint32_t loop_cnt, cnt;
    uint8_t *src_tmp;
    int16_t *dst_tmp;
    v16i8 src0, src1, src2, src3, src4, src5, src6, src7, src8;
    v8i16 filt0, filt1, filt2, filt3, filter_vec;
    v4i32 filt_h0, filt_h1, filt_h2, filt_h3;
    v16i8 mask1, mask2, mask3;
    v8u16 const_vec;
    v16i8 vec0, vec1, vec2, vec3, vec4, vec5, vec6, vec7;
    v16i8 vec8, vec9, vec10, vec11, vec12, vec13, vec14, vec15;
    v8i16 dst0, dst1, dst2, dst3, dst4, dst5, dst6, dst7, dst8;
    v4i32 dst0_r, dst0_l;
    v8i16 dst10_r, dst32_r, dst54_r, dst76_r;
    v8i16 dst10_l, dst32_l, dst54_l, dst76_l;
    v8i16 dst21_r, dst43_r, dst65_r, dst87_r;
    v8i16 dst21_l, dst43_l, dst65_l, dst87_l;
    v16i8 tmp;
    v16i8 mask0 = { 0, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 7, 8 };

    src -= ((3 * src_stride) + 3);

    filter_vec = LOAD_SH(filter_x);
    filt0 = __msa_splati_h(filter_vec, 0);
    filt1 = __msa_splati_h(filter_vec, 1);
    filt2 = __msa_splati_h(filter_vec, 2);
    filt3 = __msa_splati_h(filter_vec, 3);

    filter_vec = LOAD_SH(filter_y);
    tmp = __msa_clti_s_b((v16i8) filter_vec, 0);
    filter_vec = (v8i16) __msa_ilvr_b(tmp, (v16i8) filter_vec);

    filt_h0 = __msa_splati_w((v4i32) filter_vec, 0);
    filt_h1 = __msa_splati_w((v4i32) filter_vec, 1);
    filt_h2 = __msa_splati_w((v4i32) filter_vec, 2);
    filt_h3 = __msa_splati_w((v4i32) filter_vec, 3);

    mask1 = mask0 + 2;
    mask2 = mask0 + 4;
    mask3 = mask0 + 6;

    const_vec = (v8u16) __msa_ldi_h(128);
    const_vec <<= 6;

    for (cnt = width >> 3; cnt--;) {
        src_tmp = src;
        dst_tmp = dst;

        LOAD_7VECS_SB(src_tmp, src_stride,
                      src0, src1, src2, src3, src4, src5, src6);
        src_tmp += (7 * src_stride);

        XORI_B_7VECS_SB(src0, src1, src2, src3, src4, src5, src6,
                        src0, src1, src2, src3, src4, src5, src6, 128);

        /* Row 0 Row 1 Row 2 Row 3 */
        vec0 = __msa_vshf_b(mask0, src0, src0);
        vec1 = __msa_vshf_b(mask1, src0, src0);
        vec2 = __msa_vshf_b(mask2, src0, src0);
        vec3 = __msa_vshf_b(mask3, src0, src0);

        vec4 = __msa_vshf_b(mask0, src1, src1);
        vec5 = __msa_vshf_b(mask1, src1, src1);
        vec6 = __msa_vshf_b(mask2, src1, src1);
        vec7 = __msa_vshf_b(mask3, src1, src1);

        vec8 = __msa_vshf_b(mask0, src2, src2);
        vec9 = __msa_vshf_b(mask1, src2, src2);
        vec10 = __msa_vshf_b(mask2, src2, src2);
        vec11 = __msa_vshf_b(mask3, src2, src2);

        vec12 = __msa_vshf_b(mask0, src3, src3);
        vec13 = __msa_vshf_b(mask1, src3, src3);
        vec14 = __msa_vshf_b(mask2, src3, src3);
        vec15 = __msa_vshf_b(mask3, src3, src3);

        dst0 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        dst1 = HEVC_FILT_8TAP_DPADD_H(vec4, vec5, vec6, vec7,
                                      filt0, filt1, filt2, filt3, const_vec);

        dst2 = HEVC_FILT_8TAP_DPADD_H(vec8, vec9, vec10, vec11,
                                      filt0, filt1, filt2, filt3, const_vec);

        dst3 = HEVC_FILT_8TAP_DPADD_H(vec12, vec13, vec14, vec15,
                                      filt0, filt1, filt2, filt3, const_vec);

        /* Row 4 Row 5 Row 6 */
        vec0 = __msa_vshf_b(mask0, src4, src4);
        vec1 = __msa_vshf_b(mask1, src4, src4);
        vec2 = __msa_vshf_b(mask2, src4, src4);
        vec3 = __msa_vshf_b(mask3, src4, src4);

        vec4 = __msa_vshf_b(mask0, src5, src5);
        vec5 = __msa_vshf_b(mask1, src5, src5);
        vec6 = __msa_vshf_b(mask2, src5, src5);
        vec7 = __msa_vshf_b(mask3, src5, src5);

        vec8 = __msa_vshf_b(mask0, src6, src6);
        vec9 = __msa_vshf_b(mask1, src6, src6);
        vec10 = __msa_vshf_b(mask2, src6, src6);
        vec11 = __msa_vshf_b(mask3, src6, src6);

        dst4 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                      filt0, filt1, filt2, filt3, const_vec);

        dst5 = HEVC_FILT_8TAP_DPADD_H(vec4, vec5, vec6, vec7,
                                      filt0, filt1, filt2, filt3, const_vec);

        dst6 = HEVC_FILT_8TAP_DPADD_H(vec8, vec9, vec10, vec11,
                                      filt0, filt1, filt2, filt3, const_vec);

        ILVR_H_6VECS_SH(dst0, dst2, dst4, dst1, dst3, dst5,
                        dst1, dst3, dst5, dst2, dst4, dst6,
                        dst10_r, dst32_r, dst54_r, dst21_r, dst43_r, dst65_r);

        ILVL_H_6VECS_SH(dst0, dst2, dst4, dst1, dst3, dst5,
                        dst1, dst3, dst5, dst2, dst4, dst6,
                        dst10_l, dst32_l, dst54_l, dst21_l, dst43_l, dst65_l);

        for (loop_cnt = height >> 1; loop_cnt--;) {
            src7 = LOAD_SB(src_tmp);
            src_tmp += src_stride;

            src7 = (v16i8) __msa_xori_b((v16u8) src7, 128);

            vec0 = __msa_vshf_b(mask0, src7, src7);
            vec1 = __msa_vshf_b(mask1, src7, src7);
            vec2 = __msa_vshf_b(mask2, src7, src7);
            vec3 = __msa_vshf_b(mask3, src7, src7);

            dst7 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                          filt0, filt1, filt2, filt3,
                                          const_vec);

            dst76_r = __msa_ilvr_h(dst7, dst6);
            dst76_l = __msa_ilvl_h(dst7, dst6);

            dst0_r = HEVC_FILT_8TAP_DPADD_W(dst10_r, dst32_r, dst54_r, dst76_r,
                                            filt_h0, filt_h1, filt_h2, filt_h3);

            dst0_l = HEVC_FILT_8TAP_DPADD_W(dst10_l, dst32_l, dst54_l, dst76_l,
                                            filt_h0, filt_h1, filt_h2, filt_h3);

            dst0_r >>= 6;
            dst0_l >>= 6;

            dst0_r = (v4i32) __msa_pckev_h((v8i16) dst0_l, (v8i16) dst0_r);

            STORE_SW(dst0_r, dst_tmp);
            dst_tmp += dst_stride;

            /* Next row */
            src8 = LOAD_SB(src_tmp);
            src_tmp += src_stride;

            src8 = (v16i8) __msa_xori_b((v16u8) src8, 128);

            vec0 = __msa_vshf_b(mask0, src8, src8);
            vec1 = __msa_vshf_b(mask1, src8, src8);
            vec2 = __msa_vshf_b(mask2, src8, src8);
            vec3 = __msa_vshf_b(mask3, src8, src8);

            dst8 = HEVC_FILT_8TAP_DPADD_H(vec0, vec1, vec2, vec3,
                                          filt0, filt1, filt2, filt3,
                                          const_vec);

            dst87_r = __msa_ilvr_h(dst8, dst7);
            dst87_l = __msa_ilvl_h(dst8, dst7);

            dst6 = dst8;

            dst0_r = HEVC_FILT_8TAP_DPADD_W(dst21_r, dst43_r, dst65_r, dst87_r,
                                            filt_h0, filt_h1, filt_h2, filt_h3);

            dst0_l = HEVC_FILT_8TAP_DPADD_W(dst21_l, dst43_l, dst65_l, dst87_l,
                                            filt_h0, filt_h1, filt_h2, filt_h3);

            dst0_r >>= 6;
            dst0_l >>= 6;

            dst0_r = (v4i32) __msa_pckev_h((v8i16) dst0_l, (v8i16) dst0_r);

            STORE_SW(dst0_r, dst_tmp);
            dst_tmp += dst_stride;

            dst10_r = dst32_r;
            dst32_r = dst54_r;
            dst54_r = dst76_r;

            dst10_l = dst32_l;
            dst32_l = dst54_l;
            dst54_l = dst76_l;

            dst21_r = dst43_r;
            dst43_r = dst65_r;
            dst65_r = dst87_r;

            dst21_l = dst43_l;
            dst43_l = dst65_l;
            dst65_l = dst87_l;
        }

        src += 8;
        dst += 8;
    }
}

static void hevc_hv_8t_8w_msa(uint8_t * __restrict src, int32_t src_stride,
                              int16_t * __restrict dst, int32_t dst_stride,
                              const int8_t * __restrict filter_x,
                              const int8_t * __restrict filter_y,
                              int32_t height)
{
    hevc_hv_8t_8multx2mult_msa(src, src_stride, dst, dst_stride,
                               filter_x, filter_y, height, 8);
}

static void hevc_hv_8t_12w_msa(uint8_t * __restrict src, int32_t src_stride,
                               int16_t * __restrict dst, int32_t dst_stride,
                               const int8_t * __restrict filter_x,
                               const int8_t * __restrict filter_y,
                               int32_t height)
{
    hevc_hv_8t_8multx2mult_msa(src, src_stride, dst, dst_stride,
                               filter_x, filter_y, height, 8);

    hevc_hv_8t_4w_msa(src + 8, src_stride, dst + 8, dst_stride,
                      filter_x, filter_y, height);
}

static void hevc_hv_8t_16w_msa(uint8_t * __restrict src, int32_t src_stride,
                               int16_t * __restrict dst, int32_t dst_stride,
                               const int8_t * __restrict filter_x,
                               const int8_t * __restrict filter_y,
                               int32_t height)
{
    hevc_hv_8t_8multx2mult_msa(src, src_stride, dst, dst_stride,
                               filter_x, filter_y, height, 16);
}

static void hevc_hv_8t_24w_msa(uint8_t * __restrict src, int32_t src_stride,
                               int16_t * __restrict dst, int32_t dst_stride,
                               const int8_t * __restrict filter_x,
                               const int8_t * __restrict filter_y,
                               int32_t height)
{
    hevc_hv_8t_8multx2mult_msa(src, src_stride, dst, dst_stride,
                               filter_x, filter_y, height, 24);
}

static void hevc_hv_8t_32w_msa(uint8_t * __restrict src, int32_t src_stride,
                               int16_t * __restrict dst, int32_t dst_stride,
                               const int8_t * __restrict filter_x,
                               const int8_t * __restrict filter_y,
                               int32_t height)
{
    hevc_hv_8t_8multx2mult_msa(src, src_stride, dst, dst_stride,
                               filter_x, filter_y, height, 32);
}

static void hevc_hv_8t_48w_msa(uint8_t * __restrict src, int32_t src_stride,
                               int16_t * __restrict dst, int32_t dst_stride,
                               const int8_t * __restrict filter_x,
                               const int8_t * __restrict filter_y,
                               int32_t height)
{
    hevc_hv_8t_8multx2mult_msa(src, src_stride, dst, dst_stride,
                               filter_x, filter_y, height, 48);
}

static void hevc_hv_8t_64w_msa(uint8_t * __restrict src, int32_t src_stride,
                               int16_t * __restrict dst, int32_t dst_stride,
                               const int8_t * __restrict filter_x,
                               const int8_t * __restrict filter_y,
                               int32_t height)
{
    hevc_hv_8t_8multx2mult_msa(src, src_stride, dst, dst_stride,
                               filter_x, filter_y, height, 64);
}

static uint8_t mc_filt_mask_arr[16 * 3] = {
    /* 8 width cases */
    0, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 7, 8,
    /* 4 width cases */
    0, 1, 1, 2, 2, 3, 3, 4, 16, 17, 17, 18, 18, 19, 19, 20,
    /* 4 width cases */
    8, 9, 9, 10, 10, 11, 11, 12, 24, 25, 25, 26, 26, 27, 27, 28
};

#define FILT_8TAP_DPADD_S_H(vec0, vec1, vec2, vec3,                 \
                            filt0, filt1, filt2, filt3)             \
( {                                                                 \
    v8i16 tmp0, tmp1;                                               \
                                                                    \
    tmp0 = __msa_dotp_s_h((v16i8) (vec0), (v16i8) (filt0));         \
    tmp0 = __msa_dpadd_s_h(tmp0, (v16i8) (vec1), (v16i8) (filt1));  \
    tmp1 = __msa_dotp_s_h((v16i8) (vec2), (v16i8) (filt2));         \
    tmp1 = __msa_dpadd_s_h(tmp1, (v16i8) (vec3), (v16i8) (filt3));  \
    tmp0 = __msa_adds_s_h(tmp0, tmp1);                              \
                                                                    \
    tmp0;                                                           \
} )

#define HORIZ_8TAP_4WID_4VECS_FILT(src0, src1, src2, src3,                   \
                                   mask0, mask1, mask2, mask3,               \
                                   filt0, filt1, filt2, filt3,               \
                                   out0, out1)                               \
{                                                                            \
    v16i8 vec0_m, vec1_m, vec2_m, vec3_m,  vec4_m, vec5_m, vec6_m, vec7_m;   \
    v8i16 res0_m, res1_m, res2_m, res3_m;                                    \
                                                                             \
    vec0_m = __msa_vshf_b((v16i8) (mask0), (v16i8) (src1), (v16i8) (src0));  \
    vec1_m = __msa_vshf_b((v16i8) (mask0), (v16i8) (src3), (v16i8) (src2));  \
                                                                             \
    res0_m = __msa_dotp_s_h(vec0_m, (v16i8) (filt0));                        \
    res1_m = __msa_dotp_s_h(vec1_m, (v16i8) (filt0));                        \
                                                                             \
    vec2_m = __msa_vshf_b((v16i8) (mask1), (v16i8) (src1), (v16i8) (src0));  \
    vec3_m = __msa_vshf_b((v16i8) (mask1), (v16i8) (src3), (v16i8) (src2));  \
                                                                             \
    res0_m = __msa_dpadd_s_h(res0_m, (filt1), vec2_m);                       \
    res1_m = __msa_dpadd_s_h(res1_m, (filt1), vec3_m);                       \
                                                                             \
    vec4_m = __msa_vshf_b((v16i8) (mask2), (v16i8) (src1), (v16i8) (src0));  \
    vec5_m = __msa_vshf_b((v16i8) (mask2), (v16i8) (src3), (v16i8) (src2));  \
                                                                             \
    res2_m = __msa_dotp_s_h((v16i8) (filt2), vec4_m);                        \
    res3_m = __msa_dotp_s_h((v16i8) (filt2), vec5_m);                        \
                                                                             \
    vec6_m = __msa_vshf_b((v16i8) (mask3), (v16i8) (src1), (v16i8) (src0));  \
    vec7_m = __msa_vshf_b((v16i8) (mask3), (v16i8) (src3), (v16i8) (src2));  \
                                                                             \
    res2_m = __msa_dpadd_s_h(res2_m, (v16i8) (filt3), vec6_m);               \
    res3_m = __msa_dpadd_s_h(res3_m, (v16i8) (filt3), vec7_m);               \
                                                                             \
    out0 = __msa_adds_s_h(res0_m, res2_m);                                   \
    out1 = __msa_adds_s_h(res1_m, res3_m);                                   \
}

#define HORIZ_8TAP_8WID_4VECS_FILT(src0, src1, src2, src3,                   \
                                   mask0, mask1, mask2, mask3,               \
                                   filt0, filt1, filt2, filt3,               \
                                   out0, out1, out2, out3)                   \
{                                                                            \
    v16i8 vec0_m, vec1_m, vec2_m, vec3_m, vec4_m, vec5_m, vec6_m, vec7_m;    \
    v8i16 res0_m, res1_m, res2_m, res3_m, res4_m, res5_m, res6_m, res7_m;    \
                                                                             \
    vec0_m = __msa_vshf_b((v16i8) (mask0), (v16i8) (src0), (v16i8) (src0));  \
    vec1_m = __msa_vshf_b((v16i8) (mask0), (v16i8) (src1), (v16i8) (src1));  \
    vec2_m = __msa_vshf_b((v16i8) (mask0), (v16i8) (src2), (v16i8) (src2));  \
    vec3_m = __msa_vshf_b((v16i8) (mask0), (v16i8) (src3), (v16i8) (src3));  \
                                                                             \
    res0_m = __msa_dotp_s_h(vec0_m, (v16i8) (filt0));                        \
    res1_m = __msa_dotp_s_h(vec1_m, (v16i8) (filt0));                        \
    res2_m = __msa_dotp_s_h(vec2_m, (v16i8) (filt0));                        \
    res3_m = __msa_dotp_s_h(vec3_m, (v16i8) (filt0));                        \
                                                                             \
    vec0_m = __msa_vshf_b((v16i8) (mask2), (v16i8) (src0), (v16i8) (src0));  \
    vec1_m = __msa_vshf_b((v16i8) (mask2), (v16i8) (src1), (v16i8) (src1));  \
    vec2_m = __msa_vshf_b((v16i8) (mask2), (v16i8) (src2), (v16i8) (src2));  \
    vec3_m = __msa_vshf_b((v16i8) (mask2), (v16i8) (src3), (v16i8) (src3));  \
                                                                             \
    res4_m = __msa_dotp_s_h(vec0_m, (v16i8) (filt2));                        \
    res5_m = __msa_dotp_s_h(vec1_m, (v16i8) (filt2));                        \
    res6_m = __msa_dotp_s_h(vec2_m, (v16i8) (filt2));                        \
    res7_m = __msa_dotp_s_h(vec3_m, (v16i8) (filt2));                        \
                                                                             \
    vec4_m = __msa_vshf_b((v16i8) (mask1), (v16i8) (src0), (v16i8) (src0));  \
    vec5_m = __msa_vshf_b((v16i8) (mask1), (v16i8) (src1), (v16i8) (src1));  \
    vec6_m = __msa_vshf_b((v16i8) (mask1), (v16i8) (src2), (v16i8) (src2));  \
    vec7_m = __msa_vshf_b((v16i8) (mask1), (v16i8) (src3), (v16i8) (src3));  \
                                                                             \
    res0_m = __msa_dpadd_s_h(res0_m, (v16i8) (filt1), vec4_m);               \
    res1_m = __msa_dpadd_s_h(res1_m, (v16i8) (filt1), vec5_m);               \
    res2_m = __msa_dpadd_s_h(res2_m, (v16i8) (filt1), vec6_m);               \
    res3_m = __msa_dpadd_s_h(res3_m, (v16i8) (filt1), vec7_m);               \
                                                                             \
    vec4_m = __msa_vshf_b((v16i8) (mask3), (v16i8) (src0), (v16i8) (src0));  \
    vec5_m = __msa_vshf_b((v16i8) (mask3), (v16i8) (src1), (v16i8) (src1));  \
    vec6_m = __msa_vshf_b((v16i8) (mask3), (v16i8) (src2), (v16i8) (src2));  \
    vec7_m = __msa_vshf_b((v16i8) (mask3), (v16i8) (src3), (v16i8) (src3));  \
                                                                             \
    res4_m = __msa_dpadd_s_h(res4_m, (v16i8) (filt3), vec4_m);               \
    res5_m = __msa_dpadd_s_h(res5_m, (v16i8) (filt3), vec5_m);               \
    res6_m = __msa_dpadd_s_h(res6_m, (v16i8) (filt3), vec6_m);               \
    res7_m = __msa_dpadd_s_h(res7_m, (v16i8) (filt3), vec7_m);               \
                                                                             \
    out0 = __msa_adds_s_h(res0_m, res4_m);                                   \
    out1 = __msa_adds_s_h(res1_m, res5_m);                                   \
    out2 = __msa_adds_s_h(res2_m, res6_m);                                   \
    out3 = __msa_adds_s_h(res3_m, res7_m);                                   \
}

static void common_hz_8t_4x4_msa(uint8_t *src, int32_t src_stride,
                                 uint8_t *dst, int32_t dst_stride,
                                 const int8_t *filter, uint8_t rnd_val)
{
    v16i8 filt0, filt1, filt2, filt3;
    v16i8 src0, src1, src2, src3;
    v16u8 mask0, mask1, mask2, mask3;
    v8i16 filt, out0, out1;
    v8u16 rnd_vec;

    mask0 = LOAD_UB(&mc_filt_mask_arr[16]);

    src -= 3;

    /* rearranging filter */
    filt = LOAD_SH(filter);
    filt0 = (v16i8) __msa_splati_h(filt, 0);
    filt1 = (v16i8) __msa_splati_h(filt, 1);
    filt2 = (v16i8) __msa_splati_h(filt, 2);
    filt3 = (v16i8) __msa_splati_h(filt, 3);

    mask1 = mask0 + 2;
    mask2 = mask0 + 4;
    mask3 = mask0 + 6;

    LOAD_4VECS_SB(src, src_stride, src0, src1, src2, src3);

    XORI_B_4VECS_SB(src0, src1, src2, src3, src0, src1, src2, src3, 128);

    HORIZ_8TAP_4WID_4VECS_FILT(src0, src1, src2, src3, mask0, mask1, mask2,
                               mask3, filt0, filt1, filt2, filt3, out0, out1);

    rnd_vec = (v8u16) __msa_fill_h(rnd_val);
    out0 = SRAR_SATURATE_SIGNED_H(out0, rnd_vec, 7);
    out1 = SRAR_SATURATE_SIGNED_H(out1, rnd_vec, 7);

    PCKEV_2B_XORI128_STORE_4_BYTES_4(out0, out1, dst, dst_stride);
}

static void common_hz_8t_4x8_msa(uint8_t *src, int32_t src_stride,
                                 uint8_t *dst, int32_t dst_stride,
                                 const int8_t *filter, uint8_t rnd_val)
{
    v16i8 filt0, filt1, filt2, filt3;
    v16i8 src0, src1, src2, src3;
    v16u8 mask0, mask1, mask2, mask3;
    v8i16 filt, out0, out1, out2, out3;
    v8u16 rnd_vec;

    mask0 = LOAD_UB(&mc_filt_mask_arr[16]);

    src -= 3;

    /* rearranging filter */
    filt = LOAD_SH(filter);
    filt0 = (v16i8) __msa_splati_h(filt, 0);
    filt1 = (v16i8) __msa_splati_h(filt, 1);
    filt2 = (v16i8) __msa_splati_h(filt, 2);
    filt3 = (v16i8) __msa_splati_h(filt, 3);

    mask1 = mask0 + 2;
    mask2 = mask0 + 4;
    mask3 = mask0 + 6;

    LOAD_4VECS_SB(src, src_stride, src0, src1, src2, src3);
    src += (4 * src_stride);

    XORI_B_4VECS_SB(src0, src1, src2, src3, src0, src1, src2, src3, 128);

    HORIZ_8TAP_4WID_4VECS_FILT(src0, src1, src2, src3, mask0, mask1, mask2,
                               mask3, filt0, filt1, filt2, filt3, out0, out1);

    LOAD_4VECS_SB(src, src_stride, src0, src1, src2, src3);

    XORI_B_4VECS_SB(src0, src1, src2, src3, src0, src1, src2, src3, 128);

    HORIZ_8TAP_4WID_4VECS_FILT(src0, src1, src2, src3, mask0, mask1, mask2,
                               mask3, filt0, filt1, filt2, filt3, out2, out3);

    rnd_vec = (v8u16) __msa_fill_h(rnd_val);
    out0 = SRAR_SATURATE_SIGNED_H(out0, rnd_vec, 7);
    out1 = SRAR_SATURATE_SIGNED_H(out1, rnd_vec, 7);
    out2 = SRAR_SATURATE_SIGNED_H(out2, rnd_vec, 7);
    out3 = SRAR_SATURATE_SIGNED_H(out3, rnd_vec, 7);

    PCKEV_2B_XORI128_STORE_4_BYTES_4(out0, out1, dst, dst_stride);
    dst += (4 * dst_stride);
    PCKEV_2B_XORI128_STORE_4_BYTES_4(out2, out3, dst, dst_stride);
}

static void common_hz_8t_4x16_msa(uint8_t *src, int32_t src_stride,
                                  uint8_t *dst, int32_t dst_stride,
                                  const int8_t *filter, uint8_t rnd_val)
{
    v16i8 filt0, filt1, filt2, filt3;
    v16i8 src0, src1, src2, src3;
    v16u8 mask0, mask1, mask2, mask3;
    v8i16 filt, out0, out1, out2, out3;
    v8u16 rnd_vec;

    mask0 = LOAD_UB(&mc_filt_mask_arr[16]);

    src -= 3;

    /* rearranging filter */
    filt = LOAD_SH(filter);
    filt0 = (v16i8) __msa_splati_h(filt, 0);
    filt1 = (v16i8) __msa_splati_h(filt, 1);
    filt2 = (v16i8) __msa_splati_h(filt, 2);
    filt3 = (v16i8) __msa_splati_h(filt, 3);

    mask1 = mask0 + 2;
    mask2 = mask0 + 4;
    mask3 = mask0 + 6;

    LOAD_4VECS_SB(src, src_stride, src0, src1, src2, src3);
    src += (4 * src_stride);

    XORI_B_4VECS_SB(src0, src1, src2, src3, src0, src1, src2, src3, 128);

    HORIZ_8TAP_4WID_4VECS_FILT(src0, src1, src2, src3, mask0, mask1, mask2,
                               mask3, filt0, filt1, filt2, filt3, out0, out1);

    LOAD_4VECS_SB(src, src_stride, src0, src1, src2, src3);
    src += (4 * src_stride);

    XORI_B_4VECS_SB(src0, src1, src2, src3, src0, src1, src2, src3, 128);

    HORIZ_8TAP_4WID_4VECS_FILT(src0, src1, src2, src3, mask0, mask1, mask2,
                               mask3, filt0, filt1, filt2, filt3, out2, out3);

    rnd_vec = (v8u16) __msa_fill_h(rnd_val);
    out0 = SRAR_SATURATE_SIGNED_H(out0, rnd_vec, 7);
    out1 = SRAR_SATURATE_SIGNED_H(out1, rnd_vec, 7);
    out2 = SRAR_SATURATE_SIGNED_H(out2, rnd_vec, 7);
    out3 = SRAR_SATURATE_SIGNED_H(out3, rnd_vec, 7);

    PCKEV_2B_XORI128_STORE_4_BYTES_4(out0, out1, dst, dst_stride);
    dst += (4 * dst_stride);
    PCKEV_2B_XORI128_STORE_4_BYTES_4(out2, out3, dst, dst_stride);
    dst += (4 * dst_stride);

    LOAD_4VECS_SB(src, src_stride, src0, src1, src2, src3);
    src += (4 * src_stride);

    XORI_B_4VECS_SB(src0, src1, src2, src3, src0, src1, src2, src3, 128);

    HORIZ_8TAP_4WID_4VECS_FILT(src0, src1, src2, src3, mask0, mask1, mask2,
                               mask3, filt0, filt1, filt2, filt3, out0, out1);

    LOAD_4VECS_SB(src, src_stride, src0, src1, src2, src3);
    src += (4 * src_stride);

    XORI_B_4VECS_SB(src0, src1, src2, src3, src0, src1, src2, src3, 128);

    HORIZ_8TAP_4WID_4VECS_FILT(src0, src1, src2, src3, mask0, mask1, mask2,
                               mask3, filt0, filt1, filt2, filt3, out2, out3);

    rnd_vec = (v8u16) __msa_fill_h(rnd_val);
    out0 = SRAR_SATURATE_SIGNED_H(out0, rnd_vec, 7);
    out1 = SRAR_SATURATE_SIGNED_H(out1, rnd_vec, 7);
    out2 = SRAR_SATURATE_SIGNED_H(out2, rnd_vec, 7);
    out3 = SRAR_SATURATE_SIGNED_H(out3, rnd_vec, 7);

    PCKEV_2B_XORI128_STORE_4_BYTES_4(out0, out1, dst, dst_stride);
    dst += (4 * dst_stride);
    PCKEV_2B_XORI128_STORE_4_BYTES_4(out2, out3, dst, dst_stride);
}

static void common_hz_8t_4w_msa(uint8_t *src, int32_t src_stride,
                                uint8_t *dst, int32_t dst_stride,
                                const int8_t *filter, int32_t height,
                                uint8_t rnd_val)
{
    if (4 == height) {
        common_hz_8t_4x4_msa(src, src_stride, dst, dst_stride, filter, rnd_val);
    } else if (8 == height) {
        common_hz_8t_4x8_msa(src, src_stride, dst, dst_stride, filter, rnd_val);
    } else if (16 == height) {
        common_hz_8t_4x16_msa(src, src_stride, dst, dst_stride, filter,
                              rnd_val);
    }
}

static void common_hz_8t_8x4_msa(uint8_t *src, int32_t src_stride,
                                 uint8_t *dst, int32_t dst_stride,
                                 const int8_t *filter, uint8_t rnd_val)
{
    v16i8 filt0, filt1, filt2, filt3;
    v16i8 src0, src1, src2, src3;
    v16u8 mask0, mask1, mask2, mask3;
    v8i16 filt, out0, out1, out2, out3;
    v8u16 rnd_vec;

    mask0 = LOAD_UB(&mc_filt_mask_arr[0]);

    src -= 3;

    /* rearranging filter */
    filt = LOAD_SH(filter);
    filt0 = (v16i8) __msa_splati_h(filt, 0);
    filt1 = (v16i8) __msa_splati_h(filt, 1);
    filt2 = (v16i8) __msa_splati_h(filt, 2);
    filt3 = (v16i8) __msa_splati_h(filt, 3);

    mask1 = mask0 + 2;
    mask2 = mask0 + 4;
    mask3 = mask0 + 6;

    LOAD_4VECS_SB(src, src_stride, src0, src1, src2, src3);

    XORI_B_4VECS_SB(src0, src1, src2, src3, src0, src1, src2, src3, 128);

    HORIZ_8TAP_8WID_4VECS_FILT(src0, src1, src2, src3, mask0, mask1, mask2,
                               mask3, filt0, filt1, filt2, filt3, out0, out1,
                               out2, out3);

    rnd_vec = (v8u16) __msa_fill_h(rnd_val);
    out0 = SRAR_SATURATE_SIGNED_H(out0, rnd_vec, 7);
    out1 = SRAR_SATURATE_SIGNED_H(out1, rnd_vec, 7);
    out2 = SRAR_SATURATE_SIGNED_H(out2, rnd_vec, 7);
    out3 = SRAR_SATURATE_SIGNED_H(out3, rnd_vec, 7);

    PCKEV_B_4_XORI128_STORE_8_BYTES_4(out0, out1, out2, out3, dst, dst_stride);
}

static void common_hz_8t_8x8mult_msa(uint8_t *src, int32_t src_stride,
                                     uint8_t *dst, int32_t dst_stride,
                                     const int8_t *filter, int32_t height,
                                     uint8_t rnd_val)
{
    uint32_t loop_cnt;
    v16i8 filt0, filt1, filt2, filt3;
    v16i8 src0, src1, src2, src3;
    v16u8 mask0, mask1, mask2, mask3;
    v8i16 filt, out0, out1, out2, out3;
    v8u16 rnd_vec;

    mask0 = LOAD_UB(&mc_filt_mask_arr[0]);

    src -= 3;

    /* rearranging filter */
    filt = LOAD_SH(filter);
    filt0 = (v16i8) __msa_splati_h(filt, 0);
    filt1 = (v16i8) __msa_splati_h(filt, 1);
    filt2 = (v16i8) __msa_splati_h(filt, 2);
    filt3 = (v16i8) __msa_splati_h(filt, 3);

    mask1 = mask0 + 2;
    mask2 = mask0 + 4;
    mask3 = mask0 + 6;

    rnd_vec = (v8u16) __msa_fill_h(rnd_val);

    for (loop_cnt = (height >> 2); loop_cnt--;) {
        LOAD_4VECS_SB(src, src_stride, src0, src1, src2, src3);
        src += (4 * src_stride);

        XORI_B_4VECS_SB(src0, src1, src2, src3, src0, src1, src2, src3, 128);

        HORIZ_8TAP_8WID_4VECS_FILT(src0, src1, src2, src3, mask0, mask1, mask2,
                                   mask3, filt0, filt1, filt2, filt3, out0,
                                   out1, out2, out3);

        out0 = SRAR_SATURATE_SIGNED_H(out0, rnd_vec, 7);
        out1 = SRAR_SATURATE_SIGNED_H(out1, rnd_vec, 7);
        out2 = SRAR_SATURATE_SIGNED_H(out2, rnd_vec, 7);
        out3 = SRAR_SATURATE_SIGNED_H(out3, rnd_vec, 7);

        PCKEV_B_4_XORI128_STORE_8_BYTES_4(out0, out1, out2, out3,
                                          dst, dst_stride);
        dst += (4 * dst_stride);
    }
}

static void common_hz_8t_8w_msa(uint8_t *src, int32_t src_stride,
                                uint8_t *dst, int32_t dst_stride,
                                const int8_t *filter, int32_t height,
                                uint8_t rnd_val)
{
    if (4 == height) {
        common_hz_8t_8x4_msa(src, src_stride, dst, dst_stride, filter, rnd_val);
    } else {
        common_hz_8t_8x8mult_msa(src, src_stride, dst, dst_stride,
                                 filter, height, rnd_val);
    }
}

static void common_hz_8t_12w_msa(uint8_t *src, int32_t src_stride,
                                 uint8_t *dst, int32_t dst_stride,
                                 const int8_t *filter, int32_t height,
                                 uint8_t rnd_val)
{
    uint8_t *src1_ptr, *dst1;
    uint32_t loop_cnt;
    v16i8 filt0, filt1, filt2, filt3;
    v16i8 src0, src1, src2, src3;
    v8u16 rnd_vec;
    v16u8 mask0, mask1, mask2, mask3, mask4, mask5, mask6, mask00;
    v8i16 filt, out0, out1, out2, out3;

    mask00 = LOAD_UB(&mc_filt_mask_arr[0]);
    mask0 = LOAD_UB(&mc_filt_mask_arr[16]);

    src1_ptr = src - 3;
    dst1 = dst;

    dst = dst1 + 8;
    src = src1_ptr + 8;

    /* rearranging filter */
    filt = LOAD_SH(filter);
    filt0 = (v16i8) __msa_splati_h(filt, 0);
    filt1 = (v16i8) __msa_splati_h(filt, 1);
    filt2 = (v16i8) __msa_splati_h(filt, 2);
    filt3 = (v16i8) __msa_splati_h(filt, 3);

    mask1 = mask00 + 2;
    mask2 = mask00 + 4;
    mask3 = mask00 + 6;
    mask4 = mask0 + 2;
    mask5 = mask0 + 4;
    mask6 = mask0 + 6;

    rnd_vec = (v8u16) __msa_fill_h(rnd_val);

    for (loop_cnt = (height >> 2); loop_cnt--;) {
        /* 8 width */
        LOAD_4VECS_SB(src1_ptr, src_stride, src0, src1, src2, src3);
        src1_ptr += (4 * src_stride);

        XORI_B_4VECS_SB(src0, src1, src2, src3, src0, src1, src2, src3, 128);

        HORIZ_8TAP_8WID_4VECS_FILT(src0, src1, src2, src3, mask00, mask1, mask2,
                                   mask3, filt0, filt1, filt2, filt3, out0,
                                   out1, out2, out3);

        out0 = SRAR_SATURATE_SIGNED_H(out0, rnd_vec, 7);
        out1 = SRAR_SATURATE_SIGNED_H(out1, rnd_vec, 7);
        out2 = SRAR_SATURATE_SIGNED_H(out2, rnd_vec, 7);
        out3 = SRAR_SATURATE_SIGNED_H(out3, rnd_vec, 7);

        PCKEV_B_4_XORI128_STORE_8_BYTES_4(out0, out1, out2, out3,
                                          dst1, dst_stride);
        dst1 += (4 * dst_stride);

        /* 4 width */
        LOAD_4VECS_SB(src, src_stride, src0, src1, src2, src3);
        src += (4 * src_stride);

        XORI_B_4VECS_SB(src0, src1, src2, src3, src0, src1, src2, src3, 128);

        HORIZ_8TAP_4WID_4VECS_FILT(src0, src1, src2, src3, mask0, mask4, mask5,
                                   mask6, filt0, filt1, filt2, filt3, out0,
                                   out1);

        out0 = SRAR_SATURATE_SIGNED_H(out0, rnd_vec, 7);
        out1 = SRAR_SATURATE_SIGNED_H(out1, rnd_vec, 7);

        PCKEV_2B_XORI128_STORE_4_BYTES_4(out0, out1, dst, dst_stride);
        dst += (4 * dst_stride);
    }
}

static void common_hz_8t_16w_msa(uint8_t *src, int32_t src_stride,
                                 uint8_t *dst, int32_t dst_stride,
                                 const int8_t *filter, int32_t height,
                                 uint8_t rnd_val)
{
    uint32_t loop_cnt;
    v16i8 src0, src1, src2, src3;
    v16i8 filt0, filt1, filt2, filt3;
    v16u8 mask0, mask1, mask2, mask3;
    v8i16 filt, out0, out1, out2, out3;
    v8u16 rnd_vec;

    mask0 = LOAD_UB(&mc_filt_mask_arr[0]);

    src -= 3;

    /* rearranging filter */
    filt = LOAD_SH(filter);
    filt0 = (v16i8) __msa_splati_h(filt, 0);
    filt1 = (v16i8) __msa_splati_h(filt, 1);
    filt2 = (v16i8) __msa_splati_h(filt, 2);
    filt3 = (v16i8) __msa_splati_h(filt, 3);

    mask1 = mask0 + 2;
    mask2 = mask0 + 4;
    mask3 = mask0 + 6;

    rnd_vec = (v8u16) __msa_fill_h(rnd_val);

    for (loop_cnt = (height >> 1); loop_cnt--;) {
        src0 = LOAD_SB(src);
        src1 = LOAD_SB(src + 8);
        src += src_stride;
        src2 = LOAD_SB(src);
        src3 = LOAD_SB(src + 8);
        src += src_stride;

        XORI_B_4VECS_SB(src0, src1, src2, src3, src0, src1, src2, src3, 128);

        HORIZ_8TAP_8WID_4VECS_FILT(src0, src1, src2, src3, mask0, mask1, mask2,
                                   mask3, filt0, filt1, filt2, filt3, out0,
                                   out1, out2, out3);

        out0 = SRAR_SATURATE_SIGNED_H(out0, rnd_vec, 7);
        out1 = SRAR_SATURATE_SIGNED_H(out1, rnd_vec, 7);
        out2 = SRAR_SATURATE_SIGNED_H(out2, rnd_vec, 7);
        out3 = SRAR_SATURATE_SIGNED_H(out3, rnd_vec, 7);

        PCKEV_B_XORI128_STORE_VEC(out1, out0, dst);
        dst += dst_stride;
        PCKEV_B_XORI128_STORE_VEC(out3, out2, dst);
        dst += dst_stride;
    }
}

static void common_hz_8t_24w_msa(uint8_t *src, int32_t src_stride,
                                 uint8_t *dst, int32_t dst_stride,
                                 const int8_t *filter, int32_t height,
                                 uint8_t rnd_val)
{
    uint32_t loop_cnt;
    v16i8 src0, src1, src2, src3;
    v16i8 filt0, filt1, filt2, filt3;
    v16i8 mask0, mask1, mask2, mask3, mask4, mask5, mask6, mask7;
    v16i8 vec0, vec1, vec2, vec3, vec4, vec5;
    v16i8 vec6, vec7, vec8, vec9, vec10, vec11;
    v8i16 out0, out1, out2, out3, out4, out5;
    v8i16 out6, out7, out8, out9, out10, out11;
    v8i16 filt;
    v8u16 rnd_vec;

    mask0 = LOAD_SB(&mc_filt_mask_arr[0]);

    src -= 3;

    /* rearranging filter */
    filt = LOAD_SH(filter);
    filt0 = (v16i8) __msa_splati_h(filt, 0);
    filt1 = (v16i8) __msa_splati_h(filt, 1);
    filt2 = (v16i8) __msa_splati_h(filt, 2);
    filt3 = (v16i8) __msa_splati_h(filt, 3);

    mask1 = mask0 + 2;
    mask2 = mask0 + 4;
    mask3 = mask0 + 6;
    mask4 = mask0 + 8;
    mask5 = mask0 + 10;
    mask6 = mask0 + 12;
    mask7 = mask0 + 14;

    rnd_vec = (v8u16) __msa_fill_h(rnd_val);

    for (loop_cnt = (height >> 1); loop_cnt--;) {
        src0 = LOAD_SB(src);
        src1 = LOAD_SB(src + 16);
        src += src_stride;
        src2 = LOAD_SB(src);
        src3 = LOAD_SB(src + 16);
        src += src_stride;

        XORI_B_4VECS_SB(src0, src1, src2, src3, src0, src1, src2, src3, 128);

        vec0 = __msa_vshf_b(mask0, src0, src0);
        vec8 = __msa_vshf_b(mask0, src1, src1);
        vec2 = __msa_vshf_b(mask0, src2, src2);
        vec9 = __msa_vshf_b(mask0, src3, src3);
        vec1 = __msa_vshf_b(mask4, src1, src0);
        vec3 = __msa_vshf_b(mask4, src3, src2);

        out0 = __msa_dotp_s_h(vec0, filt0);
        out8 = __msa_dotp_s_h(vec8, filt0);
        out2 = __msa_dotp_s_h(vec2, filt0);
        out9 = __msa_dotp_s_h(vec9, filt0);
        out1 = __msa_dotp_s_h(vec1, filt0);
        out3 = __msa_dotp_s_h(vec3, filt0);

        vec0 = __msa_vshf_b(mask2, src0, src0);
        vec8 = __msa_vshf_b(mask2, src1, src1);
        vec2 = __msa_vshf_b(mask2, src2, src2);
        vec9 = __msa_vshf_b(mask2, src3, src3);
        vec1 = __msa_vshf_b(mask6, src1, src0);
        vec3 = __msa_vshf_b(mask6, src3, src2);

        out4 = __msa_dotp_s_h(vec0, filt2);
        out10 = __msa_dotp_s_h(vec8, filt2);
        out6 = __msa_dotp_s_h(vec2, filt2);
        out11 = __msa_dotp_s_h(vec9, filt2);
        out5 = __msa_dotp_s_h(vec1, filt2);
        out7 = __msa_dotp_s_h(vec3, filt2);

        vec4 = __msa_vshf_b(mask1, src0, src0);
        vec10 = __msa_vshf_b(mask1, src1, src1);
        vec6 = __msa_vshf_b(mask1, src2, src2);
        vec11 = __msa_vshf_b(mask1, src3, src3);
        vec5 = __msa_vshf_b(mask5, src1, src0);
        vec7 = __msa_vshf_b(mask5, src3, src2);

        out0 = __msa_dpadd_s_h(out0, vec4, filt1);
        out8 = __msa_dpadd_s_h(out8, vec10, filt1);
        out2 = __msa_dpadd_s_h(out2, vec6, filt1);
        out9 = __msa_dpadd_s_h(out9, vec11, filt1);
        out1 = __msa_dpadd_s_h(out1, vec5, filt1);
        out3 = __msa_dpadd_s_h(out3, vec7, filt1);

        vec4 = __msa_vshf_b(mask3, src0, src0);
        vec10 = __msa_vshf_b(mask3, src1, src1);
        vec6 = __msa_vshf_b(mask3, src2, src2);
        vec11 = __msa_vshf_b(mask3, src3, src3);
        vec5 = __msa_vshf_b(mask7, src1, src0);
        vec7 = __msa_vshf_b(mask7, src3, src2);

        out4 = __msa_dpadd_s_h(out4, vec4, filt3);
        out10 = __msa_dpadd_s_h(out10, vec10, filt3);
        out6 = __msa_dpadd_s_h(out6, vec6, filt3);
        out11 = __msa_dpadd_s_h(out11, vec11, filt3);
        out5 = __msa_dpadd_s_h(out5, vec5, filt3);
        out7 = __msa_dpadd_s_h(out7, vec7, filt3);

        out0 = __msa_adds_s_h(out0, out4);
        out8 = __msa_adds_s_h(out8, out10);
        out2 = __msa_adds_s_h(out2, out6);
        out9 = __msa_adds_s_h(out9, out11);
        out1 = __msa_adds_s_h(out1, out5);
        out3 = __msa_adds_s_h(out3, out7);

        out0 = SRAR_SATURATE_SIGNED_H(out0, rnd_vec, 7);
        out8 = SRAR_SATURATE_SIGNED_H(out8, rnd_vec, 7);
        out2 = SRAR_SATURATE_SIGNED_H(out2, rnd_vec, 7);
        out9 = SRAR_SATURATE_SIGNED_H(out9, rnd_vec, 7);
        out1 = SRAR_SATURATE_SIGNED_H(out1, rnd_vec, 7);
        out3 = SRAR_SATURATE_SIGNED_H(out3, rnd_vec, 7);

        PCKEV_B_XORI128_STORE_8_BYTES_2(out8, out9, dst + 16, dst_stride);

        PCKEV_B_XORI128_STORE_VEC(out1, out0, dst);
        dst += dst_stride;
        PCKEV_B_XORI128_STORE_VEC(out3, out2, dst);
        dst += dst_stride;
    }
}

static void common_hz_8t_32w_msa(uint8_t *src, int32_t src_stride,
                                 uint8_t *dst, int32_t dst_stride,
                                 const int8_t *filter, int32_t height,
                                 uint8_t rnd_val)
{
    uint32_t loop_cnt;
    v16i8 src0, src1, src2, src3;
    v16i8 filt0, filt1, filt2, filt3;
    v16u8 mask0, mask1, mask2, mask3;
    v8i16 filt, out0, out1, out2, out3;
    v8u16 rnd_vec;

    mask0 = LOAD_UB(&mc_filt_mask_arr[0]);

    src -= 3;

    /* rearranging filter */
    filt = LOAD_SH(filter);
    filt0 = (v16i8) __msa_splati_h(filt, 0);
    filt1 = (v16i8) __msa_splati_h(filt, 1);
    filt2 = (v16i8) __msa_splati_h(filt, 2);
    filt3 = (v16i8) __msa_splati_h(filt, 3);

    mask1 = mask0 + 2;
    mask2 = mask0 + 4;
    mask3 = mask0 + 6;

    rnd_vec = (v8u16) __msa_fill_h(rnd_val);

    for (loop_cnt = (height >> 1); loop_cnt--;) {
        src0 = LOAD_SB(src);
        src2 = LOAD_SB(src + 16);
        src3 = LOAD_SB(src + 24);
        src1 = __msa_sld_b(src2, src0, 8);
        src += src_stride;

        XORI_B_4VECS_SB(src0, src1, src2, src3, src0, src1, src2, src3, 128);

        HORIZ_8TAP_8WID_4VECS_FILT(src0, src1, src2, src3, mask0, mask1, mask2,
                                   mask3, filt0, filt1, filt2, filt3, out0,
                                   out1, out2, out3);

        out0 = SRAR_SATURATE_SIGNED_H(out0, rnd_vec, 7);
        out1 = SRAR_SATURATE_SIGNED_H(out1, rnd_vec, 7);
        out2 = SRAR_SATURATE_SIGNED_H(out2, rnd_vec, 7);
        out3 = SRAR_SATURATE_SIGNED_H(out3, rnd_vec, 7);

        src0 = LOAD_SB(src);
        src2 = LOAD_SB(src + 16);
        src3 = LOAD_SB(src + 24);
        src1 = __msa_sld_b(src2, src0, 8);
        src += src_stride;

        PCKEV_B_XORI128_STORE_VEC(out1, out0, dst);
        PCKEV_B_XORI128_STORE_VEC(out3, out2, (dst + 16));
        dst += dst_stride;

        XORI_B_4VECS_SB(src0, src1, src2, src3, src0, src1, src2, src3, 128);

        HORIZ_8TAP_8WID_4VECS_FILT(src0, src1, src2, src3, mask0, mask1, mask2,
                                   mask3, filt0, filt1, filt2, filt3, out0,
                                   out1, out2, out3);

        out0 = SRAR_SATURATE_SIGNED_H(out0, rnd_vec, 7);
        out1 = SRAR_SATURATE_SIGNED_H(out1, rnd_vec, 7);
        out2 = SRAR_SATURATE_SIGNED_H(out2, rnd_vec, 7);
        out3 = SRAR_SATURATE_SIGNED_H(out3, rnd_vec, 7);

        PCKEV_B_XORI128_STORE_VEC(out1, out0, dst);
        PCKEV_B_XORI128_STORE_VEC(out3, out2, (dst + 16));
        dst += dst_stride;
    }
}

static void common_hz_8t_48w_msa(uint8_t *src, int32_t src_stride,
                                 uint8_t *dst, int32_t dst_stride,
                                 const int8_t *filter, int32_t height,
                                 uint8_t rnd_val)
{
    uint32_t loop_cnt;
    v16i8 src0, src1, src2, src3;
    v16i8 vec0, vec1, vec2;
    v16i8 filt0, filt1, filt2, filt3;
    v16i8 mask0, mask1, mask2, mask3, mask4, mask5, mask6, mask7;
    v8i16 filt, out0, out1, out2, out3, out4, out5, out6;
    v8u16 rnd_vec;

    mask0 = LOAD_SB(&mc_filt_mask_arr[0]);

    /* rearranging filter */
    filt = LOAD_SH(filter);
    filt0 = (v16i8) __msa_splati_h(filt, 0);
    filt1 = (v16i8) __msa_splati_h(filt, 1);
    filt2 = (v16i8) __msa_splati_h(filt, 2);
    filt3 = (v16i8) __msa_splati_h(filt, 3);

    mask1 = mask0 + 2;
    mask2 = mask0 + 4;
    mask3 = mask0 + 6;
    mask4 = mask0 + 8;
    mask5 = mask0 + 10;
    mask6 = mask0 + 12;
    mask7 = mask0 + 14;

    rnd_vec = (v8u16) __msa_fill_h(rnd_val);

    src -= 3;

    for (loop_cnt = height; loop_cnt--;) {
        src0 = LOAD_SB(src);
        src2 = LOAD_SB(src + 16);
        src3 = LOAD_SB(src + 32);
        src1 = __msa_sld_b(src2, src0, 8);

        XORI_B_4VECS_SB(src0, src1, src2, src3, src0, src1, src2, src3, 128);

        vec0 = __msa_vshf_b(mask0, src0, src0);
        vec1 = __msa_vshf_b(mask0, src1, src1);
        vec2 = __msa_vshf_b(mask0, src2, src2);
        out0 = __msa_dotp_s_h(vec0, filt0);
        out1 = __msa_dotp_s_h(vec1, filt0);
        out2 = __msa_dotp_s_h(vec2, filt0);

        vec0 = __msa_vshf_b(mask1, src0, src0);
        vec1 = __msa_vshf_b(mask1, src1, src1);
        vec2 = __msa_vshf_b(mask1, src2, src2);
        out0 = __msa_dpadd_s_h(out0, vec0, filt1);
        out1 = __msa_dpadd_s_h(out1, vec1, filt1);
        out2 = __msa_dpadd_s_h(out2, vec2, filt1);

        vec0 = __msa_vshf_b(mask2, src0, src0);
        vec1 = __msa_vshf_b(mask2, src1, src1);
        vec2 = __msa_vshf_b(mask2, src2, src2);
        out3 = __msa_dotp_s_h(vec0, filt2);
        out4 = __msa_dotp_s_h(vec1, filt2);
        out5 = __msa_dotp_s_h(vec2, filt2);

        vec0 = __msa_vshf_b(mask3, src0, src0);
        vec1 = __msa_vshf_b(mask3, src1, src1);
        vec2 = __msa_vshf_b(mask3, src2, src2);
        out3 = __msa_dpadd_s_h(out3, vec0, filt3);
        out4 = __msa_dpadd_s_h(out4, vec1, filt3);
        out5 = __msa_dpadd_s_h(out5, vec2, filt3);

        out0 = __msa_adds_s_h(out0, out3);
        out1 = __msa_adds_s_h(out1, out4);
        out2 = __msa_adds_s_h(out2, out5);

        out0 = SRAR_SATURATE_SIGNED_H(out0, rnd_vec, 7);
        out1 = SRAR_SATURATE_SIGNED_H(out1, rnd_vec, 7);
        out6 = SRAR_SATURATE_SIGNED_H(out2, rnd_vec, 7);

        PCKEV_B_XORI128_STORE_VEC(out1, out0, dst);

        src1 = LOAD_SB(src + 40);
        src1 = (v16i8) __msa_xori_b((v16u8) src1, 128);

        vec0 = __msa_vshf_b(mask4, src3, src2);
        vec1 = __msa_vshf_b(mask0, src3, src3);
        vec2 = __msa_vshf_b(mask0, src1, src1);
        out0 = __msa_dotp_s_h(vec0, filt0);
        out1 = __msa_dotp_s_h(vec1, filt0);
        out2 = __msa_dotp_s_h(vec2, filt0);

        vec0 = __msa_vshf_b(mask5, src3, src2);
        vec1 = __msa_vshf_b(mask1, src3, src3);
        vec2 = __msa_vshf_b(mask1, src1, src1);
        out0 = __msa_dpadd_s_h(out0, vec0, filt1);
        out1 = __msa_dpadd_s_h(out1, vec1, filt1);
        out2 = __msa_dpadd_s_h(out2, vec2, filt1);

        vec0 = __msa_vshf_b(mask6, src3, src2);
        vec1 = __msa_vshf_b(mask2, src3, src3);
        vec2 = __msa_vshf_b(mask2, src1, src1);
        out3 = __msa_dotp_s_h(vec0, filt2);
        out4 = __msa_dotp_s_h(vec1, filt2);
        out5 = __msa_dotp_s_h(vec2, filt2);

        vec0 = __msa_vshf_b(mask7, src3, src2);
        vec1 = __msa_vshf_b(mask3, src3, src3);
        vec2 = __msa_vshf_b(mask3, src1, src1);
        out3 = __msa_dpadd_s_h(out3, vec0, filt3);
        out4 = __msa_dpadd_s_h(out4, vec1, filt3);
        out5 = __msa_dpadd_s_h(out5, vec2, filt3);

        out3 = __msa_adds_s_h(out0, out3);
        out4 = __msa_adds_s_h(out1, out4);
        out5 = __msa_adds_s_h(out2, out5);

        out3 = SRAR_SATURATE_SIGNED_H(out3, rnd_vec, 7);
        out4 = SRAR_SATURATE_SIGNED_H(out4, rnd_vec, 7);
        out5 = SRAR_SATURATE_SIGNED_H(out5, rnd_vec, 7);

        PCKEV_B_XORI128_STORE_VEC(out3, out6, (dst + 16));
        PCKEV_B_XORI128_STORE_VEC(out5, out4, (dst + 32));

        src += src_stride;
        dst += dst_stride;
    }
}

static void common_hz_8t_64w_msa(uint8_t *src, int32_t src_stride,
                                 uint8_t *dst, int32_t dst_stride,
                                 const int8_t *filter, int32_t height,
                                 uint8_t rnd_val)
{
    int32_t loop_cnt;
    v16i8 src0, src1, src2, src3;
    v16i8 filt0, filt1, filt2, filt3;
    v16u8 mask0, mask1, mask2, mask3;
    v8i16 filt, out0, out1, out2, out3;
    v8u16 rnd_vec;

    mask0 = LOAD_UB(&mc_filt_mask_arr[0]);

    src -= 3;

    /* rearranging filter */
    filt = LOAD_SH(filter);
    filt0 = (v16i8) __msa_splati_h(filt, 0);
    filt1 = (v16i8) __msa_splati_h(filt, 1);
    filt2 = (v16i8) __msa_splati_h(filt, 2);
    filt3 = (v16i8) __msa_splati_h(filt, 3);

    mask1 = mask0 + 2;
    mask2 = mask0 + 4;
    mask3 = mask0 + 6;

    rnd_vec = (v8u16) __msa_fill_h(rnd_val);

    for (loop_cnt = height; loop_cnt--;) {
        src0 = LOAD_SB(src);
        src2 = LOAD_SB(src + 16);
        src3 = LOAD_SB(src + 24);
        src1 = __msa_sld_b(src2, src0, 8);

        XORI_B_4VECS_SB(src0, src1, src2, src3, src0, src1, src2, src3, 128);

        HORIZ_8TAP_8WID_4VECS_FILT(src0, src1, src2, src3, mask0, mask1,
                                   mask2, mask3, filt0, filt1, filt2, filt3,
                                   out0, out1, out2, out3);

        out0 = SRAR_SATURATE_SIGNED_H(out0, rnd_vec, 7);
        out1 = SRAR_SATURATE_SIGNED_H(out1, rnd_vec, 7);
        out2 = SRAR_SATURATE_SIGNED_H(out2, rnd_vec, 7);
        out3 = SRAR_SATURATE_SIGNED_H(out3, rnd_vec, 7);

        PCKEV_B_XORI128_STORE_VEC(out1, out0, dst);
        PCKEV_B_XORI128_STORE_VEC(out3, out2, dst + 16);

        src0 = LOAD_SB(src + 32);
        src2 = LOAD_SB(src + 48);
        src3 = LOAD_SB(src + 56);
        src1 = __msa_sld_b(src2, src0, 8);

        XORI_B_4VECS_SB(src0, src1, src2, src3, src0, src1, src2, src3, 128);

        HORIZ_8TAP_8WID_4VECS_FILT(src0, src1, src2, src3, mask0, mask1,
                                   mask2, mask3, filt0, filt1, filt2, filt3,
                                   out0, out1, out2, out3);

        out0 = SRAR_SATURATE_SIGNED_H(out0, rnd_vec, 7);
        out1 = SRAR_SATURATE_SIGNED_H(out1, rnd_vec, 7);
        out2 = SRAR_SATURATE_SIGNED_H(out2, rnd_vec, 7);
        out3 = SRAR_SATURATE_SIGNED_H(out3, rnd_vec, 7);

        PCKEV_B_XORI128_STORE_VEC(out1, out0, dst + 32);
        PCKEV_B_XORI128_STORE_VEC(out3, out2, dst + 48);

        src += src_stride;
        dst += dst_stride;
    }
}

static void common_vt_8t_4w_msa(uint8_t *src, int32_t src_stride,
                                uint8_t *dst, int32_t dst_stride,
                                const int8_t *filter, int32_t height,
                                uint8_t rnd_val)
{
    uint32_t loop_cnt;
    v16i8 src0, src1, src2, src3, src4, src5, src6, src7, src8, src9, src10;
    v16i8 src10_r, src32_r, src54_r, src76_r, src98_r;
    v16i8 src21_r, src43_r, src65_r, src87_r, src109_r;
    v16i8 src2110, src4332, src6554, src8776, src10998;
    v16i8 filt0, filt1, filt2, filt3;
    v8i16 filt, out10, out32;
    v8u16 rnd_vec;

    src -= (3 * src_stride);

    filt = LOAD_SH(filter);
    filt0 = (v16i8) __msa_splati_h(filt, 0);
    filt1 = (v16i8) __msa_splati_h(filt, 1);
    filt2 = (v16i8) __msa_splati_h(filt, 2);
    filt3 = (v16i8) __msa_splati_h(filt, 3);

    LOAD_7VECS_SB(src, src_stride, src0, src1, src2, src3, src4, src5, src6);
    src += (7 * src_stride);

    ILVR_B_6VECS_SB(src0, src2, src4, src1, src3, src5,
                    src1, src3, src5, src2, src4, src6,
                    src10_r, src32_r, src54_r, src21_r, src43_r, src65_r);

    ILVR_D_3VECS_SB(src2110, src21_r, src10_r, src4332, src43_r, src32_r,
                    src6554, src65_r, src54_r);

    XORI_B_3VECS_SB(src2110, src4332, src6554, src2110, src4332, src6554, 128);

    rnd_vec = (v8u16) __msa_fill_h(rnd_val);

    for (loop_cnt = (height >> 2); loop_cnt--;) {
        LOAD_4VECS_SB(src, src_stride, src7, src8, src9, src10);
        src += (4 * src_stride);

        ILVR_B_4VECS_SB(src6, src7, src8, src9, src7, src8, src9, src10,
                        src76_r, src87_r, src98_r, src109_r);

        ILVR_D_2VECS_SB(src8776, src87_r, src76_r, src10998, src109_r, src98_r);

        XORI_B_2VECS_SB(src8776, src10998, src8776, src10998, 128);

        out10 = FILT_8TAP_DPADD_S_H(src2110, src4332, src6554, src8776,
                                    filt0, filt1, filt2, filt3);
        out32 = FILT_8TAP_DPADD_S_H(src4332, src6554, src8776, src10998,
                                    filt0, filt1, filt2, filt3);

        out10 = SRAR_SATURATE_SIGNED_H(out10, rnd_vec, 7);
        out32 = SRAR_SATURATE_SIGNED_H(out32, rnd_vec, 7);

        PCKEV_2B_XORI128_STORE_4_BYTES_4(out10, out32, dst, dst_stride);
        dst += (4 * dst_stride);

        src2110 = src6554;
        src4332 = src8776;
        src6554 = src10998;

        src6 = src10;
    }
}

static void common_vt_8t_8w_msa(uint8_t *src, int32_t src_stride,
                                uint8_t *dst, int32_t dst_stride,
                                const int8_t *filter, int32_t height,
                                uint8_t rnd_val)
{
    uint32_t loop_cnt;
    v16i8 src0, src1, src2, src3, src4, src5, src6, src7, src8, src9, src10;
    v16i8 src10_r, src32_r, src54_r, src76_r, src98_r;
    v16i8 src21_r, src43_r, src65_r, src87_r, src109_r;
    v16i8 filt0, filt1, filt2, filt3;
    v8i16 filt, out0_r, out1_r, out2_r, out3_r;
    v8u16 rnd_vec;

    src -= (3 * src_stride);

    rnd_vec = (v8u16) __msa_fill_h(rnd_val);

    filt = LOAD_SH(filter);
    filt0 = (v16i8) __msa_splati_h(filt, 0);
    filt1 = (v16i8) __msa_splati_h(filt, 1);
    filt2 = (v16i8) __msa_splati_h(filt, 2);
    filt3 = (v16i8) __msa_splati_h(filt, 3);

    LOAD_7VECS_SB(src, src_stride, src0, src1, src2, src3, src4, src5, src6);
    src += (7 * src_stride);

    XORI_B_7VECS_SB(src0, src1, src2, src3, src4, src5, src6,
                    src0, src1, src2, src3, src4, src5, src6, 128);

    ILVR_B_6VECS_SB(src0, src2, src4, src1, src3, src5,
                    src1, src3, src5, src2, src4, src6,
                    src10_r, src32_r, src54_r, src21_r, src43_r, src65_r);

    for (loop_cnt = (height >> 2); loop_cnt--;) {
        LOAD_4VECS_SB(src, src_stride, src7, src8, src9, src10);
        src += (4 * src_stride);

        XORI_B_4VECS_SB(src7, src8, src9, src10, src7, src8, src9, src10, 128);

        ILVR_B_4VECS_SB(src6, src7, src8, src9, src7, src8, src9, src10,
                        src76_r, src87_r, src98_r, src109_r);

        out0_r = FILT_8TAP_DPADD_S_H(src10_r, src32_r, src54_r, src76_r,
                                     filt0, filt1, filt2, filt3);
        out1_r = FILT_8TAP_DPADD_S_H(src21_r, src43_r, src65_r, src87_r,
                                     filt0, filt1, filt2, filt3);
        out2_r = FILT_8TAP_DPADD_S_H(src32_r, src54_r, src76_r, src98_r,
                                     filt0, filt1, filt2, filt3);
        out3_r = FILT_8TAP_DPADD_S_H(src43_r, src65_r, src87_r, src109_r,
                                     filt0, filt1, filt2, filt3);

        out0_r = SRAR_SATURATE_SIGNED_H(out0_r, rnd_vec, 7);
        out1_r = SRAR_SATURATE_SIGNED_H(out1_r, rnd_vec, 7);
        out2_r = SRAR_SATURATE_SIGNED_H(out2_r, rnd_vec, 7);
        out3_r = SRAR_SATURATE_SIGNED_H(out3_r, rnd_vec, 7);

        PCKEV_B_4_XORI128_STORE_8_BYTES_4(out0_r, out1_r, out2_r, out3_r,
                                          dst, dst_stride);
        dst += (4 * dst_stride);

        src10_r = src54_r;
        src32_r = src76_r;
        src54_r = src98_r;
        src21_r = src65_r;
        src43_r = src87_r;
        src65_r = src109_r;

        src6 = src10;
    }
}

static void common_vt_8t_12w_msa(uint8_t *src, int32_t src_stride,
                                 uint8_t *dst, int32_t dst_stride,
                                 const int8_t *filter, int32_t height,
                                 uint8_t rnd_val)
{
    int32_t loop_cnt;
    uint32_t out2, out3;
    uint64_t out0, out1;
    v16u8 src0, src1, src2, src3, src4, src5, src6, src7, src8;
    v16i8 vec0, vec1, vec2, vec3, vec4, vec5, vec6, vec7;
    v16i8 vec8, vec9, vec10, vec11, vec12, vec13, vec14, vec15, vec16;
    v16i8 res0, res1, res2;
    v8i16 vec01, vec23, vec45, vec67;
    v8i16 tmp0, tmp1, tmp2;
    v8i16 filt, filt0, filt1, filt2, filt3;
    v8u16 rnd_vec;
    v4i32 mask = { 2, 6, 2, 6 };

    src -= (3 * src_stride);

    LOAD_7VECS_UB(src, src_stride, src0, src1, src2, src3, src4, src5, src6);
    src += (7 * src_stride);

    XORI_B_4VECS_SB(src0, src1, src2, src3, vec0, vec1, vec2, vec3, 128);
    vec4 = (v16i8) __msa_xori_b((v16u8) src4, 128);
    vec5 = (v16i8) __msa_xori_b((v16u8) src5, 128);
    vec6 = (v16i8) __msa_xori_b((v16u8) src6, 128);

    /* 4 width */
    vec9 = (v16i8) __msa_vshf_w(mask, (v4i32) vec1, (v4i32) vec0);
    vec10 = (v16i8) __msa_vshf_w(mask, (v4i32) vec2, (v4i32) vec1);
    vec11 = (v16i8) __msa_vshf_w(mask, (v4i32) vec3, (v4i32) vec2);
    vec12 = (v16i8) __msa_vshf_w(mask, (v4i32) vec4, (v4i32) vec3);
    vec13 = (v16i8) __msa_vshf_w(mask, (v4i32) vec5, (v4i32) vec4);
    vec14 = (v16i8) __msa_vshf_w(mask, (v4i32) vec6, (v4i32) vec5);

    rnd_vec = (v8u16) __msa_fill_h(rnd_val);

    /* rearranging filter_y */
    filt = LOAD_SH(filter);
    filt0 = (v8i16) __msa_splati_h(filt, 0);
    filt1 = (v8i16) __msa_splati_h(filt, 1);
    filt2 = (v8i16) __msa_splati_h(filt, 2);
    filt3 = (v8i16) __msa_splati_h(filt, 3);

    for (loop_cnt = (height >> 1); loop_cnt--;) {
        LOAD_2VECS_UB(src, src_stride, src7, src8);
        src += (2 * src_stride);

        XORI_B_2VECS_SB(src7, src8, vec7, vec8, 128);

        ILVR_B_4VECS_SH(vec0, vec2, vec4, vec6, vec1, vec3, vec5, vec7,
                        vec01, vec23, vec45, vec67);

        tmp0 = FILT_8TAP_DPADD_S_H(vec01, vec23, vec45, vec67, filt0, filt1,
                                   filt2, filt3);

        ILVR_B_4VECS_SH(vec1, vec3, vec5, vec7, vec2, vec4, vec6, vec8,
                        vec01, vec23, vec45, vec67);

        tmp1 = FILT_8TAP_DPADD_S_H(vec01, vec23, vec45, vec67, filt0, filt1,
                                   filt2, filt3);

        /* 4 width */
        vec15 = (v16i8) __msa_vshf_w(mask, (v4i32) vec7, (v4i32) vec6);
        vec16 = (v16i8) __msa_vshf_w(mask, (v4i32) vec8, (v4i32) vec7);

        ILVR_B_4VECS_SH(vec9, vec11, vec13, vec15, vec10, vec12, vec14, vec16,
                        vec01, vec23, vec45, vec67);

        tmp2 = FILT_8TAP_DPADD_S_H(vec01, vec23, vec45, vec67, filt0, filt1,
                                   filt2, filt3);

        tmp0 = SRAR_SATURATE_SIGNED_H(tmp0, rnd_vec, 7);
        tmp1 = SRAR_SATURATE_SIGNED_H(tmp1, rnd_vec, 7);
        tmp2 = SRAR_SATURATE_SIGNED_H(tmp2, rnd_vec, 7);

        res0 = __msa_pckev_b((v16i8) tmp0, (v16i8) tmp0);
        res1 = __msa_pckev_b((v16i8) tmp1, (v16i8) tmp1);
        res2 = __msa_pckev_b((v16i8) tmp2, (v16i8) tmp2);

        XORI_B_3VECS_SB(res0, res1, res2, res0, res1, res2, 128);

        out0 = __msa_copy_u_d((v2i64) res0, 0);
        out1 = __msa_copy_u_d((v2i64) res1, 0);
        out2 = __msa_copy_u_w((v4i32) res2, 0);
        out3 = __msa_copy_u_w((v4i32) res2, 1);

        STORE_DWORD(dst, out0);
        STORE_WORD((dst + 8), out2);
        dst += dst_stride;
        STORE_DWORD(dst, out1);
        STORE_WORD((dst + 8), out3);
        dst += dst_stride;

        vec0 = vec2;
        vec1 = vec3;
        vec2 = vec4;
        vec3 = vec5;
        vec4 = vec6;
        vec5 = vec7;
        vec6 = vec8;
        vec9 = vec11;
        vec10 = vec12;
        vec11 = vec13;
        vec12 = vec14;
        vec13 = vec15;
        vec14 = vec16;
    }
}

static void common_vt_8t_16w_msa(uint8_t *src, int32_t src_stride,
                                 uint8_t *dst, int32_t dst_stride,
                                 const int8_t *filter, int32_t height,
                                 uint8_t rnd_val)
{
    uint32_t loop_cnt;
    v16i8 src0, src1, src2, src3, src4, src5, src6, src7, src8, src9, src10;
    v16i8 filt0, filt1, filt2, filt3;
    v16i8 src10_r, src32_r, src54_r, src76_r, src98_r;
    v16i8 src21_r, src43_r, src65_r, src87_r, src109_r;
    v16i8 src10_l, src32_l, src54_l, src76_l, src98_l;
    v16i8 src21_l, src43_l, src65_l, src87_l, src109_l;
    v8i16 out0_r, out1_r, out2_r, out3_r, out0_l, out1_l, out2_l, out3_l;
    v8i16 filt;
    v8u16 rnd_vec;
    v16u8 tmp0, tmp1, tmp2, tmp3;

    src -= (3 * src_stride);

    rnd_vec = (v8u16) __msa_fill_h(rnd_val);

    filt = LOAD_SH(filter);
    filt0 = (v16i8) __msa_splati_h(filt, 0);
    filt1 = (v16i8) __msa_splati_h(filt, 1);
    filt2 = (v16i8) __msa_splati_h(filt, 2);
    filt3 = (v16i8) __msa_splati_h(filt, 3);

    LOAD_7VECS_SB(src, src_stride,
                  src0, src1, src2, src3, src4, src5, src6);
    src += (7 * src_stride);

    XORI_B_7VECS_SB(src0, src1, src2, src3, src4, src5, src6,
                    src0, src1, src2, src3, src4, src5, src6, 128);

    ILVR_B_6VECS_SB(src0, src2, src4, src1, src3, src5,
                    src1, src3, src5, src2, src4, src6,
                    src10_r, src32_r, src54_r, src21_r, src43_r, src65_r);

    ILVL_B_6VECS_SB(src0, src2, src4, src1, src3, src5,
                    src1, src3, src5, src2, src4, src6,
                    src10_l, src32_l, src54_l, src21_l, src43_l, src65_l);

    for (loop_cnt = (height >> 2); loop_cnt--;) {
        LOAD_4VECS_SB(src, src_stride, src7, src8, src9, src10);
        src += (4 * src_stride);

        XORI_B_4VECS_SB(src7, src8, src9, src10,
                        src7, src8, src9, src10, 128);

        ILVR_B_4VECS_SB(src6, src7, src8, src9, src7, src8, src9, src10,
                        src76_r, src87_r, src98_r, src109_r);

        ILVL_B_4VECS_SB(src6, src7, src8, src9, src7, src8, src9, src10,
                        src76_l, src87_l, src98_l, src109_l);

        out0_r = FILT_8TAP_DPADD_S_H(src10_r, src32_r, src54_r, src76_r,
                                     filt0, filt1, filt2, filt3);
        out1_r = FILT_8TAP_DPADD_S_H(src21_r, src43_r, src65_r, src87_r,
                                     filt0, filt1, filt2, filt3);
        out2_r = FILT_8TAP_DPADD_S_H(src32_r, src54_r, src76_r, src98_r,
                                     filt0, filt1, filt2, filt3);
        out3_r = FILT_8TAP_DPADD_S_H(src43_r, src65_r, src87_r, src109_r,
                                     filt0, filt1, filt2, filt3);

        out0_l = FILT_8TAP_DPADD_S_H(src10_l, src32_l, src54_l, src76_l,
                                     filt0, filt1, filt2, filt3);
        out1_l = FILT_8TAP_DPADD_S_H(src21_l, src43_l, src65_l, src87_l,
                                     filt0, filt1, filt2, filt3);
        out2_l = FILT_8TAP_DPADD_S_H(src32_l, src54_l, src76_l, src98_l,
                                     filt0, filt1, filt2, filt3);
        out3_l = FILT_8TAP_DPADD_S_H(src43_l, src65_l, src87_l, src109_l,
                                     filt0, filt1, filt2, filt3);

        out0_r = SRAR_SATURATE_SIGNED_H(out0_r, rnd_vec, 7);
        out1_r = SRAR_SATURATE_SIGNED_H(out1_r, rnd_vec, 7);
        out2_r = SRAR_SATURATE_SIGNED_H(out2_r, rnd_vec, 7);
        out3_r = SRAR_SATURATE_SIGNED_H(out3_r, rnd_vec, 7);
        out0_l = SRAR_SATURATE_SIGNED_H(out0_l, rnd_vec, 7);
        out1_l = SRAR_SATURATE_SIGNED_H(out1_l, rnd_vec, 7);
        out2_l = SRAR_SATURATE_SIGNED_H(out2_l, rnd_vec, 7);
        out3_l = SRAR_SATURATE_SIGNED_H(out3_l, rnd_vec, 7);

        PCKEV_B_4VECS_UB(out0_l, out1_l, out2_l, out3_l, out0_r, out1_r, out2_r,
                         out3_r, tmp0, tmp1, tmp2, tmp3);

        XORI_B_4VECS_UB(tmp0, tmp1, tmp2, tmp3, tmp0, tmp1, tmp2, tmp3, 128);

        STORE_4VECS_UB(dst, dst_stride, tmp0, tmp1, tmp2, tmp3);
        dst += (4 * dst_stride);

        src10_r = src54_r;
        src32_r = src76_r;
        src54_r = src98_r;
        src21_r = src65_r;
        src43_r = src87_r;
        src65_r = src109_r;
        src10_l = src54_l;
        src32_l = src76_l;
        src54_l = src98_l;
        src21_l = src65_l;
        src43_l = src87_l;
        src65_l = src109_l;

        src6 = src10;
    }
}

static void common_vt_8t_16w_mult_msa(uint8_t *src, int32_t src_stride,
                                      uint8_t *dst, int32_t dst_stride,
                                      const int8_t *filter, int32_t height,
                                      uint8_t rnd_val, int32_t width)
{
    uint8_t *src_tmp;
    uint8_t *dst_tmp;
    uint32_t loop_cnt, cnt;
    v16i8 src0, src1, src2, src3, src4, src5, src6, src7, src8, src9, src10;
    v16i8 filt0, filt1, filt2, filt3;
    v16i8 src10_r, src32_r, src54_r, src76_r, src98_r;
    v16i8 src21_r, src43_r, src65_r, src87_r, src109_r;
    v16i8 src10_l, src32_l, src54_l, src76_l, src98_l;
    v16i8 src21_l, src43_l, src65_l, src87_l, src109_l;
    v8i16 filt, out0_r, out1_r, out2_r, out3_r, out0_l, out1_l, out2_l, out3_l;
    v16u8 tmp0, tmp1, tmp2, tmp3;
    v8u16 rnd_vec;

    src -= (3 * src_stride);

    rnd_vec = (v8u16) __msa_fill_h(rnd_val);

    filt = LOAD_SH(filter);
    filt0 = (v16i8) __msa_splati_h(filt, 0);
    filt1 = (v16i8) __msa_splati_h(filt, 1);
    filt2 = (v16i8) __msa_splati_h(filt, 2);
    filt3 = (v16i8) __msa_splati_h(filt, 3);

    for (cnt = (width >> 4); cnt--;) {
        src_tmp = src;
        dst_tmp = dst;

        LOAD_7VECS_SB(src_tmp, src_stride,
                      src0, src1, src2, src3, src4, src5, src6);
        src_tmp += (7 * src_stride);

        XORI_B_7VECS_SB(src0, src1, src2, src3, src4, src5, src6,
                        src0, src1, src2, src3, src4, src5, src6, 128);

        ILVR_B_6VECS_SB(src0, src2, src4, src1, src3, src5,
                        src1, src3, src5, src2, src4, src6,
                        src10_r, src32_r, src54_r, src21_r, src43_r, src65_r);

        ILVL_B_6VECS_SB(src0, src2, src4, src1, src3, src5,
                        src1, src3, src5, src2, src4, src6,
                        src10_l, src32_l, src54_l, src21_l, src43_l, src65_l);

        for (loop_cnt = (height >> 2); loop_cnt--;) {
            LOAD_4VECS_SB(src_tmp, src_stride, src7, src8, src9, src10);
            src_tmp += (4 * src_stride);

            XORI_B_4VECS_SB(src7, src8, src9, src10,
                            src7, src8, src9, src10, 128);

            ILVR_B_4VECS_SB(src6, src7, src8, src9, src7, src8, src9, src10,
                            src76_r, src87_r, src98_r, src109_r);

            ILVL_B_4VECS_SB(src6, src7, src8, src9, src7, src8, src9, src10,
                            src76_l, src87_l, src98_l, src109_l);

            out0_r = FILT_8TAP_DPADD_S_H(src10_r, src32_r, src54_r, src76_r,
                                         filt0, filt1, filt2, filt3);
            out1_r = FILT_8TAP_DPADD_S_H(src21_r, src43_r, src65_r, src87_r,
                                         filt0, filt1, filt2, filt3);
            out2_r = FILT_8TAP_DPADD_S_H(src32_r, src54_r, src76_r, src98_r,
                                         filt0, filt1, filt2, filt3);
            out3_r = FILT_8TAP_DPADD_S_H(src43_r, src65_r, src87_r, src109_r,
                                         filt0, filt1, filt2, filt3);

            out0_l = FILT_8TAP_DPADD_S_H(src10_l, src32_l, src54_l, src76_l,
                                         filt0, filt1, filt2, filt3);
            out1_l = FILT_8TAP_DPADD_S_H(src21_l, src43_l, src65_l, src87_l,
                                         filt0, filt1, filt2, filt3);
            out2_l = FILT_8TAP_DPADD_S_H(src32_l, src54_l, src76_l, src98_l,
                                         filt0, filt1, filt2, filt3);
            out3_l = FILT_8TAP_DPADD_S_H(src43_l, src65_l, src87_l, src109_l,
                                         filt0, filt1, filt2, filt3);

            out0_r = SRAR_SATURATE_SIGNED_H(out0_r, rnd_vec, 7);
            out1_r = SRAR_SATURATE_SIGNED_H(out1_r, rnd_vec, 7);
            out2_r = SRAR_SATURATE_SIGNED_H(out2_r, rnd_vec, 7);
            out3_r = SRAR_SATURATE_SIGNED_H(out3_r, rnd_vec, 7);
            out0_l = SRAR_SATURATE_SIGNED_H(out0_l, rnd_vec, 7);
            out1_l = SRAR_SATURATE_SIGNED_H(out1_l, rnd_vec, 7);
            out2_l = SRAR_SATURATE_SIGNED_H(out2_l, rnd_vec, 7);
            out3_l = SRAR_SATURATE_SIGNED_H(out3_l, rnd_vec, 7);

            PCKEV_B_4VECS_UB(out0_l, out1_l, out2_l, out3_l, out0_r, out1_r,
                             out2_r, out3_r, tmp0, tmp1, tmp2, tmp3);

            XORI_B_4VECS_UB(tmp0, tmp1, tmp2, tmp3,
                            tmp0, tmp1, tmp2, tmp3, 128);

            STORE_4VECS_UB(dst_tmp, dst_stride, tmp0, tmp1, tmp2, tmp3);
            dst_tmp += (4 * dst_stride);

            src10_r = src54_r;
            src32_r = src76_r;
            src54_r = src98_r;
            src21_r = src65_r;
            src43_r = src87_r;
            src65_r = src109_r;
            src10_l = src54_l;
            src32_l = src76_l;
            src54_l = src98_l;
            src21_l = src65_l;
            src43_l = src87_l;
            src65_l = src109_l;

            src6 = src10;
        }

        src += 16;
        dst += 16;
    }
}

static void common_vt_8t_24w_msa(uint8_t *src, int32_t src_stride,
                                 uint8_t *dst, int32_t dst_stride,
                                 const int8_t *filter, int32_t height,
                                 uint8_t rnd_val)
{
    common_vt_8t_16w_mult_msa(src, src_stride, dst, dst_stride, filter,
                              height, rnd_val, 16);

    common_vt_8t_8w_msa(src + 16, src_stride, dst + 16, dst_stride,
                        filter, height, rnd_val);
}

static void common_vt_8t_32w_msa(uint8_t *src, int32_t src_stride,
                                 uint8_t *dst, int32_t dst_stride,
                                 const int8_t *filter, int32_t height,
                                 uint8_t rnd_val)
{
    common_vt_8t_16w_mult_msa(src, src_stride, dst, dst_stride,
                              filter, height, rnd_val, 32);
}

static void common_vt_8t_48w_msa(uint8_t *src, int32_t src_stride,
                                 uint8_t *dst, int32_t dst_stride,
                                 const int8_t *filter, int32_t height,
                                 uint8_t rnd_val)
{
    common_vt_8t_16w_mult_msa(src, src_stride, dst, dst_stride,
                              filter, height, rnd_val, 48);
}

static void common_vt_8t_64w_msa(uint8_t *src, int32_t src_stride,
                                 uint8_t *dst, int32_t dst_stride,
                                 const int8_t *filter, int32_t height,
                                 uint8_t rnd_val)
{
    common_vt_8t_16w_mult_msa(src, src_stride, dst, dst_stride,
                              filter, height, rnd_val, 64);
}

static void copy_width8_msa(uint8_t *src, int32_t src_stride,
                            uint8_t *dst, int32_t dst_stride,
                            int32_t height)
{
    int32_t cnt;
    uint64_t out0, out1, out2, out3, out4, out5, out6, out7;
    v16u8 src0, src1, src2, src3, src4, src5, src6, src7;

    if (0 == height % 12) {
        for (cnt = (height / 12); cnt--;) {
            LOAD_8VECS_UB(src, src_stride,
                          src0, src1, src2, src3, src4, src5, src6, src7);
            src += (8 * src_stride);

            out0 = __msa_copy_u_d((v2i64) src0, 0);
            out1 = __msa_copy_u_d((v2i64) src1, 0);
            out2 = __msa_copy_u_d((v2i64) src2, 0);
            out3 = __msa_copy_u_d((v2i64) src3, 0);
            out4 = __msa_copy_u_d((v2i64) src4, 0);
            out5 = __msa_copy_u_d((v2i64) src5, 0);
            out6 = __msa_copy_u_d((v2i64) src6, 0);
            out7 = __msa_copy_u_d((v2i64) src7, 0);

            STORE_DWORD(dst, out0);
            dst += dst_stride;
            STORE_DWORD(dst, out1);
            dst += dst_stride;
            STORE_DWORD(dst, out2);
            dst += dst_stride;
            STORE_DWORD(dst, out3);
            dst += dst_stride;
            STORE_DWORD(dst, out4);
            dst += dst_stride;
            STORE_DWORD(dst, out5);
            dst += dst_stride;
            STORE_DWORD(dst, out6);
            dst += dst_stride;
            STORE_DWORD(dst, out7);
            dst += dst_stride;

            LOAD_4VECS_UB(src, src_stride, src0, src1, src2, src3);
            src += (4 * src_stride);

            out0 = __msa_copy_u_d((v2i64) src0, 0);
            out1 = __msa_copy_u_d((v2i64) src1, 0);
            out2 = __msa_copy_u_d((v2i64) src2, 0);
            out3 = __msa_copy_u_d((v2i64) src3, 0);

            STORE_DWORD(dst, out0);
            dst += dst_stride;
            STORE_DWORD(dst, out1);
            dst += dst_stride;
            STORE_DWORD(dst, out2);
            dst += dst_stride;
            STORE_DWORD(dst, out3);
            dst += dst_stride;
        }
    } else if (0 == height % 8) {
        for (cnt = height >> 3; cnt--;) {
            LOAD_8VECS_UB(src, src_stride,
                          src0, src1, src2, src3, src4, src5, src6, src7);
            src += (8 * src_stride);

            out0 = __msa_copy_u_d((v2i64) src0, 0);
            out1 = __msa_copy_u_d((v2i64) src1, 0);
            out2 = __msa_copy_u_d((v2i64) src2, 0);
            out3 = __msa_copy_u_d((v2i64) src3, 0);
            out4 = __msa_copy_u_d((v2i64) src4, 0);
            out5 = __msa_copy_u_d((v2i64) src5, 0);
            out6 = __msa_copy_u_d((v2i64) src6, 0);
            out7 = __msa_copy_u_d((v2i64) src7, 0);

            STORE_DWORD(dst, out0);
            dst += dst_stride;
            STORE_DWORD(dst, out1);
            dst += dst_stride;
            STORE_DWORD(dst, out2);
            dst += dst_stride;
            STORE_DWORD(dst, out3);
            dst += dst_stride;
            STORE_DWORD(dst, out4);
            dst += dst_stride;
            STORE_DWORD(dst, out5);
            dst += dst_stride;
            STORE_DWORD(dst, out6);
            dst += dst_stride;
            STORE_DWORD(dst, out7);
            dst += dst_stride;
        }
    } else if (0 == height % 4) {
        for (cnt = (height / 4); cnt--;) {
            LOAD_4VECS_UB(src, src_stride, src0, src1, src2, src3);
            src += (4 * src_stride);

            out0 = __msa_copy_u_d((v2i64) src0, 0);
            out1 = __msa_copy_u_d((v2i64) src1, 0);
            out2 = __msa_copy_u_d((v2i64) src2, 0);
            out3 = __msa_copy_u_d((v2i64) src3, 0);

            STORE_DWORD(dst, out0);
            dst += dst_stride;
            STORE_DWORD(dst, out1);
            dst += dst_stride;
            STORE_DWORD(dst, out2);
            dst += dst_stride;
            STORE_DWORD(dst, out3);
            dst += dst_stride;
        }
    } else if (0 == height % 2) {
        for (cnt = (height / 2); cnt--;) {
            LOAD_2VECS_UB(src, src_stride, src0, src1);
            src += (2 * src_stride);

            out0 = __msa_copy_u_d((v2i64) src0, 0);
            out1 = __msa_copy_u_d((v2i64) src1, 0);

            STORE_DWORD(dst, out0);
            dst += dst_stride;
            STORE_DWORD(dst, out1);
            dst += dst_stride;
        }
    }
}

static void copy_width12_msa(uint8_t *src, int32_t src_stride,
                             uint8_t *dst, int32_t dst_stride,
                             int32_t height)
{
    int32_t cnt;
    uint64_t out0, out1, out2, out3, out4, out5, out6, out7;
    uint32_t out8, out9, out10, out11, out12, out13, out14, out15;
    v16u8 src0, src1, src2, src3, src4, src5, src6, src7;

    for (cnt = 2; cnt--;) {
        LOAD_8VECS_UB(src, src_stride,
                      src0, src1, src2, src3, src4, src5, src6, src7);
        src += (8 * src_stride);

        out0 = __msa_copy_u_d((v2i64) src0, 0);
        out1 = __msa_copy_u_d((v2i64) src1, 0);
        out2 = __msa_copy_u_d((v2i64) src2, 0);
        out3 = __msa_copy_u_d((v2i64) src3, 0);
        out4 = __msa_copy_u_d((v2i64) src4, 0);
        out5 = __msa_copy_u_d((v2i64) src5, 0);
        out6 = __msa_copy_u_d((v2i64) src6, 0);
        out7 = __msa_copy_u_d((v2i64) src7, 0);

        out8 = __msa_copy_u_w((v4i32) src0, 2);
        out9 = __msa_copy_u_w((v4i32) src1, 2);
        out10 = __msa_copy_u_w((v4i32) src2, 2);
        out11 = __msa_copy_u_w((v4i32) src3, 2);
        out12 = __msa_copy_u_w((v4i32) src4, 2);
        out13 = __msa_copy_u_w((v4i32) src5, 2);
        out14 = __msa_copy_u_w((v4i32) src6, 2);
        out15 = __msa_copy_u_w((v4i32) src7, 2);

        STORE_DWORD(dst, out0);
        STORE_WORD(dst + 8, out8);
        dst += dst_stride;
        STORE_DWORD(dst, out1);
        STORE_WORD(dst + 8, out9);
        dst += dst_stride;
        STORE_DWORD(dst, out2);
        STORE_WORD(dst + 8, out10);
        dst += dst_stride;
        STORE_DWORD(dst, out3);
        STORE_WORD(dst + 8, out11);
        dst += dst_stride;
        STORE_DWORD(dst, out4);
        STORE_WORD(dst + 8, out12);
        dst += dst_stride;
        STORE_DWORD(dst, out5);
        STORE_WORD(dst + 8, out13);
        dst += dst_stride;
        STORE_DWORD(dst, out6);
        STORE_WORD(dst + 8, out14);
        dst += dst_stride;
        STORE_DWORD(dst, out7);
        STORE_WORD(dst + 8, out15);
        dst += dst_stride;
    }
}

static void copy_16multx8mult_msa(uint8_t *src, int32_t src_stride,
                                  uint8_t *dst, int32_t dst_stride,
                                  int32_t height, int32_t width)
{
    int32_t cnt, loop_cnt;
    uint8_t *src_tmp, *dst_tmp;
    v16u8 src0, src1, src2, src3, src4, src5, src6, src7;

    for (cnt = (width >> 4); cnt--;) {
        src_tmp = src;
        dst_tmp = dst;

        for (loop_cnt = (height >> 3); loop_cnt--;) {
            LOAD_8VECS_UB(src_tmp, src_stride,
                          src0, src1, src2, src3, src4, src5, src6, src7);
            src_tmp += (8 * src_stride);

            STORE_8VECS_UB(dst_tmp, dst_stride,
                           src0, src1, src2, src3, src4, src5, src6, src7);
            dst_tmp += (8 * dst_stride);
        }

        src += 16;
        dst += 16;
    }
}

static void copy_width16_msa(uint8_t *src, int32_t src_stride,
                             uint8_t *dst, int32_t dst_stride,
                             int32_t height)
{
    int32_t cnt;
    v16u8 src0, src1, src2, src3, src4, src5, src6, src7;

    if (0 == height % 12) {
        for (cnt = (height / 12); cnt--;) {
            LOAD_8VECS_UB(src, src_stride,
                          src0, src1, src2, src3, src4, src5, src6, src7);
            src += (8 * src_stride);

            STORE_8VECS_UB(dst, dst_stride,
                           src0, src1, src2, src3, src4, src5, src6, src7);
            dst += (8 * dst_stride);

            LOAD_4VECS_UB(src, src_stride, src0, src1, src2, src3);
            src += (4 * src_stride);

            STORE_4VECS_UB(dst, dst_stride, src0, src1, src2, src3);
            dst += (4 * dst_stride);
        }
    } else if (0 == height % 8) {
        copy_16multx8mult_msa(src, src_stride, dst, dst_stride, height, 16);
    } else if (0 == height % 4) {
        for (cnt = (height >> 2); cnt--;) {
            LOAD_4VECS_UB(src, src_stride, src0, src1, src2, src3);
            src += (4 * src_stride);

            STORE_4VECS_UB(dst, dst_stride, src0, src1, src2, src3);
            dst += (4 * dst_stride);
        }
    }
}

static void copy_width24_msa(uint8_t *src, int32_t src_stride,
                             uint8_t *dst, int32_t dst_stride,
                             int32_t height)
{
    copy_16multx8mult_msa(src, src_stride, dst, dst_stride, height, 16);

    copy_width8_msa(src + 16, src_stride, dst + 16, dst_stride, height);
}

static void copy_width32_msa(uint8_t *src, int32_t src_stride,
                             uint8_t *dst, int32_t dst_stride,
                             int32_t height)
{
    int32_t cnt;
    v16u8 src0, src1, src2, src3, src4, src5, src6, src7;

    if (0 == height % 12) {
        for (cnt = (height / 12); cnt--;) {
            LOAD_4VECS_UB(src, src_stride, src0, src1, src2, src3);
            LOAD_4VECS_UB(src + 16, src_stride, src4, src5, src6, src7);
            src += (4 * src_stride);

            STORE_4VECS_UB(dst, dst_stride, src0, src1, src2, src3);
            STORE_4VECS_UB(dst + 16, dst_stride, src4, src5, src6, src7);
            dst += (4 * dst_stride);

            LOAD_4VECS_UB(src, src_stride, src0, src1, src2, src3);
            LOAD_4VECS_UB(src + 16, src_stride, src4, src5, src6, src7);
            src += (4 * src_stride);

            STORE_4VECS_UB(dst, dst_stride, src0, src1, src2, src3);
            STORE_4VECS_UB(dst + 16, dst_stride, src4, src5, src6, src7);
            dst += (4 * dst_stride);

            LOAD_4VECS_UB(src, src_stride, src0, src1, src2, src3);
            LOAD_4VECS_UB(src + 16, src_stride, src4, src5, src6, src7);
            src += (4 * src_stride);

            STORE_4VECS_UB(dst, dst_stride, src0, src1, src2, src3);
            STORE_4VECS_UB(dst + 16, dst_stride, src4, src5, src6, src7);
            dst += (4 * dst_stride);
        }
    } else if (0 == height % 8) {
        copy_16multx8mult_msa(src, src_stride, dst, dst_stride, height, 32);
    } else if (0 == height % 4) {
        for (cnt = (height >> 2); cnt--;) {
            LOAD_4VECS_UB(src, src_stride, src0, src1, src2, src3);
            LOAD_4VECS_UB(src + 16, src_stride, src4, src5, src6, src7);
            src += (4 * src_stride);

            STORE_4VECS_UB(dst, dst_stride, src0, src1, src2, src3);
            STORE_4VECS_UB(dst + 16, dst_stride, src4, src5, src6, src7);
            dst += (4 * dst_stride);
        }
    }
}

static void copy_width48_msa(uint8_t *src, int32_t src_stride,
                             uint8_t *dst, int32_t dst_stride,
                             int32_t height)
{
    copy_16multx8mult_msa(src, src_stride, dst, dst_stride, height, 48);
}

static void copy_width64_msa(uint8_t *src, int32_t src_stride,
                             uint8_t *dst, int32_t dst_stride,
                             int32_t height)
{
    copy_16multx8mult_msa(src, src_stride, dst, dst_stride, height, 64);
}

#define MC_COPY(WIDTH)                                                    \
void ff_hevc_put_hevc_pel_pixels##WIDTH##_8_msa(int16_t *dst,             \
                                                uint8_t *src,             \
                                                ptrdiff_t src_stride,     \
                                                int height,               \
                                                intptr_t mx,              \
                                                intptr_t my,              \
                                                int width)                \
{                                                                         \
    hevc_copy_##WIDTH##w_msa(src, src_stride, dst, MAX_PB_SIZE, height);  \
}

MC_COPY(4);
MC_COPY(6);
MC_COPY(8);
MC_COPY(12);
MC_COPY(16);
MC_COPY(24);
MC_COPY(32);
MC_COPY(48);
MC_COPY(64);

#undef MC_COPY

#define MC(PEL, DIR, WIDTH, TAP, DIR1, FILT_DIR)                            \
void ff_hevc_put_hevc_##PEL##_##DIR####WIDTH##_8_msa(int16_t *dst,          \
                                                     uint8_t *src,          \
                                                     ptrdiff_t src_stride,  \
                                                     int height,            \
                                                     intptr_t mx,           \
                                                     intptr_t my,           \
                                                     int width)             \
{                                                                           \
    const int8_t *filter = ff_hevc_##PEL##_filters[FILT_DIR - 1];           \
                                                                            \
    hevc_##DIR1##_##TAP##t_##WIDTH##w_msa(src, src_stride, dst,             \
                                          MAX_PB_SIZE, filter, height);     \
}

MC(qpel, h, 4, 8, hz, mx);
MC(qpel, h, 8, 8, hz, mx);
MC(qpel, h, 12, 8, hz, mx);
MC(qpel, h, 16, 8, hz, mx);
MC(qpel, h, 24, 8, hz, mx);
MC(qpel, h, 32, 8, hz, mx);
MC(qpel, h, 48, 8, hz, mx);
MC(qpel, h, 64, 8, hz, mx);

MC(qpel, v, 4, 8, vt, my);
MC(qpel, v, 8, 8, vt, my);
MC(qpel, v, 12, 8, vt, my);
MC(qpel, v, 16, 8, vt, my);
MC(qpel, v, 24, 8, vt, my);
MC(qpel, v, 32, 8, vt, my);
MC(qpel, v, 48, 8, vt, my);
MC(qpel, v, 64, 8, vt, my);

#undef MC

#define MC_HV(PEL, DIR, WIDTH, TAP, DIR1)                                     \
void ff_hevc_put_hevc_##PEL##_##DIR####WIDTH##_8_msa(int16_t *dst,            \
                                                     uint8_t *src,            \
                                                     ptrdiff_t src_stride,    \
                                                     int height,              \
                                                     intptr_t mx,             \
                                                     intptr_t my,             \
                                                     int width)               \
{                                                                             \
    const int8_t *filter_x = ff_hevc_##PEL##_filters[mx - 1];                 \
    const int8_t *filter_y = ff_hevc_##PEL##_filters[my - 1];                 \
                                                                              \
    hevc_##DIR1##_##TAP##t_##WIDTH##w_msa(src, src_stride, dst, MAX_PB_SIZE,  \
                                          filter_x, filter_y, height);        \
}

MC_HV(qpel, hv, 4, 8, hv);
MC_HV(qpel, hv, 8, 8, hv);
MC_HV(qpel, hv, 12, 8, hv);
MC_HV(qpel, hv, 16, 8, hv);
MC_HV(qpel, hv, 24, 8, hv);
MC_HV(qpel, hv, 32, 8, hv);
MC_HV(qpel, hv, 48, 8, hv);
MC_HV(qpel, hv, 64, 8, hv);

#undef MC_HV

#define UNI_MC_COPY(WIDTH)                                                 \
void ff_hevc_put_hevc_uni_pel_pixels##WIDTH##_8_msa(uint8_t *dst,          \
                                                    ptrdiff_t dst_stride,  \
                                                    uint8_t *src,          \
                                                    ptrdiff_t src_stride,  \
                                                    int height,            \
                                                    intptr_t mx,           \
                                                    intptr_t my,           \
                                                    int width)             \
{                                                                          \
    copy_width##WIDTH##_msa(src, src_stride, dst, dst_stride, height);     \
}

UNI_MC_COPY(8);
UNI_MC_COPY(12);
UNI_MC_COPY(16);
UNI_MC_COPY(24);
UNI_MC_COPY(32);
UNI_MC_COPY(48);
UNI_MC_COPY(64);

#undef UNI_MC_COPY

#define UNI_MC(PEL, DIR, WIDTH, TAP, DIR1, FILT_DIR)                           \
void ff_hevc_put_hevc_uni_##PEL##_##DIR####WIDTH##_8_msa(uint8_t *dst,         \
                                                         ptrdiff_t             \
                                                         dst_stride,           \
                                                         uint8_t *src,         \
                                                         ptrdiff_t             \
                                                         src_stride,           \
                                                         int height,           \
                                                         intptr_t mx,          \
                                                         intptr_t my,          \
                                                         int width)            \
{                                                                              \
    const int8_t *filter = ff_hevc_##PEL##_filters[FILT_DIR - 1];              \
                                                                               \
    common_##DIR1##_##TAP##t_##WIDTH##w_msa(src, src_stride, dst, dst_stride,  \
                                            filter, height, 6);                \
}

UNI_MC(qpel, h, 4, 8, hz, mx);
UNI_MC(qpel, h, 8, 8, hz, mx);
UNI_MC(qpel, h, 12, 8, hz, mx);
UNI_MC(qpel, h, 16, 8, hz, mx);
UNI_MC(qpel, h, 24, 8, hz, mx);
UNI_MC(qpel, h, 32, 8, hz, mx);
UNI_MC(qpel, h, 48, 8, hz, mx);
UNI_MC(qpel, h, 64, 8, hz, mx);

UNI_MC(qpel, v, 4, 8, vt, my);
UNI_MC(qpel, v, 8, 8, vt, my);
UNI_MC(qpel, v, 12, 8, vt, my);
UNI_MC(qpel, v, 16, 8, vt, my);
UNI_MC(qpel, v, 24, 8, vt, my);
UNI_MC(qpel, v, 32, 8, vt, my);
UNI_MC(qpel, v, 48, 8, vt, my);
UNI_MC(qpel, v, 64, 8, vt, my);

#undef UNI_MC