avcodec/vp9dsp: fix overwrite by 1 in vert_left pred.
[ffmpeg.git] / libavcodec / vp9dsp.c
1 /*
2  * VP9 compatible video decoder
3  *
4  * Copyright (C) 2013 Ronald S. Bultje <rsbultje gmail com>
5  * Copyright (C) 2013 Clément Bœsch <u pkh me>
6  *
7  * This file is part of FFmpeg.
8  *
9  * FFmpeg is free software; you can redistribute it and/or
10  * modify it under the terms of the GNU Lesser General Public
11  * License as published by the Free Software Foundation; either
12  * version 2.1 of the License, or (at your option) any later version.
13  *
14  * FFmpeg is distributed in the hope that it will be useful,
15  * but WITHOUT ANY WARRANTY; without even the implied warranty of
16  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
17  * Lesser General Public License for more details.
18  *
19  * You should have received a copy of the GNU Lesser General Public
20  * License along with FFmpeg; if not, write to the Free Software
21  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
22  */
23
24 #include "libavutil/common.h"
25 #include "libavutil/intreadwrite.h"
26 #include "vp9dsp.h"
27 #include "rnd_avg.h"
28
29 // FIXME see whether we can merge parts of this (perhaps at least 4x4 and 8x8)
30 // back with h264pred.[ch]
31
32 static void vert_4x4_c(uint8_t *dst, ptrdiff_t stride,
33                        const uint8_t *left, const uint8_t *top)
34 {
35     unsigned p4 = AV_RN32A(top);
36
37     AV_WN32A(dst + stride * 0, p4);
38     AV_WN32A(dst + stride * 1, p4);
39     AV_WN32A(dst + stride * 2, p4);
40     AV_WN32A(dst + stride * 3, p4);
41 }
42
43 static void vert_8x8_c(uint8_t *dst, ptrdiff_t stride,
44                        const uint8_t *left, const uint8_t *top)
45 {
46     uint64_t p8 = AV_RN64A(top);
47     int y;
48
49     for (y = 0; y < 8; y++) {
50         AV_WN64A(dst, p8);
51         dst += stride;
52     }
53 }
54
55 static void vert_16x16_c(uint8_t *dst, ptrdiff_t stride,
56                          const uint8_t *left, const uint8_t *top)
57 {
58     uint64_t p8a = AV_RN64A(top + 0), p8b = AV_RN64A(top + 8);
59     int y;
60
61     for (y = 0; y < 16; y++) {
62         AV_WN64A(dst + 0, p8a);
63         AV_WN64A(dst + 8, p8b);
64         dst += stride;
65     }
66 }
67
68 static void vert_32x32_c(uint8_t *dst, ptrdiff_t stride,
69                          const uint8_t *left, const uint8_t *top)
70 {
71     uint64_t p8a = AV_RN64A(top + 0),  p8b = AV_RN64A(top + 8),
72              p8c = AV_RN64A(top + 16), p8d = AV_RN64A(top + 24);
73     int y;
74
75     for (y = 0; y < 32; y++) {
76         AV_WN64A(dst +  0, p8a);
77         AV_WN64A(dst +  8, p8b);
78         AV_WN64A(dst + 16, p8c);
79         AV_WN64A(dst + 24, p8d);
80         dst += stride;
81     }
82 }
83
84 static void hor_4x4_c(uint8_t *dst, ptrdiff_t stride,
85                       const uint8_t *left, const uint8_t *top)
86 {
87     AV_WN32A(dst + stride * 0, left[0] * 0x01010101U);
88     AV_WN32A(dst + stride * 1, left[1] * 0x01010101U);
89     AV_WN32A(dst + stride * 2, left[2] * 0x01010101U);
90     AV_WN32A(dst + stride * 3, left[3] * 0x01010101U);
91 }
92
93 static void hor_8x8_c(uint8_t *dst, ptrdiff_t stride,
94                       const uint8_t *left, const uint8_t *top)
95 {
96     int y;
97
98     for (y = 0; y < 8; y++) {
99         AV_WN64A(dst, left[y] * 0x0101010101010101ULL);
100         dst += stride;
101     }
102 }
103
104 static void hor_16x16_c(uint8_t *dst, ptrdiff_t stride,
105                         const uint8_t *left, const uint8_t *top)
106 {
107     int y;
108
109     for (y = 0; y < 16; y++) {
110         uint64_t p8 = left[y] * 0x0101010101010101ULL;
111
112         AV_WN64A(dst + 0, p8);
113         AV_WN64A(dst + 8, p8);
114         dst += stride;
115     }
116 }
117
118 static void hor_32x32_c(uint8_t *dst, ptrdiff_t stride,
119                         const uint8_t *left, const uint8_t *top)
120 {
121     int y;
122
123     for (y = 0; y < 32; y++) {
124         uint64_t p8 = left[y] * 0x0101010101010101ULL;
125
126         AV_WN64A(dst +  0, p8);
127         AV_WN64A(dst +  8, p8);
128         AV_WN64A(dst + 16, p8);
129         AV_WN64A(dst + 24, p8);
130         dst += stride;
131     }
132 }
133
134 static void tm_4x4_c(uint8_t *dst, ptrdiff_t stride,
135                      const uint8_t *left, const uint8_t *top)
136 {
137     int y, tl = top[-1];
138
139     for (y = 0; y < 4; y++) {
140         int l_m_tl = left[y] - tl;
141
142         dst[0] = av_clip_uint8(top[0] + l_m_tl);
143         dst[1] = av_clip_uint8(top[1] + l_m_tl);
144         dst[2] = av_clip_uint8(top[2] + l_m_tl);
145         dst[3] = av_clip_uint8(top[3] + l_m_tl);
146         dst += stride;
147     }
148 }
149
150 static void tm_8x8_c(uint8_t *dst, ptrdiff_t stride,
151                      const uint8_t *left, const uint8_t *top)
152 {
153     int y, tl = top[-1];
154
155     for (y = 0; y < 8; y++) {
156         int l_m_tl = left[y] - tl;
157
158         dst[0] = av_clip_uint8(top[0] + l_m_tl);
159         dst[1] = av_clip_uint8(top[1] + l_m_tl);
160         dst[2] = av_clip_uint8(top[2] + l_m_tl);
161         dst[3] = av_clip_uint8(top[3] + l_m_tl);
162         dst[4] = av_clip_uint8(top[4] + l_m_tl);
163         dst[5] = av_clip_uint8(top[5] + l_m_tl);
164         dst[6] = av_clip_uint8(top[6] + l_m_tl);
165         dst[7] = av_clip_uint8(top[7] + l_m_tl);
166         dst += stride;
167     }
168 }
169
170 static void tm_16x16_c(uint8_t *dst, ptrdiff_t stride,
171                        const uint8_t *left, const uint8_t *top)
172 {
173     int y, tl = top[-1];
174
175     for (y = 0; y < 16; y++) {
176         int l_m_tl = left[y] - tl;
177
178         dst[ 0] = av_clip_uint8(top[ 0] + l_m_tl);
179         dst[ 1] = av_clip_uint8(top[ 1] + l_m_tl);
180         dst[ 2] = av_clip_uint8(top[ 2] + l_m_tl);
181         dst[ 3] = av_clip_uint8(top[ 3] + l_m_tl);
182         dst[ 4] = av_clip_uint8(top[ 4] + l_m_tl);
183         dst[ 5] = av_clip_uint8(top[ 5] + l_m_tl);
184         dst[ 6] = av_clip_uint8(top[ 6] + l_m_tl);
185         dst[ 7] = av_clip_uint8(top[ 7] + l_m_tl);
186         dst[ 8] = av_clip_uint8(top[ 8] + l_m_tl);
187         dst[ 9] = av_clip_uint8(top[ 9] + l_m_tl);
188         dst[10] = av_clip_uint8(top[10] + l_m_tl);
189         dst[11] = av_clip_uint8(top[11] + l_m_tl);
190         dst[12] = av_clip_uint8(top[12] + l_m_tl);
191         dst[13] = av_clip_uint8(top[13] + l_m_tl);
192         dst[14] = av_clip_uint8(top[14] + l_m_tl);
193         dst[15] = av_clip_uint8(top[15] + l_m_tl);
194         dst += stride;
195     }
196 }
197
198 static void tm_32x32_c(uint8_t *dst, ptrdiff_t stride,
199                        const uint8_t *left, const uint8_t *top)
200 {
201     int y, tl = top[-1];
202
203     for (y = 0; y < 32; y++) {
204         int l_m_tl = left[y] - tl;
205
206         dst[ 0] = av_clip_uint8(top[ 0] + l_m_tl);
207         dst[ 1] = av_clip_uint8(top[ 1] + l_m_tl);
208         dst[ 2] = av_clip_uint8(top[ 2] + l_m_tl);
209         dst[ 3] = av_clip_uint8(top[ 3] + l_m_tl);
210         dst[ 4] = av_clip_uint8(top[ 4] + l_m_tl);
211         dst[ 5] = av_clip_uint8(top[ 5] + l_m_tl);
212         dst[ 6] = av_clip_uint8(top[ 6] + l_m_tl);
213         dst[ 7] = av_clip_uint8(top[ 7] + l_m_tl);
214         dst[ 8] = av_clip_uint8(top[ 8] + l_m_tl);
215         dst[ 9] = av_clip_uint8(top[ 9] + l_m_tl);
216         dst[10] = av_clip_uint8(top[10] + l_m_tl);
217         dst[11] = av_clip_uint8(top[11] + l_m_tl);
218         dst[12] = av_clip_uint8(top[12] + l_m_tl);
219         dst[13] = av_clip_uint8(top[13] + l_m_tl);
220         dst[14] = av_clip_uint8(top[14] + l_m_tl);
221         dst[15] = av_clip_uint8(top[15] + l_m_tl);
222         dst[16] = av_clip_uint8(top[16] + l_m_tl);
223         dst[17] = av_clip_uint8(top[17] + l_m_tl);
224         dst[18] = av_clip_uint8(top[18] + l_m_tl);
225         dst[19] = av_clip_uint8(top[19] + l_m_tl);
226         dst[20] = av_clip_uint8(top[20] + l_m_tl);
227         dst[21] = av_clip_uint8(top[21] + l_m_tl);
228         dst[22] = av_clip_uint8(top[22] + l_m_tl);
229         dst[23] = av_clip_uint8(top[23] + l_m_tl);
230         dst[24] = av_clip_uint8(top[24] + l_m_tl);
231         dst[25] = av_clip_uint8(top[25] + l_m_tl);
232         dst[26] = av_clip_uint8(top[26] + l_m_tl);
233         dst[27] = av_clip_uint8(top[27] + l_m_tl);
234         dst[28] = av_clip_uint8(top[28] + l_m_tl);
235         dst[29] = av_clip_uint8(top[29] + l_m_tl);
236         dst[30] = av_clip_uint8(top[30] + l_m_tl);
237         dst[31] = av_clip_uint8(top[31] + l_m_tl);
238         dst += stride;
239     }
240 }
241
242 static void dc_4x4_c(uint8_t *dst, ptrdiff_t stride,
243                      const uint8_t *left, const uint8_t *top)
244 {
245     unsigned dc = 0x01010101U * ((left[0] + left[1] + left[2] + left[3] +
246                                   top[0] + top[1] + top[2] + top[3] + 4) >> 3);
247
248     AV_WN32A(dst + stride * 0, dc);
249     AV_WN32A(dst + stride * 1, dc);
250     AV_WN32A(dst + stride * 2, dc);
251     AV_WN32A(dst + stride * 3, dc);
252 }
253
254 static void dc_8x8_c(uint8_t *dst, ptrdiff_t stride,
255                      const uint8_t *left, const uint8_t *top)
256 {
257     uint64_t dc = 0x0101010101010101ULL *
258         ((left[0] + left[1] + left[2] + left[3] + left[4] + left[5] +
259           left[6] + left[7] + top[0] + top[1] + top[2] + top[3] +
260           top[4] + top[5] + top[6] + top[7] + 8) >> 4);
261     int y;
262
263     for (y = 0; y < 8; y++) {
264         AV_WN64A(dst, dc);
265         dst += stride;
266     }
267 }
268
269 static void dc_16x16_c(uint8_t *dst, ptrdiff_t stride,
270                        const uint8_t *left, const uint8_t *top)
271 {
272     uint64_t dc = 0x0101010101010101ULL *
273         ((left[0] + left[1] + left[2] + left[3] + left[4] + left[5] + left[6] +
274           left[7] + left[8] + left[9] + left[10] + left[11] + left[12] +
275           left[13] + left[14] + left[15] + top[0] + top[1] + top[2] + top[3] +
276           top[4] + top[5] + top[6] + top[7] + top[8] + top[9] + top[10] +
277           top[11] + top[12] + top[13] + top[14] + top[15] + 16) >> 5);
278     int y;
279
280     for (y = 0; y < 16; y++) {
281         AV_WN64A(dst + 0, dc);
282         AV_WN64A(dst + 8, dc);
283         dst += stride;
284     }
285 }
286
287 static void dc_32x32_c(uint8_t *dst, ptrdiff_t stride,
288                        const uint8_t *left, const uint8_t *top)
289 {
290     uint64_t dc = 0x0101010101010101ULL *
291         ((left[0] + left[1] + left[2] + left[3] + left[4] + left[5] + left[6] +
292           left[7] + left[8] + left[9] + left[10] + left[11] + left[12] +
293           left[13] + left[14] + left[15] + left[16] + left[17] + left[18] +
294           left[19] + left[20] + left[21] + left[22] + left[23] + left[24] +
295           left[25] + left[26] + left[27] + left[28] + left[29] + left[30] +
296           left[31] + top[0] + top[1] + top[2] + top[3] + top[4] + top[5] +
297           top[6] + top[7] + top[8] + top[9] + top[10] + top[11] + top[12] +
298           top[13] + top[14] + top[15] + top[16] + top[17] + top[18] + top[19] +
299           top[20] + top[21] + top[22] + top[23] + top[24] + top[25] + top[26] +
300           top[27] + top[28] + top[29] + top[30] + top[31] + 32) >> 6);
301     int y;
302
303     for (y = 0; y < 32; y++) {
304         AV_WN64A(dst +  0, dc);
305         AV_WN64A(dst +  8, dc);
306         AV_WN64A(dst + 16, dc);
307         AV_WN64A(dst + 24, dc);
308         dst += stride;
309     }
310 }
311
312 static void dc_left_4x4_c(uint8_t *dst, ptrdiff_t stride,
313                           const uint8_t *left, const uint8_t *top)
314 {
315     unsigned dc = 0x01010101U * ((left[0] + left[1] + left[2] + left[3] + 2) >> 2);
316
317     AV_WN32A(dst + stride * 0, dc);
318     AV_WN32A(dst + stride * 1, dc);
319     AV_WN32A(dst + stride * 2, dc);
320     AV_WN32A(dst + stride * 3, dc);
321 }
322
323 static void dc_left_8x8_c(uint8_t *dst, ptrdiff_t stride,
324                           const uint8_t *left, const uint8_t *top)
325 {
326     uint64_t dc = 0x0101010101010101ULL *
327         ((left[0] + left[1] + left[2] + left[3] +
328           left[4] + left[5] + left[6] + left[7] + 4) >> 3);
329     int y;
330
331     for (y = 0; y < 8; y++) {
332         AV_WN64A(dst, dc);
333         dst += stride;
334     }
335 }
336
337 static void dc_left_16x16_c(uint8_t *dst, ptrdiff_t stride,
338                             const uint8_t *left, const uint8_t *top)
339 {
340     uint64_t dc = 0x0101010101010101ULL *
341         ((left[0] + left[1] + left[2] + left[3] + left[4] + left[5] +
342           left[6] + left[7] + left[8] + left[9] + left[10] + left[11] +
343           left[12] + left[13] + left[14] + left[15] + 8) >> 4);
344     int y;
345
346     for (y = 0; y < 16; y++) {
347         AV_WN64A(dst + 0, dc);
348         AV_WN64A(dst + 8, dc);
349         dst += stride;
350     }
351 }
352
353 static void dc_left_32x32_c(uint8_t *dst, ptrdiff_t stride,
354                             const uint8_t *left, const uint8_t *top)
355 {
356     uint64_t dc = 0x0101010101010101ULL *
357         ((left[0] + left[1] + left[2] + left[3] + left[4] + left[5] +
358           left[6] + left[7] + left[8] + left[9] + left[10] + left[11] +
359           left[12] + left[13] + left[14] + left[15] + left[16] + left[17] +
360           left[18] + left[19] + left[20] + left[21] + left[22] + left[23] +
361           left[24] + left[25] + left[26] + left[27] + left[28] + left[29] +
362           left[30] + left[31] + 16) >> 5);
363     int y;
364
365     for (y = 0; y < 32; y++) {
366         AV_WN64A(dst +  0, dc);
367         AV_WN64A(dst +  8, dc);
368         AV_WN64A(dst + 16, dc);
369         AV_WN64A(dst + 24, dc);
370         dst += stride;
371     }
372 }
373
374 static void dc_top_4x4_c(uint8_t *dst, ptrdiff_t stride,
375                          const uint8_t *left, const uint8_t *top)
376 {
377     unsigned dc = 0x01010101U * ((top[0] + top[1] + top[2] + top[3] + 2) >> 2);
378
379     AV_WN32A(dst + stride * 0, dc);
380     AV_WN32A(dst + stride * 1, dc);
381     AV_WN32A(dst + stride * 2, dc);
382     AV_WN32A(dst + stride * 3, dc);
383 }
384
385 static void dc_top_8x8_c(uint8_t *dst, ptrdiff_t stride,
386                          const uint8_t *left, const uint8_t *top)
387 {
388     uint64_t dc = 0x0101010101010101ULL *
389         ((top[0] + top[1] + top[2] + top[3] +
390           top[4] + top[5] + top[6] + top[7] + 4) >> 3);
391     int y;
392
393     for (y = 0; y < 8; y++) {
394         AV_WN64A(dst, dc);
395         dst += stride;
396     }
397 }
398
399 static void dc_top_16x16_c(uint8_t *dst, ptrdiff_t stride,
400                            const uint8_t *left, const uint8_t *top)
401 {
402     uint64_t dc = 0x0101010101010101ULL *
403         ((top[0] + top[1] + top[2] + top[3] + top[4] + top[5] +
404           top[6] + top[7] + top[8] + top[9] + top[10] + top[11] +
405           top[12] + top[13] + top[14] + top[15] + 8) >> 4);
406     int y;
407
408     for (y = 0; y < 16; y++) {
409         AV_WN64A(dst + 0, dc);
410         AV_WN64A(dst + 8, dc);
411         dst += stride;
412     }
413 }
414
415 static void dc_top_32x32_c(uint8_t *dst, ptrdiff_t stride,
416                            const uint8_t *left, const uint8_t *top)
417 {
418     uint64_t dc = 0x0101010101010101ULL *
419         ((top[0] + top[1] + top[2] + top[3] + top[4] + top[5] +
420           top[6] + top[7] + top[8] + top[9] + top[10] + top[11] +
421           top[12] + top[13] + top[14] + top[15] + top[16] + top[17] +
422           top[18] + top[19] + top[20] + top[21] + top[22] + top[23] +
423           top[24] + top[25] + top[26] + top[27] + top[28] + top[29] +
424           top[30] + top[31] + 16) >> 5);
425     int y;
426
427     for (y = 0; y < 32; y++) {
428         AV_WN64A(dst +  0, dc);
429         AV_WN64A(dst +  8, dc);
430         AV_WN64A(dst + 16, dc);
431         AV_WN64A(dst + 24, dc);
432         dst += stride;
433     }
434 }
435
436 static void dc_128_4x4_c(uint8_t *dst, ptrdiff_t stride,
437                          const uint8_t *left, const uint8_t *top)
438 {
439     AV_WN32A(dst + stride * 0, 0x80808080U);
440     AV_WN32A(dst + stride * 1, 0x80808080U);
441     AV_WN32A(dst + stride * 2, 0x80808080U);
442     AV_WN32A(dst + stride * 3, 0x80808080U);
443 }
444
445 static void dc_128_8x8_c(uint8_t *dst, ptrdiff_t stride,
446                          const uint8_t *left, const uint8_t *top)
447 {
448     int y;
449
450     for (y = 0; y < 8; y++) {
451         AV_WN64A(dst, 0x8080808080808080ULL);
452         dst += stride;
453     }
454 }
455
456 static void dc_128_16x16_c(uint8_t *dst, ptrdiff_t stride,
457                            const uint8_t *left, const uint8_t *top)
458 {
459     int y;
460
461     for (y = 0; y < 16; y++) {
462         AV_WN64A(dst + 0, 0x8080808080808080ULL);
463         AV_WN64A(dst + 8, 0x8080808080808080ULL);
464         dst += stride;
465     }
466 }
467
468 static void dc_128_32x32_c(uint8_t *dst, ptrdiff_t stride,
469                            const uint8_t *left, const uint8_t *top)
470 {
471     int y;
472
473     for (y = 0; y < 32; y++) {
474         AV_WN64A(dst +  0, 0x8080808080808080ULL);
475         AV_WN64A(dst +  8, 0x8080808080808080ULL);
476         AV_WN64A(dst + 16, 0x8080808080808080ULL);
477         AV_WN64A(dst + 24, 0x8080808080808080ULL);
478         dst += stride;
479     }
480 }
481
482 static void dc_127_4x4_c(uint8_t *dst, ptrdiff_t stride,
483                          const uint8_t *left, const uint8_t *top)
484 {
485     AV_WN32A(dst + stride * 0, 0x7F7F7F7FU);
486     AV_WN32A(dst + stride * 1, 0x7F7F7F7FU);
487     AV_WN32A(dst + stride * 2, 0x7F7F7F7FU);
488     AV_WN32A(dst + stride * 3, 0x7F7F7F7FU);
489 }
490
491 static void dc_127_8x8_c(uint8_t *dst, ptrdiff_t stride,
492                          const uint8_t *left, const uint8_t *top)
493 {
494     int y;
495
496     for (y = 0; y < 8; y++) {
497         AV_WN64A(dst, 0x7F7F7F7F7F7F7F7FULL);
498         dst += stride;
499     }
500 }
501
502 static void dc_127_16x16_c(uint8_t *dst, ptrdiff_t stride,
503                            const uint8_t *left, const uint8_t *top)
504 {
505     int y;
506
507     for (y = 0; y < 16; y++) {
508         AV_WN64A(dst + 0, 0x7F7F7F7F7F7F7F7FULL);
509         AV_WN64A(dst + 8, 0x7F7F7F7F7F7F7F7FULL);
510         dst += stride;
511     }
512 }
513
514 static void dc_127_32x32_c(uint8_t *dst, ptrdiff_t stride,
515                            const uint8_t *left, const uint8_t *top)
516 {
517     int y;
518
519     for (y = 0; y < 32; y++) {
520         AV_WN64A(dst +  0, 0x7F7F7F7F7F7F7F7FULL);
521         AV_WN64A(dst +  8, 0x7F7F7F7F7F7F7F7FULL);
522         AV_WN64A(dst + 16, 0x7F7F7F7F7F7F7F7FULL);
523         AV_WN64A(dst + 24, 0x7F7F7F7F7F7F7F7FULL);
524         dst += stride;
525     }
526 }
527
528 static void dc_129_4x4_c(uint8_t *dst, ptrdiff_t stride,
529                          const uint8_t *left, const uint8_t *top)
530 {
531     AV_WN32A(dst + stride * 0, 0x81818181U);
532     AV_WN32A(dst + stride * 1, 0x81818181U);
533     AV_WN32A(dst + stride * 2, 0x81818181U);
534     AV_WN32A(dst + stride * 3, 0x81818181U);
535 }
536
537 static void dc_129_8x8_c(uint8_t *dst, ptrdiff_t stride,
538                          const uint8_t *left, const uint8_t *top)
539 {
540     int y;
541
542     for (y = 0; y < 8; y++) {
543         AV_WN64A(dst, 0x8181818181818181ULL);
544         dst += stride;
545     }
546 }
547
548 static void dc_129_16x16_c(uint8_t *dst, ptrdiff_t stride,
549                            const uint8_t *left, const uint8_t *top)
550 {
551     int y;
552
553     for (y = 0; y < 16; y++) {
554         AV_WN64A(dst + 0, 0x8181818181818181ULL);
555         AV_WN64A(dst + 8, 0x8181818181818181ULL);
556         dst += stride;
557     }
558 }
559
560 static void dc_129_32x32_c(uint8_t *dst, ptrdiff_t stride,
561                            const uint8_t *left, const uint8_t *top)
562 {
563     int y;
564
565     for (y = 0; y < 32; y++) {
566         AV_WN64A(dst +  0, 0x8181818181818181ULL);
567         AV_WN64A(dst +  8, 0x8181818181818181ULL);
568         AV_WN64A(dst + 16, 0x8181818181818181ULL);
569         AV_WN64A(dst + 24, 0x8181818181818181ULL);
570         dst += stride;
571     }
572 }
573
574 #define DST(x, y) dst[(x) + (y) * stride]
575
576 static void diag_downleft_4x4_c(uint8_t *dst, ptrdiff_t stride,
577                                 const uint8_t *left, const uint8_t *top)
578 {
579     int a0 = top[0], a1 = top[1], a2 = top[2], a3 = top[3],
580         a4 = top[4], a5 = top[5], a6 = top[6], a7 = top[7];
581
582     DST(0,0) = (a0 + a1 * 2 + a2 + 2) >> 2;
583     DST(1,0) = DST(0,1) = (a1 + a2 * 2 + a3 + 2) >> 2;
584     DST(2,0) = DST(1,1) = DST(0,2) = (a2 + a3 * 2 + a4 + 2) >> 2;
585     DST(3,0) = DST(2,1) = DST(1,2) = DST(0,3) = (a3 + a4 * 2 + a5 + 2) >> 2;
586     DST(3,1) = DST(2,2) = DST(1,3) = (a4 + a5 * 2 + a6 + 2) >> 2;
587     DST(3,2) = DST(2,3) = (a5 + a6 * 2 + a7 + 2) >> 2;
588     DST(3,3) = a7;  // note: this is different from vp8 and such
589 }
590
591 #define def_diag_downleft(size) \
592 static void diag_downleft_##size##x##size##_c(uint8_t *dst, ptrdiff_t stride, \
593                                               const uint8_t *left, const uint8_t *top) \
594 { \
595     int i, j; \
596     uint8_t v[size - 1]; \
597 \
598     for (i = 0; i < size - 2; i++) \
599         v[i] = (top[i] + top[i + 1] * 2 + top[i + 2] + 2) >> 2; \
600     v[size - 2] = (top[size - 2] + top[size - 1] * 3 + 2) >> 2; \
601 \
602     for (j = 0; j < size; j++) { \
603         memcpy(dst + j*stride, v + j, size - 1 - j); \
604         memset(dst + j*stride + size - 1 - j, top[size - 1], j + 1); \
605     } \
606 }
607
608 def_diag_downleft(8)
609 def_diag_downleft(16)
610 def_diag_downleft(32)
611
612 static void diag_downright_4x4_c(uint8_t *dst, ptrdiff_t stride,
613                                  const uint8_t *left, const uint8_t *top)
614 {
615     int tl = top[-1], a0 = top[0], a1 = top[1], a2 = top[2], a3 = top[3],
616         l0 = left[0], l1 = left[1], l2 = left[2], l3 = left[3];
617
618     DST(0,3) = (l1 + l2 * 2 + l3 + 2) >> 2;
619     DST(0,2) = DST(1,3) = (l0 + l1 * 2 + l2 + 2) >> 2;
620     DST(0,1) = DST(1,2) = DST(2,3) = (tl + l0 * 2 + l1 + 2) >> 2;
621     DST(0,0) = DST(1,1) = DST(2,2) = DST(3,3) = (l0 + tl * 2 + a0 + 2) >> 2;
622     DST(1,0) = DST(2,1) = DST(3,2) = (tl + a0 * 2 + a1 + 2) >> 2;
623     DST(2,0) = DST(3,1) = (a0 + a1 * 2 + a2 + 2) >> 2;
624     DST(3,0) = (a1 + a2 * 2 + a3 + 2) >> 2;
625 }
626
627 #define def_diag_downright(size) \
628 static void diag_downright_##size##x##size##_c(uint8_t *dst, ptrdiff_t stride, \
629                                                const uint8_t *left, const uint8_t *top) \
630 { \
631     int i, j; \
632     uint8_t v[size + size - 1]; \
633 \
634     for (i = 0; i < size - 2; i++) { \
635         v[i           ] = (left[size - 1 - i] + left[size - 2 - i] * 2 + left[size - 3 - i] + 2) >> 2; \
636         v[size + 1 + i] = (top[i]             + top[i + 1]         * 2 + top[i + 2]         + 2) >> 2; \
637     } \
638     v[size - 2] = (left[1] + left[0] * 2 + top[-1] + 2) >> 2; \
639     v[size - 1] = (left[0] + top[-1] * 2 + top[ 0] + 2) >> 2; \
640     v[size    ] = (top[-1] + top[0]  * 2 + top[ 1] + 2) >> 2; \
641 \
642     for (j = 0; j < size; j++) \
643         memcpy(dst + j*stride, v + size - 1 - j, size); \
644 }
645
646 def_diag_downright(8)
647 def_diag_downright(16)
648 def_diag_downright(32)
649
650 static void vert_right_4x4_c(uint8_t *dst, ptrdiff_t stride,
651                              const uint8_t *left, const uint8_t *top)
652 {
653     int tl = top[-1], a0 = top[0], a1 = top[1], a2 = top[2], a3 = top[3],
654         l0 = left[0], l1 = left[1], l2 = left[2];
655
656     DST(0,3) = (l0 + l1 * 2 + l2 + 2) >> 2;
657     DST(0,2) = (tl + l0 * 2 + l1 + 2) >> 2;
658     DST(0,0) = DST(1,2) = (tl + a0 + 1) >> 1;
659     DST(0,1) = DST(1,3) = (l0 + tl * 2 + a0 + 2) >> 2;
660     DST(1,0) = DST(2,2) = (a0 + a1 + 1) >> 1;
661     DST(1,1) = DST(2,3) = (tl + a0 * 2 + a1 + 2) >> 2;
662     DST(2,0) = DST(3,2) = (a1 + a2 + 1) >> 1;
663     DST(2,1) = DST(3,3) = (a0 + a1 * 2 + a2 + 2) >> 2;
664     DST(3,0) = (a2 + a3 + 1) >> 1;
665     DST(3,1) = (a1 + a2 * 2 + a3 + 2) >> 2;
666 }
667
668 #define def_vert_right(size) \
669 static void vert_right_##size##x##size##_c(uint8_t *dst, ptrdiff_t stride, \
670                                            const uint8_t *left, const uint8_t *top) \
671 { \
672     int i, j; \
673     uint8_t ve[size + size/2 - 1], vo[size + size/2 - 1]; \
674 \
675     for (i = 0; i < size/2 - 2; i++) { \
676         vo[i] = (left[size - 4 - i*2] + left[size - 3 - i*2] * 2 + left[size - 2 - i*2] + 2) >> 2; \
677         ve[i] = (left[size - 5 - i*2] + left[size - 4 - i*2] * 2 + left[size - 3 - i*2] + 2) >> 2; \
678     } \
679     vo[size/2 - 2] = (left[0] + left[1] * 2 + left[2] + 2) >> 2; \
680     ve[size/2 - 2] = (top[-1] + left[0] * 2 + left[1] + 2) >> 2; \
681 \
682     ve[size/2 - 1] = (top[-1] + top[0] + 1) >> 1; \
683     vo[size/2 - 1] = (left[0] + top[-1] * 2 + top[0] + 2) >> 2; \
684     for (i = 0; i < size - 1; i++) { \
685         ve[size/2 + i] = (top[i] + top[i + 1] + 1) >> 1; \
686         vo[size/2 + i] = (top[i - 1] + top[i] * 2 + top[i + 1] + 2) >> 2; \
687     } \
688 \
689     for (j = 0; j < size / 2; j++) { \
690         memcpy(dst +  j*2     *stride, ve + size/2 - 1 - j, size); \
691         memcpy(dst + (j*2 + 1)*stride, vo + size/2 - 1 - j, size); \
692     } \
693 }
694
695 def_vert_right(8)
696 def_vert_right(16)
697 def_vert_right(32)
698
699 static void hor_down_4x4_c(uint8_t *dst, ptrdiff_t stride,
700                            const uint8_t *left, const uint8_t *top)
701 {
702     int l0 = left[0], l1 = left[1], l2 = left[2], l3 = left[3],
703         tl = top[-1], a0 = top[0], a1 = top[1], a2 = top[2];
704
705     DST(2,0) = (tl + a0 * 2 + a1 + 2) >> 2;
706     DST(3,0) = (a0 + a1 * 2 + a2 + 2) >> 2;
707     DST(0,0) = DST(2,1) = (tl + l0 + 1) >> 1;
708     DST(1,0) = DST(3,1) = (a0 + tl * 2 + l0 + 2) >> 2;
709     DST(0,1) = DST(2,2) = (l0 + l1 + 1) >> 1;
710     DST(1,1) = DST(3,2) = (tl + l0 * 2 + l1 + 2) >> 2;
711     DST(0,2) = DST(2,3) = (l1 + l2 + 1) >> 1;
712     DST(1,2) = DST(3,3) = (l0 + l1 * 2 + l2 + 2) >> 2;
713     DST(0,3) = (l2 + l3 + 1) >> 1;
714     DST(1,3) = (l1 + l2 * 2 + l3 + 2) >> 2;
715 }
716
717 #define def_hor_down(size) \
718 static void hor_down_##size##x##size##_c(uint8_t *dst, ptrdiff_t stride, \
719                                          const uint8_t *left, const uint8_t *top) \
720 { \
721     int i, j; \
722     uint8_t v[size * 3 - 2]; \
723 \
724     for (i = 0; i < size - 2; i++) { \
725         v[i*2       ] = (left[size - 2 - i] + left[size - 1 - i] + 1) >> 1; \
726         v[i*2    + 1] = (left[size - 3 - i] + left[size - 2 - i] * 2 + left[size - 1 - i] + 2) >> 2; \
727         v[size*2 + i] = (top[i - 1] + top[i] * 2 + top[i + 1] + 2) >> 2; \
728     } \
729     v[size*2 - 2] = (top[-1] + left[0] + 1) >> 1; \
730     v[size*2 - 4] = (left[0] + left[1] + 1) >> 1; \
731     v[size*2 - 1] = (top[0]  + top[-1] * 2 + left[0] + 2) >> 2; \
732     v[size*2 - 3] = (top[-1] + left[0] * 2 + left[1] + 2) >> 2; \
733 \
734     for (j = 0; j < size; j++) \
735         memcpy(dst + j*stride, v + size*2 - 2 - j*2, size); \
736 }
737
738 def_hor_down(8)
739 def_hor_down(16)
740 def_hor_down(32)
741
742 static void vert_left_4x4_c(uint8_t *dst, ptrdiff_t stride,
743                             const uint8_t *left, const uint8_t *top)
744 {
745     int a0 = top[0], a1 = top[1], a2 = top[2], a3 = top[3],
746         a4 = top[4], a5 = top[5], a6 = top[6];
747
748     DST(0,0) = (a0 + a1 + 1) >> 1;
749     DST(0,1) = (a0 + a1 * 2 + a2 + 2) >> 2;
750     DST(1,0) = DST(0,2) = (a1 + a2 + 1) >> 1;
751     DST(1,1) = DST(0,3) = (a1 + a2 * 2 + a3 + 2) >> 2;
752     DST(2,0) = DST(1,2) = (a2 + a3 + 1) >> 1;
753     DST(2,1) = DST(1,3) = (a2 + a3 * 2 + a4 + 2) >> 2;
754     DST(3,0) = DST(2,2) = (a3 + a4 + 1) >> 1;
755     DST(3,1) = DST(2,3) = (a3 + a4 * 2 + a5 + 2) >> 2;
756     DST(3,2) = (a4 + a5 + 1) >> 1;
757     DST(3,3) = (a4 + a5 * 2 + a6 + 2) >> 2;
758 }
759
760 #define def_vert_left(size) \
761 static void vert_left_##size##x##size##_c(uint8_t *dst, ptrdiff_t stride, \
762                                           const uint8_t *left, const uint8_t *top) \
763 { \
764     int i, j; \
765     uint8_t ve[size - 1], vo[size - 1]; \
766 \
767     for (i = 0; i < size - 2; i++) { \
768         ve[i] = (top[i] + top[i + 1] + 1) >> 1; \
769         vo[i] = (top[i] + top[i + 1] * 2 + top[i + 2] + 2) >> 2; \
770     } \
771     ve[size - 2] = (top[size - 2] + top[size - 1] + 1) >> 1; \
772     vo[size - 2] = (top[size - 2] + top[size - 1] * 3 + 2) >> 2; \
773 \
774     for (j = 0; j < size / 2; j++) { \
775         memcpy(dst +  j*2      * stride, ve + j, size - j - 1); \
776         memset(dst +  j*2      * stride + size - j - 1, top[size - 1], j + 1); \
777         memcpy(dst + (j*2 + 1) * stride, vo + j, size - j - 1); \
778         memset(dst + (j*2 + 1) * stride + size - j - 1, top[size - 1], j + 1); \
779     } \
780 }
781
782 def_vert_left(8)
783 def_vert_left(16)
784 def_vert_left(32)
785
786 static void hor_up_4x4_c(uint8_t *dst, ptrdiff_t stride,
787                          const uint8_t *left, const uint8_t *top)
788 {
789     int l0 = left[0], l1 = left[1], l2 = left[2], l3 = left[3];
790
791     DST(0,0) = (l0 + l1 + 1) >> 1;
792     DST(1,0) = (l0 + l1 * 2 + l2 + 2) >> 2;
793     DST(0,1) = DST(2,0) = (l1 + l2 + 1) >> 1;
794     DST(1,1) = DST(3,0) = (l1 + l2 * 2 + l3 + 2) >> 2;
795     DST(0,2) = DST(2,1) = (l2 + l3 + 1) >> 1;
796     DST(1,2) = DST(3,1) = (l2 + l3 * 3 + 2) >> 2;
797     DST(0,3) = DST(1,3) = DST(2,2) = DST(2,3) = DST(3,2) = DST(3,3) = l3;
798 }
799
800 #define def_hor_up(size) \
801 static void hor_up_##size##x##size##_c(uint8_t *dst, ptrdiff_t stride, \
802                                        const uint8_t *left, const uint8_t *top) \
803 { \
804     int i, j; \
805     uint8_t v[size*2 - 2]; \
806 \
807     for (i = 0; i < size - 2; i++) { \
808         v[i*2    ] = (left[i] + left[i + 1] + 1) >> 1; \
809         v[i*2 + 1] = (left[i] + left[i + 1] * 2 + left[i + 2] + 2) >> 2; \
810     } \
811     v[size*2 - 4] = (left[size - 2] + left[size - 1] + 1) >> 1; \
812     v[size*2 - 3] = (left[size - 2] + left[size - 1] * 3 + 2) >> 2; \
813 \
814     for (j = 0; j < size / 2; j++) \
815         memcpy(dst + j*stride, v + j*2, size); \
816     for (j = size / 2; j < size; j++) { \
817         memcpy(dst + j*stride, v + j*2, size*2 - 2 - j*2); \
818         memset(dst + j*stride + size*2 - 2 - j*2, left[size - 1], \
819                2 + j*2 - size); \
820     } \
821 }
822
823 def_hor_up(8)
824 def_hor_up(16)
825 def_hor_up(32)
826
827 #undef DST
828
829 static av_cold void vp9dsp_intrapred_init(VP9DSPContext *dsp)
830 {
831 #define init_intra_pred(tx, sz) \
832     dsp->intra_pred[tx][VERT_PRED]            = vert_##sz##_c; \
833     dsp->intra_pred[tx][HOR_PRED]             = hor_##sz##_c; \
834     dsp->intra_pred[tx][DC_PRED]              = dc_##sz##_c; \
835     dsp->intra_pred[tx][DIAG_DOWN_LEFT_PRED]  = diag_downleft_##sz##_c; \
836     dsp->intra_pred[tx][DIAG_DOWN_RIGHT_PRED] = diag_downright_##sz##_c; \
837     dsp->intra_pred[tx][VERT_RIGHT_PRED]      = vert_right_##sz##_c; \
838     dsp->intra_pred[tx][HOR_DOWN_PRED]        = hor_down_##sz##_c; \
839     dsp->intra_pred[tx][VERT_LEFT_PRED]       = vert_left_##sz##_c; \
840     dsp->intra_pred[tx][HOR_UP_PRED]          = hor_up_##sz##_c; \
841     dsp->intra_pred[tx][TM_VP8_PRED]          = tm_##sz##_c; \
842     dsp->intra_pred[tx][LEFT_DC_PRED]         = dc_left_##sz##_c; \
843     dsp->intra_pred[tx][TOP_DC_PRED]          = dc_top_##sz##_c; \
844     dsp->intra_pred[tx][DC_128_PRED]          = dc_128_##sz##_c; \
845     dsp->intra_pred[tx][DC_127_PRED]          = dc_127_##sz##_c; \
846     dsp->intra_pred[tx][DC_129_PRED]          = dc_129_##sz##_c
847
848     init_intra_pred(TX_4X4,   4x4);
849     init_intra_pred(TX_8X8,   8x8);
850     init_intra_pred(TX_16X16, 16x16);
851     init_intra_pred(TX_32X32, 32x32);
852
853 #undef init_intra_pred
854 }
855
856 #define itxfm_wrapper(type_a, type_b, sz, bits) \
857 static void type_a##_##type_b##_##sz##x##sz##_add_c(uint8_t *dst, \
858                                                     ptrdiff_t stride, \
859                                                     int16_t *block, int eob) \
860 { \
861     int i, j; \
862     int16_t tmp[sz * sz], out[sz]; \
863     for (i = 0; i < sz; i++) \
864         type_a##sz##_1d(block + i, sz, tmp + i * sz, 0); \
865     memset(block, 0, sz * sz * sizeof(*block)); \
866     for (i = 0; i < sz; i++) { \
867         type_b##sz##_1d(tmp + i, sz, out, 1); \
868         for (j = 0; j < sz; j++) \
869             dst[j * stride] = av_clip_uint8(dst[j * stride] + \
870                                             (bits ? \
871                                              (out[j] + (1 << (bits - 1))) >> bits : \
872                                              out[j])); \
873         dst++; \
874     } \
875 }
876
877 #define itxfm_wrap(sz, bits) \
878 itxfm_wrapper(idct,  idct,  sz, bits) \
879 itxfm_wrapper(iadst, idct,  sz, bits) \
880 itxfm_wrapper(idct,  iadst, sz, bits) \
881 itxfm_wrapper(iadst, iadst, sz, bits)
882
883 #define IN(x) in[x * stride]
884
885 static av_always_inline void idct4_1d(const int16_t *in, ptrdiff_t stride,
886                                       int16_t *out, int pass)
887 {
888     int t0, t1, t2, t3;
889
890     t0 = ((IN(0) + IN(2)) * 11585 + (1 << 13)) >> 14;
891     t1 = ((IN(0) - IN(2)) * 11585 + (1 << 13)) >> 14;
892     t2 = (IN(1) *  6270 - IN(3) * 15137 + (1 << 13)) >> 14;
893     t3 = (IN(1) * 15137 + IN(3) *  6270 + (1 << 13)) >> 14;
894
895     out[0] = t0 + t3;
896     out[1] = t1 + t2;
897     out[2] = t1 - t2;
898     out[3] = t0 - t3;
899 }
900
901 static av_always_inline void iadst4_1d(const int16_t *in, ptrdiff_t stride,
902                                        int16_t *out, int pass)
903 {
904     int t0, t1, t2, t3;
905
906     t0 =  5283 * IN(0) + 15212 * IN(2) +  9929 * IN(3);
907     t1 =  9929 * IN(0) -  5283 * IN(2) - 15212 * IN(3);
908     t2 = 13377 * (IN(0) - IN(2) + IN(3));
909     t3 = 13377 * IN(1);
910
911     out[0] = (t0 + t3      + (1 << 13)) >> 14;
912     out[1] = (t1 + t3      + (1 << 13)) >> 14;
913     out[2] = (t2           + (1 << 13)) >> 14;
914     out[3] = (t0 + t1 - t3 + (1 << 13)) >> 14;
915 }
916
917 itxfm_wrap(4, 4)
918
919 static av_always_inline void idct8_1d(const int16_t *in, ptrdiff_t stride,
920                                       int16_t *out, int pass)
921 {
922     int t0, t0a, t1, t1a, t2, t2a, t3, t3a, t4, t4a, t5, t5a, t6, t6a, t7, t7a;
923
924     t0a = ((IN(0) + IN(4)) * 11585 + (1 << 13)) >> 14;
925     t1a = ((IN(0) - IN(4)) * 11585 + (1 << 13)) >> 14;
926     t2a = (IN(2) *  6270 - IN(6) * 15137 + (1 << 13)) >> 14;
927     t3a = (IN(2) * 15137 + IN(6) *  6270 + (1 << 13)) >> 14;
928     t4a = (IN(1) *  3196 - IN(7) * 16069 + (1 << 13)) >> 14;
929     t5a = (IN(5) * 13623 - IN(3) *  9102 + (1 << 13)) >> 14;
930     t6a = (IN(5) *  9102 + IN(3) * 13623 + (1 << 13)) >> 14;
931     t7a = (IN(1) * 16069 + IN(7) *  3196 + (1 << 13)) >> 14;
932
933     t0  = t0a + t3a;
934     t1  = t1a + t2a;
935     t2  = t1a - t2a;
936     t3  = t0a - t3a;
937     t4  = t4a + t5a;
938     t5a = t4a - t5a;
939     t7  = t7a + t6a;
940     t6a = t7a - t6a;
941
942     t5  = ((t6a - t5a) * 11585 + (1 << 13)) >> 14;
943     t6  = ((t6a + t5a) * 11585 + (1 << 13)) >> 14;
944
945     out[0] = t0 + t7;
946     out[1] = t1 + t6;
947     out[2] = t2 + t5;
948     out[3] = t3 + t4;
949     out[4] = t3 - t4;
950     out[5] = t2 - t5;
951     out[6] = t1 - t6;
952     out[7] = t0 - t7;
953 }
954
955 static av_always_inline void iadst8_1d(const int16_t *in, ptrdiff_t stride,
956                                        int16_t *out, int pass)
957 {
958     int t0, t0a, t1, t1a, t2, t2a, t3, t3a, t4, t4a, t5, t5a, t6, t6a, t7, t7a;
959
960     t0a = 16305 * IN(7) +  1606 * IN(0);
961     t1a =  1606 * IN(7) - 16305 * IN(0);
962     t2a = 14449 * IN(5) +  7723 * IN(2);
963     t3a =  7723 * IN(5) - 14449 * IN(2);
964     t4a = 10394 * IN(3) + 12665 * IN(4);
965     t5a = 12665 * IN(3) - 10394 * IN(4);
966     t6a =  4756 * IN(1) + 15679 * IN(6);
967     t7a = 15679 * IN(1) -  4756 * IN(6);
968
969     t0 = (t0a + t4a + (1 << 13)) >> 14;
970     t1 = (t1a + t5a + (1 << 13)) >> 14;
971     t2 = (t2a + t6a + (1 << 13)) >> 14;
972     t3 = (t3a + t7a + (1 << 13)) >> 14;
973     t4 = (t0a - t4a + (1 << 13)) >> 14;
974     t5 = (t1a - t5a + (1 << 13)) >> 14;
975     t6 = (t2a - t6a + (1 << 13)) >> 14;
976     t7 = (t3a - t7a + (1 << 13)) >> 14;
977
978     t4a = 15137 * t4 +  6270 * t5;
979     t5a =  6270 * t4 - 15137 * t5;
980     t6a = 15137 * t7 -  6270 * t6;
981     t7a =  6270 * t7 + 15137 * t6;
982
983     out[0] =   t0 + t2;
984     out[7] = -(t1 + t3);
985     t2     =   t0 - t2;
986     t3     =   t1 - t3;
987
988     out[1] = -((t4a + t6a + (1 << 13)) >> 14);
989     out[6] =   (t5a + t7a + (1 << 13)) >> 14;
990     t6     =   (t4a - t6a + (1 << 13)) >> 14;
991     t7     =   (t5a - t7a + (1 << 13)) >> 14;
992
993     out[3] = -(((t2 + t3) * 11585 + (1 << 13)) >> 14);
994     out[4] =   ((t2 - t3) * 11585 + (1 << 13)) >> 14;
995     out[2] =   ((t6 + t7) * 11585 + (1 << 13)) >> 14;
996     out[5] = -(((t6 - t7) * 11585 + (1 << 13)) >> 14);
997 }
998
999 itxfm_wrap(8, 5)
1000
1001 static av_always_inline void idct16_1d(const int16_t *in, ptrdiff_t stride,
1002                                        int16_t *out, int pass)
1003 {
1004     int t0, t1, t2, t3, t4, t5, t6, t7, t8, t9, t10, t11, t12, t13, t14, t15;
1005     int t0a, t1a, t2a, t3a, t4a, t5a, t6a, t7a;
1006     int t8a, t9a, t10a, t11a, t12a, t13a, t14a, t15a;
1007
1008     t0a  = ((IN(0) + IN(8)) * 11585 + (1 << 13)) >> 14;
1009     t1a  = ((IN(0) - IN(8)) * 11585 + (1 << 13)) >> 14;
1010     t2a  = (IN(4)  *  6270 - IN(12) * 15137 + (1 << 13)) >> 14;
1011     t3a  = (IN(4)  * 15137 + IN(12) *  6270 + (1 << 13)) >> 14;
1012     t4a  = (IN(2)  *  3196 - IN(14) * 16069 + (1 << 13)) >> 14;
1013     t7a  = (IN(2)  * 16069 + IN(14) *  3196 + (1 << 13)) >> 14;
1014     t5a  = (IN(10) * 13623 - IN(6)  *  9102 + (1 << 13)) >> 14;
1015     t6a  = (IN(10) *  9102 + IN(6)  * 13623 + (1 << 13)) >> 14;
1016     t8a  = (IN(1)  *  1606 - IN(15) * 16305 + (1 << 13)) >> 14;
1017     t15a = (IN(1)  * 16305 + IN(15) *  1606 + (1 << 13)) >> 14;
1018     t9a  = (IN(9)  * 12665 - IN(7)  * 10394 + (1 << 13)) >> 14;
1019     t14a = (IN(9)  * 10394 + IN(7)  * 12665 + (1 << 13)) >> 14;
1020     t10a = (IN(5)  *  7723 - IN(11) * 14449 + (1 << 13)) >> 14;
1021     t13a = (IN(5)  * 14449 + IN(11) *  7723 + (1 << 13)) >> 14;
1022     t11a = (IN(13) * 15679 - IN(3)  *  4756 + (1 << 13)) >> 14;
1023     t12a = (IN(13) *  4756 + IN(3)  * 15679 + (1 << 13)) >> 14;
1024
1025     t0  = t0a  + t3a;
1026     t1  = t1a  + t2a;
1027     t2  = t1a  - t2a;
1028     t3  = t0a  - t3a;
1029     t4  = t4a  + t5a;
1030     t5  = t4a  - t5a;
1031     t6  = t7a  - t6a;
1032     t7  = t7a  + t6a;
1033     t8  = t8a  + t9a;
1034     t9  = t8a  - t9a;
1035     t10 = t11a - t10a;
1036     t11 = t11a + t10a;
1037     t12 = t12a + t13a;
1038     t13 = t12a - t13a;
1039     t14 = t15a - t14a;
1040     t15 = t15a + t14a;
1041
1042     t5a  = ((t6 - t5) * 11585 + (1 << 13)) >> 14;
1043     t6a  = ((t6 + t5) * 11585 + (1 << 13)) >> 14;
1044     t9a  = (  t14 *  6270 - t9  * 15137  + (1 << 13)) >> 14;
1045     t14a = (  t14 * 15137 + t9  *  6270  + (1 << 13)) >> 14;
1046     t10a = (-(t13 * 15137 + t10 *  6270) + (1 << 13)) >> 14;
1047     t13a = (  t13 *  6270 - t10 * 15137  + (1 << 13)) >> 14;
1048
1049     t0a  = t0   + t7;
1050     t1a  = t1   + t6a;
1051     t2a  = t2   + t5a;
1052     t3a  = t3   + t4;
1053     t4   = t3   - t4;
1054     t5   = t2   - t5a;
1055     t6   = t1   - t6a;
1056     t7   = t0   - t7;
1057     t8a  = t8   + t11;
1058     t9   = t9a  + t10a;
1059     t10  = t9a  - t10a;
1060     t11a = t8   - t11;
1061     t12a = t15  - t12;
1062     t13  = t14a - t13a;
1063     t14  = t14a + t13a;
1064     t15a = t15  + t12;
1065
1066     t10a = ((t13  - t10)  * 11585 + (1 << 13)) >> 14;
1067     t13a = ((t13  + t10)  * 11585 + (1 << 13)) >> 14;
1068     t11  = ((t12a - t11a) * 11585 + (1 << 13)) >> 14;
1069     t12  = ((t12a + t11a) * 11585 + (1 << 13)) >> 14;
1070
1071     out[ 0] = t0a + t15a;
1072     out[ 1] = t1a + t14;
1073     out[ 2] = t2a + t13a;
1074     out[ 3] = t3a + t12;
1075     out[ 4] = t4  + t11;
1076     out[ 5] = t5  + t10a;
1077     out[ 6] = t6  + t9;
1078     out[ 7] = t7  + t8a;
1079     out[ 8] = t7  - t8a;
1080     out[ 9] = t6  - t9;
1081     out[10] = t5  - t10a;
1082     out[11] = t4  - t11;
1083     out[12] = t3a - t12;
1084     out[13] = t2a - t13a;
1085     out[14] = t1a - t14;
1086     out[15] = t0a - t15a;
1087 }
1088
1089 static av_always_inline void iadst16_1d(const int16_t *in, ptrdiff_t stride,
1090                                         int16_t *out, int pass)
1091 {
1092     int t0, t1, t2, t3, t4, t5, t6, t7, t8, t9, t10, t11, t12, t13, t14, t15;
1093     int t0a, t1a, t2a, t3a, t4a, t5a, t6a, t7a;
1094     int t8a, t9a, t10a, t11a, t12a, t13a, t14a, t15a;
1095
1096     t0  = IN(15) * 16364 + IN(0)  *   804;
1097     t1  = IN(15) *   804 - IN(0)  * 16364;
1098     t2  = IN(13) * 15893 + IN(2)  *  3981;
1099     t3  = IN(13) *  3981 - IN(2)  * 15893;
1100     t4  = IN(11) * 14811 + IN(4)  *  7005;
1101     t5  = IN(11) *  7005 - IN(4)  * 14811;
1102     t6  = IN(9)  * 13160 + IN(6)  *  9760;
1103     t7  = IN(9)  *  9760 - IN(6)  * 13160;
1104     t8  = IN(7)  * 11003 + IN(8)  * 12140;
1105     t9  = IN(7)  * 12140 - IN(8)  * 11003;
1106     t10 = IN(5)  *  8423 + IN(10) * 14053;
1107     t11 = IN(5)  * 14053 - IN(10) *  8423;
1108     t12 = IN(3)  *  5520 + IN(12) * 15426;
1109     t13 = IN(3)  * 15426 - IN(12) *  5520;
1110     t14 = IN(1)  *  2404 + IN(14) * 16207;
1111     t15 = IN(1)  * 16207 - IN(14) *  2404;
1112
1113     t0a  = (t0 + t8  + (1 << 13)) >> 14;
1114     t1a  = (t1 + t9  + (1 << 13)) >> 14;
1115     t2a  = (t2 + t10 + (1 << 13)) >> 14;
1116     t3a  = (t3 + t11 + (1 << 13)) >> 14;
1117     t4a  = (t4 + t12 + (1 << 13)) >> 14;
1118     t5a  = (t5 + t13 + (1 << 13)) >> 14;
1119     t6a  = (t6 + t14 + (1 << 13)) >> 14;
1120     t7a  = (t7 + t15 + (1 << 13)) >> 14;
1121     t8a  = (t0 - t8  + (1 << 13)) >> 14;
1122     t9a  = (t1 - t9  + (1 << 13)) >> 14;
1123     t10a = (t2 - t10 + (1 << 13)) >> 14;
1124     t11a = (t3 - t11 + (1 << 13)) >> 14;
1125     t12a = (t4 - t12 + (1 << 13)) >> 14;
1126     t13a = (t5 - t13 + (1 << 13)) >> 14;
1127     t14a = (t6 - t14 + (1 << 13)) >> 14;
1128     t15a = (t7 - t15 + (1 << 13)) >> 14;
1129
1130     t8   = t8a  * 16069 + t9a  *  3196;
1131     t9   = t8a  *  3196 - t9a  * 16069;
1132     t10  = t10a *  9102 + t11a * 13623;
1133     t11  = t10a * 13623 - t11a *  9102;
1134     t12  = t13a * 16069 - t12a *  3196;
1135     t13  = t13a *  3196 + t12a * 16069;
1136     t14  = t15a *  9102 - t14a * 13623;
1137     t15  = t15a * 13623 + t14a *  9102;
1138
1139     t0   = t0a + t4a;
1140     t1   = t1a + t5a;
1141     t2   = t2a + t6a;
1142     t3   = t3a + t7a;
1143     t4   = t0a - t4a;
1144     t5   = t1a - t5a;
1145     t6   = t2a - t6a;
1146     t7   = t3a - t7a;
1147     t8a  = (t8  + t12 + (1 << 13)) >> 14;
1148     t9a  = (t9  + t13 + (1 << 13)) >> 14;
1149     t10a = (t10 + t14 + (1 << 13)) >> 14;
1150     t11a = (t11 + t15 + (1 << 13)) >> 14;
1151     t12a = (t8  - t12 + (1 << 13)) >> 14;
1152     t13a = (t9  - t13 + (1 << 13)) >> 14;
1153     t14a = (t10 - t14 + (1 << 13)) >> 14;
1154     t15a = (t11 - t15 + (1 << 13)) >> 14;
1155
1156     t4a  = t4 * 15137 + t5 *  6270;
1157     t5a  = t4 *  6270 - t5 * 15137;
1158     t6a  = t7 * 15137 - t6 *  6270;
1159     t7a  = t7 *  6270 + t6 * 15137;
1160     t12  = t12a * 15137 + t13a *  6270;
1161     t13  = t12a *  6270 - t13a * 15137;
1162     t14  = t15a * 15137 - t14a *  6270;
1163     t15  = t15a *  6270 + t14a * 15137;
1164
1165     out[ 0] =   t0 + t2;
1166     out[15] = -(t1 + t3);
1167     t2a     =   t0 - t2;
1168     t3a     =   t1 - t3;
1169     out[ 3] = -((t4a + t6a + (1 << 13)) >> 14);
1170     out[12] =   (t5a + t7a + (1 << 13)) >> 14;
1171     t6      =   (t4a - t6a + (1 << 13)) >> 14;
1172     t7      =   (t5a - t7a + (1 << 13)) >> 14;
1173     out[ 1] = -(t8a + t10a);
1174     out[14] =   t9a + t11a;
1175     t10     =   t8a - t10a;
1176     t11     =   t9a - t11a;
1177     out[ 2] =   (t12 + t14 + (1 << 13)) >> 14;
1178     out[13] = -((t13 + t15 + (1 << 13)) >> 14);
1179     t14a    =   (t12 - t14 + (1 << 13)) >> 14;
1180     t15a    =   (t13 - t15 + (1 << 13)) >> 14;
1181
1182     out[ 7] = ((t2a  + t3a)  * -11585 + (1 << 13)) >> 14;
1183     out[ 8] = ((t2a  - t3a)  *  11585 + (1 << 13)) >> 14;
1184     out[ 4] = ((t7   + t6)   *  11585 + (1 << 13)) >> 14;
1185     out[11] = ((t7   - t6)   *  11585 + (1 << 13)) >> 14;
1186     out[ 6] = ((t11  + t10)  *  11585 + (1 << 13)) >> 14;
1187     out[ 9] = ((t11  - t10)  *  11585 + (1 << 13)) >> 14;
1188     out[ 5] = ((t14a + t15a) * -11585 + (1 << 13)) >> 14;
1189     out[10] = ((t14a - t15a) *  11585 + (1 << 13)) >> 14;
1190 }
1191
1192 itxfm_wrap(16, 6)
1193
1194 static av_always_inline void idct32_1d(const int16_t *in, ptrdiff_t stride,
1195                                        int16_t *out, int pass)
1196 {
1197     int t0a  = ((IN(0) + IN(16)) * 11585 + (1 << 13)) >> 14;
1198     int t1a  = ((IN(0) - IN(16)) * 11585 + (1 << 13)) >> 14;
1199     int t2a  = (IN( 8) *  6270 - IN(24) * 15137 + (1 << 13)) >> 14;
1200     int t3a  = (IN( 8) * 15137 + IN(24) *  6270 + (1 << 13)) >> 14;
1201     int t4a  = (IN( 4) *  3196 - IN(28) * 16069 + (1 << 13)) >> 14;
1202     int t7a  = (IN( 4) * 16069 + IN(28) *  3196 + (1 << 13)) >> 14;
1203     int t5a  = (IN(20) * 13623 - IN(12) *  9102 + (1 << 13)) >> 14;
1204     int t6a  = (IN(20) *  9102 + IN(12) * 13623 + (1 << 13)) >> 14;
1205     int t8a  = (IN( 2) *  1606 - IN(30) * 16305 + (1 << 13)) >> 14;
1206     int t15a = (IN( 2) * 16305 + IN(30) *  1606 + (1 << 13)) >> 14;
1207     int t9a  = (IN(18) * 12665 - IN(14) * 10394 + (1 << 13)) >> 14;
1208     int t14a = (IN(18) * 10394 + IN(14) * 12665 + (1 << 13)) >> 14;
1209     int t10a = (IN(10) *  7723 - IN(22) * 14449 + (1 << 13)) >> 14;
1210     int t13a = (IN(10) * 14449 + IN(22) *  7723 + (1 << 13)) >> 14;
1211     int t11a = (IN(26) * 15679 - IN( 6) *  4756 + (1 << 13)) >> 14;
1212     int t12a = (IN(26) *  4756 + IN( 6) * 15679 + (1 << 13)) >> 14;
1213     int t16a = (IN( 1) *   804 - IN(31) * 16364 + (1 << 13)) >> 14;
1214     int t31a = (IN( 1) * 16364 + IN(31) *   804 + (1 << 13)) >> 14;
1215     int t17a = (IN(17) * 12140 - IN(15) * 11003 + (1 << 13)) >> 14;
1216     int t30a = (IN(17) * 11003 + IN(15) * 12140 + (1 << 13)) >> 14;
1217     int t18a = (IN( 9) *  7005 - IN(23) * 14811 + (1 << 13)) >> 14;
1218     int t29a = (IN( 9) * 14811 + IN(23) *  7005 + (1 << 13)) >> 14;
1219     int t19a = (IN(25) * 15426 - IN( 7) *  5520 + (1 << 13)) >> 14;
1220     int t28a = (IN(25) *  5520 + IN( 7) * 15426 + (1 << 13)) >> 14;
1221     int t20a = (IN( 5) *  3981 - IN(27) * 15893 + (1 << 13)) >> 14;
1222     int t27a = (IN( 5) * 15893 + IN(27) *  3981 + (1 << 13)) >> 14;
1223     int t21a = (IN(21) * 14053 - IN(11) *  8423 + (1 << 13)) >> 14;
1224     int t26a = (IN(21) *  8423 + IN(11) * 14053 + (1 << 13)) >> 14;
1225     int t22a = (IN(13) *  9760 - IN(19) * 13160 + (1 << 13)) >> 14;
1226     int t25a = (IN(13) * 13160 + IN(19) *  9760 + (1 << 13)) >> 14;
1227     int t23a = (IN(29) * 16207 - IN( 3) *  2404 + (1 << 13)) >> 14;
1228     int t24a = (IN(29) *  2404 + IN( 3) * 16207 + (1 << 13)) >> 14;
1229
1230     int t0  = t0a  + t3a;
1231     int t1  = t1a  + t2a;
1232     int t2  = t1a  - t2a;
1233     int t3  = t0a  - t3a;
1234     int t4  = t4a  + t5a;
1235     int t5  = t4a  - t5a;
1236     int t6  = t7a  - t6a;
1237     int t7  = t7a  + t6a;
1238     int t8  = t8a  + t9a;
1239     int t9  = t8a  - t9a;
1240     int t10 = t11a - t10a;
1241     int t11 = t11a + t10a;
1242     int t12 = t12a + t13a;
1243     int t13 = t12a - t13a;
1244     int t14 = t15a - t14a;
1245     int t15 = t15a + t14a;
1246     int t16 = t16a + t17a;
1247     int t17 = t16a - t17a;
1248     int t18 = t19a - t18a;
1249     int t19 = t19a + t18a;
1250     int t20 = t20a + t21a;
1251     int t21 = t20a - t21a;
1252     int t22 = t23a - t22a;
1253     int t23 = t23a + t22a;
1254     int t24 = t24a + t25a;
1255     int t25 = t24a - t25a;
1256     int t26 = t27a - t26a;
1257     int t27 = t27a + t26a;
1258     int t28 = t28a + t29a;
1259     int t29 = t28a - t29a;
1260     int t30 = t31a - t30a;
1261     int t31 = t31a + t30a;
1262
1263     t5a = ((t6 - t5) * 11585 + (1 << 13)) >> 14;
1264     t6a = ((t6 + t5) * 11585 + (1 << 13)) >> 14;
1265     t9a  = (  t14 *  6270 - t9  * 15137  + (1 << 13)) >> 14;
1266     t14a = (  t14 * 15137 + t9  *  6270  + (1 << 13)) >> 14;
1267     t10a = (-(t13 * 15137 + t10 *  6270) + (1 << 13)) >> 14;
1268     t13a = (  t13 *  6270 - t10 * 15137  + (1 << 13)) >> 14;
1269     t17a = (  t30 *  3196 - t17 * 16069  + (1 << 13)) >> 14;
1270     t30a = (  t30 * 16069 + t17 *  3196  + (1 << 13)) >> 14;
1271     t18a = (-(t29 * 16069 + t18 *  3196) + (1 << 13)) >> 14;
1272     t29a = (  t29 *  3196 - t18 * 16069  + (1 << 13)) >> 14;
1273     t21a = (  t26 * 13623 - t21 *  9102  + (1 << 13)) >> 14;
1274     t26a = (  t26 *  9102 + t21 * 13623  + (1 << 13)) >> 14;
1275     t22a = (-(t25 *  9102 + t22 * 13623) + (1 << 13)) >> 14;
1276     t25a = (  t25 * 13623 - t22 *  9102  + (1 << 13)) >> 14;
1277
1278     t0a  = t0   + t7;
1279     t1a  = t1   + t6a;
1280     t2a  = t2   + t5a;
1281     t3a  = t3   + t4;
1282     t4a  = t3   - t4;
1283     t5   = t2   - t5a;
1284     t6   = t1   - t6a;
1285     t7a  = t0   - t7;
1286     t8a  = t8   + t11;
1287     t9   = t9a  + t10a;
1288     t10  = t9a  - t10a;
1289     t11a = t8   - t11;
1290     t12a = t15  - t12;
1291     t13  = t14a - t13a;
1292     t14  = t14a + t13a;
1293     t15a = t15  + t12;
1294     t16a = t16  + t19;
1295     t17  = t17a + t18a;
1296     t18  = t17a - t18a;
1297     t19a = t16  - t19;
1298     t20a = t23  - t20;
1299     t21  = t22a - t21a;
1300     t22  = t22a + t21a;
1301     t23a = t23  + t20;
1302     t24a = t24  + t27;
1303     t25  = t25a + t26a;
1304     t26  = t25a - t26a;
1305     t27a = t24  - t27;
1306     t28a = t31  - t28;
1307     t29  = t30a - t29a;
1308     t30  = t30a + t29a;
1309     t31a = t31  + t28;
1310
1311     t10a = ((t13  - t10)  * 11585 + (1 << 13)) >> 14;
1312     t13a = ((t13  + t10)  * 11585 + (1 << 13)) >> 14;
1313     t11  = ((t12a - t11a) * 11585 + (1 << 13)) >> 14;
1314     t12  = ((t12a + t11a) * 11585 + (1 << 13)) >> 14;
1315     t18a = (  t29  *  6270 - t18  * 15137  + (1 << 13)) >> 14;
1316     t29a = (  t29  * 15137 + t18  *  6270  + (1 << 13)) >> 14;
1317     t19  = (  t28a *  6270 - t19a * 15137  + (1 << 13)) >> 14;
1318     t28  = (  t28a * 15137 + t19a *  6270  + (1 << 13)) >> 14;
1319     t20  = (-(t27a * 15137 + t20a *  6270) + (1 << 13)) >> 14;
1320     t27  = (  t27a *  6270 - t20a * 15137  + (1 << 13)) >> 14;
1321     t21a = (-(t26  * 15137 + t21  *  6270) + (1 << 13)) >> 14;
1322     t26a = (  t26  *  6270 - t21  * 15137  + (1 << 13)) >> 14;
1323
1324     t0   = t0a + t15a;
1325     t1   = t1a + t14;
1326     t2   = t2a + t13a;
1327     t3   = t3a + t12;
1328     t4   = t4a + t11;
1329     t5a  = t5  + t10a;
1330     t6a  = t6  + t9;
1331     t7   = t7a + t8a;
1332     t8   = t7a - t8a;
1333     t9a  = t6  - t9;
1334     t10  = t5  - t10a;
1335     t11a = t4a - t11;
1336     t12a = t3a - t12;
1337     t13  = t2a - t13a;
1338     t14a = t1a - t14;
1339     t15  = t0a - t15a;
1340     t16  = t16a + t23a;
1341     t17a = t17  + t22;
1342     t18  = t18a + t21a;
1343     t19a = t19  + t20;
1344     t20a = t19  - t20;
1345     t21  = t18a - t21a;
1346     t22a = t17  - t22;
1347     t23  = t16a - t23a;
1348     t24  = t31a - t24a;
1349     t25a = t30  - t25;
1350     t26  = t29a - t26a;
1351     t27a = t28  - t27;
1352     t28a = t28  + t27;
1353     t29  = t29a + t26a;
1354     t30a = t30  + t25;
1355     t31  = t31a + t24a;
1356
1357     t20  = ((t27a - t20a) * 11585 + (1 << 13)) >> 14;
1358     t27  = ((t27a + t20a) * 11585 + (1 << 13)) >> 14;
1359     t21a = ((t26  - t21 ) * 11585 + (1 << 13)) >> 14;
1360     t26a = ((t26  + t21 ) * 11585 + (1 << 13)) >> 14;
1361     t22  = ((t25a - t22a) * 11585 + (1 << 13)) >> 14;
1362     t25  = ((t25a + t22a) * 11585 + (1 << 13)) >> 14;
1363     t23a = ((t24  - t23 ) * 11585 + (1 << 13)) >> 14;
1364     t24a = ((t24  + t23 ) * 11585 + (1 << 13)) >> 14;
1365
1366     out[ 0] = t0   + t31;
1367     out[ 1] = t1   + t30a;
1368     out[ 2] = t2   + t29;
1369     out[ 3] = t3   + t28a;
1370     out[ 4] = t4   + t27;
1371     out[ 5] = t5a  + t26a;
1372     out[ 6] = t6a  + t25;
1373     out[ 7] = t7   + t24a;
1374     out[ 8] = t8   + t23a;
1375     out[ 9] = t9a  + t22;
1376     out[10] = t10  + t21a;
1377     out[11] = t11a + t20;
1378     out[12] = t12a + t19a;
1379     out[13] = t13  + t18;
1380     out[14] = t14a + t17a;
1381     out[15] = t15  + t16;
1382     out[16] = t15  - t16;
1383     out[17] = t14a - t17a;
1384     out[18] = t13  - t18;
1385     out[19] = t12a - t19a;
1386     out[20] = t11a - t20;
1387     out[21] = t10  - t21a;
1388     out[22] = t9a  - t22;
1389     out[23] = t8   - t23a;
1390     out[24] = t7   - t24a;
1391     out[25] = t6a  - t25;
1392     out[26] = t5a  - t26a;
1393     out[27] = t4   - t27;
1394     out[28] = t3   - t28a;
1395     out[29] = t2   - t29;
1396     out[30] = t1   - t30a;
1397     out[31] = t0   - t31;
1398 }
1399
1400 itxfm_wrapper(idct, idct, 32, 6)
1401
1402 static av_always_inline void iwht4_1d(const int16_t *in, ptrdiff_t stride,
1403                                       int16_t *out, int pass)
1404 {
1405     int t0, t1, t2, t3, t4;
1406
1407     if (pass == 0) {
1408         t0 = IN(0) >> 2;
1409         t1 = IN(3) >> 2;
1410         t2 = IN(1) >> 2;
1411         t3 = IN(2) >> 2;
1412     } else {
1413         t0 = IN(0);
1414         t1 = IN(3);
1415         t2 = IN(1);
1416         t3 = IN(2);
1417     }
1418
1419     t0 += t2;
1420     t3 -= t1;
1421     t4 = (t0 - t3) >> 1;
1422     t1 = t4 - t1;
1423     t2 = t4 - t2;
1424     t0 -= t1;
1425     t3 += t2;
1426
1427     out[0] = t0;
1428     out[1] = t1;
1429     out[2] = t2;
1430     out[3] = t3;
1431 }
1432
1433 itxfm_wrapper(iwht, iwht, 4, 0)
1434
1435 #undef IN
1436 #undef itxfm_wrapper
1437 #undef itxfm_wrap
1438
1439 static av_cold void vp9dsp_itxfm_init(VP9DSPContext *dsp)
1440 {
1441 #define init_itxfm(tx, sz) \
1442     dsp->itxfm_add[tx][DCT_DCT]   = idct_idct_##sz##_add_c; \
1443     dsp->itxfm_add[tx][DCT_ADST]  = iadst_idct_##sz##_add_c; \
1444     dsp->itxfm_add[tx][ADST_DCT]  = idct_iadst_##sz##_add_c; \
1445     dsp->itxfm_add[tx][ADST_ADST] = iadst_iadst_##sz##_add_c
1446
1447 #define init_idct(tx, nm) \
1448     dsp->itxfm_add[tx][DCT_DCT]   = \
1449     dsp->itxfm_add[tx][ADST_DCT]  = \
1450     dsp->itxfm_add[tx][DCT_ADST]  = \
1451     dsp->itxfm_add[tx][ADST_ADST] = nm##_add_c
1452
1453     init_itxfm(TX_4X4,   4x4);
1454     init_itxfm(TX_8X8,   8x8);
1455     init_itxfm(TX_16X16, 16x16);
1456     init_idct(TX_32X32,  idct_idct_32x32);
1457     init_idct(4 /* lossless */, iwht_iwht_4x4);
1458
1459 #undef init_itxfm
1460 #undef init_idct
1461 }
1462
1463 static av_always_inline void loop_filter(uint8_t *dst,  ptrdiff_t stride,
1464                                          int E, int I, int H,
1465                                          ptrdiff_t stridea, ptrdiff_t strideb,
1466                                          int wd)
1467 {
1468     int i;
1469
1470     for (i = 0; i < 8; i++, dst += stridea) {
1471         int p7, p6, p5, p4;
1472         int p3 = dst[strideb * -4], p2 = dst[strideb * -3];
1473         int p1 = dst[strideb * -2], p0 = dst[strideb * -1];
1474         int q0 = dst[strideb * +0], q1 = dst[strideb * +1];
1475         int q2 = dst[strideb * +2], q3 = dst[strideb * +3];
1476         int q4, q5, q6, q7;
1477         int fm = FFABS(p3 - p2) <= I && FFABS(p2 - p1) <= I &&
1478                  FFABS(p1 - p0) <= I && FFABS(q1 - q0) <= I &&
1479                  FFABS(q2 - q1) <= I && FFABS(q3 - q2) <= I &&
1480                  FFABS(p0 - q0) * 2 + (FFABS(p1 - q1) >> 1) <= E;
1481         int flat8out, flat8in;
1482
1483         if (!fm)
1484             continue;
1485
1486         if (wd >= 16) {
1487             p7 = dst[strideb * -8];
1488             p6 = dst[strideb * -7];
1489             p5 = dst[strideb * -6];
1490             p4 = dst[strideb * -5];
1491             q4 = dst[strideb * +4];
1492             q5 = dst[strideb * +5];
1493             q6 = dst[strideb * +6];
1494             q7 = dst[strideb * +7];
1495
1496             flat8out = FFABS(p7 - p0) <= 1 && FFABS(p6 - p0) <= 1 &&
1497                        FFABS(p5 - p0) <= 1 && FFABS(p4 - p0) <= 1 &&
1498                        FFABS(q4 - q0) <= 1 && FFABS(q5 - q0) <= 1 &&
1499                        FFABS(q6 - q0) <= 1 && FFABS(q7 - q0) <= 1;
1500         }
1501
1502         if (wd >= 8)
1503             flat8in = FFABS(p3 - p0) <= 1 && FFABS(p2 - p0) <= 1 &&
1504                       FFABS(p1 - p0) <= 1 && FFABS(q1 - q0) <= 1 &&
1505                       FFABS(q2 - q0) <= 1 && FFABS(q3 - q0) <= 1;
1506
1507         if (wd >= 16 && flat8out && flat8in) {
1508             dst[strideb * -7] = (p7 + p7 + p7 + p7 + p7 + p7 + p7 + p6 * 2 +
1509                                  p5 + p4 + p3 + p2 + p1 + p0 + q0 + 8) >> 4;
1510             dst[strideb * -6] = (p7 + p7 + p7 + p7 + p7 + p7 + p6 + p5 * 2 +
1511                                  p4 + p3 + p2 + p1 + p0 + q0 + q1 + 8) >> 4;
1512             dst[strideb * -5] = (p7 + p7 + p7 + p7 + p7 + p6 + p5 + p4 * 2 +
1513                                  p3 + p2 + p1 + p0 + q0 + q1 + q2 + 8) >> 4;
1514             dst[strideb * -4] = (p7 + p7 + p7 + p7 + p6 + p5 + p4 + p3 * 2 +
1515                                  p2 + p1 + p0 + q0 + q1 + q2 + q3 + 8) >> 4;
1516             dst[strideb * -3] = (p7 + p7 + p7 + p6 + p5 + p4 + p3 + p2 * 2 +
1517                                  p1 + p0 + q0 + q1 + q2 + q3 + q4 + 8) >> 4;
1518             dst[strideb * -2] = (p7 + p7 + p6 + p5 + p4 + p3 + p2 + p1 * 2 +
1519                                  p0 + q0 + q1 + q2 + q3 + q4 + q5 + 8) >> 4;
1520             dst[strideb * -1] = (p7 + p6 + p5 + p4 + p3 + p2 + p1 + p0 * 2 +
1521                                  q0 + q1 + q2 + q3 + q4 + q5 + q6 + 8) >> 4;
1522             dst[strideb * +0] = (p6 + p5 + p4 + p3 + p2 + p1 + p0 + q0 * 2 +
1523                                  q1 + q2 + q3 + q4 + q5 + q6 + q7 + 8) >> 4;
1524             dst[strideb * +1] = (p5 + p4 + p3 + p2 + p1 + p0 + q0 + q1 * 2 +
1525                                  q2 + q3 + q4 + q5 + q6 + q7 + q7 + 8) >> 4;
1526             dst[strideb * +2] = (p4 + p3 + p2 + p1 + p0 + q0 + q1 + q2 * 2 +
1527                                  q3 + q4 + q5 + q6 + q7 + q7 + q7 + 8) >> 4;
1528             dst[strideb * +3] = (p3 + p2 + p1 + p0 + q0 + q1 + q2 + q3 * 2 +
1529                                  q4 + q5 + q6 + q7 + q7 + q7 + q7 + 8) >> 4;
1530             dst[strideb * +4] = (p2 + p1 + p0 + q0 + q1 + q2 + q3 + q4 * 2 +
1531                                  q5 + q6 + q7 + q7 + q7 + q7 + q7 + 8) >> 4;
1532             dst[strideb * +5] = (p1 + p0 + q0 + q1 + q2 + q3 + q4 + q5 * 2 +
1533                                  q6 + q7 + q7 + q7 + q7 + q7 + q7 + 8) >> 4;
1534             dst[strideb * +6] = (p0 + q0 + q1 + q2 + q3 + q4 + q5 + q6 * 2 +
1535                                  q7 + q7 + q7 + q7 + q7 + q7 + q7 + 8) >> 4;
1536         } else if (wd >= 8 && flat8in) {
1537             dst[strideb * -3] = (p3 + p3 + p3 + 2 * p2 + p1 + p0 + q0 + 4) >> 3;
1538             dst[strideb * -2] = (p3 + p3 + p2 + 2 * p1 + p0 + q0 + q1 + 4) >> 3;
1539             dst[strideb * -1] = (p3 + p2 + p1 + 2 * p0 + q0 + q1 + q2 + 4) >> 3;
1540             dst[strideb * +0] = (p2 + p1 + p0 + 2 * q0 + q1 + q2 + q3 + 4) >> 3;
1541             dst[strideb * +1] = (p1 + p0 + q0 + 2 * q1 + q2 + q3 + q3 + 4) >> 3;
1542             dst[strideb * +2] = (p0 + q0 + q1 + 2 * q2 + q3 + q3 + q3 + 4) >> 3;
1543         } else {
1544             int hev = FFABS(p1 - p0) > H || FFABS(q1 - q0) > H;
1545
1546             if (hev) {
1547                 int f = av_clip_int8(3 * (q0 - p0) + av_clip_int8(p1 - q1)), f1, f2;
1548
1549                 f1 = FFMIN(f + 4, 127) >> 3;
1550                 f2 = FFMIN(f + 3, 127) >> 3;
1551
1552                 dst[strideb * -1] = av_clip_uint8(p0 + f2);
1553                 dst[strideb * +0] = av_clip_uint8(q0 - f1);
1554             } else {
1555                 int f = av_clip_int8(3 * (q0 - p0)), f1, f2;
1556
1557                 f1 = FFMIN(f + 4, 127) >> 3;
1558                 f2 = FFMIN(f + 3, 127) >> 3;
1559
1560                 dst[strideb * -1] = av_clip_uint8(p0 + f2);
1561                 dst[strideb * +0] = av_clip_uint8(q0 - f1);
1562
1563                 f = (f1 + 1) >> 1;
1564                 dst[strideb * -2] = av_clip_uint8(p1 + f);
1565                 dst[strideb * +1] = av_clip_uint8(q1 - f);
1566             }
1567         }
1568     }
1569 }
1570
1571 #define lf_8_fn(dir, wd, stridea, strideb) \
1572 static void loop_filter_##dir##_##wd##_8_c(uint8_t *dst, \
1573                                            ptrdiff_t stride, \
1574                                            int E, int I, int H) \
1575 { \
1576     loop_filter(dst, stride, E, I, H, stridea, strideb, wd); \
1577 }
1578
1579 #define lf_8_fns(wd) \
1580 lf_8_fn(h, wd, stride, 1) \
1581 lf_8_fn(v, wd, 1, stride)
1582
1583 lf_8_fns(4)
1584 lf_8_fns(8)
1585 lf_8_fns(16)
1586
1587 #undef lf_8_fn
1588 #undef lf_8_fns
1589
1590 #define lf_16_fn(dir, stridea) \
1591 static void loop_filter_##dir##_16_16_c(uint8_t *dst, \
1592                                         ptrdiff_t stride, \
1593                                         int E, int I, int H) \
1594 { \
1595     loop_filter_##dir##_16_8_c(dst, stride, E, I, H); \
1596     loop_filter_##dir##_16_8_c(dst + 8 * stridea, stride, E, I, H); \
1597 }
1598
1599 lf_16_fn(h, stride)
1600 lf_16_fn(v, 1)
1601
1602 #undef lf_16_fn
1603
1604 #define lf_mix_fn(dir, wd1, wd2, stridea) \
1605 static void loop_filter_##dir##_##wd1##wd2##_16_c(uint8_t *dst, \
1606                                                   ptrdiff_t stride, \
1607                                                   int E, int I, int H) \
1608 { \
1609     loop_filter_##dir##_##wd1##_8_c(dst, stride, E & 0xff, I & 0xff, H & 0xff); \
1610     loop_filter_##dir##_##wd2##_8_c(dst + 8 * stridea, stride, E >> 8, I >> 8, H >> 8); \
1611 }
1612
1613 #define lf_mix_fns(wd1, wd2) \
1614 lf_mix_fn(h, wd1, wd2, stride) \
1615 lf_mix_fn(v, wd1, wd2, 1)
1616
1617 lf_mix_fns(4, 4)
1618 lf_mix_fns(4, 8)
1619 lf_mix_fns(8, 4)
1620 lf_mix_fns(8, 8)
1621
1622 #undef lf_mix_fn
1623 #undef lf_mix_fns
1624
1625 static av_cold void vp9dsp_loopfilter_init(VP9DSPContext *dsp)
1626 {
1627     dsp->loop_filter_8[0][0] = loop_filter_h_4_8_c;
1628     dsp->loop_filter_8[0][1] = loop_filter_v_4_8_c;
1629     dsp->loop_filter_8[1][0] = loop_filter_h_8_8_c;
1630     dsp->loop_filter_8[1][1] = loop_filter_v_8_8_c;
1631     dsp->loop_filter_8[2][0] = loop_filter_h_16_8_c;
1632     dsp->loop_filter_8[2][1] = loop_filter_v_16_8_c;
1633
1634     dsp->loop_filter_16[0] = loop_filter_h_16_16_c;
1635     dsp->loop_filter_16[1] = loop_filter_v_16_16_c;
1636
1637     dsp->loop_filter_mix2[0][0][0] = loop_filter_h_44_16_c;
1638     dsp->loop_filter_mix2[0][0][1] = loop_filter_v_44_16_c;
1639     dsp->loop_filter_mix2[0][1][0] = loop_filter_h_48_16_c;
1640     dsp->loop_filter_mix2[0][1][1] = loop_filter_v_48_16_c;
1641     dsp->loop_filter_mix2[1][0][0] = loop_filter_h_84_16_c;
1642     dsp->loop_filter_mix2[1][0][1] = loop_filter_v_84_16_c;
1643     dsp->loop_filter_mix2[1][1][0] = loop_filter_h_88_16_c;
1644     dsp->loop_filter_mix2[1][1][1] = loop_filter_v_88_16_c;
1645 }
1646
1647 static av_always_inline void copy_c(uint8_t *dst, ptrdiff_t dst_stride,
1648                                     const uint8_t *src, ptrdiff_t src_stride,
1649                                     int w, int h)
1650 {
1651     do {
1652         memcpy(dst, src, w);
1653
1654         dst += dst_stride;
1655         src += src_stride;
1656     } while (--h);
1657 }
1658
1659 static av_always_inline void avg_c(uint8_t *dst, ptrdiff_t dst_stride,
1660                                    const uint8_t *src, ptrdiff_t src_stride,
1661                                    int w, int h)
1662 {
1663     do {
1664         int x;
1665
1666         for (x = 0; x < w; x += 4)
1667             AV_WN32A(&dst[x], rnd_avg32(AV_RN32A(&dst[x]), AV_RN32(&src[x])));
1668
1669         dst += dst_stride;
1670         src += src_stride;
1671     } while (--h);
1672 }
1673
1674 #define fpel_fn(type, sz) \
1675 static void type##sz##_c(uint8_t *dst, ptrdiff_t dst_stride, \
1676                          const uint8_t *src, ptrdiff_t src_stride, \
1677                          int h, int mx, int my) \
1678 { \
1679     type##_c(dst, dst_stride, src, src_stride, sz, h); \
1680 }
1681
1682 #define copy_avg_fn(sz) \
1683 fpel_fn(copy, sz) \
1684 fpel_fn(avg,  sz)
1685
1686 copy_avg_fn(64)
1687 copy_avg_fn(32)
1688 copy_avg_fn(16)
1689 copy_avg_fn(8)
1690 copy_avg_fn(4)
1691
1692 #undef fpel_fn
1693 #undef copy_avg_fn
1694
1695 static const int8_t vp9_subpel_filters[3][15][8] = {
1696     [FILTER_8TAP_REGULAR] = {
1697         {  0,  1,  -5, 126,   8,  -3,  1,  0 },
1698         { -1,  3, -10, 122,  18,  -6,  2,  0 },
1699         { -1,  4, -13, 118,  27,  -9,  3, -1 },
1700         { -1,  4, -16, 112,  37, -11,  4, -1 },
1701         { -1,  5, -18, 105,  48, -14,  4, -1 },
1702         { -1,  5, -19,  97,  58, -16,  5, -1 },
1703         { -1,  6, -19,  88,  68, -18,  5, -1 },
1704         { -1,  6, -19,  78,  78, -19,  6, -1 },
1705         { -1,  5, -18,  68,  88, -19,  6, -1 },
1706         { -1,  5, -16,  58,  97, -19,  5, -1 },
1707         { -1,  4, -14,  48, 105, -18,  5, -1 },
1708         { -1,  4, -11,  37, 112, -16,  4, -1 },
1709         { -1,  3,  -9,  27, 118, -13,  4, -1 },
1710         {  0,  2,  -6,  18, 122, -10,  3, -1 },
1711         {  0,  1,  -3,   8, 126,  -5,  1,  0 },
1712     }, [FILTER_8TAP_SHARP] = {
1713         { -1,  3,  -7, 127,   8,  -3,  1,  0 },
1714         { -2,  5, -13, 125,  17,  -6,  3, -1 },
1715         { -3,  7, -17, 121,  27, -10,  5, -2 },
1716         { -4,  9, -20, 115,  37, -13,  6, -2 },
1717         { -4, 10, -23, 108,  48, -16,  8, -3 },
1718         { -4, 10, -24, 100,  59, -19,  9, -3 },
1719         { -4, 11, -24,  90,  70, -21, 10, -4 },
1720         { -4, 11, -23,  80,  80, -23, 11, -4 },
1721         { -4, 10, -21,  70,  90, -24, 11, -4 },
1722         { -3,  9, -19,  59, 100, -24, 10, -4 },
1723         { -3,  8, -16,  48, 108, -23, 10, -4 },
1724         { -2,  6, -13,  37, 115, -20,  9, -4 },
1725         { -2,  5, -10,  27, 121, -17,  7, -3 },
1726         { -1,  3,  -6,  17, 125, -13,  5, -2 },
1727         {  0,  1,  -3,   8, 127,  -7,  3, -1 },
1728     }, [FILTER_8TAP_SMOOTH] = {
1729         { -3, -1,  32,  64,  38,   1, -3,  0 },
1730         { -2, -2,  29,  63,  41,   2, -3,  0 },
1731         { -2, -2,  26,  63,  43,   4, -4,  0 },
1732         { -2, -3,  24,  62,  46,   5, -4,  0 },
1733         { -2, -3,  21,  60,  49,   7, -4,  0 },
1734         { -1, -4,  18,  59,  51,   9, -4,  0 },
1735         { -1, -4,  16,  57,  53,  12, -4, -1 },
1736         { -1, -4,  14,  55,  55,  14, -4, -1 },
1737         { -1, -4,  12,  53,  57,  16, -4, -1 },
1738         {  0, -4,   9,  51,  59,  18, -4, -1 },
1739         {  0, -4,   7,  49,  60,  21, -3, -2 },
1740         {  0, -4,   5,  46,  62,  24, -3, -2 },
1741         {  0, -4,   4,  43,  63,  26, -2, -2 },
1742         {  0, -3,   2,  41,  63,  29, -2, -2 },
1743         {  0, -3,   1,  38,  64,  32, -1, -3 },
1744     }
1745 };
1746
1747 #define FILTER_8TAP(src, x, F, stride) \
1748     av_clip_uint8((F[0] * src[x + -3 * stride] + \
1749                    F[1] * src[x + -2 * stride] + \
1750                    F[2] * src[x + -1 * stride] + \
1751                    F[3] * src[x + +0 * stride] + \
1752                    F[4] * src[x + +1 * stride] + \
1753                    F[5] * src[x + +2 * stride] + \
1754                    F[6] * src[x + +3 * stride] + \
1755                    F[7] * src[x + +4 * stride] + 64) >> 7)
1756
1757 static av_always_inline void do_8tap_1d_c(uint8_t *dst, ptrdiff_t dst_stride,
1758                                           const uint8_t *src, ptrdiff_t src_stride,
1759                                           int w, int h, ptrdiff_t ds,
1760                                           const int8_t *filter, int avg)
1761 {
1762     do {
1763         int x;
1764
1765         for (x = 0; x < w; x++)
1766             if (avg) {
1767                 dst[x] = (dst[x] + FILTER_8TAP(src, x, filter, ds) + 1) >> 1;
1768             } else {
1769                 dst[x] = FILTER_8TAP(src, x, filter, ds);
1770             }
1771
1772         dst += dst_stride;
1773         src += src_stride;
1774     } while (--h);
1775 }
1776
1777 #define filter_8tap_1d_fn(opn, opa, dir, ds) \
1778 static av_noinline void opn##_8tap_1d_##dir##_c(uint8_t *dst, ptrdiff_t dst_stride, \
1779                                                 const uint8_t *src, ptrdiff_t src_stride, \
1780                                                 int w, int h, const int8_t *filter) \
1781 { \
1782     do_8tap_1d_c(dst, dst_stride, src, src_stride, w, h, ds, filter, opa); \
1783 }
1784
1785 filter_8tap_1d_fn(put, 0, v, src_stride)
1786 filter_8tap_1d_fn(put, 0, h, 1)
1787 filter_8tap_1d_fn(avg, 1, v, src_stride)
1788 filter_8tap_1d_fn(avg, 1, h, 1)
1789
1790 #undef filter_8tap_1d_fn
1791
1792 static av_always_inline void do_8tap_2d_c(uint8_t *dst, ptrdiff_t dst_stride,
1793                                           const uint8_t *src, ptrdiff_t src_stride,
1794                                           int w, int h, const int8_t *filterx,
1795                                           const int8_t *filtery, int avg)
1796 {
1797     int tmp_h = h + 7;
1798     uint8_t tmp[64 * 71], *tmp_ptr = tmp;
1799
1800     src -= src_stride * 3;
1801     do {
1802         int x;
1803
1804         for (x = 0; x < w; x++)
1805             tmp_ptr[x] = FILTER_8TAP(src, x, filterx, 1);
1806
1807         tmp_ptr += 64;
1808         src += src_stride;
1809     } while (--tmp_h);
1810
1811     tmp_ptr = tmp + 64 * 3;
1812     do {
1813         int x;
1814
1815         for (x = 0; x < w; x++)
1816             if (avg) {
1817                 dst[x] = (dst[x] + FILTER_8TAP(tmp_ptr, x, filtery, 64) + 1) >> 1;
1818             } else {
1819                 dst[x] = FILTER_8TAP(tmp_ptr, x, filtery, 64);
1820             }
1821
1822         tmp_ptr += 64;
1823         dst += dst_stride;
1824     } while (--h);
1825 }
1826
1827 #define filter_8tap_2d_fn(opn, opa) \
1828 static av_noinline void opn##_8tap_2d_hv_c(uint8_t *dst, ptrdiff_t dst_stride, \
1829                                            const uint8_t *src, ptrdiff_t src_stride, \
1830                                            int w, int h, const int8_t *filterx, \
1831                                            const int8_t *filtery) \
1832 { \
1833     do_8tap_2d_c(dst, dst_stride, src, src_stride, w, h, filterx, filtery, opa); \
1834 }
1835
1836 filter_8tap_2d_fn(put, 0)
1837 filter_8tap_2d_fn(avg, 1)
1838
1839 #undef filter_8tap_2d_fn
1840
1841 #undef FILTER_8TAP
1842
1843 #define filter_fn_1d(sz, dir, dir_m, type, type_idx, avg) \
1844 static void avg##_8tap_##type##_##sz##dir##_c(uint8_t *dst, ptrdiff_t dst_stride, \
1845                                               const uint8_t *src, ptrdiff_t src_stride, \
1846                                               int h, int mx, int my) \
1847 { \
1848     avg##_8tap_1d_##dir##_c(dst, dst_stride, src, src_stride, sz, h, \
1849                             vp9_subpel_filters[type_idx][dir_m - 1]); \
1850 }
1851
1852 #define filter_fn_2d(sz, type, type_idx, avg) \
1853 static void avg##_8tap_##type##_##sz##hv_c(uint8_t *dst, ptrdiff_t dst_stride, \
1854                                            const uint8_t *src, ptrdiff_t src_stride, \
1855                                            int h, int mx, int my) \
1856 { \
1857     avg##_8tap_2d_hv_c(dst, dst_stride, src, src_stride, sz, h, \
1858                        vp9_subpel_filters[type_idx][mx - 1], \
1859                        vp9_subpel_filters[type_idx][my - 1]); \
1860 }
1861
1862 #define FILTER_BILIN(src, x, mxy, stride) \
1863     (src[x] + ((mxy * (src[x + stride] - src[x]) + 8) >> 4))
1864
1865 static av_always_inline void do_bilin_1d_c(uint8_t *dst, ptrdiff_t dst_stride,
1866                                            const uint8_t *src, ptrdiff_t src_stride,
1867                                            int w, int h, ptrdiff_t ds, int mxy, int avg)
1868 {
1869     do {
1870         int x;
1871
1872         for (x = 0; x < w; x++)
1873             if (avg) {
1874                 dst[x] = (dst[x] + FILTER_BILIN(src, x, mxy, ds) + 1) >> 1;
1875             } else {
1876                 dst[x] = FILTER_BILIN(src, x, mxy, ds);
1877             }
1878
1879         dst += dst_stride;
1880         src += src_stride;
1881     } while (--h);
1882 }
1883
1884 #define bilin_1d_fn(opn, opa, dir, ds) \
1885 static av_noinline void opn##_bilin_1d_##dir##_c(uint8_t *dst, ptrdiff_t dst_stride, \
1886                                                  const uint8_t *src, ptrdiff_t src_stride, \
1887                                                  int w, int h, int mxy) \
1888 { \
1889     do_bilin_1d_c(dst, dst_stride, src, src_stride, w, h, ds, mxy, opa); \
1890 }
1891
1892 bilin_1d_fn(put, 0, v, src_stride)
1893 bilin_1d_fn(put, 0, h, 1)
1894 bilin_1d_fn(avg, 1, v, src_stride)
1895 bilin_1d_fn(avg, 1, h, 1)
1896
1897 #undef bilin_1d_fn
1898
1899 static av_always_inline void do_bilin_2d_c(uint8_t *dst, ptrdiff_t dst_stride,
1900                                            const uint8_t *src, ptrdiff_t src_stride,
1901                                            int w, int h, int mx, int my, int avg)
1902 {
1903     uint8_t tmp[64 * 65], *tmp_ptr = tmp;
1904     int tmp_h = h + 1;
1905
1906     do {
1907         int x;
1908
1909         for (x = 0; x < w; x++)
1910             tmp_ptr[x] = FILTER_BILIN(src, x, mx, 1);
1911
1912         tmp_ptr += 64;
1913         src += src_stride;
1914     } while (--tmp_h);
1915
1916     tmp_ptr = tmp;
1917     do {
1918         int x;
1919
1920         for (x = 0; x < w; x++)
1921             if (avg) {
1922                 dst[x] = (dst[x] + FILTER_BILIN(tmp_ptr, x, my, 64) + 1) >> 1;
1923             } else {
1924                 dst[x] = FILTER_BILIN(tmp_ptr, x, my, 64);
1925             }
1926
1927         tmp_ptr += 64;
1928         dst += dst_stride;
1929     } while (--h);
1930 }
1931
1932 #define bilin_2d_fn(opn, opa) \
1933 static av_noinline void opn##_bilin_2d_hv_c(uint8_t *dst, ptrdiff_t dst_stride, \
1934                                             const uint8_t *src, ptrdiff_t src_stride, \
1935                                             int w, int h, int mx, int my) \
1936 { \
1937     do_bilin_2d_c(dst, dst_stride, src, src_stride, w, h, mx, my, opa); \
1938 }
1939
1940 bilin_2d_fn(put, 0)
1941 bilin_2d_fn(avg, 1)
1942
1943 #undef bilin_2d_fn
1944
1945 #undef FILTER_BILIN
1946
1947 #define bilinf_fn_1d(sz, dir, dir_m, avg) \
1948 static void avg##_bilin_##sz##dir##_c(uint8_t *dst, ptrdiff_t dst_stride, \
1949                                       const uint8_t *src, ptrdiff_t src_stride, \
1950                                       int h, int mx, int my) \
1951 { \
1952     avg##_bilin_1d_##dir##_c(dst, dst_stride, src, src_stride, sz, h, dir_m); \
1953 }
1954
1955 #define bilinf_fn_2d(sz, avg) \
1956 static void avg##_bilin_##sz##hv_c(uint8_t *dst, ptrdiff_t dst_stride, \
1957                                    const uint8_t *src, ptrdiff_t src_stride, \
1958                                    int h, int mx, int my) \
1959 { \
1960     avg##_bilin_2d_hv_c(dst, dst_stride, src, src_stride, sz, h, mx, my); \
1961 }
1962
1963 #define filter_fn(sz, avg) \
1964 filter_fn_1d(sz, h, mx, regular, FILTER_8TAP_REGULAR, avg) \
1965 filter_fn_1d(sz, v, my, regular, FILTER_8TAP_REGULAR, avg) \
1966 filter_fn_2d(sz,        regular, FILTER_8TAP_REGULAR, avg) \
1967 filter_fn_1d(sz, h, mx, smooth,  FILTER_8TAP_SMOOTH,  avg) \
1968 filter_fn_1d(sz, v, my, smooth,  FILTER_8TAP_SMOOTH,  avg) \
1969 filter_fn_2d(sz,        smooth,  FILTER_8TAP_SMOOTH,  avg) \
1970 filter_fn_1d(sz, h, mx, sharp,   FILTER_8TAP_SHARP,   avg) \
1971 filter_fn_1d(sz, v, my, sharp,   FILTER_8TAP_SHARP,   avg) \
1972 filter_fn_2d(sz,        sharp,   FILTER_8TAP_SHARP,   avg) \
1973 bilinf_fn_1d(sz, h, mx,                               avg) \
1974 bilinf_fn_1d(sz, v, my,                               avg) \
1975 bilinf_fn_2d(sz,                                      avg)
1976
1977 #define filter_fn_set(avg) \
1978 filter_fn(64, avg) \
1979 filter_fn(32, avg) \
1980 filter_fn(16, avg) \
1981 filter_fn(8,  avg) \
1982 filter_fn(4,  avg)
1983
1984 filter_fn_set(put)
1985 filter_fn_set(avg)
1986
1987 #undef filter_fn
1988 #undef filter_fn_set
1989 #undef filter_fn_1d
1990 #undef filter_fn_2d
1991 #undef bilinf_fn_1d
1992 #undef bilinf_fn_2d
1993
1994 static av_cold void vp9dsp_mc_init(VP9DSPContext *dsp)
1995 {
1996 #define init_fpel(idx1, idx2, sz, type) \
1997     dsp->mc[idx1][FILTER_8TAP_SMOOTH ][idx2][0][0] = type##sz##_c; \
1998     dsp->mc[idx1][FILTER_8TAP_REGULAR][idx2][0][0] = type##sz##_c; \
1999     dsp->mc[idx1][FILTER_8TAP_SHARP  ][idx2][0][0] = type##sz##_c; \
2000     dsp->mc[idx1][FILTER_BILINEAR    ][idx2][0][0] = type##sz##_c
2001
2002 #define init_copy_avg(idx, sz) \
2003     init_fpel(idx, 0, sz, copy); \
2004     init_fpel(idx, 1, sz, avg)
2005
2006     init_copy_avg(0, 64);
2007     init_copy_avg(1, 32);
2008     init_copy_avg(2, 16);
2009     init_copy_avg(3,  8);
2010     init_copy_avg(4,  4);
2011
2012 #undef init_copy_avg
2013 #undef init_fpel
2014
2015 #define init_subpel1(idx1, idx2, idxh, idxv, sz, dir, type) \
2016     dsp->mc[idx1][FILTER_8TAP_SMOOTH ][idx2][idxh][idxv] = type##_8tap_smooth_##sz##dir##_c; \
2017     dsp->mc[idx1][FILTER_8TAP_REGULAR][idx2][idxh][idxv] = type##_8tap_regular_##sz##dir##_c; \
2018     dsp->mc[idx1][FILTER_8TAP_SHARP  ][idx2][idxh][idxv] = type##_8tap_sharp_##sz##dir##_c; \
2019     dsp->mc[idx1][FILTER_BILINEAR    ][idx2][idxh][idxv] = type##_bilin_##sz##dir##_c
2020
2021 #define init_subpel2(idx, idxh, idxv, dir, type) \
2022     init_subpel1(0, idx, idxh, idxv, 64, dir, type); \
2023     init_subpel1(1, idx, idxh, idxv, 32, dir, type); \
2024     init_subpel1(2, idx, idxh, idxv, 16, dir, type); \
2025     init_subpel1(3, idx, idxh, idxv,  8, dir, type); \
2026     init_subpel1(4, idx, idxh, idxv,  4, dir, type)
2027
2028 #define init_subpel3(idx, type) \
2029     init_subpel2(idx, 1, 1, hv, type); \
2030     init_subpel2(idx, 0, 1, v, type); \
2031     init_subpel2(idx, 1, 0, h, type)
2032
2033     init_subpel3(0, put);
2034     init_subpel3(1, avg);
2035
2036 #undef init_subpel1
2037 #undef init_subpel2
2038 #undef init_subpel3
2039 }
2040
2041 av_cold void ff_vp9dsp_init(VP9DSPContext *dsp)
2042 {
2043     vp9dsp_intrapred_init(dsp);
2044     vp9dsp_itxfm_init(dsp);
2045     vp9dsp_loopfilter_init(dsp);
2046     vp9dsp_mc_init(dsp);
2047
2048     if (ARCH_X86) ff_vp9dsp_init_x86(dsp);
2049 }