checkasm: add vp9dsp.itxfm_add tests.
[ffmpeg.git] / tests / checkasm / vp9dsp.c
1 /*
2  * Copyright (c) 2015 Ronald S. Bultje <rsbultje@gmail.com>
3  *
4  * This file is part of FFmpeg.
5  *
6  * FFmpeg is free software; you can redistribute it and/or modify
7  * it under the terms of the GNU General Public License as published by
8  * the Free Software Foundation; either version 2 of the License, or
9  * (at your option) any later version.
10  *
11  * FFmpeg is distributed in the hope that it will be useful,
12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
14  * GNU General Public License for more details.
15  *
16  * You should have received a copy of the GNU General Public License along
17  * with FFmpeg; if not, write to the Free Software Foundation, Inc.,
18  * 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
19  */
20
21 #include <math.h>
22 #include <string.h>
23 #include "checkasm.h"
24 #include "libavcodec/vp9data.h"
25 #include "libavcodec/vp9dsp.h"
26 #include "libavutil/common.h"
27 #include "libavutil/internal.h"
28 #include "libavutil/intreadwrite.h"
29 #include "libavutil/mathematics.h"
30
31 static const uint32_t pixel_mask[3] = { 0xffffffff, 0x03ff03ff, 0x0fff0fff };
32 #define SIZEOF_PIXEL ((bit_depth + 7) / 8)
33
34 #define randomize_buffers()                                        \
35     do {                                                           \
36         uint32_t mask = pixel_mask[(bit_depth - 8) >> 1];          \
37         int k;                                                     \
38         for (k = -4;  k < SIZEOF_PIXEL * FFMAX(8, size); k += 4) { \
39             uint32_t r = rnd() & mask;                             \
40             AV_WN32A(a + k, r);                                    \
41         }                                                          \
42         for (k = 0; k < size * SIZEOF_PIXEL; k += 4) {             \
43             uint32_t r = rnd() & mask;                             \
44             AV_WN32A(l + k, r);                                    \
45         }                                                          \
46     } while (0)
47
48 static void check_ipred(void)
49 {
50     LOCAL_ALIGNED_32(uint8_t, a_buf, [64 * 2]);
51     uint8_t *a = &a_buf[32 * 2];
52     LOCAL_ALIGNED_32(uint8_t, l, [32 * 2]);
53     LOCAL_ALIGNED_32(uint8_t, dst0, [32 * 32 * 2]);
54     LOCAL_ALIGNED_32(uint8_t, dst1, [32 * 32 * 2]);
55     VP9DSPContext dsp;
56     int tx, mode, bit_depth;
57     declare_func(void, uint8_t *dst, ptrdiff_t stride,
58                  const uint8_t *left, const uint8_t *top);
59     static const char *const mode_names[N_INTRA_PRED_MODES] = {
60         [VERT_PRED] = "vert",
61         [HOR_PRED] = "hor",
62         [DC_PRED] = "dc",
63         [DIAG_DOWN_LEFT_PRED] = "diag_downleft",
64         [DIAG_DOWN_RIGHT_PRED] = "diag_downright",
65         [VERT_RIGHT_PRED] = "vert_right",
66         [HOR_DOWN_PRED] = "hor_down",
67         [VERT_LEFT_PRED] = "vert_left",
68         [HOR_UP_PRED] = "hor_up",
69         [TM_VP8_PRED] = "tm",
70         [LEFT_DC_PRED] = "dc_left",
71         [TOP_DC_PRED] = "dc_top",
72         [DC_128_PRED] = "dc_128",
73         [DC_127_PRED] = "dc_127",
74         [DC_129_PRED] = "dc_129",
75     };
76
77     for (bit_depth = 8; bit_depth <= 12; bit_depth += 2) {
78         ff_vp9dsp_init(&dsp, bit_depth, 0);
79         for (tx = 0; tx < 4; tx++) {
80             int size = 4 << tx;
81
82             for (mode = 0; mode < N_INTRA_PRED_MODES; mode++) {
83                 if (check_func(dsp.intra_pred[tx][mode], "vp9_%s_%dx%d_%dbpp",
84                                mode_names[mode], size, size, bit_depth)) {
85                     randomize_buffers();
86                     call_ref(dst0, size * SIZEOF_PIXEL, l, a);
87                     call_new(dst1, size * SIZEOF_PIXEL, l, a);
88                     if (memcmp(dst0, dst1, size * size * SIZEOF_PIXEL))
89                         fail();
90                     bench_new(dst1, size * SIZEOF_PIXEL,l, a);
91                 }
92             }
93         }
94     }
95     report("ipred");
96 }
97
98 #undef randomize_buffers
99
100 #define randomize_buffers() \
101     do { \
102         uint32_t mask = pixel_mask[(bit_depth - 8) >> 1];                  \
103         for (y = 0; y < sz; y++) {                                         \
104             for (x = 0; x < sz * SIZEOF_PIXEL; x += 4) {                   \
105                 uint32_t r = rnd() & mask;                                 \
106                 AV_WN32A(dst + y * sz * SIZEOF_PIXEL + x, r);              \
107                 AV_WN32A(src + y * sz * SIZEOF_PIXEL + x, rnd() & mask);   \
108             }                                                              \
109             for (x = 0; x < sz; x++) {                                     \
110                 if (bit_depth == 8) {                                      \
111                     coef[y * sz + x] = src[y * sz + x] - dst[y * sz + x];  \
112                 } else {                                                   \
113                     ((int32_t *) coef)[y * sz + x] =                       \
114                         ((uint16_t *) src)[y * sz + x] -                   \
115                         ((uint16_t *) dst)[y * sz + x];                    \
116                 }                                                          \
117             }                                                              \
118         }                                                                  \
119     } while(0)
120
121 // wht function copied from libvpx
122 static void fwht_1d(double *out, const double *in, int sz)
123 {
124     double t0 = in[0] + in[1];
125     double t3 = in[3] - in[2];
126     double t4 = trunc((t0 - t3) * 0.5);
127     double t1 = t4 - in[1];
128     double t2 = t4 - in[2];
129
130     out[0] = t0 - t2;
131     out[1] = t2;
132     out[2] = t3 + t1;
133     out[3] = t1;
134 }
135
136 // standard DCT-II
137 static void fdct_1d(double *out, const double *in, int sz)
138 {
139     int k, n;
140
141     for (k = 0; k < sz; k++) {
142         out[k] = 0.0;
143         for (n = 0; n < sz; n++)
144             out[k] += in[n] * cos(M_PI * (2 * n + 1) * k / (sz * 2.0));
145     }
146     out[0] *= M_SQRT1_2;
147 }
148
149 // see "Towards jointly optimal spatial prediction and adaptive transform in
150 // video/image coding", by J. Han, A. Saxena, and K. Rose
151 // IEEE Proc. ICASSP, pp. 726-729, Mar. 2010.
152 static void fadst4_1d(double *out, const double *in, int sz)
153 {
154     int k, n;
155
156     for (k = 0; k < sz; k++) {
157         out[k] = 0.0;
158         for (n = 0; n < sz; n++)
159             out[k] += in[n] * sin(M_PI * (n + 1) * (2 * k + 1) / (sz * 2.0 + 1.0));
160     }
161 }
162
163 // see "A Butterfly Structured Design of The Hybrid Transform Coding Scheme",
164 // by Jingning Han, Yaowu Xu, and Debargha Mukherjee
165 // http://static.googleusercontent.com/media/research.google.com/en//pubs/archive/41418.pdf
166 static void fadst_1d(double *out, const double *in, int sz)
167 {
168     int k, n;
169
170     for (k = 0; k < sz; k++) {
171         out[k] = 0.0;
172         for (n = 0; n < sz; n++)
173             out[k] += in[n] * sin(M_PI * (2 * n + 1) * (2 * k + 1) / (sz * 4.0));
174     }
175 }
176
177 typedef void (*ftx1d_fn)(double *out, const double *in, int sz);
178 static void ftx_2d(double *out, const double *in, enum TxfmMode tx,
179                    enum TxfmType txtp, int sz)
180 {
181     static const double scaling_factors[5][4] = {
182         { 4.0, 16.0 * M_SQRT1_2 / 3.0, 16.0 * M_SQRT1_2 / 3.0, 32.0 / 9.0 },
183         { 2.0, 2.0, 2.0, 2.0 },
184         { 1.0, 1.0, 1.0, 1.0 },
185         { 0.25 },
186         { 4.0 }
187     };
188     static const ftx1d_fn ftx1d_tbl[5][4][2] = {
189         {
190             { fdct_1d, fdct_1d },
191             { fadst4_1d, fdct_1d },
192             { fdct_1d, fadst4_1d },
193             { fadst4_1d, fadst4_1d },
194         }, {
195             { fdct_1d, fdct_1d },
196             { fadst_1d, fdct_1d },
197             { fdct_1d, fadst_1d },
198             { fadst_1d, fadst_1d },
199         }, {
200             { fdct_1d, fdct_1d },
201             { fadst_1d, fdct_1d },
202             { fdct_1d, fadst_1d },
203             { fadst_1d, fadst_1d },
204         }, {
205             { fdct_1d, fdct_1d },
206         }, {
207             { fwht_1d, fwht_1d },
208         },
209     };
210     double temp[1024];
211     double scaling_factor = scaling_factors[tx][txtp];
212     int i, j;
213
214     // cols
215     for (i = 0; i < sz; ++i) {
216         double temp_out[32];
217
218         ftx1d_tbl[tx][txtp][0](temp_out, &in[i * sz], sz);
219         // scale and transpose
220         for (j = 0; j < sz; ++j)
221             temp[j * sz + i] = temp_out[j] * scaling_factor;
222     }
223
224     // rows
225     for (i = 0; i < sz; i++)
226         ftx1d_tbl[tx][txtp][1](&out[i * sz], &temp[i * sz], sz);
227 }
228
229 static void ftx(int16_t *buf, enum TxfmMode tx,
230                 enum TxfmType txtp, int sz, int bit_depth)
231 {
232     double ind[1024], outd[1024];
233     int n;
234
235     emms_c();
236     for (n = 0; n < sz * sz; n++) {
237         if (bit_depth == 8)
238             ind[n] = buf[n];
239         else
240             ind[n] = ((int32_t *) buf)[n];
241     }
242     ftx_2d(outd, ind, tx, txtp, sz);
243     for (n = 0; n < sz * sz; n++) {
244         if (bit_depth == 8)
245             buf[n] = lrint(outd[n]);
246         else
247             ((int32_t *) buf)[n] = lrint(outd[n]);
248     }
249 }
250
251 static int copy_subcoefs(int16_t *out, const int16_t *in, enum TxfmMode tx,
252                          enum TxfmType txtp, int sz, int sub, int bit_depth)
253 {
254     // copy the topleft coefficients such that the return value (being the
255     // coefficient scantable index for the eob token) guarantees that only
256     // the topleft $sub out of $sz (where $sz >= $sub) coefficients in both
257     // dimensions are non-zero. This leads to braching to specific optimized
258     // simd versions (e.g. dc-only) so that we get full asm coverage in this
259     // test
260
261     int n;
262     const int16_t *scan = vp9_scans[tx][txtp];
263     int eob;
264
265     for (n = 0; n < sz * sz; n++) {
266         int rc = scan[n], rcx = rc % sz, rcy = rc / sz;
267
268         // find eob for this sub-idct
269         if (rcx >= sub || rcy >= sub)
270             break;
271
272         // copy coef
273         if (bit_depth == 8) {
274             out[rc] = in[rc];
275         } else {
276             AV_COPY32(&out[rc * 2], &in[rc * 2]);
277         }
278     }
279
280     eob = n;
281
282     for (; n < sz * sz; n++) {
283         int rc = scan[n];
284
285         // zero
286         if (bit_depth == 8) {
287             out[rc] = 0;
288         } else {
289             AV_ZERO32(&out[rc * 2]);
290         }
291     }
292
293     return eob;
294 }
295
296 static int iszero(const int16_t *c, int sz)
297 {
298     int n;
299
300     for (n = 0; n < sz; n += 4)
301         if (AV_RN32A(&c[n]))
302             return 0;
303
304     return 1;
305 }
306
307 #define SIZEOF_COEF (2 * ((bit_depth + 7) / 8))
308
309 static void check_itxfm(void)
310 {
311     LOCAL_ALIGNED_32(uint8_t, src, [32 * 32 * 2]);
312     LOCAL_ALIGNED_32(uint8_t, dst, [32 * 32 * 2]);
313     LOCAL_ALIGNED_32(uint8_t, dst0, [32 * 32 * 2]);
314     LOCAL_ALIGNED_32(uint8_t, dst1, [32 * 32 * 2]);
315     LOCAL_ALIGNED_32(int16_t, coef, [32 * 32 * 2]);
316     LOCAL_ALIGNED_32(int16_t, subcoef0, [32 * 32 * 2]);
317     LOCAL_ALIGNED_32(int16_t, subcoef1, [32 * 32 * 2]);
318     declare_func(void, uint8_t *dst, ptrdiff_t stride, int16_t *block, int eob);
319     VP9DSPContext dsp;
320     int y, x, tx, txtp, bit_depth, sub;
321     static const char *const txtp_types[N_TXFM_TYPES] = {
322         [DCT_DCT] = "dct_dct", [DCT_ADST] = "adst_dct",
323         [ADST_DCT] = "dct_adst", [ADST_ADST] = "adst_adst"
324     };
325
326     for (bit_depth = 8; bit_depth <= 12; bit_depth += 2) {
327         ff_vp9dsp_init(&dsp, bit_depth, 0);
328
329         for (tx = TX_4X4; tx <= N_TXFM_SIZES /* 4 = lossless */; tx++) {
330             int sz = 4 << (tx & 3);
331             int n_txtps = tx < TX_32X32 ? N_TXFM_TYPES : 1;
332
333             for (txtp = 0; txtp < n_txtps; txtp++) {
334                 if (check_func(dsp.itxfm_add[tx][txtp], "vp9_inv_%s_%dx%d_add_%d",
335                                tx == 4 ? "wht_wht" : txtp_types[txtp], sz, sz,
336                                bit_depth)) {
337                     randomize_buffers();
338                     ftx(coef, tx, txtp, sz, bit_depth);
339
340                     for (sub = (txtp == 0) ? 1 : sz; sub <= sz; sub <<= 1) {
341                         int eob;
342
343                         if (sub < sz) {
344                             eob = copy_subcoefs(subcoef0, coef, tx, txtp,
345                                                 sz, sub, bit_depth);
346                         } else {
347                             eob = sz * sz;
348                             memcpy(subcoef0, coef, sz * sz * SIZEOF_COEF);
349                         }
350
351                         memcpy(dst0, dst, sz * sz * SIZEOF_PIXEL);
352                         memcpy(dst1, dst, sz * sz * SIZEOF_PIXEL);
353                         memcpy(subcoef1, subcoef0, sz * sz * SIZEOF_COEF);
354                         call_ref(dst0, sz * SIZEOF_PIXEL, subcoef0, eob);
355                         call_new(dst1, sz * SIZEOF_PIXEL, subcoef1, eob);
356                         if (memcmp(dst0, dst1, sz * sz * SIZEOF_PIXEL) ||
357                             !iszero(subcoef0, sz * sz * SIZEOF_COEF) ||
358                             !iszero(subcoef1, sz * sz * SIZEOF_COEF))
359                             fail();
360                     }
361                     bench_new(dst, sz * SIZEOF_PIXEL, coef, sz * sz);
362                 }
363             }
364         }
365     }
366     report("itxfm");
367 }
368
369 #undef randomize_buffers
370
371 #define setpx(a,b,c) \
372     do { \
373         if (SIZEOF_PIXEL == 1) { \
374             buf0[(a) + (b) * jstride] = av_clip_uint8(c); \
375         } else { \
376             ((uint16_t *)buf0)[(a) + (b) * jstride] = av_clip_uintp2(c, bit_depth); \
377         } \
378     } while (0)
379
380 // c can be an assignment and must not be put under ()
381 #define setdx(a,b,c,d) setpx(a,b,c-(d)+(rnd()%((d)*2+1)))
382 #define setsx(a,b,c,d) setdx(a,b,c,(d) << (bit_depth - 8))
383 static void randomize_loopfilter_buffers(int bidx, int lineoff, int str,
384                                          int bit_depth, int dir, const int *E,
385                                          const int *F, const int *H, const int *I,
386                                          uint8_t *buf0, uint8_t *buf1)
387 {
388     uint32_t mask = (1 << bit_depth) - 1;
389     int off = dir ? lineoff : lineoff * 16;
390     int istride = dir ? 1 : 16;
391     int jstride = dir ? str : 1;
392     int i, j;
393     for (i = 0; i < 2; i++) /* flat16 */ {
394         int idx = off + i * istride, p0, q0;
395         setpx(idx,  0, q0 = rnd() & mask);
396         setsx(idx, -1, p0 = q0, E[bidx] >> 2);
397         for (j = 1; j < 8; j++) {
398             setsx(idx, -1 - j, p0, F[bidx]);
399             setsx(idx, j, q0, F[bidx]);
400         }
401     }
402     for (i = 2; i < 4; i++) /* flat8 */ {
403         int idx = off + i * istride, p0, q0;
404         setpx(idx,  0, q0 = rnd() & mask);
405         setsx(idx, -1, p0 = q0, E[bidx] >> 2);
406         for (j = 1; j < 4; j++) {
407             setsx(idx, -1 - j, p0, F[bidx]);
408             setsx(idx, j, q0, F[bidx]);
409         }
410         for (j = 4; j < 8; j++) {
411             setpx(idx, -1 - j, rnd() & mask);
412             setpx(idx, j, rnd() & mask);
413         }
414     }
415     for (i = 4; i < 6; i++) /* regular */ {
416         int idx = off + i * istride, p2, p1, p0, q0, q1, q2;
417         setpx(idx,  0, q0 = rnd() & mask);
418         setsx(idx,  1, q1 = q0, I[bidx]);
419         setsx(idx,  2, q2 = q1, I[bidx]);
420         setsx(idx,  3, q2,      I[bidx]);
421         setsx(idx, -1, p0 = q0, E[bidx] >> 2);
422         setsx(idx, -2, p1 = p0, I[bidx]);
423         setsx(idx, -3, p2 = p1, I[bidx]);
424         setsx(idx, -4, p2,      I[bidx]);
425         for (j = 4; j < 8; j++) {
426             setpx(idx, -1 - j, rnd() & mask);
427             setpx(idx, j, rnd() & mask);
428         }
429     }
430     for (i = 6; i < 8; i++) /* off */ {
431         int idx = off + i * istride;
432         for (j = 0; j < 8; j++) {
433             setpx(idx, -1 - j, rnd() & mask);
434             setpx(idx, j, rnd() & mask);
435         }
436     }
437 }
438 #define randomize_buffers(bidx, lineoff, str) \
439         randomize_loopfilter_buffers(bidx, lineoff, str, bit_depth, dir, \
440                                      E, F, H, I, buf0, buf1)
441
442 static void check_loopfilter(void)
443 {
444     LOCAL_ALIGNED_32(uint8_t, base0, [32 + 16 * 16 * 2]);
445     LOCAL_ALIGNED_32(uint8_t, base1, [32 + 16 * 16 * 2]);
446     VP9DSPContext dsp;
447     int dir, wd, wd2, bit_depth;
448     static const char *const dir_name[2] = { "h", "v" };
449     static const int E[2] = { 20, 28 }, I[2] = { 10, 16 };
450     static const int H[2] = { 7, 11 }, F[2] = { 1, 1 };
451     declare_func(void, uint8_t *dst, ptrdiff_t stride, int E, int I, int H);
452
453     for (bit_depth = 8; bit_depth <= 12; bit_depth += 2) {
454         ff_vp9dsp_init(&dsp, bit_depth, 0);
455
456         for (dir = 0; dir < 2; dir++) {
457             int midoff = (dir ? 8 * 8 : 8) * SIZEOF_PIXEL;
458             int midoff_aligned = (dir ? 8 * 8 : 16) * SIZEOF_PIXEL;
459             uint8_t *buf0 = base0 + midoff_aligned;
460             uint8_t *buf1 = base1 + midoff_aligned;
461
462             for (wd = 0; wd < 3; wd++) {
463                 // 4/8/16wd_8px
464                 if (check_func(dsp.loop_filter_8[wd][dir],
465                                "vp9_loop_filter_%s_%d_8_%dbpp",
466                                dir_name[dir], 4 << wd, bit_depth)) {
467                     randomize_buffers(0, 0, 8);
468                     memcpy(buf1 - midoff, buf0 - midoff,
469                            16 * 8 * SIZEOF_PIXEL);
470                     call_ref(buf0, 16 * SIZEOF_PIXEL >> dir, E[0], I[0], H[0]);
471                     call_new(buf1, 16 * SIZEOF_PIXEL >> dir, E[0], I[0], H[0]);
472                     if (memcmp(buf0 - midoff, buf1 - midoff, 16 * 8 * SIZEOF_PIXEL))
473                         fail();
474                     bench_new(buf1, 16 * SIZEOF_PIXEL >> dir, E[0], I[0], H[0]);
475                 }
476             }
477
478             midoff = (dir ? 16 * 8 : 8) * SIZEOF_PIXEL;
479             midoff_aligned = (dir ? 16 * 8 : 16) * SIZEOF_PIXEL;
480
481             buf0 = base0 + midoff_aligned;
482             buf1 = base1 + midoff_aligned;
483
484             // 16wd_16px loopfilter
485             if (check_func(dsp.loop_filter_16[dir],
486                            "vp9_loop_filter_%s_16_16_%dbpp",
487                            dir_name[dir], bit_depth)) {
488                 randomize_buffers(0, 0, 16);
489                 randomize_buffers(0, 8, 16);
490                 memcpy(buf1 - midoff, buf0 - midoff, 16 * 16 * SIZEOF_PIXEL);
491                 call_ref(buf0, 16 * SIZEOF_PIXEL, E[0], I[0], H[0]);
492                 call_new(buf1, 16 * SIZEOF_PIXEL, E[0], I[0], H[0]);
493                 if (memcmp(buf0 - midoff, buf1 - midoff, 16 * 16 * SIZEOF_PIXEL))
494                     fail();
495                 bench_new(buf1, 16 * SIZEOF_PIXEL, E[0], I[0], H[0]);
496             }
497
498             for (wd = 0; wd < 2; wd++) {
499                 for (wd2 = 0; wd2 < 2; wd2++) {
500                     // mix2 loopfilter
501                     if (check_func(dsp.loop_filter_mix2[wd][wd2][dir],
502                                    "vp9_loop_filter_mix2_%s_%d%d_16_%dbpp",
503                                    dir_name[dir], 4 << wd, 4 << wd2, bit_depth)) {
504                         randomize_buffers(0, 0, 16);
505                         randomize_buffers(1, 8, 16);
506                         memcpy(buf1 - midoff, buf0 - midoff, 16 * 16 * SIZEOF_PIXEL);
507 #define M(a) (((a)[1] << 8) | (a)[0])
508                         call_ref(buf0, 16 * SIZEOF_PIXEL, M(E), M(I), M(H));
509                         call_new(buf1, 16 * SIZEOF_PIXEL, M(E), M(I), M(H));
510                         if (memcmp(buf0 - midoff, buf1 - midoff, 16 * 16 * SIZEOF_PIXEL))
511                             fail();
512                         bench_new(buf1, 16 * SIZEOF_PIXEL, M(E), M(I), M(H));
513 #undef M
514                     }
515                 }
516             }
517         }
518     }
519     report("loopfilter");
520 }
521
522 #undef setsx
523 #undef setpx
524 #undef setdx
525 #undef randomize_buffers
526
527 #define DST_BUF_SIZE (size * size * SIZEOF_PIXEL)
528 #define SRC_BUF_STRIDE 72
529 #define SRC_BUF_SIZE ((size + 7) * SRC_BUF_STRIDE * SIZEOF_PIXEL)
530 #define src (buf + 3 * SIZEOF_PIXEL * (SRC_BUF_STRIDE + 1))
531
532 #define randomize_buffers()                               \
533     do {                                                  \
534         uint32_t mask = pixel_mask[(bit_depth - 8) >> 1]; \
535         int k;                                            \
536         for (k = 0; k < SRC_BUF_SIZE; k += 4) {           \
537             uint32_t r = rnd() & mask;                    \
538             AV_WN32A(buf + k, r);                         \
539         }                                                 \
540         if (op == 1) {                                    \
541             for (k = 0; k < DST_BUF_SIZE; k += 4) {       \
542                 uint32_t r = rnd() & mask;                \
543                 AV_WN32A(dst0 + k, r);                    \
544                 AV_WN32A(dst1 + k, r);                    \
545             }                                             \
546         }                                                 \
547     } while (0)
548
549 static void check_mc(void)
550 {
551     LOCAL_ALIGNED_32(uint8_t, buf, [72 * 72 * 2]);
552     LOCAL_ALIGNED_32(uint8_t, dst0, [64 * 64 * 2]);
553     LOCAL_ALIGNED_32(uint8_t, dst1, [64 * 64 * 2]);
554     VP9DSPContext dsp;
555     int op, hsize, bit_depth, filter, dx, dy;
556     declare_func(void, uint8_t *dst, ptrdiff_t dst_stride,
557                  const uint8_t *ref, ptrdiff_t ref_stride,
558                  int h, int mx, int my);
559     static const char *const filter_names[4] = {
560         "8tap_smooth", "8tap_regular", "8tap_sharp", "bilin"
561     };
562     static const char *const subpel_names[2][2] = { { "", "h" }, { "v", "hv" } };
563     static const char *const op_names[2] = { "put", "avg" };
564     char str[256];
565
566     for (op = 0; op < 2; op++) {
567         for (bit_depth = 8; bit_depth <= 12; bit_depth += 2) {
568             ff_vp9dsp_init(&dsp, bit_depth, 0);
569             for (hsize = 0; hsize < 5; hsize++) {
570                 int size = 64 >> hsize;
571
572                 for (filter = 0; filter < 4; filter++) {
573                     for (dx = 0; dx < 2; dx++) {
574                         for (dy = 0; dy < 2; dy++) {
575                             if (dx || dy) {
576                                 snprintf(str, sizeof(str),
577                                          "%s_%s_%d%s", op_names[op],
578                                          filter_names[filter], size,
579                                          subpel_names[dy][dx]);
580                             } else {
581                                 snprintf(str, sizeof(str),
582                                          "%s%d", op_names[op], size);
583                             }
584                             if (check_func(dsp.mc[hsize][filter][op][dx][dy],
585                                            "vp9_%s_%dbpp", str, bit_depth)) {
586                                 int mx = dx ? 1 + (rnd() % 14) : 0;
587                                 int my = dy ? 1 + (rnd() % 14) : 0;
588                                 randomize_buffers();
589                                 call_ref(dst0, size * SIZEOF_PIXEL,
590                                          src, SRC_BUF_STRIDE * SIZEOF_PIXEL,
591                                          size, mx, my);
592                                 call_new(dst1, size * SIZEOF_PIXEL,
593                                          src, SRC_BUF_STRIDE * SIZEOF_PIXEL,
594                                          size, mx, my);
595                                 if (memcmp(dst0, dst1, DST_BUF_SIZE))
596                                     fail();
597
598                                 // simd implementations for each filter of subpel
599                                 // functions are identical
600                                 if (filter >= 1 && filter <= 2) continue;
601                                 // 10/12 bpp for bilin are identical
602                                 if (bit_depth == 12 && filter == 3) continue;
603
604                                 bench_new(dst1, size * SIZEOF_PIXEL,
605                                           src, SRC_BUF_STRIDE * SIZEOF_PIXEL,
606                                           size, mx, my);
607                             }
608                         }
609                     }
610                 }
611             }
612         }
613     }
614     report("mc");
615 }
616
617 void checkasm_check_vp9dsp(void)
618 {
619     check_ipred();
620     check_itxfm();
621     check_loopfilter();
622     check_mc();
623 }