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