1 /*
2  *  Copyright (c) 2015 The WebM project authors. All Rights Reserved.
3  *
4  *  Use of this source code is governed by a BSD-style license
5  *  that can be found in the LICENSE file in the root of the source
6  *  tree. An additional intellectual property rights grant can be found
7  *  in the file PATENTS.  All contributing project authors may
8  *  be found in the AUTHORS file in the root of the source tree.
9  */
10 
11 #include <assert.h>
12 #include "./vpx_dsp_rtcd.h"
13 #include "vpx_dsp/mips/vpx_convolve_msa.h"
14 
common_vt_8t_4w_msa(const uint8_t * src,int32_t src_stride,uint8_t * dst,int32_t dst_stride,int8_t * filter,int32_t height)15 static void common_vt_8t_4w_msa(const uint8_t *src, int32_t src_stride,
16                                 uint8_t *dst, int32_t dst_stride,
17                                 int8_t *filter, int32_t height) {
18   uint32_t loop_cnt;
19   v16i8 src0, src1, src2, src3, src4, src5, src6, src7, src8, src9, src10;
20   v16i8 src10_r, src32_r, src54_r, src76_r, src98_r, src21_r, src43_r;
21   v16i8 src65_r, src87_r, src109_r, src2110, src4332, src6554, src8776;
22   v16i8 src10998, filt0, filt1, filt2, filt3;
23   v16u8 out;
24   v8i16 filt, out10, out32;
25 
26   src -= (3 * src_stride);
27 
28   filt = LD_SH(filter);
29   SPLATI_H4_SB(filt, 0, 1, 2, 3, filt0, filt1, filt2, filt3);
30 
31   LD_SB7(src, src_stride, src0, src1, src2, src3, src4, src5, src6);
32   src += (7 * src_stride);
33 
34   ILVR_B4_SB(src1, src0, src3, src2, src5, src4, src2, src1, src10_r, src32_r,
35              src54_r, src21_r);
36   ILVR_B2_SB(src4, src3, src6, src5, src43_r, src65_r);
37   ILVR_D3_SB(src21_r, src10_r, src43_r, src32_r, src65_r, src54_r, src2110,
38              src4332, src6554);
39   XORI_B3_128_SB(src2110, src4332, src6554);
40 
41   for (loop_cnt = (height >> 2); loop_cnt--;) {
42     LD_SB4(src, src_stride, src7, src8, src9, src10);
43     src += (4 * src_stride);
44 
45     ILVR_B4_SB(src7, src6, src8, src7, src9, src8, src10, src9, src76_r,
46                src87_r, src98_r, src109_r);
47     ILVR_D2_SB(src87_r, src76_r, src109_r, src98_r, src8776, src10998);
48     XORI_B2_128_SB(src8776, src10998);
49     out10 = FILT_8TAP_DPADD_S_H(src2110, src4332, src6554, src8776, filt0,
50                                 filt1, filt2, filt3);
51     out32 = FILT_8TAP_DPADD_S_H(src4332, src6554, src8776, src10998, filt0,
52                                 filt1, filt2, filt3);
53     SRARI_H2_SH(out10, out32, FILTER_BITS);
54     SAT_SH2_SH(out10, out32, 7);
55     out = PCKEV_XORI128_UB(out10, out32);
56     ST4x4_UB(out, out, 0, 1, 2, 3, dst, dst_stride);
57     dst += (4 * dst_stride);
58 
59     src2110 = src6554;
60     src4332 = src8776;
61     src6554 = src10998;
62     src6 = src10;
63   }
64 }
65 
common_vt_8t_8w_msa(const uint8_t * src,int32_t src_stride,uint8_t * dst,int32_t dst_stride,int8_t * filter,int32_t height)66 static void common_vt_8t_8w_msa(const uint8_t *src, int32_t src_stride,
67                                 uint8_t *dst, int32_t dst_stride,
68                                 int8_t *filter, int32_t height) {
69   uint32_t loop_cnt;
70   v16i8 src0, src1, src2, src3, src4, src5, src6, src7, src8, src9, src10;
71   v16i8 src10_r, src32_r, src54_r, src76_r, src98_r, src21_r, src43_r;
72   v16i8 src65_r, src87_r, src109_r, filt0, filt1, filt2, filt3;
73   v16u8 tmp0, tmp1;
74   v8i16 filt, out0_r, out1_r, out2_r, out3_r;
75 
76   src -= (3 * src_stride);
77 
78   filt = LD_SH(filter);
79   SPLATI_H4_SB(filt, 0, 1, 2, 3, filt0, filt1, filt2, filt3);
80 
81   LD_SB7(src, src_stride, src0, src1, src2, src3, src4, src5, src6);
82   XORI_B7_128_SB(src0, src1, src2, src3, src4, src5, src6);
83   src += (7 * src_stride);
84   ILVR_B4_SB(src1, src0, src3, src2, src5, src4, src2, src1, src10_r, src32_r,
85              src54_r, src21_r);
86   ILVR_B2_SB(src4, src3, src6, src5, src43_r, src65_r);
87 
88   for (loop_cnt = (height >> 2); loop_cnt--;) {
89     LD_SB4(src, src_stride, src7, src8, src9, src10);
90     XORI_B4_128_SB(src7, src8, src9, src10);
91     src += (4 * src_stride);
92 
93     ILVR_B4_SB(src7, src6, src8, src7, src9, src8, src10, src9, src76_r,
94                src87_r, src98_r, src109_r);
95     out0_r = FILT_8TAP_DPADD_S_H(src10_r, src32_r, src54_r, src76_r, filt0,
96                                  filt1, filt2, filt3);
97     out1_r = FILT_8TAP_DPADD_S_H(src21_r, src43_r, src65_r, src87_r, filt0,
98                                  filt1, filt2, filt3);
99     out2_r = FILT_8TAP_DPADD_S_H(src32_r, src54_r, src76_r, src98_r, filt0,
100                                  filt1, filt2, filt3);
101     out3_r = FILT_8TAP_DPADD_S_H(src43_r, src65_r, src87_r, src109_r, filt0,
102                                  filt1, filt2, filt3);
103     SRARI_H4_SH(out0_r, out1_r, out2_r, out3_r, FILTER_BITS);
104     SAT_SH4_SH(out0_r, out1_r, out2_r, out3_r, 7);
105     tmp0 = PCKEV_XORI128_UB(out0_r, out1_r);
106     tmp1 = PCKEV_XORI128_UB(out2_r, out3_r);
107     ST8x4_UB(tmp0, tmp1, dst, dst_stride);
108     dst += (4 * dst_stride);
109 
110     src10_r = src54_r;
111     src32_r = src76_r;
112     src54_r = src98_r;
113     src21_r = src65_r;
114     src43_r = src87_r;
115     src65_r = src109_r;
116     src6 = src10;
117   }
118 }
119 
common_vt_8t_16w_msa(const uint8_t * src,int32_t src_stride,uint8_t * dst,int32_t dst_stride,int8_t * filter,int32_t height)120 static void common_vt_8t_16w_msa(const uint8_t *src, int32_t src_stride,
121                                  uint8_t *dst, int32_t dst_stride,
122                                  int8_t *filter, int32_t height) {
123   uint32_t loop_cnt;
124   v16i8 src0, src1, src2, src3, src4, src5, src6, src7, src8, src9, src10;
125   v16i8 filt0, filt1, filt2, filt3;
126   v16i8 src10_r, src32_r, src54_r, src76_r, src98_r, src21_r, src43_r;
127   v16i8 src65_r, src87_r, src109_r, src10_l, src32_l, src54_l, src76_l;
128   v16i8 src98_l, src21_l, src43_l, src65_l, src87_l, src109_l;
129   v16u8 tmp0, tmp1, tmp2, tmp3;
130   v8i16 filt, out0_r, out1_r, out2_r, out3_r, out0_l, out1_l, out2_l, out3_l;
131 
132   src -= (3 * src_stride);
133 
134   filt = LD_SH(filter);
135   SPLATI_H4_SB(filt, 0, 1, 2, 3, filt0, filt1, filt2, filt3);
136 
137   LD_SB7(src, src_stride, src0, src1, src2, src3, src4, src5, src6);
138   XORI_B7_128_SB(src0, src1, src2, src3, src4, src5, src6);
139   src += (7 * src_stride);
140   ILVR_B4_SB(src1, src0, src3, src2, src5, src4, src2, src1, src10_r, src32_r,
141              src54_r, src21_r);
142   ILVR_B2_SB(src4, src3, src6, src5, src43_r, src65_r);
143   ILVL_B4_SB(src1, src0, src3, src2, src5, src4, src2, src1, src10_l, src32_l,
144              src54_l, src21_l);
145   ILVL_B2_SB(src4, src3, src6, src5, src43_l, src65_l);
146 
147   for (loop_cnt = (height >> 2); loop_cnt--;) {
148     LD_SB4(src, src_stride, src7, src8, src9, src10);
149     XORI_B4_128_SB(src7, src8, src9, src10);
150     src += (4 * src_stride);
151 
152     ILVR_B4_SB(src7, src6, src8, src7, src9, src8, src10, src9, src76_r,
153                src87_r, src98_r, src109_r);
154     ILVL_B4_SB(src7, src6, src8, src7, src9, src8, src10, src9, src76_l,
155                src87_l, src98_l, src109_l);
156     out0_r = FILT_8TAP_DPADD_S_H(src10_r, src32_r, src54_r, src76_r, filt0,
157                                  filt1, filt2, filt3);
158     out1_r = FILT_8TAP_DPADD_S_H(src21_r, src43_r, src65_r, src87_r, filt0,
159                                  filt1, filt2, filt3);
160     out2_r = FILT_8TAP_DPADD_S_H(src32_r, src54_r, src76_r, src98_r, filt0,
161                                  filt1, filt2, filt3);
162     out3_r = FILT_8TAP_DPADD_S_H(src43_r, src65_r, src87_r, src109_r, filt0,
163                                  filt1, filt2, filt3);
164     out0_l = FILT_8TAP_DPADD_S_H(src10_l, src32_l, src54_l, src76_l, filt0,
165                                  filt1, filt2, filt3);
166     out1_l = FILT_8TAP_DPADD_S_H(src21_l, src43_l, src65_l, src87_l, filt0,
167                                  filt1, filt2, filt3);
168     out2_l = FILT_8TAP_DPADD_S_H(src32_l, src54_l, src76_l, src98_l, filt0,
169                                  filt1, filt2, filt3);
170     out3_l = FILT_8TAP_DPADD_S_H(src43_l, src65_l, src87_l, src109_l, filt0,
171                                  filt1, filt2, filt3);
172     SRARI_H4_SH(out0_r, out1_r, out2_r, out3_r, FILTER_BITS);
173     SRARI_H4_SH(out0_l, out1_l, out2_l, out3_l, FILTER_BITS);
174     SAT_SH4_SH(out0_r, out1_r, out2_r, out3_r, 7);
175     SAT_SH4_SH(out0_l, out1_l, out2_l, out3_l, 7);
176     PCKEV_B4_UB(out0_l, out0_r, out1_l, out1_r, out2_l, out2_r, out3_l, out3_r,
177                 tmp0, tmp1, tmp2, tmp3);
178     XORI_B4_128_UB(tmp0, tmp1, tmp2, tmp3);
179     ST_UB4(tmp0, tmp1, tmp2, tmp3, dst, dst_stride);
180     dst += (4 * dst_stride);
181 
182     src10_r = src54_r;
183     src32_r = src76_r;
184     src54_r = src98_r;
185     src21_r = src65_r;
186     src43_r = src87_r;
187     src65_r = src109_r;
188     src10_l = src54_l;
189     src32_l = src76_l;
190     src54_l = src98_l;
191     src21_l = src65_l;
192     src43_l = src87_l;
193     src65_l = src109_l;
194     src6 = src10;
195   }
196 }
197 
common_vt_8t_16w_mult_msa(const uint8_t * src,int32_t src_stride,uint8_t * dst,int32_t dst_stride,int8_t * filter,int32_t height,int32_t width)198 static void common_vt_8t_16w_mult_msa(const uint8_t *src, int32_t src_stride,
199                                       uint8_t *dst, int32_t dst_stride,
200                                       int8_t *filter, int32_t height,
201                                       int32_t width) {
202   const uint8_t *src_tmp;
203   uint8_t *dst_tmp;
204   uint32_t loop_cnt, cnt;
205   v16i8 src0, src1, src2, src3, src4, src5, src6, src7, src8, src9, src10;
206   v16i8 filt0, filt1, filt2, filt3;
207   v16i8 src10_r, src32_r, src54_r, src76_r, src98_r, src21_r, src43_r;
208   v16i8 src65_r, src87_r, src109_r, src10_l, src32_l, src54_l, src76_l;
209   v16i8 src98_l, src21_l, src43_l, src65_l, src87_l, src109_l;
210   v16u8 tmp0, tmp1, tmp2, tmp3;
211   v8i16 filt, out0_r, out1_r, out2_r, out3_r, out0_l, out1_l, out2_l, out3_l;
212 
213   src -= (3 * src_stride);
214 
215   filt = LD_SH(filter);
216   SPLATI_H4_SB(filt, 0, 1, 2, 3, filt0, filt1, filt2, filt3);
217 
218   for (cnt = (width >> 4); cnt--;) {
219     src_tmp = src;
220     dst_tmp = dst;
221 
222     LD_SB7(src_tmp, src_stride, src0, src1, src2, src3, src4, src5, src6);
223     XORI_B7_128_SB(src0, src1, src2, src3, src4, src5, src6);
224     src_tmp += (7 * src_stride);
225     ILVR_B4_SB(src1, src0, src3, src2, src5, src4, src2, src1, src10_r, src32_r,
226                src54_r, src21_r);
227     ILVR_B2_SB(src4, src3, src6, src5, src43_r, src65_r);
228     ILVL_B4_SB(src1, src0, src3, src2, src5, src4, src2, src1, src10_l, src32_l,
229                src54_l, src21_l);
230     ILVL_B2_SB(src4, src3, src6, src5, src43_l, src65_l);
231 
232     for (loop_cnt = (height >> 2); loop_cnt--;) {
233       LD_SB4(src_tmp, src_stride, src7, src8, src9, src10);
234       XORI_B4_128_SB(src7, src8, src9, src10);
235       src_tmp += (4 * src_stride);
236       ILVR_B4_SB(src7, src6, src8, src7, src9, src8, src10, src9, src76_r,
237                  src87_r, src98_r, src109_r);
238       ILVL_B4_SB(src7, src6, src8, src7, src9, src8, src10, src9, src76_l,
239                  src87_l, src98_l, src109_l);
240       out0_r = FILT_8TAP_DPADD_S_H(src10_r, src32_r, src54_r, src76_r, filt0,
241                                    filt1, filt2, filt3);
242       out1_r = FILT_8TAP_DPADD_S_H(src21_r, src43_r, src65_r, src87_r, filt0,
243                                    filt1, filt2, filt3);
244       out2_r = FILT_8TAP_DPADD_S_H(src32_r, src54_r, src76_r, src98_r, filt0,
245                                    filt1, filt2, filt3);
246       out3_r = FILT_8TAP_DPADD_S_H(src43_r, src65_r, src87_r, src109_r, filt0,
247                                    filt1, filt2, filt3);
248       out0_l = FILT_8TAP_DPADD_S_H(src10_l, src32_l, src54_l, src76_l, filt0,
249                                    filt1, filt2, filt3);
250       out1_l = FILT_8TAP_DPADD_S_H(src21_l, src43_l, src65_l, src87_l, filt0,
251                                    filt1, filt2, filt3);
252       out2_l = FILT_8TAP_DPADD_S_H(src32_l, src54_l, src76_l, src98_l, filt0,
253                                    filt1, filt2, filt3);
254       out3_l = FILT_8TAP_DPADD_S_H(src43_l, src65_l, src87_l, src109_l, filt0,
255                                    filt1, filt2, filt3);
256       SRARI_H4_SH(out0_r, out1_r, out2_r, out3_r, FILTER_BITS);
257       SRARI_H4_SH(out0_l, out1_l, out2_l, out3_l, FILTER_BITS);
258       SAT_SH4_SH(out0_r, out1_r, out2_r, out3_r, 7);
259       SAT_SH4_SH(out0_l, out1_l, out2_l, out3_l, 7);
260       PCKEV_B4_UB(out0_l, out0_r, out1_l, out1_r, out2_l, out2_r, out3_l,
261                   out3_r, tmp0, tmp1, tmp2, tmp3);
262       XORI_B4_128_UB(tmp0, tmp1, tmp2, tmp3);
263       ST_UB4(tmp0, tmp1, tmp2, tmp3, dst_tmp, dst_stride);
264       dst_tmp += (4 * dst_stride);
265 
266       src10_r = src54_r;
267       src32_r = src76_r;
268       src54_r = src98_r;
269       src21_r = src65_r;
270       src43_r = src87_r;
271       src65_r = src109_r;
272       src10_l = src54_l;
273       src32_l = src76_l;
274       src54_l = src98_l;
275       src21_l = src65_l;
276       src43_l = src87_l;
277       src65_l = src109_l;
278       src6 = src10;
279     }
280 
281     src += 16;
282     dst += 16;
283   }
284 }
285 
common_vt_8t_32w_msa(const uint8_t * src,int32_t src_stride,uint8_t * dst,int32_t dst_stride,int8_t * filter,int32_t height)286 static void common_vt_8t_32w_msa(const uint8_t *src, int32_t src_stride,
287                                  uint8_t *dst, int32_t dst_stride,
288                                  int8_t *filter, int32_t height) {
289   common_vt_8t_16w_mult_msa(src, src_stride, dst, dst_stride, filter, height,
290                             32);
291 }
292 
common_vt_8t_64w_msa(const uint8_t * src,int32_t src_stride,uint8_t * dst,int32_t dst_stride,int8_t * filter,int32_t height)293 static void common_vt_8t_64w_msa(const uint8_t *src, int32_t src_stride,
294                                  uint8_t *dst, int32_t dst_stride,
295                                  int8_t *filter, int32_t height) {
296   common_vt_8t_16w_mult_msa(src, src_stride, dst, dst_stride, filter, height,
297                             64);
298 }
299 
common_vt_2t_4x4_msa(const uint8_t * src,int32_t src_stride,uint8_t * dst,int32_t dst_stride,int8_t * filter)300 static void common_vt_2t_4x4_msa(const uint8_t *src, int32_t src_stride,
301                                  uint8_t *dst, int32_t dst_stride,
302                                  int8_t *filter) {
303   v16i8 src0, src1, src2, src3, src4;
304   v16i8 src10_r, src32_r, src21_r, src43_r, src2110, src4332;
305   v16u8 filt0;
306   v8i16 filt;
307   v8u16 tmp0, tmp1;
308 
309   filt = LD_SH(filter);
310   filt0 = (v16u8)__msa_splati_h(filt, 0);
311 
312   LD_SB5(src, src_stride, src0, src1, src2, src3, src4);
313   src += (5 * src_stride);
314 
315   ILVR_B4_SB(src1, src0, src2, src1, src3, src2, src4, src3, src10_r, src21_r,
316              src32_r, src43_r);
317   ILVR_D2_SB(src21_r, src10_r, src43_r, src32_r, src2110, src4332);
318   DOTP_UB2_UH(src2110, src4332, filt0, filt0, tmp0, tmp1);
319   SRARI_H2_UH(tmp0, tmp1, FILTER_BITS);
320   src2110 = __msa_pckev_b((v16i8)tmp1, (v16i8)tmp0);
321   ST4x4_UB(src2110, src2110, 0, 1, 2, 3, dst, dst_stride);
322 }
323 
common_vt_2t_4x8_msa(const uint8_t * src,int32_t src_stride,uint8_t * dst,int32_t dst_stride,int8_t * filter)324 static void common_vt_2t_4x8_msa(const uint8_t *src, int32_t src_stride,
325                                  uint8_t *dst, int32_t dst_stride,
326                                  int8_t *filter) {
327   v16i8 src0, src1, src2, src3, src4, src5, src6, src7, src8;
328   v16i8 src10_r, src32_r, src54_r, src76_r, src21_r, src43_r;
329   v16i8 src65_r, src87_r, src2110, src4332, src6554, src8776;
330   v8u16 tmp0, tmp1, tmp2, tmp3;
331   v16u8 filt0;
332   v8i16 filt;
333 
334   filt = LD_SH(filter);
335   filt0 = (v16u8)__msa_splati_h(filt, 0);
336 
337   LD_SB8(src, src_stride, src0, src1, src2, src3, src4, src5, src6, src7);
338   src += (8 * src_stride);
339 
340   src8 = LD_SB(src);
341   src += src_stride;
342 
343   ILVR_B4_SB(src1, src0, src2, src1, src3, src2, src4, src3, src10_r, src21_r,
344              src32_r, src43_r);
345   ILVR_B4_SB(src5, src4, src6, src5, src7, src6, src8, src7, src54_r, src65_r,
346              src76_r, src87_r);
347   ILVR_D4_SB(src21_r, src10_r, src43_r, src32_r, src65_r, src54_r, src87_r,
348              src76_r, src2110, src4332, src6554, src8776);
349   DOTP_UB4_UH(src2110, src4332, src6554, src8776, filt0, filt0, filt0, filt0,
350               tmp0, tmp1, tmp2, tmp3);
351   SRARI_H4_UH(tmp0, tmp1, tmp2, tmp3, FILTER_BITS);
352   PCKEV_B2_SB(tmp1, tmp0, tmp3, tmp2, src2110, src4332);
353   ST4x4_UB(src2110, src2110, 0, 1, 2, 3, dst, dst_stride);
354   ST4x4_UB(src4332, src4332, 0, 1, 2, 3, dst + 4 * dst_stride, dst_stride);
355 }
356 
common_vt_2t_4w_msa(const uint8_t * src,int32_t src_stride,uint8_t * dst,int32_t dst_stride,int8_t * filter,int32_t height)357 static void common_vt_2t_4w_msa(const uint8_t *src, int32_t src_stride,
358                                 uint8_t *dst, int32_t dst_stride,
359                                 int8_t *filter, int32_t height) {
360   if (4 == height) {
361     common_vt_2t_4x4_msa(src, src_stride, dst, dst_stride, filter);
362   } else if (8 == height) {
363     common_vt_2t_4x8_msa(src, src_stride, dst, dst_stride, filter);
364   }
365 }
366 
common_vt_2t_8x4_msa(const uint8_t * src,int32_t src_stride,uint8_t * dst,int32_t dst_stride,int8_t * filter)367 static void common_vt_2t_8x4_msa(const uint8_t *src, int32_t src_stride,
368                                  uint8_t *dst, int32_t dst_stride,
369                                  int8_t *filter) {
370   v16u8 src0, src1, src2, src3, src4, vec0, vec1, vec2, vec3, filt0;
371   v16i8 out0, out1;
372   v8u16 tmp0, tmp1, tmp2, tmp3;
373   v8i16 filt;
374 
375   /* rearranging filter_y */
376   filt = LD_SH(filter);
377   filt0 = (v16u8)__msa_splati_h(filt, 0);
378 
379   LD_UB5(src, src_stride, src0, src1, src2, src3, src4);
380   ILVR_B2_UB(src1, src0, src2, src1, vec0, vec1);
381   ILVR_B2_UB(src3, src2, src4, src3, vec2, vec3);
382   DOTP_UB4_UH(vec0, vec1, vec2, vec3, filt0, filt0, filt0, filt0, tmp0, tmp1,
383               tmp2, tmp3);
384   SRARI_H4_UH(tmp0, tmp1, tmp2, tmp3, FILTER_BITS);
385   PCKEV_B2_SB(tmp1, tmp0, tmp3, tmp2, out0, out1);
386   ST8x4_UB(out0, out1, dst, dst_stride);
387 }
388 
common_vt_2t_8x8mult_msa(const uint8_t * src,int32_t src_stride,uint8_t * dst,int32_t dst_stride,int8_t * filter,int32_t height)389 static void common_vt_2t_8x8mult_msa(const uint8_t *src, int32_t src_stride,
390                                      uint8_t *dst, int32_t dst_stride,
391                                      int8_t *filter, int32_t height) {
392   uint32_t loop_cnt;
393   v16u8 src0, src1, src2, src3, src4, src5, src6, src7, src8;
394   v16u8 vec0, vec1, vec2, vec3, vec4, vec5, vec6, vec7, filt0;
395   v16i8 out0, out1;
396   v8u16 tmp0, tmp1, tmp2, tmp3;
397   v8i16 filt;
398 
399   /* rearranging filter_y */
400   filt = LD_SH(filter);
401   filt0 = (v16u8)__msa_splati_h(filt, 0);
402 
403   src0 = LD_UB(src);
404   src += src_stride;
405 
406   for (loop_cnt = (height >> 3); loop_cnt--;) {
407     LD_UB8(src, src_stride, src1, src2, src3, src4, src5, src6, src7, src8);
408     src += (8 * src_stride);
409 
410     ILVR_B4_UB(src1, src0, src2, src1, src3, src2, src4, src3, vec0, vec1, vec2,
411                vec3);
412     ILVR_B4_UB(src5, src4, src6, src5, src7, src6, src8, src7, vec4, vec5, vec6,
413                vec7);
414     DOTP_UB4_UH(vec0, vec1, vec2, vec3, filt0, filt0, filt0, filt0, tmp0, tmp1,
415                 tmp2, tmp3);
416     SRARI_H4_UH(tmp0, tmp1, tmp2, tmp3, FILTER_BITS);
417     PCKEV_B2_SB(tmp1, tmp0, tmp3, tmp2, out0, out1);
418     ST8x4_UB(out0, out1, dst, dst_stride);
419     dst += (4 * dst_stride);
420 
421     DOTP_UB4_UH(vec4, vec5, vec6, vec7, filt0, filt0, filt0, filt0, tmp0, tmp1,
422                 tmp2, tmp3);
423     SRARI_H4_UH(tmp0, tmp1, tmp2, tmp3, FILTER_BITS);
424     PCKEV_B2_SB(tmp1, tmp0, tmp3, tmp2, out0, out1);
425     ST8x4_UB(out0, out1, dst, dst_stride);
426     dst += (4 * dst_stride);
427 
428     src0 = src8;
429   }
430 }
431 
common_vt_2t_8w_msa(const uint8_t * src,int32_t src_stride,uint8_t * dst,int32_t dst_stride,int8_t * filter,int32_t height)432 static void common_vt_2t_8w_msa(const uint8_t *src, int32_t src_stride,
433                                 uint8_t *dst, int32_t dst_stride,
434                                 int8_t *filter, int32_t height) {
435   if (4 == height) {
436     common_vt_2t_8x4_msa(src, src_stride, dst, dst_stride, filter);
437   } else {
438     common_vt_2t_8x8mult_msa(src, src_stride, dst, dst_stride, filter, height);
439   }
440 }
441 
common_vt_2t_16w_msa(const uint8_t * src,int32_t src_stride,uint8_t * dst,int32_t dst_stride,int8_t * filter,int32_t height)442 static void common_vt_2t_16w_msa(const uint8_t *src, int32_t src_stride,
443                                  uint8_t *dst, int32_t dst_stride,
444                                  int8_t *filter, int32_t height) {
445   uint32_t loop_cnt;
446   v16u8 src0, src1, src2, src3, src4;
447   v16u8 vec0, vec1, vec2, vec3, vec4, vec5, vec6, vec7, filt0;
448   v8u16 tmp0, tmp1, tmp2, tmp3;
449   v8i16 filt;
450 
451   /* rearranging filter_y */
452   filt = LD_SH(filter);
453   filt0 = (v16u8)__msa_splati_h(filt, 0);
454 
455   src0 = LD_UB(src);
456   src += src_stride;
457 
458   for (loop_cnt = (height >> 2); loop_cnt--;) {
459     LD_UB4(src, src_stride, src1, src2, src3, src4);
460     src += (4 * src_stride);
461 
462     ILVR_B2_UB(src1, src0, src2, src1, vec0, vec2);
463     ILVL_B2_UB(src1, src0, src2, src1, vec1, vec3);
464     DOTP_UB2_UH(vec0, vec1, filt0, filt0, tmp0, tmp1);
465     SRARI_H2_UH(tmp0, tmp1, FILTER_BITS);
466     PCKEV_ST_SB(tmp0, tmp1, dst);
467     dst += dst_stride;
468 
469     ILVR_B2_UB(src3, src2, src4, src3, vec4, vec6);
470     ILVL_B2_UB(src3, src2, src4, src3, vec5, vec7);
471     DOTP_UB2_UH(vec2, vec3, filt0, filt0, tmp2, tmp3);
472     SRARI_H2_UH(tmp2, tmp3, FILTER_BITS);
473     PCKEV_ST_SB(tmp2, tmp3, dst);
474     dst += dst_stride;
475 
476     DOTP_UB2_UH(vec4, vec5, filt0, filt0, tmp0, tmp1);
477     SRARI_H2_UH(tmp0, tmp1, FILTER_BITS);
478     PCKEV_ST_SB(tmp0, tmp1, dst);
479     dst += dst_stride;
480 
481     DOTP_UB2_UH(vec6, vec7, filt0, filt0, tmp2, tmp3);
482     SRARI_H2_UH(tmp2, tmp3, FILTER_BITS);
483     PCKEV_ST_SB(tmp2, tmp3, dst);
484     dst += dst_stride;
485 
486     src0 = src4;
487   }
488 }
489 
common_vt_2t_32w_msa(const uint8_t * src,int32_t src_stride,uint8_t * dst,int32_t dst_stride,int8_t * filter,int32_t height)490 static void common_vt_2t_32w_msa(const uint8_t *src, int32_t src_stride,
491                                  uint8_t *dst, int32_t dst_stride,
492                                  int8_t *filter, int32_t height) {
493   uint32_t loop_cnt;
494   v16u8 src0, src1, src2, src3, src4, src5, src6, src7, src8, src9;
495   v16u8 vec0, vec1, vec2, vec3, vec4, vec5, vec6, vec7, filt0;
496   v8u16 tmp0, tmp1, tmp2, tmp3;
497   v8i16 filt;
498 
499   /* rearranging filter_y */
500   filt = LD_SH(filter);
501   filt0 = (v16u8)__msa_splati_h(filt, 0);
502 
503   src0 = LD_UB(src);
504   src5 = LD_UB(src + 16);
505   src += src_stride;
506 
507   for (loop_cnt = (height >> 2); loop_cnt--;) {
508     LD_UB4(src, src_stride, src1, src2, src3, src4);
509     ILVR_B2_UB(src1, src0, src2, src1, vec0, vec2);
510     ILVL_B2_UB(src1, src0, src2, src1, vec1, vec3);
511 
512     LD_UB4(src + 16, src_stride, src6, src7, src8, src9);
513     src += (4 * src_stride);
514 
515     DOTP_UB2_UH(vec0, vec1, filt0, filt0, tmp0, tmp1);
516     SRARI_H2_UH(tmp0, tmp1, FILTER_BITS);
517     PCKEV_ST_SB(tmp0, tmp1, dst);
518     DOTP_UB2_UH(vec2, vec3, filt0, filt0, tmp2, tmp3);
519     SRARI_H2_UH(tmp2, tmp3, FILTER_BITS);
520     PCKEV_ST_SB(tmp2, tmp3, dst + dst_stride);
521 
522     ILVR_B2_UB(src3, src2, src4, src3, vec4, vec6);
523     ILVL_B2_UB(src3, src2, src4, src3, vec5, vec7);
524     DOTP_UB2_UH(vec4, vec5, filt0, filt0, tmp0, tmp1);
525     SRARI_H2_UH(tmp0, tmp1, FILTER_BITS);
526     PCKEV_ST_SB(tmp0, tmp1, dst + 2 * dst_stride);
527 
528     DOTP_UB2_UH(vec6, vec7, filt0, filt0, tmp2, tmp3);
529     SRARI_H2_UH(tmp2, tmp3, FILTER_BITS);
530     PCKEV_ST_SB(tmp2, tmp3, dst + 3 * dst_stride);
531 
532     ILVR_B2_UB(src6, src5, src7, src6, vec0, vec2);
533     ILVL_B2_UB(src6, src5, src7, src6, vec1, vec3);
534     DOTP_UB2_UH(vec0, vec1, filt0, filt0, tmp0, tmp1);
535     SRARI_H2_UH(tmp0, tmp1, FILTER_BITS);
536     PCKEV_ST_SB(tmp0, tmp1, dst + 16);
537 
538     DOTP_UB2_UH(vec2, vec3, filt0, filt0, tmp2, tmp3);
539     SRARI_H2_UH(tmp2, tmp3, FILTER_BITS);
540     PCKEV_ST_SB(tmp2, tmp3, dst + 16 + dst_stride);
541 
542     ILVR_B2_UB(src8, src7, src9, src8, vec4, vec6);
543     ILVL_B2_UB(src8, src7, src9, src8, vec5, vec7);
544     DOTP_UB2_UH(vec4, vec5, filt0, filt0, tmp0, tmp1);
545     SRARI_H2_UH(tmp0, tmp1, FILTER_BITS);
546     PCKEV_ST_SB(tmp0, tmp1, dst + 16 + 2 * dst_stride);
547 
548     DOTP_UB2_UH(vec6, vec7, filt0, filt0, tmp2, tmp3);
549     SRARI_H2_UH(tmp2, tmp3, FILTER_BITS);
550     PCKEV_ST_SB(tmp2, tmp3, dst + 16 + 3 * dst_stride);
551     dst += (4 * dst_stride);
552 
553     src0 = src4;
554     src5 = src9;
555   }
556 }
557 
common_vt_2t_64w_msa(const uint8_t * src,int32_t src_stride,uint8_t * dst,int32_t dst_stride,int8_t * filter,int32_t height)558 static void common_vt_2t_64w_msa(const uint8_t *src, int32_t src_stride,
559                                  uint8_t *dst, int32_t dst_stride,
560                                  int8_t *filter, int32_t height) {
561   uint32_t loop_cnt;
562   v16u8 src0, src1, src2, src3, src4, src5, src6, src7, src8, src9, src10;
563   v16u8 src11, vec0, vec1, vec2, vec3, vec4, vec5, vec6, vec7, filt0;
564   v8u16 tmp0, tmp1, tmp2, tmp3, tmp4, tmp5, tmp6, tmp7;
565   v8i16 filt;
566 
567   /* rearranging filter_y */
568   filt = LD_SH(filter);
569   filt0 = (v16u8)__msa_splati_h(filt, 0);
570 
571   LD_UB4(src, 16, src0, src3, src6, src9);
572   src += src_stride;
573 
574   for (loop_cnt = (height >> 1); loop_cnt--;) {
575     LD_UB2(src, src_stride, src1, src2);
576     LD_UB2(src + 16, src_stride, src4, src5);
577     LD_UB2(src + 32, src_stride, src7, src8);
578     LD_UB2(src + 48, src_stride, src10, src11);
579     src += (2 * src_stride);
580 
581     ILVR_B2_UB(src1, src0, src2, src1, vec0, vec2);
582     ILVL_B2_UB(src1, src0, src2, src1, vec1, vec3);
583     DOTP_UB2_UH(vec0, vec1, filt0, filt0, tmp0, tmp1);
584     SRARI_H2_UH(tmp0, tmp1, FILTER_BITS);
585     PCKEV_ST_SB(tmp0, tmp1, dst);
586 
587     DOTP_UB2_UH(vec2, vec3, filt0, filt0, tmp2, tmp3);
588     SRARI_H2_UH(tmp2, tmp3, FILTER_BITS);
589     PCKEV_ST_SB(tmp2, tmp3, dst + dst_stride);
590 
591     ILVR_B2_UB(src4, src3, src5, src4, vec4, vec6);
592     ILVL_B2_UB(src4, src3, src5, src4, vec5, vec7);
593     DOTP_UB2_UH(vec4, vec5, filt0, filt0, tmp4, tmp5);
594     SRARI_H2_UH(tmp4, tmp5, FILTER_BITS);
595     PCKEV_ST_SB(tmp4, tmp5, dst + 16);
596 
597     DOTP_UB2_UH(vec6, vec7, filt0, filt0, tmp6, tmp7);
598     SRARI_H2_UH(tmp6, tmp7, FILTER_BITS);
599     PCKEV_ST_SB(tmp6, tmp7, dst + 16 + dst_stride);
600 
601     ILVR_B2_UB(src7, src6, src8, src7, vec0, vec2);
602     ILVL_B2_UB(src7, src6, src8, src7, vec1, vec3);
603     DOTP_UB2_UH(vec0, vec1, filt0, filt0, tmp0, tmp1);
604     SRARI_H2_UH(tmp0, tmp1, FILTER_BITS);
605     PCKEV_ST_SB(tmp0, tmp1, dst + 32);
606 
607     DOTP_UB2_UH(vec2, vec3, filt0, filt0, tmp2, tmp3);
608     SRARI_H2_UH(tmp2, tmp3, FILTER_BITS);
609     PCKEV_ST_SB(tmp2, tmp3, dst + 32 + dst_stride);
610 
611     ILVR_B2_UB(src10, src9, src11, src10, vec4, vec6);
612     ILVL_B2_UB(src10, src9, src11, src10, vec5, vec7);
613     DOTP_UB2_UH(vec4, vec5, filt0, filt0, tmp4, tmp5);
614     SRARI_H2_UH(tmp4, tmp5, FILTER_BITS);
615     PCKEV_ST_SB(tmp4, tmp5, dst + 48);
616 
617     DOTP_UB2_UH(vec6, vec7, filt0, filt0, tmp6, tmp7);
618     SRARI_H2_UH(tmp6, tmp7, FILTER_BITS);
619     PCKEV_ST_SB(tmp6, tmp7, dst + 48 + dst_stride);
620     dst += (2 * dst_stride);
621 
622     src0 = src2;
623     src3 = src5;
624     src6 = src8;
625     src9 = src11;
626   }
627 }
628 
vpx_convolve8_vert_msa(const uint8_t * src,ptrdiff_t src_stride,uint8_t * dst,ptrdiff_t dst_stride,const InterpKernel * filter,int x0_q4,int32_t x_step_q4,int y0_q4,int y_step_q4,int w,int h)629 void vpx_convolve8_vert_msa(const uint8_t *src, ptrdiff_t src_stride,
630                             uint8_t *dst, ptrdiff_t dst_stride,
631                             const InterpKernel *filter, int x0_q4,
632                             int32_t x_step_q4, int y0_q4, int y_step_q4, int w,
633                             int h) {
634   const int16_t *const filter_y = filter[y0_q4];
635   int8_t cnt, filt_ver[8];
636 
637   assert(y_step_q4 == 16);
638   assert(((const int32_t *)filter_y)[1] != 0x800000);
639 
640   for (cnt = 8; cnt--;) {
641     filt_ver[cnt] = filter_y[cnt];
642   }
643 
644   if (vpx_get_filter_taps(filter_y) == 2) {
645     switch (w) {
646       case 4:
647         common_vt_2t_4w_msa(src, (int32_t)src_stride, dst, (int32_t)dst_stride,
648                             &filt_ver[3], h);
649         break;
650       case 8:
651         common_vt_2t_8w_msa(src, (int32_t)src_stride, dst, (int32_t)dst_stride,
652                             &filt_ver[3], h);
653         break;
654       case 16:
655         common_vt_2t_16w_msa(src, (int32_t)src_stride, dst, (int32_t)dst_stride,
656                              &filt_ver[3], h);
657         break;
658       case 32:
659         common_vt_2t_32w_msa(src, (int32_t)src_stride, dst, (int32_t)dst_stride,
660                              &filt_ver[3], h);
661         break;
662       case 64:
663         common_vt_2t_64w_msa(src, (int32_t)src_stride, dst, (int32_t)dst_stride,
664                              &filt_ver[3], h);
665         break;
666       default:
667         vpx_convolve8_vert_c(src, src_stride, dst, dst_stride, filter, x0_q4,
668                              x_step_q4, y0_q4, y_step_q4, w, h);
669         break;
670     }
671   } else {
672     switch (w) {
673       case 4:
674         common_vt_8t_4w_msa(src, (int32_t)src_stride, dst, (int32_t)dst_stride,
675                             filt_ver, h);
676         break;
677       case 8:
678         common_vt_8t_8w_msa(src, (int32_t)src_stride, dst, (int32_t)dst_stride,
679                             filt_ver, h);
680         break;
681       case 16:
682         common_vt_8t_16w_msa(src, (int32_t)src_stride, dst, (int32_t)dst_stride,
683                              filt_ver, h);
684         break;
685       case 32:
686         common_vt_8t_32w_msa(src, (int32_t)src_stride, dst, (int32_t)dst_stride,
687                              filt_ver, h);
688         break;
689       case 64:
690         common_vt_8t_64w_msa(src, (int32_t)src_stride, dst, (int32_t)dst_stride,
691                              filt_ver, h);
692         break;
693       default:
694         vpx_convolve8_vert_c(src, src_stride, dst, dst_stride, filter, x0_q4,
695                              x_step_q4, y0_q4, y_step_q4, w, h);
696         break;
697     }
698   }
699 }
700