1/*
2 * Copyright (c) 2016 Google Inc.
3 *
4 * This file is part of FFmpeg.
5 *
6 * FFmpeg is free software; you can redistribute it and/or
7 * modify it under the terms of the GNU Lesser General Public
8 * License as published by the Free Software Foundation; either
9 * version 2.1 of the License, or (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 GNU
14 * Lesser General Public License for more details.
15 *
16 * You should have received a copy of the GNU Lesser General Public
17 * License along with FFmpeg; if not, write to the Free Software
18 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
19 */
20
21#include <stdint.h>
22
23#include "libavutil/attributes.h"
24#include "libavutil/internal.h"
25#include "libavutil/mem_internal.h"
26#include "libavutil/aarch64/cpu.h"
27#include "libavcodec/vp9dsp.h"
28#include "vp9dsp_init.h"
29
30#define declare_fpel(type, sz)                                          \
31void ff_vp9_##type##sz##_neon(uint8_t *dst, ptrdiff_t dst_stride,       \
32                              const uint8_t *src, ptrdiff_t src_stride, \
33                              int h, int mx, int my)
34
35#define declare_copy_avg(sz) \
36    declare_fpel(copy, sz);  \
37    declare_fpel(avg , sz)
38
39#define decl_mc_func(op, filter, dir, sz)                                                \
40void ff_vp9_##op##_##filter##sz##_##dir##_neon(uint8_t *dst, ptrdiff_t dst_stride,       \
41                                               const uint8_t *src, ptrdiff_t src_stride, \
42                                               int h, int mx, int my)
43
44#define define_8tap_2d_fn(op, filter, sz)                                         \
45static void op##_##filter##sz##_hv_neon(uint8_t *dst, ptrdiff_t dst_stride,       \
46                                        const uint8_t *src, ptrdiff_t src_stride, \
47                                        int h, int mx, int my)                    \
48{                                                                                 \
49    LOCAL_ALIGNED_16(uint8_t, temp, [((1 + (sz < 64)) * sz + 8) * sz]);           \
50    /* We only need h + 7 lines, but the horizontal filter assumes an             \
51     * even number of rows, so filter h + 8 lines here. */                        \
52    ff_vp9_put_##filter##sz##_h_neon(temp, sz,                                    \
53                                     src - 3 * src_stride, src_stride,            \
54                                     h + 8, mx, 0);                               \
55    ff_vp9_##op##_##filter##sz##_v_neon(dst, dst_stride,                          \
56                                        temp + 3 * sz, sz,                        \
57                                        h, 0, my);                                \
58}
59
60#define decl_filter_funcs(op, dir, sz)  \
61    decl_mc_func(op, regular, dir, sz); \
62    decl_mc_func(op, sharp,   dir, sz); \
63    decl_mc_func(op, smooth,  dir, sz)
64
65#define decl_mc_funcs(sz)           \
66    decl_filter_funcs(put, h,  sz); \
67    decl_filter_funcs(avg, h,  sz); \
68    decl_filter_funcs(put, v,  sz); \
69    decl_filter_funcs(avg, v,  sz); \
70    decl_filter_funcs(put, hv, sz); \
71    decl_filter_funcs(avg, hv, sz)
72
73#define ff_vp9_copy32_neon ff_vp9_copy32_aarch64
74#define ff_vp9_copy64_neon ff_vp9_copy64_aarch64
75
76declare_copy_avg(64);
77declare_copy_avg(32);
78declare_copy_avg(16);
79declare_copy_avg(8);
80declare_copy_avg(4);
81
82decl_mc_funcs(64);
83decl_mc_funcs(32);
84decl_mc_funcs(16);
85decl_mc_funcs(8);
86decl_mc_funcs(4);
87
88#define define_8tap_2d_funcs(sz)        \
89    define_8tap_2d_fn(put, regular, sz) \
90    define_8tap_2d_fn(put, sharp,   sz) \
91    define_8tap_2d_fn(put, smooth,  sz) \
92    define_8tap_2d_fn(avg, regular, sz) \
93    define_8tap_2d_fn(avg, sharp,   sz) \
94    define_8tap_2d_fn(avg, smooth,  sz)
95
96define_8tap_2d_funcs(64)
97define_8tap_2d_funcs(32)
98define_8tap_2d_funcs(16)
99define_8tap_2d_funcs(8)
100define_8tap_2d_funcs(4)
101
102static av_cold void vp9dsp_mc_init_aarch64(VP9DSPContext *dsp)
103{
104    int cpu_flags = av_get_cpu_flags();
105
106#define init_fpel(idx1, idx2, sz, type, suffix)      \
107    dsp->mc[idx1][FILTER_8TAP_SMOOTH ][idx2][0][0] = \
108    dsp->mc[idx1][FILTER_8TAP_REGULAR][idx2][0][0] = \
109    dsp->mc[idx1][FILTER_8TAP_SHARP  ][idx2][0][0] = \
110    dsp->mc[idx1][FILTER_BILINEAR    ][idx2][0][0] = ff_vp9_##type##sz##suffix
111
112#define init_copy(idx, sz, suffix) \
113    init_fpel(idx, 0, sz, copy, suffix)
114
115#define init_avg(idx, sz, suffix) \
116    init_fpel(idx, 1, sz, avg,  suffix)
117
118#define init_copy_avg(idx, sz) \
119    init_copy(idx, sz, _neon); \
120    init_avg (idx, sz, _neon)
121
122    if (have_armv8(cpu_flags)) {
123        init_copy(0, 64, _aarch64);
124        init_copy(1, 32, _aarch64);
125    }
126
127    if (have_neon(cpu_flags)) {
128#define init_mc_func(idx1, idx2, op, filter, fname, dir, mx, my, sz, pfx) \
129    dsp->mc[idx1][filter][idx2][mx][my] = pfx##op##_##fname##sz##_##dir##_neon
130
131#define init_mc_funcs(idx, dir, mx, my, sz, pfx)                                   \
132    init_mc_func(idx, 0, put, FILTER_8TAP_REGULAR, regular, dir, mx, my, sz, pfx); \
133    init_mc_func(idx, 0, put, FILTER_8TAP_SHARP,   sharp,   dir, mx, my, sz, pfx); \
134    init_mc_func(idx, 0, put, FILTER_8TAP_SMOOTH,  smooth,  dir, mx, my, sz, pfx); \
135    init_mc_func(idx, 1, avg, FILTER_8TAP_REGULAR, regular, dir, mx, my, sz, pfx); \
136    init_mc_func(idx, 1, avg, FILTER_8TAP_SHARP,   sharp,   dir, mx, my, sz, pfx); \
137    init_mc_func(idx, 1, avg, FILTER_8TAP_SMOOTH,  smooth,  dir, mx, my, sz, pfx)
138
139#define init_mc_funcs_dirs(idx, sz)            \
140    init_mc_funcs(idx, h,  1, 0, sz, ff_vp9_); \
141    init_mc_funcs(idx, v,  0, 1, sz, ff_vp9_); \
142    init_mc_funcs(idx, hv, 1, 1, sz,)
143
144        init_avg(0, 64, _neon);
145        init_avg(1, 32, _neon);
146        init_copy_avg(2, 16);
147        init_copy_avg(3, 8);
148        init_copy_avg(4, 4);
149
150        init_mc_funcs_dirs(0, 64);
151        init_mc_funcs_dirs(1, 32);
152        init_mc_funcs_dirs(2, 16);
153        init_mc_funcs_dirs(3, 8);
154        init_mc_funcs_dirs(4, 4);
155    }
156}
157
158#define define_itxfm(type_a, type_b, sz)                                   \
159void ff_vp9_##type_a##_##type_b##_##sz##x##sz##_add_neon(uint8_t *_dst,    \
160                                                         ptrdiff_t stride, \
161                                                         int16_t *_block, int eob)
162
163#define define_itxfm_funcs(sz)      \
164    define_itxfm(idct,  idct,  sz); \
165    define_itxfm(iadst, idct,  sz); \
166    define_itxfm(idct,  iadst, sz); \
167    define_itxfm(iadst, iadst, sz)
168
169define_itxfm_funcs(4);
170define_itxfm_funcs(8);
171define_itxfm_funcs(16);
172define_itxfm(idct, idct, 32);
173define_itxfm(iwht, iwht, 4);
174
175
176static av_cold void vp9dsp_itxfm_init_aarch64(VP9DSPContext *dsp)
177{
178    int cpu_flags = av_get_cpu_flags();
179
180    if (have_neon(cpu_flags)) {
181#define init_itxfm(tx, sz)                                             \
182    dsp->itxfm_add[tx][DCT_DCT]   = ff_vp9_idct_idct_##sz##_add_neon;  \
183    dsp->itxfm_add[tx][DCT_ADST]  = ff_vp9_iadst_idct_##sz##_add_neon; \
184    dsp->itxfm_add[tx][ADST_DCT]  = ff_vp9_idct_iadst_##sz##_add_neon; \
185    dsp->itxfm_add[tx][ADST_ADST] = ff_vp9_iadst_iadst_##sz##_add_neon
186
187#define init_idct(tx, nm)           \
188    dsp->itxfm_add[tx][DCT_DCT]   = \
189    dsp->itxfm_add[tx][ADST_DCT]  = \
190    dsp->itxfm_add[tx][DCT_ADST]  = \
191    dsp->itxfm_add[tx][ADST_ADST] = ff_vp9_##nm##_add_neon
192
193        init_itxfm(TX_4X4, 4x4);
194        init_itxfm(TX_8X8, 8x8);
195        init_itxfm(TX_16X16, 16x16);
196        init_idct(TX_32X32, idct_idct_32x32);
197        init_idct(4, iwht_iwht_4x4);
198    }
199}
200
201#define define_loop_filter(dir, wd, len) \
202void ff_vp9_loop_filter_##dir##_##wd##_##len##_neon(uint8_t *dst, ptrdiff_t stride, int E, int I, int H)
203
204#define define_loop_filters(wd, len) \
205    define_loop_filter(h, wd, len);  \
206    define_loop_filter(v, wd, len)
207
208define_loop_filters(4, 8);
209define_loop_filters(8, 8);
210define_loop_filters(16, 8);
211
212define_loop_filters(16, 16);
213
214define_loop_filters(44, 16);
215define_loop_filters(48, 16);
216define_loop_filters(84, 16);
217define_loop_filters(88, 16);
218
219static av_cold void vp9dsp_loopfilter_init_aarch64(VP9DSPContext *dsp)
220{
221    int cpu_flags = av_get_cpu_flags();
222
223    if (have_neon(cpu_flags)) {
224        dsp->loop_filter_8[0][1] = ff_vp9_loop_filter_v_4_8_neon;
225        dsp->loop_filter_8[0][0] = ff_vp9_loop_filter_h_4_8_neon;
226        dsp->loop_filter_8[1][1] = ff_vp9_loop_filter_v_8_8_neon;
227        dsp->loop_filter_8[1][0] = ff_vp9_loop_filter_h_8_8_neon;
228        dsp->loop_filter_8[2][1] = ff_vp9_loop_filter_v_16_8_neon;
229        dsp->loop_filter_8[2][0] = ff_vp9_loop_filter_h_16_8_neon;
230
231        dsp->loop_filter_16[0] = ff_vp9_loop_filter_h_16_16_neon;
232        dsp->loop_filter_16[1] = ff_vp9_loop_filter_v_16_16_neon;
233
234        dsp->loop_filter_mix2[0][0][0] = ff_vp9_loop_filter_h_44_16_neon;
235        dsp->loop_filter_mix2[0][0][1] = ff_vp9_loop_filter_v_44_16_neon;
236        dsp->loop_filter_mix2[0][1][0] = ff_vp9_loop_filter_h_48_16_neon;
237        dsp->loop_filter_mix2[0][1][1] = ff_vp9_loop_filter_v_48_16_neon;
238        dsp->loop_filter_mix2[1][0][0] = ff_vp9_loop_filter_h_84_16_neon;
239        dsp->loop_filter_mix2[1][0][1] = ff_vp9_loop_filter_v_84_16_neon;
240        dsp->loop_filter_mix2[1][1][0] = ff_vp9_loop_filter_h_88_16_neon;
241        dsp->loop_filter_mix2[1][1][1] = ff_vp9_loop_filter_v_88_16_neon;
242    }
243}
244
245av_cold void ff_vp9dsp_init_aarch64(VP9DSPContext *dsp, int bpp)
246{
247    if (bpp == 10) {
248        ff_vp9dsp_init_10bpp_aarch64(dsp);
249        return;
250    } else if (bpp == 12) {
251        ff_vp9dsp_init_12bpp_aarch64(dsp);
252        return;
253    } else if (bpp != 8)
254        return;
255
256    vp9dsp_mc_init_aarch64(dsp);
257    vp9dsp_loopfilter_init_aarch64(dsp);
258    vp9dsp_itxfm_init_aarch64(dsp);
259}
260