1/*
2 * Copyright (c) 2016 Alexandra Hájková
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 <string.h>
22
23#include "libavutil/intreadwrite.h"
24#include "libavutil/mem_internal.h"
25
26#include "libavcodec/hevcdsp.h"
27
28#include "checkasm.h"
29
30#define randomize_buffers(buf, size)            \
31    do {                                        \
32        int j;                                  \
33        for (j = 0; j < size; j++) {            \
34            int16_t r = rnd();                  \
35            AV_WN16A(buf + j, r >> 3);          \
36        }                                       \
37    } while (0)
38
39#define randomize_buffers2(buf, size)             \
40    do {                                          \
41        int j;                                    \
42        for (j = 0; j < size; j++)                \
43            AV_WN16A(buf + j * 2, rnd() & 0x3FF); \
44    } while (0)
45
46static void compare_add_res(int size, ptrdiff_t stride, int overflow_test)
47{
48    LOCAL_ALIGNED_32(int16_t, res0, [32 * 32]);
49    LOCAL_ALIGNED_32(int16_t, res1, [32 * 32]);
50    LOCAL_ALIGNED_32(uint8_t, dst0, [32 * 32 * 2]);
51    LOCAL_ALIGNED_32(uint8_t, dst1, [32 * 32 * 2]);
52
53    declare_func_emms(AV_CPU_FLAG_MMX, void, uint8_t *dst, int16_t *res, ptrdiff_t stride);
54
55    randomize_buffers(res0, size);
56    randomize_buffers2(dst0, size);
57    if (overflow_test)
58        res0[0] = 0x8000;
59    memcpy(res1, res0, sizeof(*res0) * size);
60    memcpy(dst1, dst0, sizeof(int16_t) * size);
61
62    call_ref(dst0, res0, stride);
63    call_new(dst1, res1, stride);
64    if (memcmp(dst0, dst1, size))
65        fail();
66    bench_new(dst1, res1, stride);
67}
68
69static void check_add_res(HEVCDSPContext h, int bit_depth)
70{
71    int i;
72
73    for (i = 2; i <= 5; i++) {
74        int block_size = 1 << i;
75        int size = block_size * block_size;
76        ptrdiff_t stride = block_size << (bit_depth > 8);
77
78        if (check_func(h.add_residual[i - 2], "hevc_add_res_%dx%d_%d", block_size, block_size, bit_depth)) {
79            compare_add_res(size, stride, 0);
80            // overflow test for res = -32768
81            compare_add_res(size, stride, 1);
82        }
83    }
84}
85
86void checkasm_check_hevc_add_res(void)
87{
88    int bit_depth;
89
90    for (bit_depth = 8; bit_depth <= 10; bit_depth++) {
91        HEVCDSPContext h;
92
93        ff_hevc_dsp_init(&h, bit_depth);
94        check_add_res(h, bit_depth);
95    }
96    report("add_residual");
97}
98