Coverage Report

Created: 2022-08-24 06:17

/src/aom/av1/common/av1_inv_txfm1d.c
Line
Count
Source
1
/*
2
 * Copyright (c) 2016, Alliance for Open Media. All rights reserved
3
 *
4
 * This source code is subject to the terms of the BSD 2 Clause License and
5
 * the Alliance for Open Media Patent License 1.0. If the BSD 2 Clause License
6
 * was not distributed with this source code in the LICENSE file, you can
7
 * obtain it at www.aomedia.org/license/software. If the Alliance for Open
8
 * Media Patent License 1.0 was not distributed with this source code in the
9
 * PATENTS file, you can obtain it at www.aomedia.org/license/patent.
10
 */
11
12
#include <stdlib.h>
13
#include "av1/common/av1_inv_txfm1d.h"
14
#include "av1/common/av1_txfm.h"
15
16
void av1_idct4(const int32_t *input, int32_t *output, int8_t cos_bit,
17
1.36M
               const int8_t *stage_range) {
18
1.36M
  assert(output != input);
19
1.36M
  const int32_t size = 4;
20
1.36M
  const int32_t *cospi = cospi_arr(cos_bit);
21
22
1.36M
  int32_t stage = 0;
23
1.36M
  int32_t *bf0, *bf1;
24
1.36M
  int32_t step[4];
25
26
  // stage 0;
27
28
  // stage 1;
29
1.36M
  stage++;
30
1.36M
  bf1 = output;
31
1.36M
  bf1[0] = input[0];
32
1.36M
  bf1[1] = input[2];
33
1.36M
  bf1[2] = input[1];
34
1.36M
  bf1[3] = input[3];
35
1.36M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
36
37
  // stage 2
38
1.36M
  stage++;
39
1.36M
  bf0 = output;
40
1.36M
  bf1 = step;
41
1.36M
  bf1[0] = half_btf(cospi[32], bf0[0], cospi[32], bf0[1], cos_bit);
42
1.36M
  bf1[1] = half_btf(cospi[32], bf0[0], -cospi[32], bf0[1], cos_bit);
43
1.36M
  bf1[2] = half_btf(cospi[48], bf0[2], -cospi[16], bf0[3], cos_bit);
44
1.36M
  bf1[3] = half_btf(cospi[16], bf0[2], cospi[48], bf0[3], cos_bit);
45
1.36M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
46
47
  // stage 3
48
1.36M
  stage++;
49
1.36M
  bf0 = step;
50
1.36M
  bf1 = output;
51
1.36M
  bf1[0] = clamp_value(bf0[0] + bf0[3], stage_range[stage]);
52
1.36M
  bf1[1] = clamp_value(bf0[1] + bf0[2], stage_range[stage]);
53
1.36M
  bf1[2] = clamp_value(bf0[1] - bf0[2], stage_range[stage]);
54
1.36M
  bf1[3] = clamp_value(bf0[0] - bf0[3], stage_range[stage]);
55
1.36M
}
56
57
void av1_idct8(const int32_t *input, int32_t *output, int8_t cos_bit,
58
3.41M
               const int8_t *stage_range) {
59
3.41M
  assert(output != input);
60
3.41M
  const int32_t size = 8;
61
3.41M
  const int32_t *cospi = cospi_arr(cos_bit);
62
63
3.41M
  int32_t stage = 0;
64
3.41M
  int32_t *bf0, *bf1;
65
3.41M
  int32_t step[8];
66
67
  // stage 0;
68
69
  // stage 1;
70
3.41M
  stage++;
71
3.41M
  bf1 = output;
72
3.41M
  bf1[0] = input[0];
73
3.41M
  bf1[1] = input[4];
74
3.41M
  bf1[2] = input[2];
75
3.41M
  bf1[3] = input[6];
76
3.41M
  bf1[4] = input[1];
77
3.41M
  bf1[5] = input[5];
78
3.41M
  bf1[6] = input[3];
79
3.41M
  bf1[7] = input[7];
80
3.41M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
81
82
  // stage 2
83
3.41M
  stage++;
84
3.41M
  bf0 = output;
85
3.41M
  bf1 = step;
86
3.41M
  bf1[0] = bf0[0];
87
3.41M
  bf1[1] = bf0[1];
88
3.41M
  bf1[2] = bf0[2];
89
3.41M
  bf1[3] = bf0[3];
90
3.41M
  bf1[4] = half_btf(cospi[56], bf0[4], -cospi[8], bf0[7], cos_bit);
91
3.41M
  bf1[5] = half_btf(cospi[24], bf0[5], -cospi[40], bf0[6], cos_bit);
92
3.41M
  bf1[6] = half_btf(cospi[40], bf0[5], cospi[24], bf0[6], cos_bit);
93
3.41M
  bf1[7] = half_btf(cospi[8], bf0[4], cospi[56], bf0[7], cos_bit);
94
3.41M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
95
96
  // stage 3
97
3.41M
  stage++;
98
3.41M
  bf0 = step;
99
3.41M
  bf1 = output;
100
3.41M
  bf1[0] = half_btf(cospi[32], bf0[0], cospi[32], bf0[1], cos_bit);
101
3.41M
  bf1[1] = half_btf(cospi[32], bf0[0], -cospi[32], bf0[1], cos_bit);
102
3.41M
  bf1[2] = half_btf(cospi[48], bf0[2], -cospi[16], bf0[3], cos_bit);
103
3.41M
  bf1[3] = half_btf(cospi[16], bf0[2], cospi[48], bf0[3], cos_bit);
104
3.41M
  bf1[4] = clamp_value(bf0[4] + bf0[5], stage_range[stage]);
105
3.41M
  bf1[5] = clamp_value(bf0[4] - bf0[5], stage_range[stage]);
106
3.41M
  bf1[6] = clamp_value(-bf0[6] + bf0[7], stage_range[stage]);
107
3.41M
  bf1[7] = clamp_value(bf0[6] + bf0[7], stage_range[stage]);
108
3.41M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
109
110
  // stage 4
111
3.41M
  stage++;
112
3.41M
  bf0 = output;
113
3.41M
  bf1 = step;
114
3.41M
  bf1[0] = clamp_value(bf0[0] + bf0[3], stage_range[stage]);
115
3.41M
  bf1[1] = clamp_value(bf0[1] + bf0[2], stage_range[stage]);
116
3.41M
  bf1[2] = clamp_value(bf0[1] - bf0[2], stage_range[stage]);
117
3.41M
  bf1[3] = clamp_value(bf0[0] - bf0[3], stage_range[stage]);
118
3.41M
  bf1[4] = bf0[4];
119
3.41M
  bf1[5] = half_btf(-cospi[32], bf0[5], cospi[32], bf0[6], cos_bit);
120
3.41M
  bf1[6] = half_btf(cospi[32], bf0[5], cospi[32], bf0[6], cos_bit);
121
3.41M
  bf1[7] = bf0[7];
122
3.41M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
123
124
  // stage 5
125
3.41M
  stage++;
126
3.41M
  bf0 = step;
127
3.41M
  bf1 = output;
128
3.41M
  bf1[0] = clamp_value(bf0[0] + bf0[7], stage_range[stage]);
129
3.41M
  bf1[1] = clamp_value(bf0[1] + bf0[6], stage_range[stage]);
130
3.41M
  bf1[2] = clamp_value(bf0[2] + bf0[5], stage_range[stage]);
131
3.41M
  bf1[3] = clamp_value(bf0[3] + bf0[4], stage_range[stage]);
132
3.41M
  bf1[4] = clamp_value(bf0[3] - bf0[4], stage_range[stage]);
133
3.41M
  bf1[5] = clamp_value(bf0[2] - bf0[5], stage_range[stage]);
134
3.41M
  bf1[6] = clamp_value(bf0[1] - bf0[6], stage_range[stage]);
135
3.41M
  bf1[7] = clamp_value(bf0[0] - bf0[7], stage_range[stage]);
136
3.41M
}
137
138
void av1_idct16(const int32_t *input, int32_t *output, int8_t cos_bit,
139
5.97M
                const int8_t *stage_range) {
140
5.97M
  assert(output != input);
141
5.97M
  const int32_t size = 16;
142
5.97M
  const int32_t *cospi = cospi_arr(cos_bit);
143
144
5.97M
  int32_t stage = 0;
145
5.97M
  int32_t *bf0, *bf1;
146
5.97M
  int32_t step[16];
147
148
  // stage 0;
149
150
  // stage 1;
151
5.97M
  stage++;
152
5.97M
  bf1 = output;
153
5.97M
  bf1[0] = input[0];
154
5.97M
  bf1[1] = input[8];
155
5.97M
  bf1[2] = input[4];
156
5.97M
  bf1[3] = input[12];
157
5.97M
  bf1[4] = input[2];
158
5.97M
  bf1[5] = input[10];
159
5.97M
  bf1[6] = input[6];
160
5.97M
  bf1[7] = input[14];
161
5.97M
  bf1[8] = input[1];
162
5.97M
  bf1[9] = input[9];
163
5.97M
  bf1[10] = input[5];
164
5.97M
  bf1[11] = input[13];
165
5.97M
  bf1[12] = input[3];
166
5.97M
  bf1[13] = input[11];
167
5.97M
  bf1[14] = input[7];
168
5.97M
  bf1[15] = input[15];
169
5.97M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
170
171
  // stage 2
172
5.97M
  stage++;
173
5.97M
  bf0 = output;
174
5.97M
  bf1 = step;
175
5.97M
  bf1[0] = bf0[0];
176
5.97M
  bf1[1] = bf0[1];
177
5.97M
  bf1[2] = bf0[2];
178
5.97M
  bf1[3] = bf0[3];
179
5.97M
  bf1[4] = bf0[4];
180
5.97M
  bf1[5] = bf0[5];
181
5.97M
  bf1[6] = bf0[6];
182
5.97M
  bf1[7] = bf0[7];
183
5.97M
  bf1[8] = half_btf(cospi[60], bf0[8], -cospi[4], bf0[15], cos_bit);
184
5.97M
  bf1[9] = half_btf(cospi[28], bf0[9], -cospi[36], bf0[14], cos_bit);
185
5.97M
  bf1[10] = half_btf(cospi[44], bf0[10], -cospi[20], bf0[13], cos_bit);
186
5.97M
  bf1[11] = half_btf(cospi[12], bf0[11], -cospi[52], bf0[12], cos_bit);
187
5.97M
  bf1[12] = half_btf(cospi[52], bf0[11], cospi[12], bf0[12], cos_bit);
188
5.97M
  bf1[13] = half_btf(cospi[20], bf0[10], cospi[44], bf0[13], cos_bit);
189
5.97M
  bf1[14] = half_btf(cospi[36], bf0[9], cospi[28], bf0[14], cos_bit);
190
5.97M
  bf1[15] = half_btf(cospi[4], bf0[8], cospi[60], bf0[15], cos_bit);
191
5.97M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
192
193
  // stage 3
194
5.97M
  stage++;
195
5.97M
  bf0 = step;
196
5.97M
  bf1 = output;
197
5.97M
  bf1[0] = bf0[0];
198
5.97M
  bf1[1] = bf0[1];
199
5.97M
  bf1[2] = bf0[2];
200
5.97M
  bf1[3] = bf0[3];
201
5.97M
  bf1[4] = half_btf(cospi[56], bf0[4], -cospi[8], bf0[7], cos_bit);
202
5.97M
  bf1[5] = half_btf(cospi[24], bf0[5], -cospi[40], bf0[6], cos_bit);
203
5.97M
  bf1[6] = half_btf(cospi[40], bf0[5], cospi[24], bf0[6], cos_bit);
204
5.97M
  bf1[7] = half_btf(cospi[8], bf0[4], cospi[56], bf0[7], cos_bit);
205
5.97M
  bf1[8] = clamp_value(bf0[8] + bf0[9], stage_range[stage]);
206
5.97M
  bf1[9] = clamp_value(bf0[8] - bf0[9], stage_range[stage]);
207
5.97M
  bf1[10] = clamp_value(-bf0[10] + bf0[11], stage_range[stage]);
208
5.97M
  bf1[11] = clamp_value(bf0[10] + bf0[11], stage_range[stage]);
209
5.97M
  bf1[12] = clamp_value(bf0[12] + bf0[13], stage_range[stage]);
210
5.97M
  bf1[13] = clamp_value(bf0[12] - bf0[13], stage_range[stage]);
211
5.97M
  bf1[14] = clamp_value(-bf0[14] + bf0[15], stage_range[stage]);
212
5.97M
  bf1[15] = clamp_value(bf0[14] + bf0[15], stage_range[stage]);
213
5.97M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
214
215
  // stage 4
216
5.97M
  stage++;
217
5.97M
  bf0 = output;
218
5.97M
  bf1 = step;
219
5.97M
  bf1[0] = half_btf(cospi[32], bf0[0], cospi[32], bf0[1], cos_bit);
220
5.97M
  bf1[1] = half_btf(cospi[32], bf0[0], -cospi[32], bf0[1], cos_bit);
221
5.97M
  bf1[2] = half_btf(cospi[48], bf0[2], -cospi[16], bf0[3], cos_bit);
222
5.97M
  bf1[3] = half_btf(cospi[16], bf0[2], cospi[48], bf0[3], cos_bit);
223
5.97M
  bf1[4] = clamp_value(bf0[4] + bf0[5], stage_range[stage]);
224
5.97M
  bf1[5] = clamp_value(bf0[4] - bf0[5], stage_range[stage]);
225
5.97M
  bf1[6] = clamp_value(-bf0[6] + bf0[7], stage_range[stage]);
226
5.97M
  bf1[7] = clamp_value(bf0[6] + bf0[7], stage_range[stage]);
227
5.97M
  bf1[8] = bf0[8];
228
5.97M
  bf1[9] = half_btf(-cospi[16], bf0[9], cospi[48], bf0[14], cos_bit);
229
5.97M
  bf1[10] = half_btf(-cospi[48], bf0[10], -cospi[16], bf0[13], cos_bit);
230
5.97M
  bf1[11] = bf0[11];
231
5.97M
  bf1[12] = bf0[12];
232
5.97M
  bf1[13] = half_btf(-cospi[16], bf0[10], cospi[48], bf0[13], cos_bit);
233
5.97M
  bf1[14] = half_btf(cospi[48], bf0[9], cospi[16], bf0[14], cos_bit);
234
5.97M
  bf1[15] = bf0[15];
235
5.97M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
236
237
  // stage 5
238
5.97M
  stage++;
239
5.97M
  bf0 = step;
240
5.97M
  bf1 = output;
241
5.97M
  bf1[0] = clamp_value(bf0[0] + bf0[3], stage_range[stage]);
242
5.97M
  bf1[1] = clamp_value(bf0[1] + bf0[2], stage_range[stage]);
243
5.97M
  bf1[2] = clamp_value(bf0[1] - bf0[2], stage_range[stage]);
244
5.97M
  bf1[3] = clamp_value(bf0[0] - bf0[3], stage_range[stage]);
245
5.97M
  bf1[4] = bf0[4];
246
5.97M
  bf1[5] = half_btf(-cospi[32], bf0[5], cospi[32], bf0[6], cos_bit);
247
5.97M
  bf1[6] = half_btf(cospi[32], bf0[5], cospi[32], bf0[6], cos_bit);
248
5.97M
  bf1[7] = bf0[7];
249
5.97M
  bf1[8] = clamp_value(bf0[8] + bf0[11], stage_range[stage]);
250
5.97M
  bf1[9] = clamp_value(bf0[9] + bf0[10], stage_range[stage]);
251
5.97M
  bf1[10] = clamp_value(bf0[9] - bf0[10], stage_range[stage]);
252
5.97M
  bf1[11] = clamp_value(bf0[8] - bf0[11], stage_range[stage]);
253
5.97M
  bf1[12] = clamp_value(-bf0[12] + bf0[15], stage_range[stage]);
254
5.97M
  bf1[13] = clamp_value(-bf0[13] + bf0[14], stage_range[stage]);
255
5.97M
  bf1[14] = clamp_value(bf0[13] + bf0[14], stage_range[stage]);
256
5.97M
  bf1[15] = clamp_value(bf0[12] + bf0[15], stage_range[stage]);
257
5.97M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
258
259
  // stage 6
260
5.97M
  stage++;
261
5.97M
  bf0 = output;
262
5.97M
  bf1 = step;
263
5.97M
  bf1[0] = clamp_value(bf0[0] + bf0[7], stage_range[stage]);
264
5.97M
  bf1[1] = clamp_value(bf0[1] + bf0[6], stage_range[stage]);
265
5.97M
  bf1[2] = clamp_value(bf0[2] + bf0[5], stage_range[stage]);
266
5.97M
  bf1[3] = clamp_value(bf0[3] + bf0[4], stage_range[stage]);
267
5.97M
  bf1[4] = clamp_value(bf0[3] - bf0[4], stage_range[stage]);
268
5.97M
  bf1[5] = clamp_value(bf0[2] - bf0[5], stage_range[stage]);
269
5.97M
  bf1[6] = clamp_value(bf0[1] - bf0[6], stage_range[stage]);
270
5.97M
  bf1[7] = clamp_value(bf0[0] - bf0[7], stage_range[stage]);
271
5.97M
  bf1[8] = bf0[8];
272
5.97M
  bf1[9] = bf0[9];
273
5.97M
  bf1[10] = half_btf(-cospi[32], bf0[10], cospi[32], bf0[13], cos_bit);
274
5.97M
  bf1[11] = half_btf(-cospi[32], bf0[11], cospi[32], bf0[12], cos_bit);
275
5.97M
  bf1[12] = half_btf(cospi[32], bf0[11], cospi[32], bf0[12], cos_bit);
276
5.97M
  bf1[13] = half_btf(cospi[32], bf0[10], cospi[32], bf0[13], cos_bit);
277
5.97M
  bf1[14] = bf0[14];
278
5.97M
  bf1[15] = bf0[15];
279
5.97M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
280
281
  // stage 7
282
5.97M
  stage++;
283
5.97M
  bf0 = step;
284
5.97M
  bf1 = output;
285
5.97M
  bf1[0] = clamp_value(bf0[0] + bf0[15], stage_range[stage]);
286
5.97M
  bf1[1] = clamp_value(bf0[1] + bf0[14], stage_range[stage]);
287
5.97M
  bf1[2] = clamp_value(bf0[2] + bf0[13], stage_range[stage]);
288
5.97M
  bf1[3] = clamp_value(bf0[3] + bf0[12], stage_range[stage]);
289
5.97M
  bf1[4] = clamp_value(bf0[4] + bf0[11], stage_range[stage]);
290
5.97M
  bf1[5] = clamp_value(bf0[5] + bf0[10], stage_range[stage]);
291
5.97M
  bf1[6] = clamp_value(bf0[6] + bf0[9], stage_range[stage]);
292
5.97M
  bf1[7] = clamp_value(bf0[7] + bf0[8], stage_range[stage]);
293
5.97M
  bf1[8] = clamp_value(bf0[7] - bf0[8], stage_range[stage]);
294
5.97M
  bf1[9] = clamp_value(bf0[6] - bf0[9], stage_range[stage]);
295
5.97M
  bf1[10] = clamp_value(bf0[5] - bf0[10], stage_range[stage]);
296
5.97M
  bf1[11] = clamp_value(bf0[4] - bf0[11], stage_range[stage]);
297
5.97M
  bf1[12] = clamp_value(bf0[3] - bf0[12], stage_range[stage]);
298
5.97M
  bf1[13] = clamp_value(bf0[2] - bf0[13], stage_range[stage]);
299
5.97M
  bf1[14] = clamp_value(bf0[1] - bf0[14], stage_range[stage]);
300
5.97M
  bf1[15] = clamp_value(bf0[0] - bf0[15], stage_range[stage]);
301
5.97M
}
302
303
void av1_idct32(const int32_t *input, int32_t *output, int8_t cos_bit,
304
7.37M
                const int8_t *stage_range) {
305
7.37M
  assert(output != input);
306
7.37M
  const int32_t size = 32;
307
7.37M
  const int32_t *cospi = cospi_arr(cos_bit);
308
309
7.37M
  int32_t stage = 0;
310
7.37M
  int32_t *bf0, *bf1;
311
7.37M
  int32_t step[32];
312
313
  // stage 0;
314
315
  // stage 1;
316
7.37M
  stage++;
317
7.37M
  bf1 = output;
318
7.37M
  bf1[0] = input[0];
319
7.37M
  bf1[1] = input[16];
320
7.37M
  bf1[2] = input[8];
321
7.37M
  bf1[3] = input[24];
322
7.37M
  bf1[4] = input[4];
323
7.37M
  bf1[5] = input[20];
324
7.37M
  bf1[6] = input[12];
325
7.37M
  bf1[7] = input[28];
326
7.37M
  bf1[8] = input[2];
327
7.37M
  bf1[9] = input[18];
328
7.37M
  bf1[10] = input[10];
329
7.37M
  bf1[11] = input[26];
330
7.37M
  bf1[12] = input[6];
331
7.37M
  bf1[13] = input[22];
332
7.37M
  bf1[14] = input[14];
333
7.37M
  bf1[15] = input[30];
334
7.37M
  bf1[16] = input[1];
335
7.37M
  bf1[17] = input[17];
336
7.37M
  bf1[18] = input[9];
337
7.37M
  bf1[19] = input[25];
338
7.37M
  bf1[20] = input[5];
339
7.37M
  bf1[21] = input[21];
340
7.37M
  bf1[22] = input[13];
341
7.37M
  bf1[23] = input[29];
342
7.37M
  bf1[24] = input[3];
343
7.37M
  bf1[25] = input[19];
344
7.37M
  bf1[26] = input[11];
345
7.37M
  bf1[27] = input[27];
346
7.37M
  bf1[28] = input[7];
347
7.37M
  bf1[29] = input[23];
348
7.37M
  bf1[30] = input[15];
349
7.37M
  bf1[31] = input[31];
350
7.37M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
351
352
  // stage 2
353
7.37M
  stage++;
354
7.37M
  bf0 = output;
355
7.37M
  bf1 = step;
356
7.37M
  bf1[0] = bf0[0];
357
7.37M
  bf1[1] = bf0[1];
358
7.37M
  bf1[2] = bf0[2];
359
7.37M
  bf1[3] = bf0[3];
360
7.37M
  bf1[4] = bf0[4];
361
7.37M
  bf1[5] = bf0[5];
362
7.37M
  bf1[6] = bf0[6];
363
7.37M
  bf1[7] = bf0[7];
364
7.37M
  bf1[8] = bf0[8];
365
7.37M
  bf1[9] = bf0[9];
366
7.37M
  bf1[10] = bf0[10];
367
7.37M
  bf1[11] = bf0[11];
368
7.37M
  bf1[12] = bf0[12];
369
7.37M
  bf1[13] = bf0[13];
370
7.37M
  bf1[14] = bf0[14];
371
7.37M
  bf1[15] = bf0[15];
372
7.37M
  bf1[16] = half_btf(cospi[62], bf0[16], -cospi[2], bf0[31], cos_bit);
373
7.37M
  bf1[17] = half_btf(cospi[30], bf0[17], -cospi[34], bf0[30], cos_bit);
374
7.37M
  bf1[18] = half_btf(cospi[46], bf0[18], -cospi[18], bf0[29], cos_bit);
375
7.37M
  bf1[19] = half_btf(cospi[14], bf0[19], -cospi[50], bf0[28], cos_bit);
376
7.37M
  bf1[20] = half_btf(cospi[54], bf0[20], -cospi[10], bf0[27], cos_bit);
377
7.37M
  bf1[21] = half_btf(cospi[22], bf0[21], -cospi[42], bf0[26], cos_bit);
378
7.37M
  bf1[22] = half_btf(cospi[38], bf0[22], -cospi[26], bf0[25], cos_bit);
379
7.37M
  bf1[23] = half_btf(cospi[6], bf0[23], -cospi[58], bf0[24], cos_bit);
380
7.37M
  bf1[24] = half_btf(cospi[58], bf0[23], cospi[6], bf0[24], cos_bit);
381
7.37M
  bf1[25] = half_btf(cospi[26], bf0[22], cospi[38], bf0[25], cos_bit);
382
7.37M
  bf1[26] = half_btf(cospi[42], bf0[21], cospi[22], bf0[26], cos_bit);
383
7.37M
  bf1[27] = half_btf(cospi[10], bf0[20], cospi[54], bf0[27], cos_bit);
384
7.37M
  bf1[28] = half_btf(cospi[50], bf0[19], cospi[14], bf0[28], cos_bit);
385
7.37M
  bf1[29] = half_btf(cospi[18], bf0[18], cospi[46], bf0[29], cos_bit);
386
7.37M
  bf1[30] = half_btf(cospi[34], bf0[17], cospi[30], bf0[30], cos_bit);
387
7.37M
  bf1[31] = half_btf(cospi[2], bf0[16], cospi[62], bf0[31], cos_bit);
388
7.37M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
389
390
  // stage 3
391
7.37M
  stage++;
392
7.37M
  bf0 = step;
393
7.37M
  bf1 = output;
394
7.37M
  bf1[0] = bf0[0];
395
7.37M
  bf1[1] = bf0[1];
396
7.37M
  bf1[2] = bf0[2];
397
7.37M
  bf1[3] = bf0[3];
398
7.37M
  bf1[4] = bf0[4];
399
7.37M
  bf1[5] = bf0[5];
400
7.37M
  bf1[6] = bf0[6];
401
7.37M
  bf1[7] = bf0[7];
402
7.37M
  bf1[8] = half_btf(cospi[60], bf0[8], -cospi[4], bf0[15], cos_bit);
403
7.37M
  bf1[9] = half_btf(cospi[28], bf0[9], -cospi[36], bf0[14], cos_bit);
404
7.37M
  bf1[10] = half_btf(cospi[44], bf0[10], -cospi[20], bf0[13], cos_bit);
405
7.37M
  bf1[11] = half_btf(cospi[12], bf0[11], -cospi[52], bf0[12], cos_bit);
406
7.37M
  bf1[12] = half_btf(cospi[52], bf0[11], cospi[12], bf0[12], cos_bit);
407
7.37M
  bf1[13] = half_btf(cospi[20], bf0[10], cospi[44], bf0[13], cos_bit);
408
7.37M
  bf1[14] = half_btf(cospi[36], bf0[9], cospi[28], bf0[14], cos_bit);
409
7.37M
  bf1[15] = half_btf(cospi[4], bf0[8], cospi[60], bf0[15], cos_bit);
410
7.37M
  bf1[16] = clamp_value(bf0[16] + bf0[17], stage_range[stage]);
411
7.37M
  bf1[17] = clamp_value(bf0[16] - bf0[17], stage_range[stage]);
412
7.37M
  bf1[18] = clamp_value(-bf0[18] + bf0[19], stage_range[stage]);
413
7.37M
  bf1[19] = clamp_value(bf0[18] + bf0[19], stage_range[stage]);
414
7.37M
  bf1[20] = clamp_value(bf0[20] + bf0[21], stage_range[stage]);
415
7.37M
  bf1[21] = clamp_value(bf0[20] - bf0[21], stage_range[stage]);
416
7.37M
  bf1[22] = clamp_value(-bf0[22] + bf0[23], stage_range[stage]);
417
7.37M
  bf1[23] = clamp_value(bf0[22] + bf0[23], stage_range[stage]);
418
7.37M
  bf1[24] = clamp_value(bf0[24] + bf0[25], stage_range[stage]);
419
7.37M
  bf1[25] = clamp_value(bf0[24] - bf0[25], stage_range[stage]);
420
7.37M
  bf1[26] = clamp_value(-bf0[26] + bf0[27], stage_range[stage]);
421
7.37M
  bf1[27] = clamp_value(bf0[26] + bf0[27], stage_range[stage]);
422
7.37M
  bf1[28] = clamp_value(bf0[28] + bf0[29], stage_range[stage]);
423
7.37M
  bf1[29] = clamp_value(bf0[28] - bf0[29], stage_range[stage]);
424
7.37M
  bf1[30] = clamp_value(-bf0[30] + bf0[31], stage_range[stage]);
425
7.37M
  bf1[31] = clamp_value(bf0[30] + bf0[31], stage_range[stage]);
426
7.37M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
427
428
  // stage 4
429
7.37M
  stage++;
430
7.37M
  bf0 = output;
431
7.37M
  bf1 = step;
432
7.37M
  bf1[0] = bf0[0];
433
7.37M
  bf1[1] = bf0[1];
434
7.37M
  bf1[2] = bf0[2];
435
7.37M
  bf1[3] = bf0[3];
436
7.37M
  bf1[4] = half_btf(cospi[56], bf0[4], -cospi[8], bf0[7], cos_bit);
437
7.37M
  bf1[5] = half_btf(cospi[24], bf0[5], -cospi[40], bf0[6], cos_bit);
438
7.37M
  bf1[6] = half_btf(cospi[40], bf0[5], cospi[24], bf0[6], cos_bit);
439
7.37M
  bf1[7] = half_btf(cospi[8], bf0[4], cospi[56], bf0[7], cos_bit);
440
7.37M
  bf1[8] = clamp_value(bf0[8] + bf0[9], stage_range[stage]);
441
7.37M
  bf1[9] = clamp_value(bf0[8] - bf0[9], stage_range[stage]);
442
7.37M
  bf1[10] = clamp_value(-bf0[10] + bf0[11], stage_range[stage]);
443
7.37M
  bf1[11] = clamp_value(bf0[10] + bf0[11], stage_range[stage]);
444
7.37M
  bf1[12] = clamp_value(bf0[12] + bf0[13], stage_range[stage]);
445
7.37M
  bf1[13] = clamp_value(bf0[12] - bf0[13], stage_range[stage]);
446
7.37M
  bf1[14] = clamp_value(-bf0[14] + bf0[15], stage_range[stage]);
447
7.37M
  bf1[15] = clamp_value(bf0[14] + bf0[15], stage_range[stage]);
448
7.37M
  bf1[16] = bf0[16];
449
7.37M
  bf1[17] = half_btf(-cospi[8], bf0[17], cospi[56], bf0[30], cos_bit);
450
7.37M
  bf1[18] = half_btf(-cospi[56], bf0[18], -cospi[8], bf0[29], cos_bit);
451
7.37M
  bf1[19] = bf0[19];
452
7.37M
  bf1[20] = bf0[20];
453
7.37M
  bf1[21] = half_btf(-cospi[40], bf0[21], cospi[24], bf0[26], cos_bit);
454
7.37M
  bf1[22] = half_btf(-cospi[24], bf0[22], -cospi[40], bf0[25], cos_bit);
455
7.37M
  bf1[23] = bf0[23];
456
7.37M
  bf1[24] = bf0[24];
457
7.37M
  bf1[25] = half_btf(-cospi[40], bf0[22], cospi[24], bf0[25], cos_bit);
458
7.37M
  bf1[26] = half_btf(cospi[24], bf0[21], cospi[40], bf0[26], cos_bit);
459
7.37M
  bf1[27] = bf0[27];
460
7.37M
  bf1[28] = bf0[28];
461
7.37M
  bf1[29] = half_btf(-cospi[8], bf0[18], cospi[56], bf0[29], cos_bit);
462
7.37M
  bf1[30] = half_btf(cospi[56], bf0[17], cospi[8], bf0[30], cos_bit);
463
7.37M
  bf1[31] = bf0[31];
464
7.37M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
465
466
  // stage 5
467
7.37M
  stage++;
468
7.37M
  bf0 = step;
469
7.37M
  bf1 = output;
470
7.37M
  bf1[0] = half_btf(cospi[32], bf0[0], cospi[32], bf0[1], cos_bit);
471
7.37M
  bf1[1] = half_btf(cospi[32], bf0[0], -cospi[32], bf0[1], cos_bit);
472
7.37M
  bf1[2] = half_btf(cospi[48], bf0[2], -cospi[16], bf0[3], cos_bit);
473
7.37M
  bf1[3] = half_btf(cospi[16], bf0[2], cospi[48], bf0[3], cos_bit);
474
7.37M
  bf1[4] = clamp_value(bf0[4] + bf0[5], stage_range[stage]);
475
7.37M
  bf1[5] = clamp_value(bf0[4] - bf0[5], stage_range[stage]);
476
7.37M
  bf1[6] = clamp_value(-bf0[6] + bf0[7], stage_range[stage]);
477
7.37M
  bf1[7] = clamp_value(bf0[6] + bf0[7], stage_range[stage]);
478
7.37M
  bf1[8] = bf0[8];
479
7.37M
  bf1[9] = half_btf(-cospi[16], bf0[9], cospi[48], bf0[14], cos_bit);
480
7.37M
  bf1[10] = half_btf(-cospi[48], bf0[10], -cospi[16], bf0[13], cos_bit);
481
7.37M
  bf1[11] = bf0[11];
482
7.37M
  bf1[12] = bf0[12];
483
7.37M
  bf1[13] = half_btf(-cospi[16], bf0[10], cospi[48], bf0[13], cos_bit);
484
7.37M
  bf1[14] = half_btf(cospi[48], bf0[9], cospi[16], bf0[14], cos_bit);
485
7.37M
  bf1[15] = bf0[15];
486
7.37M
  bf1[16] = clamp_value(bf0[16] + bf0[19], stage_range[stage]);
487
7.37M
  bf1[17] = clamp_value(bf0[17] + bf0[18], stage_range[stage]);
488
7.37M
  bf1[18] = clamp_value(bf0[17] - bf0[18], stage_range[stage]);
489
7.37M
  bf1[19] = clamp_value(bf0[16] - bf0[19], stage_range[stage]);
490
7.37M
  bf1[20] = clamp_value(-bf0[20] + bf0[23], stage_range[stage]);
491
7.37M
  bf1[21] = clamp_value(-bf0[21] + bf0[22], stage_range[stage]);
492
7.37M
  bf1[22] = clamp_value(bf0[21] + bf0[22], stage_range[stage]);
493
7.37M
  bf1[23] = clamp_value(bf0[20] + bf0[23], stage_range[stage]);
494
7.37M
  bf1[24] = clamp_value(bf0[24] + bf0[27], stage_range[stage]);
495
7.37M
  bf1[25] = clamp_value(bf0[25] + bf0[26], stage_range[stage]);
496
7.37M
  bf1[26] = clamp_value(bf0[25] - bf0[26], stage_range[stage]);
497
7.37M
  bf1[27] = clamp_value(bf0[24] - bf0[27], stage_range[stage]);
498
7.37M
  bf1[28] = clamp_value(-bf0[28] + bf0[31], stage_range[stage]);
499
7.37M
  bf1[29] = clamp_value(-bf0[29] + bf0[30], stage_range[stage]);
500
7.37M
  bf1[30] = clamp_value(bf0[29] + bf0[30], stage_range[stage]);
501
7.37M
  bf1[31] = clamp_value(bf0[28] + bf0[31], stage_range[stage]);
502
7.37M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
503
504
  // stage 6
505
7.37M
  stage++;
506
7.37M
  bf0 = output;
507
7.37M
  bf1 = step;
508
7.37M
  bf1[0] = clamp_value(bf0[0] + bf0[3], stage_range[stage]);
509
7.37M
  bf1[1] = clamp_value(bf0[1] + bf0[2], stage_range[stage]);
510
7.37M
  bf1[2] = clamp_value(bf0[1] - bf0[2], stage_range[stage]);
511
7.37M
  bf1[3] = clamp_value(bf0[0] - bf0[3], stage_range[stage]);
512
7.37M
  bf1[4] = bf0[4];
513
7.37M
  bf1[5] = half_btf(-cospi[32], bf0[5], cospi[32], bf0[6], cos_bit);
514
7.37M
  bf1[6] = half_btf(cospi[32], bf0[5], cospi[32], bf0[6], cos_bit);
515
7.37M
  bf1[7] = bf0[7];
516
7.37M
  bf1[8] = clamp_value(bf0[8] + bf0[11], stage_range[stage]);
517
7.37M
  bf1[9] = clamp_value(bf0[9] + bf0[10], stage_range[stage]);
518
7.37M
  bf1[10] = clamp_value(bf0[9] - bf0[10], stage_range[stage]);
519
7.37M
  bf1[11] = clamp_value(bf0[8] - bf0[11], stage_range[stage]);
520
7.37M
  bf1[12] = clamp_value(-bf0[12] + bf0[15], stage_range[stage]);
521
7.37M
  bf1[13] = clamp_value(-bf0[13] + bf0[14], stage_range[stage]);
522
7.37M
  bf1[14] = clamp_value(bf0[13] + bf0[14], stage_range[stage]);
523
7.37M
  bf1[15] = clamp_value(bf0[12] + bf0[15], stage_range[stage]);
524
7.37M
  bf1[16] = bf0[16];
525
7.37M
  bf1[17] = bf0[17];
526
7.37M
  bf1[18] = half_btf(-cospi[16], bf0[18], cospi[48], bf0[29], cos_bit);
527
7.37M
  bf1[19] = half_btf(-cospi[16], bf0[19], cospi[48], bf0[28], cos_bit);
528
7.37M
  bf1[20] = half_btf(-cospi[48], bf0[20], -cospi[16], bf0[27], cos_bit);
529
7.37M
  bf1[21] = half_btf(-cospi[48], bf0[21], -cospi[16], bf0[26], cos_bit);
530
7.37M
  bf1[22] = bf0[22];
531
7.37M
  bf1[23] = bf0[23];
532
7.37M
  bf1[24] = bf0[24];
533
7.37M
  bf1[25] = bf0[25];
534
7.37M
  bf1[26] = half_btf(-cospi[16], bf0[21], cospi[48], bf0[26], cos_bit);
535
7.37M
  bf1[27] = half_btf(-cospi[16], bf0[20], cospi[48], bf0[27], cos_bit);
536
7.37M
  bf1[28] = half_btf(cospi[48], bf0[19], cospi[16], bf0[28], cos_bit);
537
7.37M
  bf1[29] = half_btf(cospi[48], bf0[18], cospi[16], bf0[29], cos_bit);
538
7.37M
  bf1[30] = bf0[30];
539
7.37M
  bf1[31] = bf0[31];
540
7.37M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
541
542
  // stage 7
543
7.37M
  stage++;
544
7.37M
  bf0 = step;
545
7.37M
  bf1 = output;
546
7.37M
  bf1[0] = clamp_value(bf0[0] + bf0[7], stage_range[stage]);
547
7.37M
  bf1[1] = clamp_value(bf0[1] + bf0[6], stage_range[stage]);
548
7.37M
  bf1[2] = clamp_value(bf0[2] + bf0[5], stage_range[stage]);
549
7.37M
  bf1[3] = clamp_value(bf0[3] + bf0[4], stage_range[stage]);
550
7.37M
  bf1[4] = clamp_value(bf0[3] - bf0[4], stage_range[stage]);
551
7.37M
  bf1[5] = clamp_value(bf0[2] - bf0[5], stage_range[stage]);
552
7.37M
  bf1[6] = clamp_value(bf0[1] - bf0[6], stage_range[stage]);
553
7.37M
  bf1[7] = clamp_value(bf0[0] - bf0[7], stage_range[stage]);
554
7.37M
  bf1[8] = bf0[8];
555
7.37M
  bf1[9] = bf0[9];
556
7.37M
  bf1[10] = half_btf(-cospi[32], bf0[10], cospi[32], bf0[13], cos_bit);
557
7.37M
  bf1[11] = half_btf(-cospi[32], bf0[11], cospi[32], bf0[12], cos_bit);
558
7.37M
  bf1[12] = half_btf(cospi[32], bf0[11], cospi[32], bf0[12], cos_bit);
559
7.37M
  bf1[13] = half_btf(cospi[32], bf0[10], cospi[32], bf0[13], cos_bit);
560
7.37M
  bf1[14] = bf0[14];
561
7.37M
  bf1[15] = bf0[15];
562
7.37M
  bf1[16] = clamp_value(bf0[16] + bf0[23], stage_range[stage]);
563
7.37M
  bf1[17] = clamp_value(bf0[17] + bf0[22], stage_range[stage]);
564
7.37M
  bf1[18] = clamp_value(bf0[18] + bf0[21], stage_range[stage]);
565
7.37M
  bf1[19] = clamp_value(bf0[19] + bf0[20], stage_range[stage]);
566
7.37M
  bf1[20] = clamp_value(bf0[19] - bf0[20], stage_range[stage]);
567
7.37M
  bf1[21] = clamp_value(bf0[18] - bf0[21], stage_range[stage]);
568
7.37M
  bf1[22] = clamp_value(bf0[17] - bf0[22], stage_range[stage]);
569
7.37M
  bf1[23] = clamp_value(bf0[16] - bf0[23], stage_range[stage]);
570
7.37M
  bf1[24] = clamp_value(-bf0[24] + bf0[31], stage_range[stage]);
571
7.37M
  bf1[25] = clamp_value(-bf0[25] + bf0[30], stage_range[stage]);
572
7.37M
  bf1[26] = clamp_value(-bf0[26] + bf0[29], stage_range[stage]);
573
7.37M
  bf1[27] = clamp_value(-bf0[27] + bf0[28], stage_range[stage]);
574
7.37M
  bf1[28] = clamp_value(bf0[27] + bf0[28], stage_range[stage]);
575
7.37M
  bf1[29] = clamp_value(bf0[26] + bf0[29], stage_range[stage]);
576
7.37M
  bf1[30] = clamp_value(bf0[25] + bf0[30], stage_range[stage]);
577
7.37M
  bf1[31] = clamp_value(bf0[24] + bf0[31], stage_range[stage]);
578
7.37M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
579
580
  // stage 8
581
7.37M
  stage++;
582
7.37M
  bf0 = output;
583
7.37M
  bf1 = step;
584
7.37M
  bf1[0] = clamp_value(bf0[0] + bf0[15], stage_range[stage]);
585
7.37M
  bf1[1] = clamp_value(bf0[1] + bf0[14], stage_range[stage]);
586
7.37M
  bf1[2] = clamp_value(bf0[2] + bf0[13], stage_range[stage]);
587
7.37M
  bf1[3] = clamp_value(bf0[3] + bf0[12], stage_range[stage]);
588
7.37M
  bf1[4] = clamp_value(bf0[4] + bf0[11], stage_range[stage]);
589
7.37M
  bf1[5] = clamp_value(bf0[5] + bf0[10], stage_range[stage]);
590
7.37M
  bf1[6] = clamp_value(bf0[6] + bf0[9], stage_range[stage]);
591
7.37M
  bf1[7] = clamp_value(bf0[7] + bf0[8], stage_range[stage]);
592
7.37M
  bf1[8] = clamp_value(bf0[7] - bf0[8], stage_range[stage]);
593
7.37M
  bf1[9] = clamp_value(bf0[6] - bf0[9], stage_range[stage]);
594
7.37M
  bf1[10] = clamp_value(bf0[5] - bf0[10], stage_range[stage]);
595
7.37M
  bf1[11] = clamp_value(bf0[4] - bf0[11], stage_range[stage]);
596
7.37M
  bf1[12] = clamp_value(bf0[3] - bf0[12], stage_range[stage]);
597
7.37M
  bf1[13] = clamp_value(bf0[2] - bf0[13], stage_range[stage]);
598
7.37M
  bf1[14] = clamp_value(bf0[1] - bf0[14], stage_range[stage]);
599
7.37M
  bf1[15] = clamp_value(bf0[0] - bf0[15], stage_range[stage]);
600
7.37M
  bf1[16] = bf0[16];
601
7.37M
  bf1[17] = bf0[17];
602
7.37M
  bf1[18] = bf0[18];
603
7.37M
  bf1[19] = bf0[19];
604
7.37M
  bf1[20] = half_btf(-cospi[32], bf0[20], cospi[32], bf0[27], cos_bit);
605
7.37M
  bf1[21] = half_btf(-cospi[32], bf0[21], cospi[32], bf0[26], cos_bit);
606
7.37M
  bf1[22] = half_btf(-cospi[32], bf0[22], cospi[32], bf0[25], cos_bit);
607
7.37M
  bf1[23] = half_btf(-cospi[32], bf0[23], cospi[32], bf0[24], cos_bit);
608
7.37M
  bf1[24] = half_btf(cospi[32], bf0[23], cospi[32], bf0[24], cos_bit);
609
7.37M
  bf1[25] = half_btf(cospi[32], bf0[22], cospi[32], bf0[25], cos_bit);
610
7.37M
  bf1[26] = half_btf(cospi[32], bf0[21], cospi[32], bf0[26], cos_bit);
611
7.37M
  bf1[27] = half_btf(cospi[32], bf0[20], cospi[32], bf0[27], cos_bit);
612
7.37M
  bf1[28] = bf0[28];
613
7.37M
  bf1[29] = bf0[29];
614
7.37M
  bf1[30] = bf0[30];
615
7.37M
  bf1[31] = bf0[31];
616
7.37M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
617
618
  // stage 9
619
7.37M
  stage++;
620
7.37M
  bf0 = step;
621
7.37M
  bf1 = output;
622
7.37M
  bf1[0] = clamp_value(bf0[0] + bf0[31], stage_range[stage]);
623
7.37M
  bf1[1] = clamp_value(bf0[1] + bf0[30], stage_range[stage]);
624
7.37M
  bf1[2] = clamp_value(bf0[2] + bf0[29], stage_range[stage]);
625
7.37M
  bf1[3] = clamp_value(bf0[3] + bf0[28], stage_range[stage]);
626
7.37M
  bf1[4] = clamp_value(bf0[4] + bf0[27], stage_range[stage]);
627
7.37M
  bf1[5] = clamp_value(bf0[5] + bf0[26], stage_range[stage]);
628
7.37M
  bf1[6] = clamp_value(bf0[6] + bf0[25], stage_range[stage]);
629
7.37M
  bf1[7] = clamp_value(bf0[7] + bf0[24], stage_range[stage]);
630
7.37M
  bf1[8] = clamp_value(bf0[8] + bf0[23], stage_range[stage]);
631
7.37M
  bf1[9] = clamp_value(bf0[9] + bf0[22], stage_range[stage]);
632
7.37M
  bf1[10] = clamp_value(bf0[10] + bf0[21], stage_range[stage]);
633
7.37M
  bf1[11] = clamp_value(bf0[11] + bf0[20], stage_range[stage]);
634
7.37M
  bf1[12] = clamp_value(bf0[12] + bf0[19], stage_range[stage]);
635
7.37M
  bf1[13] = clamp_value(bf0[13] + bf0[18], stage_range[stage]);
636
7.37M
  bf1[14] = clamp_value(bf0[14] + bf0[17], stage_range[stage]);
637
7.37M
  bf1[15] = clamp_value(bf0[15] + bf0[16], stage_range[stage]);
638
7.37M
  bf1[16] = clamp_value(bf0[15] - bf0[16], stage_range[stage]);
639
7.37M
  bf1[17] = clamp_value(bf0[14] - bf0[17], stage_range[stage]);
640
7.37M
  bf1[18] = clamp_value(bf0[13] - bf0[18], stage_range[stage]);
641
7.37M
  bf1[19] = clamp_value(bf0[12] - bf0[19], stage_range[stage]);
642
7.37M
  bf1[20] = clamp_value(bf0[11] - bf0[20], stage_range[stage]);
643
7.37M
  bf1[21] = clamp_value(bf0[10] - bf0[21], stage_range[stage]);
644
7.37M
  bf1[22] = clamp_value(bf0[9] - bf0[22], stage_range[stage]);
645
7.37M
  bf1[23] = clamp_value(bf0[8] - bf0[23], stage_range[stage]);
646
7.37M
  bf1[24] = clamp_value(bf0[7] - bf0[24], stage_range[stage]);
647
7.37M
  bf1[25] = clamp_value(bf0[6] - bf0[25], stage_range[stage]);
648
7.37M
  bf1[26] = clamp_value(bf0[5] - bf0[26], stage_range[stage]);
649
7.37M
  bf1[27] = clamp_value(bf0[4] - bf0[27], stage_range[stage]);
650
7.37M
  bf1[28] = clamp_value(bf0[3] - bf0[28], stage_range[stage]);
651
7.37M
  bf1[29] = clamp_value(bf0[2] - bf0[29], stage_range[stage]);
652
7.37M
  bf1[30] = clamp_value(bf0[1] - bf0[30], stage_range[stage]);
653
7.37M
  bf1[31] = clamp_value(bf0[0] - bf0[31], stage_range[stage]);
654
7.37M
}
655
656
void av1_iadst4(const int32_t *input, int32_t *output, int8_t cos_bit,
657
901k
                const int8_t *stage_range) {
658
901k
  int bit = cos_bit;
659
901k
  const int32_t *sinpi = sinpi_arr(bit);
660
901k
  int32_t s0, s1, s2, s3, s4, s5, s6, s7;
661
662
901k
  int32_t x0 = input[0];
663
901k
  int32_t x1 = input[1];
664
901k
  int32_t x2 = input[2];
665
901k
  int32_t x3 = input[3];
666
667
901k
  if (!(x0 | x1 | x2 | x3)) {
668
284k
    output[0] = output[1] = output[2] = output[3] = 0;
669
284k
    return;
670
284k
  }
671
672
617k
  assert(sinpi[1] + sinpi[2] == sinpi[4]);
673
674
  // stage 1
675
617k
  s0 = range_check_value(sinpi[1] * x0, stage_range[1] + bit);
676
617k
  s1 = range_check_value(sinpi[2] * x0, stage_range[1] + bit);
677
617k
  s2 = range_check_value(sinpi[3] * x1, stage_range[1] + bit);
678
617k
  s3 = range_check_value(sinpi[4] * x2, stage_range[1] + bit);
679
617k
  s4 = range_check_value(sinpi[1] * x2, stage_range[1] + bit);
680
617k
  s5 = range_check_value(sinpi[2] * x3, stage_range[1] + bit);
681
617k
  s6 = range_check_value(sinpi[4] * x3, stage_range[1] + bit);
682
683
  // stage 2
684
  // NOTICE: (x0 - x2) here may use one extra bit compared to the
685
  // opt_range_row/col specified in av1_gen_inv_stage_range()
686
617k
  s7 = range_check_value((x0 - x2) + x3, stage_range[2]);
687
688
  // stage 3
689
617k
  s0 = range_check_value(s0 + s3, stage_range[3] + bit);
690
617k
  s1 = range_check_value(s1 - s4, stage_range[3] + bit);
691
617k
  s3 = range_check_value(s2, stage_range[3] + bit);
692
617k
  s2 = range_check_value(sinpi[3] * s7, stage_range[3] + bit);
693
694
  // stage 4
695
617k
  s0 = range_check_value(s0 + s5, stage_range[4] + bit);
696
617k
  s1 = range_check_value(s1 - s6, stage_range[4] + bit);
697
698
  // stage 5
699
617k
  x0 = range_check_value(s0 + s3, stage_range[5] + bit);
700
617k
  x1 = range_check_value(s1 + s3, stage_range[5] + bit);
701
617k
  x2 = range_check_value(s2, stage_range[5] + bit);
702
617k
  x3 = range_check_value(s0 + s1, stage_range[5] + bit);
703
704
  // stage 6
705
617k
  x3 = range_check_value(x3 - s3, stage_range[6] + bit);
706
707
617k
  output[0] = round_shift(x0, bit);
708
617k
  output[1] = round_shift(x1, bit);
709
617k
  output[2] = round_shift(x2, bit);
710
617k
  output[3] = round_shift(x3, bit);
711
617k
}
712
713
void av1_iadst8(const int32_t *input, int32_t *output, int8_t cos_bit,
714
1.24M
                const int8_t *stage_range) {
715
1.24M
  assert(output != input);
716
1.24M
  const int32_t size = 8;
717
1.24M
  const int32_t *cospi = cospi_arr(cos_bit);
718
719
1.24M
  int32_t stage = 0;
720
1.24M
  int32_t *bf0, *bf1;
721
1.24M
  int32_t step[8];
722
723
  // stage 0;
724
725
  // stage 1;
726
1.24M
  stage++;
727
1.24M
  bf1 = output;
728
1.24M
  bf1[0] = input[7];
729
1.24M
  bf1[1] = input[0];
730
1.24M
  bf1[2] = input[5];
731
1.24M
  bf1[3] = input[2];
732
1.24M
  bf1[4] = input[3];
733
1.24M
  bf1[5] = input[4];
734
1.24M
  bf1[6] = input[1];
735
1.24M
  bf1[7] = input[6];
736
1.24M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
737
738
  // stage 2
739
1.24M
  stage++;
740
1.24M
  bf0 = output;
741
1.24M
  bf1 = step;
742
1.24M
  bf1[0] = half_btf(cospi[4], bf0[0], cospi[60], bf0[1], cos_bit);
743
1.24M
  bf1[1] = half_btf(cospi[60], bf0[0], -cospi[4], bf0[1], cos_bit);
744
1.24M
  bf1[2] = half_btf(cospi[20], bf0[2], cospi[44], bf0[3], cos_bit);
745
1.24M
  bf1[3] = half_btf(cospi[44], bf0[2], -cospi[20], bf0[3], cos_bit);
746
1.24M
  bf1[4] = half_btf(cospi[36], bf0[4], cospi[28], bf0[5], cos_bit);
747
1.24M
  bf1[5] = half_btf(cospi[28], bf0[4], -cospi[36], bf0[5], cos_bit);
748
1.24M
  bf1[6] = half_btf(cospi[52], bf0[6], cospi[12], bf0[7], cos_bit);
749
1.24M
  bf1[7] = half_btf(cospi[12], bf0[6], -cospi[52], bf0[7], cos_bit);
750
1.24M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
751
752
  // stage 3
753
1.24M
  stage++;
754
1.24M
  bf0 = step;
755
1.24M
  bf1 = output;
756
1.24M
  bf1[0] = clamp_value(bf0[0] + bf0[4], stage_range[stage]);
757
1.24M
  bf1[1] = clamp_value(bf0[1] + bf0[5], stage_range[stage]);
758
1.24M
  bf1[2] = clamp_value(bf0[2] + bf0[6], stage_range[stage]);
759
1.24M
  bf1[3] = clamp_value(bf0[3] + bf0[7], stage_range[stage]);
760
1.24M
  bf1[4] = clamp_value(bf0[0] - bf0[4], stage_range[stage]);
761
1.24M
  bf1[5] = clamp_value(bf0[1] - bf0[5], stage_range[stage]);
762
1.24M
  bf1[6] = clamp_value(bf0[2] - bf0[6], stage_range[stage]);
763
1.24M
  bf1[7] = clamp_value(bf0[3] - bf0[7], stage_range[stage]);
764
1.24M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
765
766
  // stage 4
767
1.24M
  stage++;
768
1.24M
  bf0 = output;
769
1.24M
  bf1 = step;
770
1.24M
  bf1[0] = bf0[0];
771
1.24M
  bf1[1] = bf0[1];
772
1.24M
  bf1[2] = bf0[2];
773
1.24M
  bf1[3] = bf0[3];
774
1.24M
  bf1[4] = half_btf(cospi[16], bf0[4], cospi[48], bf0[5], cos_bit);
775
1.24M
  bf1[5] = half_btf(cospi[48], bf0[4], -cospi[16], bf0[5], cos_bit);
776
1.24M
  bf1[6] = half_btf(-cospi[48], bf0[6], cospi[16], bf0[7], cos_bit);
777
1.24M
  bf1[7] = half_btf(cospi[16], bf0[6], cospi[48], bf0[7], cos_bit);
778
1.24M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
779
780
  // stage 5
781
1.24M
  stage++;
782
1.24M
  bf0 = step;
783
1.24M
  bf1 = output;
784
1.24M
  bf1[0] = clamp_value(bf0[0] + bf0[2], stage_range[stage]);
785
1.24M
  bf1[1] = clamp_value(bf0[1] + bf0[3], stage_range[stage]);
786
1.24M
  bf1[2] = clamp_value(bf0[0] - bf0[2], stage_range[stage]);
787
1.24M
  bf1[3] = clamp_value(bf0[1] - bf0[3], stage_range[stage]);
788
1.24M
  bf1[4] = clamp_value(bf0[4] + bf0[6], stage_range[stage]);
789
1.24M
  bf1[5] = clamp_value(bf0[5] + bf0[7], stage_range[stage]);
790
1.24M
  bf1[6] = clamp_value(bf0[4] - bf0[6], stage_range[stage]);
791
1.24M
  bf1[7] = clamp_value(bf0[5] - bf0[7], stage_range[stage]);
792
1.24M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
793
794
  // stage 6
795
1.24M
  stage++;
796
1.24M
  bf0 = output;
797
1.24M
  bf1 = step;
798
1.24M
  bf1[0] = bf0[0];
799
1.24M
  bf1[1] = bf0[1];
800
1.24M
  bf1[2] = half_btf(cospi[32], bf0[2], cospi[32], bf0[3], cos_bit);
801
1.24M
  bf1[3] = half_btf(cospi[32], bf0[2], -cospi[32], bf0[3], cos_bit);
802
1.24M
  bf1[4] = bf0[4];
803
1.24M
  bf1[5] = bf0[5];
804
1.24M
  bf1[6] = half_btf(cospi[32], bf0[6], cospi[32], bf0[7], cos_bit);
805
1.24M
  bf1[7] = half_btf(cospi[32], bf0[6], -cospi[32], bf0[7], cos_bit);
806
1.24M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
807
808
  // stage 7
809
1.24M
  bf0 = step;
810
1.24M
  bf1 = output;
811
1.24M
  bf1[0] = bf0[0];
812
1.24M
  bf1[1] = -bf0[4];
813
1.24M
  bf1[2] = bf0[6];
814
1.24M
  bf1[3] = -bf0[2];
815
1.24M
  bf1[4] = bf0[3];
816
1.24M
  bf1[5] = -bf0[7];
817
1.24M
  bf1[6] = bf0[5];
818
1.24M
  bf1[7] = -bf0[1];
819
1.24M
}
820
821
void av1_iadst16(const int32_t *input, int32_t *output, int8_t cos_bit,
822
1.99M
                 const int8_t *stage_range) {
823
1.99M
  assert(output != input);
824
1.99M
  const int32_t size = 16;
825
1.99M
  const int32_t *cospi = cospi_arr(cos_bit);
826
827
1.99M
  int32_t stage = 0;
828
1.99M
  int32_t *bf0, *bf1;
829
1.99M
  int32_t step[16];
830
831
  // stage 0;
832
833
  // stage 1;
834
1.99M
  stage++;
835
1.99M
  bf1 = output;
836
1.99M
  bf1[0] = input[15];
837
1.99M
  bf1[1] = input[0];
838
1.99M
  bf1[2] = input[13];
839
1.99M
  bf1[3] = input[2];
840
1.99M
  bf1[4] = input[11];
841
1.99M
  bf1[5] = input[4];
842
1.99M
  bf1[6] = input[9];
843
1.99M
  bf1[7] = input[6];
844
1.99M
  bf1[8] = input[7];
845
1.99M
  bf1[9] = input[8];
846
1.99M
  bf1[10] = input[5];
847
1.99M
  bf1[11] = input[10];
848
1.99M
  bf1[12] = input[3];
849
1.99M
  bf1[13] = input[12];
850
1.99M
  bf1[14] = input[1];
851
1.99M
  bf1[15] = input[14];
852
1.99M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
853
854
  // stage 2
855
1.99M
  stage++;
856
1.99M
  bf0 = output;
857
1.99M
  bf1 = step;
858
1.99M
  bf1[0] = half_btf(cospi[2], bf0[0], cospi[62], bf0[1], cos_bit);
859
1.99M
  bf1[1] = half_btf(cospi[62], bf0[0], -cospi[2], bf0[1], cos_bit);
860
1.99M
  bf1[2] = half_btf(cospi[10], bf0[2], cospi[54], bf0[3], cos_bit);
861
1.99M
  bf1[3] = half_btf(cospi[54], bf0[2], -cospi[10], bf0[3], cos_bit);
862
1.99M
  bf1[4] = half_btf(cospi[18], bf0[4], cospi[46], bf0[5], cos_bit);
863
1.99M
  bf1[5] = half_btf(cospi[46], bf0[4], -cospi[18], bf0[5], cos_bit);
864
1.99M
  bf1[6] = half_btf(cospi[26], bf0[6], cospi[38], bf0[7], cos_bit);
865
1.99M
  bf1[7] = half_btf(cospi[38], bf0[6], -cospi[26], bf0[7], cos_bit);
866
1.99M
  bf1[8] = half_btf(cospi[34], bf0[8], cospi[30], bf0[9], cos_bit);
867
1.99M
  bf1[9] = half_btf(cospi[30], bf0[8], -cospi[34], bf0[9], cos_bit);
868
1.99M
  bf1[10] = half_btf(cospi[42], bf0[10], cospi[22], bf0[11], cos_bit);
869
1.99M
  bf1[11] = half_btf(cospi[22], bf0[10], -cospi[42], bf0[11], cos_bit);
870
1.99M
  bf1[12] = half_btf(cospi[50], bf0[12], cospi[14], bf0[13], cos_bit);
871
1.99M
  bf1[13] = half_btf(cospi[14], bf0[12], -cospi[50], bf0[13], cos_bit);
872
1.99M
  bf1[14] = half_btf(cospi[58], bf0[14], cospi[6], bf0[15], cos_bit);
873
1.99M
  bf1[15] = half_btf(cospi[6], bf0[14], -cospi[58], bf0[15], cos_bit);
874
1.99M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
875
876
  // stage 3
877
1.99M
  stage++;
878
1.99M
  bf0 = step;
879
1.99M
  bf1 = output;
880
1.99M
  bf1[0] = clamp_value(bf0[0] + bf0[8], stage_range[stage]);
881
1.99M
  bf1[1] = clamp_value(bf0[1] + bf0[9], stage_range[stage]);
882
1.99M
  bf1[2] = clamp_value(bf0[2] + bf0[10], stage_range[stage]);
883
1.99M
  bf1[3] = clamp_value(bf0[3] + bf0[11], stage_range[stage]);
884
1.99M
  bf1[4] = clamp_value(bf0[4] + bf0[12], stage_range[stage]);
885
1.99M
  bf1[5] = clamp_value(bf0[5] + bf0[13], stage_range[stage]);
886
1.99M
  bf1[6] = clamp_value(bf0[6] + bf0[14], stage_range[stage]);
887
1.99M
  bf1[7] = clamp_value(bf0[7] + bf0[15], stage_range[stage]);
888
1.99M
  bf1[8] = clamp_value(bf0[0] - bf0[8], stage_range[stage]);
889
1.99M
  bf1[9] = clamp_value(bf0[1] - bf0[9], stage_range[stage]);
890
1.99M
  bf1[10] = clamp_value(bf0[2] - bf0[10], stage_range[stage]);
891
1.99M
  bf1[11] = clamp_value(bf0[3] - bf0[11], stage_range[stage]);
892
1.99M
  bf1[12] = clamp_value(bf0[4] - bf0[12], stage_range[stage]);
893
1.99M
  bf1[13] = clamp_value(bf0[5] - bf0[13], stage_range[stage]);
894
1.99M
  bf1[14] = clamp_value(bf0[6] - bf0[14], stage_range[stage]);
895
1.99M
  bf1[15] = clamp_value(bf0[7] - bf0[15], stage_range[stage]);
896
1.99M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
897
898
  // stage 4
899
1.99M
  stage++;
900
1.99M
  bf0 = output;
901
1.99M
  bf1 = step;
902
1.99M
  bf1[0] = bf0[0];
903
1.99M
  bf1[1] = bf0[1];
904
1.99M
  bf1[2] = bf0[2];
905
1.99M
  bf1[3] = bf0[3];
906
1.99M
  bf1[4] = bf0[4];
907
1.99M
  bf1[5] = bf0[5];
908
1.99M
  bf1[6] = bf0[6];
909
1.99M
  bf1[7] = bf0[7];
910
1.99M
  bf1[8] = half_btf(cospi[8], bf0[8], cospi[56], bf0[9], cos_bit);
911
1.99M
  bf1[9] = half_btf(cospi[56], bf0[8], -cospi[8], bf0[9], cos_bit);
912
1.99M
  bf1[10] = half_btf(cospi[40], bf0[10], cospi[24], bf0[11], cos_bit);
913
1.99M
  bf1[11] = half_btf(cospi[24], bf0[10], -cospi[40], bf0[11], cos_bit);
914
1.99M
  bf1[12] = half_btf(-cospi[56], bf0[12], cospi[8], bf0[13], cos_bit);
915
1.99M
  bf1[13] = half_btf(cospi[8], bf0[12], cospi[56], bf0[13], cos_bit);
916
1.99M
  bf1[14] = half_btf(-cospi[24], bf0[14], cospi[40], bf0[15], cos_bit);
917
1.99M
  bf1[15] = half_btf(cospi[40], bf0[14], cospi[24], bf0[15], cos_bit);
918
1.99M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
919
920
  // stage 5
921
1.99M
  stage++;
922
1.99M
  bf0 = step;
923
1.99M
  bf1 = output;
924
1.99M
  bf1[0] = clamp_value(bf0[0] + bf0[4], stage_range[stage]);
925
1.99M
  bf1[1] = clamp_value(bf0[1] + bf0[5], stage_range[stage]);
926
1.99M
  bf1[2] = clamp_value(bf0[2] + bf0[6], stage_range[stage]);
927
1.99M
  bf1[3] = clamp_value(bf0[3] + bf0[7], stage_range[stage]);
928
1.99M
  bf1[4] = clamp_value(bf0[0] - bf0[4], stage_range[stage]);
929
1.99M
  bf1[5] = clamp_value(bf0[1] - bf0[5], stage_range[stage]);
930
1.99M
  bf1[6] = clamp_value(bf0[2] - bf0[6], stage_range[stage]);
931
1.99M
  bf1[7] = clamp_value(bf0[3] - bf0[7], stage_range[stage]);
932
1.99M
  bf1[8] = clamp_value(bf0[8] + bf0[12], stage_range[stage]);
933
1.99M
  bf1[9] = clamp_value(bf0[9] + bf0[13], stage_range[stage]);
934
1.99M
  bf1[10] = clamp_value(bf0[10] + bf0[14], stage_range[stage]);
935
1.99M
  bf1[11] = clamp_value(bf0[11] + bf0[15], stage_range[stage]);
936
1.99M
  bf1[12] = clamp_value(bf0[8] - bf0[12], stage_range[stage]);
937
1.99M
  bf1[13] = clamp_value(bf0[9] - bf0[13], stage_range[stage]);
938
1.99M
  bf1[14] = clamp_value(bf0[10] - bf0[14], stage_range[stage]);
939
1.99M
  bf1[15] = clamp_value(bf0[11] - bf0[15], stage_range[stage]);
940
1.99M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
941
942
  // stage 6
943
1.99M
  stage++;
944
1.99M
  bf0 = output;
945
1.99M
  bf1 = step;
946
1.99M
  bf1[0] = bf0[0];
947
1.99M
  bf1[1] = bf0[1];
948
1.99M
  bf1[2] = bf0[2];
949
1.99M
  bf1[3] = bf0[3];
950
1.99M
  bf1[4] = half_btf(cospi[16], bf0[4], cospi[48], bf0[5], cos_bit);
951
1.99M
  bf1[5] = half_btf(cospi[48], bf0[4], -cospi[16], bf0[5], cos_bit);
952
1.99M
  bf1[6] = half_btf(-cospi[48], bf0[6], cospi[16], bf0[7], cos_bit);
953
1.99M
  bf1[7] = half_btf(cospi[16], bf0[6], cospi[48], bf0[7], cos_bit);
954
1.99M
  bf1[8] = bf0[8];
955
1.99M
  bf1[9] = bf0[9];
956
1.99M
  bf1[10] = bf0[10];
957
1.99M
  bf1[11] = bf0[11];
958
1.99M
  bf1[12] = half_btf(cospi[16], bf0[12], cospi[48], bf0[13], cos_bit);
959
1.99M
  bf1[13] = half_btf(cospi[48], bf0[12], -cospi[16], bf0[13], cos_bit);
960
1.99M
  bf1[14] = half_btf(-cospi[48], bf0[14], cospi[16], bf0[15], cos_bit);
961
1.99M
  bf1[15] = half_btf(cospi[16], bf0[14], cospi[48], bf0[15], cos_bit);
962
1.99M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
963
964
  // stage 7
965
1.99M
  stage++;
966
1.99M
  bf0 = step;
967
1.99M
  bf1 = output;
968
1.99M
  bf1[0] = clamp_value(bf0[0] + bf0[2], stage_range[stage]);
969
1.99M
  bf1[1] = clamp_value(bf0[1] + bf0[3], stage_range[stage]);
970
1.99M
  bf1[2] = clamp_value(bf0[0] - bf0[2], stage_range[stage]);
971
1.99M
  bf1[3] = clamp_value(bf0[1] - bf0[3], stage_range[stage]);
972
1.99M
  bf1[4] = clamp_value(bf0[4] + bf0[6], stage_range[stage]);
973
1.99M
  bf1[5] = clamp_value(bf0[5] + bf0[7], stage_range[stage]);
974
1.99M
  bf1[6] = clamp_value(bf0[4] - bf0[6], stage_range[stage]);
975
1.99M
  bf1[7] = clamp_value(bf0[5] - bf0[7], stage_range[stage]);
976
1.99M
  bf1[8] = clamp_value(bf0[8] + bf0[10], stage_range[stage]);
977
1.99M
  bf1[9] = clamp_value(bf0[9] + bf0[11], stage_range[stage]);
978
1.99M
  bf1[10] = clamp_value(bf0[8] - bf0[10], stage_range[stage]);
979
1.99M
  bf1[11] = clamp_value(bf0[9] - bf0[11], stage_range[stage]);
980
1.99M
  bf1[12] = clamp_value(bf0[12] + bf0[14], stage_range[stage]);
981
1.99M
  bf1[13] = clamp_value(bf0[13] + bf0[15], stage_range[stage]);
982
1.99M
  bf1[14] = clamp_value(bf0[12] - bf0[14], stage_range[stage]);
983
1.99M
  bf1[15] = clamp_value(bf0[13] - bf0[15], stage_range[stage]);
984
1.99M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
985
986
  // stage 8
987
1.99M
  stage++;
988
1.99M
  bf0 = output;
989
1.99M
  bf1 = step;
990
1.99M
  bf1[0] = bf0[0];
991
1.99M
  bf1[1] = bf0[1];
992
1.99M
  bf1[2] = half_btf(cospi[32], bf0[2], cospi[32], bf0[3], cos_bit);
993
1.99M
  bf1[3] = half_btf(cospi[32], bf0[2], -cospi[32], bf0[3], cos_bit);
994
1.99M
  bf1[4] = bf0[4];
995
1.99M
  bf1[5] = bf0[5];
996
1.99M
  bf1[6] = half_btf(cospi[32], bf0[6], cospi[32], bf0[7], cos_bit);
997
1.99M
  bf1[7] = half_btf(cospi[32], bf0[6], -cospi[32], bf0[7], cos_bit);
998
1.99M
  bf1[8] = bf0[8];
999
1.99M
  bf1[9] = bf0[9];
1000
1.99M
  bf1[10] = half_btf(cospi[32], bf0[10], cospi[32], bf0[11], cos_bit);
1001
1.99M
  bf1[11] = half_btf(cospi[32], bf0[10], -cospi[32], bf0[11], cos_bit);
1002
1.99M
  bf1[12] = bf0[12];
1003
1.99M
  bf1[13] = bf0[13];
1004
1.99M
  bf1[14] = half_btf(cospi[32], bf0[14], cospi[32], bf0[15], cos_bit);
1005
1.99M
  bf1[15] = half_btf(cospi[32], bf0[14], -cospi[32], bf0[15], cos_bit);
1006
1.99M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
1007
1008
  // stage 9
1009
1.99M
  bf0 = step;
1010
1.99M
  bf1 = output;
1011
1.99M
  bf1[0] = bf0[0];
1012
1.99M
  bf1[1] = -bf0[8];
1013
1.99M
  bf1[2] = bf0[12];
1014
1.99M
  bf1[3] = -bf0[4];
1015
1.99M
  bf1[4] = bf0[6];
1016
1.99M
  bf1[5] = -bf0[14];
1017
1.99M
  bf1[6] = bf0[10];
1018
1.99M
  bf1[7] = -bf0[2];
1019
1.99M
  bf1[8] = bf0[3];
1020
1.99M
  bf1[9] = -bf0[11];
1021
1.99M
  bf1[10] = bf0[15];
1022
1.99M
  bf1[11] = -bf0[7];
1023
1.99M
  bf1[12] = bf0[5];
1024
1.99M
  bf1[13] = -bf0[13];
1025
1.99M
  bf1[14] = bf0[9];
1026
1.99M
  bf1[15] = -bf0[1];
1027
1.99M
}
1028
1029
void av1_iidentity4_c(const int32_t *input, int32_t *output, int8_t cos_bit,
1030
295k
                      const int8_t *stage_range) {
1031
295k
  (void)cos_bit;
1032
295k
  (void)stage_range;
1033
1.47M
  for (int i = 0; i < 4; ++i) {
1034
1.18M
    output[i] = round_shift((int64_t)NewSqrt2 * input[i], NewSqrt2Bits);
1035
1.18M
  }
1036
295k
  assert(stage_range[0] + NewSqrt2Bits <= 32);
1037
295k
}
1038
1039
void av1_iidentity8_c(const int32_t *input, int32_t *output, int8_t cos_bit,
1040
397k
                      const int8_t *stage_range) {
1041
397k
  (void)cos_bit;
1042
397k
  (void)stage_range;
1043
3.57M
  for (int i = 0; i < 8; ++i) output[i] = (int32_t)((int64_t)input[i] * 2);
1044
397k
}
1045
1046
void av1_iidentity16_c(const int32_t *input, int32_t *output, int8_t cos_bit,
1047
216k
                       const int8_t *stage_range) {
1048
216k
  (void)cos_bit;
1049
216k
  (void)stage_range;
1050
3.67M
  for (int i = 0; i < 16; ++i)
1051
3.45M
    output[i] = round_shift((int64_t)NewSqrt2 * 2 * input[i], NewSqrt2Bits);
1052
216k
  assert(stage_range[0] + NewSqrt2Bits <= 32);
1053
216k
}
1054
1055
void av1_iidentity32_c(const int32_t *input, int32_t *output, int8_t cos_bit,
1056
1.87k
                       const int8_t *stage_range) {
1057
1.87k
  (void)cos_bit;
1058
1.87k
  (void)stage_range;
1059
61.7k
  for (int i = 0; i < 32; ++i) output[i] = (int32_t)((int64_t)input[i] * 4);
1060
1.87k
}
1061
1062
void av1_idct64(const int32_t *input, int32_t *output, int8_t cos_bit,
1063
10.4M
                const int8_t *stage_range) {
1064
10.4M
  assert(output != input);
1065
10.4M
  const int32_t size = 64;
1066
10.4M
  const int32_t *cospi = cospi_arr(cos_bit);
1067
1068
10.4M
  int32_t stage = 0;
1069
10.4M
  int32_t *bf0, *bf1;
1070
10.4M
  int32_t step[64];
1071
1072
  // stage 0;
1073
1074
  // stage 1;
1075
10.4M
  stage++;
1076
10.4M
  bf1 = output;
1077
10.4M
  bf1[0] = input[0];
1078
10.4M
  bf1[1] = input[32];
1079
10.4M
  bf1[2] = input[16];
1080
10.4M
  bf1[3] = input[48];
1081
10.4M
  bf1[4] = input[8];
1082
10.4M
  bf1[5] = input[40];
1083
10.4M
  bf1[6] = input[24];
1084
10.4M
  bf1[7] = input[56];
1085
10.4M
  bf1[8] = input[4];
1086
10.4M
  bf1[9] = input[36];
1087
10.4M
  bf1[10] = input[20];
1088
10.4M
  bf1[11] = input[52];
1089
10.4M
  bf1[12] = input[12];
1090
10.4M
  bf1[13] = input[44];
1091
10.4M
  bf1[14] = input[28];
1092
10.4M
  bf1[15] = input[60];
1093
10.4M
  bf1[16] = input[2];
1094
10.4M
  bf1[17] = input[34];
1095
10.4M
  bf1[18] = input[18];
1096
10.4M
  bf1[19] = input[50];
1097
10.4M
  bf1[20] = input[10];
1098
10.4M
  bf1[21] = input[42];
1099
10.4M
  bf1[22] = input[26];
1100
10.4M
  bf1[23] = input[58];
1101
10.4M
  bf1[24] = input[6];
1102
10.4M
  bf1[25] = input[38];
1103
10.4M
  bf1[26] = input[22];
1104
10.4M
  bf1[27] = input[54];
1105
10.4M
  bf1[28] = input[14];
1106
10.4M
  bf1[29] = input[46];
1107
10.4M
  bf1[30] = input[30];
1108
10.4M
  bf1[31] = input[62];
1109
10.4M
  bf1[32] = input[1];
1110
10.4M
  bf1[33] = input[33];
1111
10.4M
  bf1[34] = input[17];
1112
10.4M
  bf1[35] = input[49];
1113
10.4M
  bf1[36] = input[9];
1114
10.4M
  bf1[37] = input[41];
1115
10.4M
  bf1[38] = input[25];
1116
10.4M
  bf1[39] = input[57];
1117
10.4M
  bf1[40] = input[5];
1118
10.4M
  bf1[41] = input[37];
1119
10.4M
  bf1[42] = input[21];
1120
10.4M
  bf1[43] = input[53];
1121
10.4M
  bf1[44] = input[13];
1122
10.4M
  bf1[45] = input[45];
1123
10.4M
  bf1[46] = input[29];
1124
10.4M
  bf1[47] = input[61];
1125
10.4M
  bf1[48] = input[3];
1126
10.4M
  bf1[49] = input[35];
1127
10.4M
  bf1[50] = input[19];
1128
10.4M
  bf1[51] = input[51];
1129
10.4M
  bf1[52] = input[11];
1130
10.4M
  bf1[53] = input[43];
1131
10.4M
  bf1[54] = input[27];
1132
10.4M
  bf1[55] = input[59];
1133
10.4M
  bf1[56] = input[7];
1134
10.4M
  bf1[57] = input[39];
1135
10.4M
  bf1[58] = input[23];
1136
10.4M
  bf1[59] = input[55];
1137
10.4M
  bf1[60] = input[15];
1138
10.4M
  bf1[61] = input[47];
1139
10.4M
  bf1[62] = input[31];
1140
10.4M
  bf1[63] = input[63];
1141
10.4M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
1142
1143
  // stage 2
1144
10.4M
  stage++;
1145
10.4M
  bf0 = output;
1146
10.4M
  bf1 = step;
1147
10.4M
  bf1[0] = bf0[0];
1148
10.4M
  bf1[1] = bf0[1];
1149
10.4M
  bf1[2] = bf0[2];
1150
10.4M
  bf1[3] = bf0[3];
1151
10.4M
  bf1[4] = bf0[4];
1152
10.4M
  bf1[5] = bf0[5];
1153
10.4M
  bf1[6] = bf0[6];
1154
10.4M
  bf1[7] = bf0[7];
1155
10.4M
  bf1[8] = bf0[8];
1156
10.4M
  bf1[9] = bf0[9];
1157
10.4M
  bf1[10] = bf0[10];
1158
10.4M
  bf1[11] = bf0[11];
1159
10.4M
  bf1[12] = bf0[12];
1160
10.4M
  bf1[13] = bf0[13];
1161
10.4M
  bf1[14] = bf0[14];
1162
10.4M
  bf1[15] = bf0[15];
1163
10.4M
  bf1[16] = bf0[16];
1164
10.4M
  bf1[17] = bf0[17];
1165
10.4M
  bf1[18] = bf0[18];
1166
10.4M
  bf1[19] = bf0[19];
1167
10.4M
  bf1[20] = bf0[20];
1168
10.4M
  bf1[21] = bf0[21];
1169
10.4M
  bf1[22] = bf0[22];
1170
10.4M
  bf1[23] = bf0[23];
1171
10.4M
  bf1[24] = bf0[24];
1172
10.4M
  bf1[25] = bf0[25];
1173
10.4M
  bf1[26] = bf0[26];
1174
10.4M
  bf1[27] = bf0[27];
1175
10.4M
  bf1[28] = bf0[28];
1176
10.4M
  bf1[29] = bf0[29];
1177
10.4M
  bf1[30] = bf0[30];
1178
10.4M
  bf1[31] = bf0[31];
1179
10.4M
  bf1[32] = half_btf(cospi[63], bf0[32], -cospi[1], bf0[63], cos_bit);
1180
10.4M
  bf1[33] = half_btf(cospi[31], bf0[33], -cospi[33], bf0[62], cos_bit);
1181
10.4M
  bf1[34] = half_btf(cospi[47], bf0[34], -cospi[17], bf0[61], cos_bit);
1182
10.4M
  bf1[35] = half_btf(cospi[15], bf0[35], -cospi[49], bf0[60], cos_bit);
1183
10.4M
  bf1[36] = half_btf(cospi[55], bf0[36], -cospi[9], bf0[59], cos_bit);
1184
10.4M
  bf1[37] = half_btf(cospi[23], bf0[37], -cospi[41], bf0[58], cos_bit);
1185
10.4M
  bf1[38] = half_btf(cospi[39], bf0[38], -cospi[25], bf0[57], cos_bit);
1186
10.4M
  bf1[39] = half_btf(cospi[7], bf0[39], -cospi[57], bf0[56], cos_bit);
1187
10.4M
  bf1[40] = half_btf(cospi[59], bf0[40], -cospi[5], bf0[55], cos_bit);
1188
10.4M
  bf1[41] = half_btf(cospi[27], bf0[41], -cospi[37], bf0[54], cos_bit);
1189
10.4M
  bf1[42] = half_btf(cospi[43], bf0[42], -cospi[21], bf0[53], cos_bit);
1190
10.4M
  bf1[43] = half_btf(cospi[11], bf0[43], -cospi[53], bf0[52], cos_bit);
1191
10.4M
  bf1[44] = half_btf(cospi[51], bf0[44], -cospi[13], bf0[51], cos_bit);
1192
10.4M
  bf1[45] = half_btf(cospi[19], bf0[45], -cospi[45], bf0[50], cos_bit);
1193
10.4M
  bf1[46] = half_btf(cospi[35], bf0[46], -cospi[29], bf0[49], cos_bit);
1194
10.4M
  bf1[47] = half_btf(cospi[3], bf0[47], -cospi[61], bf0[48], cos_bit);
1195
10.4M
  bf1[48] = half_btf(cospi[61], bf0[47], cospi[3], bf0[48], cos_bit);
1196
10.4M
  bf1[49] = half_btf(cospi[29], bf0[46], cospi[35], bf0[49], cos_bit);
1197
10.4M
  bf1[50] = half_btf(cospi[45], bf0[45], cospi[19], bf0[50], cos_bit);
1198
10.4M
  bf1[51] = half_btf(cospi[13], bf0[44], cospi[51], bf0[51], cos_bit);
1199
10.4M
  bf1[52] = half_btf(cospi[53], bf0[43], cospi[11], bf0[52], cos_bit);
1200
10.4M
  bf1[53] = half_btf(cospi[21], bf0[42], cospi[43], bf0[53], cos_bit);
1201
10.4M
  bf1[54] = half_btf(cospi[37], bf0[41], cospi[27], bf0[54], cos_bit);
1202
10.4M
  bf1[55] = half_btf(cospi[5], bf0[40], cospi[59], bf0[55], cos_bit);
1203
10.4M
  bf1[56] = half_btf(cospi[57], bf0[39], cospi[7], bf0[56], cos_bit);
1204
10.4M
  bf1[57] = half_btf(cospi[25], bf0[38], cospi[39], bf0[57], cos_bit);
1205
10.4M
  bf1[58] = half_btf(cospi[41], bf0[37], cospi[23], bf0[58], cos_bit);
1206
10.4M
  bf1[59] = half_btf(cospi[9], bf0[36], cospi[55], bf0[59], cos_bit);
1207
10.4M
  bf1[60] = half_btf(cospi[49], bf0[35], cospi[15], bf0[60], cos_bit);
1208
10.4M
  bf1[61] = half_btf(cospi[17], bf0[34], cospi[47], bf0[61], cos_bit);
1209
10.4M
  bf1[62] = half_btf(cospi[33], bf0[33], cospi[31], bf0[62], cos_bit);
1210
10.4M
  bf1[63] = half_btf(cospi[1], bf0[32], cospi[63], bf0[63], cos_bit);
1211
10.4M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
1212
1213
  // stage 3
1214
10.4M
  stage++;
1215
10.4M
  bf0 = step;
1216
10.4M
  bf1 = output;
1217
10.4M
  bf1[0] = bf0[0];
1218
10.4M
  bf1[1] = bf0[1];
1219
10.4M
  bf1[2] = bf0[2];
1220
10.4M
  bf1[3] = bf0[3];
1221
10.4M
  bf1[4] = bf0[4];
1222
10.4M
  bf1[5] = bf0[5];
1223
10.4M
  bf1[6] = bf0[6];
1224
10.4M
  bf1[7] = bf0[7];
1225
10.4M
  bf1[8] = bf0[8];
1226
10.4M
  bf1[9] = bf0[9];
1227
10.4M
  bf1[10] = bf0[10];
1228
10.4M
  bf1[11] = bf0[11];
1229
10.4M
  bf1[12] = bf0[12];
1230
10.4M
  bf1[13] = bf0[13];
1231
10.4M
  bf1[14] = bf0[14];
1232
10.4M
  bf1[15] = bf0[15];
1233
10.4M
  bf1[16] = half_btf(cospi[62], bf0[16], -cospi[2], bf0[31], cos_bit);
1234
10.4M
  bf1[17] = half_btf(cospi[30], bf0[17], -cospi[34], bf0[30], cos_bit);
1235
10.4M
  bf1[18] = half_btf(cospi[46], bf0[18], -cospi[18], bf0[29], cos_bit);
1236
10.4M
  bf1[19] = half_btf(cospi[14], bf0[19], -cospi[50], bf0[28], cos_bit);
1237
10.4M
  bf1[20] = half_btf(cospi[54], bf0[20], -cospi[10], bf0[27], cos_bit);
1238
10.4M
  bf1[21] = half_btf(cospi[22], bf0[21], -cospi[42], bf0[26], cos_bit);
1239
10.4M
  bf1[22] = half_btf(cospi[38], bf0[22], -cospi[26], bf0[25], cos_bit);
1240
10.4M
  bf1[23] = half_btf(cospi[6], bf0[23], -cospi[58], bf0[24], cos_bit);
1241
10.4M
  bf1[24] = half_btf(cospi[58], bf0[23], cospi[6], bf0[24], cos_bit);
1242
10.4M
  bf1[25] = half_btf(cospi[26], bf0[22], cospi[38], bf0[25], cos_bit);
1243
10.4M
  bf1[26] = half_btf(cospi[42], bf0[21], cospi[22], bf0[26], cos_bit);
1244
10.4M
  bf1[27] = half_btf(cospi[10], bf0[20], cospi[54], bf0[27], cos_bit);
1245
10.4M
  bf1[28] = half_btf(cospi[50], bf0[19], cospi[14], bf0[28], cos_bit);
1246
10.4M
  bf1[29] = half_btf(cospi[18], bf0[18], cospi[46], bf0[29], cos_bit);
1247
10.4M
  bf1[30] = half_btf(cospi[34], bf0[17], cospi[30], bf0[30], cos_bit);
1248
10.4M
  bf1[31] = half_btf(cospi[2], bf0[16], cospi[62], bf0[31], cos_bit);
1249
10.4M
  bf1[32] = clamp_value(bf0[32] + bf0[33], stage_range[stage]);
1250
10.4M
  bf1[33] = clamp_value(bf0[32] - bf0[33], stage_range[stage]);
1251
10.4M
  bf1[34] = clamp_value(-bf0[34] + bf0[35], stage_range[stage]);
1252
10.4M
  bf1[35] = clamp_value(bf0[34] + bf0[35], stage_range[stage]);
1253
10.4M
  bf1[36] = clamp_value(bf0[36] + bf0[37], stage_range[stage]);
1254
10.4M
  bf1[37] = clamp_value(bf0[36] - bf0[37], stage_range[stage]);
1255
10.4M
  bf1[38] = clamp_value(-bf0[38] + bf0[39], stage_range[stage]);
1256
10.4M
  bf1[39] = clamp_value(bf0[38] + bf0[39], stage_range[stage]);
1257
10.4M
  bf1[40] = clamp_value(bf0[40] + bf0[41], stage_range[stage]);
1258
10.4M
  bf1[41] = clamp_value(bf0[40] - bf0[41], stage_range[stage]);
1259
10.4M
  bf1[42] = clamp_value(-bf0[42] + bf0[43], stage_range[stage]);
1260
10.4M
  bf1[43] = clamp_value(bf0[42] + bf0[43], stage_range[stage]);
1261
10.4M
  bf1[44] = clamp_value(bf0[44] + bf0[45], stage_range[stage]);
1262
10.4M
  bf1[45] = clamp_value(bf0[44] - bf0[45], stage_range[stage]);
1263
10.4M
  bf1[46] = clamp_value(-bf0[46] + bf0[47], stage_range[stage]);
1264
10.4M
  bf1[47] = clamp_value(bf0[46] + bf0[47], stage_range[stage]);
1265
10.4M
  bf1[48] = clamp_value(bf0[48] + bf0[49], stage_range[stage]);
1266
10.4M
  bf1[49] = clamp_value(bf0[48] - bf0[49], stage_range[stage]);
1267
10.4M
  bf1[50] = clamp_value(-bf0[50] + bf0[51], stage_range[stage]);
1268
10.4M
  bf1[51] = clamp_value(bf0[50] + bf0[51], stage_range[stage]);
1269
10.4M
  bf1[52] = clamp_value(bf0[52] + bf0[53], stage_range[stage]);
1270
10.4M
  bf1[53] = clamp_value(bf0[52] - bf0[53], stage_range[stage]);
1271
10.4M
  bf1[54] = clamp_value(-bf0[54] + bf0[55], stage_range[stage]);
1272
10.4M
  bf1[55] = clamp_value(bf0[54] + bf0[55], stage_range[stage]);
1273
10.4M
  bf1[56] = clamp_value(bf0[56] + bf0[57], stage_range[stage]);
1274
10.4M
  bf1[57] = clamp_value(bf0[56] - bf0[57], stage_range[stage]);
1275
10.4M
  bf1[58] = clamp_value(-bf0[58] + bf0[59], stage_range[stage]);
1276
10.4M
  bf1[59] = clamp_value(bf0[58] + bf0[59], stage_range[stage]);
1277
10.4M
  bf1[60] = clamp_value(bf0[60] + bf0[61], stage_range[stage]);
1278
10.4M
  bf1[61] = clamp_value(bf0[60] - bf0[61], stage_range[stage]);
1279
10.4M
  bf1[62] = clamp_value(-bf0[62] + bf0[63], stage_range[stage]);
1280
10.4M
  bf1[63] = clamp_value(bf0[62] + bf0[63], stage_range[stage]);
1281
10.4M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
1282
1283
  // stage 4
1284
10.4M
  stage++;
1285
10.4M
  bf0 = output;
1286
10.4M
  bf1 = step;
1287
10.4M
  bf1[0] = bf0[0];
1288
10.4M
  bf1[1] = bf0[1];
1289
10.4M
  bf1[2] = bf0[2];
1290
10.4M
  bf1[3] = bf0[3];
1291
10.4M
  bf1[4] = bf0[4];
1292
10.4M
  bf1[5] = bf0[5];
1293
10.4M
  bf1[6] = bf0[6];
1294
10.4M
  bf1[7] = bf0[7];
1295
10.4M
  bf1[8] = half_btf(cospi[60], bf0[8], -cospi[4], bf0[15], cos_bit);
1296
10.4M
  bf1[9] = half_btf(cospi[28], bf0[9], -cospi[36], bf0[14], cos_bit);
1297
10.4M
  bf1[10] = half_btf(cospi[44], bf0[10], -cospi[20], bf0[13], cos_bit);
1298
10.4M
  bf1[11] = half_btf(cospi[12], bf0[11], -cospi[52], bf0[12], cos_bit);
1299
10.4M
  bf1[12] = half_btf(cospi[52], bf0[11], cospi[12], bf0[12], cos_bit);
1300
10.4M
  bf1[13] = half_btf(cospi[20], bf0[10], cospi[44], bf0[13], cos_bit);
1301
10.4M
  bf1[14] = half_btf(cospi[36], bf0[9], cospi[28], bf0[14], cos_bit);
1302
10.4M
  bf1[15] = half_btf(cospi[4], bf0[8], cospi[60], bf0[15], cos_bit);
1303
10.4M
  bf1[16] = clamp_value(bf0[16] + bf0[17], stage_range[stage]);
1304
10.4M
  bf1[17] = clamp_value(bf0[16] - bf0[17], stage_range[stage]);
1305
10.4M
  bf1[18] = clamp_value(-bf0[18] + bf0[19], stage_range[stage]);
1306
10.4M
  bf1[19] = clamp_value(bf0[18] + bf0[19], stage_range[stage]);
1307
10.4M
  bf1[20] = clamp_value(bf0[20] + bf0[21], stage_range[stage]);
1308
10.4M
  bf1[21] = clamp_value(bf0[20] - bf0[21], stage_range[stage]);
1309
10.4M
  bf1[22] = clamp_value(-bf0[22] + bf0[23], stage_range[stage]);
1310
10.4M
  bf1[23] = clamp_value(bf0[22] + bf0[23], stage_range[stage]);
1311
10.4M
  bf1[24] = clamp_value(bf0[24] + bf0[25], stage_range[stage]);
1312
10.4M
  bf1[25] = clamp_value(bf0[24] - bf0[25], stage_range[stage]);
1313
10.4M
  bf1[26] = clamp_value(-bf0[26] + bf0[27], stage_range[stage]);
1314
10.4M
  bf1[27] = clamp_value(bf0[26] + bf0[27], stage_range[stage]);
1315
10.4M
  bf1[28] = clamp_value(bf0[28] + bf0[29], stage_range[stage]);
1316
10.4M
  bf1[29] = clamp_value(bf0[28] - bf0[29], stage_range[stage]);
1317
10.4M
  bf1[30] = clamp_value(-bf0[30] + bf0[31], stage_range[stage]);
1318
10.4M
  bf1[31] = clamp_value(bf0[30] + bf0[31], stage_range[stage]);
1319
10.4M
  bf1[32] = bf0[32];
1320
10.4M
  bf1[33] = half_btf(-cospi[4], bf0[33], cospi[60], bf0[62], cos_bit);
1321
10.4M
  bf1[34] = half_btf(-cospi[60], bf0[34], -cospi[4], bf0[61], cos_bit);
1322
10.4M
  bf1[35] = bf0[35];
1323
10.4M
  bf1[36] = bf0[36];
1324
10.4M
  bf1[37] = half_btf(-cospi[36], bf0[37], cospi[28], bf0[58], cos_bit);
1325
10.4M
  bf1[38] = half_btf(-cospi[28], bf0[38], -cospi[36], bf0[57], cos_bit);
1326
10.4M
  bf1[39] = bf0[39];
1327
10.4M
  bf1[40] = bf0[40];
1328
10.4M
  bf1[41] = half_btf(-cospi[20], bf0[41], cospi[44], bf0[54], cos_bit);
1329
10.4M
  bf1[42] = half_btf(-cospi[44], bf0[42], -cospi[20], bf0[53], cos_bit);
1330
10.4M
  bf1[43] = bf0[43];
1331
10.4M
  bf1[44] = bf0[44];
1332
10.4M
  bf1[45] = half_btf(-cospi[52], bf0[45], cospi[12], bf0[50], cos_bit);
1333
10.4M
  bf1[46] = half_btf(-cospi[12], bf0[46], -cospi[52], bf0[49], cos_bit);
1334
10.4M
  bf1[47] = bf0[47];
1335
10.4M
  bf1[48] = bf0[48];
1336
10.4M
  bf1[49] = half_btf(-cospi[52], bf0[46], cospi[12], bf0[49], cos_bit);
1337
10.4M
  bf1[50] = half_btf(cospi[12], bf0[45], cospi[52], bf0[50], cos_bit);
1338
10.4M
  bf1[51] = bf0[51];
1339
10.4M
  bf1[52] = bf0[52];
1340
10.4M
  bf1[53] = half_btf(-cospi[20], bf0[42], cospi[44], bf0[53], cos_bit);
1341
10.4M
  bf1[54] = half_btf(cospi[44], bf0[41], cospi[20], bf0[54], cos_bit);
1342
10.4M
  bf1[55] = bf0[55];
1343
10.4M
  bf1[56] = bf0[56];
1344
10.4M
  bf1[57] = half_btf(-cospi[36], bf0[38], cospi[28], bf0[57], cos_bit);
1345
10.4M
  bf1[58] = half_btf(cospi[28], bf0[37], cospi[36], bf0[58], cos_bit);
1346
10.4M
  bf1[59] = bf0[59];
1347
10.4M
  bf1[60] = bf0[60];
1348
10.4M
  bf1[61] = half_btf(-cospi[4], bf0[34], cospi[60], bf0[61], cos_bit);
1349
10.4M
  bf1[62] = half_btf(cospi[60], bf0[33], cospi[4], bf0[62], cos_bit);
1350
10.4M
  bf1[63] = bf0[63];
1351
10.4M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
1352
1353
  // stage 5
1354
10.4M
  stage++;
1355
10.4M
  bf0 = step;
1356
10.4M
  bf1 = output;
1357
10.4M
  bf1[0] = bf0[0];
1358
10.4M
  bf1[1] = bf0[1];
1359
10.4M
  bf1[2] = bf0[2];
1360
10.4M
  bf1[3] = bf0[3];
1361
10.4M
  bf1[4] = half_btf(cospi[56], bf0[4], -cospi[8], bf0[7], cos_bit);
1362
10.4M
  bf1[5] = half_btf(cospi[24], bf0[5], -cospi[40], bf0[6], cos_bit);
1363
10.4M
  bf1[6] = half_btf(cospi[40], bf0[5], cospi[24], bf0[6], cos_bit);
1364
10.4M
  bf1[7] = half_btf(cospi[8], bf0[4], cospi[56], bf0[7], cos_bit);
1365
10.4M
  bf1[8] = clamp_value(bf0[8] + bf0[9], stage_range[stage]);
1366
10.4M
  bf1[9] = clamp_value(bf0[8] - bf0[9], stage_range[stage]);
1367
10.4M
  bf1[10] = clamp_value(-bf0[10] + bf0[11], stage_range[stage]);
1368
10.4M
  bf1[11] = clamp_value(bf0[10] + bf0[11], stage_range[stage]);
1369
10.4M
  bf1[12] = clamp_value(bf0[12] + bf0[13], stage_range[stage]);
1370
10.4M
  bf1[13] = clamp_value(bf0[12] - bf0[13], stage_range[stage]);
1371
10.4M
  bf1[14] = clamp_value(-bf0[14] + bf0[15], stage_range[stage]);
1372
10.4M
  bf1[15] = clamp_value(bf0[14] + bf0[15], stage_range[stage]);
1373
10.4M
  bf1[16] = bf0[16];
1374
10.4M
  bf1[17] = half_btf(-cospi[8], bf0[17], cospi[56], bf0[30], cos_bit);
1375
10.4M
  bf1[18] = half_btf(-cospi[56], bf0[18], -cospi[8], bf0[29], cos_bit);
1376
10.4M
  bf1[19] = bf0[19];
1377
10.4M
  bf1[20] = bf0[20];
1378
10.4M
  bf1[21] = half_btf(-cospi[40], bf0[21], cospi[24], bf0[26], cos_bit);
1379
10.4M
  bf1[22] = half_btf(-cospi[24], bf0[22], -cospi[40], bf0[25], cos_bit);
1380
10.4M
  bf1[23] = bf0[23];
1381
10.4M
  bf1[24] = bf0[24];
1382
10.4M
  bf1[25] = half_btf(-cospi[40], bf0[22], cospi[24], bf0[25], cos_bit);
1383
10.4M
  bf1[26] = half_btf(cospi[24], bf0[21], cospi[40], bf0[26], cos_bit);
1384
10.4M
  bf1[27] = bf0[27];
1385
10.4M
  bf1[28] = bf0[28];
1386
10.4M
  bf1[29] = half_btf(-cospi[8], bf0[18], cospi[56], bf0[29], cos_bit);
1387
10.4M
  bf1[30] = half_btf(cospi[56], bf0[17], cospi[8], bf0[30], cos_bit);
1388
10.4M
  bf1[31] = bf0[31];
1389
10.4M
  bf1[32] = clamp_value(bf0[32] + bf0[35], stage_range[stage]);
1390
10.4M
  bf1[33] = clamp_value(bf0[33] + bf0[34], stage_range[stage]);
1391
10.4M
  bf1[34] = clamp_value(bf0[33] - bf0[34], stage_range[stage]);
1392
10.4M
  bf1[35] = clamp_value(bf0[32] - bf0[35], stage_range[stage]);
1393
10.4M
  bf1[36] = clamp_value(-bf0[36] + bf0[39], stage_range[stage]);
1394
10.4M
  bf1[37] = clamp_value(-bf0[37] + bf0[38], stage_range[stage]);
1395
10.4M
  bf1[38] = clamp_value(bf0[37] + bf0[38], stage_range[stage]);
1396
10.4M
  bf1[39] = clamp_value(bf0[36] + bf0[39], stage_range[stage]);
1397
10.4M
  bf1[40] = clamp_value(bf0[40] + bf0[43], stage_range[stage]);
1398
10.4M
  bf1[41] = clamp_value(bf0[41] + bf0[42], stage_range[stage]);
1399
10.4M
  bf1[42] = clamp_value(bf0[41] - bf0[42], stage_range[stage]);
1400
10.4M
  bf1[43] = clamp_value(bf0[40] - bf0[43], stage_range[stage]);
1401
10.4M
  bf1[44] = clamp_value(-bf0[44] + bf0[47], stage_range[stage]);
1402
10.4M
  bf1[45] = clamp_value(-bf0[45] + bf0[46], stage_range[stage]);
1403
10.4M
  bf1[46] = clamp_value(bf0[45] + bf0[46], stage_range[stage]);
1404
10.4M
  bf1[47] = clamp_value(bf0[44] + bf0[47], stage_range[stage]);
1405
10.4M
  bf1[48] = clamp_value(bf0[48] + bf0[51], stage_range[stage]);
1406
10.4M
  bf1[49] = clamp_value(bf0[49] + bf0[50], stage_range[stage]);
1407
10.4M
  bf1[50] = clamp_value(bf0[49] - bf0[50], stage_range[stage]);
1408
10.4M
  bf1[51] = clamp_value(bf0[48] - bf0[51], stage_range[stage]);
1409
10.4M
  bf1[52] = clamp_value(-bf0[52] + bf0[55], stage_range[stage]);
1410
10.4M
  bf1[53] = clamp_value(-bf0[53] + bf0[54], stage_range[stage]);
1411
10.4M
  bf1[54] = clamp_value(bf0[53] + bf0[54], stage_range[stage]);
1412
10.4M
  bf1[55] = clamp_value(bf0[52] + bf0[55], stage_range[stage]);
1413
10.4M
  bf1[56] = clamp_value(bf0[56] + bf0[59], stage_range[stage]);
1414
10.4M
  bf1[57] = clamp_value(bf0[57] + bf0[58], stage_range[stage]);
1415
10.4M
  bf1[58] = clamp_value(bf0[57] - bf0[58], stage_range[stage]);
1416
10.4M
  bf1[59] = clamp_value(bf0[56] - bf0[59], stage_range[stage]);
1417
10.4M
  bf1[60] = clamp_value(-bf0[60] + bf0[63], stage_range[stage]);
1418
10.4M
  bf1[61] = clamp_value(-bf0[61] + bf0[62], stage_range[stage]);
1419
10.4M
  bf1[62] = clamp_value(bf0[61] + bf0[62], stage_range[stage]);
1420
10.4M
  bf1[63] = clamp_value(bf0[60] + bf0[63], stage_range[stage]);
1421
10.4M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
1422
1423
  // stage 6
1424
10.4M
  stage++;
1425
10.4M
  bf0 = output;
1426
10.4M
  bf1 = step;
1427
10.4M
  bf1[0] = half_btf(cospi[32], bf0[0], cospi[32], bf0[1], cos_bit);
1428
10.4M
  bf1[1] = half_btf(cospi[32], bf0[0], -cospi[32], bf0[1], cos_bit);
1429
10.4M
  bf1[2] = half_btf(cospi[48], bf0[2], -cospi[16], bf0[3], cos_bit);
1430
10.4M
  bf1[3] = half_btf(cospi[16], bf0[2], cospi[48], bf0[3], cos_bit);
1431
10.4M
  bf1[4] = clamp_value(bf0[4] + bf0[5], stage_range[stage]);
1432
10.4M
  bf1[5] = clamp_value(bf0[4] - bf0[5], stage_range[stage]);
1433
10.4M
  bf1[6] = clamp_value(-bf0[6] + bf0[7], stage_range[stage]);
1434
10.4M
  bf1[7] = clamp_value(bf0[6] + bf0[7], stage_range[stage]);
1435
10.4M
  bf1[8] = bf0[8];
1436
10.4M
  bf1[9] = half_btf(-cospi[16], bf0[9], cospi[48], bf0[14], cos_bit);
1437
10.4M
  bf1[10] = half_btf(-cospi[48], bf0[10], -cospi[16], bf0[13], cos_bit);
1438
10.4M
  bf1[11] = bf0[11];
1439
10.4M
  bf1[12] = bf0[12];
1440
10.4M
  bf1[13] = half_btf(-cospi[16], bf0[10], cospi[48], bf0[13], cos_bit);
1441
10.4M
  bf1[14] = half_btf(cospi[48], bf0[9], cospi[16], bf0[14], cos_bit);
1442
10.4M
  bf1[15] = bf0[15];
1443
10.4M
  bf1[16] = clamp_value(bf0[16] + bf0[19], stage_range[stage]);
1444
10.4M
  bf1[17] = clamp_value(bf0[17] + bf0[18], stage_range[stage]);
1445
10.4M
  bf1[18] = clamp_value(bf0[17] - bf0[18], stage_range[stage]);
1446
10.4M
  bf1[19] = clamp_value(bf0[16] - bf0[19], stage_range[stage]);
1447
10.4M
  bf1[20] = clamp_value(-bf0[20] + bf0[23], stage_range[stage]);
1448
10.4M
  bf1[21] = clamp_value(-bf0[21] + bf0[22], stage_range[stage]);
1449
10.4M
  bf1[22] = clamp_value(bf0[21] + bf0[22], stage_range[stage]);
1450
10.4M
  bf1[23] = clamp_value(bf0[20] + bf0[23], stage_range[stage]);
1451
10.4M
  bf1[24] = clamp_value(bf0[24] + bf0[27], stage_range[stage]);
1452
10.4M
  bf1[25] = clamp_value(bf0[25] + bf0[26], stage_range[stage]);
1453
10.4M
  bf1[26] = clamp_value(bf0[25] - bf0[26], stage_range[stage]);
1454
10.4M
  bf1[27] = clamp_value(bf0[24] - bf0[27], stage_range[stage]);
1455
10.4M
  bf1[28] = clamp_value(-bf0[28] + bf0[31], stage_range[stage]);
1456
10.4M
  bf1[29] = clamp_value(-bf0[29] + bf0[30], stage_range[stage]);
1457
10.4M
  bf1[30] = clamp_value(bf0[29] + bf0[30], stage_range[stage]);
1458
10.4M
  bf1[31] = clamp_value(bf0[28] + bf0[31], stage_range[stage]);
1459
10.4M
  bf1[32] = bf0[32];
1460
10.4M
  bf1[33] = bf0[33];
1461
10.4M
  bf1[34] = half_btf(-cospi[8], bf0[34], cospi[56], bf0[61], cos_bit);
1462
10.4M
  bf1[35] = half_btf(-cospi[8], bf0[35], cospi[56], bf0[60], cos_bit);
1463
10.4M
  bf1[36] = half_btf(-cospi[56], bf0[36], -cospi[8], bf0[59], cos_bit);
1464
10.4M
  bf1[37] = half_btf(-cospi[56], bf0[37], -cospi[8], bf0[58], cos_bit);
1465
10.4M
  bf1[38] = bf0[38];
1466
10.4M
  bf1[39] = bf0[39];
1467
10.4M
  bf1[40] = bf0[40];
1468
10.4M
  bf1[41] = bf0[41];
1469
10.4M
  bf1[42] = half_btf(-cospi[40], bf0[42], cospi[24], bf0[53], cos_bit);
1470
10.4M
  bf1[43] = half_btf(-cospi[40], bf0[43], cospi[24], bf0[52], cos_bit);
1471
10.4M
  bf1[44] = half_btf(-cospi[24], bf0[44], -cospi[40], bf0[51], cos_bit);
1472
10.4M
  bf1[45] = half_btf(-cospi[24], bf0[45], -cospi[40], bf0[50], cos_bit);
1473
10.4M
  bf1[46] = bf0[46];
1474
10.4M
  bf1[47] = bf0[47];
1475
10.4M
  bf1[48] = bf0[48];
1476
10.4M
  bf1[49] = bf0[49];
1477
10.4M
  bf1[50] = half_btf(-cospi[40], bf0[45], cospi[24], bf0[50], cos_bit);
1478
10.4M
  bf1[51] = half_btf(-cospi[40], bf0[44], cospi[24], bf0[51], cos_bit);
1479
10.4M
  bf1[52] = half_btf(cospi[24], bf0[43], cospi[40], bf0[52], cos_bit);
1480
10.4M
  bf1[53] = half_btf(cospi[24], bf0[42], cospi[40], bf0[53], cos_bit);
1481
10.4M
  bf1[54] = bf0[54];
1482
10.4M
  bf1[55] = bf0[55];
1483
10.4M
  bf1[56] = bf0[56];
1484
10.4M
  bf1[57] = bf0[57];
1485
10.4M
  bf1[58] = half_btf(-cospi[8], bf0[37], cospi[56], bf0[58], cos_bit);
1486
10.4M
  bf1[59] = half_btf(-cospi[8], bf0[36], cospi[56], bf0[59], cos_bit);
1487
10.4M
  bf1[60] = half_btf(cospi[56], bf0[35], cospi[8], bf0[60], cos_bit);
1488
10.4M
  bf1[61] = half_btf(cospi[56], bf0[34], cospi[8], bf0[61], cos_bit);
1489
10.4M
  bf1[62] = bf0[62];
1490
10.4M
  bf1[63] = bf0[63];
1491
10.4M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
1492
1493
  // stage 7
1494
10.4M
  stage++;
1495
10.4M
  bf0 = step;
1496
10.4M
  bf1 = output;
1497
10.4M
  bf1[0] = clamp_value(bf0[0] + bf0[3], stage_range[stage]);
1498
10.4M
  bf1[1] = clamp_value(bf0[1] + bf0[2], stage_range[stage]);
1499
10.4M
  bf1[2] = clamp_value(bf0[1] - bf0[2], stage_range[stage]);
1500
10.4M
  bf1[3] = clamp_value(bf0[0] - bf0[3], stage_range[stage]);
1501
10.4M
  bf1[4] = bf0[4];
1502
10.4M
  bf1[5] = half_btf(-cospi[32], bf0[5], cospi[32], bf0[6], cos_bit);
1503
10.4M
  bf1[6] = half_btf(cospi[32], bf0[5], cospi[32], bf0[6], cos_bit);
1504
10.4M
  bf1[7] = bf0[7];
1505
10.4M
  bf1[8] = clamp_value(bf0[8] + bf0[11], stage_range[stage]);
1506
10.4M
  bf1[9] = clamp_value(bf0[9] + bf0[10], stage_range[stage]);
1507
10.4M
  bf1[10] = clamp_value(bf0[9] - bf0[10], stage_range[stage]);
1508
10.4M
  bf1[11] = clamp_value(bf0[8] - bf0[11], stage_range[stage]);
1509
10.4M
  bf1[12] = clamp_value(-bf0[12] + bf0[15], stage_range[stage]);
1510
10.4M
  bf1[13] = clamp_value(-bf0[13] + bf0[14], stage_range[stage]);
1511
10.4M
  bf1[14] = clamp_value(bf0[13] + bf0[14], stage_range[stage]);
1512
10.4M
  bf1[15] = clamp_value(bf0[12] + bf0[15], stage_range[stage]);
1513
10.4M
  bf1[16] = bf0[16];
1514
10.4M
  bf1[17] = bf0[17];
1515
10.4M
  bf1[18] = half_btf(-cospi[16], bf0[18], cospi[48], bf0[29], cos_bit);
1516
10.4M
  bf1[19] = half_btf(-cospi[16], bf0[19], cospi[48], bf0[28], cos_bit);
1517
10.4M
  bf1[20] = half_btf(-cospi[48], bf0[20], -cospi[16], bf0[27], cos_bit);
1518
10.4M
  bf1[21] = half_btf(-cospi[48], bf0[21], -cospi[16], bf0[26], cos_bit);
1519
10.4M
  bf1[22] = bf0[22];
1520
10.4M
  bf1[23] = bf0[23];
1521
10.4M
  bf1[24] = bf0[24];
1522
10.4M
  bf1[25] = bf0[25];
1523
10.4M
  bf1[26] = half_btf(-cospi[16], bf0[21], cospi[48], bf0[26], cos_bit);
1524
10.4M
  bf1[27] = half_btf(-cospi[16], bf0[20], cospi[48], bf0[27], cos_bit);
1525
10.4M
  bf1[28] = half_btf(cospi[48], bf0[19], cospi[16], bf0[28], cos_bit);
1526
10.4M
  bf1[29] = half_btf(cospi[48], bf0[18], cospi[16], bf0[29], cos_bit);
1527
10.4M
  bf1[30] = bf0[30];
1528
10.4M
  bf1[31] = bf0[31];
1529
10.4M
  bf1[32] = clamp_value(bf0[32] + bf0[39], stage_range[stage]);
1530
10.4M
  bf1[33] = clamp_value(bf0[33] + bf0[38], stage_range[stage]);
1531
10.4M
  bf1[34] = clamp_value(bf0[34] + bf0[37], stage_range[stage]);
1532
10.4M
  bf1[35] = clamp_value(bf0[35] + bf0[36], stage_range[stage]);
1533
10.4M
  bf1[36] = clamp_value(bf0[35] - bf0[36], stage_range[stage]);
1534
10.4M
  bf1[37] = clamp_value(bf0[34] - bf0[37], stage_range[stage]);
1535
10.4M
  bf1[38] = clamp_value(bf0[33] - bf0[38], stage_range[stage]);
1536
10.4M
  bf1[39] = clamp_value(bf0[32] - bf0[39], stage_range[stage]);
1537
10.4M
  bf1[40] = clamp_value(-bf0[40] + bf0[47], stage_range[stage]);
1538
10.4M
  bf1[41] = clamp_value(-bf0[41] + bf0[46], stage_range[stage]);
1539
10.4M
  bf1[42] = clamp_value(-bf0[42] + bf0[45], stage_range[stage]);
1540
10.4M
  bf1[43] = clamp_value(-bf0[43] + bf0[44], stage_range[stage]);
1541
10.4M
  bf1[44] = clamp_value(bf0[43] + bf0[44], stage_range[stage]);
1542
10.4M
  bf1[45] = clamp_value(bf0[42] + bf0[45], stage_range[stage]);
1543
10.4M
  bf1[46] = clamp_value(bf0[41] + bf0[46], stage_range[stage]);
1544
10.4M
  bf1[47] = clamp_value(bf0[40] + bf0[47], stage_range[stage]);
1545
10.4M
  bf1[48] = clamp_value(bf0[48] + bf0[55], stage_range[stage]);
1546
10.4M
  bf1[49] = clamp_value(bf0[49] + bf0[54], stage_range[stage]);
1547
10.4M
  bf1[50] = clamp_value(bf0[50] + bf0[53], stage_range[stage]);
1548
10.4M
  bf1[51] = clamp_value(bf0[51] + bf0[52], stage_range[stage]);
1549
10.4M
  bf1[52] = clamp_value(bf0[51] - bf0[52], stage_range[stage]);
1550
10.4M
  bf1[53] = clamp_value(bf0[50] - bf0[53], stage_range[stage]);
1551
10.4M
  bf1[54] = clamp_value(bf0[49] - bf0[54], stage_range[stage]);
1552
10.4M
  bf1[55] = clamp_value(bf0[48] - bf0[55], stage_range[stage]);
1553
10.4M
  bf1[56] = clamp_value(-bf0[56] + bf0[63], stage_range[stage]);
1554
10.4M
  bf1[57] = clamp_value(-bf0[57] + bf0[62], stage_range[stage]);
1555
10.4M
  bf1[58] = clamp_value(-bf0[58] + bf0[61], stage_range[stage]);
1556
10.4M
  bf1[59] = clamp_value(-bf0[59] + bf0[60], stage_range[stage]);
1557
10.4M
  bf1[60] = clamp_value(bf0[59] + bf0[60], stage_range[stage]);
1558
10.4M
  bf1[61] = clamp_value(bf0[58] + bf0[61], stage_range[stage]);
1559
10.4M
  bf1[62] = clamp_value(bf0[57] + bf0[62], stage_range[stage]);
1560
10.4M
  bf1[63] = clamp_value(bf0[56] + bf0[63], stage_range[stage]);
1561
10.4M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
1562
1563
  // stage 8
1564
10.4M
  stage++;
1565
10.4M
  bf0 = output;
1566
10.4M
  bf1 = step;
1567
10.4M
  bf1[0] = clamp_value(bf0[0] + bf0[7], stage_range[stage]);
1568
10.4M
  bf1[1] = clamp_value(bf0[1] + bf0[6], stage_range[stage]);
1569
10.4M
  bf1[2] = clamp_value(bf0[2] + bf0[5], stage_range[stage]);
1570
10.4M
  bf1[3] = clamp_value(bf0[3] + bf0[4], stage_range[stage]);
1571
10.4M
  bf1[4] = clamp_value(bf0[3] - bf0[4], stage_range[stage]);
1572
10.4M
  bf1[5] = clamp_value(bf0[2] - bf0[5], stage_range[stage]);
1573
10.4M
  bf1[6] = clamp_value(bf0[1] - bf0[6], stage_range[stage]);
1574
10.4M
  bf1[7] = clamp_value(bf0[0] - bf0[7], stage_range[stage]);
1575
10.4M
  bf1[8] = bf0[8];
1576
10.4M
  bf1[9] = bf0[9];
1577
10.4M
  bf1[10] = half_btf(-cospi[32], bf0[10], cospi[32], bf0[13], cos_bit);
1578
10.4M
  bf1[11] = half_btf(-cospi[32], bf0[11], cospi[32], bf0[12], cos_bit);
1579
10.4M
  bf1[12] = half_btf(cospi[32], bf0[11], cospi[32], bf0[12], cos_bit);
1580
10.4M
  bf1[13] = half_btf(cospi[32], bf0[10], cospi[32], bf0[13], cos_bit);
1581
10.4M
  bf1[14] = bf0[14];
1582
10.4M
  bf1[15] = bf0[15];
1583
10.4M
  bf1[16] = clamp_value(bf0[16] + bf0[23], stage_range[stage]);
1584
10.4M
  bf1[17] = clamp_value(bf0[17] + bf0[22], stage_range[stage]);
1585
10.4M
  bf1[18] = clamp_value(bf0[18] + bf0[21], stage_range[stage]);
1586
10.4M
  bf1[19] = clamp_value(bf0[19] + bf0[20], stage_range[stage]);
1587
10.4M
  bf1[20] = clamp_value(bf0[19] - bf0[20], stage_range[stage]);
1588
10.4M
  bf1[21] = clamp_value(bf0[18] - bf0[21], stage_range[stage]);
1589
10.4M
  bf1[22] = clamp_value(bf0[17] - bf0[22], stage_range[stage]);
1590
10.4M
  bf1[23] = clamp_value(bf0[16] - bf0[23], stage_range[stage]);
1591
10.4M
  bf1[24] = clamp_value(-bf0[24] + bf0[31], stage_range[stage]);
1592
10.4M
  bf1[25] = clamp_value(-bf0[25] + bf0[30], stage_range[stage]);
1593
10.4M
  bf1[26] = clamp_value(-bf0[26] + bf0[29], stage_range[stage]);
1594
10.4M
  bf1[27] = clamp_value(-bf0[27] + bf0[28], stage_range[stage]);
1595
10.4M
  bf1[28] = clamp_value(bf0[27] + bf0[28], stage_range[stage]);
1596
10.4M
  bf1[29] = clamp_value(bf0[26] + bf0[29], stage_range[stage]);
1597
10.4M
  bf1[30] = clamp_value(bf0[25] + bf0[30], stage_range[stage]);
1598
10.4M
  bf1[31] = clamp_value(bf0[24] + bf0[31], stage_range[stage]);
1599
10.4M
  bf1[32] = bf0[32];
1600
10.4M
  bf1[33] = bf0[33];
1601
10.4M
  bf1[34] = bf0[34];
1602
10.4M
  bf1[35] = bf0[35];
1603
10.4M
  bf1[36] = half_btf(-cospi[16], bf0[36], cospi[48], bf0[59], cos_bit);
1604
10.4M
  bf1[37] = half_btf(-cospi[16], bf0[37], cospi[48], bf0[58], cos_bit);
1605
10.4M
  bf1[38] = half_btf(-cospi[16], bf0[38], cospi[48], bf0[57], cos_bit);
1606
10.4M
  bf1[39] = half_btf(-cospi[16], bf0[39], cospi[48], bf0[56], cos_bit);
1607
10.4M
  bf1[40] = half_btf(-cospi[48], bf0[40], -cospi[16], bf0[55], cos_bit);
1608
10.4M
  bf1[41] = half_btf(-cospi[48], bf0[41], -cospi[16], bf0[54], cos_bit);
1609
10.4M
  bf1[42] = half_btf(-cospi[48], bf0[42], -cospi[16], bf0[53], cos_bit);
1610
10.4M
  bf1[43] = half_btf(-cospi[48], bf0[43], -cospi[16], bf0[52], cos_bit);
1611
10.4M
  bf1[44] = bf0[44];
1612
10.4M
  bf1[45] = bf0[45];
1613
10.4M
  bf1[46] = bf0[46];
1614
10.4M
  bf1[47] = bf0[47];
1615
10.4M
  bf1[48] = bf0[48];
1616
10.4M
  bf1[49] = bf0[49];
1617
10.4M
  bf1[50] = bf0[50];
1618
10.4M
  bf1[51] = bf0[51];
1619
10.4M
  bf1[52] = half_btf(-cospi[16], bf0[43], cospi[48], bf0[52], cos_bit);
1620
10.4M
  bf1[53] = half_btf(-cospi[16], bf0[42], cospi[48], bf0[53], cos_bit);
1621
10.4M
  bf1[54] = half_btf(-cospi[16], bf0[41], cospi[48], bf0[54], cos_bit);
1622
10.4M
  bf1[55] = half_btf(-cospi[16], bf0[40], cospi[48], bf0[55], cos_bit);
1623
10.4M
  bf1[56] = half_btf(cospi[48], bf0[39], cospi[16], bf0[56], cos_bit);
1624
10.4M
  bf1[57] = half_btf(cospi[48], bf0[38], cospi[16], bf0[57], cos_bit);
1625
10.4M
  bf1[58] = half_btf(cospi[48], bf0[37], cospi[16], bf0[58], cos_bit);
1626
10.4M
  bf1[59] = half_btf(cospi[48], bf0[36], cospi[16], bf0[59], cos_bit);
1627
10.4M
  bf1[60] = bf0[60];
1628
10.4M
  bf1[61] = bf0[61];
1629
10.4M
  bf1[62] = bf0[62];
1630
10.4M
  bf1[63] = bf0[63];
1631
10.4M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
1632
1633
  // stage 9
1634
10.4M
  stage++;
1635
10.4M
  bf0 = step;
1636
10.4M
  bf1 = output;
1637
10.4M
  bf1[0] = clamp_value(bf0[0] + bf0[15], stage_range[stage]);
1638
10.4M
  bf1[1] = clamp_value(bf0[1] + bf0[14], stage_range[stage]);
1639
10.4M
  bf1[2] = clamp_value(bf0[2] + bf0[13], stage_range[stage]);
1640
10.4M
  bf1[3] = clamp_value(bf0[3] + bf0[12], stage_range[stage]);
1641
10.4M
  bf1[4] = clamp_value(bf0[4] + bf0[11], stage_range[stage]);
1642
10.4M
  bf1[5] = clamp_value(bf0[5] + bf0[10], stage_range[stage]);
1643
10.4M
  bf1[6] = clamp_value(bf0[6] + bf0[9], stage_range[stage]);
1644
10.4M
  bf1[7] = clamp_value(bf0[7] + bf0[8], stage_range[stage]);
1645
10.4M
  bf1[8] = clamp_value(bf0[7] - bf0[8], stage_range[stage]);
1646
10.4M
  bf1[9] = clamp_value(bf0[6] - bf0[9], stage_range[stage]);
1647
10.4M
  bf1[10] = clamp_value(bf0[5] - bf0[10], stage_range[stage]);
1648
10.4M
  bf1[11] = clamp_value(bf0[4] - bf0[11], stage_range[stage]);
1649
10.4M
  bf1[12] = clamp_value(bf0[3] - bf0[12], stage_range[stage]);
1650
10.4M
  bf1[13] = clamp_value(bf0[2] - bf0[13], stage_range[stage]);
1651
10.4M
  bf1[14] = clamp_value(bf0[1] - bf0[14], stage_range[stage]);
1652
10.4M
  bf1[15] = clamp_value(bf0[0] - bf0[15], stage_range[stage]);
1653
10.4M
  bf1[16] = bf0[16];
1654
10.4M
  bf1[17] = bf0[17];
1655
10.4M
  bf1[18] = bf0[18];
1656
10.4M
  bf1[19] = bf0[19];
1657
10.4M
  bf1[20] = half_btf(-cospi[32], bf0[20], cospi[32], bf0[27], cos_bit);
1658
10.4M
  bf1[21] = half_btf(-cospi[32], bf0[21], cospi[32], bf0[26], cos_bit);
1659
10.4M
  bf1[22] = half_btf(-cospi[32], bf0[22], cospi[32], bf0[25], cos_bit);
1660
10.4M
  bf1[23] = half_btf(-cospi[32], bf0[23], cospi[32], bf0[24], cos_bit);
1661
10.4M
  bf1[24] = half_btf(cospi[32], bf0[23], cospi[32], bf0[24], cos_bit);
1662
10.4M
  bf1[25] = half_btf(cospi[32], bf0[22], cospi[32], bf0[25], cos_bit);
1663
10.4M
  bf1[26] = half_btf(cospi[32], bf0[21], cospi[32], bf0[26], cos_bit);
1664
10.4M
  bf1[27] = half_btf(cospi[32], bf0[20], cospi[32], bf0[27], cos_bit);
1665
10.4M
  bf1[28] = bf0[28];
1666
10.4M
  bf1[29] = bf0[29];
1667
10.4M
  bf1[30] = bf0[30];
1668
10.4M
  bf1[31] = bf0[31];
1669
10.4M
  bf1[32] = clamp_value(bf0[32] + bf0[47], stage_range[stage]);
1670
10.4M
  bf1[33] = clamp_value(bf0[33] + bf0[46], stage_range[stage]);
1671
10.4M
  bf1[34] = clamp_value(bf0[34] + bf0[45], stage_range[stage]);
1672
10.4M
  bf1[35] = clamp_value(bf0[35] + bf0[44], stage_range[stage]);
1673
10.4M
  bf1[36] = clamp_value(bf0[36] + bf0[43], stage_range[stage]);
1674
10.4M
  bf1[37] = clamp_value(bf0[37] + bf0[42], stage_range[stage]);
1675
10.4M
  bf1[38] = clamp_value(bf0[38] + bf0[41], stage_range[stage]);
1676
10.4M
  bf1[39] = clamp_value(bf0[39] + bf0[40], stage_range[stage]);
1677
10.4M
  bf1[40] = clamp_value(bf0[39] - bf0[40], stage_range[stage]);
1678
10.4M
  bf1[41] = clamp_value(bf0[38] - bf0[41], stage_range[stage]);
1679
10.4M
  bf1[42] = clamp_value(bf0[37] - bf0[42], stage_range[stage]);
1680
10.4M
  bf1[43] = clamp_value(bf0[36] - bf0[43], stage_range[stage]);
1681
10.4M
  bf1[44] = clamp_value(bf0[35] - bf0[44], stage_range[stage]);
1682
10.4M
  bf1[45] = clamp_value(bf0[34] - bf0[45], stage_range[stage]);
1683
10.4M
  bf1[46] = clamp_value(bf0[33] - bf0[46], stage_range[stage]);
1684
10.4M
  bf1[47] = clamp_value(bf0[32] - bf0[47], stage_range[stage]);
1685
10.4M
  bf1[48] = clamp_value(-bf0[48] + bf0[63], stage_range[stage]);
1686
10.4M
  bf1[49] = clamp_value(-bf0[49] + bf0[62], stage_range[stage]);
1687
10.4M
  bf1[50] = clamp_value(-bf0[50] + bf0[61], stage_range[stage]);
1688
10.4M
  bf1[51] = clamp_value(-bf0[51] + bf0[60], stage_range[stage]);
1689
10.4M
  bf1[52] = clamp_value(-bf0[52] + bf0[59], stage_range[stage]);
1690
10.4M
  bf1[53] = clamp_value(-bf0[53] + bf0[58], stage_range[stage]);
1691
10.4M
  bf1[54] = clamp_value(-bf0[54] + bf0[57], stage_range[stage]);
1692
10.4M
  bf1[55] = clamp_value(-bf0[55] + bf0[56], stage_range[stage]);
1693
10.4M
  bf1[56] = clamp_value(bf0[55] + bf0[56], stage_range[stage]);
1694
10.4M
  bf1[57] = clamp_value(bf0[54] + bf0[57], stage_range[stage]);
1695
10.4M
  bf1[58] = clamp_value(bf0[53] + bf0[58], stage_range[stage]);
1696
10.4M
  bf1[59] = clamp_value(bf0[52] + bf0[59], stage_range[stage]);
1697
10.4M
  bf1[60] = clamp_value(bf0[51] + bf0[60], stage_range[stage]);
1698
10.4M
  bf1[61] = clamp_value(bf0[50] + bf0[61], stage_range[stage]);
1699
10.4M
  bf1[62] = clamp_value(bf0[49] + bf0[62], stage_range[stage]);
1700
10.4M
  bf1[63] = clamp_value(bf0[48] + bf0[63], stage_range[stage]);
1701
10.4M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
1702
1703
  // stage 10
1704
10.4M
  stage++;
1705
10.4M
  bf0 = output;
1706
10.4M
  bf1 = step;
1707
10.4M
  bf1[0] = clamp_value(bf0[0] + bf0[31], stage_range[stage]);
1708
10.4M
  bf1[1] = clamp_value(bf0[1] + bf0[30], stage_range[stage]);
1709
10.4M
  bf1[2] = clamp_value(bf0[2] + bf0[29], stage_range[stage]);
1710
10.4M
  bf1[3] = clamp_value(bf0[3] + bf0[28], stage_range[stage]);
1711
10.4M
  bf1[4] = clamp_value(bf0[4] + bf0[27], stage_range[stage]);
1712
10.4M
  bf1[5] = clamp_value(bf0[5] + bf0[26], stage_range[stage]);
1713
10.4M
  bf1[6] = clamp_value(bf0[6] + bf0[25], stage_range[stage]);
1714
10.4M
  bf1[7] = clamp_value(bf0[7] + bf0[24], stage_range[stage]);
1715
10.4M
  bf1[8] = clamp_value(bf0[8] + bf0[23], stage_range[stage]);
1716
10.4M
  bf1[9] = clamp_value(bf0[9] + bf0[22], stage_range[stage]);
1717
10.4M
  bf1[10] = clamp_value(bf0[10] + bf0[21], stage_range[stage]);
1718
10.4M
  bf1[11] = clamp_value(bf0[11] + bf0[20], stage_range[stage]);
1719
10.4M
  bf1[12] = clamp_value(bf0[12] + bf0[19], stage_range[stage]);
1720
10.4M
  bf1[13] = clamp_value(bf0[13] + bf0[18], stage_range[stage]);
1721
10.4M
  bf1[14] = clamp_value(bf0[14] + bf0[17], stage_range[stage]);
1722
10.4M
  bf1[15] = clamp_value(bf0[15] + bf0[16], stage_range[stage]);
1723
10.4M
  bf1[16] = clamp_value(bf0[15] - bf0[16], stage_range[stage]);
1724
10.4M
  bf1[17] = clamp_value(bf0[14] - bf0[17], stage_range[stage]);
1725
10.4M
  bf1[18] = clamp_value(bf0[13] - bf0[18], stage_range[stage]);
1726
10.4M
  bf1[19] = clamp_value(bf0[12] - bf0[19], stage_range[stage]);
1727
10.4M
  bf1[20] = clamp_value(bf0[11] - bf0[20], stage_range[stage]);
1728
10.4M
  bf1[21] = clamp_value(bf0[10] - bf0[21], stage_range[stage]);
1729
10.4M
  bf1[22] = clamp_value(bf0[9] - bf0[22], stage_range[stage]);
1730
10.4M
  bf1[23] = clamp_value(bf0[8] - bf0[23], stage_range[stage]);
1731
10.4M
  bf1[24] = clamp_value(bf0[7] - bf0[24], stage_range[stage]);
1732
10.4M
  bf1[25] = clamp_value(bf0[6] - bf0[25], stage_range[stage]);
1733
10.4M
  bf1[26] = clamp_value(bf0[5] - bf0[26], stage_range[stage]);
1734
10.4M
  bf1[27] = clamp_value(bf0[4] - bf0[27], stage_range[stage]);
1735
10.4M
  bf1[28] = clamp_value(bf0[3] - bf0[28], stage_range[stage]);
1736
10.4M
  bf1[29] = clamp_value(bf0[2] - bf0[29], stage_range[stage]);
1737
10.4M
  bf1[30] = clamp_value(bf0[1] - bf0[30], stage_range[stage]);
1738
10.4M
  bf1[31] = clamp_value(bf0[0] - bf0[31], stage_range[stage]);
1739
10.4M
  bf1[32] = bf0[32];
1740
10.4M
  bf1[33] = bf0[33];
1741
10.4M
  bf1[34] = bf0[34];
1742
10.4M
  bf1[35] = bf0[35];
1743
10.4M
  bf1[36] = bf0[36];
1744
10.4M
  bf1[37] = bf0[37];
1745
10.4M
  bf1[38] = bf0[38];
1746
10.4M
  bf1[39] = bf0[39];
1747
10.4M
  bf1[40] = half_btf(-cospi[32], bf0[40], cospi[32], bf0[55], cos_bit);
1748
10.4M
  bf1[41] = half_btf(-cospi[32], bf0[41], cospi[32], bf0[54], cos_bit);
1749
10.4M
  bf1[42] = half_btf(-cospi[32], bf0[42], cospi[32], bf0[53], cos_bit);
1750
10.4M
  bf1[43] = half_btf(-cospi[32], bf0[43], cospi[32], bf0[52], cos_bit);
1751
10.4M
  bf1[44] = half_btf(-cospi[32], bf0[44], cospi[32], bf0[51], cos_bit);
1752
10.4M
  bf1[45] = half_btf(-cospi[32], bf0[45], cospi[32], bf0[50], cos_bit);
1753
10.4M
  bf1[46] = half_btf(-cospi[32], bf0[46], cospi[32], bf0[49], cos_bit);
1754
10.4M
  bf1[47] = half_btf(-cospi[32], bf0[47], cospi[32], bf0[48], cos_bit);
1755
10.4M
  bf1[48] = half_btf(cospi[32], bf0[47], cospi[32], bf0[48], cos_bit);
1756
10.4M
  bf1[49] = half_btf(cospi[32], bf0[46], cospi[32], bf0[49], cos_bit);
1757
10.4M
  bf1[50] = half_btf(cospi[32], bf0[45], cospi[32], bf0[50], cos_bit);
1758
10.4M
  bf1[51] = half_btf(cospi[32], bf0[44], cospi[32], bf0[51], cos_bit);
1759
10.4M
  bf1[52] = half_btf(cospi[32], bf0[43], cospi[32], bf0[52], cos_bit);
1760
10.4M
  bf1[53] = half_btf(cospi[32], bf0[42], cospi[32], bf0[53], cos_bit);
1761
10.4M
  bf1[54] = half_btf(cospi[32], bf0[41], cospi[32], bf0[54], cos_bit);
1762
10.4M
  bf1[55] = half_btf(cospi[32], bf0[40], cospi[32], bf0[55], cos_bit);
1763
10.4M
  bf1[56] = bf0[56];
1764
10.4M
  bf1[57] = bf0[57];
1765
10.4M
  bf1[58] = bf0[58];
1766
10.4M
  bf1[59] = bf0[59];
1767
10.4M
  bf1[60] = bf0[60];
1768
10.4M
  bf1[61] = bf0[61];
1769
10.4M
  bf1[62] = bf0[62];
1770
10.4M
  bf1[63] = bf0[63];
1771
10.4M
  av1_range_check_buf(stage, input, bf1, size, stage_range[stage]);
1772
1773
  // stage 11
1774
10.4M
  stage++;
1775
10.4M
  bf0 = step;
1776
10.4M
  bf1 = output;
1777
10.4M
  bf1[0] = clamp_value(bf0[0] + bf0[63], stage_range[stage]);
1778
10.4M
  bf1[1] = clamp_value(bf0[1] + bf0[62], stage_range[stage]);
1779
10.4M
  bf1[2] = clamp_value(bf0[2] + bf0[61], stage_range[stage]);
1780
10.4M
  bf1[3] = clamp_value(bf0[3] + bf0[60], stage_range[stage]);
1781
10.4M
  bf1[4] = clamp_value(bf0[4] + bf0[59], stage_range[stage]);
1782
10.4M
  bf1[5] = clamp_value(bf0[5] + bf0[58], stage_range[stage]);
1783
10.4M
  bf1[6] = clamp_value(bf0[6] + bf0[57], stage_range[stage]);
1784
10.4M
  bf1[7] = clamp_value(bf0[7] + bf0[56], stage_range[stage]);
1785
10.4M
  bf1[8] = clamp_value(bf0[8] + bf0[55], stage_range[stage]);
1786
10.4M
  bf1[9] = clamp_value(bf0[9] + bf0[54], stage_range[stage]);
1787
10.4M
  bf1[10] = clamp_value(bf0[10] + bf0[53], stage_range[stage]);
1788
10.4M
  bf1[11] = clamp_value(bf0[11] + bf0[52], stage_range[stage]);
1789
10.4M
  bf1[12] = clamp_value(bf0[12] + bf0[51], stage_range[stage]);
1790
10.4M
  bf1[13] = clamp_value(bf0[13] + bf0[50], stage_range[stage]);
1791
10.4M
  bf1[14] = clamp_value(bf0[14] + bf0[49], stage_range[stage]);
1792
10.4M
  bf1[15] = clamp_value(bf0[15] + bf0[48], stage_range[stage]);
1793
10.4M
  bf1[16] = clamp_value(bf0[16] + bf0[47], stage_range[stage]);
1794
10.4M
  bf1[17] = clamp_value(bf0[17] + bf0[46], stage_range[stage]);
1795
10.4M
  bf1[18] = clamp_value(bf0[18] + bf0[45], stage_range[stage]);
1796
10.4M
  bf1[19] = clamp_value(bf0[19] + bf0[44], stage_range[stage]);
1797
10.4M
  bf1[20] = clamp_value(bf0[20] + bf0[43], stage_range[stage]);
1798
10.4M
  bf1[21] = clamp_value(bf0[21] + bf0[42], stage_range[stage]);
1799
10.4M
  bf1[22] = clamp_value(bf0[22] + bf0[41], stage_range[stage]);
1800
10.4M
  bf1[23] = clamp_value(bf0[23] + bf0[40], stage_range[stage]);
1801
10.4M
  bf1[24] = clamp_value(bf0[24] + bf0[39], stage_range[stage]);
1802
10.4M
  bf1[25] = clamp_value(bf0[25] + bf0[38], stage_range[stage]);
1803
10.4M
  bf1[26] = clamp_value(bf0[26] + bf0[37], stage_range[stage]);
1804
10.4M
  bf1[27] = clamp_value(bf0[27] + bf0[36], stage_range[stage]);
1805
10.4M
  bf1[28] = clamp_value(bf0[28] + bf0[35], stage_range[stage]);
1806
10.4M
  bf1[29] = clamp_value(bf0[29] + bf0[34], stage_range[stage]);
1807
10.4M
  bf1[30] = clamp_value(bf0[30] + bf0[33], stage_range[stage]);
1808
10.4M
  bf1[31] = clamp_value(bf0[31] + bf0[32], stage_range[stage]);
1809
10.4M
  bf1[32] = clamp_value(bf0[31] - bf0[32], stage_range[stage]);
1810
10.4M
  bf1[33] = clamp_value(bf0[30] - bf0[33], stage_range[stage]);
1811
10.4M
  bf1[34] = clamp_value(bf0[29] - bf0[34], stage_range[stage]);
1812
10.4M
  bf1[35] = clamp_value(bf0[28] - bf0[35], stage_range[stage]);
1813
10.4M
  bf1[36] = clamp_value(bf0[27] - bf0[36], stage_range[stage]);
1814
10.4M
  bf1[37] = clamp_value(bf0[26] - bf0[37], stage_range[stage]);
1815
10.4M
  bf1[38] = clamp_value(bf0[25] - bf0[38], stage_range[stage]);
1816
10.4M
  bf1[39] = clamp_value(bf0[24] - bf0[39], stage_range[stage]);
1817
10.4M
  bf1[40] = clamp_value(bf0[23] - bf0[40], stage_range[stage]);
1818
10.4M
  bf1[41] = clamp_value(bf0[22] - bf0[41], stage_range[stage]);
1819
10.4M
  bf1[42] = clamp_value(bf0[21] - bf0[42], stage_range[stage]);
1820
10.4M
  bf1[43] = clamp_value(bf0[20] - bf0[43], stage_range[stage]);
1821
10.4M
  bf1[44] = clamp_value(bf0[19] - bf0[44], stage_range[stage]);
1822
10.4M
  bf1[45] = clamp_value(bf0[18] - bf0[45], stage_range[stage]);
1823
10.4M
  bf1[46] = clamp_value(bf0[17] - bf0[46], stage_range[stage]);
1824
10.4M
  bf1[47] = clamp_value(bf0[16] - bf0[47], stage_range[stage]);
1825
10.4M
  bf1[48] = clamp_value(bf0[15] - bf0[48], stage_range[stage]);
1826
10.4M
  bf1[49] = clamp_value(bf0[14] - bf0[49], stage_range[stage]);
1827
10.4M
  bf1[50] = clamp_value(bf0[13] - bf0[50], stage_range[stage]);
1828
10.4M
  bf1[51] = clamp_value(bf0[12] - bf0[51], stage_range[stage]);
1829
10.4M
  bf1[52] = clamp_value(bf0[11] - bf0[52], stage_range[stage]);
1830
10.4M
  bf1[53] = clamp_value(bf0[10] - bf0[53], stage_range[stage]);
1831
10.4M
  bf1[54] = clamp_value(bf0[9] - bf0[54], stage_range[stage]);
1832
10.4M
  bf1[55] = clamp_value(bf0[8] - bf0[55], stage_range[stage]);
1833
10.4M
  bf1[56] = clamp_value(bf0[7] - bf0[56], stage_range[stage]);
1834
10.4M
  bf1[57] = clamp_value(bf0[6] - bf0[57], stage_range[stage]);
1835
10.4M
  bf1[58] = clamp_value(bf0[5] - bf0[58], stage_range[stage]);
1836
10.4M
  bf1[59] = clamp_value(bf0[4] - bf0[59], stage_range[stage]);
1837
10.4M
  bf1[60] = clamp_value(bf0[3] - bf0[60], stage_range[stage]);
1838
10.4M
  bf1[61] = clamp_value(bf0[2] - bf0[61], stage_range[stage]);
1839
10.4M
  bf1[62] = clamp_value(bf0[1] - bf0[62], stage_range[stage]);
1840
10.4M
  bf1[63] = clamp_value(bf0[0] - bf0[63], stage_range[stage]);
1841
10.4M
}