/* ------------------------------------------------------------------ * Copyright (C) 1998-2009 PacketVideo * * Licensed under the Apache License, Version 2.0 (the "License"); * you may not use this file except in compliance with the License. * You may obtain a copy of the License at * * http://www.apache.org/licenses/LICENSE-2.0 * * Unless required by applicable law or agreed to in writing, software * distributed under the License is distributed on an "AS IS" BASIS, * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either * express or implied. * See the License for the specific language governing permissions * and limitations under the License. * ------------------------------------------------------------------- */ /* Filename: ps_stereo_processing.c ------------------------------------------------------------------------------ REVISION HISTORY Who: Date: MM/DD/YYYY Description: ------------------------------------------------------------------------------ INPUT AND OUTPUT DEFINITIONS ------------------------------------------------------------------------------ FUNCTION DESCRIPTION Stereo Process or reconstruction l_k(n) = H11(k,n)*s_k(n) + H21(k,n)*d_k(n) r_k(n) = H12(k,n)*s_k(n) + H22(k,n)*d_k(n) _______ ________ | | _______ | | ->|Hybrid | LF ---- | |->| Hybrid |--> | Anal. | | | | | Synth | QMF -> L ------- o----------------------->| | -------- Synth QMF | s_k(n) |Stereo |--------------> Anal. -------------------------->| | _______ | | | | ________ | | HF --o | ----------- |Process| | | ->| Delay | | ->| |-------->| |->| Hybrid |--> ------- | |decorrelate| d_k(n) | | | Synth | QMF -> R ---->| |-------->| | -------- Synth ----------- |_______|--------------> ------------------------------------------------------------------------------ REQUIREMENTS ------------------------------------------------------------------------------ REFERENCES SC 29 Software Copyright Licencing Disclaimer: This software module was originally developed by Coding Technologies and edited by - in the course of development of the ISO/IEC 13818-7 and ISO/IEC 14496-3 standards for reference purposes and its performance may not have been optimized. This software module is an implementation of one or more tools as specified by the ISO/IEC 13818-7 and ISO/IEC 14496-3 standards. ISO/IEC gives users free license to this software module or modifications thereof for use in products claiming conformance to audiovisual and image-coding related ITU Recommendations and/or ISO/IEC International Standards. ISO/IEC gives users the same free license to this software module or modifications thereof for research purposes and further ISO/IEC standardisation. Those intending to use this software module in products are advised that its use may infringe existing patents. ISO/IEC have no liability for use of this software module or modifications thereof. Copyright is not released for products that do not conform to audiovisual and image-coding related ITU Recommendations and/or ISO/IEC International Standards. The original developer retains full right to modify and use the code for its own purpose, assign or donate the code to a third party and to inhibit third parties from using the code for products that do not conform to audiovisual and image-coding related ITU Recommendations and/or ISO/IEC International Standards. This copyright notice must be included in all copies or derivative works. Copyright (c) ISO/IEC 2003. ------------------------------------------------------------------------------ PSEUDO-CODE ------------------------------------------------------------------------------ */ /*---------------------------------------------------------------------------- ; INCLUDES ----------------------------------------------------------------------------*/ #ifdef AAC_PLUS #ifdef PARAMETRICSTEREO #include "pv_audio_type_defs.h" #include "ps_stereo_processing.h" #include "fxp_mul32.h" #include "ps_all_pass_filter_coeff.h" /*---------------------------------------------------------------------------- ; MACROS ; Define module specific macros here ----------------------------------------------------------------------------*/ #ifndef min #define min(a, b) ((a) < (b) ? (a) : (b)) #endif /*---------------------------------------------------------------------------- ; DEFINES ; Include all pre-processor statements here. Include conditional ; compile variables also. ----------------------------------------------------------------------------*/ /*---------------------------------------------------------------------------- ; LOCAL FUNCTION DEFINITIONS ; Function Prototype declaration ----------------------------------------------------------------------------*/ /*---------------------------------------------------------------------------- ; LOCAL STORE/BUFFER/POINTER DEFINITIONS ; Variable declaration - defined here and used outside this module ----------------------------------------------------------------------------*/ /*---------------------------------------------------------------------------- ; EXTERNAL FUNCTION REFERENCES ; Declare functions defined elsewhere and referenced in this module ----------------------------------------------------------------------------*/ /*---------------------------------------------------------------------------- ; EXTERNAL GLOBAL STORE/BUFFER/POINTER REFERENCES ; Declare variables used in this module but defined elsewhere ----------------------------------------------------------------------------*/ /*---------------------------------------------------------------------------- ; FUNCTION CODE ----------------------------------------------------------------------------*/ void ps_stereo_processing(STRUCT_PS_DEC *pms, Int32 *qmfLeftReal, Int32 *qmfLeftImag, Int32 *qmfRightReal, Int32 *qmfRightImag) { Int32 group; Int32 subband; Int32 maxSubband; Int32 usb; Char index; Int32 *hybrLeftReal; Int32 *hybrLeftImag; Int32 *hybrRightReal; Int32 *hybrRightImag; Int32 *ptr_hybrLeftReal; Int32 *ptr_hybrLeftImag; Int32 *ptr_hybrRightReal; Int32 *ptr_hybrRightImag; Int16 h11; Int16 h12; Int16 h21; Int16 h22; Int32 temp1; Int32 temp2; Int32 temp3; usb = pms->usb; /* * Complete Linear interpolation */ hybrLeftReal = pms->mHybridRealLeft; hybrLeftImag = pms->mHybridImagLeft; hybrRightReal = pms->mHybridRealRight; hybrRightImag = pms->mHybridImagRight; for (group = 0; group < SUBQMF_GROUPS; group++) /* SUBQMF_GROUPS == 10 */ { temp1 = pms->deltaH11[group]; temp2 = pms->deltaH12[group]; pms->H11[group] += temp1; h11 = (Int16)(pms->H11[group] >> 16); pms->H12[group] += temp2; h12 = (Int16)(pms->H12[group] >> 16); temp1 = pms->deltaH21[group]; temp2 = pms->deltaH22[group]; pms->H21[group] += temp1; h21 = (Int16)(pms->H21[group] >> 16); pms->H22[group] += temp2; h22 = (Int16)(pms->H22[group] >> 16); index = groupBorders[group]; /* * Reconstruction of Stereo sub-band signal * * l_k(n) = H11(k,n)*s_k(n) + H21(k,n)*d_k(n) * * r_k(n) = H12(k,n)*s_k(n) + H22(k,n)*d_k(n) */ ptr_hybrLeftReal = &hybrLeftReal[ index]; ptr_hybrRightReal = &hybrRightReal[ index]; temp1 = *(ptr_hybrLeftReal) << 1; temp2 = *(ptr_hybrRightReal) << 1; temp3 = fxp_mul32_by_16(temp1, h11); *(ptr_hybrLeftReal) = fxp_mac32_by_16(temp2, h21, temp3) << 1; temp3 = fxp_mul32_by_16(temp1, h12); *(ptr_hybrRightReal) = fxp_mac32_by_16(temp2, h22, temp3) << 1; ptr_hybrLeftImag = &hybrLeftImag[ index]; ptr_hybrRightImag = &hybrRightImag[ index]; temp1 = *(ptr_hybrLeftImag) << 1; temp2 = *(ptr_hybrRightImag) << 1; temp3 = fxp_mul32_by_16(temp1, h11); *(ptr_hybrLeftImag) = fxp_mac32_by_16(temp2, h21, temp3) << 1; temp3 = fxp_mul32_by_16(temp1, h12); *(ptr_hybrRightImag) = fxp_mac32_by_16(temp2, h22, temp3) << 1; } /* groups loop */ temp1 = pms->deltaH11[SUBQMF_GROUPS]; temp2 = pms->deltaH12[SUBQMF_GROUPS]; pms->H11[SUBQMF_GROUPS] += temp1; h11 = (Int16)(pms->H11[SUBQMF_GROUPS] >> 16); pms->H12[SUBQMF_GROUPS] += temp2; h12 = (Int16)(pms->H12[SUBQMF_GROUPS] >> 16); temp1 = pms->deltaH21[SUBQMF_GROUPS]; temp2 = pms->deltaH22[SUBQMF_GROUPS]; pms->H21[SUBQMF_GROUPS] += temp1; h21 = (Int16)(pms->H21[SUBQMF_GROUPS] >> 16); pms->H22[SUBQMF_GROUPS] += temp2; h22 = (Int16)(pms->H22[SUBQMF_GROUPS] >> 16); ptr_hybrLeftReal = &qmfLeftReal[ 3]; ptr_hybrRightReal = &qmfRightReal[ 3]; /* * Reconstruction of Stereo sub-band signal * * l_k(n) = H11(k,n)*s_k(n) + H21(k,n)*d_k(n) * * r_k(n) = H12(k,n)*s_k(n) + H22(k,n)*d_k(n) */ temp1 = *(ptr_hybrLeftReal) << 1; temp2 = *(ptr_hybrRightReal) << 1; temp3 = fxp_mul32_by_16(temp1, h11); *(ptr_hybrLeftReal) = fxp_mac32_by_16(temp2, h21, temp3) << 1; temp3 = fxp_mul32_by_16(temp1, h12); *(ptr_hybrRightReal) = fxp_mac32_by_16(temp2, h22, temp3) << 1; ptr_hybrLeftImag = &qmfLeftImag[ 3]; ptr_hybrRightImag = &qmfRightImag[ 3]; temp1 = *(ptr_hybrLeftImag) << 1; temp2 = *(ptr_hybrRightImag) << 1; temp3 = fxp_mul32_by_16(temp1, h11); *(ptr_hybrLeftImag) = fxp_mac32_by_16(temp2, h21, temp3) << 1; temp3 = fxp_mul32_by_16(temp1, h12); *(ptr_hybrRightImag) = fxp_mac32_by_16(temp2, h22, temp3) << 1; for (group = SUBQMF_GROUPS + 1; group < NO_IID_GROUPS; group++) /* 11 to NO_IID_GROUPS == 22 */ { temp1 = pms->deltaH11[group]; temp2 = pms->deltaH12[group]; pms->H11[group] += temp1; h11 = (Int16)(pms->H11[group] >> 16); pms->H12[group] += temp2; h12 = (Int16)(pms->H12[group] >> 16); temp1 = pms->deltaH21[group]; temp2 = pms->deltaH22[group]; pms->H21[group] += temp1; h21 = (Int16)(pms->H21[group] >> 16); pms->H22[group] += temp2; h22 = (Int16)(pms->H22[group] >> 16); index = groupBorders[group]; maxSubband = groupBorders[group + 1]; maxSubband = min(usb, maxSubband); /* * Reconstruction of Stereo sub-band signal * * l_k(n) = H11(k,n)*s_k(n) + H21(k,n)*d_k(n) * * r_k(n) = H12(k,n)*s_k(n) + H22(k,n)*d_k(n) */ ptr_hybrLeftReal = &qmfLeftReal[ index]; ptr_hybrRightReal = &qmfRightReal[ index]; for (subband = index; subband < maxSubband; subband++) { temp1 = *(ptr_hybrLeftReal) << 1; temp2 = *(ptr_hybrRightReal) << 1; temp3 = fxp_mul32_by_16(temp1, h11); *(ptr_hybrLeftReal++) = fxp_mac32_by_16(temp2, h21, temp3) << 1; temp3 = fxp_mul32_by_16(temp1, h12); *(ptr_hybrRightReal++) = fxp_mac32_by_16(temp2, h22, temp3) << 1; } ptr_hybrLeftImag = &qmfLeftImag[ index]; ptr_hybrRightImag = &qmfRightImag[ index]; for (subband = index; subband < maxSubband; subband++) { temp1 = *(ptr_hybrLeftImag) << 1; temp2 = *(ptr_hybrRightImag) << 1; temp3 = fxp_mul32_by_16(temp1, h11); *(ptr_hybrLeftImag++) = fxp_mac32_by_16(temp2, h21, temp3) << 1; temp3 = fxp_mul32_by_16(temp1, h12); *(ptr_hybrRightImag++) = fxp_mac32_by_16(temp2, h22, temp3) << 1; } /* subband loop */ } /* groups loop */ } /* END ps_stereo_processing */ #endif #endif