/* ------------------------------------------------------------------ * 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. * ------------------------------------------------------------------- */ /* ------------------------------------------------------------------------------ PacketVideo Corp. MP3 Decoder Library Filename: pvmp3_polyphase_filter_window.cpp Date: 09/21/2007 ------------------------------------------------------------------------------ REVISION HISTORY Description: ------------------------------------------------------------------------------ INPUT AND OUTPUT DEFINITIONS Input int32 *synth_buffer, synthesis input buffer int16 *outPcm, generated output ( 32 values) int32 numChannels number of channels Returns int16 *outPcm ------------------------------------------------------------------------------ FUNCTION DESCRIPTION apply polyphase filter window Input 32 subband samples Calculate 64 values ------------------------------------------------------------------------------ REQUIREMENTS ------------------------------------------------------------------------------ REFERENCES [1] ISO MPEG Audio Subgroup Software Simulation Group (1996) ISO 13818-3 MPEG-2 Audio Decoder - Lower Sampling Frequency Extension ------------------------------------------------------------------------------ PSEUDO-CODE ------------------------------------------------------------------------------ */ #if ( !defined(PV_ARM_GCC_V5) && !defined(PV_ARM_GCC_V4) && !defined(PV_ARM_V5) && !defined(PV_ARM_V4) ) /*---------------------------------------------------------------------------- ; INCLUDES ----------------------------------------------------------------------------*/ #include "pvmp3_polyphase_filter_window.h" #include "pv_mp3dec_fxd_op.h" #include "pvmp3_dec_defs.h" #include "pvmp3_tables.h" /*---------------------------------------------------------------------------- ; MACROS ; Define module1 specific macros here ----------------------------------------------------------------------------*/ /*---------------------------------------------------------------------------- ; 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 module1 ----------------------------------------------------------------------------*/ /*---------------------------------------------------------------------------- ; EXTERNAL FUNCTION REFERENCES ; Declare functions defined elsewhere and referenced in this module_x ----------------------------------------------------------------------------*/ /*---------------------------------------------------------------------------- ; EXTERNAL GLOBAL STORE/BUFFER/POINTER REFERENCES ; Declare variables used in this module_x but defined elsewhere ----------------------------------------------------------------------------*/ /*---------------------------------------------------------------------------- ; FUNCTION CODE ----------------------------------------------------------------------------*/ void pvmp3_polyphase_filter_window(int32 *synth_buffer, int16 *outPcm, int32 numChannels) { int32 sum1; int32 sum2; const int32 *winPtr = pqmfSynthWin; int32 i; for (int16 j = 1; j < SUBBANDS_NUMBER / 2; j++) { sum1 = 0x00000020; sum2 = 0x00000020; for (i = (SUBBANDS_NUMBER >> 1); i < HAN_SIZE + (SUBBANDS_NUMBER >> 1); i += SUBBANDS_NUMBER << 4) { int32 *pt_1 = &synth_buffer[ i+j]; int32 *pt_2 = &synth_buffer[ i-j]; int32 temp1 = pt_1[ 0]; int32 temp3 = pt_2[ SUBBANDS_NUMBER*15 ]; int32 temp2 = pt_2[ SUBBANDS_NUMBER* 1 ]; int32 temp4 = pt_1[ SUBBANDS_NUMBER*14 ]; sum1 = fxp_mac32_Q32(sum1, temp1, winPtr[ 0]); sum2 = fxp_mac32_Q32(sum2, temp3, winPtr[ 0]); sum2 = fxp_mac32_Q32(sum2, temp1, winPtr[ 1]); sum1 = fxp_msb32_Q32(sum1, temp3, winPtr[ 1]); sum1 = fxp_mac32_Q32(sum1, temp2, winPtr[ 2]); sum2 = fxp_msb32_Q32(sum2, temp4, winPtr[ 2]); sum2 = fxp_mac32_Q32(sum2, temp2, winPtr[ 3]); sum1 = fxp_mac32_Q32(sum1, temp4, winPtr[ 3]); temp1 = pt_1[ SUBBANDS_NUMBER* 2]; temp3 = pt_2[ SUBBANDS_NUMBER*13]; temp2 = pt_2[ SUBBANDS_NUMBER* 3]; temp4 = pt_1[ SUBBANDS_NUMBER*12]; sum1 = fxp_mac32_Q32(sum1, temp1, winPtr[ 4]); sum2 = fxp_mac32_Q32(sum2, temp3, winPtr[ 4]); sum2 = fxp_mac32_Q32(sum2, temp1, winPtr[ 5]); sum1 = fxp_msb32_Q32(sum1, temp3, winPtr[ 5]); sum1 = fxp_mac32_Q32(sum1, temp2, winPtr[ 6]); sum2 = fxp_msb32_Q32(sum2, temp4, winPtr[ 6]); sum2 = fxp_mac32_Q32(sum2, temp2, winPtr[ 7]); sum1 = fxp_mac32_Q32(sum1, temp4, winPtr[ 7]); temp1 = pt_1[ SUBBANDS_NUMBER* 4 ]; temp3 = pt_2[ SUBBANDS_NUMBER*11 ]; temp2 = pt_2[ SUBBANDS_NUMBER* 5 ]; temp4 = pt_1[ SUBBANDS_NUMBER*10 ]; sum1 = fxp_mac32_Q32(sum1, temp1, winPtr[ 8]); sum2 = fxp_mac32_Q32(sum2, temp3, winPtr[ 8]); sum2 = fxp_mac32_Q32(sum2, temp1, winPtr[ 9]); sum1 = fxp_msb32_Q32(sum1, temp3, winPtr[ 9]); sum1 = fxp_mac32_Q32(sum1, temp2, winPtr[10]); sum2 = fxp_msb32_Q32(sum2, temp4, winPtr[10]); sum2 = fxp_mac32_Q32(sum2, temp2, winPtr[11]); sum1 = fxp_mac32_Q32(sum1, temp4, winPtr[11]); temp1 = pt_1[ SUBBANDS_NUMBER*6 ]; temp3 = pt_2[ SUBBANDS_NUMBER*9 ]; temp2 = pt_2[ SUBBANDS_NUMBER*7 ]; temp4 = pt_1[ SUBBANDS_NUMBER*8 ]; sum1 = fxp_mac32_Q32(sum1, temp1, winPtr[12]); sum2 = fxp_mac32_Q32(sum2, temp3, winPtr[12]); sum2 = fxp_mac32_Q32(sum2, temp1, winPtr[13]); sum1 = fxp_msb32_Q32(sum1, temp3, winPtr[13]); sum1 = fxp_mac32_Q32(sum1, temp2, winPtr[14]); sum2 = fxp_msb32_Q32(sum2, temp4, winPtr[14]); sum2 = fxp_mac32_Q32(sum2, temp2, winPtr[15]); sum1 = fxp_mac32_Q32(sum1, temp4, winPtr[15]); winPtr += 16; } int32 k = j << (numChannels - 1); outPcm[k] = saturate16(sum1 >> 6); outPcm[(numChannels<<5) - k] = saturate16(sum2 >> 6); } sum1 = 0x00000020; sum2 = 0x00000020; for (i = 16; i < HAN_SIZE + 16; i += (SUBBANDS_NUMBER << 2)) { int32 *pt_synth = &synth_buffer[i]; int32 temp1 = pt_synth[ 0 ]; int32 temp2 = pt_synth[ SUBBANDS_NUMBER ]; int32 temp3 = pt_synth[ SUBBANDS_NUMBER/2]; sum1 = fxp_mac32_Q32(sum1, temp1, winPtr[0]) ; sum1 = fxp_mac32_Q32(sum1, temp2, winPtr[1]) ; sum2 = fxp_mac32_Q32(sum2, temp3, winPtr[2]) ; temp1 = pt_synth[ SUBBANDS_NUMBER<<1 ]; temp2 = pt_synth[ 3*SUBBANDS_NUMBER ]; temp3 = pt_synth[ SUBBANDS_NUMBER*5/2]; sum1 = fxp_mac32_Q32(sum1, temp1, winPtr[3]) ; sum1 = fxp_mac32_Q32(sum1, temp2, winPtr[4]) ; sum2 = fxp_mac32_Q32(sum2, temp3, winPtr[5]) ; winPtr += 6; } outPcm[0] = saturate16(sum1 >> 6); outPcm[(SUBBANDS_NUMBER/2)<<(numChannels-1)] = saturate16(sum2 >> 6); } #endif // If not assembly