Snapshot 2bda038c163298531d47394bc2c09e1409c5d0db

Change-Id: If584e579464f28b97d50e51fc76ba654a5536c54

1 files changed, 709 insertions, 0 deletions

diff --git a/libFDK/src/FDK_hybrid.cpp b/libFDK/src/FDK_hybrid.cpp
new file mode 100644
index 0000000..3d3ab7b
--- /dev/null
+++ b/libFDK/src/FDK_hybrid.cpp
@@ -0,0 +1,709 @@
+
+/* -----------------------------------------------------------------------------------------------------------
+Software License for The Fraunhofer FDK AAC Codec Library for Android
+
+© Copyright  1995 - 2012 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V.
+  All rights reserved.
+
+ 1.    INTRODUCTION
+The Fraunhofer FDK AAC Codec Library for Android ("FDK AAC Codec") is software that implements
+the MPEG Advanced Audio Coding ("AAC") encoding and decoding scheme for digital audio.
+This FDK AAC Codec software is intended to be used on a wide variety of Android devices.
+
+AAC's HE-AAC and HE-AAC v2 versions are regarded as today's most efficient general perceptual
+audio codecs. AAC-ELD is considered the best-performing full-bandwidth communications codec by
+independent studies and is widely deployed. AAC has been standardized by ISO and IEC as part
+of the MPEG specifications.
+
+Patent licenses for necessary patent claims for the FDK AAC Codec (including those of Fraunhofer)
+may be obtained through Via Licensing (www.vialicensing.com) or through the respective patent owners
+individually for the purpose of encoding or decoding bit streams in products that are compliant with
+the ISO/IEC MPEG audio standards. Please note that most manufacturers of Android devices already license
+these patent claims through Via Licensing or directly from the patent owners, and therefore FDK AAC Codec
+software may already be covered under those patent licenses when it is used for those licensed purposes only.
+
+Commercially-licensed AAC software libraries, including floating-point versions with enhanced sound quality,
+are also available from Fraunhofer. Users are encouraged to check the Fraunhofer website for additional
+applications information and documentation.
+
+2.    COPYRIGHT LICENSE
+
+Redistribution and use in source and binary forms, with or without modification, are permitted without
+payment of copyright license fees provided that you satisfy the following conditions:
+
+You must retain the complete text of this software license in redistributions of the FDK AAC Codec or
+your modifications thereto in source code form.
+
+You must retain the complete text of this software license in the documentation and/or other materials
+provided with redistributions of the FDK AAC Codec or your modifications thereto in binary form.
+You must make available free of charge copies of the complete source code of the FDK AAC Codec and your
+modifications thereto to recipients of copies in binary form.
+
+The name of Fraunhofer may not be used to endorse or promote products derived from this library without
+prior written permission.
+
+You may not charge copyright license fees for anyone to use, copy or distribute the FDK AAC Codec
+software or your modifications thereto.
+
+Your modified versions of the FDK AAC Codec must carry prominent notices stating that you changed the software
+and the date of any change. For modified versions of the FDK AAC Codec, the term
+"Fraunhofer FDK AAC Codec Library for Android" must be replaced by the term
+"Third-Party Modified Version of the Fraunhofer FDK AAC Codec Library for Android."
+
+3.    NO PATENT LICENSE
+
+NO EXPRESS OR IMPLIED LICENSES TO ANY PATENT CLAIMS, including without limitation the patents of Fraunhofer,
+ARE GRANTED BY THIS SOFTWARE LICENSE. Fraunhofer provides no warranty of patent non-infringement with
+respect to this software.
+
+You may use this FDK AAC Codec software or modifications thereto only for purposes that are authorized
+by appropriate patent licenses.
+
+4.    DISCLAIMER
+
+This FDK AAC Codec software is provided by Fraunhofer on behalf of the copyright holders and contributors
+"AS IS" and WITHOUT ANY EXPRESS OR IMPLIED WARRANTIES, including but not limited to the implied warranties
+of merchantability and fitness for a particular purpose. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR
+CONTRIBUTORS BE LIABLE for any direct, indirect, incidental, special, exemplary, or consequential damages,
+including but not limited to procurement of substitute goods or services; loss of use, data, or profits,
+or business interruption, however caused and on any theory of liability, whether in contract, strict
+liability, or tort (including negligence), arising in any way out of the use of this software, even if
+advised of the possibility of such damage.
+
+5.    CONTACT INFORMATION
+
+Fraunhofer Institute for Integrated Circuits IIS
+Attention: Audio and Multimedia Departments - FDK AAC LL
+Am Wolfsmantel 33
+91058 Erlangen, Germany
+
+www.iis.fraunhofer.de/amm
+amm-info@iis.fraunhofer.de
+----------------------------------------------------------------------------------------------------------- */
+
+/***************************  Fraunhofer IIS FDK Tools  **********************
+
+   Author(s): Markus Lohwasser
+   Description: FDK Tools Hybrid Filterbank
+
+******************************************************************************/
+
+#include "FDK_hybrid.h"
+
+
+#include "fft.h"
+
+/*--------------- defines -----------------------------*/
+#define FFT_IDX_R(a)  (2*a)
+#define FFT_IDX_I(a)  (2*a+1)
+
+#define HYB_COEF8_0  ( 0.00746082949812f )
+#define HYB_COEF8_1  ( 0.02270420949825f )
+#define HYB_COEF8_2  ( 0.04546865930473f )
+#define HYB_COEF8_3  ( 0.07266113929591f )
+#define HYB_COEF8_4  ( 0.09885108575264f )
+#define HYB_COEF8_5  ( 0.11793710567217f )
+#define HYB_COEF8_6  ( 0.12500000000000f )
+#define HYB_COEF8_7  ( HYB_COEF8_5 )
+#define HYB_COEF8_8  ( HYB_COEF8_4 )
+#define HYB_COEF8_9  ( HYB_COEF8_3 )
+#define HYB_COEF8_10 ( HYB_COEF8_2 )
+#define HYB_COEF8_11 ( HYB_COEF8_1 )
+#define HYB_COEF8_12 ( HYB_COEF8_0 )
+
+
+/*--------------- structure definitions ---------------*/
+
+#if defined(ARCH_PREFER_MULT_32x16)
+  #define FIXP_HTB FIXP_SGL               /* SGL data type. */
+  #define FIXP_HTP FIXP_SPK               /* Packed SGL data type. */
+  #define HTC(a) (FX_DBL2FXCONST_SGL(a))  /* Cast to SGL */
+  #define FL2FXCONST_HTB FL2FXCONST_SGL
+#else
+  #define FIXP_HTB FIXP_DBL               /* SGL data type. */
+  #define FIXP_HTP FIXP_DPK               /* Packed DBL data type. */
+  #define HTC(a) ((FIXP_DBL)(LONG)(a))    /* Cast to DBL */
+  #define FL2FXCONST_HTB FL2FXCONST_DBL
+#endif
+
+#define HTCP(real,imag) { { HTC(real), HTC(imag) } } /* How to arrange the packed values. */
+
+
+struct FDK_HYBRID_SETUP
+{
+    UCHAR               nrQmfBands;          /*!< Number of QMF bands to be converted to hybrid. */
+    UCHAR               nHybBands[3];        /*!< Number of Hybrid bands generated by nrQmfBands. */
+    SCHAR               kHybrid[3];          /*!< Filter configuration of each QMF band. */
+    UCHAR               protoLen;            /*!< Prototype filter length. */
+    UCHAR               filterDelay;         /*!< Delay caused by hybrid filter. */
+    const INT          *pReadIdxTable;       /*!< Helper table to access input data ringbuffer. */
+
+};
+
+/*--------------- constants ---------------------------*/
+static const INT ringbuffIdxTab[2*13] = { 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12 };
+
+static const FDK_HYBRID_SETUP setup_3_16 = { 3, { 8, 4, 4}, {  8,  4,  4}, 13, (13-1)/2, ringbuffIdxTab};
+static const FDK_HYBRID_SETUP setup_3_12 = { 3, { 8, 2, 2}, {  8,  2,  2}, 13, (13-1)/2, ringbuffIdxTab};
+static const FDK_HYBRID_SETUP setup_3_10 = { 3, { 6, 2, 2}, { -8, -2,  2}, 13, (13-1)/2, ringbuffIdxTab};
+
+
+static const FIXP_HTP HybFilterCoef8[] = {
+  HTCP(0x10000000, 0x00000000), HTCP(0x0df26407, 0xfa391882), HTCP(0xff532109, 0x00acdef7), HTCP(0x08f26d36, 0xf70d92ca),
+  HTCP(0xfee34b5f, 0x02af570f), HTCP(0x038f276e, 0xf7684793), HTCP(0x00000000, 0x05d1eac2), HTCP(0x00000000, 0x05d1eac2),
+  HTCP(0x038f276e, 0x0897b86d), HTCP(0xfee34b5f, 0xfd50a8f1), HTCP(0x08f26d36, 0x08f26d36), HTCP(0xff532109, 0xff532109),
+  HTCP(0x0df26407, 0x05c6e77e)
+};
+
+static const FIXP_HTB HybFilterCoef2[13] = {
+  FL2FXCONST_HTB( 0.00000000000000f), FL2FXCONST_HTB( 0.01899487526049f), FL2FXCONST_HTB( 0.00000000000000f), FL2FXCONST_HTB(-0.07293139167538f), FL2FXCONST_HTB( 0.00000000000000f), FL2FXCONST_HTB( 0.30596630545168f),
+  FL2FXCONST_HTB( 0.50000000000000f), FL2FXCONST_HTB( 0.30596630545168f), FL2FXCONST_HTB( 0.00000000000000f), FL2FXCONST_HTB(-0.07293139167538f), FL2FXCONST_HTB( 0.00000000000000f), FL2FXCONST_HTB( 0.01899487526049f),
+  FL2FXCONST_HTB( 0.00000000000000f)
+};
+
+static const FIXP_HTB HybFilterCoef4[13] = {
+  FL2FXCONST_HTB(-0.00305151927305f), FL2FXCONST_HTB(-0.00794862316203f), FL2FXCONST_HTB(              0.0f), FL2FXCONST_HTB( 0.04318924038756f), FL2FXCONST_HTB( 0.12542448210445f), FL2FXCONST_HTB( 0.21227807049160f),
+  FL2FXCONST_HTB(             0.25f), FL2FXCONST_HTB( 0.21227807049160f), FL2FXCONST_HTB( 0.12542448210445f), FL2FXCONST_HTB( 0.04318924038756f), FL2FXCONST_HTB(              0.0f), FL2FXCONST_HTB(-0.00794862316203f),
+  FL2FXCONST_HTB(-0.00305151927305f)
+};
+
+/*--------------- function declarations ---------------*/
+static INT kChannelFiltering(
+        const FIXP_DBL *const      pQmfReal,
+        const FIXP_DBL *const      pQmfImag,
+        const INT *const           pReadIdx,
+        FIXP_DBL *const            mHybridReal,
+        FIXP_DBL *const            mHybridImag,
+        const SCHAR                hybridConfig
+        );
+
+
+/*--------------- function definitions ----------------*/
+
+INT FDKhybridAnalysisOpen(
+        HANDLE_FDK_ANA_HYB_FILTER  hAnalysisHybFilter,
+        FIXP_DBL *const            pLFmemory,
+        const UINT                 LFmemorySize,
+        FIXP_DBL *const            pHFmemory,
+        const UINT                 HFmemorySize
+        )
+{
+  INT err = 0;
+
+  /* Save pointer to extern memory. */
+  hAnalysisHybFilter->pLFmemory    = pLFmemory;
+  hAnalysisHybFilter->LFmemorySize = LFmemorySize;
+
+  hAnalysisHybFilter->pHFmemory    = pHFmemory;
+  hAnalysisHybFilter->HFmemorySize = HFmemorySize;
+
+  return err;
+}
+
+INT FDKhybridAnalysisInit(
+        HANDLE_FDK_ANA_HYB_FILTER  hAnalysisHybFilter,
+        const FDK_HYBRID_MODE      mode,
+        const INT                  qmfBands,
+        const INT                  cplxBands,
+        const INT                  initStatesFlag
+        )
+{
+    int k;
+    INT err = 0;
+    FIXP_DBL *pMem = NULL;
+    HANDLE_FDK_HYBRID_SETUP setup = NULL;
+
+    switch (mode) {
+      case THREE_TO_TEN:     setup = (HANDLE_FDK_HYBRID_SETUP)&setup_3_10; break;
+      case THREE_TO_TWELVE:  setup = (HANDLE_FDK_HYBRID_SETUP)&setup_3_12; break;
+      case THREE_TO_SIXTEEN: setup = (HANDLE_FDK_HYBRID_SETUP)&setup_3_16; break;
+      default:               err = -1; goto bail;
+    }
+
+    /* Initialize handle. */
+    hAnalysisHybFilter->pSetup      = setup;
+    hAnalysisHybFilter->bufferLFpos = setup->protoLen-1;
+    hAnalysisHybFilter->bufferHFpos = 0;
+    hAnalysisHybFilter->nrBands     = qmfBands;
+    hAnalysisHybFilter->cplxBands   = cplxBands;
+    hAnalysisHybFilter->hfMode      = 0;
+
+    /* Check available memory. */
+    if ( ((2*setup->nrQmfBands*setup->protoLen*sizeof(FIXP_DBL)) > hAnalysisHybFilter->LFmemorySize)
+      || ((setup->filterDelay*((qmfBands-setup->nrQmfBands)+(cplxBands-setup->nrQmfBands))*sizeof(FIXP_DBL)) > hAnalysisHybFilter->HFmemorySize) )
+    {
+      err = -2;
+      goto bail;
+    }
+
+    /* Distribut LF memory. */
+    pMem = hAnalysisHybFilter->pLFmemory;
+    for (k=0; k<setup->nrQmfBands; k++) {
+      hAnalysisHybFilter->bufferLFReal[k] = pMem; pMem += setup->protoLen;
+      hAnalysisHybFilter->bufferLFImag[k] = pMem; pMem += setup->protoLen;
+    }
+
+    /* Distribut HF memory. */
+    pMem = hAnalysisHybFilter->pHFmemory;
+    for (k=0; k<setup->filterDelay; k++) {
+      hAnalysisHybFilter->bufferHFReal[k] = pMem; pMem += (qmfBands-setup->nrQmfBands);
+      hAnalysisHybFilter->bufferHFImag[k] = pMem; pMem += (cplxBands-setup->nrQmfBands);
+    }
+
+    if (initStatesFlag) {
+      /* Clear LF buffer */
+      for (k=0; k<setup->nrQmfBands; k++) {
+        FDKmemclear(hAnalysisHybFilter->bufferLFReal[k], setup->protoLen*sizeof(FIXP_DBL));
+        FDKmemclear(hAnalysisHybFilter->bufferLFImag[k], setup->protoLen*sizeof(FIXP_DBL));
+      }
+
+      if (qmfBands > setup->nrQmfBands) {
+      /* Clear HF buffer */
+      for (k=0; k<setup->filterDelay; k++) {
+        FDKmemclear(hAnalysisHybFilter->bufferHFReal[k], (qmfBands-setup->nrQmfBands)*sizeof(FIXP_DBL));
+        FDKmemclear(hAnalysisHybFilter->bufferHFImag[k], (cplxBands-setup->nrQmfBands)*sizeof(FIXP_DBL));
+      }
+    }
+    }
+
+bail:
+    return err;
+}
+
+INT FDKhybridAnalysisScaleStates(
+        HANDLE_FDK_ANA_HYB_FILTER  hAnalysisHybFilter,
+        const INT                  scalingValue
+        )
+{
+    INT err = 0;
+
+    if (hAnalysisHybFilter==NULL) {
+      err = 1; /* invalid handle */
+    }
+    else {
+      int k;
+      HANDLE_FDK_HYBRID_SETUP setup = hAnalysisHybFilter->pSetup;
+
+      /* Scale LF buffer */
+      for (k=0; k<setup->nrQmfBands; k++) {
+        scaleValues(hAnalysisHybFilter->bufferLFReal[k], setup->protoLen, scalingValue);
+        scaleValues(hAnalysisHybFilter->bufferLFImag[k], setup->protoLen, scalingValue);
+      }
+      if (hAnalysisHybFilter->nrBands > setup->nrQmfBands) {
+        /* Scale HF buffer */
+        for (k=0; k<setup->filterDelay; k++) {
+          scaleValues(hAnalysisHybFilter->bufferHFReal[k], (hAnalysisHybFilter->nrBands-setup->nrQmfBands), scalingValue);
+          scaleValues(hAnalysisHybFilter->bufferHFImag[k], (hAnalysisHybFilter->cplxBands-setup->nrQmfBands), scalingValue);
+        }
+      }
+    }
+    return err;
+}
+
+INT FDKhybridAnalysisApply(
+        HANDLE_FDK_ANA_HYB_FILTER  hAnalysisHybFilter,
+        const FIXP_DBL *const      pQmfReal,
+        const FIXP_DBL *const      pQmfImag,
+        FIXP_DBL *const            pHybridReal,
+        FIXP_DBL *const            pHybridImag)
+{
+    int k, hybOffset = 0;
+    INT err = 0;
+    const int nrQmfBandsLF = hAnalysisHybFilter->pSetup->nrQmfBands; /* number of QMF bands to be converted to hybrid */
+
+    const int writIndex = hAnalysisHybFilter->bufferLFpos;
+    int readIndex = hAnalysisHybFilter->bufferLFpos;
+
+    if (++readIndex>=hAnalysisHybFilter->pSetup->protoLen) readIndex = 0;
+    const INT* pBufferLFreadIdx = &hAnalysisHybFilter->pSetup->pReadIdxTable[readIndex];
+
+    /*
+     * LF buffer.
+     */
+    for (k=0; k<nrQmfBandsLF; k++) {
+        /* New input sample. */
+        hAnalysisHybFilter->bufferLFReal[k][writIndex] = pQmfReal[k];
+        hAnalysisHybFilter->bufferLFImag[k][writIndex] = pQmfImag[k];
+
+    /* Perform hybrid filtering. */
+        kChannelFiltering(
+                hAnalysisHybFilter->bufferLFReal[k],
+                hAnalysisHybFilter->bufferLFImag[k],
+                pBufferLFreadIdx,
+                pHybridReal+hybOffset,
+                pHybridImag+hybOffset,
+                hAnalysisHybFilter->pSetup->kHybrid[k]);
+
+        hybOffset += hAnalysisHybFilter->pSetup->nHybBands[k];
+    }
+
+    hAnalysisHybFilter->bufferLFpos = readIndex; /* Index where to write next input sample. */
+
+    if (hAnalysisHybFilter->nrBands > nrQmfBandsLF) {
+    /*
+     * HF buffer.
+     */
+    if (hAnalysisHybFilter->hfMode!=0) {
+        /* HF delay compensation was applied outside. */
+        FDKmemcpy(pHybridReal+hybOffset, &pQmfReal[nrQmfBandsLF], (hAnalysisHybFilter->nrBands-nrQmfBandsLF)*sizeof(FIXP_DBL));
+        FDKmemcpy(pHybridImag+hybOffset, &pQmfImag[nrQmfBandsLF], (hAnalysisHybFilter->cplxBands-nrQmfBandsLF)*sizeof(FIXP_DBL));
+    }
+    else {
+        /* HF delay compensation, filterlength/2. */
+    FDKmemcpy(pHybridReal+hybOffset, hAnalysisHybFilter->bufferHFReal[hAnalysisHybFilter->bufferHFpos], (hAnalysisHybFilter->nrBands-nrQmfBandsLF)*sizeof(FIXP_DBL));
+    FDKmemcpy(pHybridImag+hybOffset, hAnalysisHybFilter->bufferHFImag[hAnalysisHybFilter->bufferHFpos], (hAnalysisHybFilter->cplxBands-nrQmfBandsLF)*sizeof(FIXP_DBL));
+
+    FDKmemcpy(hAnalysisHybFilter->bufferHFReal[hAnalysisHybFilter->bufferHFpos], &pQmfReal[nrQmfBandsLF], (hAnalysisHybFilter->nrBands-nrQmfBandsLF)*sizeof(FIXP_DBL));
+    FDKmemcpy(hAnalysisHybFilter->bufferHFImag[hAnalysisHybFilter->bufferHFpos], &pQmfImag[nrQmfBandsLF], (hAnalysisHybFilter->cplxBands-nrQmfBandsLF)*sizeof(FIXP_DBL));
+
+    if (++hAnalysisHybFilter->bufferHFpos>=hAnalysisHybFilter->pSetup->filterDelay) hAnalysisHybFilter->bufferHFpos = 0;
+    }
+    } /* process HF part*/
+
+    return err;
+}
+
+INT FDKhybridAnalysisClose(
+        HANDLE_FDK_ANA_HYB_FILTER hAnalysisHybFilter
+        )
+{
+  INT err = 0;
+
+  if (hAnalysisHybFilter != NULL) {
+    hAnalysisHybFilter->pLFmemory    = NULL;
+    hAnalysisHybFilter->pHFmemory    = NULL;
+    hAnalysisHybFilter->LFmemorySize = 0;
+    hAnalysisHybFilter->HFmemorySize = 0;
+  }
+
+  return err;
+}
+
+INT FDKhybridSynthesisInit(
+        HANDLE_FDK_SYN_HYB_FILTER  hSynthesisHybFilter,
+        const FDK_HYBRID_MODE      mode,
+        const INT                  qmfBands,
+        const INT                  cplxBands
+        )
+{
+    INT err = 0;
+    HANDLE_FDK_HYBRID_SETUP setup = NULL;
+
+    switch (mode) {
+      case THREE_TO_TEN:     setup = (HANDLE_FDK_HYBRID_SETUP)&setup_3_10; break;
+      case THREE_TO_TWELVE:  setup = (HANDLE_FDK_HYBRID_SETUP)&setup_3_12; break;
+      case THREE_TO_SIXTEEN: setup = (HANDLE_FDK_HYBRID_SETUP)&setup_3_16; break;
+      default:               err = -1; goto bail;
+    }
+
+    hSynthesisHybFilter->pSetup      = setup;
+    hSynthesisHybFilter->nrBands     = qmfBands;
+    hSynthesisHybFilter->cplxBands   = cplxBands;
+
+bail:
+    return err;
+}
+
+
+INT FDKhybridSynthesisApply(
+        HANDLE_FDK_SYN_HYB_FILTER  hSynthesisHybFilter,
+        const FIXP_DBL *const      pHybridReal,
+        const FIXP_DBL *const      pHybridImag,
+        FIXP_DBL *const            pQmfReal,
+        FIXP_DBL *const            pQmfImag
+        )
+{
+    int k, n, hybOffset=0;
+    INT err = 0;
+    const INT nrQmfBandsLF = hSynthesisHybFilter->pSetup->nrQmfBands;
+
+    /*
+     * LF buffer.
+     */
+    for (k=0; k<nrQmfBandsLF; k++) {
+      const int nHybBands = hSynthesisHybFilter->pSetup->nHybBands[k];
+
+      FIXP_DBL accu1 = FL2FXCONST_DBL(0.f);
+      FIXP_DBL accu2 = FL2FXCONST_DBL(0.f);
+
+      /* Perform hybrid filtering. */
+      for (n=0; n<nHybBands; n++) {
+          accu1 += pHybridReal[hybOffset+n];
+          accu2 += pHybridImag[hybOffset+n];
+      }
+      pQmfReal[k] = accu1;
+      pQmfImag[k] = accu2;
+
+      hybOffset += nHybBands;
+    }
+
+    if (hSynthesisHybFilter->nrBands > nrQmfBandsLF) {
+      /*
+       * HF buffer.
+       */
+    FDKmemcpy(&pQmfReal[nrQmfBandsLF], &pHybridReal[hybOffset], (hSynthesisHybFilter->nrBands-nrQmfBandsLF)*sizeof(FIXP_DBL));
+    FDKmemcpy(&pQmfImag[nrQmfBandsLF], &pHybridImag[hybOffset], (hSynthesisHybFilter->cplxBands-nrQmfBandsLF)*sizeof(FIXP_DBL));
+    }
+
+    return err;
+}
+
+static void dualChannelFiltering(
+        const FIXP_DBL *const      pQmfReal,
+        const FIXP_DBL *const      pQmfImag,
+        const INT *const           pReadIdx,
+        FIXP_DBL *const            mHybridReal,
+        FIXP_DBL *const            mHybridImag,
+        const INT                  invert
+        )
+{
+    const FIXP_HTB *p = HybFilterCoef2;
+
+    FIXP_DBL  r1, r6;
+    FIXP_DBL  i1, i6;
+
+    /* symmetric filter coefficients */
+    r1  = fMultDiv2(p[1], pQmfReal[pReadIdx[1]]) + fMultDiv2(p[1], pQmfReal[pReadIdx[11]]) ;
+    i1  = fMultDiv2(p[1], pQmfImag[pReadIdx[1]]) + fMultDiv2(p[1], pQmfImag[pReadIdx[11]]) ;
+    r1 += fMultDiv2(p[3], pQmfReal[pReadIdx[3]]) + fMultDiv2(p[3], pQmfReal[pReadIdx[ 9]]) ;
+    i1 += fMultDiv2(p[3], pQmfImag[pReadIdx[3]]) + fMultDiv2(p[3], pQmfImag[pReadIdx[ 9]]) ;
+    r1 += fMultDiv2(p[5], pQmfReal[pReadIdx[5]]) + fMultDiv2(p[5], pQmfReal[pReadIdx[ 7]]) ;
+    i1 += fMultDiv2(p[5], pQmfImag[pReadIdx[5]]) + fMultDiv2(p[5], pQmfImag[pReadIdx[ 7]]) ;
+    r6  = fMultDiv2(p[6], pQmfReal[pReadIdx[6]]) ;
+    i6  = fMultDiv2(p[6], pQmfImag[pReadIdx[6]]) ;
+
+    if (invert) {
+      mHybridReal[1] = (r1 + r6) << 1;
+      mHybridImag[1] = (i1 + i6) << 1;
+
+      mHybridReal[0] = (r6 - r1) << 1;
+      mHybridImag[0] = (i6 - i1) << 1;
+    }
+    else {
+      mHybridReal[0] = (r1 + r6) << 1;
+      mHybridImag[0] = (i1 + i6) << 1;
+
+      mHybridReal[1] = (r6 - r1) << 1;
+      mHybridImag[1] = (i6 - i1) << 1;
+    }
+}
+
+static void fourChannelFiltering(
+        const FIXP_DBL *const      pQmfReal,
+        const FIXP_DBL *const      pQmfImag,
+        const INT *const           pReadIdx,
+        FIXP_DBL *const            mHybridReal,
+        FIXP_DBL *const            mHybridImag,
+        const INT                  invert
+        )
+{
+    const FIXP_HTB *p = HybFilterCoef4;
+
+    FIXP_DBL fft[8];
+
+    static const FIXP_DBL  cr[13] = {
+      FL2FXCONST_DBL(               0.f), FL2FXCONST_DBL(-0.70710678118655f), FL2FXCONST_DBL(              -1.f),
+      FL2FXCONST_DBL(-0.70710678118655f), FL2FXCONST_DBL(               0.f), FL2FXCONST_DBL( 0.70710678118655f),
+      FL2FXCONST_DBL(               1.f),
+      FL2FXCONST_DBL( 0.70710678118655f), FL2FXCONST_DBL(               0.f), FL2FXCONST_DBL(-0.70710678118655f),
+      FL2FXCONST_DBL(              -1.f), FL2FXCONST_DBL(-0.70710678118655f), FL2FXCONST_DBL(               0.f)
+    };
+    static const FIXP_DBL  ci[13] = {
+      FL2FXCONST_DBL(              -1.f), FL2FXCONST_DBL(-0.70710678118655f), FL2FXCONST_DBL(               0.f),
+      FL2FXCONST_DBL( 0.70710678118655f), FL2FXCONST_DBL(               1.f), FL2FXCONST_DBL( 0.70710678118655f),
+      FL2FXCONST_DBL(               0.f),
+      FL2FXCONST_DBL(-0.70710678118655f), FL2FXCONST_DBL(              -1.f), FL2FXCONST_DBL(-0.70710678118655f),
+      FL2FXCONST_DBL(               0.f), FL2FXCONST_DBL( 0.70710678118655f), FL2FXCONST_DBL(               1.f)
+    };
+
+
+    /* FIR filter. */
+    /* pre twiddeling with pre-twiddling coefficients c[n]  */
+    /* multiplication with filter coefficients p[n]         */
+    /* hint: (a + ib)*(c + id) = (a*c - b*d) + i(a*d + b*c) */
+    /* write to fft coefficient n'                          */
+    fft[FFT_IDX_R(0)] =  ( fMult(p[10], ( fMultSub(fMultDiv2(cr[ 2], pQmfReal[pReadIdx[ 2]]), ci[ 2], pQmfImag[pReadIdx[ 2]]))) +
+                           fMult(p[ 6], ( fMultSub(fMultDiv2(cr[ 6], pQmfReal[pReadIdx[ 6]]), ci[ 6], pQmfImag[pReadIdx[ 6]]))) +
+                           fMult(p[ 2], ( fMultSub(fMultDiv2(cr[10], pQmfReal[pReadIdx[10]]), ci[10], pQmfImag[pReadIdx[10]]))) );
+    fft[FFT_IDX_I(0)] =  ( fMult(p[10], ( fMultAdd(fMultDiv2(ci[ 2], pQmfReal[pReadIdx[ 2]]), cr[ 2], pQmfImag[pReadIdx[ 2]]))) +
+                           fMult(p[ 6], ( fMultAdd(fMultDiv2(ci[ 6], pQmfReal[pReadIdx[ 6]]), cr[ 6], pQmfImag[pReadIdx[ 6]]))) +
+                           fMult(p[ 2], ( fMultAdd(fMultDiv2(ci[10], pQmfReal[pReadIdx[10]]), cr[10], pQmfImag[pReadIdx[10]]))) );
+
+    /* twiddle dee dum */
+    fft[FFT_IDX_R(1)] =  ( fMult(p[ 9], ( fMultSub(fMultDiv2(cr[ 3], pQmfReal[pReadIdx[ 3]]), ci[ 3], pQmfImag[pReadIdx[ 3]]))) +
+                           fMult(p[ 5], ( fMultSub(fMultDiv2(cr[ 7], pQmfReal[pReadIdx[ 7]]), ci[ 7], pQmfImag[pReadIdx[ 7]]))) +
+                           fMult(p[ 1], ( fMultSub(fMultDiv2(cr[11], pQmfReal[pReadIdx[11]]), ci[11], pQmfImag[pReadIdx[11]]))) );
+    fft[FFT_IDX_I(1)] =  ( fMult(p[ 9], ( fMultAdd(fMultDiv2(ci[ 3], pQmfReal[pReadIdx[ 3]]), cr[ 3], pQmfImag[pReadIdx[ 3]]))) +
+                           fMult(p[ 5], ( fMultAdd(fMultDiv2(ci[ 7], pQmfReal[pReadIdx[ 7]]), cr[ 7], pQmfImag[pReadIdx[ 7]]))) +
+                           fMult(p[ 1], ( fMultAdd(fMultDiv2(ci[11], pQmfReal[pReadIdx[11]]), cr[11], pQmfImag[pReadIdx[11]]))) );
+
+    /* twiddle dee dee */
+    fft[FFT_IDX_R(2)] =  ( fMult(p[12], ( fMultSub(fMultDiv2(cr[ 0], pQmfReal[pReadIdx[ 0]]), ci[ 0], pQmfImag[pReadIdx[ 0]]))) +
+                           fMult(p[ 8], ( fMultSub(fMultDiv2(cr[ 4], pQmfReal[pReadIdx[ 4]]), ci[ 4], pQmfImag[pReadIdx[ 4]]))) +
+                           fMult(p[ 4], ( fMultSub(fMultDiv2(cr[ 8], pQmfReal[pReadIdx[ 8]]), ci[ 8], pQmfImag[pReadIdx[ 8]]))) +
+                           fMult(p[ 0], ( fMultSub(fMultDiv2(cr[12], pQmfReal[pReadIdx[12]]), ci[12], pQmfImag[pReadIdx[12]]))) );
+    fft[FFT_IDX_I(2)] =  ( fMult(p[12], ( fMultAdd(fMultDiv2(ci[ 0], pQmfReal[pReadIdx[ 0]]), cr[ 0], pQmfImag[pReadIdx[ 0]]))) +
+                           fMult(p[ 8], ( fMultAdd(fMultDiv2(ci[ 4], pQmfReal[pReadIdx[ 4]]), cr[ 4], pQmfImag[pReadIdx[ 4]]))) +
+                           fMult(p[ 4], ( fMultAdd(fMultDiv2(ci[ 8], pQmfReal[pReadIdx[ 8]]), cr[ 8], pQmfImag[pReadIdx[ 8]]))) +
+                           fMult(p[ 0], ( fMultAdd(fMultDiv2(ci[12], pQmfReal[pReadIdx[12]]), cr[12], pQmfImag[pReadIdx[12]]))) );
+
+    fft[FFT_IDX_R(3)] =  ( fMult(p[11], ( fMultSub(fMultDiv2(cr[ 1], pQmfReal[pReadIdx[ 1]]), ci[ 1], pQmfImag[pReadIdx[ 1]]))) +
+                           fMult(p[ 7], ( fMultSub(fMultDiv2(cr[ 5], pQmfReal[pReadIdx[ 5]]), ci[ 5], pQmfImag[pReadIdx[ 5]]))) +
+                           fMult(p[ 3], ( fMultSub(fMultDiv2(cr[ 9], pQmfReal[pReadIdx[ 9]]), ci[ 9], pQmfImag[pReadIdx[ 9]]))) );
+    fft[FFT_IDX_I(3)] =  ( fMult(p[11], ( fMultAdd(fMultDiv2(ci[ 1], pQmfReal[pReadIdx[ 1]]), cr[ 1], pQmfImag[pReadIdx[ 1]]))) +
+                           fMult(p[ 7], ( fMultAdd(fMultDiv2(ci[ 5], pQmfReal[pReadIdx[ 5]]), cr[ 5], pQmfImag[pReadIdx[ 5]]))) +
+                           fMult(p[ 3], ( fMultAdd(fMultDiv2(ci[ 9], pQmfReal[pReadIdx[ 9]]), cr[ 9], pQmfImag[pReadIdx[ 9]]))) );
+
+    /* fft modulation                                                    */
+    /* here: fast manual fft modulation for a fft of length M=4          */
+    /* fft_4{x[n]} = x[0]*exp(-i*2*pi/4*m*0) + x[1]*exp(-i*2*pi/4*m*1) +
+    x[2]*exp(-i*2*pi/4*m*2) + x[3]*exp(-i*2*pi/4*m*3)   */
+
+    /*
+    fft bin m=0:
+    X[0, n] = x[0] +   x[1] + x[2] +   x[3]
+    */
+    mHybridReal[0] = fft[FFT_IDX_R(0)] + fft[FFT_IDX_R(1)] + fft[FFT_IDX_R(2)] + fft[FFT_IDX_R(3)];
+    mHybridImag[0] = fft[FFT_IDX_I(0)] + fft[FFT_IDX_I(1)] + fft[FFT_IDX_I(2)] + fft[FFT_IDX_I(3)];
+
+    /*
+    fft bin m=1:
+    X[1, n] = x[0] - i*x[1] - x[2] + i*x[3]
+    */
+    mHybridReal[1] = fft[FFT_IDX_R(0)] + fft[FFT_IDX_I(1)] - fft[FFT_IDX_R(2)] - fft[FFT_IDX_I(3)];
+    mHybridImag[1] = fft[FFT_IDX_I(0)] - fft[FFT_IDX_R(1)] - fft[FFT_IDX_I(2)] + fft[FFT_IDX_R(3)];
+
+    /*
+    fft bin m=2:
+    X[2, n] = x[0] -   x[1] + x[2] -   x[3]
+    */
+    mHybridReal[2] = fft[FFT_IDX_R(0)] - fft[FFT_IDX_R(1)] + fft[FFT_IDX_R(2)] - fft[FFT_IDX_R(3)];
+    mHybridImag[2] = fft[FFT_IDX_I(0)] - fft[FFT_IDX_I(1)] + fft[FFT_IDX_I(2)] - fft[FFT_IDX_I(3)];
+
+    /*
+    fft bin m=3:
+    X[3, n] = x[0] + j*x[1] - x[2] - j*x[3]
+    */
+    mHybridReal[3] = fft[FFT_IDX_R(0)] - fft[FFT_IDX_I(1)] - fft[FFT_IDX_R(2)] + fft[FFT_IDX_I(3)];
+    mHybridImag[3] = fft[FFT_IDX_I(0)] + fft[FFT_IDX_R(1)] - fft[FFT_IDX_I(2)] - fft[FFT_IDX_R(3)];
+}
+
+
+static void eightChannelFiltering(
+        const FIXP_DBL *const      pQmfReal,
+        const FIXP_DBL *const      pQmfImag,
+        const INT *const           pReadIdx,
+        FIXP_DBL *const            mHybridReal,
+        FIXP_DBL *const            mHybridImag,
+        const INT                  invert
+        )
+{
+    const FIXP_HTP *p = HybFilterCoef8;
+    INT k, sc;
+
+    FIXP_DBL mfft[16+ALIGNMENT_DEFAULT];
+    FIXP_DBL *pfft = (FIXP_DBL*)ALIGN_PTR(mfft);
+
+    FIXP_DBL accu1, accu2, accu3, accu4;
+
+    /* pre twiddeling */
+    pfft[FFT_IDX_R(0)] = fMultDiv2(p[0].v.re, pQmfReal[pReadIdx[6]]);
+    pfft[FFT_IDX_I(0)] = fMultDiv2(p[0].v.re, pQmfImag[pReadIdx[6]]);
+
+    cplxMultDiv2(&accu1, &accu2, pQmfReal[pReadIdx[7]], pQmfImag[pReadIdx[7]], p[1]);
+    pfft[FFT_IDX_R(1)] = accu1;
+    pfft[FFT_IDX_I(1)] = accu2;
+
+    cplxMultDiv2(&accu1, &accu2, pQmfReal[pReadIdx[0]], pQmfImag[pReadIdx[0]], p[2]);
+    cplxMultDiv2(&accu3, &accu4, pQmfReal[pReadIdx[8]], pQmfImag[pReadIdx[8]], p[3]);
+    pfft[FFT_IDX_R(2)] = accu1 + accu3;
+    pfft[FFT_IDX_I(2)] = accu2 + accu4;
+
+    cplxMultDiv2(&accu1, &accu2, pQmfReal[pReadIdx[1]], pQmfImag[pReadIdx[1]], p[4]);
+    cplxMultDiv2(&accu3, &accu4, pQmfReal[pReadIdx[9]], pQmfImag[pReadIdx[9]], p[5]);
+    pfft[FFT_IDX_R(3)] = accu1 + accu3;
+    pfft[FFT_IDX_I(3)] = accu2 + accu4;
+
+    pfft[FFT_IDX_R(4)] = fMultDiv2(pQmfImag[pReadIdx[10]], p[7].v.im) - fMultDiv2(pQmfImag[pReadIdx[ 2]], p[6].v.im);
+    pfft[FFT_IDX_I(4)] = fMultDiv2(pQmfReal[pReadIdx[ 2]], p[6].v.im) - fMultDiv2(pQmfReal[pReadIdx[10]], p[7].v.im);
+
+    cplxMultDiv2(&accu1, &accu2, pQmfReal[pReadIdx[ 3]], pQmfImag[pReadIdx[ 3]], p[8]);
+    cplxMultDiv2(&accu3, &accu4, pQmfReal[pReadIdx[11]], pQmfImag[pReadIdx[11]], p[9]);
+    pfft[FFT_IDX_R(5)] = accu1 + accu3;
+    pfft[FFT_IDX_I(5)] = accu2 + accu4;
+
+    cplxMultDiv2(&accu1, &accu2, pQmfReal[pReadIdx[ 4]], pQmfImag[pReadIdx[ 4]], p[10]);
+    cplxMultDiv2(&accu3, &accu4, pQmfReal[pReadIdx[12]], pQmfImag[pReadIdx[12]], p[11]);
+    pfft[FFT_IDX_R(6)] = accu1 + accu3;
+    pfft[FFT_IDX_I(6)] = accu2 + accu4;
+
+    cplxMultDiv2(&accu1, &accu2, pQmfReal[pReadIdx[ 5]], pQmfImag[pReadIdx[ 5]], p[12]);
+    pfft[FFT_IDX_R(7)] = accu1;
+    pfft[FFT_IDX_I(7)] = accu2;
+
+    /* fft modulation */
+    fft_8 (pfft);
+    sc = 1 + 2;
+
+    if (invert) {
+      mHybridReal[0]  = pfft[FFT_IDX_R(7)] << sc;
+      mHybridImag[0]  = pfft[FFT_IDX_I(7)] << sc;
+      mHybridReal[1]  = pfft[FFT_IDX_R(0)] << sc;
+      mHybridImag[1]  = pfft[FFT_IDX_I(0)] << sc;
+
+      mHybridReal[2]  = pfft[FFT_IDX_R(6)] << sc;
+      mHybridImag[2]  = pfft[FFT_IDX_I(6)] << sc;
+      mHybridReal[3]  = pfft[FFT_IDX_R(1)] << sc;
+      mHybridImag[3]  = pfft[FFT_IDX_I(1)] << sc;
+
+      mHybridReal[4]  = pfft[FFT_IDX_R(2)] << sc;
+      mHybridReal[4] += pfft[FFT_IDX_R(5)] << sc;
+      mHybridImag[4]  = pfft[FFT_IDX_I(2)] << sc;
+      mHybridImag[4] += pfft[FFT_IDX_I(5)] << sc;
+
+      mHybridReal[5]  = pfft[FFT_IDX_R(3)] << sc;
+      mHybridReal[5] += pfft[FFT_IDX_R(4)] << sc;
+      mHybridImag[5]  = pfft[FFT_IDX_I(3)] << sc;
+      mHybridImag[5] += pfft[FFT_IDX_I(4)] << sc;
+    }
+    else {
+      for(k=0; k<8;k++ ) {
+        mHybridReal[k] = pfft[FFT_IDX_R(k)] << sc;
+        mHybridImag[k] = pfft[FFT_IDX_I(k)] << sc;
+      }
+    }
+}
+
+static INT kChannelFiltering(
+        const FIXP_DBL *const      pQmfReal,
+        const FIXP_DBL *const      pQmfImag,
+        const INT *const           pReadIdx,
+        FIXP_DBL *const            mHybridReal,
+        FIXP_DBL *const            mHybridImag,
+        const SCHAR                hybridConfig
+        )
+{
+    INT err = 0;
+
+    switch (hybridConfig) {
+      case  2:
+      case -2:
+        dualChannelFiltering(pQmfReal, pQmfImag, pReadIdx, mHybridReal, mHybridImag, (hybridConfig<0) ? 1 : 0 );
+        break;
+      case  4:
+      case -4:
+        fourChannelFiltering(pQmfReal, pQmfImag, pReadIdx, mHybridReal, mHybridImag, (hybridConfig<0) ? 1 : 0 );
+        break;
+      case  8:
+      case -8:
+        eightChannelFiltering(pQmfReal, pQmfImag, pReadIdx, mHybridReal, mHybridImag, (hybridConfig<0) ? 1 : 0 );
+        break;
+      default:
+        err = -1;
+    }
+
+    return err;
+}
+
+
+