• Home
  • Line#
  • Scopes#
  • Navigate#
  • Raw
  • Download
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 "vpx_dsp/mips/inv_txfm_msa.h"
12 
vpx_idct16_1d_rows_msa(const int16_t * input,int16_t * output)13 void vpx_idct16_1d_rows_msa(const int16_t *input, int16_t *output) {
14   v8i16 loc0, loc1, loc2, loc3;
15   v8i16 reg0, reg2, reg4, reg6, reg8, reg10, reg12, reg14;
16   v8i16 reg3, reg13, reg11, reg5, reg7, reg9, reg1, reg15;
17   v8i16 tmp5, tmp6, tmp7;
18 
19   LD_SH8(input, 16, reg0, reg1, reg2, reg3, reg4, reg5, reg6, reg7);
20   input += 8;
21   LD_SH8(input, 16, reg8, reg9, reg10, reg11, reg12, reg13, reg14, reg15);
22 
23   TRANSPOSE8x8_SH_SH(reg0, reg1, reg2, reg3, reg4, reg5, reg6, reg7, reg0, reg1,
24                      reg2, reg3, reg4, reg5, reg6, reg7);
25   TRANSPOSE8x8_SH_SH(reg8, reg9, reg10, reg11, reg12, reg13, reg14, reg15, reg8,
26                      reg9, reg10, reg11, reg12, reg13, reg14, reg15);
27   DOTP_CONST_PAIR(reg2, reg14, cospi_28_64, cospi_4_64, reg2, reg14);
28   DOTP_CONST_PAIR(reg10, reg6, cospi_12_64, cospi_20_64, reg10, reg6);
29   BUTTERFLY_4(reg2, reg14, reg6, reg10, loc0, loc1, reg14, reg2);
30   DOTP_CONST_PAIR(reg14, reg2, cospi_16_64, cospi_16_64, loc2, loc3);
31   DOTP_CONST_PAIR(reg0, reg8, cospi_16_64, cospi_16_64, reg0, reg8);
32   DOTP_CONST_PAIR(reg4, reg12, cospi_24_64, cospi_8_64, reg4, reg12);
33   BUTTERFLY_4(reg8, reg0, reg4, reg12, reg2, reg6, reg10, reg14);
34   SUB4(reg2, loc1, reg14, loc0, reg6, loc3, reg10, loc2, reg0, reg12, reg4,
35        reg8);
36   ADD4(reg2, loc1, reg14, loc0, reg6, loc3, reg10, loc2, reg2, reg14, reg6,
37        reg10);
38 
39   /* stage 2 */
40   DOTP_CONST_PAIR(reg1, reg15, cospi_30_64, cospi_2_64, reg1, reg15);
41   DOTP_CONST_PAIR(reg9, reg7, cospi_14_64, cospi_18_64, loc2, loc3);
42 
43   reg9 = reg1 - loc2;
44   reg1 = reg1 + loc2;
45   reg7 = reg15 - loc3;
46   reg15 = reg15 + loc3;
47 
48   DOTP_CONST_PAIR(reg5, reg11, cospi_22_64, cospi_10_64, reg5, reg11);
49   DOTP_CONST_PAIR(reg13, reg3, cospi_6_64, cospi_26_64, loc0, loc1);
50   BUTTERFLY_4(loc0, loc1, reg11, reg5, reg13, reg3, reg11, reg5);
51 
52   loc1 = reg15 + reg3;
53   reg3 = reg15 - reg3;
54   loc2 = reg2 + loc1;
55   reg15 = reg2 - loc1;
56 
57   loc1 = reg1 + reg13;
58   reg13 = reg1 - reg13;
59   loc0 = reg0 + loc1;
60   loc1 = reg0 - loc1;
61   tmp6 = loc0;
62   tmp7 = loc1;
63   reg0 = loc2;
64 
65   DOTP_CONST_PAIR(reg7, reg9, cospi_24_64, cospi_8_64, reg7, reg9);
66   DOTP_CONST_PAIR((-reg5), (-reg11), cospi_8_64, cospi_24_64, reg5, reg11);
67 
68   loc0 = reg9 + reg5;
69   reg5 = reg9 - reg5;
70   reg2 = reg6 + loc0;
71   reg1 = reg6 - loc0;
72 
73   loc0 = reg7 + reg11;
74   reg11 = reg7 - reg11;
75   loc1 = reg4 + loc0;
76   loc2 = reg4 - loc0;
77   tmp5 = loc1;
78 
79   DOTP_CONST_PAIR(reg5, reg11, cospi_16_64, cospi_16_64, reg5, reg11);
80   BUTTERFLY_4(reg8, reg10, reg11, reg5, loc0, reg4, reg9, loc1);
81 
82   reg10 = loc0;
83   reg11 = loc1;
84 
85   DOTP_CONST_PAIR(reg3, reg13, cospi_16_64, cospi_16_64, reg3, reg13);
86   BUTTERFLY_4(reg12, reg14, reg13, reg3, reg8, reg6, reg7, reg5);
87 
88   reg13 = loc2;
89 
90   /* Transpose and store the output */
91   reg12 = tmp5;
92   reg14 = tmp6;
93   reg3 = tmp7;
94 
95   /* transpose block */
96   TRANSPOSE8x8_SH_SH(reg0, reg2, reg4, reg6, reg8, reg10, reg12, reg14, reg0,
97                      reg2, reg4, reg6, reg8, reg10, reg12, reg14);
98   ST_SH8(reg0, reg2, reg4, reg6, reg8, reg10, reg12, reg14, output, 16);
99 
100   /* transpose block */
101   TRANSPOSE8x8_SH_SH(reg3, reg13, reg11, reg5, reg7, reg9, reg1, reg15, reg3,
102                      reg13, reg11, reg5, reg7, reg9, reg1, reg15);
103   ST_SH8(reg3, reg13, reg11, reg5, reg7, reg9, reg1, reg15, (output + 8), 16);
104 }
105 
vpx_idct16_1d_columns_addblk_msa(int16_t * input,uint8_t * dst,int32_t dst_stride)106 void vpx_idct16_1d_columns_addblk_msa(int16_t *input, uint8_t *dst,
107                                       int32_t dst_stride) {
108   v8i16 loc0, loc1, loc2, loc3;
109   v8i16 reg0, reg2, reg4, reg6, reg8, reg10, reg12, reg14;
110   v8i16 reg3, reg13, reg11, reg5, reg7, reg9, reg1, reg15;
111   v8i16 tmp5, tmp6, tmp7;
112 
113   /* load up 8x8 */
114   LD_SH8(input, 16, reg0, reg1, reg2, reg3, reg4, reg5, reg6, reg7);
115   input += 8 * 16;
116   /* load bottom 8x8 */
117   LD_SH8(input, 16, reg8, reg9, reg10, reg11, reg12, reg13, reg14, reg15);
118 
119   DOTP_CONST_PAIR(reg2, reg14, cospi_28_64, cospi_4_64, reg2, reg14);
120   DOTP_CONST_PAIR(reg10, reg6, cospi_12_64, cospi_20_64, reg10, reg6);
121   BUTTERFLY_4(reg2, reg14, reg6, reg10, loc0, loc1, reg14, reg2);
122   DOTP_CONST_PAIR(reg14, reg2, cospi_16_64, cospi_16_64, loc2, loc3);
123   DOTP_CONST_PAIR(reg0, reg8, cospi_16_64, cospi_16_64, reg0, reg8);
124   DOTP_CONST_PAIR(reg4, reg12, cospi_24_64, cospi_8_64, reg4, reg12);
125   BUTTERFLY_4(reg8, reg0, reg4, reg12, reg2, reg6, reg10, reg14);
126 
127   reg0 = reg2 - loc1;
128   reg2 = reg2 + loc1;
129   reg12 = reg14 - loc0;
130   reg14 = reg14 + loc0;
131   reg4 = reg6 - loc3;
132   reg6 = reg6 + loc3;
133   reg8 = reg10 - loc2;
134   reg10 = reg10 + loc2;
135 
136   /* stage 2 */
137   DOTP_CONST_PAIR(reg1, reg15, cospi_30_64, cospi_2_64, reg1, reg15);
138   DOTP_CONST_PAIR(reg9, reg7, cospi_14_64, cospi_18_64, loc2, loc3);
139 
140   reg9 = reg1 - loc2;
141   reg1 = reg1 + loc2;
142   reg7 = reg15 - loc3;
143   reg15 = reg15 + loc3;
144 
145   DOTP_CONST_PAIR(reg5, reg11, cospi_22_64, cospi_10_64, reg5, reg11);
146   DOTP_CONST_PAIR(reg13, reg3, cospi_6_64, cospi_26_64, loc0, loc1);
147   BUTTERFLY_4(loc0, loc1, reg11, reg5, reg13, reg3, reg11, reg5);
148 
149   loc1 = reg15 + reg3;
150   reg3 = reg15 - reg3;
151   loc2 = reg2 + loc1;
152   reg15 = reg2 - loc1;
153 
154   loc1 = reg1 + reg13;
155   reg13 = reg1 - reg13;
156   loc0 = reg0 + loc1;
157   loc1 = reg0 - loc1;
158   tmp6 = loc0;
159   tmp7 = loc1;
160   reg0 = loc2;
161 
162   DOTP_CONST_PAIR(reg7, reg9, cospi_24_64, cospi_8_64, reg7, reg9);
163   DOTP_CONST_PAIR((-reg5), (-reg11), cospi_8_64, cospi_24_64, reg5, reg11);
164 
165   loc0 = reg9 + reg5;
166   reg5 = reg9 - reg5;
167   reg2 = reg6 + loc0;
168   reg1 = reg6 - loc0;
169 
170   loc0 = reg7 + reg11;
171   reg11 = reg7 - reg11;
172   loc1 = reg4 + loc0;
173   loc2 = reg4 - loc0;
174   tmp5 = loc1;
175 
176   DOTP_CONST_PAIR(reg5, reg11, cospi_16_64, cospi_16_64, reg5, reg11);
177   BUTTERFLY_4(reg8, reg10, reg11, reg5, loc0, reg4, reg9, loc1);
178 
179   reg10 = loc0;
180   reg11 = loc1;
181 
182   DOTP_CONST_PAIR(reg3, reg13, cospi_16_64, cospi_16_64, reg3, reg13);
183   BUTTERFLY_4(reg12, reg14, reg13, reg3, reg8, reg6, reg7, reg5);
184   reg13 = loc2;
185 
186   /* Transpose and store the output */
187   reg12 = tmp5;
188   reg14 = tmp6;
189   reg3 = tmp7;
190 
191   SRARI_H4_SH(reg0, reg2, reg4, reg6, 6);
192   VP9_ADDBLK_ST8x4_UB(dst, dst_stride, reg0, reg2, reg4, reg6);
193   dst += (4 * dst_stride);
194   SRARI_H4_SH(reg8, reg10, reg12, reg14, 6);
195   VP9_ADDBLK_ST8x4_UB(dst, dst_stride, reg8, reg10, reg12, reg14);
196   dst += (4 * dst_stride);
197   SRARI_H4_SH(reg3, reg13, reg11, reg5, 6);
198   VP9_ADDBLK_ST8x4_UB(dst, dst_stride, reg3, reg13, reg11, reg5);
199   dst += (4 * dst_stride);
200   SRARI_H4_SH(reg7, reg9, reg1, reg15, 6);
201   VP9_ADDBLK_ST8x4_UB(dst, dst_stride, reg7, reg9, reg1, reg15);
202 }
203 
vpx_idct16x16_256_add_msa(const int16_t * input,uint8_t * dst,int32_t dst_stride)204 void vpx_idct16x16_256_add_msa(const int16_t *input, uint8_t *dst,
205                                int32_t dst_stride) {
206   int32_t i;
207   DECLARE_ALIGNED(32, int16_t, out_arr[16 * 16]);
208   int16_t *out = out_arr;
209 
210   /* transform rows */
211   for (i = 0; i < 2; ++i) {
212     /* process 16 * 8 block */
213     vpx_idct16_1d_rows_msa((input + (i << 7)), (out + (i << 7)));
214   }
215 
216   /* transform columns */
217   for (i = 0; i < 2; ++i) {
218     /* process 8 * 16 block */
219     vpx_idct16_1d_columns_addblk_msa((out + (i << 3)), (dst + (i << 3)),
220                                      dst_stride);
221   }
222 }
223 
vpx_idct16x16_10_add_msa(const int16_t * input,uint8_t * dst,int32_t dst_stride)224 void vpx_idct16x16_10_add_msa(const int16_t *input, uint8_t *dst,
225                               int32_t dst_stride) {
226   uint8_t i;
227   DECLARE_ALIGNED(32, int16_t, out_arr[16 * 16]);
228   int16_t *out = out_arr;
229 
230   /* process 16 * 8 block */
231   vpx_idct16_1d_rows_msa(input, out);
232 
233   /* short case just considers top 4 rows as valid output */
234   out += 4 * 16;
235   for (i = 12; i--;) {
236     __asm__ __volatile__(
237         "sw     $zero,   0(%[out])     \n\t"
238         "sw     $zero,   4(%[out])     \n\t"
239         "sw     $zero,   8(%[out])     \n\t"
240         "sw     $zero,  12(%[out])     \n\t"
241         "sw     $zero,  16(%[out])     \n\t"
242         "sw     $zero,  20(%[out])     \n\t"
243         "sw     $zero,  24(%[out])     \n\t"
244         "sw     $zero,  28(%[out])     \n\t"
245 
246         :
247         : [out] "r"(out));
248 
249     out += 16;
250   }
251 
252   out = out_arr;
253 
254   /* transform columns */
255   for (i = 0; i < 2; ++i) {
256     /* process 8 * 16 block */
257     vpx_idct16_1d_columns_addblk_msa((out + (i << 3)), (dst + (i << 3)),
258                                      dst_stride);
259   }
260 }
261 
vpx_idct16x16_1_add_msa(const int16_t * input,uint8_t * dst,int32_t dst_stride)262 void vpx_idct16x16_1_add_msa(const int16_t *input, uint8_t *dst,
263                              int32_t dst_stride) {
264   uint8_t i;
265   int16_t out;
266   v8i16 vec, res0, res1, res2, res3, res4, res5, res6, res7;
267   v16u8 dst0, dst1, dst2, dst3, tmp0, tmp1, tmp2, tmp3;
268 
269   out = ROUND_POWER_OF_TWO((input[0] * cospi_16_64), DCT_CONST_BITS);
270   out = ROUND_POWER_OF_TWO((out * cospi_16_64), DCT_CONST_BITS);
271   out = ROUND_POWER_OF_TWO(out, 6);
272 
273   vec = __msa_fill_h(out);
274 
275   for (i = 4; i--;) {
276     LD_UB4(dst, dst_stride, dst0, dst1, dst2, dst3);
277     UNPCK_UB_SH(dst0, res0, res4);
278     UNPCK_UB_SH(dst1, res1, res5);
279     UNPCK_UB_SH(dst2, res2, res6);
280     UNPCK_UB_SH(dst3, res3, res7);
281     ADD4(res0, vec, res1, vec, res2, vec, res3, vec, res0, res1, res2, res3);
282     ADD4(res4, vec, res5, vec, res6, vec, res7, vec, res4, res5, res6, res7);
283     CLIP_SH4_0_255(res0, res1, res2, res3);
284     CLIP_SH4_0_255(res4, res5, res6, res7);
285     PCKEV_B4_UB(res4, res0, res5, res1, res6, res2, res7, res3, tmp0, tmp1,
286                 tmp2, tmp3);
287     ST_UB4(tmp0, tmp1, tmp2, tmp3, dst, dst_stride);
288     dst += (4 * dst_stride);
289   }
290 }
291 
vpx_iadst16_1d_rows_msa(const int16_t * input,int16_t * output)292 void vpx_iadst16_1d_rows_msa(const int16_t *input, int16_t *output) {
293   v8i16 r0, r1, r2, r3, r4, r5, r6, r7, r8, r9, r10, r11, r12, r13, r14, r15;
294   v8i16 l0, l1, l2, l3, l4, l5, l6, l7, l8, l9, l10, l11, l12, l13, l14, l15;
295 
296   /* load input data */
297   LD_SH16(input, 8, l0, l8, l1, l9, l2, l10, l3, l11, l4, l12, l5, l13, l6, l14,
298           l7, l15);
299   TRANSPOSE8x8_SH_SH(l0, l1, l2, l3, l4, l5, l6, l7, l0, l1, l2, l3, l4, l5, l6,
300                      l7);
301   TRANSPOSE8x8_SH_SH(l8, l9, l10, l11, l12, l13, l14, l15, l8, l9, l10, l11,
302                      l12, l13, l14, l15);
303 
304   /* ADST in horizontal */
305   VP9_IADST8x16_1D(l0, l1, l2, l3, l4, l5, l6, l7, l8, l9, l10, l11, l12, l13,
306                    l14, l15, r0, r1, r2, r3, r4, r5, r6, r7, r8, r9, r10, r11,
307                    r12, r13, r14, r15);
308 
309   l1 = -r8;
310   l3 = -r4;
311   l13 = -r13;
312   l15 = -r1;
313 
314   TRANSPOSE8x8_SH_SH(r0, l1, r12, l3, r6, r14, r10, r2, l0, l1, l2, l3, l4, l5,
315                      l6, l7);
316   ST_SH8(l0, l1, l2, l3, l4, l5, l6, l7, output, 16);
317   TRANSPOSE8x8_SH_SH(r3, r11, r15, r7, r5, l13, r9, l15, l8, l9, l10, l11, l12,
318                      l13, l14, l15);
319   ST_SH8(l8, l9, l10, l11, l12, l13, l14, l15, (output + 8), 16);
320 }
321 
vpx_iadst16_1d_columns_addblk_msa(int16_t * input,uint8_t * dst,int32_t dst_stride)322 void vpx_iadst16_1d_columns_addblk_msa(int16_t *input, uint8_t *dst,
323                                        int32_t dst_stride) {
324   v8i16 v0, v2, v4, v6, k0, k1, k2, k3;
325   v8i16 r0, r1, r2, r3, r4, r5, r6, r7, r8, r9, r10, r11, r12, r13, r14, r15;
326   v8i16 out0, out1, out2, out3, out4, out5, out6, out7;
327   v8i16 out8, out9, out10, out11, out12, out13, out14, out15;
328   v8i16 g0, g1, g2, g3, g4, g5, g6, g7, g8, g9, g10, g11, g12, g13, g14, g15;
329   v8i16 h0, h1, h2, h3, h4, h5, h6, h7, h8, h9, h10, h11;
330   v8i16 res0, res1, res2, res3, res4, res5, res6, res7;
331   v8i16 res8, res9, res10, res11, res12, res13, res14, res15;
332   v16u8 dst0, dst1, dst2, dst3, dst4, dst5, dst6, dst7;
333   v16u8 dst8, dst9, dst10, dst11, dst12, dst13, dst14, dst15;
334   v16i8 zero = { 0 };
335 
336   r0 = LD_SH(input + 0 * 16);
337   r3 = LD_SH(input + 3 * 16);
338   r4 = LD_SH(input + 4 * 16);
339   r7 = LD_SH(input + 7 * 16);
340   r8 = LD_SH(input + 8 * 16);
341   r11 = LD_SH(input + 11 * 16);
342   r12 = LD_SH(input + 12 * 16);
343   r15 = LD_SH(input + 15 * 16);
344 
345   /* stage 1 */
346   k0 = VP9_SET_COSPI_PAIR(cospi_1_64, cospi_31_64);
347   k1 = VP9_SET_COSPI_PAIR(cospi_31_64, -cospi_1_64);
348   k2 = VP9_SET_COSPI_PAIR(cospi_17_64, cospi_15_64);
349   k3 = VP9_SET_COSPI_PAIR(cospi_15_64, -cospi_17_64);
350   MADD_BF(r15, r0, r7, r8, k0, k1, k2, k3, g0, g1, g2, g3);
351   k0 = VP9_SET_COSPI_PAIR(cospi_9_64, cospi_23_64);
352   k1 = VP9_SET_COSPI_PAIR(cospi_23_64, -cospi_9_64);
353   k2 = VP9_SET_COSPI_PAIR(cospi_25_64, cospi_7_64);
354   k3 = VP9_SET_COSPI_PAIR(cospi_7_64, -cospi_25_64);
355   MADD_BF(r11, r4, r3, r12, k0, k1, k2, k3, g8, g9, g10, g11);
356   BUTTERFLY_4(g0, g2, g10, g8, h8, h9, v2, v0);
357   k0 = VP9_SET_COSPI_PAIR(cospi_4_64, cospi_28_64);
358   k1 = VP9_SET_COSPI_PAIR(cospi_28_64, -cospi_4_64);
359   k2 = VP9_SET_COSPI_PAIR(-cospi_28_64, cospi_4_64);
360   MADD_BF(g1, g3, g9, g11, k0, k1, k2, k0, h0, h1, h2, h3);
361 
362   r1 = LD_SH(input + 1 * 16);
363   r2 = LD_SH(input + 2 * 16);
364   r5 = LD_SH(input + 5 * 16);
365   r6 = LD_SH(input + 6 * 16);
366   r9 = LD_SH(input + 9 * 16);
367   r10 = LD_SH(input + 10 * 16);
368   r13 = LD_SH(input + 13 * 16);
369   r14 = LD_SH(input + 14 * 16);
370 
371   k0 = VP9_SET_COSPI_PAIR(cospi_5_64, cospi_27_64);
372   k1 = VP9_SET_COSPI_PAIR(cospi_27_64, -cospi_5_64);
373   k2 = VP9_SET_COSPI_PAIR(cospi_21_64, cospi_11_64);
374   k3 = VP9_SET_COSPI_PAIR(cospi_11_64, -cospi_21_64);
375   MADD_BF(r13, r2, r5, r10, k0, k1, k2, k3, g4, g5, g6, g7);
376   k0 = VP9_SET_COSPI_PAIR(cospi_13_64, cospi_19_64);
377   k1 = VP9_SET_COSPI_PAIR(cospi_19_64, -cospi_13_64);
378   k2 = VP9_SET_COSPI_PAIR(cospi_29_64, cospi_3_64);
379   k3 = VP9_SET_COSPI_PAIR(cospi_3_64, -cospi_29_64);
380   MADD_BF(r9, r6, r1, r14, k0, k1, k2, k3, g12, g13, g14, g15);
381   BUTTERFLY_4(g4, g6, g14, g12, h10, h11, v6, v4);
382   BUTTERFLY_4(h8, h9, h11, h10, out0, out1, h11, h10);
383   out1 = -out1;
384   SRARI_H2_SH(out0, out1, 6);
385   dst0 = LD_UB(dst + 0 * dst_stride);
386   dst1 = LD_UB(dst + 15 * dst_stride);
387   ILVR_B2_SH(zero, dst0, zero, dst1, res0, res1);
388   ADD2(res0, out0, res1, out1, res0, res1);
389   CLIP_SH2_0_255(res0, res1);
390   PCKEV_B2_SH(res0, res0, res1, res1, res0, res1);
391   ST8x1_UB(res0, dst);
392   ST8x1_UB(res1, dst + 15 * dst_stride);
393 
394   k0 = VP9_SET_COSPI_PAIR(cospi_12_64, cospi_20_64);
395   k1 = VP9_SET_COSPI_PAIR(-cospi_20_64, cospi_12_64);
396   k2 = VP9_SET_COSPI_PAIR(cospi_20_64, -cospi_12_64);
397   MADD_BF(g7, g5, g15, g13, k0, k1, k2, k0, h4, h5, h6, h7);
398   BUTTERFLY_4(h0, h2, h6, h4, out8, out9, out11, out10);
399   out8 = -out8;
400 
401   SRARI_H2_SH(out8, out9, 6);
402   dst8 = LD_UB(dst + 1 * dst_stride);
403   dst9 = LD_UB(dst + 14 * dst_stride);
404   ILVR_B2_SH(zero, dst8, zero, dst9, res8, res9);
405   ADD2(res8, out8, res9, out9, res8, res9);
406   CLIP_SH2_0_255(res8, res9);
407   PCKEV_B2_SH(res8, res8, res9, res9, res8, res9);
408   ST8x1_UB(res8, dst + dst_stride);
409   ST8x1_UB(res9, dst + 14 * dst_stride);
410 
411   k0 = VP9_SET_COSPI_PAIR(cospi_8_64, cospi_24_64);
412   k1 = VP9_SET_COSPI_PAIR(cospi_24_64, -cospi_8_64);
413   k2 = VP9_SET_COSPI_PAIR(-cospi_24_64, cospi_8_64);
414   MADD_BF(v0, v2, v4, v6, k0, k1, k2, k0, out4, out6, out5, out7);
415   out4 = -out4;
416   SRARI_H2_SH(out4, out5, 6);
417   dst4 = LD_UB(dst + 3 * dst_stride);
418   dst5 = LD_UB(dst + 12 * dst_stride);
419   ILVR_B2_SH(zero, dst4, zero, dst5, res4, res5);
420   ADD2(res4, out4, res5, out5, res4, res5);
421   CLIP_SH2_0_255(res4, res5);
422   PCKEV_B2_SH(res4, res4, res5, res5, res4, res5);
423   ST8x1_UB(res4, dst + 3 * dst_stride);
424   ST8x1_UB(res5, dst + 12 * dst_stride);
425 
426   MADD_BF(h1, h3, h5, h7, k0, k1, k2, k0, out12, out14, out13, out15);
427   out13 = -out13;
428   SRARI_H2_SH(out12, out13, 6);
429   dst12 = LD_UB(dst + 2 * dst_stride);
430   dst13 = LD_UB(dst + 13 * dst_stride);
431   ILVR_B2_SH(zero, dst12, zero, dst13, res12, res13);
432   ADD2(res12, out12, res13, out13, res12, res13);
433   CLIP_SH2_0_255(res12, res13);
434   PCKEV_B2_SH(res12, res12, res13, res13, res12, res13);
435   ST8x1_UB(res12, dst + 2 * dst_stride);
436   ST8x1_UB(res13, dst + 13 * dst_stride);
437 
438   k0 = VP9_SET_COSPI_PAIR(cospi_16_64, cospi_16_64);
439   k3 = VP9_SET_COSPI_PAIR(-cospi_16_64, cospi_16_64);
440   MADD_SHORT(out6, out7, k0, k3, out6, out7);
441   SRARI_H2_SH(out6, out7, 6);
442   dst6 = LD_UB(dst + 4 * dst_stride);
443   dst7 = LD_UB(dst + 11 * dst_stride);
444   ILVR_B2_SH(zero, dst6, zero, dst7, res6, res7);
445   ADD2(res6, out6, res7, out7, res6, res7);
446   CLIP_SH2_0_255(res6, res7);
447   PCKEV_B2_SH(res6, res6, res7, res7, res6, res7);
448   ST8x1_UB(res6, dst + 4 * dst_stride);
449   ST8x1_UB(res7, dst + 11 * dst_stride);
450 
451   MADD_SHORT(out10, out11, k0, k3, out10, out11);
452   SRARI_H2_SH(out10, out11, 6);
453   dst10 = LD_UB(dst + 6 * dst_stride);
454   dst11 = LD_UB(dst + 9 * dst_stride);
455   ILVR_B2_SH(zero, dst10, zero, dst11, res10, res11);
456   ADD2(res10, out10, res11, out11, res10, res11);
457   CLIP_SH2_0_255(res10, res11);
458   PCKEV_B2_SH(res10, res10, res11, res11, res10, res11);
459   ST8x1_UB(res10, dst + 6 * dst_stride);
460   ST8x1_UB(res11, dst + 9 * dst_stride);
461 
462   k1 = VP9_SET_COSPI_PAIR(-cospi_16_64, -cospi_16_64);
463   k2 = VP9_SET_COSPI_PAIR(cospi_16_64, -cospi_16_64);
464   MADD_SHORT(h10, h11, k1, k2, out2, out3);
465   SRARI_H2_SH(out2, out3, 6);
466   dst2 = LD_UB(dst + 7 * dst_stride);
467   dst3 = LD_UB(dst + 8 * dst_stride);
468   ILVR_B2_SH(zero, dst2, zero, dst3, res2, res3);
469   ADD2(res2, out2, res3, out3, res2, res3);
470   CLIP_SH2_0_255(res2, res3);
471   PCKEV_B2_SH(res2, res2, res3, res3, res2, res3);
472   ST8x1_UB(res2, dst + 7 * dst_stride);
473   ST8x1_UB(res3, dst + 8 * dst_stride);
474 
475   MADD_SHORT(out14, out15, k1, k2, out14, out15);
476   SRARI_H2_SH(out14, out15, 6);
477   dst14 = LD_UB(dst + 5 * dst_stride);
478   dst15 = LD_UB(dst + 10 * dst_stride);
479   ILVR_B2_SH(zero, dst14, zero, dst15, res14, res15);
480   ADD2(res14, out14, res15, out15, res14, res15);
481   CLIP_SH2_0_255(res14, res15);
482   PCKEV_B2_SH(res14, res14, res15, res15, res14, res15);
483   ST8x1_UB(res14, dst + 5 * dst_stride);
484   ST8x1_UB(res15, dst + 10 * dst_stride);
485 }
486