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/*! 852228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project \file 862228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project \brief FDK Fixed Point Arithmetic Library Interface 872228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project*/ 882228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 892228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#ifndef __TRANSCENDENT_H 902228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#define __TRANSCENDENT_H 912228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 922228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#include "sbrdecoder.h" 932228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#include "sbr_rom.h" 942228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 952228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/************************************************************************/ 962228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/*! 972228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project \brief Get number of octaves between frequencies a and b 982228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 992228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project The Result is scaled with 1/8. 1002228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project The valid range for a and b is 1 to LOG_DUALIS_TABLE_SIZE. 1012228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1022228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project \return ld(a/b) / 8 1032228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project*/ 1042228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/************************************************************************/ 1052228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectstatic inline FIXP_SGL FDK_getNumOctavesDiv8(INT a, /*!< lower band */ 1062228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT b) /*!< upper band */ 1072228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 1082228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project return ( (SHORT)((LONG)(CalcLdInt(b) - CalcLdInt(a))>>(FRACT_BITS-3)) ); 1092228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} 1102228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1112228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1122228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/************************************************************************/ 1132228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/*! 1142228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project \brief Add two values given by mantissa and exponent. 1152228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1162228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project Mantissas are in fract format with values between 0 and 1. <br> 1172228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project The base for exponents is 2. Example: \f$ a = a\_m * 2^{a\_e} \f$<br> 1182228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project*/ 1192228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/************************************************************************/ 1202228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectinline void FDK_add_MantExp(FIXP_SGL a_m, /*!< Mantissa of 1st operand a */ 1212228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project SCHAR a_e, /*!< Exponent of 1st operand a */ 1222228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_SGL b_m, /*!< Mantissa of 2nd operand b */ 1232228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project SCHAR b_e, /*!< Exponent of 2nd operand b */ 1242228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_SGL *ptrSum_m, /*!< Mantissa of result */ 1252228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project SCHAR *ptrSum_e) /*!< Exponent of result */ 1262228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 1272228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL accu; 1282228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int shift; 1292228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int shiftAbs; 1302228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1312228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL shiftedMantissa; 1322228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL otherMantissa; 1332228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1342228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* Equalize exponents of the summands. 1352228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project For the smaller summand, the exponent is adapted and 1362228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for compensation, the mantissa is shifted right. */ 1372228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1382228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project shift = (int)(a_e - b_e); 1392228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1402228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project shiftAbs = (shift>0)? shift : -shift; 1412228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project shiftAbs = (shiftAbs < DFRACT_BITS-1)? shiftAbs : DFRACT_BITS-1; 1422228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project shiftedMantissa = (shift>0)? (FX_SGL2FX_DBL(b_m) >> shiftAbs) : (FX_SGL2FX_DBL(a_m) >> shiftAbs); 1432228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project otherMantissa = (shift>0)? FX_SGL2FX_DBL(a_m) : FX_SGL2FX_DBL(b_m); 1442228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *ptrSum_e = (shift>0)? a_e : b_e; 1452228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1462228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu = (shiftedMantissa >> 1) + (otherMantissa >> 1); 1472228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* shift by 1 bit to avoid overflow */ 1482228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1492228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if ( (accu >= (FL2FXCONST_DBL(0.5f) - (FIXP_DBL)1)) || (accu <= FL2FXCONST_DBL(-0.5f)) ) 1502228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *ptrSum_e += 1; 1512228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project else 1522228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu = (shiftedMantissa + otherMantissa); 1532228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1542228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *ptrSum_m = FX_DBL2FX_SGL(accu); 1552228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1562228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} 1572228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1582228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectinline void FDK_add_MantExp(FIXP_DBL a, /*!< Mantissa of 1st operand a */ 1592228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project SCHAR a_e, /*!< Exponent of 1st operand a */ 1602228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL b, /*!< Mantissa of 2nd operand b */ 1612228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project SCHAR b_e, /*!< Exponent of 2nd operand b */ 1622228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *ptrSum, /*!< Mantissa of result */ 1632228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project SCHAR *ptrSum_e) /*!< Exponent of result */ 1642228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 1652228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL accu; 1662228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int shift; 1672228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int shiftAbs; 1682228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1692228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL shiftedMantissa; 1702228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL otherMantissa; 1712228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1722228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* Equalize exponents of the summands. 1732228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project For the smaller summand, the exponent is adapted and 1742228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for compensation, the mantissa is shifted right. */ 1752228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1762228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project shift = (int)(a_e - b_e); 1772228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1782228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project shiftAbs = (shift>0)? shift : -shift; 1792228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project shiftAbs = (shiftAbs < DFRACT_BITS-1)? shiftAbs : DFRACT_BITS-1; 1802228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project shiftedMantissa = (shift>0)? (b >> shiftAbs) : (a >> shiftAbs); 1812228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project otherMantissa = (shift>0)? a : b; 1822228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *ptrSum_e = (shift>0)? a_e : b_e; 1832228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1842228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu = (shiftedMantissa >> 1) + (otherMantissa >> 1); 1852228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* shift by 1 bit to avoid overflow */ 1862228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1872228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if ( (accu >= (FL2FXCONST_DBL(0.5f) - (FIXP_DBL)1)) || (accu <= FL2FXCONST_DBL(-0.5f)) ) 1882228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *ptrSum_e += 1; 1892228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project else 1902228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu = (shiftedMantissa + otherMantissa); 1912228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1922228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *ptrSum = accu; 1932228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1942228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} 1952228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1962228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/************************************************************************/ 1972228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/*! 1982228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project \brief Divide two values given by mantissa and exponent. 1992228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2002228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project Mantissas are in fract format with values between 0 and 1. <br> 2012228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project The base for exponents is 2. Example: \f$ a = a\_m * 2^{a\_e} \f$<br> 2022228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2032228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project For performance reasons, the division is based on a table lookup 2042228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project which limits accuracy. 2052228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project*/ 2062228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/************************************************************************/ 2072228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectstatic inline void FDK_divide_MantExp(FIXP_SGL a_m, /*!< Mantissa of dividend a */ 2082228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project SCHAR a_e, /*!< Exponent of dividend a */ 2092228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_SGL b_m, /*!< Mantissa of divisor b */ 2102228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project SCHAR b_e, /*!< Exponent of divisor b */ 2112228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_SGL *ptrResult_m, /*!< Mantissa of quotient a/b */ 2122228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project SCHAR *ptrResult_e) /*!< Exponent of quotient a/b */ 2132228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2142228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 2152228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int preShift, postShift, index, shift; 2162228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL ratio_m; 2172228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_SGL bInv_m = FL2FXCONST_SGL(0.0f); 2182228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2192228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project preShift = CntLeadingZeros(FX_SGL2FX_DBL(b_m)); 2202228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2212228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* 2222228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project Shift b into the range from 0..INV_TABLE_SIZE-1, 2232228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2242228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project E.g. 10 bits must be skipped for INV_TABLE_BITS 8: 2252228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project - leave 8 bits as index for table 2262228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project - skip sign bit, 2272228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project - skip first bit of mantissa, because this is always the same (>0.5) 2282228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2292228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project We are dealing with energies, so we need not care 2302228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project about negative numbers 2312228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project */ 2322228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2332228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* 2342228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project The first interval has half width so the lowest bit of the index is 2352228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project needed for a doubled resolution. 2362228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project */ 2372228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project shift = (FRACT_BITS - 2 - INV_TABLE_BITS - preShift); 2382228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2392228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project index = (shift<0)? (LONG)b_m << (-shift) : (LONG)b_m >> shift; 2402228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2412228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2422228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* The index has INV_TABLE_BITS +1 valid bits here. Clear the other bits. */ 2432228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project index &= (1 << (INV_TABLE_BITS+1)) - 1; 2442228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2452228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* Remove offset of half an interval */ 2462228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project index--; 2472228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2482228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* Now the lowest bit is shifted out */ 2492228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project index = index >> 1; 2502228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2512228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* Fetch inversed mantissa from table: */ 2522228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project bInv_m = (index<0)? bInv_m : FDK_sbrDecoder_invTable[index]; 2532228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2542228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* Multiply a with the inverse of b: */ 2552228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project ratio_m = (index<0)? FX_SGL2FX_DBL(a_m >> 1) : fMultDiv2(bInv_m,a_m); 2562228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2572228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project postShift = CntLeadingZeros(ratio_m)-1; 2582228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2592228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *ptrResult_m = FX_DBL2FX_SGL(ratio_m << postShift); 2602228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *ptrResult_e = a_e - b_e + 1 + preShift - postShift; 2612228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} 2622228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2632228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectstatic inline void FDK_divide_MantExp(FIXP_DBL a_m, /*!< Mantissa of dividend a */ 2642228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project SCHAR a_e, /*!< Exponent of dividend a */ 2652228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL b_m, /*!< Mantissa of divisor b */ 2662228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project SCHAR b_e, /*!< Exponent of divisor b */ 2672228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *ptrResult_m, /*!< Mantissa of quotient a/b */ 2682228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project SCHAR *ptrResult_e) /*!< Exponent of quotient a/b */ 2692228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2702228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 2712228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int preShift, postShift, index, shift; 2722228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL ratio_m; 2732228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_SGL bInv_m = FL2FXCONST_SGL(0.0f); 2742228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2752228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project preShift = CntLeadingZeros(b_m); 2762228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2772228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* 2782228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project Shift b into the range from 0..INV_TABLE_SIZE-1, 2792228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2802228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project E.g. 10 bits must be skipped for INV_TABLE_BITS 8: 2812228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project - leave 8 bits as index for table 2822228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project - skip sign bit, 2832228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project - skip first bit of mantissa, because this is always the same (>0.5) 2842228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2852228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project We are dealing with energies, so we need not care 2862228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project about negative numbers 2872228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project */ 2882228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2892228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* 2902228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project The first interval has half width so the lowest bit of the index is 2912228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project needed for a doubled resolution. 2922228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project */ 2932228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project shift = (DFRACT_BITS - 2 - INV_TABLE_BITS - preShift); 2942228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2952228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project index = (shift<0)? (LONG)b_m << (-shift) : (LONG)b_m >> shift; 2962228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2972228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2982228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* The index has INV_TABLE_BITS +1 valid bits here. Clear the other bits. */ 2992228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project index &= (1 << (INV_TABLE_BITS+1)) - 1; 3002228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3012228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* Remove offset of half an interval */ 3022228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project index--; 3032228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3042228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* Now the lowest bit is shifted out */ 3052228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project index = index >> 1; 3062228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3072228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* Fetch inversed mantissa from table: */ 3082228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project bInv_m = (index<0)? bInv_m : FDK_sbrDecoder_invTable[index]; 3092228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3102228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* Multiply a with the inverse of b: */ 3112228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project ratio_m = (index<0)? (a_m >> 1) : fMultDiv2(bInv_m,a_m); 3122228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3132228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project postShift = CntLeadingZeros(ratio_m)-1; 3142228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3152228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *ptrResult_m = ratio_m << postShift; 3162228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *ptrResult_e = a_e - b_e + 1 + preShift - postShift; 3172228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} 3182228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3192228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/*! 3202228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project \brief Calculate the squareroot of a number given by mantissa and exponent 3212228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3222228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project Mantissa is in fract format with values between 0 and 1. <br> 3232228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project The base for the exponent is 2. Example: \f$ a = a\_m * 2^{a\_e} \f$<br> 3242228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project The operand is addressed via pointers and will be overwritten with the result. 3252228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3262228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project For performance reasons, the square root is based on a table lookup 3272228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project which limits accuracy. 3282228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project*/ 3292228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectstatic inline void FDK_sqrt_MantExp(FIXP_DBL *mantissa, /*!< Pointer to mantissa */ 3302228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project SCHAR *exponent, 3312228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const SCHAR *destScale) 3322228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 3332228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL input_m = *mantissa; 3342228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int input_e = (int) *exponent; 3352228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL result = FL2FXCONST_DBL(0.0f); 3362228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int result_e = -FRACT_BITS; 3372228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3382228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* Call lookup square root, which does internally normalization. */ 3392228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project result = sqrtFixp_lookup(input_m, &input_e); 3402228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project result_e = input_e; 3412228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3422228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* Write result */ 3432228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (exponent==destScale) { 3442228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *mantissa = result; 3452228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *exponent = result_e; 3462228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } else { 3472228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int shift = result_e - *destScale; 3482228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *mantissa = (shift>=0) ? result << (INT)fixMin(DFRACT_BITS-1,shift) 3492228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project : result >> (INT)fixMin(DFRACT_BITS-1,-shift); 3502228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *exponent = *destScale; 3512228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 3522228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} 3532228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3542228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3552228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#endif 356