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: A. Horndasch (code originally from lwr) / Josef Hoepfl (FDK) 872228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project contents/description: intensity stereo processing 882228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 892228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project******************************************************************************/ 902228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 912228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#include "intensity.h" 922228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#include "interface.h" 932228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#include "psy_configuration.h" 942228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#include "psy_const.h" 952228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#include "qc_main.h" 962228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#include "bit_cnt.h" 972228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 982228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/* only set an IS seed it left/right channel correlation is above IS_CORR_THRESH */ 992228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#define IS_CORR_THRESH FL2FXCONST_DBL(0.95f) 1002228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1012228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/* when expanding the IS region to more SFBs only accept an error that is 1022228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * not more than IS_TOTAL_ERROR_THRESH overall and 1032228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * not more than IS_LOCAL_ERROR_THRESH for the current SFB */ 1042228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#define IS_TOTAL_ERROR_THRESH FL2FXCONST_DBL(0.04f) 1052228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#define IS_LOCAL_ERROR_THRESH FL2FXCONST_DBL(0.01f) 1062228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1072228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/* the maximum allowed change of the intensity direction (unit: IS scale) - scaled with factor 0.25 - */ 1082228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#define IS_DIRECTION_DEVIATION_THRESH_SF 2 1092228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#define IS_DIRECTION_DEVIATION_THRESH FL2FXCONST_DBL(2.0f/(1<<IS_DIRECTION_DEVIATION_THRESH_SF)) 1102228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1112228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/* IS regions need to have a minimal percentage of the overall loudness, e.g. 0.06 == 6% */ 1122228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#define IS_REGION_MIN_LOUDNESS FL2FXCONST_DBL(0.1f) 1132228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1142228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/* only perform IS if IS_MIN_SFBS neighboring SFBs can be processed */ 1152228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#define IS_MIN_SFBS 6 1162228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1172228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/* only do IS if 1182228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * if IS_LEFT_RIGHT_RATIO_THRESH < sfbEnergyLeft[sfb]/sfbEnergyRight[sfb] < 1 / IS_LEFT_RIGHT_RATIO_THRESH 1192228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * -> no IS if the panning angle is not far from the middle, MS will do */ 1202228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/* this is equivalent to a scale of +/-1.02914634566 */ 1212228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#define IS_LEFT_RIGHT_RATIO_THRESH FL2FXCONST_DBL(0.7f) 1222228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1232228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/* scalefactor of realScale */ 1242228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#define REAL_SCALE_SF 1 1252228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1262228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/* scalefactor overallLoudness */ 1272228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#define OVERALL_LOUDNESS_SF 6 1282228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1292228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/* scalefactor for sum over max samples per goup */ 1302228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#define MAX_SFB_PER_GROUP_SF 6 1312228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1322228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/* scalefactor for sum of mdct spectrum */ 1332228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#define MDCT_SPEC_SF 6 1342228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1352228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1362228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projecttypedef struct 1372228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 1382228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1392228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL corr_thresh; /*!< Only set an IS seed it left/right channel correlation is above corr_thresh */ 1402228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1412228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL total_error_thresh; /*!< When expanding the IS region to more SFBs only accept an error that is 1422228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project not more than 'total_error_thresh' overall. */ 1432228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1442228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL local_error_thresh; /*!< When expanding the IS region to more SFBs only accept an error that is 1452228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project not more than 'local_error_thresh' for the current SFB. */ 1462228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1472228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL direction_deviation_thresh; /*!< The maximum allowed change of the intensity direction (unit: IS scale) */ 1482228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1492228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL is_region_min_loudness; /*!< IS regions need to have a minimal percentage of the overall loudness, e.g. 0.06 == 6% */ 1502228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1512228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT min_is_sfbs; /*!< Only perform IS if 'min_is_sfbs' neighboring SFBs can be processed */ 1522228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1532228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL left_right_ratio_threshold; /*!< No IS if the panning angle is not far from the middle, MS will do */ 1542228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1552228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} INTENSITY_PARAMETERS; 1562228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1572228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1582228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/***************************************************************************** 1592228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1602228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project functionname: calcSfbMaxScale 1612228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1622228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project description: Calc max value in scalefactor band 1632228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1642228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project input: *mdctSpectrum 1652228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project l1 1662228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project l2 1672228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1682228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project output: none 1692228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1702228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project returns: scalefactor 1712228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1722228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project*****************************************************************************/ 1732228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectstatic INT 1742228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectcalcSfbMaxScale(const FIXP_DBL *mdctSpectrum, 1752228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const INT l1, 1762228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const INT l2) 1772228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 1782228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT i; 1792228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT sfbMaxScale; 1802228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL maxSpc; 1812228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1822228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project maxSpc = FL2FXCONST_DBL(0.0); 1832228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for (i=l1; i<l2; i++) { 1842228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL tmp = fixp_abs((FIXP_DBL)mdctSpectrum[i]); 1852228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project maxSpc = fixMax(maxSpc, tmp); 1862228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 1872228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbMaxScale = (maxSpc==FL2FXCONST_DBL(0.0)) ? (DFRACT_BITS-2) : CntLeadingZeros(maxSpc)-1; 1882228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1892228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project return sfbMaxScale; 1902228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 1912228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1922228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1932228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/***************************************************************************** 1942228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1952228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project functionname: FDKaacEnc_initIsParams 1962228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1972228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project description: Initialization of intensity parameters 1982228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1992228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project input: isParams 2002228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2012228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project output: isParams 2022228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2032228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project returns: none 2042228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2052228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project*****************************************************************************/ 2062228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectstatic void 2072228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectFDKaacEnc_initIsParams(INTENSITY_PARAMETERS *isParams) 2082228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 2092228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project isParams->corr_thresh = IS_CORR_THRESH; 2102228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project isParams->total_error_thresh = IS_TOTAL_ERROR_THRESH; 2112228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project isParams->local_error_thresh = IS_LOCAL_ERROR_THRESH; 2122228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project isParams->direction_deviation_thresh = IS_DIRECTION_DEVIATION_THRESH; 2132228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project isParams->is_region_min_loudness = IS_REGION_MIN_LOUDNESS; 2142228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project isParams->min_is_sfbs = IS_MIN_SFBS; 2152228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project isParams->left_right_ratio_threshold = IS_LEFT_RIGHT_RATIO_THRESH; 2162228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} 2172228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2182228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2192228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/***************************************************************************** 2202228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2212228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project functionname: FDKaacEnc_prepareIntensityDecision 2222228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2232228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project description: Prepares intensity decision 2242228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2252228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project input: sfbEnergyLeft 2262228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbEnergyRight 2272228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbEnergyLdDataLeft 2282228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbEnergyLdDataRight 2292228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project mdctSpectrumLeft 2302228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbEnergyLdDataRight 2312228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project isParams 2322228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2332228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project output: hrrErr scale: none 2342228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project isMask scale: none 2352228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project realScale scale: LD_DATA_SHIFT + REAL_SCALE_SF 2362228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project normSfbLoudness scale: none 2372228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2382228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project returns: none 2392228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2402228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project*****************************************************************************/ 2412228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectstatic void 2422228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectFDKaacEnc_prepareIntensityDecision(const FIXP_DBL *sfbEnergyLeft, 2432228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const FIXP_DBL *sfbEnergyRight, 2442228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const FIXP_DBL *sfbEnergyLdDataLeft, 2452228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const FIXP_DBL *sfbEnergyLdDataRight, 2462228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const FIXP_DBL *mdctSpectrumLeft, 2472228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const FIXP_DBL *mdctSpectrumRight, 2482228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const INTENSITY_PARAMETERS *isParams, 2492228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *hrrErr, 2502228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT *isMask, 2512228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *realScale, 2522228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *normSfbLoudness, 2532228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const INT sfbCnt, 2542228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const INT sfbPerGroup, 2552228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const INT maxSfbPerGroup, 2562228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const INT *sfbOffset) 2572228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 2582228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT j,sfb,sfboffs; 2592228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT grpCounter; 2602228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2612228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* temporary variables to compute loudness */ 2622228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL overallLoudness[MAX_NO_OF_GROUPS]; 2632228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2642228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* temporary variables to compute correlation */ 2652228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL channelCorr[MAX_GROUPED_SFB]; 2662228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL ml, mr; 2672228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL prod_lr; 2682228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL square_l, square_r; 2692228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL tmp_l, tmp_r; 2702228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL inv_n; 2712228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2722228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDKmemclear(channelCorr, MAX_GROUPED_SFB*sizeof(FIXP_DBL)); 2732228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDKmemclear(normSfbLoudness, MAX_GROUPED_SFB*sizeof(FIXP_DBL)); 2742228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDKmemclear(overallLoudness, MAX_NO_OF_GROUPS*sizeof(FIXP_DBL)); 2752228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDKmemclear(realScale, MAX_GROUPED_SFB*sizeof(FIXP_DBL)); 2762228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2772228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for (grpCounter = 0, sfboffs = 0; sfboffs < sfbCnt; sfboffs += sfbPerGroup, grpCounter++) { 2782228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project overallLoudness[grpCounter] = FL2FXCONST_DBL(0.0f); 2792228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for (sfb = 0; sfb < maxSfbPerGroup; sfb++) { 2802228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT sL,sR,s; 2812228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL isValue = sfbEnergyLdDataLeft[sfb+sfboffs]-sfbEnergyLdDataRight[sfb+sfboffs]; 2822228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2832228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* delimitate intensity scale value to representable range */ 2842228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project realScale[sfb + sfboffs] = fixMin(FL2FXCONST_DBL(60.f/(1<<(REAL_SCALE_SF+LD_DATA_SHIFT))), fixMax(FL2FXCONST_DBL(-60.f/(1<<(REAL_SCALE_SF+LD_DATA_SHIFT))), isValue)); 2852228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2862228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sL = fixMax(0,(CntLeadingZeros(sfbEnergyLeft[sfb + sfboffs])-1)); 2872228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sR = fixMax(0,(CntLeadingZeros(sfbEnergyRight[sfb + sfboffs])-1)); 2882228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project s = (fixMin(sL,sR)>>2)<<2; 2892228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project normSfbLoudness[sfb + sfboffs] = sqrtFixp(sqrtFixp(((sfbEnergyLeft[sfb + sfboffs]<<s) >> 1) + ((sfbEnergyRight[sfb + sfboffs]<<s) >> 1))) >> (s>>2); 2902228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2912228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project overallLoudness[grpCounter] += normSfbLoudness[sfb + sfboffs] >> OVERALL_LOUDNESS_SF; 2922228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* don't do intensity if 2932228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * - panning angle is too close to the middle or 2942228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * - one channel is non-existent or 2952228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * - if it is dual mono */ 2962228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if( (sfbEnergyLeft[sfb + sfboffs] >= fMult(isParams->left_right_ratio_threshold,sfbEnergyRight[sfb + sfboffs])) 2972228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project && (fMult(isParams->left_right_ratio_threshold,sfbEnergyLeft[sfb + sfboffs]) <= sfbEnergyRight[sfb + sfboffs]) ) { 2982228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2992228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* this will prevent post processing from considering this SFB for merging */ 3002228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project hrrErr[sfb + sfboffs] = FL2FXCONST_DBL(1.0/8.0); 3012228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 3022228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 3032228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 3042228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3052228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for (grpCounter = 0, sfboffs = 0; sfboffs < sfbCnt; sfboffs += sfbPerGroup, grpCounter++) { 3062228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT invOverallLoudnessSF; 3072228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL invOverallLoudness; 3082228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3092228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (overallLoudness[grpCounter] == FL2FXCONST_DBL(0.0)) { 3102228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project invOverallLoudness = FL2FXCONST_DBL(0.0); 3112228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project invOverallLoudnessSF = 0; 3122228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 3132228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project else { 3142228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project invOverallLoudness = fDivNorm((FIXP_DBL)MAXVAL_DBL, overallLoudness[grpCounter],&invOverallLoudnessSF); 3152228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project invOverallLoudnessSF = invOverallLoudnessSF - OVERALL_LOUDNESS_SF + 1; /* +1: compensate fMultDiv2() in subsequent loop */ 3162228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 3172228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project invOverallLoudnessSF = fixMin(fixMax(invOverallLoudnessSF,-(DFRACT_BITS-1)),DFRACT_BITS-1); 3182228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3192228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for (sfb = 0; sfb < maxSfbPerGroup; sfb++) { 3202228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL tmp; 3212228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3222228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project tmp = fMultDiv2((normSfbLoudness[sfb + sfboffs]>>OVERALL_LOUDNESS_SF)<<OVERALL_LOUDNESS_SF,invOverallLoudness); 3232228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3242228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project normSfbLoudness[sfb + sfboffs] = scaleValue(tmp, invOverallLoudnessSF); 3252228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3262228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project channelCorr[sfb + sfboffs] = FL2FXCONST_DBL(0.0f); 3272228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3282228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDK_ASSERT(50 >= 49); 3292228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* max width of scalefactorband is 96; width's are always even */ 3302228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* inv_n is scaled with factor 2 to compensate fMultDiv2() in subsequent loops */ 3312228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project inv_n = GetInvInt((sfbOffset[sfb + sfboffs + 1] - sfbOffset[sfb + sfboffs])>>1); 3322228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3332228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (inv_n > FL2FXCONST_DBL(0.0f)) { 3342228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT s,sL,sR; 3352228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3362228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* correlation := Pearson's product-moment coefficient */ 3372228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* compute correlation between channels and check if it is over threshold */ 3382228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project ml = FL2FXCONST_DBL(0.0f); 3392228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project mr = FL2FXCONST_DBL(0.0f); 3402228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project prod_lr = FL2FXCONST_DBL(0.0f); 3412228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project square_l = FL2FXCONST_DBL(0.0f); 3422228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project square_r = FL2FXCONST_DBL(0.0f); 3432228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3442228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sL = calcSfbMaxScale(mdctSpectrumLeft,sfbOffset[sfb+sfboffs],sfbOffset[sfb+sfboffs+1]); 3452228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sR = calcSfbMaxScale(mdctSpectrumRight,sfbOffset[sfb+sfboffs],sfbOffset[sfb+sfboffs+1]); 3462228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project s = fixMin(sL,sR); 3472228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3482228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for (j = sfbOffset[sfb + sfboffs]; j < sfbOffset[sfb + sfboffs + 1]; j++) { 3492228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project ml += fMultDiv2((mdctSpectrumLeft[j] << s),inv_n); // scaled with mdctScale - s + inv_n 3502228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project mr += fMultDiv2((mdctSpectrumRight[j] << s),inv_n); // scaled with mdctScale - s + inv_n 3512228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 3522228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project ml = fMultDiv2(ml,inv_n); // scaled with mdctScale - s + inv_n 3532228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project mr = fMultDiv2(mr,inv_n); // scaled with mdctScale - s + inv_n 3542228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3552228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for (j = sfbOffset[sfb + sfboffs]; j < sfbOffset[sfb + sfboffs + 1]; j++) { 3562228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project tmp_l = fMultDiv2((mdctSpectrumLeft[j] << s),inv_n) - ml; // scaled with mdctScale - s + inv_n 3572228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project tmp_r = fMultDiv2((mdctSpectrumRight[j] << s),inv_n) - mr; // scaled with mdctScale - s + inv_n 3582228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3592228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project prod_lr += fMultDiv2(tmp_l,tmp_r); // scaled with 2*(mdctScale - s + inv_n) + 1 3602228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project square_l += fPow2Div2(tmp_l); // scaled with 2*(mdctScale - s + inv_n) + 1 3612228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project square_r += fPow2Div2(tmp_r); // scaled with 2*(mdctScale - s + inv_n) + 1 3622228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 3632228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project prod_lr = prod_lr << 1; // scaled with 2*(mdctScale - s + inv_n) 3642228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project square_l = square_l << 1; // scaled with 2*(mdctScale - s + inv_n) 3652228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project square_r = square_r << 1; // scaled with 2*(mdctScale - s + inv_n) 3662228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3672228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (square_l > FL2FXCONST_DBL(0.0f) && square_r > FL2FXCONST_DBL(0.0f)) { 3682228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT channelCorrSF = 0; 3692228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3702228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* local scaling of square_l and square_r is compensated after sqrt calculation */ 3712228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sL = fixMax(0,(CntLeadingZeros(square_l)-1)); 3722228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sR = fixMax(0,(CntLeadingZeros(square_r)-1)); 3732228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project s = ((sL + sR)>>1)<<1; 3742228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sL = fixMin(sL,s); 3752228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sR = s-sL; 3762228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project tmp = fMult(square_l<<sL,square_r<<sR); 3772228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project tmp = sqrtFixp(tmp); 3782228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3792228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDK_ASSERT(tmp > FL2FXCONST_DBL(0.0f)); 3802228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3812228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* numerator and denominator have the same scaling */ 3822228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (prod_lr < FL2FXCONST_DBL(0.0f) ) { 3832228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project channelCorr[sfb + sfboffs] = -(fDivNorm(-prod_lr,tmp,&channelCorrSF)); 3842228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3852228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 3862228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project else { 3872228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project channelCorr[sfb + sfboffs] = (fDivNorm( prod_lr,tmp,&channelCorrSF)); 3882228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 3892228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project channelCorrSF = fixMin(fixMax(( channelCorrSF + ((sL+sR)>>1)),-(DFRACT_BITS-1)),DFRACT_BITS-1); 3902228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3912228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (channelCorrSF < 0) { 3922228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project channelCorr[sfb + sfboffs] = channelCorr[sfb + sfboffs] >> (-channelCorrSF); 3932228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 3942228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project else { 3952228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* avoid overflows due to limited computational accuracy */ 3962228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if ( fAbs(channelCorr[sfb + sfboffs]) > (((FIXP_DBL)MAXVAL_DBL)>>channelCorrSF) ) { 3972228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (channelCorr[sfb + sfboffs] < FL2FXCONST_DBL(0.0f)) 3982228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project channelCorr[sfb + sfboffs] = -(FIXP_DBL) MAXVAL_DBL; 3992228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project else 4002228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project channelCorr[sfb + sfboffs] = (FIXP_DBL) MAXVAL_DBL; 4012228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 4022228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project else { 4032228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project channelCorr[sfb + sfboffs] = channelCorr[sfb + sfboffs] << channelCorrSF; 4042228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 4052228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 4062228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 4072228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 4082228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4092228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* for post processing: hrrErr is the error in terms of (too little) correlation 4102228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * weighted with the loudness of the SFB; SFBs with small hrrErr can be merged */ 4112228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (hrrErr[sfb + sfboffs] == FL2FXCONST_DBL(1.0/8.0)) { 4122228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project continue; 4132228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 4142228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4152228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project hrrErr[sfb + sfboffs] = fMultDiv2((FL2FXCONST_DBL(0.25f)-(channelCorr[sfb + sfboffs]>>2)),normSfbLoudness[sfb + sfboffs]); 4162228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4172228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* set IS mask/vector to 1, if correlation is high enough */ 4182228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (fAbs(channelCorr[sfb + sfboffs]) >= isParams->corr_thresh) { 4192228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project isMask[sfb + sfboffs] = 1; 4202228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 4212228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 4222228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 4232228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} 4242228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4252228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4262228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/***************************************************************************** 4272228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4282228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project functionname: FDKaacEnc_finalizeIntensityDecision 4292228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4302228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project description: Finalizes intensity decision 4312228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4322228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project input: isParams scale: none 4332228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project hrrErr scale: none 4342228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project realIsScale scale: LD_DATA_SHIFT + REAL_SCALE_SF 4352228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project normSfbLoudness scale: none 4362228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4372228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project output: isMask scale: none 4382228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4392228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project returns: none 4402228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4412228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project*****************************************************************************/ 4422228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectstatic void 4432228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectFDKaacEnc_finalizeIntensityDecision(const FIXP_DBL *hrrErr, 4442228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT *isMask, 4452228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const FIXP_DBL *realIsScale, 4462228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const FIXP_DBL *normSfbLoudness, 4472228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const INTENSITY_PARAMETERS *isParams, 4482228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const INT sfbCnt, 4492228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const INT sfbPerGroup, 4502228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const INT maxSfbPerGroup) 4512228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 4522228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT sfb,sfboffs, j; 4532228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL isScaleLast = FL2FXCONST_DBL(0.0f); 45460a27b7fe5bf2c3b07ae7f679ca67c0c5cba3b00Jean-Michel Trivi INT isStartValueFound = 0; 4552228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4562228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for (sfboffs = 0; sfboffs < sfbCnt; sfboffs += sfbPerGroup) { 45760a27b7fe5bf2c3b07ae7f679ca67c0c5cba3b00Jean-Michel Trivi INT startIsSfb = 0; 45860a27b7fe5bf2c3b07ae7f679ca67c0c5cba3b00Jean-Michel Trivi INT inIsBlock = 0; 45960a27b7fe5bf2c3b07ae7f679ca67c0c5cba3b00Jean-Michel Trivi INT currentIsSfbCount = 0; 46060a27b7fe5bf2c3b07ae7f679ca67c0c5cba3b00Jean-Michel Trivi FIXP_DBL overallHrrError = FL2FXCONST_DBL(0.0f); 46160a27b7fe5bf2c3b07ae7f679ca67c0c5cba3b00Jean-Michel Trivi FIXP_DBL isRegionLoudness = FL2FXCONST_DBL(0.0f); 46260a27b7fe5bf2c3b07ae7f679ca67c0c5cba3b00Jean-Michel Trivi 4632228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for (sfb = 0; sfb < maxSfbPerGroup; sfb++) { 4642228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (isMask[sfboffs + sfb] == 1) { 4652228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (currentIsSfbCount == 0) { 4662228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project startIsSfb = sfboffs + sfb; 46760a27b7fe5bf2c3b07ae7f679ca67c0c5cba3b00Jean-Michel Trivi } 46860a27b7fe5bf2c3b07ae7f679ca67c0c5cba3b00Jean-Michel Trivi if (isStartValueFound==0) { 4692228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project isScaleLast = realIsScale[sfboffs + sfb]; 47060a27b7fe5bf2c3b07ae7f679ca67c0c5cba3b00Jean-Michel Trivi isStartValueFound = 1; 4712228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 4722228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project inIsBlock = 1; 4732228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project currentIsSfbCount++; 4742228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project overallHrrError += hrrErr[sfboffs + sfb] >> (MAX_SFB_PER_GROUP_SF-3); 4752228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project isRegionLoudness += normSfbLoudness[sfboffs + sfb] >> MAX_SFB_PER_GROUP_SF; 4762228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 4772228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project else { 4782228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* based on correlation, IS should not be used 4792228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * -> use it anyway, if overall error is below threshold 4802228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * and if local error does not exceed threshold 4812228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * otherwise: check if there are enough IS SFBs 4822228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project */ 4832228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (inIsBlock) { 4842228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project overallHrrError += hrrErr[sfboffs + sfb] >> (MAX_SFB_PER_GROUP_SF-3); 4852228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project isRegionLoudness += normSfbLoudness[sfboffs + sfb] >> MAX_SFB_PER_GROUP_SF; 4862228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4872228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if ( (hrrErr[sfboffs + sfb] < (isParams->local_error_thresh>>3)) && (overallHrrError < (isParams->total_error_thresh>>MAX_SFB_PER_GROUP_SF)) ) { 4882228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project currentIsSfbCount++; 4892228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* overwrite correlation based decision */ 4902228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project isMask[sfboffs + sfb] = 1; 4912228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } else { 4922228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project inIsBlock = 0; 4932228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 4942228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 4952228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 4962228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* check for large direction deviation */ 4972228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (inIsBlock) { 4982228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if( fAbs(isScaleLast-realIsScale[sfboffs + sfb]) < (isParams->direction_deviation_thresh>>(REAL_SCALE_SF+LD_DATA_SHIFT-IS_DIRECTION_DEVIATION_THRESH_SF)) ) { 4992228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project isScaleLast = realIsScale[sfboffs + sfb]; 5002228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 5012228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project else{ 5022228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project isMask[sfboffs + sfb] = 0; 5032228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project inIsBlock = 0; 5042228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project currentIsSfbCount--; 5052228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 5062228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 5072228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 5082228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (currentIsSfbCount > 0 && (!inIsBlock || sfb == maxSfbPerGroup - 1)) { 5092228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* not enough SFBs -> do not use IS */ 5102228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (currentIsSfbCount < isParams->min_is_sfbs || (isRegionLoudness < isParams->is_region_min_loudness>>MAX_SFB_PER_GROUP_SF)) { 5112228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for(j = startIsSfb; j <= sfboffs + sfb; j++) { 5122228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project isMask[j] = 0; 5132228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 51460a27b7fe5bf2c3b07ae7f679ca67c0c5cba3b00Jean-Michel Trivi isScaleLast = FL2FXCONST_DBL(0.0f); 51560a27b7fe5bf2c3b07ae7f679ca67c0c5cba3b00Jean-Michel Trivi isStartValueFound = 0; 51660a27b7fe5bf2c3b07ae7f679ca67c0c5cba3b00Jean-Michel Trivi for (j=0; j < startIsSfb; j++) { 51760a27b7fe5bf2c3b07ae7f679ca67c0c5cba3b00Jean-Michel Trivi if (isMask[j]!=0) { 51860a27b7fe5bf2c3b07ae7f679ca67c0c5cba3b00Jean-Michel Trivi isScaleLast = realIsScale[j]; 51960a27b7fe5bf2c3b07ae7f679ca67c0c5cba3b00Jean-Michel Trivi isStartValueFound = 1; 52060a27b7fe5bf2c3b07ae7f679ca67c0c5cba3b00Jean-Michel Trivi } 52160a27b7fe5bf2c3b07ae7f679ca67c0c5cba3b00Jean-Michel Trivi } 5222228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 5232228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project currentIsSfbCount = 0; 5242228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project overallHrrError = FL2FXCONST_DBL(0.0f); 5252228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project isRegionLoudness = FL2FXCONST_DBL(0.0f); 5262228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 5272228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 5282228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 5292228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} 5302228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 5312228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 5322228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/***************************************************************************** 5332228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 5342228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project functionname: FDKaacEnc_IntensityStereoProcessing 5352228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 5362228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project description: Intensity stereo processing tool 5372228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 5382228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project input: sfbEnergyLeft 5392228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbEnergyRight 5402228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project mdctSpectrumLeft 5412228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project mdctSpectrumRight 5422228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbThresholdLeft 5432228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbThresholdRight 5442228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbSpreadEnLeft 5452228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbSpreadEnRight 5462228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbEnergyLdDataLeft 5472228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbEnergyLdDataRight 5482228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 5492228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project output: isBook 5502228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project isScale 5512228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pnsData->pnsFlag 5522228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project msDigest zeroed from start to sfbCnt 5532228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project msMask zeroed from start to sfbCnt 5542228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project mdctSpectrumRight zeroed where isBook!=0 5552228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbEnergyRight zeroed where isBook!=0 5562228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbSpreadEnRight zeroed where isBook!=0 5572228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbThresholdRight zeroed where isBook!=0 5582228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbEnergyLdDataRight FL2FXCONST_DBL(-1.0) where isBook!=0 5592228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbThresholdLdDataRight FL2FXCONST_DBL(-0.515625f) where isBook!=0 5602228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 5612228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project returns: none 5622228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 5632228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project*****************************************************************************/ 5642228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectvoid FDKaacEnc_IntensityStereoProcessing( 5652228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *sfbEnergyLeft, 5662228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *sfbEnergyRight, 5672228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *mdctSpectrumLeft, 5682228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *mdctSpectrumRight, 5692228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *sfbThresholdLeft, 5702228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *sfbThresholdRight, 5712228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *sfbThresholdLdDataRight, 5722228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *sfbSpreadEnLeft, 5732228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *sfbSpreadEnRight, 5742228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *sfbEnergyLdDataLeft, 5752228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *sfbEnergyLdDataRight, 5762228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT *msDigest, 5772228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT *msMask, 5782228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const INT sfbCnt, 5792228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const INT sfbPerGroup, 5802228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const INT maxSfbPerGroup, 5812228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const INT *sfbOffset, 5822228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const INT allowIS, 5832228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT *isBook, 5842228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT *isScale, 5852228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project PNS_DATA *RESTRICT pnsData[2] 5862228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project ) 5872228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 5882228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT sfb,sfboffs, j; 5892228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL scale; 5902228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL lr; 5912228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL hrrErr[MAX_GROUPED_SFB]; 5922228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL normSfbLoudness[MAX_GROUPED_SFB]; 5932228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL realIsScale[MAX_GROUPED_SFB]; 5942228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INTENSITY_PARAMETERS isParams; 5952228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT isMask[MAX_GROUPED_SFB]; 5962228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 5972228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDKmemclear((void*)isBook,sfbCnt*sizeof(INT)); 5982228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDKmemclear((void*)isMask,sfbCnt*sizeof(INT)); 5992228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDKmemclear((void*)realIsScale,sfbCnt*sizeof(FIXP_DBL)); 6002228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDKmemclear((void*)isScale,sfbCnt*sizeof(INT)); 6012228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDKmemclear((void*)hrrErr,sfbCnt*sizeof(FIXP_DBL)); 6022228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 6032228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (!allowIS) 6042228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project return; 6052228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 6062228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDKaacEnc_initIsParams(&isParams); 6072228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 6082228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* compute / set the following values per SFB: 6092228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * - left/right ratio between channels 6102228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * - normalized loudness 6112228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * + loudness == average of energy in channels to 0.25 6122228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * + normalization: division by sum of all SFB loudnesses 6132228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * - isMask (is set to 0 if channels are the same or one is 0) 6142228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project */ 6152228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDKaacEnc_prepareIntensityDecision(sfbEnergyLeft, 6162228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbEnergyRight, 6172228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbEnergyLdDataLeft, 6182228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbEnergyLdDataRight, 6192228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project mdctSpectrumLeft, 6202228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project mdctSpectrumRight, 6212228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project &isParams, 6222228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project hrrErr, 6232228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project isMask, 6242228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project realIsScale, 6252228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project normSfbLoudness, 6262228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbCnt, 6272228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbPerGroup, 6282228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project maxSfbPerGroup, 6292228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbOffset); 6302228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 6312228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDKaacEnc_finalizeIntensityDecision(hrrErr, 6322228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project isMask, 6332228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project realIsScale, 6342228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project normSfbLoudness, 6352228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project &isParams, 6362228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbCnt, 6372228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbPerGroup, 6382228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project maxSfbPerGroup); 6392228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 6402228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for (sfb=0; sfb<sfbCnt; sfb+=sfbPerGroup) { 6412228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for (sfboffs=0; sfboffs<maxSfbPerGroup; sfboffs++) { 6422228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT sL, sR; 6432228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL inv_n; 6442228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 6452228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project msMask[sfb+sfboffs] = 0; 6462228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (isMask[sfb+sfboffs] == 0) { 6472228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project continue; 6482228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 6492228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 6502228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if ( (sfbEnergyLeft[sfb+sfboffs] < sfbThresholdLeft[sfb+sfboffs]) 6512228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project &&(fMult(FL2FXCONST_DBL(1.0f/1.5f),sfbEnergyRight[sfb+sfboffs]) > sfbThresholdRight[sfb+sfboffs]) ) { 6522228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project continue; 6532228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 6542228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* NEW: if there is a big-enough IS region, switch off PNS */ 6552228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (pnsData[0]) { 6562228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if(pnsData[0]->pnsFlag[sfb+sfboffs]) { 6572228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pnsData[0]->pnsFlag[sfb+sfboffs] = 0; 6582228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 6592228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if(pnsData[1]->pnsFlag[sfb+sfboffs]) { 6602228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pnsData[1]->pnsFlag[sfb+sfboffs] = 0; 6612228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 6622228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 6632228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 6642228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project inv_n = GetInvInt((sfbOffset[sfb + sfboffs + 1] - sfbOffset[sfb + sfboffs])>>1); // scaled with 2 to compensate fMultDiv2() in subsequent loop 6652228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sL = calcSfbMaxScale(mdctSpectrumLeft,sfbOffset[sfb+sfboffs],sfbOffset[sfb+sfboffs+1]); 6662228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sR = calcSfbMaxScale(mdctSpectrumRight,sfbOffset[sfb+sfboffs],sfbOffset[sfb+sfboffs+1]); 6672228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 6682228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project lr = FL2FXCONST_DBL(0.0f); 6692228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for (j=sfbOffset[sfb+sfboffs]; j<sfbOffset[sfb+sfboffs+1]; j++) 6702228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project lr += fMultDiv2(fMultDiv2(mdctSpectrumLeft[j]<<sL,mdctSpectrumRight[j]<<sR),inv_n); 6712228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project lr = lr<<1; 6722228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 6732228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (lr < FL2FXCONST_DBL(0.0f)) { 6742228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* This means OUT OF phase intensity stereo, cf. standard */ 6752228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT s0, s1, s2; 6762228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL tmp, d, ed = FL2FXCONST_DBL(0.0f); 6772228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 6782228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project s0 = fixMin(sL,sR); 6792228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for (j=sfbOffset[sfb+sfboffs]; j<sfbOffset[sfb+sfboffs+1]; j++) { 6802228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project d = ((mdctSpectrumLeft[j]<<s0)>>1) - ((mdctSpectrumRight[j]<<s0)>>1); 6812228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project ed += fMultDiv2(d,d)>>(MDCT_SPEC_SF-1); 6822228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 6832228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project msMask[sfb+sfboffs] = 1; 6842228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project tmp = fDivNorm(sfbEnergyLeft[sfb+sfboffs],ed,&s1); 6852228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project s2 = (s1) + (2*s0) - 2 - MDCT_SPEC_SF; 6862228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (s2 & 1) { 6872228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project tmp = tmp>>1; 6882228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project s2 = s2+1; 6892228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 6902228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project s2 = (s2>>1) + 1; // +1 compensate fMultDiv2() in subsequent loop 6912228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project s2 = fixMin(fixMax(s2,-(DFRACT_BITS-1)),(DFRACT_BITS-1)); 6922228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project scale = sqrtFixp(tmp); 6932228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (s2 < 0) { 6942228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project s2 = -s2; 6952228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for (j=sfbOffset[sfb+sfboffs]; j<sfbOffset[sfb+sfboffs+1]; j++) { 6962228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project mdctSpectrumLeft[j] = (fMultDiv2(mdctSpectrumLeft[j],scale) - fMultDiv2(mdctSpectrumRight[j],scale)) >> s2; 6972228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project mdctSpectrumRight[j] = FL2FXCONST_DBL(0.0f); 6982228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 6992228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 7002228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project else { 7012228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for (j=sfbOffset[sfb+sfboffs]; j<sfbOffset[sfb+sfboffs+1]; j++) { 7022228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project mdctSpectrumLeft[j] = (fMultDiv2(mdctSpectrumLeft[j],scale) - fMultDiv2(mdctSpectrumRight[j],scale)) << s2; 7032228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project mdctSpectrumRight[j] = FL2FXCONST_DBL(0.0f); 7042228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 7052228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 7062228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 7072228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project else { 7082228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* This means IN phase intensity stereo, cf. standard */ 7092228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT s0,s1,s2; 7102228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL tmp, s, es = FL2FXCONST_DBL(0.0f); 7112228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 7122228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project s0 = fixMin(sL,sR); 7132228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for (j=sfbOffset[sfb+sfboffs]; j<sfbOffset[sfb+sfboffs+1]; j++) { 7142228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project s = ((mdctSpectrumLeft[j]<<s0)>>1) + ((mdctSpectrumRight[j]<<s0)>>1); 7152228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project es += fMultDiv2(s,s)>>(MDCT_SPEC_SF-1); // scaled 2*(mdctScale - s0 + 1) + MDCT_SPEC_SF 7162228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 7172228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project msMask[sfb+sfboffs] = 0; 7182228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project tmp = fDivNorm(sfbEnergyLeft[sfb+sfboffs],es,&s1); 7192228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project s2 = (s1) + (2*s0) - 2 - MDCT_SPEC_SF; 7202228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (s2 & 1) { 7212228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project tmp = tmp>>1; 7222228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project s2 = s2 + 1; 7232228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 7242228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project s2 = (s2>>1) + 1; // +1 compensate fMultDiv2() in subsequent loop 7252228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project s2 = fixMin(fixMax(s2,-(DFRACT_BITS-1)),(DFRACT_BITS-1)); 7262228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project scale = sqrtFixp(tmp); 7272228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (s2 < 0) { 7282228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project s2 = -s2; 7292228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for (j=sfbOffset[sfb+sfboffs]; j<sfbOffset[sfb+sfboffs+1]; j++) { 7302228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project mdctSpectrumLeft[j] = (fMultDiv2(mdctSpectrumLeft[j],scale) + fMultDiv2(mdctSpectrumRight[j],scale)) >> s2; 7312228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project mdctSpectrumRight[j] = FL2FXCONST_DBL(0.0f); 7322228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 7332228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 7342228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project else { 7352228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for (j=sfbOffset[sfb+sfboffs]; j<sfbOffset[sfb+sfboffs+1]; j++) { 7362228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project mdctSpectrumLeft[j] = (fMultDiv2(mdctSpectrumLeft[j],scale) + fMultDiv2(mdctSpectrumRight[j],scale)) << s2; 7372228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project mdctSpectrumRight[j] = FL2FXCONST_DBL(0.0f); 7382228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 7392228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 7402228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 7412228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 7422228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project isBook[sfb+sfboffs] = CODE_BOOK_IS_IN_PHASE_NO; 7432228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 7442228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if ( realIsScale[sfb+sfboffs] < FL2FXCONST_DBL(0.0f) ) { 7452228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project isScale[sfb+sfboffs] = (INT)(((realIsScale[sfb+sfboffs]>>1)-FL2FXCONST_DBL(0.5f/(1<<(REAL_SCALE_SF+LD_DATA_SHIFT+1))))>>(DFRACT_BITS-1-REAL_SCALE_SF-LD_DATA_SHIFT-1)) + 1; 7462228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 7472228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project else { 7482228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project isScale[sfb+sfboffs] = (INT)(((realIsScale[sfb+sfboffs]>>1)+FL2FXCONST_DBL(0.5f/(1<<(REAL_SCALE_SF+LD_DATA_SHIFT+1))))>>(DFRACT_BITS-1-REAL_SCALE_SF-LD_DATA_SHIFT-1)); 7492228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 7502228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 7512228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbEnergyRight[sfb+sfboffs] = FL2FXCONST_DBL(0.0f); 7522228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbEnergyLdDataRight[sfb+sfboffs] = FL2FXCONST_DBL(-1.0f); 7532228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbThresholdRight[sfb+sfboffs] = FL2FXCONST_DBL(0.0f); 7542228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbThresholdLdDataRight[sfb+sfboffs] = FL2FXCONST_DBL(-0.515625f); 7552228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sfbSpreadEnRight[sfb+sfboffs] = FL2FXCONST_DBL(0.0f); 7562228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 7572228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *msDigest = MS_SOME; 7582228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 7592228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 7602228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} 7612228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 762