spstproc.cpp revision b841f14f8e51f2365945281fbfa54ef6a1b1b5a6
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/spstproc.c 35b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Functions: subframePostProc 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 Changed to accept the pOverflow flag for EPOC compatibility. 44b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 45b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: 46b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1. Eliminated unused include files. 47b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 2. Replaced array addressing by pointers 48b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 3. Eliminated math operations that unnecessary checked for 49b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber saturation 50b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 4. Replaced loop counter with decrement loops 51b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 52b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: Added casting to eliminate warnings 53b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 54b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: Replaced "int" and/or "char" with OSCL defined types. 55b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 56b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: 57b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 58b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 59b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber MODULE DESCRIPTION 60b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 61b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Subframe post processing 62b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 63b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber*/ 64b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 65b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*---------------------------------------------------------------------------- 66b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; INCLUDES 67b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/ 68b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "spstproc.h" 69b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "syn_filt.h" 70b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "cnst.h" 71b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 72b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*---------------------------------------------------------------------------- 73b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; MACROS 74b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; Define module specific macros here 75b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/ 76b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 77b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*---------------------------------------------------------------------------- 78b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; DEFINES 79b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; Include all pre-processor statements here. Include conditional 80b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; compile variables also. 81b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/ 82b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 83b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*---------------------------------------------------------------------------- 84b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; LOCAL FUNCTION DEFINITIONS 85b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; Function Prototype declaration 86b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/ 87b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 88b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*---------------------------------------------------------------------------- 89b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; LOCAL VARIABLE DEFINITIONS 90b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; Variable declaration - defined here and used outside this module 91b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/ 92b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 93b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/* 94b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 95b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber FUNCTION NAME: subframePostProc 96b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 97b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber INPUT AND OUTPUT DEFINITIONS 98b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 99b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 100b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Inputs: 101b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber speech -- Pointer to Word16 -- speech segment 102b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber mode -- enum Mode -- coder mode 103b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber i_subfr -- Word16 -- Subframe nr 104b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber gain_pit -- Word16 -- Pitch gain Q14 105b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber gain_code -- Word16 -- Decoded innovation gain 106b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Aq -- Pointer to Word16 -- A(z) quantized for the 4 subframes 107b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber synth -- Word16 Array -- Local synthesis 108b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber xn -- Word16 Array -- Target vector for pitch search 109b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber code -- Word16 Array -- Fixed codebook exitation 110b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber y1 -- Word16 Array -- Filtered adaptive exitation 111b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber y2 -- Word16 Array -- Filtered fixed codebook excitation 112b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber mem_syn -- Pointer to Word16 -- memory of synthesis filter 113b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 114b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Outputs: 115b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber mem_syn -- Pointer to Word16 -- memory of synthesis filter 116b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber mem_err -- Pointer to Word16 -- pointer to error signal 117b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber mem_w0 -- Pointer to Word16 -- memory of weighting filter 118b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exc -- Pointer to Word16 -- long term prediction residual 119b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sharp -- Pointer to Word16 -- pitch sharpening value 120b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber pOverflow -- Pointer to Flag -- overflow indicator 121b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 122b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Returns: 123b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber None 124b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 125b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Global Variables Used: 126b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber None 127b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 128b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Local Variables Needed: 129b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber None 130b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 131b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 132b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber FUNCTION DESCRIPTION 133b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 134b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 135b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 136b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber REQUIREMENTS 137b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 138b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber None 139b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 140b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 141b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber REFERENCES 142b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 143b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber spstproc.c, UMTS GSM AMR speech codec, R99 - Version 3.2.0, March 2, 2001 144b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 145b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 146b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber PSEUDO-CODE 147b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 148b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 149b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 150b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber RESOURCES USED [optional] 151b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 152b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber When the code is written for a specific target processor the 153b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber the resources used should be documented below. 154b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 155b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber HEAP MEMORY USED: x bytes 156b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 157b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber STACK MEMORY USED: x bytes 158b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 159b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber CLOCK CYCLES: (cycle count equation for this function) + (variable 160b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber used to represent cycle count for each subroutine 161b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber called) 162b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber where: (cycle count variable) = cycle count for [subroutine 163b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber name] 164b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 165b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 166b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber CAUTION [optional] 167b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber [State any special notes, constraints or cautions for users of this function] 168b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 169b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 170b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber*/ 171b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 172b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Hubervoid subframePostProc( 173b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *speech, /* i : speech segment */ 174b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber enum Mode mode, /* i : coder mode */ 175b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 i_subfr, /* i : Subframe nr */ 176b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 gain_pit, /* i : Pitch gain Q14 */ 177b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 gain_code, /* i : Decoded innovation gain */ 178b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *Aq, /* i : A(z) quantized for the 4 subframes */ 179b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 synth[], /* i : Local snthesis */ 180b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 xn[], /* i : Target vector for pitch search */ 181b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 code[], /* i : Fixed codebook exitation */ 182b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 y1[], /* i : Filtered adaptive exitation */ 183b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 y2[], /* i : Filtered fixed codebook excitation */ 184b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *mem_syn, /* i/o : memory of synthesis filter */ 185b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *mem_err, /* o : pointer to error signal */ 186b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *mem_w0, /* o : memory of weighting filter */ 187b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *exc, /* o : long term prediction residual */ 188b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *sharp, /* o : pitch sharpening value */ 189b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Flag *pOverflow /* o : overflow indicator */ 190b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber) 191b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber{ 192b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 i; 193b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 j; 194b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 temp; 195b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word32 L_temp; 196b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word32 L_temp2; 197b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 tempShift; 198b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 kShift; 199b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 pitch_fac; 200b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *p_exc; 201b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *p_code; 202b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 203b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber OSCL_UNUSED_ARG(pOverflow); 204b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 205b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber if (mode != MR122) 206b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 207b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tempShift = 1; 208b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber kShift = 16 - 2 - 1; 209b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber pitch_fac = gain_pit; 210b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 211b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber else 212b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 213b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tempShift = 2; 214b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber kShift = 16 - 4 - 1; 215b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber pitch_fac = gain_pit >> 1; 216b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 217b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 218b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /*------------------------------------------------------------* 219b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * - Update pitch sharpening "sharp" with quantized gain_pit * 220b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *------------------------------------------------------------*/ 221b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 222b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber if (gain_pit < SHARPMAX) 223b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 224b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *sharp = gain_pit; 225b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 226b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber else 227b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 228b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *sharp = SHARPMAX; 229b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 230b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 231b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /*------------------------------------------------------* 232b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * - Find the total excitation * 233b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * - find synthesis speech corresponding to exc[] * 234b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * - update filters memories for finding the target * 235b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * vector in the next subframe * 236b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * (update error[-m..-1] and mem_w_err[]) * 237b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *------------------------------------------------------*/ 238b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 239b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber p_exc = &exc[ i_subfr]; 240b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber p_code = &code[0]; 241b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 242b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (i = L_SUBFR >> 1; i != 0 ; i--) 243b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 244b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* exc[i] = gain_pit*exc[i] + gain_code*code[i]; */ 245b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 246b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* 247b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * 12k2 others 248b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * --------------------------------- 249b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * exc Q0 Q0 250b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * gain_pit Q14 Q14 251b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * pitch_fac Q13 Q14 252b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * product: Q14 Q15 253b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * 254b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * code Q12 Q13 255b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * gain_code Q1 Q1 256b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * product Q14 Q15 257b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * sum Q14 Q15 258b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * 259b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * tempShift 2 1 260b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * sum<<tempShift Q16 Q16 261b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * result -> exc Q0 Q0 262b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber */ 263b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_temp = ((Word32) * (p_exc++) * pitch_fac) << 1; 264b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_temp2 = ((Word32) * (p_exc--) * pitch_fac) << 1; 265b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_temp += ((Word32) * (p_code++) * gain_code) << 1; 266b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_temp2 += ((Word32) * (p_code++) * gain_code) << 1; 267b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_temp <<= tempShift; 268b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_temp2 <<= tempShift; 269b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *(p_exc++) = (Word16)((L_temp + 0x08000L) >> 16); 270b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *(p_exc++) = (Word16)((L_temp2 + 0x08000L) >> 16); 271b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 272b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 273b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 274b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Syn_filt( 275b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Aq, 276b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber &exc[i_subfr], 277b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber &synth[i_subfr], 278b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_SUBFR, 279b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber mem_syn, 280b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1); 281b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 282b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (i = L_SUBFR - M, j = 0; i < L_SUBFR; i++, j++) 283b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 284b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber mem_err[j] = speech[i_subfr + i] - synth[i_subfr + i]; 285b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 286b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* 287b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * 12k2 others 288b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * --------------------------------- 289b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * y1 Q0 Q0 290b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * gain_pit Q14 Q14 291b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * product Q15 Q15 292b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * shifted prod. Q16 Q16 293b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * temp Q0 Q0 294b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * 295b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * y2 Q10 Q12 296b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * gain_code Q1 Q1 297b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * product Q12 Q14 298b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * kshift 4 2 299b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * shifted prod. Q16 Q16 300b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * k Q0 Q0 301b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * mem_w0,xn,sum Q0 Q0 302b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber */ 303b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 304b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_temp = ((Word32)y1[i] * gain_pit); 305b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber temp = (Word16)(L_temp >> 14); 306b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 307b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_temp = ((Word32)y2[i] * gain_code); 308b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber temp += (Word16)(L_temp >> kShift); 309b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 310b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber mem_w0[j] = xn[i] - temp; 311b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 312b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 313b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber return; 314b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber} 315