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_9.cpp 25d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 26d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber Date: 09/21/2007 27d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 28d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber------------------------------------------------------------------------------ 29d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber REVISION HISTORY 30d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 31d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 32d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber Description: 33d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 34d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber------------------------------------------------------------------------------ 35d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber INPUT AND OUTPUT DEFINITIONS 36d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 37d9618f23226f46c752e56f712bc4b505117d8b4bAndreas HuberInput 38d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 vec[] vector of 9 32-bit integers 39d9618f23226f46c752e56f712bc4b505117d8b4bAndreas HuberReturns 40d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 vec[] dct computation in-place 41d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 42d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 43d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber------------------------------------------------------------------------------ 44d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber FUNCTION DESCRIPTION 45d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 46d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber Returns the dct of length 9 of the input vector 47d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 48d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber------------------------------------------------------------------------------ 49d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber REQUIREMENTS 50d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 51d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 52d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber------------------------------------------------------------------------------ 53d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber REFERENCES 54d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 55d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber------------------------------------------------------------------------------ 56d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber PSEUDO-CODE 57d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 58d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber------------------------------------------------------------------------------ 59d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber*/ 60d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 61d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber#if ( !defined(PV_ARM_GCC_V5) && !defined(PV_ARM_GCC_V4) && !defined(PV_ARM_V5) && !defined(PV_ARM_V4) ) 62d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber/*---------------------------------------------------------------------------- 63d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; INCLUDES 64d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber----------------------------------------------------------------------------*/ 65d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber#include "pvmp3_audio_type_defs.h" 66d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber#include "pv_mp3dec_fxd_op.h" 67d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber#include "pvmp3_mdct_18.h" 68d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 69d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber/*---------------------------------------------------------------------------- 70d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; MACROS 71d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; Define module specific macros here 72d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber----------------------------------------------------------------------------*/ 73d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 74d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 75d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber/*---------------------------------------------------------------------------- 76d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; DEFINES 77d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; Include all pre-processor statements here. Include conditional 78d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; compile variables also. 79d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber----------------------------------------------------------------------------*/ 80d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber#define Qfmt31(a) (int32)(a*(0x7FFFFFFF)) 81d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 82d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber#define cos_pi_9 Qfmt31( 0.93969262078591f) 83d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber#define cos_2pi_9 Qfmt31( 0.76604444311898f) 84d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber#define cos_4pi_9 Qfmt31( 0.17364817766693f) 85d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber#define cos_5pi_9 Qfmt31(-0.17364817766693f) 86d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber#define cos_7pi_9 Qfmt31(-0.76604444311898f) 87d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber#define cos_8pi_9 Qfmt31(-0.93969262078591f) 88d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber#define cos_pi_6 Qfmt31( 0.86602540378444f) 89d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber#define cos_5pi_6 Qfmt31(-0.86602540378444f) 90d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber#define cos_5pi_18 Qfmt31( 0.64278760968654f) 91d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber#define cos_7pi_18 Qfmt31( 0.34202014332567f) 92d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber#define cos_11pi_18 Qfmt31(-0.34202014332567f) 93d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber#define cos_13pi_18 Qfmt31(-0.64278760968654f) 94d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber#define cos_17pi_18 Qfmt31(-0.98480775301221f) 95d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 96d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber/*---------------------------------------------------------------------------- 97d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; LOCAL FUNCTION DEFINITIONS 98d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; Function Prototype declaration 99d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber----------------------------------------------------------------------------*/ 100d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 101d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber/*---------------------------------------------------------------------------- 102d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; LOCAL STORE/BUFFER/POINTER DEFINITIONS 103d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; Variable declaration - defined here and used outside this module 104d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber----------------------------------------------------------------------------*/ 105d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 106d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber/*---------------------------------------------------------------------------- 107d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; EXTERNAL FUNCTION REFERENCES 108d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; Declare functions defined elsewhere and referenced in this module 109d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber----------------------------------------------------------------------------*/ 110d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 111d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber/*---------------------------------------------------------------------------- 112d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; EXTERNAL GLOBAL STORE/BUFFER/POINTER REFERENCES 113d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; Declare variables used in this module but defined elsewhere 114d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber----------------------------------------------------------------------------*/ 115d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 116d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber/*---------------------------------------------------------------------------- 117d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber; FUNCTION CODE 118d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber----------------------------------------------------------------------------*/ 119d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 120d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Hubervoid pvmp3_dct_9(int32 vec[]) 121d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber{ 122d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 123d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber /* split input vector */ 124d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 125d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp0 = vec[8] + vec[0]; 126d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp8 = vec[8] - vec[0]; 127d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp1 = vec[7] + vec[1]; 128d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp7 = vec[7] - vec[1]; 129d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp2 = vec[6] + vec[2]; 130d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp6 = vec[6] - vec[2]; 131d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp3 = vec[5] + vec[3]; 132d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber int32 tmp5 = vec[5] - vec[3]; 133d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 134d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[0] = (tmp0 + tmp2 + tmp3) + (tmp1 + vec[4]); 135d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[6] = ((tmp0 + tmp2 + tmp3) >> 1) - (tmp1 + vec[4]); 136d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[2] = (tmp1 >> 1) - vec[4]; 137d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[4] = -vec[2]; 138d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[8] = -vec[2]; 139d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[4] = fxp_mac32_Q32(vec[4], tmp0 << 1, cos_2pi_9); 140d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[8] = fxp_mac32_Q32(vec[8], tmp0 << 1, cos_4pi_9); 141d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[2] = fxp_mac32_Q32(vec[2], tmp0 << 1, cos_pi_9); 142d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[2] = fxp_mac32_Q32(vec[2], tmp2 << 1, cos_5pi_9); 143d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[4] = fxp_mac32_Q32(vec[4], tmp2 << 1, cos_8pi_9); 144d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[8] = fxp_mac32_Q32(vec[8], tmp2 << 1, cos_2pi_9); 145d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[8] = fxp_mac32_Q32(vec[8], tmp3 << 1, cos_8pi_9); 146d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[4] = fxp_mac32_Q32(vec[4], tmp3 << 1, cos_4pi_9); 147d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[2] = fxp_mac32_Q32(vec[2], tmp3 << 1, cos_7pi_9); 148d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 149d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[1] = fxp_mul32_Q32(tmp5 << 1, cos_11pi_18); 150d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[1] = fxp_mac32_Q32(vec[1], tmp6 << 1, cos_13pi_18); 151d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[1] = fxp_mac32_Q32(vec[1], tmp7 << 1, cos_5pi_6); 152d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[1] = fxp_mac32_Q32(vec[1], tmp8 << 1, cos_17pi_18); 153d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[3] = fxp_mul32_Q32((tmp5 + tmp6 - tmp8) << 1, cos_pi_6); 154d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[5] = fxp_mul32_Q32(tmp5 << 1, cos_17pi_18); 155d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[5] = fxp_mac32_Q32(vec[5], tmp6 << 1, cos_7pi_18); 156d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[5] = fxp_mac32_Q32(vec[5], tmp7 << 1, cos_pi_6); 157d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[5] = fxp_mac32_Q32(vec[5], tmp8 << 1, cos_13pi_18); 158d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[7] = fxp_mul32_Q32(tmp5 << 1, cos_5pi_18); 159d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[7] = fxp_mac32_Q32(vec[7], tmp6 << 1, cos_17pi_18); 160d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[7] = fxp_mac32_Q32(vec[7], tmp7 << 1, cos_pi_6); 161d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber vec[7] = fxp_mac32_Q32(vec[7], tmp8 << 1, cos_11pi_18); 162d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 163d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber} 164d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 165d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 166d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber 167d9618f23226f46c752e56f712bc4b505117d8b4bAndreas Huber#endif // If not assembly 168