Home | History | Annotate | Download | only in src
      1 
      2 /* -----------------------------------------------------------------------------------------------------------
      3 Software License for The Fraunhofer FDK AAC Codec Library for Android
      4 
      5  Copyright  1995 - 2013 Fraunhofer-Gesellschaft zur Frderung 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 (at) iis.fraunhofer.de
     82 ----------------------------------------------------------------------------------------------------------- */
     83 
     84 /*****************************  MPEG-4 AAC Decoder  **************************
     85 
     86    Author(s):   Josef Hoepfl
     87    Description: long/short-block decoding
     88 
     89 ******************************************************************************/
     90 
     91 #include "block.h"
     92 
     93 #include "aac_rom.h"
     94 #include "FDK_bitstream.h"
     95 #include "FDK_tools_rom.h"
     96 
     97 
     98 
     99 
    100 #include "aacdec_hcr.h"
    101 #include "rvlc.h"
    102 
    103 
    104 #if defined(__arm__)
    105 #include "arm/block_arm.cpp"
    106 #endif
    107 
    108 /*!
    109   \brief Read escape sequence of codeword
    110 
    111   The function reads the escape sequence from the bitstream,
    112   if the absolute value of the quantized coefficient has the
    113   value 16.
    114 
    115   \return  quantized coefficient
    116 */
    117 LONG CBlock_GetEscape(HANDLE_FDK_BITSTREAM bs, /*!< pointer to bitstream */
    118                      const LONG q)        /*!< quantized coefficient */
    119 {
    120   LONG i, off, neg ;
    121 
    122   if (q < 0)
    123   {
    124     if (q != -16) return q;
    125     neg = 1;
    126   }
    127   else
    128   {
    129     if (q != +16) return q;
    130     neg = 0;
    131   }
    132 
    133   for (i=4; ; i++)
    134   {
    135     if (FDKreadBits(bs,1) == 0)
    136       break;
    137   }
    138 
    139   if (i > 16)
    140   {
    141     if (i - 16 > CACHE_BITS) { /* cannot read more than "CACHE_BITS" bits at once in the function FDKreadBits() */
    142       return (MAX_QUANTIZED_VALUE + 1); /* returning invalid value that will be captured later */
    143     }
    144 
    145     off = FDKreadBits(bs,i-16) << 16;
    146     off |= FDKreadBits(bs,16);
    147   }
    148   else
    149   {
    150     off = FDKreadBits(bs,i);
    151   }
    152 
    153   i = off + (1 << i);
    154 
    155   if (neg) i = -i;
    156 
    157   return i;
    158 }
    159 
    160 AAC_DECODER_ERROR CBlock_ReadScaleFactorData(
    161         CAacDecoderChannelInfo *pAacDecoderChannelInfo,
    162         HANDLE_FDK_BITSTREAM bs,
    163         UINT flags
    164         )
    165 {
    166   int temp;
    167   int band;
    168   int group;
    169   int position = 0; /* accu for intensity delta coding */
    170   int factor = pAacDecoderChannelInfo->pDynData->RawDataInfo.GlobalGain; /* accu for scale factor delta coding */
    171   UCHAR *pCodeBook = pAacDecoderChannelInfo->pDynData->aCodeBook;
    172   SHORT *pScaleFactor = pAacDecoderChannelInfo->pDynData->aScaleFactor;
    173   const CodeBookDescription *hcb =&AACcodeBookDescriptionTable[BOOKSCL];
    174 
    175   int ScaleFactorBandsTransmitted = GetScaleFactorBandsTransmitted(&pAacDecoderChannelInfo->icsInfo);
    176   for (group=0; group < GetWindowGroups(&pAacDecoderChannelInfo->icsInfo); group++)
    177   {
    178     for (band=0; band < ScaleFactorBandsTransmitted; band++)
    179     {
    180       switch (pCodeBook[group*16+band]) {
    181 
    182       case ZERO_HCB: /* zero book */
    183         pScaleFactor[group*16+band] = 0;
    184         break;
    185 
    186       default: /* decode scale factor */
    187         {
    188           temp = CBlock_DecodeHuffmanWord(bs,hcb);
    189           factor += temp - 60; /* MIDFAC 1.5 dB */
    190         }
    191         pScaleFactor[group*16+band] = factor - 100;
    192         break;
    193 
    194       case INTENSITY_HCB: /* intensity steering */
    195       case INTENSITY_HCB2:
    196         temp = CBlock_DecodeHuffmanWord(bs,hcb);
    197         position += temp - 60;
    198         pScaleFactor[group*16+band] = position - 100;
    199         break;
    200 
    201       case NOISE_HCB: /* PNS */
    202         if (flags & (AC_MPS_RES|AC_USAC|AC_RSVD50)) {
    203           return AAC_DEC_PARSE_ERROR;
    204         }
    205         CPns_Read( &pAacDecoderChannelInfo->data.aac.PnsData, bs, hcb, pAacDecoderChannelInfo->pDynData->aScaleFactor, pAacDecoderChannelInfo->pDynData->RawDataInfo.GlobalGain, band, group);
    206         break;
    207       }
    208     }
    209   }
    210 
    211   return AAC_DEC_OK;
    212 }
    213 
    214 void CBlock_ScaleSpectralData(CAacDecoderChannelInfo *pAacDecoderChannelInfo, SamplingRateInfo *pSamplingRateInfo)
    215 {
    216   int band;
    217   int window;
    218   const SHORT * RESTRICT pSfbScale  = pAacDecoderChannelInfo->pDynData->aSfbScale;
    219   SHORT * RESTRICT pSpecScale = pAacDecoderChannelInfo->specScale;
    220   int groupwin,group;
    221   const SHORT * RESTRICT BandOffsets = GetScaleFactorBandOffsets(&pAacDecoderChannelInfo->icsInfo, pSamplingRateInfo);
    222   SPECTRAL_PTR RESTRICT pSpectralCoefficient = pAacDecoderChannelInfo->pSpectralCoefficient;
    223 
    224 
    225   FDKmemclear(pSpecScale, 8*sizeof(SHORT));
    226 
    227   int max_band = GetScaleFactorBandsTransmitted(&pAacDecoderChannelInfo->icsInfo);
    228   for (window=0, group=0; group < GetWindowGroups(&pAacDecoderChannelInfo->icsInfo); group++)
    229   {
    230     for (groupwin=0; groupwin < GetWindowGroupLength(&pAacDecoderChannelInfo->icsInfo,group); groupwin++, window++)
    231     {
    232       int SpecScale_window = pSpecScale[window];
    233       FIXP_DBL *pSpectrum = SPEC(pSpectralCoefficient, window,  pAacDecoderChannelInfo->granuleLength);
    234 
    235       /* find scaling for current window */
    236       for (band=0; band < max_band; band++)
    237       {
    238         SpecScale_window = fMax(SpecScale_window, (int)pSfbScale[window*16+band]);
    239       }
    240 
    241       if (pAacDecoderChannelInfo->pDynData->TnsData.Active) {
    242         SpecScale_window += TNS_SCALE;
    243       }
    244 
    245       /* store scaling of current window */
    246       pSpecScale[window] = SpecScale_window;
    247 
    248 #ifdef FUNCTION_CBlock_ScaleSpectralData_func1
    249 
    250       CBlock_ScaleSpectralData_func1(pSpectrum, max_band, BandOffsets, SpecScale_window, pSfbScale, window);
    251 
    252 #else /* FUNCTION_CBlock_ScaleSpectralData_func1 */
    253       for (band=0; band < max_band; band++)
    254       {
    255         int scale = SpecScale_window - pSfbScale[window*16+band];
    256         if (scale)
    257         {
    258           /* following relation can be used for optimizations: (BandOffsets[i]%4) == 0 for all i */
    259           int max_index = BandOffsets[band+1];
    260           for (int index = BandOffsets[band]; index < max_index; index++)
    261           {
    262             pSpectrum[index] >>= scale;
    263           }
    264         }
    265       }
    266 #endif  /* FUNCTION_CBlock_ScaleSpectralData_func1 */
    267     }
    268   }
    269 
    270 }
    271 
    272 AAC_DECODER_ERROR CBlock_ReadSectionData(HANDLE_FDK_BITSTREAM bs,
    273                                          CAacDecoderChannelInfo *pAacDecoderChannelInfo,
    274                                          const SamplingRateInfo *pSamplingRateInfo,
    275                                          const UINT  flags)
    276 {
    277   int top, band;
    278   int sect_len, sect_len_incr;
    279   int group;
    280   UCHAR sect_cb;
    281   UCHAR *pCodeBook = pAacDecoderChannelInfo->pDynData->aCodeBook;
    282   /* HCR input (long) */
    283   SHORT *pNumLinesInSec    = pAacDecoderChannelInfo->pDynData->specificTo.aac.aNumLineInSec4Hcr;
    284   int    numLinesInSecIdx  = 0;
    285   UCHAR *pHcrCodeBook      = pAacDecoderChannelInfo->pDynData->specificTo.aac.aCodeBooks4Hcr;
    286   const SHORT *BandOffsets = GetScaleFactorBandOffsets(&pAacDecoderChannelInfo->icsInfo, pSamplingRateInfo);
    287   pAacDecoderChannelInfo->pDynData->specificTo.aac.numberSection = 0;
    288   AAC_DECODER_ERROR ErrorStatus = AAC_DEC_OK;
    289 
    290   FDKmemclear(pCodeBook, sizeof(UCHAR)*(8*16));
    291 
    292   const int nbits = (IsLongBlock(&pAacDecoderChannelInfo->icsInfo) == 1) ? 5 : 3;
    293 
    294   int sect_esc_val = (1 << nbits) - 1 ;
    295 
    296   UCHAR ScaleFactorBandsTransmitted = GetScaleFactorBandsTransmitted(&pAacDecoderChannelInfo->icsInfo);
    297   for (group=0; group<GetWindowGroups(&pAacDecoderChannelInfo->icsInfo); group++)
    298   {
    299     for (band=0; band < ScaleFactorBandsTransmitted; )
    300     {
    301       sect_len = 0;
    302       if ( flags & AC_ER_VCB11 )  {
    303         sect_cb = (UCHAR) FDKreadBits(bs,5);
    304       }
    305       else
    306         sect_cb = (UCHAR) FDKreadBits(bs,4);
    307 
    308       if ( ((flags & AC_ER_VCB11) == 0) || ( sect_cb < 11 ) || ((sect_cb > 11) && (sect_cb < 16)) ) {
    309         sect_len_incr = FDKreadBits(bs, nbits);
    310         while (sect_len_incr == sect_esc_val)
    311         {
    312           sect_len += sect_esc_val;
    313           sect_len_incr = FDKreadBits(bs, nbits);
    314         }
    315       }
    316       else {
    317         sect_len_incr = 1;
    318       }
    319 
    320       sect_len += sect_len_incr;
    321 
    322 
    323       top = band + sect_len;
    324 
    325       if (flags & AC_ER_HCR) {
    326         /* HCR input (long) -- collecting sideinfo (for HCR-_long_ only) */
    327         pNumLinesInSec[numLinesInSecIdx] = BandOffsets[top] - BandOffsets[band];
    328         numLinesInSecIdx++;
    329         if (numLinesInSecIdx >= MAX_SFB_HCR) {
    330           return AAC_DEC_PARSE_ERROR;
    331         }
    332         if (
    333              (sect_cb == BOOKSCL) )
    334         {
    335           return AAC_DEC_INVALID_CODE_BOOK;
    336         } else {
    337           *pHcrCodeBook++ = sect_cb;
    338         }
    339         pAacDecoderChannelInfo->pDynData->specificTo.aac.numberSection++;
    340       }
    341 
    342       /* Check spectral line limits */
    343       if (IsLongBlock( &(pAacDecoderChannelInfo->icsInfo) ))
    344       {
    345         if (top > 64) {
    346           return AAC_DEC_DECODE_FRAME_ERROR;
    347         }
    348       } else { /* short block */
    349         if (top + group*16 > (8 * 16)) {
    350           return AAC_DEC_DECODE_FRAME_ERROR;
    351         }
    352       }
    353 
    354       /* Check if decoded codebook index is feasible */
    355       if ( (sect_cb == BOOKSCL)
    356        || ( (sect_cb == INTENSITY_HCB || sect_cb == INTENSITY_HCB2) && pAacDecoderChannelInfo->pDynData->RawDataInfo.CommonWindow == 0)
    357          )
    358       {
    359         return AAC_DEC_INVALID_CODE_BOOK;
    360       }
    361 
    362       /* Store codebook index */
    363       for (; band < top; band++)
    364       {
    365         pCodeBook[group*16+band] = sect_cb;
    366       }
    367     }
    368   }
    369 
    370 
    371   return ErrorStatus;
    372 }
    373 
    374 /* mso: provides a faster way to i-quantize a whole band in one go */
    375 
    376 /**
    377  * \brief inverse quantize one sfb. Each value of the sfb is processed according to the
    378  *        formula: spectrum[i] = Sign(spectrum[i]) * Matissa(spectrum[i])^(4/3) * 2^(lsb/4).
    379  * \param spectrum pointer to first line of the sfb to be inverse quantized.
    380  * \param noLines number of lines belonging to the sfb.
    381  * \param lsb last 2 bits of the scale factor of the sfb.
    382  * \param scale max allowed shift scale for the sfb.
    383  */
    384 static
    385 void InverseQuantizeBand( FIXP_DBL * RESTRICT spectrum,
    386                               INT noLines,
    387                               INT lsb,
    388                               INT scale )
    389 {
    390     const FIXP_DBL * RESTRICT InverseQuantTabler=(FIXP_DBL *)InverseQuantTable;
    391     const FIXP_DBL * RESTRICT MantissaTabler=(FIXP_DBL *)MantissaTable[lsb];
    392     const SCHAR* RESTRICT ExponentTabler=(SCHAR*)ExponentTable[lsb];
    393 
    394     FIXP_DBL *ptr = spectrum;
    395     FIXP_DBL signedValue;
    396 
    397     FDK_ASSERT(noLines>2);
    398     for (INT i=noLines; i--; )
    399     {
    400         if ((signedValue = *ptr++) != FL2FXCONST_DBL(0))
    401         {
    402           FIXP_DBL value = fAbs(signedValue);
    403           UINT freeBits = CntLeadingZeros(value);
    404           UINT exponent = 32 - freeBits;
    405 
    406           UINT x = (UINT) (LONG)value << (INT) freeBits;
    407           x <<= 1;                                  /* shift out sign bit to avoid masking later on */
    408           UINT tableIndex = x >> 24;
    409           x = (x >> 20) &  0x0F;
    410 
    411           UINT r0=(UINT)(LONG)InverseQuantTabler[tableIndex+0];
    412           UINT r1=(UINT)(LONG)InverseQuantTabler[tableIndex+1];
    413           UINT temp= (r1 - r0)*x + (r0 << 4);
    414 
    415           value = fMultDiv2((FIXP_DBL)temp, MantissaTabler[exponent]);
    416 
    417           /* + 1 compensates fMultDiv2() */
    418           scaleValueInPlace(&value, scale + ExponentTabler[exponent] + 1);
    419 
    420           signedValue = (signedValue < (FIXP_DBL)0) ? -value : value;
    421           ptr[-1] = signedValue;
    422         }
    423     }
    424 }
    425 
    426 AAC_DECODER_ERROR CBlock_InverseQuantizeSpectralData(CAacDecoderChannelInfo *pAacDecoderChannelInfo, SamplingRateInfo *pSamplingRateInfo)
    427 {
    428   int window, group, groupwin, band;
    429   int ScaleFactorBandsTransmitted = GetScaleFactorBandsTransmitted(&pAacDecoderChannelInfo->icsInfo);
    430   UCHAR *RESTRICT pCodeBook = pAacDecoderChannelInfo->pDynData->aCodeBook;
    431   SHORT *RESTRICT pSfbScale = pAacDecoderChannelInfo->pDynData->aSfbScale;
    432   SHORT *RESTRICT pScaleFactor = pAacDecoderChannelInfo->pDynData->aScaleFactor;
    433   const SHORT *RESTRICT BandOffsets = GetScaleFactorBandOffsets(&pAacDecoderChannelInfo->icsInfo, pSamplingRateInfo);
    434 
    435   FDKmemclear(pAacDecoderChannelInfo->pDynData->aSfbScale, (8*16)*sizeof(SHORT));
    436 
    437   for (window=0, group=0; group < GetWindowGroups(&pAacDecoderChannelInfo->icsInfo); group++)
    438   {
    439     for (groupwin=0; groupwin < GetWindowGroupLength(&pAacDecoderChannelInfo->icsInfo,group); groupwin++, window++)
    440     {
    441       /* inverse quantization */
    442       for (band=0; band < ScaleFactorBandsTransmitted; band++)
    443       {
    444         FIXP_DBL *pSpectralCoefficient = SPEC(pAacDecoderChannelInfo->pSpectralCoefficient, window, pAacDecoderChannelInfo->granuleLength) + BandOffsets[band];
    445 
    446         int noLines = BandOffsets[band+1] - BandOffsets[band];
    447         int bnds = group*16+band;
    448         int i;
    449 
    450         if ((pCodeBook[bnds] == ZERO_HCB)
    451          || (pCodeBook[bnds] == INTENSITY_HCB)
    452          || (pCodeBook[bnds] == INTENSITY_HCB2)
    453            )
    454           continue;
    455 
    456         if (pCodeBook[bnds] == NOISE_HCB)
    457         {
    458           /* Leave headroom for PNS values. + 1 because ceil(log2(2^(0.25*3))) = 1,
    459              worst case of additional headroom required because of the scalefactor. */
    460           pSfbScale[window*16+band] = (pScaleFactor [bnds] >> 2) + 1 ;
    461           continue;
    462         }
    463 
    464         /* Find max spectral line value of the current sfb */
    465         FIXP_DBL locMax = (FIXP_DBL)0;
    466 
    467         for (i = noLines; i-- ; ) {
    468           /* Expensive memory access */
    469           locMax = fMax(fixp_abs(pSpectralCoefficient[i]), locMax);
    470         }
    471 
    472         /* Cheap robustness improvement - Do not remove!!! */
    473         if (fixp_abs(locMax) > (FIXP_DBL)MAX_QUANTIZED_VALUE) {
    474           return AAC_DEC_DECODE_FRAME_ERROR;
    475         }
    476 
    477         /*
    478            The inverse quantized spectral lines are defined by:
    479         pSpectralCoefficient[i] = Sign(pSpectralCoefficient[i]) * 2^(0.25*pScaleFactor[bnds]) * pSpectralCoefficient[i]^(4/3)
    480            This is equivalent to:
    481         pSpectralCoefficient[i]    = Sign(pSpectralCoefficient[i]) * (2^(pScaleFactor[bnds] % 4) * pSpectralCoefficient[i]^(4/3))
    482         pSpectralCoefficient_e[i] += pScaleFactor[bnds]/4
    483         */
    484         {
    485           int msb = pScaleFactor [bnds] >> 2 ;
    486           int lsb = pScaleFactor [bnds] & 0x03 ;
    487 
    488           int scale = GetScaleFromValue(locMax, lsb);
    489 
    490           pSfbScale[window*16+band] = msb - scale;
    491           InverseQuantizeBand(pSpectralCoefficient, noLines, lsb, scale);
    492         }
    493       }
    494     }
    495   }
    496 
    497 
    498   return AAC_DEC_OK;
    499 }
    500 
    501 
    502 AAC_DECODER_ERROR  CBlock_ReadSpectralData(HANDLE_FDK_BITSTREAM bs,
    503                                            CAacDecoderChannelInfo *pAacDecoderChannelInfo,
    504                                            const SamplingRateInfo *pSamplingRateInfo,
    505                                            const UINT  flags)
    506 {
    507   int i,index;
    508   int window,group,groupwin,groupoffset,band;
    509   UCHAR *RESTRICT pCodeBook = pAacDecoderChannelInfo->pDynData->aCodeBook;
    510   const SHORT *RESTRICT BandOffsets = GetScaleFactorBandOffsets(&pAacDecoderChannelInfo->icsInfo, pSamplingRateInfo);
    511 
    512   SPECTRAL_PTR pSpectralCoefficient = pAacDecoderChannelInfo->pSpectralCoefficient;
    513   FIXP_DBL locMax;
    514 
    515   int ScaleFactorBandsTransmitted = GetScaleFactorBandsTransmitted(&pAacDecoderChannelInfo->icsInfo);
    516 
    517   FDK_ASSERT(BandOffsets != NULL);
    518 
    519   FDKmemclear(pSpectralCoefficient, sizeof(SPECTRUM));
    520 
    521   if ( (flags & AC_ER_HCR) == 0 )
    522   {
    523     groupoffset = 0;
    524 
    525     /* plain huffman decoder  short */
    526     for (group=0; group < GetWindowGroups(&pAacDecoderChannelInfo->icsInfo); group++)
    527     {
    528       for (band=0; band < ScaleFactorBandsTransmitted; band++)
    529       {
    530         int bnds = group*16+band;
    531         UCHAR currentCB = pCodeBook[bnds];
    532 
    533         /* patch to run plain-huffman-decoder with vcb11 input codebooks (LAV-checking might be possible below using the virtual cb and a LAV-table) */
    534         if ((currentCB >= 16) && (currentCB <= 31)) {
    535           pCodeBook[bnds] = currentCB = 11;
    536         }
    537         if ( !((currentCB == ZERO_HCB)
    538             || (currentCB == NOISE_HCB)
    539             || (currentCB == INTENSITY_HCB)
    540             || (currentCB == INTENSITY_HCB2)) )
    541         {
    542           const CodeBookDescription *hcb = &AACcodeBookDescriptionTable[currentCB];
    543           int step = hcb->Dimension;
    544           int offset = hcb->Offset;
    545           int bits = hcb->numBits;
    546           int mask = (1<<bits)-1;
    547 
    548           for (groupwin=0; groupwin < GetWindowGroupLength(&pAacDecoderChannelInfo->icsInfo,group); groupwin++)
    549           {
    550             window = groupoffset + groupwin;
    551 
    552             FIXP_DBL *mdctSpectrum = SPEC(pSpectralCoefficient, window, pAacDecoderChannelInfo->granuleLength);
    553 
    554             locMax = (FIXP_DBL)0 ;
    555 
    556             for (index=BandOffsets[band]; index < BandOffsets[band+1]; index+=step)
    557             {
    558               int idx = CBlock_DecodeHuffmanWord(bs,hcb);
    559 
    560               for (i=0; i<step; i++) {
    561                 FIXP_DBL tmp;
    562 
    563                 tmp = (FIXP_DBL)((idx & mask)-offset);
    564                 idx >>= bits;
    565 
    566                 if (offset == 0) {
    567                   if (tmp != FIXP_DBL(0))
    568                     tmp = (FDKreadBits(bs,1))? -tmp : tmp;
    569                 }
    570                 mdctSpectrum[index+i] = tmp;
    571               }
    572 
    573               if (currentCB == ESCBOOK)
    574               {
    575                 mdctSpectrum[index+0] = (FIXP_DBL)CBlock_GetEscape(bs, (LONG)mdctSpectrum[index+0]);
    576                 mdctSpectrum[index+1] = (FIXP_DBL)CBlock_GetEscape(bs, (LONG)mdctSpectrum[index+1]);
    577 
    578               }
    579             }
    580           }
    581         }
    582       }
    583       groupoffset += GetWindowGroupLength(&pAacDecoderChannelInfo->icsInfo,group);
    584     }
    585     /* plain huffman decoding (short) finished */
    586   }
    587   /* HCR - Huffman Codeword Reordering  short */
    588   else  /* if ( flags & AC_ER_HCR ) */
    589   {
    590     H_HCR_INFO hHcr = &pAacDecoderChannelInfo->pComData->overlay.aac.erHcrInfo;
    591     int hcrStatus = 0;
    592 
    593     /* advanced Huffman decoding starts here (HCR decoding :) */
    594     if ( pAacDecoderChannelInfo->pDynData->specificTo.aac.lenOfReorderedSpectralData != 0 ) {
    595 
    596       /* HCR initialization short */
    597       hcrStatus = HcrInit(hHcr, pAacDecoderChannelInfo, pSamplingRateInfo, bs);
    598 
    599       if (hcrStatus != 0) {
    600         return AAC_DEC_DECODE_FRAME_ERROR;
    601       }
    602 
    603       /* HCR decoding short */
    604       hcrStatus = HcrDecoder(hHcr, pAacDecoderChannelInfo, pSamplingRateInfo, bs);
    605 
    606       if (hcrStatus != 0) {
    607 #if HCR_ERROR_CONCEALMENT
    608         HcrMuteErroneousLines(hHcr);
    609 #else
    610         return AAC_DEC_DECODE_FRAME_ERROR;
    611 #endif /* HCR_ERROR_CONCEALMENT */
    612       }
    613 
    614       FDKpushFor (bs, pAacDecoderChannelInfo->pDynData->specificTo.aac.lenOfReorderedSpectralData);
    615     }
    616   }
    617   /* HCR - Huffman Codeword Reordering short finished */
    618 
    619 
    620 
    621   if ( IsLongBlock(&pAacDecoderChannelInfo->icsInfo) && !(flags & (AC_ELD|AC_SCALABLE)) )
    622   {
    623     /* apply pulse data */
    624     CPulseData_Apply(&pAacDecoderChannelInfo->pDynData->specificTo.aac.PulseData,
    625                       GetScaleFactorBandOffsets(&pAacDecoderChannelInfo->icsInfo, pSamplingRateInfo),
    626                       SPEC_LONG(pSpectralCoefficient));
    627   }
    628 
    629 
    630   return AAC_DEC_OK;
    631 }
    632 
    633 
    634 
    635 void ApplyTools ( CAacDecoderChannelInfo *pAacDecoderChannelInfo[],
    636                   const SamplingRateInfo *pSamplingRateInfo,
    637                   const UINT flags,
    638                   const int channel )
    639 {
    640 
    641   if ( !(flags & (AC_USAC|AC_RSVD50|AC_MPS_RES)) ) {
    642     CPns_Apply(
    643            &pAacDecoderChannelInfo[channel]->data.aac.PnsData,
    644            &pAacDecoderChannelInfo[channel]->icsInfo,
    645             pAacDecoderChannelInfo[channel]->pSpectralCoefficient,
    646             pAacDecoderChannelInfo[channel]->specScale,
    647             pAacDecoderChannelInfo[channel]->pDynData->aScaleFactor,
    648             pSamplingRateInfo,
    649             pAacDecoderChannelInfo[channel]->granuleLength,
    650             channel
    651             );
    652   }
    653 
    654   CTns_Apply (
    655          &pAacDecoderChannelInfo[channel]->pDynData->TnsData,
    656          &pAacDecoderChannelInfo[channel]->icsInfo,
    657           pAacDecoderChannelInfo[channel]->pSpectralCoefficient,
    658           pSamplingRateInfo,
    659           pAacDecoderChannelInfo[channel]->granuleLength
    660           );
    661 }
    662 
    663 static
    664 int getWindow2Nr(int length, int shape)
    665 {
    666   int nr = 0;
    667 
    668   if (shape == 2) {
    669     /* Low Overlap, 3/4 zeroed */
    670     nr = (length * 3)>>2;
    671   }
    672 
    673   return nr;
    674 }
    675 
    676 void CBlock_FrequencyToTime(CAacDecoderStaticChannelInfo *pAacDecoderStaticChannelInfo,
    677                             CAacDecoderChannelInfo *pAacDecoderChannelInfo,
    678                             INT_PCM outSamples[],
    679                             const SHORT frameLen,
    680                             const int stride,
    681                             const int frameOk,
    682                             FIXP_DBL *pWorkBuffer1 )
    683 {
    684   int fr, fl, tl, nSamples, nSpec;
    685 
    686   /* Determine left slope length (fl), right slope length (fr) and transform length (tl).
    687      USAC: The slope length may mismatch with the previous frame in case of LPD / FD
    688            transitions. The adjustment is handled by the imdct implementation.
    689   */
    690   tl = frameLen;
    691   nSpec = 1;
    692 
    693   switch( pAacDecoderChannelInfo->icsInfo.WindowSequence ) {
    694     default:
    695     case OnlyLongSequence:
    696       fl = frameLen;
    697       fr = frameLen - getWindow2Nr(frameLen, GetWindowShape(&pAacDecoderChannelInfo->icsInfo));
    698       break;
    699     case LongStopSequence:
    700       fl = frameLen >> 3;
    701       fr = frameLen;
    702       break;
    703     case LongStartSequence: /* or StopStartSequence */
    704       fl = frameLen;
    705       fr = frameLen >> 3;
    706       break;
    707     case EightShortSequence:
    708       fl = fr = frameLen >> 3;
    709       tl >>= 3;
    710       nSpec = 8;
    711       break;
    712   }
    713 
    714   {
    715     int i;
    716 
    717     {
    718       FIXP_DBL *tmp = pAacDecoderChannelInfo->pComData->workBufferCore1->mdctOutTemp;
    719 
    720       nSamples = imdct_block(
    721              &pAacDecoderStaticChannelInfo->IMdct,
    722               tmp,
    723               SPEC_LONG(pAacDecoderChannelInfo->pSpectralCoefficient),
    724               pAacDecoderChannelInfo->specScale,
    725               nSpec,
    726               frameLen,
    727               tl,
    728               FDKgetWindowSlope(fl, GetWindowShape(&pAacDecoderChannelInfo->icsInfo)),
    729               fl,
    730               FDKgetWindowSlope(fr, GetWindowShape(&pAacDecoderChannelInfo->icsInfo)),
    731               fr,
    732               (FIXP_DBL)0 );
    733 
    734       for (i=0; i<frameLen; i++) {
    735         outSamples[i*stride] = IMDCT_SCALE(tmp[i]);
    736       }
    737     }
    738   }
    739 
    740   FDK_ASSERT(nSamples == frameLen);
    741 
    742 }
    743 
    744 #include "ldfiltbank.h"
    745 void CBlock_FrequencyToTimeLowDelay( CAacDecoderStaticChannelInfo *pAacDecoderStaticChannelInfo,
    746                                      CAacDecoderChannelInfo *pAacDecoderChannelInfo,
    747                                      INT_PCM outSamples[],
    748                                      const short frameLen,
    749                                      const char stride )
    750 {
    751   InvMdctTransformLowDelay_fdk (
    752           SPEC_LONG(pAacDecoderChannelInfo->pSpectralCoefficient),
    753           pAacDecoderChannelInfo->specScale[0],
    754           outSamples,
    755           pAacDecoderStaticChannelInfo->pOverlapBuffer,
    756           stride,
    757           frameLen
    758           );
    759 }
    760