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 5e1c78ed73faa51f2c7fcb0c4a17a92be9cc747f5Jean-Michel Trivi� Copyright 1995 - 2015 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/*************************** Fraunhofer IIS FDK Tools ********************** 852228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 862228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project Author(s): M. Gayer 872228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project Description: Fixed point specific mathematical functions 882228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 892228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project******************************************************************************/ 902228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 912228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#ifndef __fixpoint_math_H 922228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#define __fixpoint_math_H 932228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 942228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 952228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#include "common_fix.h" 962228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 9746ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi#if !defined(FUNCTION_fIsLessThan) 9846ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi/** 9946ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi * \brief Compares two fixpoint values incl. scaling. 10046ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi * \param a_m mantissa of the first input value. 10146ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi * \param a_e exponent of the first input value. 10246ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi * \param b_m mantissa of the second input value. 10346ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi * \param b_e exponent of the second input value. 10446ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi * \return non-zero if (a_m*2^a_e) < (b_m*2^b_e), 0 otherwise 10546ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi */ 10646ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel TriviFDK_INLINE INT fIsLessThan(FIXP_DBL a_m, INT a_e, FIXP_DBL b_m, INT b_e) 10746ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi{ 10846ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi if (a_e > b_e) { 10946ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi return (b_m >> fMin(a_e-b_e, DFRACT_BITS-1) > a_m); 11046ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi } else { 11146ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi return (a_m >> fMin(b_e-a_e, DFRACT_BITS-1) < b_m); 11246ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi } 11346ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi} 11446ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi 11546ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel TriviFDK_INLINE INT fIsLessThan(FIXP_SGL a_m, INT a_e, FIXP_SGL b_m, INT b_e) 11646ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi{ 11746ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi if (a_e > b_e) { 11846ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi return (b_m >> fMin(a_e-b_e, FRACT_BITS-1) > a_m); 11946ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi } else { 12046ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi return (a_m >> fMin(b_e-a_e, FRACT_BITS-1) < b_m); 12146ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi } 12246ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi} 12346ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi#endif 12446ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi 12546ba3676b854acbc69a4c7845f578d4c2886377bJean-Michel Trivi 1262228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1272228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#define LD_DATA_SCALING (64.0f) 1282228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#define LD_DATA_SHIFT 6 /* pow(2, LD_DATA_SHIFT) = LD_DATA_SCALING */ 1292228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1302228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/** 1312228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \brief deprecated. Use fLog2() instead. 1322228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project */ 1332228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectFIXP_DBL CalcLdData(FIXP_DBL op); 1342228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1352228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectvoid LdDataVector(FIXP_DBL *srcVector, FIXP_DBL *destVector, INT number); 1362228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1372228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectFIXP_DBL CalcInvLdData(FIXP_DBL op); 1382228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1392228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1402228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectvoid InitLdInt(); 1412228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectFIXP_DBL CalcLdInt(INT i); 1422228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1432228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectextern const USHORT sqrt_tab[49]; 1442228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1452228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectinline FIXP_DBL sqrtFixp_lookup(FIXP_DBL x) 1462228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 1472228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project UINT y = (INT)x; 1482228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project UCHAR is_zero=(y==0); 1492228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT zeros=fixnormz_D(y) & 0x1e; 1502228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project y<<=zeros; 1512228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project UINT idx=(y>>26)-16; 1522228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project USHORT frac=(y>>10)&0xffff; 1532228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project USHORT nfrac=0xffff^frac; 1542228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project UINT t=nfrac*sqrt_tab[idx]+frac*sqrt_tab[idx+1]; 1552228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project t=t>>(zeros>>1); 1562228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project return(is_zero ? 0 : t); 1572228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} 1582228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1592228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectinline FIXP_DBL sqrtFixp_lookup(FIXP_DBL x, INT *x_e) 1602228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 1612228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project UINT y = (INT)x; 1622228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT e; 1632228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1642228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (x == (FIXP_DBL)0) { 1652228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project return x; 1662228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 1672228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1682228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* Normalize */ 1692228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project e=fixnormz_D(y); 1702228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project y<<=e; 1712228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project e = *x_e - e + 2; 1722228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1732228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* Correct odd exponent. */ 1742228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (e & 1) { 1752228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project y >>= 1; 1762228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project e ++; 1772228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 1782228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* Get square root */ 1792228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project UINT idx=(y>>26)-16; 1802228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project USHORT frac=(y>>10)&0xffff; 1812228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project USHORT nfrac=0xffff^frac; 1822228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project UINT t=nfrac*sqrt_tab[idx]+frac*sqrt_tab[idx+1]; 1832228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1842228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* Write back exponent */ 1852228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *x_e = e >> 1; 1862228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project return (FIXP_DBL)(LONG)(t>>1); 1872228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} 1882228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1892228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1902228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1912228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectFIXP_DBL sqrtFixp(FIXP_DBL op); 1922228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1932228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectvoid InitInvSqrtTab(); 1942228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1952228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectFIXP_DBL invSqrtNorm2(FIXP_DBL op, INT *shift); 1962228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1972228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/***************************************************************************** 1982228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1992228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project functionname: invFixp 2002228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project description: delivers 1/(op) 2012228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2022228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project*****************************************************************************/ 2032228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectinline FIXP_DBL invFixp(FIXP_DBL op) 2042228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 2052228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT tmp_exp ; 2062228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL tmp_inv = invSqrtNorm2(op, &tmp_exp) ; 2072228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDK_ASSERT((31-(2*tmp_exp+1))>=0) ; 2082228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project return ( fPow2Div2( (FIXP_DBL)tmp_inv ) >> (31-(2*tmp_exp+1)) ) ; 2092228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} 2102228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2112228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2122228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2132228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#if defined(__mips__) && (__GNUC__==2) 2142228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2152228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#define FUNCTION_schur_div 2162228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectinline FIXP_DBL schur_div(FIXP_DBL num,FIXP_DBL denum, INT count) 2172228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 2182228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT result, tmp ; 2192228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project __asm__ ("srl %1, %2, 15\n" 2202228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project "div %3, %1\n" : "=lo" (result) 2212228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project : "%d" (tmp), "d" (denum) , "d" (num) 2222228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project : "hi" ) ; 2232228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project return result<<16 ; 2242228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} 2252228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2262228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/*###########################################################################################*/ 2272228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#elif defined(__mips__) && (__GNUC__==3) 2282228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2292228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#define FUNCTION_schur_div 2302228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectinline FIXP_DBL schur_div(FIXP_DBL num,FIXP_DBL denum, INT count) 2312228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 2322228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT result, tmp; 2332228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2342228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project __asm__ ("srl %[tmp], %[denum], 15\n" 2352228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project "div %[result], %[num], %[tmp]\n" 2362228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project : [tmp] "+r" (tmp), [result]"=r"(result) 2372228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project : [denum]"r"(denum), [num]"r"(num) 2382228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project : "hi", "lo"); 2392228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project return result << (DFRACT_BITS-16); 2402228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} 2412228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2422228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/*###########################################################################################*/ 2432228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#elif defined(SIMULATE_MIPS_DIV) 2442228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2452228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#define FUNCTION_schur_div 2462228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectinline FIXP_DBL schur_div(FIXP_DBL num, FIXP_DBL denum, INT count) 2472228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 2482228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDK_ASSERT (count<=DFRACT_BITS-1); 2492228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDK_ASSERT (num>=(FIXP_DBL)0); 2502228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDK_ASSERT (denum>(FIXP_DBL)0); 2512228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDK_ASSERT (num <= denum); 2522228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2532228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT tmp = denum >> (count-1); 2542228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT result = 0; 2552228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2562228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project while (num > tmp) 2572228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project { 2582228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project num -= tmp; 2592228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project result++; 2602228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 2612228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2622228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project return result << (DFRACT_BITS-count); 2632228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} 2642228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2652228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/*###########################################################################################*/ 2662228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#endif /* target architecture selector */ 2672228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2682228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#if !defined(FUNCTION_schur_div) 2692228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/** 2702228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \brief Divide two FIXP_DBL values with given precision. 2712228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param num dividend 2722228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param denum divisor 2732228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param count amount of significant bits of the result (starting to the MSB) 2742228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \return num/divisor 2752228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project */ 2762228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectFIXP_DBL schur_div(FIXP_DBL num,FIXP_DBL denum, INT count); 2772228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#endif 2782228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2792228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2802228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2812228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectFIXP_DBL mul_dbl_sgl_rnd (const FIXP_DBL op1, 2822228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const FIXP_SGL op2); 2832228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2842228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/** 2852228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \brief multiply two values with normalization, thus max precision. 2862228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * Author: Robert Weidner 2872228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * 2882228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param f1 first factor 2892228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param f2 secod factor 2902228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param result_e pointer to an INT where the exponent of the result is stored into 2912228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \return mantissa of the product f1*f2 2922228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project */ 2932228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectFIXP_DBL fMultNorm( 2942228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL f1, 2952228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL f2, 2962228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT *result_e 2972228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project ); 2982228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2992228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectinline FIXP_DBL fMultNorm(FIXP_DBL f1, FIXP_DBL f2) 3002228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 3012228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL m; 3022228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT e; 3032228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3042228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project m = fMultNorm(f1, f2, &e); 3052228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3062228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project m = scaleValueSaturate(m, e); 3072228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3082228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project return m; 3092228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} 3102228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3112228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/** 3122228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \brief Divide 2 FIXP_DBL values with normalization of input values. 3132228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param num numerator 3142228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param denum denomintator 3152228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \return num/denum with exponent = 0 3162228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project */ 3172228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectFIXP_DBL fDivNorm(FIXP_DBL num, FIXP_DBL denom, INT *result_e); 3182228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3192228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/** 3202228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \brief Divide 2 FIXP_DBL values with normalization of input values. 3212228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param num numerator 3222228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param denum denomintator 3232228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param result_e pointer to an INT where the exponent of the result is stored into 3242228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \return num/denum with exponent = *result_e 3252228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project */ 3262228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectFIXP_DBL fDivNorm(FIXP_DBL num, FIXP_DBL denom); 3272228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3282228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/** 3292228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \brief Divide 2 FIXP_DBL values with normalization of input values. 3302228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param num numerator 3312228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param denum denomintator 3322228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \return num/denum with exponent = 0 3332228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project */ 3342228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectFIXP_DBL fDivNormHighPrec(FIXP_DBL L_num, FIXP_DBL L_denum, INT *result_e); 3352228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3362228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/** 3372228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \brief Calculate log(argument)/log(2) (logarithm with base 2). deprecated. Use fLog2() instead. 3382228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param arg mantissa of the argument 3392228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param arg_e exponent of the argument 3402228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param result_e pointer to an INT to store the exponent of the result 3412228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \return the mantissa of the result. 3422228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param 3432228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project */ 3442228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectFIXP_DBL CalcLog2(FIXP_DBL arg, INT arg_e, INT *result_e); 3452228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3462228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/** 3472228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \brief return 2 ^ (exp * 2^exp_e) 3482228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param exp_m mantissa of the exponent to 2.0f 3492228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param exp_e exponent of the exponent to 2.0f 3502228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param result_e pointer to a INT where the exponent of the result will be stored into 3512228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \return mantissa of the result 3522228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project */ 3532228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectFIXP_DBL f2Pow(const FIXP_DBL exp_m, const INT exp_e, INT *result_e); 3542228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3552228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/** 3562228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \brief return 2 ^ (exp_m * 2^exp_e). This version returns only the mantissa with implicit exponent of zero. 3572228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param exp_m mantissa of the exponent to 2.0f 3582228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param exp_e exponent of the exponent to 2.0f 3592228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \return mantissa of the result 3602228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project */ 3612228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectFIXP_DBL f2Pow(const FIXP_DBL exp_m, const INT exp_e); 3622228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3632228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/** 3642228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \brief return x ^ (exp * 2^exp_e), where log2(x) = baseLd_m * 2^(baseLd_e). This saves 3652228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * the need to compute log2() of constant values (when x is a constant). 3662228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param ldx_m mantissa of log2() of x. 3672228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param ldx_e exponent of log2() of x. 3682228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param exp_m mantissa of the exponent to 2.0f 3692228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param exp_e exponent of the exponent to 2.0f 3702228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param result_e pointer to a INT where the exponent of the result will be stored into 3712228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \return mantissa of the result 3722228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project */ 3732228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectFIXP_DBL fLdPow( 3742228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL baseLd_m, 3752228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT baseLd_e, 3762228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL exp_m, INT exp_e, 3772228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project INT *result_e 3782228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project ); 3792228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3802228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/** 3812228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \brief return x ^ (exp * 2^exp_e), where log2(x) = baseLd_m * 2^(baseLd_e). This saves 3822228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * the need to compute log2() of constant values (when x is a constant). This version 3832228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * does not return an exponent, which is implicitly 0. 3842228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param ldx_m mantissa of log2() of x. 3852228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param ldx_e exponent of log2() of x. 3862228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param exp_m mantissa of the exponent to 2.0f 3872228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param exp_e exponent of the exponent to 2.0f 3882228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \return mantissa of the result 3892228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project */ 3902228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectFIXP_DBL fLdPow( 3912228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL baseLd_m, INT baseLd_e, 3922228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL exp_m, INT exp_e 3932228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project ); 3942228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3952228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/** 3962228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \brief return (base * 2^base_e) ^ (exp * 2^exp_e). Use fLdPow() instead whenever possible. 3972228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param base_m mantissa of the base. 3982228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param base_e exponent of the base. 3992228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param exp_m mantissa of power to be calculated of the base. 4002228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param exp_e exponent of power to be calculated of the base. 4012228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param result_e pointer to a INT where the exponent of the result will be stored into. 4022228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \return mantissa of the result. 4032228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project */ 4042228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectFIXP_DBL fPow(FIXP_DBL base_m, INT base_e, FIXP_DBL exp_m, INT exp_e, INT *result_e); 4052228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4062228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/** 4072228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \brief return (base * 2^base_e) ^ N 4082228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param base mantissa of the base 4092228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param base_e exponent of the base 4102228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param power to be calculated of the base 4112228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param result_e pointer to a INT where the exponent of the result will be stored into 4122228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \return mantissa of the result 4132228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project */ 4142228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectFIXP_DBL fPowInt(FIXP_DBL base_m, INT base_e, INT N, INT *result_e); 4152228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4162228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/** 4172228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \brief calculate logarithm of base 2 of x_m * 2^(x_e) 4182228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param x_m mantissa of the input value. 4192228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param x_e exponent of the input value. 4202228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param pointer to an INT where the exponent of the result is returned into. 4212228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \return mantissa of the result. 4222228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project */ 4232228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectFIXP_DBL fLog2(FIXP_DBL x_m, INT x_e, INT *result_e); 4242228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4252228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/** 4262228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \brief calculate logarithm of base 2 of x_m * 2^(x_e) 4272228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param x_m mantissa of the input value. 4282228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param x_e exponent of the input value. 4292228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \return mantissa of the result with implicit exponent of LD_DATA_SHIFT. 4302228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project */ 4312228e360595641dd906bf1773307f43d304f5b2The Android Open Source ProjectFIXP_DBL fLog2(FIXP_DBL x_m, INT x_e); 4322228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4332228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/** 4342228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \brief Add with saturation of the result. 4352228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param a first summand 4362228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param b second summand 4372228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \return saturated sum of a and b. 4382228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project */ 4392228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectinline FIXP_SGL fAddSaturate(const FIXP_SGL a, const FIXP_SGL b) 4402228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 4412228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project LONG sum; 4422228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4432228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sum = (LONG)(SHORT)a + (LONG)(SHORT)b; 4442228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sum = fMax(fMin((INT)sum, (INT)MAXVAL_SGL), (INT)MINVAL_SGL); 4452228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project return (FIXP_SGL)(SHORT)sum; 4462228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} 4472228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4482228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/** 4492228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \brief Add with saturation of the result. 4502228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param a first summand 4512228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param b second summand 4522228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \return saturated sum of a and b. 4532228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project */ 4542228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectinline FIXP_DBL fAddSaturate(const FIXP_DBL a, const FIXP_DBL b) 4552228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 4562228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project LONG sum; 4572228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4582228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sum = (LONG)(a>>1) + (LONG)(b>>1); 4592228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project sum = fMax(fMin((INT)sum, (INT)(MAXVAL_DBL>>1)), (INT)(MINVAL_DBL>>1)); 4602228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project return (FIXP_DBL)(LONG)(sum<<1); 4612228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} 4622228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4632228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project//#define TEST_ROUNDING 4642228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4652228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4662228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4672228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4682228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/***************************************************************************** 4692228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 470e1c78ed73faa51f2c7fcb0c4a17a92be9cc747f5Jean-Michel Trivi array for 1/n, n=1..80 4712228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4722228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project****************************************************************************/ 4732228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 474e1c78ed73faa51f2c7fcb0c4a17a92be9cc747f5Jean-Michel Trivi extern const FIXP_DBL invCount[80]; 4752228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4762228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project LNK_SECTION_INITCODE 4772228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project inline void InitInvInt(void) {} 4782228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4792228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4802228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project/** 4812228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \brief Calculate the value of 1/i where i is a integer value. It supports 48229189063770fbd7d00f04ed1fd16ec5eefee7ec9Jean-Michel Trivi * input values from 0 upto 79. 4832228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param intValue Integer input value. 4842228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project * \param FIXP_DBL representation of 1/intValue 4852228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project */ 4862228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectinline FIXP_DBL GetInvInt(int intValue) 4872228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 48829189063770fbd7d00f04ed1fd16ec5eefee7ec9Jean-Michel Trivi FDK_ASSERT((intValue >= 0) && (intValue < 80)); 48929189063770fbd7d00f04ed1fd16ec5eefee7ec9Jean-Michel Trivi if (intValue > 79) 49029189063770fbd7d00f04ed1fd16ec5eefee7ec9Jean-Michel Trivi return invCount[79]; 49129189063770fbd7d00f04ed1fd16ec5eefee7ec9Jean-Michel Trivi else if (intValue < 0) 49229189063770fbd7d00f04ed1fd16ec5eefee7ec9Jean-Michel Trivi return invCount[0]; 49329189063770fbd7d00f04ed1fd16ec5eefee7ec9Jean-Michel Trivi else 49429189063770fbd7d00f04ed1fd16ec5eefee7ec9Jean-Michel Trivi return invCount[intValue]; 4952228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} 4962228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4972228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4982228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#endif 4992228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 500