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 dct.cpp 862228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project \brief DCT Implementations 872228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project Library functions to calculate standard DCTs. This will most likely be replaced by hand-optimized 882228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project functions for the specific target processor. 892228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 902228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project Three different implementations of the dct type II and the dct type III transforms are provided. 912228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 922228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project By default implementations which are based on a single, standard complex FFT-kernel are used (dctII_f() and dctIII_f()). 932228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project These are specifically helpful in cases where optimized FFT libraries are already available. The FFT used in these 942228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project implementation is FFT rad2 from FDK_tools. 952228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 962228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project Of course, one might also use DCT-libraries should they be available. The DCT and DST 972228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project type IV implementations are only available in a version based on a complex FFT kernel. 982228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project*/ 992228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1002228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#include "dct.h" 1012228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1022228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1032228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#include "FDK_tools_rom.h" 1042228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#include "fft.h" 1052228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1062228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1072228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#if defined(__arm__) 1082228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#include "arm/dct_arm.cpp" 1092228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#endif 1102228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1112228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1122228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#if !defined(FUNCTION_dct_III) 1132228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectvoid dct_III(FIXP_DBL *pDat, /*!< pointer to input/output */ 1142228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *tmp, /*!< pointer to temporal working buffer */ 1152228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int L, /*!< lenght of transform */ 1162228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int *pDat_e 1172228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project ) 1182228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 1192228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDK_ASSERT(L == 64 || L == 32); 1202228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int i; 1212228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL xr, accu1, accu2; 1222228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int inc; 1232228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int M = L>>1; 1242228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int ld_M; 1252228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1262228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (L == 64) ld_M = 5; 1272228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project else ld_M = 4; 1282228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1292228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* This loop performs multiplication for index i (i*inc) */ 1302228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project inc = (64/2) >> ld_M; /* 64/L */ 1312228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1322228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *pTmp_0 = &tmp[2]; 1332228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *pTmp_1 = &tmp[(M-1)*2]; 1342228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1352228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for(i=1; i<M>>1; i++,pTmp_0+=2,pTmp_1-=2) { 1362228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1372228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL accu3,accu4,accu5,accu6; 1382228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1392228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project cplxMultDiv2(&accu2, &accu1, pDat[L - i], pDat[i], sin_twiddle_L64[i*inc]); 1402228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project cplxMultDiv2(&accu4, &accu3, pDat[M+i], pDat[M-i], sin_twiddle_L64[(M-i)*inc]); 1412228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu3 >>= 1; accu4 >>= 1; 1422228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1432228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* This method is better for ARM926, that uses operand2 shifted right by 1 always */ 1442228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project cplxMultDiv2(&accu6, &accu5, (accu3 - (accu1>>1)), ((accu2>>1) + accu4), sin_twiddle_L64[(4*i)*inc]); 1452228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project xr = (accu1>>1) + accu3; 1462228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pTmp_0[0] = (xr>>1) - accu5; 1472228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pTmp_1[0] = (xr>>1) + accu5; 1482228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1492228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project xr = (accu2>>1) - accu4; 1502228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pTmp_0[1] = (xr>>1) - accu6; 1512228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pTmp_1[1] = -((xr>>1) + accu6); 1522228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1532228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 1542228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1552228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project xr = fMultDiv2(pDat[M], sin_twiddle_L64[64/2].v.re );/* cos((PI/(2*L))*M); */ 1562228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project tmp[0] = ((pDat[0]>>1) + xr)>>1; 1572228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project tmp[1] = ((pDat[0]>>1) - xr)>>1; 1582228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1592228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project cplxMultDiv2(&accu2, &accu1, pDat[L - (M/2)], pDat[M/2], sin_twiddle_L64[64/4]); 1602228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project tmp[M] = accu1>>1; 1612228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project tmp[M+1] = accu2>>1; 1622228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1632228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* dit_fft expects 1 bit scaled input values */ 1642228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project fft(M, tmp, pDat_e); 1652228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1662228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* ARM926: 12 cycles per 2-iteration, no overhead code by compiler */ 1672228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pTmp_1 = &tmp[L]; 1682228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for (i = M>>1; i--;) 1692228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project { 1702228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL tmp1, tmp2, tmp3, tmp4; 1712228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project tmp1 = *tmp++; 1722228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project tmp2 = *tmp++; 1732228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project tmp3 = *--pTmp_1; 1742228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project tmp4 = *--pTmp_1; 1752228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *pDat++ = tmp1; 1762228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *pDat++ = tmp3; 1772228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *pDat++ = tmp2; 1782228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *pDat++ = tmp4; 1792228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 1802228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1812228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *pDat_e += 2; 1822228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} 1832228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#endif 1842228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1852228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#if !defined(FUNCTION_dct_II) 1862228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectvoid dct_II(FIXP_DBL *pDat, /*!< pointer to input/output */ 1872228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *tmp, /*!< pointer to temporal working buffer */ 1882228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int L, /*!< lenght of transform */ 1892228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int *pDat_e 1902228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project ) 1912228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 1922228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDK_ASSERT(L == 64 || L == 32); 1932228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL accu1,accu2; 1942228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *pTmp_0, *pTmp_1; 1952228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 1962228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int i; 1972228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int inc; 1982228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int M = L>>1; 1992228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int ld_M; 2002228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2012228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDK_ASSERT(L == 64 || L == 32); 2022228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project ld_M = 4 + (L >> 6); /* L=64: 5, L=32: 4 */ 2032228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2042228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project inc = (64/2) >> ld_M; /* L=64: 1, L=32: 2 */ 2052228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2062228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *pdat = &pDat[0]; 2072228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL accu3, accu4; 2082228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pTmp_0 = &tmp[0]; 2092228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pTmp_1 = &tmp[L-1]; 2102228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for (i = M>>1; i--; ) 2112228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project { 2122228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu1 = *pdat++; 2132228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu2 = *pdat++; 2142228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu3 = *pdat++; 2152228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu4 = *pdat++; 2162228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu1 >>= 1; 2172228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu2 >>= 1; 2182228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu3 >>= 1; 2192228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu4 >>= 1; 2202228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *pTmp_0++ = accu1; 2212228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *pTmp_0++ = accu3; 2222228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *pTmp_1-- = accu2; 2232228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *pTmp_1-- = accu4; 2242228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 2252228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2262228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2272228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project fft(M, tmp, pDat_e); 2282228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2292228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pTmp_0 = &tmp[2]; 2302228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pTmp_1 = &tmp[(M-1)*2]; 2312228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2322228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for (i=1; i<M>>1; i++,pTmp_0+=2,pTmp_1-=2) { 2332228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2342228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL a1,a2; 2352228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL accu3, accu4; 2362228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2372228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project a1 = ((pTmp_0[1]>>1) + (pTmp_1[1]>>1)); 2382228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project a2 = ((pTmp_1[0]>>1) - (pTmp_0[0]>>1)); 2392228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2402228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project cplxMultDiv2(&accu1, &accu2, a2, a1, sin_twiddle_L64[(4*i)*inc]); 2412228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu1<<=1; accu2<<=1; 2422228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2432228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project a1 = ((pTmp_0[0]>>1) + (pTmp_1[0]>>1)); 2442228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project a2 = ((pTmp_0[1]>>1) - (pTmp_1[1]>>1)); 2452228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2462228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project cplxMultDiv2(&accu3, &accu4, (a1 + accu2), -(accu1 + a2), sin_twiddle_L64[i*inc]); 2472228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat[L - i] = accu4; 2482228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat[i] = accu3; 2492228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2502228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project cplxMultDiv2(&accu3, &accu4, (a1 - accu2), -(accu1 - a2), sin_twiddle_L64[(M-i)*inc]); 2512228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat[M + i] = accu4; 2522228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat[M - i] = accu3; 2532228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2542228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 2552228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2562228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project cplxMultDiv2(&accu1, &accu2, tmp[M], tmp[M+1], sin_twiddle_L64[(M/2)*inc]); 2572228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat[L - (M/2)] = accu2; 2582228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat[M/2] = accu1; 2592228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2602228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat[0] = (tmp[0]>>1)+(tmp[1]>>1); 2612228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat[M] = fMult(((tmp[0]>>1)-(tmp[1]>>1)), sin_twiddle_L64[64/2].v.re);/* cos((PI/(2*L))*M); */ 2622228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2632228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *pDat_e += 2; 2642228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} 2652228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#endif 2662228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2672228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectstatic 2682228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectvoid getTables(const FIXP_WTP **twiddle, const FIXP_STP **sin_twiddle, int *sin_step, int length) 2692228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 2702228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int ld2_length; 2712228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2722228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* Get ld2 of length - 2 + 1 2732228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project -2: because first table entry is window of size 4 2742228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project +1: because we already include +1 because of ceil(log2(length)) */ 2752228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project ld2_length = DFRACT_BITS-1-fNormz((FIXP_DBL)length) - 1; 2762228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2772228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* Extract sort of "eigenvalue" (the 4 left most bits) of length. */ 2782228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project switch ( (length) >> (ld2_length-1) ) { 2792228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project case 0x4: /* radix 2 */ 2802228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *sin_twiddle = SineTable512; 2812228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *sin_step = 1<<(9 - ld2_length); 2822228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *twiddle = windowSlopes[0][0][ld2_length-1]; 2832228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project break; 2842228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project case 0x7: /* 10 ms */ 2852228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *sin_twiddle = SineTable480; 2862228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *sin_step = 1<<(8 - ld2_length); 2872228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *twiddle = windowSlopes[0][1][ld2_length]; 2882228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project break; 2892228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project default: 2902228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *sin_twiddle = NULL; 2912228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *sin_step = 0; 2922228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *twiddle = NULL; 2932228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project break; 2942228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 2952228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2962228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDK_ASSERT(*twiddle != NULL); 2972228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 2982228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDK_ASSERT(*sin_step > 0); 2992228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3002228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} 3012228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3022228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#if !defined(FUNCTION_dct_IV) 3032228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3042228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectvoid dct_IV(FIXP_DBL *pDat, 3052228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int L, 3062228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int *pDat_e) 3072228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 3082228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int sin_step = 0; 3092228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int M = L >> 1; 3102228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3112228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const FIXP_WTP *twiddle; 3122228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const FIXP_STP *sin_twiddle; 3132228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3142228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDK_ASSERT(L >= 4); 3152228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3162228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project getTables(&twiddle, &sin_twiddle, &sin_step, L); 3172228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3182228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#ifdef FUNCTION_dct_IV_func1 3192228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (M>=4 && (M&3) == 0) { 3202228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* ARM926: 44 cycles for 2 iterations = 22 cycles/iteration */ 3212228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project dct_IV_func1(M>>2, twiddle, &pDat[0], &pDat[L-1]); 3222228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } else 3232228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#endif /* FUNCTION_dct_IV_func1 */ 3242228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project { 3252228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *RESTRICT pDat_0 = &pDat[0]; 3262228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *RESTRICT pDat_1 = &pDat[L - 2]; 3272228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project register int i; 3282228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3292228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* 29 cycles on ARM926 */ 3302228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for (i = 0; i < M-1; i+=2,pDat_0+=2,pDat_1-=2) 3312228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project { 3322228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project register FIXP_DBL accu1,accu2,accu3,accu4; 3332228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3342228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu1 = pDat_1[1]; accu2 = pDat_0[0]; 3352228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu3 = pDat_0[1]; accu4 = pDat_1[0]; 3362228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3372228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project cplxMultDiv2(&accu1, &accu2, accu1, accu2, twiddle[i]); 3382228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project cplxMultDiv2(&accu3, &accu4, accu4, accu3, twiddle[i+1]); 3392228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3402228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_0[0] = accu2; pDat_0[1] = accu1; 3412228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_1[0] = accu4; pDat_1[1] = -accu3; 3422228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 3432228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (M&1) 3442228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project { 3452228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project register FIXP_DBL accu1,accu2; 3462228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3472228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu1 = pDat_1[1]; accu2 = pDat_0[0]; 3482228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3492228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project cplxMultDiv2(&accu1, &accu2, accu1, accu2, twiddle[i]); 3502228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3512228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_0[0] = accu2; pDat_0[1] = accu1; 3522228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 3532228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 3542228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3552228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project fft(M, pDat, pDat_e); 3562228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3572228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#ifdef FUNCTION_dct_IV_func2 3582228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (M>=4 && (M&3) == 0) { 3592228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* ARM926: 42 cycles for 2 iterations = 21 cycles/iteration */ 3602228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project dct_IV_func2(M>>2, sin_twiddle, &pDat[0], &pDat[L], sin_step); 3612228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } else 3622228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#endif /* FUNCTION_dct_IV_func2 */ 3632228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project { 3642228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *RESTRICT pDat_0 = &pDat[0]; 3652228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *RESTRICT pDat_1 = &pDat[L - 2]; 3662228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project register FIXP_DBL accu1,accu2,accu3,accu4; 3672228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int idx, i; 3682228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3692228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* Sin and Cos values are 0.0f and 1.0f */ 3702228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu1 = pDat_1[0]; 3712228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu2 = pDat_1[1]; 3722228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3732228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_1[1] = -(pDat_0[1]>>1); 3742228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_0[0] = (pDat_0[0]>>1); 3752228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3762228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3772228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* 28 cycles for ARM926 */ 3782228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for (idx = sin_step,i=1; i<(M+1)>>1; i++, idx+=sin_step) 3792228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project { 3802228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_STP twd = sin_twiddle[idx]; 3812228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project cplxMultDiv2(&accu3, &accu4, accu1, accu2, twd); 3822228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_0[1] = accu3; 3832228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_1[0] = accu4; 3842228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3852228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_0+=2; 3862228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_1-=2; 3872228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3882228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project cplxMultDiv2(&accu3, &accu4, pDat_0[1], pDat_0[0], twd); 3892228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3902228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu1 = pDat_1[0]; 3912228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu2 = pDat_1[1]; 3922228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3932228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_1[1] = -accu3; 3942228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_0[0] = accu4; 3952228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 3962228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 3972228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if ( (M&1) == 0 ) 3982228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project { 3992228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* Last Sin and Cos value pair are the same */ 4002228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu1 = fMultDiv2(accu1, WTC(0x5a82799a)); 4012228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu2 = fMultDiv2(accu2, WTC(0x5a82799a)); 4022228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4032228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_1[0] = accu1 + accu2; 4042228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_0[1] = accu1 - accu2; 4052228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 4062228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 4072228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4082228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* Add twiddeling scale. */ 4092228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *pDat_e += 2; 4102228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} 4112228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#endif /* defined (FUNCTION_dct_IV) */ 4122228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4132228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#if !defined(FUNCTION_dst_IV) 4142228e360595641dd906bf1773307f43d304f5b2The Android Open Source Projectvoid dst_IV(FIXP_DBL *pDat, 4152228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int L, 4162228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int *pDat_e ) 4172228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project{ 4182228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int sin_step = 0; 4192228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int M = L >> 1; 4202228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4212228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const FIXP_WTP *twiddle; 4222228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const FIXP_STP *sin_twiddle; 4232228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4242228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#ifdef DSTIV2_ENABLE 4252228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (L == 2) { 4262228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project const FIXP_STP tab = STCP(0x7641AF3D, 0x30FB9452); 4272228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL tmp1, tmp2; 4282228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4292228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project cplxMultDiv2(&tmp2, &tmp1, pDat[0], pDat[1], tab); 4302228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4312228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat[0] = tmp1; 4322228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat[1] = tmp2; 4332228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4342228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *pDat_e += 1; 4352228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4362228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project return; 4372228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 4382228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#else 4392228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FDK_ASSERT(L >= 4); 4402228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#endif 4412228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4422228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project getTables(&twiddle, &sin_twiddle, &sin_step, L); 4432228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4442228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#ifdef FUNCTION_dst_IV_func1 4452228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if ( (M>=4) && ((M&3) == 0) ) { 4462228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project dst_IV_func1(M, twiddle, &pDat[0], &pDat[L]); 4472228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } else 4482228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#endif 4492228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project { 4502228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *RESTRICT pDat_0 = &pDat[0]; 4512228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *RESTRICT pDat_1 = &pDat[L - 2]; 4522228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4532228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project register int i; 4542228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4552228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* 34 cycles on ARM926 */ 4562228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for (i = 0; i < M-1; i+=2,pDat_0+=2,pDat_1-=2) 4572228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project { 4582228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project register FIXP_DBL accu1,accu2,accu3,accu4; 4592228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4602228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu1 = pDat_1[1]; accu2 = -pDat_0[0]; 4612228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu3 = pDat_0[1]; accu4 = -pDat_1[0]; 4622228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4632228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project cplxMultDiv2(&accu1, &accu2, accu1, accu2, twiddle[i]); 4642228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project cplxMultDiv2(&accu3, &accu4, accu4, accu3, twiddle[i+1]); 4652228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4662228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_0[0] = accu2; pDat_0[1] = accu1; 4672228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_1[0] = accu4; pDat_1[1] = -accu3; 4682228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 4692228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if (M&1) 4702228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project { 4712228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project register FIXP_DBL accu1,accu2; 4722228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4732228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu1 = pDat_1[1]; accu2 = -pDat_0[0]; 4742228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4752228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project cplxMultDiv2(&accu1, &accu2, accu1, accu2, twiddle[i]); 4762228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4772228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_0[0] = accu2; pDat_0[1] = accu1; 4782228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 4792228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 4802228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4812228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project fft(M, pDat, pDat_e); 4822228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4832228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#ifdef FUNCTION_dst_IV_func2 4842228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if ( (M>=4) && ((M&3) == 0) ) { 4852228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project dst_IV_func2(M>>2, sin_twiddle + sin_step, &pDat[0], &pDat[L - 1], sin_step); 4862228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } else 4872228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#endif /* FUNCTION_dst_IV_func2 */ 4882228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project { 4892228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *RESTRICT pDat_0; 4902228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_DBL *RESTRICT pDat_1; 4912228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project register FIXP_DBL accu1,accu2,accu3,accu4; 4922228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project int idx, i; 4932228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4942228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_0 = &pDat[0]; 4952228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_1 = &pDat[L - 2]; 4962228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 4972228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* Sin and Cos values are 0.0f and 1.0f */ 4982228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu1 = pDat_1[0]; 4992228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu2 = pDat_1[1]; 5002228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 5012228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_1[1] = -(pDat_0[0]>>1); 5022228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_0[0] = (pDat_0[1]>>1); 5032228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 5042228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project for (idx = sin_step,i=1; i<(M+1)>>1; i++, idx+=sin_step) 5052228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project { 5062228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project FIXP_STP twd = sin_twiddle[idx]; 5072228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 5082228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project cplxMultDiv2(&accu3, &accu4, accu1, accu2, twd); 5092228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_1[0] = -accu3; 5102228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_0[1] = -accu4; 5112228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 5122228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_0+=2; 5132228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_1-=2; 5142228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 5152228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project cplxMultDiv2(&accu3, &accu4, pDat_0[1], pDat_0[0], twd); 5162228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 5172228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu1 = pDat_1[0]; 5182228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu2 = pDat_1[1]; 5192228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 5202228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_0[0] = accu3; 5212228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_1[1] = -accu4; 5222228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 5232228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 5242228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project if ( (M&1) == 0 ) 5252228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project { 5262228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* Last Sin and Cos value pair are the same */ 5272228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu1 = fMultDiv2(accu1, WTC(0x5a82799a)); 5282228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project accu2 = fMultDiv2(accu2, WTC(0x5a82799a)); 5292228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 5302228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_0[1] = - accu1 - accu2; 5312228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project pDat_1[0] = accu2 - accu1; 5322228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 5332228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project } 5342228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 5352228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project /* Add twiddeling scale. */ 5362228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project *pDat_e += 2; 5372228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project} 5382228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project#endif /* !defined(FUNCTION_dst_IV) */ 5392228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 5402228e360595641dd906bf1773307f43d304f5b2The Android Open Source Project 541