1
2 /* -----------------------------------------------------------------------------------------------------------
3 Software License for The Fraunhofer FDK AAC Codec Library for Android
4
5 � Copyright 1995 - 2012 Fraunhofer-Gesellschaft zur F�rderung der angewandten Forschung e.V.
6 All rights reserved.
7
8 1. INTRODUCTION
9 The Fraunhofer FDK AAC Codec Library for Android ("FDK AAC Codec") is software that implements
10 the MPEG Advanced Audio Coding ("AAC") encoding and decoding scheme for digital audio.
11 This FDK AAC Codec software is intended to be used on a wide variety of Android devices.
12
13 AAC's HE-AAC and HE-AAC v2 versions are regarded as today's most efficient general perceptual
14 audio codecs. AAC-ELD is considered the best-performing full-bandwidth communications codec by
15 independent studies and is widely deployed. AAC has been standardized by ISO and IEC as part
16 of the MPEG specifications.
17
18 Patent licenses for necessary patent claims for the FDK AAC Codec (including those of Fraunhofer)
19 may be obtained through Via Licensing (www.vialicensing.com) or through the respective patent owners
20 individually for the purpose of encoding or decoding bit streams in products that are compliant with
21 the ISO/IEC MPEG audio standards. Please note that most manufacturers of Android devices already license
22 these patent claims through Via Licensing or directly from the patent owners, and therefore FDK AAC Codec
23 software may already be covered under those patent licenses when it is used for those licensed purposes only.
24
25 Commercially-licensed AAC software libraries, including floating-point versions with enhanced sound quality,
26 are also available from Fraunhofer. Users are encouraged to check the Fraunhofer website for additional
27 applications information and documentation.
28
29 2. COPYRIGHT LICENSE
30
31 Redistribution and use in source and binary forms, with or without modification, are permitted without
32 payment of copyright license fees provided that you satisfy the following conditions:
33
34 You must retain the complete text of this software license in redistributions of the FDK AAC Codec or
35 your modifications thereto in source code form.
36
37 You must retain the complete text of this software license in the documentation and/or other materials
38 provided with redistributions of the FDK AAC Codec or your modifications thereto in binary form.
39 You must make available free of charge copies of the complete source code of the FDK AAC Codec and your
40 modifications thereto to recipients of copies in binary form.
41
42 The name of Fraunhofer may not be used to endorse or promote products derived from this library without
43 prior written permission.
44
45 You may not charge copyright license fees for anyone to use, copy or distribute the FDK AAC Codec
46 software or your modifications thereto.
47
48 Your modified versions of the FDK AAC Codec must carry prominent notices stating that you changed the software
49 and the date of any change. For modified versions of the FDK AAC Codec, the term
50 "Fraunhofer FDK AAC Codec Library for Android" must be replaced by the term
51 "Third-Party Modified Version of the Fraunhofer FDK AAC Codec Library for Android."
52
53 3. NO PATENT LICENSE
54
55 NO EXPRESS OR IMPLIED LICENSES TO ANY PATENT CLAIMS, including without limitation the patents of Fraunhofer,
56 ARE GRANTED BY THIS SOFTWARE LICENSE. Fraunhofer provides no warranty of patent non-infringement with
57 respect to this software.
58
59 You may use this FDK AAC Codec software or modifications thereto only for purposes that are authorized
60 by appropriate patent licenses.
61
62 4. DISCLAIMER
63
64 This FDK AAC Codec software is provided by Fraunhofer on behalf of the copyright holders and contributors
65 "AS IS" and WITHOUT ANY EXPRESS OR IMPLIED WARRANTIES, including but not limited to the implied warranties
66 of merchantability and fitness for a particular purpose. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR
67 CONTRIBUTORS BE LIABLE for any direct, indirect, incidental, special, exemplary, or consequential damages,
68 including but not limited to procurement of substitute goods or services; loss of use, data, or profits,
69 or business interruption, however caused and on any theory of liability, whether in contract, strict
70 liability, or tort (including negligence), arising in any way out of the use of this software, even if
71 advised of the possibility of such damage.
72
73 5. CONTACT INFORMATION
74
75 Fraunhofer Institute for Integrated Circuits IIS
76 Attention: Audio and Multimedia Departments - FDK AAC LL
77 Am Wolfsmantel 33
78 91058 Erlangen, Germany
79
80 www.iis.fraunhofer.de/amm
81 amm-info@iis.fraunhofer.de
82 ----------------------------------------------------------------------------------------------------------- */
83
84 /***************************** MPEG-4 AAC Decoder **************************
85
86 Author(s): Christian Griebel
87 Description: Dynamic range control (DRC) decoder tool for SBR
88
89 ******************************************************************************/
90
91 #include "sbrdec_drc.h"
92
93
94 /* DRC - Offset table for QMF interpolation. */
95 static const int offsetTab[2][16] =
96 {
97 { 0, 4, 8, 12, 16, 20, 24, 28, 0, 0, 0, 0, 0, 0, 0, 0 }, /* 1024 framing */
98 { 0, 4, 8, 12, 16, 19, 22, 26, 0, 0, 0, 0, 0, 0, 0, 0 } /* 960 framing */
99 };
100
101 /*!
102 \brief Initialize DRC QMF factors
103
104 \hDrcData Handle to DRC channel data.
105
106 \return none
107 */
sbrDecoder_drcInitChannel(HANDLE_SBR_DRC_CHANNEL hDrcData)108 void sbrDecoder_drcInitChannel (
109 HANDLE_SBR_DRC_CHANNEL hDrcData )
110 {
111 int band;
112
113 if (hDrcData == NULL) {
114 return;
115 }
116
117 for (band = 0; band < (64); band++) {
118 hDrcData->prevFact_mag[band] = (FIXP_DBL)MAXVAL_DBL /*FL2FXCONST_DBL(1.0f)*/;
119 }
120
121 for (band = 0; band < SBRDEC_MAX_DRC_BANDS; band++) {
122 hDrcData->currFact_mag[band] = (FIXP_DBL)MAXVAL_DBL /*FL2FXCONST_DBL(1.0f)*/;
123 hDrcData->nextFact_mag[band] = (FIXP_DBL)MAXVAL_DBL /*FL2FXCONST_DBL(1.0f)*/;
124 }
125
126 hDrcData->prevFact_exp = 0;
127 hDrcData->currFact_exp = 0;
128 hDrcData->nextFact_exp = 0;
129
130 hDrcData->numBandsCurr = 0;
131 hDrcData->numBandsNext = 0;
132
133 hDrcData->winSequenceCurr = 0;
134 hDrcData->winSequenceNext = 0;
135
136 hDrcData->drcInterpolationSchemeCurr = 0;
137 hDrcData->drcInterpolationSchemeNext = 0;
138
139 hDrcData->enable = 0;
140 }
141
142
143 /*!
144 \brief Swap DRC QMF scaling factors after they have been applied.
145
146 \hDrcData Handle to DRC channel data.
147
148 \return none
149 */
sbrDecoder_drcUpdateChannel(HANDLE_SBR_DRC_CHANNEL hDrcData)150 void sbrDecoder_drcUpdateChannel (
151 HANDLE_SBR_DRC_CHANNEL hDrcData )
152 {
153 if (hDrcData == NULL) {
154 return;
155 }
156 if (hDrcData->enable != 1) {
157 return;
158 }
159
160 /* swap previous data */
161 FDKmemcpy( hDrcData->currFact_mag,
162 hDrcData->nextFact_mag,
163 SBRDEC_MAX_DRC_BANDS * sizeof(FIXP_DBL) );
164
165 hDrcData->currFact_exp = hDrcData->nextFact_exp;
166
167 hDrcData->numBandsCurr = hDrcData->numBandsNext;
168
169 FDKmemcpy( hDrcData->bandTopCurr,
170 hDrcData->bandTopNext,
171 SBRDEC_MAX_DRC_BANDS * sizeof(USHORT) );
172
173 hDrcData->drcInterpolationSchemeCurr = hDrcData->drcInterpolationSchemeNext;
174
175 hDrcData->winSequenceCurr = hDrcData->winSequenceNext;
176 }
177
178
179 /*!
180 \brief Apply DRC factors slot based.
181
182 \hDrcData Handle to DRC channel data.
183 \qmfRealSlot Pointer to real valued QMF data of one time slot.
184 \qmfImagSlot Pointer to the imaginary QMF data of one time slot.
185 \col Number of the time slot.
186 \numQmfSubSamples Total number of time slots for one frame.
187 \scaleFactor Pointer to the out scale factor of the time slot.
188
189 \return None.
190 */
sbrDecoder_drcApplySlot(HANDLE_SBR_DRC_CHANNEL hDrcData,FIXP_DBL * qmfRealSlot,FIXP_DBL * qmfImagSlot,int col,int numQmfSubSamples,int maxShift)191 void sbrDecoder_drcApplySlot (
192 HANDLE_SBR_DRC_CHANNEL hDrcData,
193 FIXP_DBL *qmfRealSlot,
194 FIXP_DBL *qmfImagSlot,
195 int col,
196 int numQmfSubSamples,
197 int maxShift
198 )
199 {
200 const int *offset;
201
202 int band, bottomMdct, topMdct, bin, useLP;
203 int indx = numQmfSubSamples - (numQmfSubSamples >> 1) - 10; /* l_border */
204 int frameLenFlag = (numQmfSubSamples == 30) ? 1 : 0;
205
206 const FIXP_DBL *fact_mag = NULL;
207 INT fact_exp = 0;
208 UINT numBands = 0;
209 USHORT *bandTop = NULL;
210 int shortDrc = 0;
211
212 FIXP_DBL alphaValue = FL2FXCONST_DBL(0.0f);
213
214 if (hDrcData == NULL) {
215 return;
216 }
217 if (hDrcData->enable != 1) {
218 return;
219 }
220
221 offset = offsetTab[frameLenFlag];
222
223 useLP = (qmfImagSlot == NULL) ? 1 : 0;
224
225 col += indx;
226 bottomMdct = 0;
227 bin = 0;
228
229 /* get respective data and calc interpolation factor */
230 if (col < (numQmfSubSamples>>1)) { /* first half of current frame */
231 if (hDrcData->winSequenceCurr != 2) { /* long window */
232 int j = col + (numQmfSubSamples>>1);
233
234 if (hDrcData->drcInterpolationSchemeCurr == 0) {
235 INT k = (frameLenFlag) ? 0x4444444 : 0x4000000;
236
237 alphaValue = (FIXP_DBL)(j * k);
238 }
239 else {
240 if (j >= offset[hDrcData->drcInterpolationSchemeCurr - 1]) {
241 alphaValue = FL2FXCONST_DBL(1.0f);
242 }
243 }
244 }
245 else { /* short windows */
246 shortDrc = 1;
247 }
248
249 fact_mag = hDrcData->currFact_mag;
250 fact_exp = hDrcData->currFact_exp;
251 numBands = hDrcData->numBandsCurr;
252 bandTop = hDrcData->bandTopCurr;
253 }
254 else if (col < numQmfSubSamples) { /* second half of current frame */
255 if (hDrcData->winSequenceNext != 2) { /* next: long window */
256 int j = col - (numQmfSubSamples>>1);
257
258 if (hDrcData->drcInterpolationSchemeNext == 0) {
259 INT k = (frameLenFlag) ? 0x4444444 : 0x4000000;
260
261 alphaValue = (FIXP_DBL)(j * k);
262 }
263 else {
264 if (j >= offset[hDrcData->drcInterpolationSchemeNext - 1]) {
265 alphaValue = FL2FXCONST_DBL(1.0f);
266 }
267 }
268
269 fact_mag = hDrcData->nextFact_mag;
270 fact_exp = hDrcData->nextFact_exp;
271 numBands = hDrcData->numBandsNext;
272 bandTop = hDrcData->bandTopNext;
273 }
274 else { /* next: short windows */
275 if (hDrcData->winSequenceCurr != 2) { /* current: long window */
276 alphaValue = (FIXP_DBL)0;
277
278 fact_mag = hDrcData->nextFact_mag;
279 fact_exp = hDrcData->nextFact_exp;
280 numBands = hDrcData->numBandsNext;
281 bandTop = hDrcData->bandTopNext;
282 }
283 else { /* current: short windows */
284 shortDrc = 1;
285
286 fact_mag = hDrcData->currFact_mag;
287 fact_exp = hDrcData->currFact_exp;
288 numBands = hDrcData->numBandsCurr;
289 bandTop = hDrcData->bandTopCurr;
290 }
291 }
292 }
293 else { /* first half of next frame */
294 if (hDrcData->winSequenceNext != 2) { /* long window */
295 int j = col - (numQmfSubSamples>>1);
296
297 if (hDrcData->drcInterpolationSchemeNext == 0) {
298 INT k = (frameLenFlag) ? 0x4444444 : 0x4000000;
299
300 alphaValue = (FIXP_DBL)(j * k);
301 }
302 else {
303 if (j >= offset[hDrcData->drcInterpolationSchemeNext - 1]) {
304 alphaValue = FL2FXCONST_DBL(1.0f);
305 }
306 }
307 }
308 else { /* short windows */
309 shortDrc = 1;
310 }
311
312 fact_mag = hDrcData->nextFact_mag;
313 fact_exp = hDrcData->nextFact_exp;
314 numBands = hDrcData->numBandsNext;
315 bandTop = hDrcData->bandTopNext;
316
317 col -= numQmfSubSamples;
318 }
319
320
321 /* process bands */
322 for (band = 0; band < (int)numBands; band++) {
323 int bottomQmf, topQmf;
324
325 FIXP_DBL drcFact_mag = FL2FXCONST_DBL(1.0f);
326
327 topMdct = (bandTop[band]+1) << 2;
328
329 if (!shortDrc) { /* long window */
330 if (frameLenFlag) {
331 /* 960 framing */
332 bottomMdct = 30 * (bottomMdct / 30);
333 topMdct = 30 * (topMdct / 30);
334
335 bottomQmf = fMultIfloor((FIXP_DBL)0x4444444, bottomMdct);
336 topQmf = fMultIfloor((FIXP_DBL)0x4444444, topMdct);
337 }
338 else {
339 /* 1024 framing */
340 bottomMdct &= ~0x1f;
341 topMdct &= ~0x1f;
342
343 bottomQmf = bottomMdct >> 5;
344 topQmf = topMdct >> 5;
345 }
346
347 if (band == ((int)numBands-1)) {
348 topQmf = (64);
349 }
350
351 for (bin = bottomQmf; bin < topQmf; bin++) {
352 FIXP_DBL drcFact1_mag = hDrcData->prevFact_mag[bin];
353 FIXP_DBL drcFact2_mag = fact_mag[band];
354
355 /* normalize scale factors */
356 if (hDrcData->prevFact_exp < maxShift) {
357 drcFact1_mag >>= maxShift - hDrcData->prevFact_exp;
358 }
359 if (fact_exp < maxShift) {
360 drcFact2_mag >>= maxShift - fact_exp;
361 }
362
363 /* interpolate */
364 drcFact_mag = fMult(alphaValue, drcFact2_mag) + fMult((FL2FXCONST_DBL(1.0f) - alphaValue), drcFact1_mag);
365
366 /* apply scaling */
367 qmfRealSlot[bin] = fMult(qmfRealSlot[bin], drcFact_mag);
368 if (!useLP) {
369 qmfImagSlot[bin] = fMult(qmfImagSlot[bin], drcFact_mag);
370 }
371
372 /* save previous factors */
373 if (col == (numQmfSubSamples>>1)-1) {
374 hDrcData->prevFact_mag[bin] = fact_mag[band];
375 }
376 }
377 }
378 else { /* short windows */
379 int startSample, stopSample;
380 FIXP_DBL invFrameSizeDiv8 = (frameLenFlag) ? (FIXP_DBL)0x1111111 : (FIXP_DBL)0x1000000;
381
382 if (frameLenFlag) {
383 /* 960 framing */
384 bottomMdct = 30/8 * (bottomMdct*8/30);
385 topMdct = 30/8 * (topMdct*8/30);
386 }
387 else {
388 /* 1024 framing */
389 bottomMdct &= ~0x03;
390 topMdct &= ~0x03;
391 }
392
393 /* startSample is truncated to the nearest corresponding start subsample in
394 the QMF of the short window bottom is present in:*/
395 startSample = ((fMultIfloor( invFrameSizeDiv8, bottomMdct ) & 0x7) * numQmfSubSamples) >> 3;
396
397 /* stopSample is rounded upwards to the nearest corresponding stop subsample
398 in the QMF of the short window top is present in. */
399 stopSample = ((fMultIceil( invFrameSizeDiv8, topMdct ) & 0xf) * numQmfSubSamples) >> 3;
400
401 bottomQmf = fMultIfloor( invFrameSizeDiv8, ((bottomMdct%(numQmfSubSamples<<2)) << 5) );
402 topQmf = fMultIfloor( invFrameSizeDiv8, ((topMdct%(numQmfSubSamples<<2)) << 5) );
403
404 /* extend last band */
405 if (band == ((int)numBands-1)) {
406 topQmf = (64);
407 stopSample = numQmfSubSamples;
408 }
409
410 if (topQmf == 0) {
411 topQmf = (64);
412 }
413
414 /* save previous factors */
415 if (stopSample == numQmfSubSamples) {
416 int tmpBottom = bottomQmf;
417
418 if (((numQmfSubSamples-1) & ~0x03) > startSample) {
419 tmpBottom = 0; /* band starts in previous short window */
420 }
421
422 for (bin = tmpBottom; bin < topQmf; bin++) {
423 hDrcData->prevFact_mag[bin] = fact_mag[band];
424 }
425 }
426
427 /* apply */
428 if ((col >= startSample) && (col < stopSample)) {
429 if ((col & ~0x03) > startSample) {
430 bottomQmf = 0; /* band starts in previous short window */
431 }
432 if (col < ((stopSample-1) & ~0x03)) {
433 topQmf = (64); /* band ends in next short window */
434 }
435
436 drcFact_mag = fact_mag[band];
437
438 /* normalize scale factor */
439 if (fact_exp < maxShift) {
440 drcFact_mag >>= maxShift - fact_exp;
441 }
442
443 /* apply scaling */
444 for (bin = bottomQmf; bin < topQmf; bin++) {
445 qmfRealSlot[bin] = fMult(qmfRealSlot[bin], drcFact_mag);
446 if (!useLP) {
447 qmfImagSlot[bin] = fMult(qmfImagSlot[bin], drcFact_mag);
448 }
449 }
450 }
451 }
452
453 bottomMdct = topMdct;
454 } /* end of bands loop */
455
456 if (col == (numQmfSubSamples>>1)-1) {
457 hDrcData->prevFact_exp = fact_exp;
458 }
459 }
460
461
462 /*!
463 \brief Apply DRC factors frame based.
464
465 \hDrcData Handle to DRC channel data.
466 \qmfRealSlot Pointer to real valued QMF data of the whole frame.
467 \qmfImagSlot Pointer to the imaginary QMF data of the whole frame.
468 \numQmfSubSamples Total number of time slots for one frame.
469 \scaleFactor Pointer to the out scale factor of the frame.
470
471 \return None.
472 */
sbrDecoder_drcApply(HANDLE_SBR_DRC_CHANNEL hDrcData,FIXP_DBL ** QmfBufferReal,FIXP_DBL ** QmfBufferImag,int numQmfSubSamples,int * scaleFactor)473 void sbrDecoder_drcApply (
474 HANDLE_SBR_DRC_CHANNEL hDrcData,
475 FIXP_DBL **QmfBufferReal,
476 FIXP_DBL **QmfBufferImag,
477 int numQmfSubSamples,
478 int *scaleFactor
479 )
480 {
481 int col;
482 int maxShift = 0;
483
484 /* get max scale factor */
485 if (hDrcData->prevFact_exp > maxShift) {
486 maxShift = hDrcData->prevFact_exp;
487 }
488 if (hDrcData->currFact_exp > maxShift) {
489 maxShift = hDrcData->currFact_exp;
490 }
491 if (hDrcData->nextFact_exp > maxShift) {
492 maxShift = hDrcData->nextFact_exp;
493 }
494
495 for (col = 0; col < numQmfSubSamples; col++)
496 {
497 FIXP_DBL *qmfSlotReal = QmfBufferReal[col];
498 FIXP_DBL *qmfSlotImag = (QmfBufferImag == NULL) ? NULL : QmfBufferImag[col];
499
500 sbrDecoder_drcApplySlot (
501 hDrcData,
502 qmfSlotReal,
503 qmfSlotImag,
504 col,
505 numQmfSubSamples,
506 maxShift
507 );
508 }
509
510 *scaleFactor += maxShift;
511 }
512
513