summaryrefslogtreecommitdiff
path: root/libavcodec/dct32_template.c (plain)
blob: bb3f1553acec555b5f3b5f9075d2ebf358252516
1/*
2 * Template for the Discrete Cosine Transform for 32 samples
3 * Copyright (c) 2001, 2002 Fabrice Bellard
4 *
5 * This file is part of FFmpeg.
6 *
7 * FFmpeg is free software; you can redistribute it and/or
8 * modify it under the terms of the GNU Lesser General Public
9 * License as published by the Free Software Foundation; either
10 * version 2.1 of the License, or (at your option) any later version.
11 *
12 * FFmpeg is distributed in the hope that it will be useful,
13 * but WITHOUT ANY WARRANTY; without even the implied warranty of
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
15 * Lesser General Public License for more details.
16 *
17 * You should have received a copy of the GNU Lesser General Public
18 * License along with FFmpeg; if not, write to the Free Software
19 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
20 */
21
22#include "dct32.h"
23#include "mathops.h"
24#include "libavutil/internal.h"
25
26#ifdef CHECKED
27#define SUINT int
28#define SUINT32 int32_t
29#else
30#define SUINT unsigned
31#define SUINT32 uint32_t
32#endif
33
34#if DCT32_FLOAT
35# define dct32 ff_dct32_float
36# define FIXHR(x) ((float)(x))
37# define MULH3(x, y, s) ((s)*(y)*(x))
38# define INTFLOAT float
39# define SUINTFLOAT float
40#else
41# define dct32 ff_dct32_fixed
42# define FIXHR(a) ((int)((a) * (1LL<<32) + 0.5))
43# define MULH3(x, y, s) MULH((s)*(x), y)
44# define INTFLOAT int
45# define SUINTFLOAT SUINT
46#endif
47
48
49/* tab[i][j] = 1.0 / (2.0 * cos(pi*(2*k+1) / 2^(6 - j))) */
50
51/* cos(i*pi/64) */
52
53#define COS0_0 FIXHR(0.50060299823519630134/2)
54#define COS0_1 FIXHR(0.50547095989754365998/2)
55#define COS0_2 FIXHR(0.51544730992262454697/2)
56#define COS0_3 FIXHR(0.53104259108978417447/2)
57#define COS0_4 FIXHR(0.55310389603444452782/2)
58#define COS0_5 FIXHR(0.58293496820613387367/2)
59#define COS0_6 FIXHR(0.62250412303566481615/2)
60#define COS0_7 FIXHR(0.67480834145500574602/2)
61#define COS0_8 FIXHR(0.74453627100229844977/2)
62#define COS0_9 FIXHR(0.83934964541552703873/2)
63#define COS0_10 FIXHR(0.97256823786196069369/2)
64#define COS0_11 FIXHR(1.16943993343288495515/4)
65#define COS0_12 FIXHR(1.48416461631416627724/4)
66#define COS0_13 FIXHR(2.05778100995341155085/8)
67#define COS0_14 FIXHR(3.40760841846871878570/8)
68#define COS0_15 FIXHR(10.19000812354805681150/32)
69
70#define COS1_0 FIXHR(0.50241928618815570551/2)
71#define COS1_1 FIXHR(0.52249861493968888062/2)
72#define COS1_2 FIXHR(0.56694403481635770368/2)
73#define COS1_3 FIXHR(0.64682178335999012954/2)
74#define COS1_4 FIXHR(0.78815462345125022473/2)
75#define COS1_5 FIXHR(1.06067768599034747134/4)
76#define COS1_6 FIXHR(1.72244709823833392782/4)
77#define COS1_7 FIXHR(5.10114861868916385802/16)
78
79#define COS2_0 FIXHR(0.50979557910415916894/2)
80#define COS2_1 FIXHR(0.60134488693504528054/2)
81#define COS2_2 FIXHR(0.89997622313641570463/2)
82#define COS2_3 FIXHR(2.56291544774150617881/8)
83
84#define COS3_0 FIXHR(0.54119610014619698439/2)
85#define COS3_1 FIXHR(1.30656296487637652785/4)
86
87#define COS4_0 FIXHR(M_SQRT1_2/2)
88
89/* butterfly operator */
90#define BF(a, b, c, s)\
91{\
92 tmp0 = val##a + val##b;\
93 tmp1 = val##a - val##b;\
94 val##a = tmp0;\
95 val##b = MULH3(tmp1, c, 1<<(s));\
96}
97
98#define BF0(a, b, c, s)\
99{\
100 tmp0 = tab[a] + tab[b];\
101 tmp1 = tab[a] - tab[b];\
102 val##a = tmp0;\
103 val##b = MULH3(tmp1, c, 1<<(s));\
104}
105
106#define BF1(a, b, c, d)\
107{\
108 BF(a, b, COS4_0, 1);\
109 BF(c, d,-COS4_0, 1);\
110 val##c += val##d;\
111}
112
113#define BF2(a, b, c, d)\
114{\
115 BF(a, b, COS4_0, 1);\
116 BF(c, d,-COS4_0, 1);\
117 val##c += val##d;\
118 val##a += val##c;\
119 val##c += val##b;\
120 val##b += val##d;\
121}
122
123#define ADD(a, b) val##a += val##b
124
125/* DCT32 without 1/sqrt(2) coef zero scaling. */
126void dct32(INTFLOAT *out, const INTFLOAT *tab)
127{
128 SUINTFLOAT tmp0, tmp1;
129
130 SUINTFLOAT val0 , val1 , val2 , val3 , val4 , val5 , val6 , val7 ,
131 val8 , val9 , val10, val11, val12, val13, val14, val15,
132 val16, val17, val18, val19, val20, val21, val22, val23,
133 val24, val25, val26, val27, val28, val29, val30, val31;
134
135 /* pass 1 */
136 BF0( 0, 31, COS0_0 , 1);
137 BF0(15, 16, COS0_15, 5);
138 /* pass 2 */
139 BF( 0, 15, COS1_0 , 1);
140 BF(16, 31,-COS1_0 , 1);
141 /* pass 1 */
142 BF0( 7, 24, COS0_7 , 1);
143 BF0( 8, 23, COS0_8 , 1);
144 /* pass 2 */
145 BF( 7, 8, COS1_7 , 4);
146 BF(23, 24,-COS1_7 , 4);
147 /* pass 3 */
148 BF( 0, 7, COS2_0 , 1);
149 BF( 8, 15,-COS2_0 , 1);
150 BF(16, 23, COS2_0 , 1);
151 BF(24, 31,-COS2_0 , 1);
152 /* pass 1 */
153 BF0( 3, 28, COS0_3 , 1);
154 BF0(12, 19, COS0_12, 2);
155 /* pass 2 */
156 BF( 3, 12, COS1_3 , 1);
157 BF(19, 28,-COS1_3 , 1);
158 /* pass 1 */
159 BF0( 4, 27, COS0_4 , 1);
160 BF0(11, 20, COS0_11, 2);
161 /* pass 2 */
162 BF( 4, 11, COS1_4 , 1);
163 BF(20, 27,-COS1_4 , 1);
164 /* pass 3 */
165 BF( 3, 4, COS2_3 , 3);
166 BF(11, 12,-COS2_3 , 3);
167 BF(19, 20, COS2_3 , 3);
168 BF(27, 28,-COS2_3 , 3);
169 /* pass 4 */
170 BF( 0, 3, COS3_0 , 1);
171 BF( 4, 7,-COS3_0 , 1);
172 BF( 8, 11, COS3_0 , 1);
173 BF(12, 15,-COS3_0 , 1);
174 BF(16, 19, COS3_0 , 1);
175 BF(20, 23,-COS3_0 , 1);
176 BF(24, 27, COS3_0 , 1);
177 BF(28, 31,-COS3_0 , 1);
178
179
180
181 /* pass 1 */
182 BF0( 1, 30, COS0_1 , 1);
183 BF0(14, 17, COS0_14, 3);
184 /* pass 2 */
185 BF( 1, 14, COS1_1 , 1);
186 BF(17, 30,-COS1_1 , 1);
187 /* pass 1 */
188 BF0( 6, 25, COS0_6 , 1);
189 BF0( 9, 22, COS0_9 , 1);
190 /* pass 2 */
191 BF( 6, 9, COS1_6 , 2);
192 BF(22, 25,-COS1_6 , 2);
193 /* pass 3 */
194 BF( 1, 6, COS2_1 , 1);
195 BF( 9, 14,-COS2_1 , 1);
196 BF(17, 22, COS2_1 , 1);
197 BF(25, 30,-COS2_1 , 1);
198
199 /* pass 1 */
200 BF0( 2, 29, COS0_2 , 1);
201 BF0(13, 18, COS0_13, 3);
202 /* pass 2 */
203 BF( 2, 13, COS1_2 , 1);
204 BF(18, 29,-COS1_2 , 1);
205 /* pass 1 */
206 BF0( 5, 26, COS0_5 , 1);
207 BF0(10, 21, COS0_10, 1);
208 /* pass 2 */
209 BF( 5, 10, COS1_5 , 2);
210 BF(21, 26,-COS1_5 , 2);
211 /* pass 3 */
212 BF( 2, 5, COS2_2 , 1);
213 BF(10, 13,-COS2_2 , 1);
214 BF(18, 21, COS2_2 , 1);
215 BF(26, 29,-COS2_2 , 1);
216 /* pass 4 */
217 BF( 1, 2, COS3_1 , 2);
218 BF( 5, 6,-COS3_1 , 2);
219 BF( 9, 10, COS3_1 , 2);
220 BF(13, 14,-COS3_1 , 2);
221 BF(17, 18, COS3_1 , 2);
222 BF(21, 22,-COS3_1 , 2);
223 BF(25, 26, COS3_1 , 2);
224 BF(29, 30,-COS3_1 , 2);
225
226 /* pass 5 */
227 BF1( 0, 1, 2, 3);
228 BF2( 4, 5, 6, 7);
229 BF1( 8, 9, 10, 11);
230 BF2(12, 13, 14, 15);
231 BF1(16, 17, 18, 19);
232 BF2(20, 21, 22, 23);
233 BF1(24, 25, 26, 27);
234 BF2(28, 29, 30, 31);
235
236 /* pass 6 */
237
238 ADD( 8, 12);
239 ADD(12, 10);
240 ADD(10, 14);
241 ADD(14, 9);
242 ADD( 9, 13);
243 ADD(13, 11);
244 ADD(11, 15);
245
246 out[ 0] = val0;
247 out[16] = val1;
248 out[ 8] = val2;
249 out[24] = val3;
250 out[ 4] = val4;
251 out[20] = val5;
252 out[12] = val6;
253 out[28] = val7;
254 out[ 2] = val8;
255 out[18] = val9;
256 out[10] = val10;
257 out[26] = val11;
258 out[ 6] = val12;
259 out[22] = val13;
260 out[14] = val14;
261 out[30] = val15;
262
263 ADD(24, 28);
264 ADD(28, 26);
265 ADD(26, 30);
266 ADD(30, 25);
267 ADD(25, 29);
268 ADD(29, 27);
269 ADD(27, 31);
270
271 out[ 1] = val16 + val24;
272 out[17] = val17 + val25;
273 out[ 9] = val18 + val26;
274 out[25] = val19 + val27;
275 out[ 5] = val20 + val28;
276 out[21] = val21 + val29;
277 out[13] = val22 + val30;
278 out[29] = val23 + val31;
279 out[ 3] = val24 + val20;
280 out[19] = val25 + val21;
281 out[11] = val26 + val22;
282 out[27] = val27 + val23;
283 out[ 7] = val28 + val18;
284 out[23] = val29 + val19;
285 out[15] = val30 + val17;
286 out[31] = val31;
287}
288
generated by cgit at 2024-05-17 04:45:56 (GMT)