1d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber/* ------------------------------------------------------------------ 2d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber * Copyright (C) 1998-2009 PacketVideo 3d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber * 4d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber * Licensed under the Apache License, Version 2.0 (the "License"); 5d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber * you may not use this file except in compliance with the License. 6d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber * You may obtain a copy of the License at 7d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber * 8d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber * http://www.apache.org/licenses/LICENSE-2.0 9d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber * 10d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber * Unless required by applicable law or agreed to in writing, software 11d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber * distributed under the License is distributed on an "AS IS" BASIS, 12d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either 13d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber * express or implied. 14d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber * See the License for the specific language governing permissions 15d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber * and limitations under the License. 16d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber * ------------------------------------------------------------------- 17d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber */ 18d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber/* 19d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber------------------------------------------------------------------------------ 20d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 21d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber PacketVideo Corp. 22d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber MP3 Decoder Library 23d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 24d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber Filename: pvmp3_dct_16.cpp 25d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 26d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber Functions: 27d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber dct_16 28d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber pv_merge_in_place_N32 29d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber pv_split 30d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 31d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber Date: 09/21/2007 32d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 33d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber------------------------------------------------------------------------------ 34d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber REVISION HISTORY 35d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 36d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 37d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber Description: 38d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 39d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber------------------------------------------------------------------------------ 40d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber INPUT AND OUTPUT DEFINITIONS 41d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 42d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber dct_16 43d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 44d9618f23226f46c752e56f712bc4b505117d8b4bAndreas HuberInput 45d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 vec[], input vector length 16 46d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber Int flag processing direction: forward (1), backward ( 0) 47d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber Returns 48d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 49d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 vec[], dct length 16 50d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 51d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber------------------------------------------------------------------------------ 52d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber INPUT AND OUTPUT DEFINITIONS 53d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 54d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber pv_merge_in_place_N32 55d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 56d9618f23226f46c752e56f712bc4b505117d8b4bAndreas HuberInput 57d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 vec[], input vector length 16 58d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 59d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber Returns 60d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 61d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 vec[], merged output of two dct 16 to create a dct 32 62d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 63d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber------------------------------------------------------------------------------ 64d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber INPUT AND OUTPUT DEFINITIONS 65d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 66d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber pv_split 67d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 68d9618f23226f46c752e56f712bc4b505117d8b4bAndreas HuberInput 69d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 vec[], input vector length 16 70d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 71d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber Returns 72d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 73d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 vec[], splitted even/odd and pre processing rotation 74d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 75d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber------------------------------------------------------------------------------ 76d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber FUNCTION DESCRIPTION 77d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 78d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber dct 16 and tools to assemble a dct32 output 79d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 80d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber------------------------------------------------------------------------------ 81d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber REQUIREMENTS 82d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 83d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 84d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber------------------------------------------------------------------------------ 85d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber REFERENCES 86d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 87d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber------------------------------------------------------------------------------ 88d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber PSEUDO-CODE 89d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 90d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber------------------------------------------------------------------------------ 91d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber*/ 92d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 93d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber#if ( !defined(PV_ARM_GCC_V5) && !defined(PV_ARM_GCC_V4) ) 94d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber/*---------------------------------------------------------------------------- 95d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; INCLUDES 96d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber----------------------------------------------------------------------------*/ 97d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 98d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber#include "pvmp3_dct_16.h" 99d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber#include "pv_mp3dec_fxd_op.h" 100d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 101d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber/*---------------------------------------------------------------------------- 102d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; MACROS 103d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; Define module specific macros here 104d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber----------------------------------------------------------------------------*/ 105d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 106d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 107d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber/*---------------------------------------------------------------------------- 108d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; DEFINES 109d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; Include all pre-processor statements here. Include conditional 110d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; compile variables also. 111d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber----------------------------------------------------------------------------*/ 112d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber#define Qfmt(a) (int32)(a*((int32)1<<27)) 113d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 114d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber/*---------------------------------------------------------------------------- 115d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; LOCAL FUNCTION DEFINITIONS 116d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; Function Prototype declaration 117d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber----------------------------------------------------------------------------*/ 118d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 119d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber/*---------------------------------------------------------------------------- 120d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; LOCAL STORE/BUFFER/POINTER DEFINITIONS 121d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; Variable declaration - defined here and used outside this module 122d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber----------------------------------------------------------------------------*/ 123d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huberconst int32 CosTable_dct32[16] = 124d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber{ 125d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber Qfmt_31(0.50060299823520F) , Qfmt_31(0.50547095989754F) , 126d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber Qfmt_31(0.51544730992262F) , Qfmt_31(0.53104259108978F) , 127d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber Qfmt_31(0.55310389603444F) , Qfmt_31(0.58293496820613F) , 128d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber Qfmt_31(0.62250412303566F) , Qfmt_31(0.67480834145501F) , 129d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber Qfmt_31(0.74453627100230F) , Qfmt_31(0.83934964541553F) , 130d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 131d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber Qfmt(0.97256823786196F) , Qfmt(1.16943993343288F) , 132d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber Qfmt(1.48416461631417F) , Qfmt(2.05778100995341F) , 133d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber Qfmt(3.40760841846872F) , Qfmt(10.19000812354803F) 134d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber}; 135d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 136d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 137d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber/*---------------------------------------------------------------------------- 138d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; EXTERNAL FUNCTION REFERENCES 139d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; Declare functions defined elsewhere and referenced in this module 140d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber----------------------------------------------------------------------------*/ 141d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 142d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber/*---------------------------------------------------------------------------- 143d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; EXTERNAL GLOBAL STORE/BUFFER/POINTER REFERENCES 144d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; Declare variables used in this module but defined elsewhere 145d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber----------------------------------------------------------------------------*/ 146d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 147d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 148d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber/*---------------------------------------------------------------------------- 149d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; FUNCTION CODE 150d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber----------------------------------------------------------------------------*/ 151d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 152d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Hubervoid pvmp3_dct_16(int32 vec[], int32 flag) 153d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber{ 154d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp0; 155d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp1; 156d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp2; 157d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp3; 158d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp4; 159d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp5; 160d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp6; 161d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp7; 162d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp_o0; 163d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp_o1; 164d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp_o2; 165d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp_o3; 166d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp_o4; 167d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp_o5; 168d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp_o6; 169d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp_o7; 170d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 itmp_e0; 171d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 itmp_e1; 172d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 itmp_e2; 173d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 174d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber /* split input vector */ 175d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 176d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp_o0 = fxp_mul32_Q32((vec[ 0] - vec[15]), Qfmt_31(0.50241928618816F)); 177d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp0 = vec[ 0] + vec[15]; 178d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 179d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp_o7 = fxp_mul32_Q32((vec[ 7] - vec[ 8]) << 3, Qfmt_31(0.63764357733614F)); 180d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp7 = vec[ 7] + vec[ 8]; 181d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 182d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber itmp_e0 = fxp_mul32_Q32((tmp0 - tmp7), Qfmt_31(0.50979557910416F)); 183d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp7 = (tmp0 + tmp7); 184d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 185d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp_o1 = fxp_mul32_Q32((vec[ 1] - vec[14]), Qfmt_31(0.52249861493969F)); 186d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp1 = vec[ 1] + vec[14]; 187d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 188d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp_o6 = fxp_mul32_Q32((vec[ 6] - vec[ 9]) << 1, Qfmt_31(0.86122354911916F)); 189d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp6 = vec[ 6] + vec[ 9]; 190d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 191d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 192d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 193d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber itmp_e1 = (tmp1 + tmp6); 194d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp6 = fxp_mul32_Q32((tmp1 - tmp6), Qfmt_31(0.60134488693505F)); 195d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 196d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 197d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 198d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp_o2 = fxp_mul32_Q32((vec[ 2] - vec[13]), Qfmt_31(0.56694403481636F)); 199d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp2 = vec[ 2] + vec[13]; 200d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp_o5 = fxp_mul32_Q32((vec[ 5] - vec[10]) << 1, Qfmt_31(0.53033884299517F)); 201d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp5 = vec[ 5] + vec[10]; 202d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 203d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber itmp_e2 = (tmp2 + tmp5); 204d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp5 = fxp_mul32_Q32((tmp2 - tmp5), Qfmt_31(0.89997622313642F)); 205d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 206d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp_o3 = fxp_mul32_Q32((vec[ 3] - vec[12]), Qfmt_31(0.64682178335999F)); 207d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp3 = vec[ 3] + vec[12]; 208d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp_o4 = fxp_mul32_Q32((vec[ 4] - vec[11]), Qfmt_31(0.78815462345125F)); 209d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp4 = vec[ 4] + vec[11]; 210d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 211d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp1 = (tmp3 + tmp4); 212d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp4 = fxp_mul32_Q32((tmp3 - tmp4) << 2, Qfmt_31(0.64072886193538F)); 213d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 214d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber /* split even part of tmp_e */ 215d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 216d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp0 = (tmp7 + tmp1); 217d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp1 = fxp_mul32_Q32((tmp7 - tmp1), Qfmt_31(0.54119610014620F)); 218d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 219d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp3 = fxp_mul32_Q32((itmp_e1 - itmp_e2) << 1, Qfmt_31(0.65328148243819F)); 220d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp7 = (itmp_e1 + itmp_e2); 221d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 222d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[ 0] = (tmp0 + tmp7) >> 1; 223d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[ 8] = fxp_mul32_Q32((tmp0 - tmp7), Qfmt_31(0.70710678118655F)); 224d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp0 = fxp_mul32_Q32((tmp1 - tmp3) << 1, Qfmt_31(0.70710678118655F)); 225d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[ 4] = tmp1 + tmp3 + tmp0; 226d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[12] = tmp0; 227d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 228d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber /* split odd part of tmp_e */ 229d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 230d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp1 = fxp_mul32_Q32((itmp_e0 - tmp4) << 1, Qfmt_31(0.54119610014620F)); 231d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp7 = itmp_e0 + tmp4; 232d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 233d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp3 = fxp_mul32_Q32((tmp6 - tmp5) << 2, Qfmt_31(0.65328148243819F)); 234d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp6 += tmp5; 235d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 236d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp4 = fxp_mul32_Q32((tmp7 - tmp6) << 1, Qfmt_31(0.70710678118655F)); 237d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp6 += tmp7; 238d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp7 = fxp_mul32_Q32((tmp1 - tmp3) << 1, Qfmt_31(0.70710678118655F)); 239d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 240d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp1 += tmp3 + tmp7; 241d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[ 2] = tmp1 + tmp6; 242d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[ 6] = tmp1 + tmp4; 243d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[10] = tmp7 + tmp4; 244d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[14] = tmp7; 245d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 246d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 247d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber // dct8; 248d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 249d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp1 = fxp_mul32_Q32((tmp_o0 - tmp_o7) << 1, Qfmt_31(0.50979557910416F)); 250d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp7 = tmp_o0 + tmp_o7; 251d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 252d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp6 = tmp_o1 + tmp_o6; 253d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp_o1 = fxp_mul32_Q32((tmp_o1 - tmp_o6) << 1, Qfmt_31(0.60134488693505F)); 254d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 255d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp5 = tmp_o2 + tmp_o5; 256d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp_o5 = fxp_mul32_Q32((tmp_o2 - tmp_o5) << 1, Qfmt_31(0.89997622313642F)); 257d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 258d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp0 = fxp_mul32_Q32((tmp_o3 - tmp_o4) << 3, Qfmt_31(0.6407288619354F)); 259d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp4 = tmp_o3 + tmp_o4; 260d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 261d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber if (!flag) 262d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber { 263d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp7 = -tmp7; 264d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp1 = -tmp1; 265d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp6 = -tmp6; 266d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp_o1 = -tmp_o1; 267d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp5 = -tmp5; 268d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp_o5 = -tmp_o5; 269d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp4 = -tmp4; 270d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp0 = -tmp0; 271d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber } 272d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 273d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 274d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp2 = fxp_mul32_Q32((tmp1 - tmp0) << 1, Qfmt_31(0.54119610014620F)); 275d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp0 += tmp1; 276d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp1 = fxp_mul32_Q32((tmp7 - tmp4) << 1, Qfmt_31(0.54119610014620F)); 277d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp7 += tmp4; 278d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp4 = fxp_mul32_Q32((tmp6 - tmp5) << 2, Qfmt_31(0.65328148243819F)); 279d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp6 += tmp5; 280d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp5 = fxp_mul32_Q32((tmp_o1 - tmp_o5) << 2, Qfmt_31(0.65328148243819F)); 281d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp_o1 += tmp_o5; 282d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 283d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 284d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[13] = fxp_mul32_Q32((tmp1 - tmp4) << 1, Qfmt_31(0.70710678118655F)); 285d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[ 5] = tmp1 + tmp4 + vec[13]; 286d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 287d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[ 9] = fxp_mul32_Q32((tmp7 - tmp6) << 1, Qfmt_31(0.70710678118655F)); 288d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[ 1] = tmp7 + tmp6; 289d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 290d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp4 = fxp_mul32_Q32((tmp0 - tmp_o1) << 1, Qfmt_31(0.70710678118655F)); 291d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp0 += tmp_o1; 292d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp6 = fxp_mul32_Q32((tmp2 - tmp5) << 1, Qfmt_31(0.70710678118655F)); 293d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp2 += tmp5 + tmp6; 294d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp0 += tmp2; 295d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 296d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[ 1] += tmp0; 297d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[ 3] = tmp0 + vec[ 5]; 298d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp2 += tmp4; 299d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[ 5] = tmp2 + vec[ 5]; 300d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[ 7] = tmp2 + vec[ 9]; 301d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp4 += tmp6; 302d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[ 9] = tmp4 + vec[ 9]; 303d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[11] = tmp4 + vec[13]; 304d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[13] = tmp6 + vec[13]; 305d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[15] = tmp6; 306d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 307d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber} 308d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber/*---------------------------------------------------------------------------- 309d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; FUNCTION CODE 310d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber----------------------------------------------------------------------------*/ 311d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Hubervoid pvmp3_merge_in_place_N32(int32 vec[]) 312d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber{ 313d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 314d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 315d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 temp0; 316d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 temp1; 317d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 temp2; 318d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 temp3; 319d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 320d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber temp0 = vec[14]; 321d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[14] = vec[ 7]; 322d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber temp1 = vec[12]; 323d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[12] = vec[ 6]; 324d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber temp2 = vec[10]; 325d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[10] = vec[ 5]; 326d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber temp3 = vec[ 8]; 327d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[ 8] = vec[ 4]; 328d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[ 6] = vec[ 3]; 329d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[ 4] = vec[ 2]; 330d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[ 2] = vec[ 1]; 331d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 332d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[ 1] = (vec[16] + vec[17]); 333d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[16] = temp3; 334d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[ 3] = (vec[18] + vec[17]); 335d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[ 5] = (vec[19] + vec[18]); 336d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[18] = vec[9]; 337d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 338d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[ 7] = (vec[20] + vec[19]); 339d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[ 9] = (vec[21] + vec[20]); 340d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[20] = temp2; 341d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber temp2 = vec[13]; 342d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber temp3 = vec[11]; 343d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[11] = (vec[22] + vec[21]); 344d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[13] = (vec[23] + vec[22]); 345d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[22] = temp3; 346d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber temp3 = vec[15]; 347d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 348d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[15] = (vec[24] + vec[23]); 349d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[17] = (vec[25] + vec[24]); 350d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[19] = (vec[26] + vec[25]); 351d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[21] = (vec[27] + vec[26]); 352d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[23] = (vec[28] + vec[27]); 353d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[24] = temp1; 354d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[25] = (vec[29] + vec[28]); 355d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[26] = temp2; 356d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[27] = (vec[30] + vec[29]); 357d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[28] = temp0; 358d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[29] = (vec[30] + vec[31]); 359d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[30] = temp3; 360d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber} 361d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 362d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 363d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber/*---------------------------------------------------------------------------- 364d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; FUNCTION CODE 365d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber----------------------------------------------------------------------------*/ 366d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 367d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 368d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 369d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Hubervoid pvmp3_split(int32 *vect) 370d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber{ 371d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 372d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 i; 373d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber const int32 *pt_cosTerms = &CosTable_dct32[15]; 374d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 *pt_vect = vect; 375d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 *pt_vect_2 = pt_vect - 1; 376d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 377d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber for (i = 3; i != 0; i--) 378d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber { 379d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp2 = *(pt_vect); 380d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp1 = *(pt_vect_2); 381d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 cosx = *(pt_cosTerms--); 382d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber *(pt_vect_2--) = (tmp1 + tmp2); 383d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber *(pt_vect++) = fxp_mul32_Q27((tmp1 - tmp2), cosx); 384d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 385d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp2 = *(pt_vect); 386d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp1 = *(pt_vect_2); 387d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber cosx = *(pt_cosTerms--); 388d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber *(pt_vect_2--) = (tmp1 + tmp2); 389d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber *(pt_vect++) = fxp_mul32_Q27((tmp1 - tmp2), cosx); 390d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 391d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber } 392d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 393d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber for (i = 5; i != 0; i--) 394d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber { 395d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp2 = *(pt_vect); 396d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp1 = *(pt_vect_2); 397d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 cosx = *(pt_cosTerms--); 398d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber *(pt_vect_2--) = (tmp1 + tmp2); 399d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber *(pt_vect++) = fxp_mul32_Q32((tmp1 - tmp2) << 1, cosx); 400d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 401d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp2 = *(pt_vect); 402d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber tmp1 = *(pt_vect_2); 403d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber cosx = *(pt_cosTerms--); 404d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber *(pt_vect_2--) = (tmp1 + tmp2); 405d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber *(pt_vect++) = fxp_mul32_Q32((tmp1 - tmp2) << 1, cosx); 406d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber } 407d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 408d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber} 409d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 410d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber#endif 411