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/convolve.c 35b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 36b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Date: 06/19/2000 37b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 38b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 39b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber REVISION HISTORY 40b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 41b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: Optimize for speed. Update to code template. 42b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 43b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: Added author name and date, fixed tabs, and added missing 44b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sections. Updated Input/Output section. 45b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 46b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: Optimized code by calculating two convolution sums per iteration 47b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber of the outer loop, thereby, decreasing outer loop count by 2. 48b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Updated input/output definitions to be the same as the assembly 49b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber file (convolve.asm). Left Pseudo-code section blank. 50b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 51b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: Deleted semi-colon in the Pointers modified section. 52b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 53b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: Synchronized file with UMTS version 3.2.0. Updated coding 54b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber template. Removed unnecessary include files. 55b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 56b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: Made the following changes per comments from Phase 2/3 review: 57b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1. Fixed typecasting issue with TI C compiler. 58b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 2. Modified FOR loop to count down, wherever applicable. 59b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 60b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: Made the following changes 61b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1. Unrolled the correlation loop. 62b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 2. Performed 2 correlation per pass per sample to avoid recalling 63b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber the same data twice. 64b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 3. Eliminated math operations that check for saturation. 65b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 66b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: 67b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1. Modified loop counter, extra unrolling did speed up code 68b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 69b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: Replaced "int" and/or "char" with OSCL defined types. 70b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 71b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: Using inlines from fxp_arithmetic.h . 72b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 73b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: Replacing fxp_arithmetic.h with basic_op.h. 74b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 75b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: 76b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 77b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 78b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber*/ 79b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 80b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*---------------------------------------------------------------------------- 81b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; INCLUDES 82b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/ 83b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "typedef.h" 84b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "convolve.h" 85b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "basic_op.h" 86b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 87b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*---------------------------------------------------------------------------- 88b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; MACROS 89b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; Define module specific macros here 90b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/ 91b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 92b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 93b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*---------------------------------------------------------------------------- 94b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; DEFINES 95b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; Include all pre-processor statements here. Include conditional 96b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; compile variables also. 97b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/ 98b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 99b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 100b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*---------------------------------------------------------------------------- 101b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; LOCAL FUNCTION DEFINITIONS 102b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; Function Prototype declaration 103b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/ 104b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 105b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*---------------------------------------------------------------------------- 106b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; LOCAL STORE/BUFFER/POINTER DEFINITIONS 107b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; Variable declaration - defined here and used outside this module 108b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/ 109b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 110b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 111b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/* 112b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 113b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber FUNCTION NAME: Convolve 114b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 115b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber INPUT AND OUTPUT DEFINITIONS 116b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 117b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Inputs: 118b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber x = pointer to input vector of L elements of type Word16 119b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber h = pointer to the filter's impulse response vector of L elements 120b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber of type Word16 121b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber y = pointer to the output vector of L elements of type Word16 used for 122b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber storing the convolution of x and h; 123b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L = Length of the convolution; type definition is Word16 124b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 125b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Outputs: 126b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber y buffer contains the new convolution output 127b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 128b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Returns: 129b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber None 130b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 131b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Global Variables Used: 132b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber None 133b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 134b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Local Variables Needed: 135b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber None 136b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 137b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 138b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber FUNCTION DESCRIPTION 139b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 140b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Perform the convolution between two vectors x[] and h[] and write the result 141b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber in the vector y[]. All vectors are of length L and only the first L samples 142b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber of the convolution are computed. 143b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 144b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber The convolution is given by: 145b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 146b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber y[n] = sum_{i=0}^{n} x[i] h[n-i], n=0,...,L-1 147b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 148b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 149b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber REQUIREMENTS 150b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 151b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber None 152b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 153b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 154b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber REFERENCES 155b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 156b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber convolve.c, UMTS GSM AMR speech codec, R99 - Version 3.2.0, March 2, 2001 157b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 158b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 159b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber PSEUDO-CODE 160b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 161b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Hubervoid Convolve ( 162b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 x[], // (i) : input vector 163b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 h[], // (i) : impulse response 164b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 y[], // (o) : output vector 165b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 L // (i) : vector size 166b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber) 167b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber{ 168b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 i, n; 169b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word32 s; 170b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 171b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (n = 0; n < L; n++) 172b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 173b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s = 0; move32 (); 174b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (i = 0; i <= n; i++) 175b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 176b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s = L_mac (s, x[i], h[n - i]); 177b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 178b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s = L_shl (s, 3); 179b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber y[n] = extract_h (s); move16 (); 180b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 181b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 182b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber return; 183b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber} 184b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 185b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 186b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber RESOURCES USED [optional] 187b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 188b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber When the code is written for a specific target processor the 189b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber the resources used should be documented below. 190b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 191b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber HEAP MEMORY USED: x bytes 192b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 193b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber STACK MEMORY USED: x bytes 194b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 195b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber CLOCK CYCLES: (cycle count equation for this function) + (variable 196b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber used to represent cycle count for each subroutine 197b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber called) 198b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber where: (cycle count variable) = cycle count for [subroutine 199b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber name] 200b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 201b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 202b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber CAUTION [optional] 203b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber [State any special notes, constraints or cautions for users of this function] 204b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 205b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 206b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber*/ 207b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 208b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Hubervoid Convolve( 209b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 x[], /* (i) : input vector */ 210b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 h[], /* (i) : impulse response */ 211b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 y[], /* (o) : output vector */ 212b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 L /* (i) : vector size */ 213b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber) 214b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber{ 2152d0ac425564ff9882ebaac5267d1a04d4af67d00Bernhard Rosenkränzer Word16 i, n; 216b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word32 s1, s2; 217b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 218b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 219b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (n = 1; n < L; n = n + 2) 220b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 221b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 222b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber h = h + n; 223b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 224b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s2 = ((Word32) * (x)) * *(h--); 225b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s1 = ((Word32) * (x++)) * *(h); 226b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 227b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (i = (n - 1) >> 1; i != 0; i--) 228b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 229b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s2 = amrnb_fxp_mac_16_by_16bb((Word32) * (x), (Word32) * (h--), s2); 230b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s1 = amrnb_fxp_mac_16_by_16bb((Word32) * (x++), (Word32) * (h), s1); 231b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s2 = amrnb_fxp_mac_16_by_16bb((Word32) * (x), (Word32) * (h--), s2); 232b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s1 = amrnb_fxp_mac_16_by_16bb((Word32) * (x++), (Word32) * (h), s1); 233b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 234b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 235b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s2 = amrnb_fxp_mac_16_by_16bb((Word32) * (x), (Word32) * (h), s2); 236b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 237b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *(y++) = (Word16)(s1 >> 12); 238b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *(y++) = (Word16)(s2 >> 12); 239b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 240b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber x = x - n; 241b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 242b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 243b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 244b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber return; 245b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber} 246