12228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
22228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/* -----------------------------------------------------------------------------------------------------------
32228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectSoftware License for The Fraunhofer FDK AAC Codec Library for Android
42228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
54f0d97057c5c640b25518358886f8c47da9fc052Jean-Michel Trivi� Copyright  1995 - 2013 Fraunhofer-Gesellschaft zur F�rderung der angewandten Forschung e.V.
62228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  All rights reserved.
72228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
82228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1.    INTRODUCTION
92228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectThe Fraunhofer FDK AAC Codec Library for Android ("FDK AAC Codec") is software that implements
102228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectthe MPEG Advanced Audio Coding ("AAC") encoding and decoding scheme for digital audio.
112228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectThis FDK AAC Codec software is intended to be used on a wide variety of Android devices.
122228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
132228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectAAC's HE-AAC and HE-AAC v2 versions are regarded as today's most efficient general perceptual
142228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectaudio codecs. AAC-ELD is considered the best-performing full-bandwidth communications codec by
152228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectindependent studies and is widely deployed. AAC has been standardized by ISO and IEC as part
162228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectof the MPEG specifications.
172228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
182228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectPatent licenses for necessary patent claims for the FDK AAC Codec (including those of Fraunhofer)
192228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectmay be obtained through Via Licensing (www.vialicensing.com) or through the respective patent owners
202228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectindividually for the purpose of encoding or decoding bit streams in products that are compliant with
212228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectthe ISO/IEC MPEG audio standards. Please note that most manufacturers of Android devices already license
222228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectthese patent claims through Via Licensing or directly from the patent owners, and therefore FDK AAC Codec
232228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectsoftware may already be covered under those patent licenses when it is used for those licensed purposes only.
242228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
252228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectCommercially-licensed AAC software libraries, including floating-point versions with enhanced sound quality,
262228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectare also available from Fraunhofer. Users are encouraged to check the Fraunhofer website for additional
272228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectapplications information and documentation.
282228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
292228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project2.    COPYRIGHT LICENSE
302228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
312228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectRedistribution and use in source and binary forms, with or without modification, are permitted without
322228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectpayment of copyright license fees provided that you satisfy the following conditions:
332228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
342228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectYou must retain the complete text of this software license in redistributions of the FDK AAC Codec or
352228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectyour modifications thereto in source code form.
362228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
372228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectYou must retain the complete text of this software license in the documentation and/or other materials
382228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectprovided with redistributions of the FDK AAC Codec or your modifications thereto in binary form.
392228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectYou must make available free of charge copies of the complete source code of the FDK AAC Codec and your
402228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectmodifications thereto to recipients of copies in binary form.
412228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
422228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectThe name of Fraunhofer may not be used to endorse or promote products derived from this library without
432228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectprior written permission.
442228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
452228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectYou may not charge copyright license fees for anyone to use, copy or distribute the FDK AAC Codec
462228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectsoftware or your modifications thereto.
472228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
482228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectYour modified versions of the FDK AAC Codec must carry prominent notices stating that you changed the software
492228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectand the date of any change. For modified versions of the FDK AAC Codec, the term
502228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project"Fraunhofer FDK AAC Codec Library for Android" must be replaced by the term
512228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project"Third-Party Modified Version of the Fraunhofer FDK AAC Codec Library for Android."
522228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
532228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project3.    NO PATENT LICENSE
542228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
552228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectNO EXPRESS OR IMPLIED LICENSES TO ANY PATENT CLAIMS, including without limitation the patents of Fraunhofer,
562228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectARE GRANTED BY THIS SOFTWARE LICENSE. Fraunhofer provides no warranty of patent non-infringement with
572228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectrespect to this software.
582228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
592228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectYou may use this FDK AAC Codec software or modifications thereto only for purposes that are authorized
602228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectby appropriate patent licenses.
612228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
622228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project4.    DISCLAIMER
632228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
642228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectThis FDK AAC Codec software is provided by Fraunhofer on behalf of the copyright holders and contributors
652228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project"AS IS" and WITHOUT ANY EXPRESS OR IMPLIED WARRANTIES, including but not limited to the implied warranties
662228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectof merchantability and fitness for a particular purpose. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR
672228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectCONTRIBUTORS BE LIABLE for any direct, indirect, incidental, special, exemplary, or consequential damages,
682228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectincluding but not limited to procurement of substitute goods or services; loss of use, data, or profits,
692228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projector business interruption, however caused and on any theory of liability, whether in contract, strict
702228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectliability, or tort (including negligence), arising in any way out of the use of this software, even if
712228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectadvised of the possibility of such damage.
722228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
732228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project5.    CONTACT INFORMATION
742228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
752228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectFraunhofer Institute for Integrated Circuits IIS
762228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectAttention: Audio and Multimedia Departments - FDK AAC LL
772228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectAm Wolfsmantel 33
782228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project91058 Erlangen, Germany
792228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
802228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectwww.iis.fraunhofer.de/amm
812228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectamm-info@iis.fraunhofer.de
822228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project----------------------------------------------------------------------------------------------------------- */
832228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
842228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/******************************** MPEG Audio Encoder **************************
852228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
862228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project   Initial author:       M.Werner
872228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project   contents/description: Quantization
882228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
892228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project******************************************************************************/
902228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
912228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#include "quantize.h"
922228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
932228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#include "aacEnc_rom.h"
942228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
952228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/*****************************************************************************
962228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
972228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    functionname: FDKaacEnc_quantizeLines
982228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    description: quantizes spectrum lines
992228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    returns:
1002228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    input: global gain, number of lines to process, spectral data
1012228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    output: quantized spectrum
1022228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
1032228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project*****************************************************************************/
1042228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectstatic void FDKaacEnc_quantizeLines(INT      gain,
1052228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                          INT      noOfLines,
1062228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                          FIXP_DBL *mdctSpectrum,
1072228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                          SHORT      *quaSpectrum)
1082228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{
1092228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  int   line;
1102228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  FIXP_DBL k = FL2FXCONST_DBL(-0.0946f + 0.5f)>>16;
1112228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  FIXP_QTD quantizer = FDKaacEnc_quantTableQ[(-gain)&3];
1122228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  INT      quantizershift = ((-gain)>>2)+1;
1132228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
1142228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
1152228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  for (line = 0; line < noOfLines; line++)
1162228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  {
1172228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    FIXP_DBL accu = fMultDiv2(mdctSpectrum[line],quantizer);
1182228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
1192228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    if (accu < FL2FXCONST_DBL(0.0f))
1202228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    {
1212228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      accu=-accu;
1222228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      /* normalize */
1232228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      INT   accuShift = CntLeadingZeros(accu) - 1;  /* CountLeadingBits() is not necessary here since test value is always > 0 */
1242228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      accu <<= accuShift;
1252228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      INT tabIndex = (INT)(accu>>(DFRACT_BITS-2-MANT_DIGITS))&(~MANT_SIZE);
1262228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      INT totalShift = quantizershift-accuShift+1;
1272228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      accu = fMultDiv2(FDKaacEnc_mTab_3_4[tabIndex],FDKaacEnc_quantTableE[totalShift&3]);
1282228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      totalShift = (16-4)-(3*(totalShift>>2));
1292228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      FDK_ASSERT(totalShift >=0); /* MAX_QUANT_VIOLATION */
130ed247dfa54234c1115c35a05327955b29b161e8aJean-Michel Trivi      accu >>= fixMin(totalShift,DFRACT_BITS-1);
1312228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      quaSpectrum[line] = (SHORT)(-((LONG)(k + accu) >> (DFRACT_BITS-1-16)));
1322228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    }
1332228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    else if(accu > FL2FXCONST_DBL(0.0f))
1342228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    {
1352228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      /* normalize */
1362228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      INT   accuShift = CntLeadingZeros(accu) - 1;  /* CountLeadingBits() is not necessary here since test value is always > 0 */
1372228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      accu <<= accuShift;
1382228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      INT tabIndex = (INT)(accu>>(DFRACT_BITS-2-MANT_DIGITS))&(~MANT_SIZE);
1392228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      INT totalShift = quantizershift-accuShift+1;
1402228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      accu = fMultDiv2(FDKaacEnc_mTab_3_4[tabIndex],FDKaacEnc_quantTableE[totalShift&3]);
1412228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      totalShift = (16-4)-(3*(totalShift>>2));
1422228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      FDK_ASSERT(totalShift >=0); /* MAX_QUANT_VIOLATION */
143ed247dfa54234c1115c35a05327955b29b161e8aJean-Michel Trivi      accu >>= fixMin(totalShift,DFRACT_BITS-1);
1442228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      quaSpectrum[line] = (SHORT)((LONG)(k + accu) >> (DFRACT_BITS-1-16));
1452228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    }
1462228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    else
1472228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      quaSpectrum[line]=0;
1482228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  }
1492228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project}
1502228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
1512228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
1522228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/*****************************************************************************
1532228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
1542228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    functionname:iFDKaacEnc_quantizeLines
1552228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    description: iquantizes spectrum lines
1562228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                 mdctSpectrum = iquaSpectrum^4/3 *2^(0.25*gain)
1572228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    input: global gain, number of lines to process,quantized spectrum
1582228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    output: spectral data
1592228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
1602228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project*****************************************************************************/
1612228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectstatic void FDKaacEnc_invQuantizeLines(INT  gain,
1622228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                             INT  noOfLines,
1632228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                             SHORT *quantSpectrum,
1642228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                             FIXP_DBL *mdctSpectrum)
1652228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
1662228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{
1672228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  INT iquantizermod;
1682228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  INT iquantizershift;
1692228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  INT line;
1702228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
1712228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  iquantizermod = gain&3;
1722228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  iquantizershift = gain>>2;
1732228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
1742228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  for (line = 0; line < noOfLines; line++) {
1752228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
1762228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    if(quantSpectrum[line] < 0) {
1772228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      FIXP_DBL accu;
1782228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      INT ex,specExp,tabIndex;
1792228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      FIXP_DBL s,t;
1802228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
1812228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      accu = (FIXP_DBL) -quantSpectrum[line];
1822228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
1832228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      ex = CountLeadingBits(accu);
1842228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      accu <<= ex;
1852228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      specExp = (DFRACT_BITS-1) - ex;
1862228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
1872228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      FDK_ASSERT(specExp < 14);       /* this fails if abs(value) > 8191 */
1882228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
1892228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      tabIndex = (INT)(accu>>(DFRACT_BITS-2-MANT_DIGITS))&(~MANT_SIZE);
1902228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
1912228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      /* calculate "mantissa" ^4/3 */
1922228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      s = FDKaacEnc_mTab_4_3Elc[tabIndex];
1932228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
1942228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      /* get approperiate exponent multiplier for specExp^3/4 combined with scfMod */
1952228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      t = FDKaacEnc_specExpMantTableCombElc[iquantizermod][specExp];
1962228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
1972228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      /* multiply "mantissa" ^4/3 with exponent multiplier */
1982228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      accu = fMult(s,t);
1992228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
2002228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      /* get approperiate exponent shifter */
2012228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      specExp = FDKaacEnc_specExpTableComb[iquantizermod][specExp]-1; /* -1 to avoid overflows in accu */
2022228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
2032228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      if ((-iquantizershift-specExp) < 0)
2042228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project        accu <<= -(-iquantizershift-specExp);
2052228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      else
2062228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project        accu >>= -iquantizershift-specExp;
2072228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
2082228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      mdctSpectrum[line] = -accu;
2092228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    }
2102228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    else if (quantSpectrum[line] > 0) {
2112228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      FIXP_DBL accu;
2122228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      INT ex,specExp,tabIndex;
2132228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      FIXP_DBL s,t;
2142228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
2152228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      accu = (FIXP_DBL)(INT)quantSpectrum[line];
2162228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
2172228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      ex = CountLeadingBits(accu);
2182228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      accu <<= ex;
2192228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      specExp = (DFRACT_BITS-1) - ex;
2202228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
2212228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      FDK_ASSERT(specExp < 14);       /* this fails if abs(value) > 8191 */
2222228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
2232228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      tabIndex = (INT)(accu>>(DFRACT_BITS-2-MANT_DIGITS))&(~MANT_SIZE);
2242228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
2252228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      /* calculate "mantissa" ^4/3 */
2262228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      s = FDKaacEnc_mTab_4_3Elc[tabIndex];
2272228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
2282228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      /* get approperiate exponent multiplier for specExp^3/4 combined with scfMod */
2292228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      t = FDKaacEnc_specExpMantTableCombElc[iquantizermod][specExp];
2302228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
2312228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      /* multiply "mantissa" ^4/3 with exponent multiplier */
2322228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      accu = fMult(s,t);
2332228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
2342228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      /* get approperiate exponent shifter */
2352228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      specExp = FDKaacEnc_specExpTableComb[iquantizermod][specExp]-1; /* -1 to avoid overflows in accu */
2362228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
2372228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      if (( -iquantizershift-specExp) < 0)
2382228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project        accu <<= -(-iquantizershift-specExp);
2392228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      else
2402228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project        accu >>= -iquantizershift-specExp;
2412228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
2422228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      mdctSpectrum[line] = accu;
2432228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    }
2442228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    else {
2452228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project      mdctSpectrum[line] = FL2FXCONST_DBL(0.0f);
2462228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    }
2472228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  }
2482228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project}
2492228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
2502228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/*****************************************************************************
2512228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
2522228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    functionname: FDKaacEnc_QuantizeSpectrum
2532228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    description: quantizes the entire spectrum
2542228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    returns:
2552228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    input: number of scalefactor bands to be quantized, ...
2562228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    output: quantized spectrum
2572228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
2582228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project*****************************************************************************/
2592228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectvoid FDKaacEnc_QuantizeSpectrum(INT sfbCnt,
2602228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                      INT maxSfbPerGroup,
2612228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                      INT sfbPerGroup,
2622228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                      INT *sfbOffset,
2632228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                      FIXP_DBL *mdctSpectrum,
2642228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                      INT globalGain,
2652228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                      INT *scalefactors,
2662228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                      SHORT *quantizedSpectrum)
2672228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{
2682228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  INT sfbOffs,sfb;
2692228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
2702228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  /* in FDKaacEnc_quantizeLines quaSpectrum is calculated with:
2712228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project        spec^(3/4) * 2^(-3/16*QSS) * 2^(3/4*scale) + k
2722228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project     simplify scaling calculation and reduce QSS before:
2732228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project        spec^(3/4) * 2^(-3/16*(QSS - 4*scale)) */
2742228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
2752228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  for(sfbOffs=0;sfbOffs<sfbCnt;sfbOffs+=sfbPerGroup)
2762228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  for (sfb = 0; sfb < maxSfbPerGroup; sfb++)
2772228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  {
2782228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    INT scalefactor = scalefactors[sfbOffs+sfb] ;
2792228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
2802228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    FDKaacEnc_quantizeLines(globalGain - scalefactor, /* QSS */
2812228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                  sfbOffset[sfbOffs+sfb+1] - sfbOffset[sfbOffs+sfb],
2822228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                  mdctSpectrum + sfbOffset[sfbOffs+sfb],
2832228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                  quantizedSpectrum + sfbOffset[sfbOffs+sfb]);
2842228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  }
2852228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project}
2862228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
2872228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/*****************************************************************************
2882228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
2892228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    functionname: FDKaacEnc_calcSfbDist
2902228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    description: calculates distortion of quantized values
2912228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    returns: distortion
2922228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    input: gain, number of lines to process, spectral data
2932228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    output:
2942228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
2952228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project*****************************************************************************/
2962228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectFIXP_DBL FDKaacEnc_calcSfbDist(FIXP_DBL *mdctSpectrum,
2972228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                     SHORT *quantSpectrum,
2982228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                     INT noOfLines,
2992228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                     INT gain
3002228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                     )
3012228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{
3022228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  INT i,scale;
3032228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  FIXP_DBL xfsf;
3042228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  FIXP_DBL diff;
3052228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  FIXP_DBL invQuantSpec;
3062228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
3072228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  xfsf = FL2FXCONST_DBL(0.0f);
3082228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
3092228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  for (i=0; i<noOfLines; i++) {
3102228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    /* quantization */
3112228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    FDKaacEnc_quantizeLines(gain,
3122228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                  1,
3132228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                 &mdctSpectrum[i],
3142228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                 &quantSpectrum[i]);
3152228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
316381d69840ad3af2259f0b7ef49236f9ee9c76b76Jean-Michel Trivi    if (fAbs(quantSpectrum[i])>MAX_QUANT) {
317381d69840ad3af2259f0b7ef49236f9ee9c76b76Jean-Michel Trivi      return FL2FXCONST_DBL(0.0f);
318381d69840ad3af2259f0b7ef49236f9ee9c76b76Jean-Michel Trivi    }
3192228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    /* inverse quantization */
3202228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    FDKaacEnc_invQuantizeLines(gain,1,&quantSpectrum[i],&invQuantSpec);
3212228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
3222228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    /* dist */
3232228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    diff = fixp_abs(fixp_abs(invQuantSpec) - fixp_abs(mdctSpectrum[i]>>1));
3242228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
3252228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    scale = CountLeadingBits(diff);
3262228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    diff = scaleValue(diff, scale);
3272228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    diff = fPow2(diff);
3282228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    scale = fixMin(2*(scale-1), DFRACT_BITS-1);
3292228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
3302228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    diff = scaleValue(diff, -scale);
3312228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
3322228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    xfsf = xfsf + diff;
3332228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  }
3342228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
3352228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  xfsf = CalcLdData(xfsf);
3362228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
3372228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  return xfsf;
3382228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project}
3392228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
3402228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/*****************************************************************************
3412228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
3422228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    functionname: FDKaacEnc_calcSfbQuantEnergyAndDist
3432228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    description: calculates energy and distortion of quantized values
3442228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    returns:
3452228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    input: gain, number of lines to process, quantized spectral data,
3462228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project           spectral data
3472228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    output: energy, distortion
3482228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
3492228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project*****************************************************************************/
3502228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectvoid FDKaacEnc_calcSfbQuantEnergyAndDist(FIXP_DBL *mdctSpectrum,
3512228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                               SHORT *quantSpectrum,
3522228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                               INT noOfLines,
3532228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                               INT gain,
3542228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                               FIXP_DBL *en,
3552228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project                               FIXP_DBL *dist)
3562228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{
3572228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  INT i,scale;
3582228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  FIXP_DBL invQuantSpec;
3592228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  FIXP_DBL diff;
3602228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
361381d69840ad3af2259f0b7ef49236f9ee9c76b76Jean-Michel Trivi  FIXP_DBL energy = FL2FXCONST_DBL(0.0f);
362381d69840ad3af2259f0b7ef49236f9ee9c76b76Jean-Michel Trivi  FIXP_DBL distortion = FL2FXCONST_DBL(0.0f);
3632228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
3642228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  for (i=0; i<noOfLines; i++) {
365381d69840ad3af2259f0b7ef49236f9ee9c76b76Jean-Michel Trivi
366381d69840ad3af2259f0b7ef49236f9ee9c76b76Jean-Michel Trivi    if (fAbs(quantSpectrum[i])>MAX_QUANT) {
367381d69840ad3af2259f0b7ef49236f9ee9c76b76Jean-Michel Trivi      *en   = FL2FXCONST_DBL(0.0f);
368381d69840ad3af2259f0b7ef49236f9ee9c76b76Jean-Michel Trivi      *dist = FL2FXCONST_DBL(0.0f);
369381d69840ad3af2259f0b7ef49236f9ee9c76b76Jean-Michel Trivi      return;
370381d69840ad3af2259f0b7ef49236f9ee9c76b76Jean-Michel Trivi    }
371381d69840ad3af2259f0b7ef49236f9ee9c76b76Jean-Michel Trivi
3722228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    /* inverse quantization */
3732228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    FDKaacEnc_invQuantizeLines(gain,1,&quantSpectrum[i],&invQuantSpec);
3742228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
3752228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    /* energy */
376381d69840ad3af2259f0b7ef49236f9ee9c76b76Jean-Michel Trivi    energy += fPow2(invQuantSpec);
3772228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
3782228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    /* dist */
3792228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    diff = fixp_abs(fixp_abs(invQuantSpec) - fixp_abs(mdctSpectrum[i]>>1));
3802228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
3812228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    scale = CountLeadingBits(diff);
3822228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    diff = scaleValue(diff, scale);
3832228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    diff = fPow2(diff);
3842228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
3852228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    scale = fixMin(2*(scale-1), DFRACT_BITS-1);
3862228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
3872228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project    diff = scaleValue(diff, -scale);
3882228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
389381d69840ad3af2259f0b7ef49236f9ee9c76b76Jean-Michel Trivi    distortion += diff;
3902228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project  }
3912228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
392381d69840ad3af2259f0b7ef49236f9ee9c76b76Jean-Michel Trivi  *en   = CalcLdData(energy)+FL2FXCONST_DBL(0.03125f);
393381d69840ad3af2259f0b7ef49236f9ee9c76b76Jean-Michel Trivi  *dist = CalcLdData(distortion);
3942228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project}
3952228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project
396