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