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/spreproc.c
35b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Functions: subframePreProc
36b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
37b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber     Date: 02/06/2002
38b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
39b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------
40b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber REVISION HISTORY
41b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
42b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: Updated template used to PV coding template.
43b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Eliminated unnecessary use of the sub() function.
44b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
45b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description:
46b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber              1. Replaced copy() and for-loop with more efficient memcpy().
47b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber              2. Eliminated unused include file copy.h.
48b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
49b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description:  Replaced OSCL mem type functions and eliminated include
50b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber               files that now are chosen by OSCL definitions
51b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
52b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description:  Replaced "int" and/or "char" with OSCL defined types.
53b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
54b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description:
55b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
56b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------
57b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber MODULE DESCRIPTION
58b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
59b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
60b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------
61b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber*/
62b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
63b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*----------------------------------------------------------------------------
64b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; INCLUDES
65b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/
66b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include <string.h>
67b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
68b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "spreproc.h"
69b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "typedef.h"
70b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "weight_a.h"
71b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "syn_filt.h"
72b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "residu.h"
73b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
74b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*----------------------------------------------------------------------------
75b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; MACROS
76b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; Define module specific macros here
77b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/
78b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
79b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*----------------------------------------------------------------------------
80b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; DEFINES
81b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; Include all pre-processor statements here. Include conditional
82b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; compile variables also.
83b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/
84b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
85b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*----------------------------------------------------------------------------
86b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; LOCAL FUNCTION DEFINITIONS
87b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; Function Prototype declaration
88b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/
89b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
90b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*----------------------------------------------------------------------------
91b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; LOCAL VARIABLE DEFINITIONS
92b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; Variable declaration - defined here and used outside this module
93b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/
94b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
95b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*
96b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------
97b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber FUNCTION NAME: subframePreProc
98b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------
99b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber INPUT AND OUTPUT DEFINITIONS
100b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
101b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Inputs:
102b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    mode        -- enum Mode          -- coder mode
103b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    gamma1      -- const Word16 array -- spectral exp. factor 1
104b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    gamma1_12k2 -- const Word16 array -- spectral exp. factor 1 for EFR
105b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    gamma2      -- const Word16 array -- spectral exp. factor 2
106b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    A           -- Pointer to Word16  -- A(z) unquantized for the 4 subframes
107b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Aq          -- Pointer to Word16  -- A(z)   quantized for the 4 subframes
108b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    speech      -- Pointer to Word16  -- speech segment
109b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    mem_err     -- Pointer to Word16  -- pointer to error signal
110b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    mem_w0      -- Pointer to Word16  -- memory of weighting filter
111b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    zero        -- Pointer to Word16  -- pointer to zero vector
112b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
113b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Outputs:
114b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    ai_zero -- Word16 array -- history of weighted synth. filter
115b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    exc     -- Word16 array -- long term prediction residual
116b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    h1      -- Word16 array -- impulse response
117b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    xn      -- Word16 array -- target vector for pitch search
118b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    res2    -- Word16 array -- long term prediction residual
119b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    error   -- Word16 array -- error of LPC synthesis filter
120b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
121b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Returns:
122b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Zero
123b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
124b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Global Variables Used:
125b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    None
126b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
127b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Local Variables Needed:
128b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    None
129b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
130b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------
131b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber FUNCTION DESCRIPTION
132b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
133b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
134b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------
135b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber REQUIREMENTS
136b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
137b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber None
138b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
139b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------
140b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber REFERENCES
141b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
142b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber spreproc.c, UMTS GSM AMR speech codec, R99 - Version 3.2.0, March 2, 2001
143b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
144b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------
145b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber PSEUDO-CODE
146b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
147b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
148b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------
149b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber RESOURCES USED [optional]
150b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
151b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber When the code is written for a specific target processor the
152b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber the resources used should be documented below.
153b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
154b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber HEAP MEMORY USED: x bytes
155b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
156b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber STACK MEMORY USED: x bytes
157b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
158b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber CLOCK CYCLES: (cycle count equation for this function) + (variable
159b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber                used to represent cycle count for each subroutine
160b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber                called)
161b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber     where: (cycle count variable) = cycle count for [subroutine
162b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber                                     name]
163b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
164b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------
165b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber CAUTION [optional]
166b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber [State any special notes, constraints or cautions for users of this function]
167b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
168b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------
169b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber*/
170b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
171b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Hubervoid subframePreProc(
172b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    enum Mode mode,            /* i  : coder mode                            */
173b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    const Word16 gamma1[],     /* i  : spectral exp. factor 1                */
174b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    const Word16 gamma1_12k2[],/* i  : spectral exp. factor 1 for EFR        */
175b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    const Word16 gamma2[],     /* i  : spectral exp. factor 2                */
176b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Word16 *A,                 /* i  : A(z) unquantized for the 4 subframes  */
177b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Word16 *Aq,                /* i  : A(z)   quantized for the 4 subframes  */
178b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Word16 *speech,            /* i  : speech segment                        */
179b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Word16 *mem_err,           /* i  : pointer to error signal               */
180b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Word16 *mem_w0,            /* i  : memory of weighting filter            */
181b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Word16 *zero,              /* i  : pointer to zero vector                */
182b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Word16 ai_zero[],          /* o  : history of weighted synth. filter     */
183b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Word16 exc[],              /* o  : long term prediction residual         */
184b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Word16 h1[],               /* o  : impulse response                      */
185b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Word16 xn[],               /* o  : target vector for pitch search        */
186b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Word16 res2[],             /* o  : long term prediction residual         */
187b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Word16 error[]             /* o  : error of LPC synthesis filter         */
188b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber)
189b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber{
190b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Word16 Ap1[MP1];              /* A(z) with spectral expansion         */
191b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Word16 Ap2[MP1];              /* A(z) with spectral expansion         */
192b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    const Word16 *g1;             /* Pointer to correct gammma1 vector    */
193b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
194b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    /* mode specific pointer to gamma1 values */
195b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    if (mode == MR122 || mode == MR102)
196b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    {
197b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber        g1 = gamma1_12k2;
198b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    }
199b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    else
200b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    {
201b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber        g1 = gamma1;
202b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    }
203b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
204b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    /* Find the weighted LPC coefficients for the weighting filter. */
205b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Weight_Ai(A, g1, Ap1);
206b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Weight_Ai(A, gamma2, Ap2);
207b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
208b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    memcpy(ai_zero, Ap1, (M + 1)*sizeof(Word16));
209b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
210b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
211b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Syn_filt(Aq, ai_zero, h1, L_SUBFR, zero, 0);
212b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Syn_filt(Ap2, h1, h1, L_SUBFR, zero, 0);
213b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
214b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    /*
215b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber     *
216b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber     *          Find the target vector for pitch search:
217b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber     *
218b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber     */
219b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
220b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    /* LPC residual */
221b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Residu(Aq, speech, res2, L_SUBFR);
222b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
223b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    memcpy(exc, res2, L_SUBFR*sizeof(Word16));
224b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
225b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Syn_filt(Aq, exc, error, L_SUBFR, mem_err, 0);
226b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
227b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Residu(Ap1, error, xn, L_SUBFR);
228b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
229b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    /* target signal xn[]*/
230b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    Syn_filt(Ap2, xn, xn, L_SUBFR, mem_w0, 0);
231b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
232b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber    return;
233b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber
234b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber}
235