1b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/* ------------------------------------------------------------------
2b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * Copyright (C) 1998-2009 PacketVideo
3b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *
4b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * Licensed under the Apache License, Version 2.0 (the "License");
5b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * you may not use this file except in compliance with the License.
6b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * You may obtain a copy of the License at
7b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *
8b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *      http://www.apache.org/licenses/LICENSE-2.0
9b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *
10b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * Unless required by applicable law or agreed to in writing, software
11b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * distributed under the License is distributed on an "AS IS" BASIS,
12b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either
13b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * express or implied.
14b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * See the License for the specific language governing permissions
15b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * and limitations under the License.
16b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * -------------------------------------------------------------------
17b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber */
18b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/****************************************************************************************
19b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas HuberPortions of this file are derived from the following 3GPP standard:
20b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
21b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    3GPP TS 26.073
22b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    ANSI-C code for the Adaptive Multi-Rate (AMR) speech codec
23b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Available from http://www.3gpp.org
24b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
25b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber(C) 2004, 3GPP Organizational Partners (ARIB, ATIS, CCSA, ETSI, TTA, TTC)
26b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas HuberPermission to distribute, modify and use this file under the standard license
27b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huberterms listed above has been obtained from the copyright holder.
28b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber****************************************************************************************/
29b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*
30b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------
31b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
32b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
33b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
34b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Pathname: ./audio/gsm-amr/c/src/lag_wind.c
35b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
36b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber     Date: 01/31/2002
37b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
38b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------
39b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber REVISION HISTORY
40b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
41b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description:
42b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber              1. Eliminated unused include files.
43b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber              2. Replaced array addressing by pointers
44b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber              3. Eliminated l_extract() function call
45b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
46b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description:  Added casting to eliminate warnings
47b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
48b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description:  Replaced "int" and/or "char" with OSCL defined types.
49b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
50b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description:
51b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
52b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------
53b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber*/
54b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
55b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*----------------------------------------------------------------------------
56b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; INCLUDES
57b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/
58b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "lag_wind.h"
59b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "lag_wind_tab.h"
60b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "basic_op.h"
61b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
62b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*----------------------------------------------------------------------------
63b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; MACROS
64b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; Define module specific macros here
65b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/
66b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
67b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
68b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*----------------------------------------------------------------------------
69b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; DEFINES
70b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; Include all pre-processor statements here. Include conditional
71b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; compile variables also.
72b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/
73b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
74b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
75b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*----------------------------------------------------------------------------
76b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; LOCAL FUNCTION DEFINITIONS
77b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; Function Prototype declaration
78b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/
79b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
80b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
81b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*----------------------------------------------------------------------------
82b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; LOCAL STORE/BUFFER/POINTER DEFINITIONS
83b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; Variable declaration - defined here and used outside this module
84b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/
85b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
86b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
87b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*
88b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------
89b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber FUNCTION NAME: lag_wind
90b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------
91b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber INPUT AND OUTPUT DEFINITIONS
92b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
93b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Inputs:
94b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    m = LPC order of type Word16
95b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    r_h[] = pointer to autocorrelations (msb) of type Word16
96b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    r_l[] = pointer to autocorrelations (lsb) of type Word16
97b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    pOverflow = pointer to overflow flag
98b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
99b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Outputs:
100b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    None
101b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
102b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Returns:
103b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    None
104b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
105b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Global Variables Used:
106b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    None.
107b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
108b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Local Variables Needed:
109b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    None.
110b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
111b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------
112b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber FUNCTION DESCRIPTION
113b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
114b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber      File             : lag_wind.c
115b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber      Purpose          : Lag windowing of autocorrelations.
116b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
117b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    FUNCTION:  Lag_window()
118b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
119b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    PURPOSE:  Lag windowing of autocorrelations.
120b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
121b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    DESCRIPTION:
122b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber          r[i] = r[i]*lag_wind[i],   i=1,...,10
123b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
124b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber      r[i] and lag_wind[i] are in special double precision format.
125b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber      See "oper_32b.c" for the format.
126b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
127b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------
128b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber REQUIREMENTS
129b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
130b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber None.
131b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
132b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------
133b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber REFERENCES
134b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
135b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber lag_wind.c, UMTS GSM AMR speech codec, R99 - Version 3.2.0, March 2, 2001
136b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
137b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------
138b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber PSEUDO-CODE
139b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
140b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Word16 i;
141b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Word32 x;
142b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
143b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    for (i = 1; i <= m; i++)
144b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    {
145b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber        x = Mpy_32 (r_h[i], r_l[i], lag_h[i - 1], lag_l[i - 1], pOverflow);
146b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber        L_Extract (x, &r_h[i], &r_l[i], pOverflow);
147b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    }
148b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    return;
149b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
150b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------
151b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber RESOURCES USED [optional]
152b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
153b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber When the code is written for a specific target processor the
154b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber the resources used should be documented below.
155b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
156b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber HEAP MEMORY USED: x bytes
157b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
158b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber STACK MEMORY USED: x bytes
159b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
160b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber CLOCK CYCLES: (cycle count equation for this function) + (variable
161b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber                used to represent cycle count for each subroutine
162b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber                called)
163b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber     where: (cycle count variable) = cycle count for [subroutine
164b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber                                     name]
165b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
166b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------
167b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber CAUTION [optional]
168b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber [State any special notes, constraints or cautions for users of this function]
169b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
170b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------
171b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber*/
172b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Hubervoid Lag_window(
173b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Word16 m,           /* (i)     : LPC order                        */
174b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Word16 r_h[],       /* (i/o)   : Autocorrelations  (msb)          */
175b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Word16 r_l[],       /* (i/o)   : Autocorrelations  (lsb)          */
176b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Flag   *pOverflow
177b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber)
178b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber{
179b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Word16 i;
180b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Word32 x;
181b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    const Word16 *p_lag_h = &lag_h[0];
182b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    const Word16 *p_lag_l = &lag_l[0];
183b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Word16 *p_r_h = &r_h[1];
184b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Word16 *p_r_l = &r_l[1];
185b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
186b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    for (i = m; i != 0 ; i--)
187b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    {
188b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber        x = Mpy_32(*(p_r_h), *(p_r_l), *(p_lag_h++), *(p_lag_l++), pOverflow);
189b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber        *(p_r_h) = (Word16)(x >> 16);
190b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber        *(p_r_l++) = (x >> 1) - (*(p_r_h++) << 15);
191b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    }
192b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
193b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    return;
194b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber}
195b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
196