Coverage Report

Created: 2022-08-24 06:15

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