1 /* ------------------------------------------------------------------
2 * Copyright (C) 1998-2009 PacketVideo
3 *
4 * Licensed under the Apache License, Version 2.0 (the "License");
5 * you may not use this file except in compliance with the License.
6 * You may obtain a copy of the License at
7 *
8 * http://www.apache.org/licenses/LICENSE-2.0
9 *
10 * Unless required by applicable law or agreed to in writing, software
11 * distributed under the License is distributed on an "AS IS" BASIS,
12 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either
13 * express or implied.
14 * See the License for the specific language governing permissions
15 * and limitations under the License.
16 * -------------------------------------------------------------------
17 */
18 /*
19 ------------------------------------------------------------------------------
20
21 PacketVideo Corp.
22 MP3 Decoder Library
23
24 Filename: mdct_18.cpp
25
26 Date: 09/21/2007
27
28 ------------------------------------------------------------------------------
29 REVISION HISTORY
30
31
32 Description:
33
34 ------------------------------------------------------------------------------
35 INPUT AND OUTPUT DEFINITIONS
36
37 Input
38 int32 vec[], input vector of length 18
39 int32 *history input for overlap and add, vector updated with
40 next overlap and add values
41 const int32 *window sine window used in the mdct, three types are allowed
42 noraml, start and stop
43 Returns
44 none mdct computation in-place
45
46
47 ------------------------------------------------------------------------------
48 FUNCTION DESCRIPTION
49
50 Returns the mdct of length 18 of the input vector, as well as the overlap
51 vector for next iteration ( on history[])
52
53 ------------------------------------------------------------------------------
54 REQUIREMENTS
55
56
57 ------------------------------------------------------------------------------
58 REFERENCES
59
60 ------------------------------------------------------------------------------
61 PSEUDO-CODE
62
63 ------------------------------------------------------------------------------
64 */
65
66 #if ( !defined(PV_ARM_GCC_V5) && !defined(PV_ARM_GCC_V4) && !defined(PV_ARM_V5) && !defined(PV_ARM_V4) )
67 /*----------------------------------------------------------------------------
68 ; INCLUDES
69 ----------------------------------------------------------------------------*/
70
71 #include "pv_mp3dec_fxd_op.h"
72 #include "pvmp3_mdct_18.h"
73
74
75 /*----------------------------------------------------------------------------
76 ; MACROS
77 ; Define module specific macros here
78 ----------------------------------------------------------------------------*/
79
80
81 /*----------------------------------------------------------------------------
82 ; DEFINES
83 ; Include all pre-processor statements here. Include conditional
84 ; compile variables also.
85 ----------------------------------------------------------------------------*/
86
87 /*----------------------------------------------------------------------------
88 ; LOCAL FUNCTION DEFINITIONS
89 ; Function Prototype declaration
90 ----------------------------------------------------------------------------*/
91
92 /*----------------------------------------------------------------------------
93 ; LOCAL STORE/BUFFER/POINTER DEFINITIONS
94 ; Variable declaration - defined here and used outside this module
95 ----------------------------------------------------------------------------*/
96 const int32 cosTerms_dct18[9] =
97 {
98 Qfmt(0.50190991877167f), Qfmt(0.51763809020504f), Qfmt(0.55168895948125f),
99 Qfmt(0.61038729438073f), Qfmt(0.70710678118655f), Qfmt(0.87172339781055f),
100 Qfmt(1.18310079157625f), Qfmt(1.93185165257814f), Qfmt(5.73685662283493f)
101 };
102
103
104 const int32 cosTerms_1_ov_cos_phi[18] =
105 {
106
107 Qfmt1(0.50047634258166f), Qfmt1(0.50431448029008f), Qfmt1(0.51213975715725f),
108 Qfmt1(0.52426456257041f), Qfmt1(0.54119610014620f), Qfmt1(0.56369097343317f),
109 Qfmt1(0.59284452371708f), Qfmt1(0.63023620700513f), Qfmt1(0.67817085245463f),
110
111 Qfmt2(0.74009361646113f), Qfmt2(0.82133981585229f), Qfmt2(0.93057949835179f),
112 Qfmt2(1.08284028510010f), Qfmt2(1.30656296487638f), Qfmt2(1.66275476171152f),
113 Qfmt2(2.31011315767265f), Qfmt2(3.83064878777019f), Qfmt2(11.46279281302667f)
114 };
115
116 /*----------------------------------------------------------------------------
117 ; EXTERNAL FUNCTION REFERENCES
118 ; Declare functions defined elsewhere and referenced in this module
119 ----------------------------------------------------------------------------*/
120
121 /*----------------------------------------------------------------------------
122 ; EXTERNAL GLOBAL STORE/BUFFER/POINTER REFERENCES
123 ; Declare variables used in this module but defined elsewhere
124 ----------------------------------------------------------------------------*/
125
126 /*----------------------------------------------------------------------------
127 ; FUNCTION CODE
128 ----------------------------------------------------------------------------*/
129
130
131
132 __attribute__((no_sanitize("integer")))
pvmp3_mdct_18(int32 vec[],int32 * history,const int32 * window)133 void pvmp3_mdct_18(int32 vec[], int32 *history, const int32 *window)
134 {
135 int32 i;
136 int32 tmp;
137 int32 tmp1;
138 int32 tmp2;
139 int32 tmp3;
140 int32 tmp4;
141
142
143
144 const int32 *pt_cos_split = cosTerms_dct18;
145 const int32 *pt_cos = cosTerms_1_ov_cos_phi;
146 const int32 *pt_cos_x = &cosTerms_1_ov_cos_phi[17];
147 int32 *pt_vec = vec;
148 int32 *pt_vec_o = &vec[17];
149
150
151 for (i = 9; i != 0; i--)
152 {
153 tmp = *(pt_vec);
154 tmp1 = *(pt_vec_o);
155 tmp = fxp_mul32_Q32(tmp << 1, *(pt_cos++));
156 tmp1 = fxp_mul32_Q27(tmp1, *(pt_cos_x--));
157 *(pt_vec++) = tmp + tmp1 ;
158 *(pt_vec_o--) = fxp_mul32_Q28((tmp - tmp1), *(pt_cos_split++));
159 }
160
161
162 pvmp3_dct_9(vec); // Even terms
163 pvmp3_dct_9(&vec[9]); // Odd terms
164
165
166 tmp3 = vec[16]; //
167 vec[16] = vec[ 8];
168 tmp4 = vec[14]; //
169 vec[14] = vec[ 7];
170 tmp = vec[12];
171 vec[12] = vec[ 6];
172 tmp2 = vec[10]; // vec[10]
173 vec[10] = vec[ 5];
174 vec[ 8] = vec[ 4];
175 vec[ 6] = vec[ 3];
176 vec[ 4] = vec[ 2];
177 vec[ 2] = vec[ 1];
178 vec[ 1] = vec[ 9] - tmp2; // vec[9] + vec[10]
179 vec[ 3] = vec[11] - tmp2;
180 vec[ 5] = vec[11] - tmp;
181 vec[ 7] = vec[13] - tmp;
182 vec[ 9] = vec[13] - tmp4;
183 vec[11] = vec[15] - tmp4;
184 vec[13] = vec[15] - tmp3;
185 vec[15] = vec[17] - tmp3;
186
187
188 /* overlap and add */
189
190 tmp2 = vec[0];
191 tmp3 = vec[9];
192
193 for (i = 0; i < 6; i++)
194 {
195 tmp = history[ i];
196 tmp4 = vec[i+10];
197 vec[i+10] = tmp3 + tmp4;
198 tmp1 = vec[i+1];
199 vec[ i] = fxp_mac32_Q32(tmp, (vec[i+10]), window[ i]);
200 tmp3 = tmp4;
201 history[i ] = -(tmp2 + tmp1);
202 tmp2 = tmp1;
203 }
204
205 tmp = history[ 6];
206 tmp4 = vec[16];
207 vec[16] = tmp3 + tmp4;
208 tmp1 = vec[7];
209 vec[ 6] = fxp_mac32_Q32(tmp, vec[16] << 1, window[ i]);
210 tmp = history[ 7];
211 history[6] = -(tmp2 + tmp1);
212 history[7] = -(tmp1 + vec[8]);
213
214 tmp1 = history[ 8];
215 tmp4 = vec[17] + tmp4;
216 vec[ 7] = fxp_mac32_Q32(tmp, tmp4 << 1, window[ 7]);
217 history[8] = -(vec[8] + vec[9]);
218 vec[ 8] = fxp_mac32_Q32(tmp1, vec[17] << 1, window[ 8]);
219
220 tmp = history[9];
221 tmp1 = history[17];
222 tmp2 = history[16];
223 vec[ 9] = fxp_mac32_Q32(tmp, vec[17] << 1, window[ 9]);
224
225 vec[17] = fxp_mac32_Q32(tmp1, vec[10] << 1, window[17]);
226 vec[10] = -vec[ 16];
227 vec[16] = fxp_mac32_Q32(tmp2, vec[11] << 1, window[16]);
228 tmp1 = history[15];
229 tmp2 = history[14];
230 vec[11] = -vec[ 15];
231 vec[15] = fxp_mac32_Q32(tmp1, vec[12] << 1, window[15]);
232 vec[12] = -vec[ 14];
233 vec[14] = fxp_mac32_Q32(tmp2, vec[13] << 1, window[14]);
234
235 tmp = history[13];
236 tmp1 = history[12];
237 tmp2 = history[11];
238 tmp3 = history[10];
239 vec[13] = fxp_mac32_Q32(tmp, vec[12] << 1, window[13]);
240 vec[12] = fxp_mac32_Q32(tmp1, vec[11] << 1, window[12]);
241 vec[11] = fxp_mac32_Q32(tmp2, vec[10] << 1, window[11]);
242 vec[10] = fxp_mac32_Q32(tmp3, tmp4 << 1, window[10]);
243
244
245 /* next iteration overlap */
246
247 tmp1 = history[ 8];
248 tmp3 = history[ 7];
249 tmp2 = history[ 1];
250 tmp = history[ 0];
251 tmp1 <<= 1;
252 tmp3 <<= 1;
253
254 history[ 0] = fxp_mul32_Q32(tmp1, window[18]);
255 history[17] = fxp_mul32_Q32(tmp1, window[35]);
256 history[ 1] = fxp_mul32_Q32(tmp3, window[19]);
257 history[16] = fxp_mul32_Q32(tmp3, window[34]);
258
259 tmp2 <<= 1;
260 tmp <<= 1;
261 history[ 7] = fxp_mul32_Q32(tmp2, window[25]);
262 history[10] = fxp_mul32_Q32(tmp2, window[28]);
263 history[ 8] = fxp_mul32_Q32(tmp, window[26]);
264 history[ 9] = fxp_mul32_Q32(tmp, window[27]);
265
266 tmp1 = history[ 6];
267 tmp3 = history[ 5];
268 tmp4 = history[ 4];
269 tmp2 = history[ 3];
270 tmp = history[ 2];
271
272 tmp1 <<= 1;
273 tmp3 <<= 1;
274 tmp4 <<= 1;
275
276 history[ 2] = fxp_mul32_Q32(tmp1, window[20]);
277 history[15] = fxp_mul32_Q32(tmp1, window[33]);
278 history[ 3] = fxp_mul32_Q32(tmp3, window[21]);
279 history[14] = fxp_mul32_Q32(tmp3, window[32]);
280 history[ 4] = fxp_mul32_Q32(tmp4, window[22]);
281 history[13] = fxp_mul32_Q32(tmp4, window[31]);
282 tmp2 <<= 1;
283 tmp <<= 1;
284 history[ 5] = fxp_mul32_Q32(tmp2, window[23]);
285 history[12] = fxp_mul32_Q32(tmp2, window[30]);
286 history[ 6] = fxp_mul32_Q32(tmp, window[24]);
287 history[11] = fxp_mul32_Q32(tmp, window[29]);
288 }
289
290 #endif // If not assembly
291