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/cor_h.c 35b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 36b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Date: 06/12/2000 37b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 38b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 39b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber REVISION HISTORY 40b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 41b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: Updated template used to PV coding template. First attempt at 42b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber optimizing C code. 43b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 44b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: Used MAX_16 and MIN_16 when checking the result of Inv_sqrt. 45b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Synced up to the new template. 46b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 47b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: Added setting of Overflow flag in inlined code. 48b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 49b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: Took out cor_h_x function and put it in its own file. Sync'ed 50b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber up with the single_func_template.c template. Delete version 51b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber ID variable. 52b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 53b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: Synchronized file with UTMS version 3.2.0. Updated coding 54b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber template. Removed unnecessary include files. 55b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 56b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: Fixed portion of the code that builds the rr[] matrix. There 57b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber was an error in the original inlining of code that caused 58b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber the code to be not bit-exact with UMTS version 3.2.0. 59b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 60b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: Added calls to L_add() and mult() in the code to handle overflow 61b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber scenario. Moved cor_h.h after cnst.h in the Include section. 62b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Doing this allows the unit test to build using the cnst.h in the 63b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /test/include directory. Fixed initialization of the accumulator 64b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber in the first calculation of the sum of squares. 65b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 66b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: Made the following changes per comments from Phase 2/3 review: 67b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1. Used #define value instead of hard-coded numbers in the code. 68b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 2. Fixed typecasting issue with TI C compiler. 69b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 3. Removed typecasting of 0x00008000L in the call to L_add. 70b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 71b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: Changed pOverflow from a global variable into a function 72b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber parameter. 73b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 74b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: 75b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1. Added pointer to avoid adding offsets in every pass 76b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 2. Eliminate variables defined as registers 77b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 3. Removed extra check for overflow by doing scaling right 78b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber after overflow is detected. 79b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 4. Eliminated calls to basic operations (like extract) not 80b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber needed because of the nature of the number (all bounded) 81b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 5. Eliminated duplicate loop accessing same data 82b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 6. Simplified matrix addressing by use of pointers 83b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 84b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: 85b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1. Eliminated unused include files. 86b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 2. Access twice the number of points when delaing with matrices 87b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber and in the process only 3 pointers (instead of 4) are needed 88b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 3. Replaced array addressing (array sign[]) by pointers 89b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 90b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: Changed round function name to pv_round to avoid conflict with 91b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber round function in C standard library. 92b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 93b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: Using inlines from fxp_arithmetic.h . 94b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 95b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: Replacing fxp_arithmetic.h with basic_op.h. 96b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 97b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Description: 98b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 99b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 100b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber*/ 101b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 102b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*---------------------------------------------------------------------------- 103b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; INCLUDES 104b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/ 105b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "cnst.h" 106b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "cor_h.h" 107b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "basicop_malloc.h" 108b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "inv_sqrt.h" 109b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "basic_op.h" 110b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 111b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*---------------------------------------------------------------------------- 112b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; MACROS 113b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; Define module specific macros here 114b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/ 115b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 116b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*---------------------------------------------------------------------------- 117b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; DEFINES 118b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; Include all pre-processor statements here. Include conditional 119b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; compile variables also. 120b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/ 121b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 122b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*---------------------------------------------------------------------------- 123b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; LOCAL FUNCTION DEFINITIONS 124b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; Function Prototype declaration 125b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/ 126b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 127b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*---------------------------------------------------------------------------- 128b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; LOCAL STORE/BUFFER/POINTER DEFINITIONS 129b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; Variable declaration - defined here and used outside this module 130b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/ 131b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 132b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/* 133b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 134b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber FUNCTION NAME: cor_h 135b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 136b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber INPUT AND OUTPUT DEFINITIONS 137b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 138b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Inputs: 139b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber h = vector containing the impulse response of the weighted synthesis 140b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber filter; vector contents are of type Word16; vector length is 141b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 2 * L_SUBFR 142b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sign = vector containing the sign information for the correlation 143b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber values; vector contents are of type Word16; vector length is 144b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_CODE 145b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber rr = autocorrelation matrix; matrix contents are of type Word16; 146b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber matrix dimension is L_CODE by L_CODE 147b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 148b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Outputs: 149b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber rr contents are the newly calculated autocorrelation values 150b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 151b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Returns: 152b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber None 153b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 154b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Global Variables Used: 155b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber None 156b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 157b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Local Variables Needed: 158b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber None 159b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 160b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 161b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber FUNCTION DESCRIPTION 162b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 163b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber This function computes correlations of the impulse response (h) needed for 164b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber the codebook search, and includes the sign information into the correlations. 165b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 166b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber The correlations are given by: 167b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber rr[i][j] = sum_{n=i}^{L-1} h[n-i] h[n-j]; i>=j; i,j=0,...,L-1 168b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 169b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber The sign information is included by: 170b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber rr[i][j] = rr[i][j]*sign[i]*sign[j] 171b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 172b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 173b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber REQUIREMENTS 174b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 175b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber None 176b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 177b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 178b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber REFERENCES 179b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 180b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber cor_h.c, UMTS GSM AMR speech codec, R99 - Version 3.2.0, March 2, 2001 181b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 182b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 183b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber PSEUDO-CODE 184b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 185b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Hubervoid cor_h ( 186b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 h[], // (i) : impulse response of weighted synthesis 187b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber filter 188b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sign[], // (i) : sign of d[n] 189b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 rr[][L_CODE] // (o) : matrix of autocorrelation 190b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber) 191b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber{ 192b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 i, j, k, dec, h2[L_CODE]; 193b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word32 s; 194b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 195b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // Scaling for maximum precision 196b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 197b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s = 2; 198b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (i = 0; i < L_CODE; i++) 199b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s = L_mac (s, h[i], h[i]); 200b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 201b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber j = sub (extract_h (s), 32767); 202b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber if (j == 0) 203b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 204b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (i = 0; i < L_CODE; i++) 205b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 206b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber h2[i] = shr (h[i], 1); 207b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 208b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 209b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber else 210b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 211b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s = L_shr (s, 1); 212b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber k = extract_h (L_shl (Inv_sqrt (s), 7)); 213b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber k = mult (k, 32440); // k = 0.99*k 214b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 215b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (i = 0; i < L_CODE; i++) 216b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 217b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber h2[i] = pv_round (L_shl (L_mult (h[i], k), 9)); 218b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 219b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 220b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 221b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // build matrix rr[] 222b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s = 0; 223b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber i = L_CODE - 1; 224b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (k = 0; k < L_CODE; k++, i--) 225b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 226b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s = L_mac (s, h2[k], h2[k]); 227b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber rr[i][i] = pv_round (s); 228b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 229b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 230b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (dec = 1; dec < L_CODE; dec++) 231b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 232b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s = 0; 233b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber j = L_CODE - 1; 234b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber i = sub (j, dec); 235b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (k = 0; k < (L_CODE - dec); k++, i--, j--) 236b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 237b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s = L_mac (s, h2[k], h2[k + dec]); 238b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber rr[j][i] = mult (pv_round (s), mult (sign[i], sign[j])); 239b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber rr[i][j] = rr[j][i]; 240b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 241b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 242b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber} 243b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 244b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber--------------------------------------------------------------------------- 245b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber RESOURCES USED [optional] 246b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 247b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber When the code is written for a specific target processor the 248b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber the resources used should be documented below. 249b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 250b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber HEAP MEMORY USED: x bytes 251b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 252b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber STACK MEMORY USED: x bytes 253b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 254b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber CLOCK CYCLES: (cycle count equation for this function) + (variable 255b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber used to represent cycle count for each subroutine 256b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber called) 257b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber where: (cycle count variable) = cycle count for [subroutine 258b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber name] 259b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 260b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 261b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber CAUTION [optional] 262b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber [State any special notes, constraints or cautions for users of this function] 263b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 264b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 265b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber*/ 266b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 267b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Hubervoid cor_h( 268b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 h[], /* (i) : impulse response of weighted synthesis 269b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber filter */ 270b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sign[], /* (i) : sign of d[n] */ 271b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 rr[][L_CODE], /* (o) : matrix of autocorrelation */ 272b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Flag *pOverflow 273b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber) 274b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber{ 2752d0ac425564ff9882ebaac5267d1a04d4af67d00Bernhard Rosenkränzer Word16 i; 2762d0ac425564ff9882ebaac5267d1a04d4af67d00Bernhard Rosenkränzer Word16 dec; 277b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 278b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 h2[L_CODE]; 279b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word32 s; 280b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word32 s2; 281b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 tmp1; 282b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 tmp2; 283b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 tmp11; 284b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 tmp22; 285b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 286b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *p_h; 287b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *p_h2; 288b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *rr1; 289b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *rr2; 290b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *rr3; 291b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *p_rr_ref1; 292b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *p_sign1; 293b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *p_sign2; 294b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 295b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* Scaling for maximum precision */ 296b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 297b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* Initialize accumulator to 1 since left shift happens */ 298b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* after the accumulation of the sum of squares (original */ 299b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* code initialized s to 2) */ 300b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s = 1; 301b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber p_h = h; 302b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 303b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (i = (L_CODE >> 1); i != 0 ; i--) 304b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 305b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp1 = *(p_h++); 306b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s = amrnb_fxp_mac_16_by_16bb((Word32) tmp1, (Word32) tmp1, s); 307b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp1 = *(p_h++); 308b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s = amrnb_fxp_mac_16_by_16bb((Word32) tmp1, (Word32) tmp1, s); 309b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 310b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 311b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 312b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s <<= 1; 313b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 314b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber if (s & MIN_32) 315b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 316b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber p_h2 = h2; 317b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber p_h = h; 318b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 319b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (i = (L_CODE >> 1); i != 0; i--) 320b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 321b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *(p_h2++) = *(p_h++) >> 1; 322b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *(p_h2++) = *(p_h++) >> 1; 323b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 324b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 325b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber else 326b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 327b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 328b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s >>= 1; 329b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 330b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s = Inv_sqrt(s, pOverflow); 331b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 332b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber if (s < (Word32) 0x00ffffffL) 333b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 334b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* k = 0.99*k */ 335b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber dec = (Word16)(((s >> 9) * 32440) >> 15); 336b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 337b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber else 338b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 339b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber dec = 32440; /* 0.99 */ 340b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 341b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 342b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber p_h = h; 343b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber p_h2 = h2; 344b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 345b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (i = (L_CODE >> 1); i != 0; i--) 346b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 347b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *(p_h2++) = (Word16)((amrnb_fxp_mac_16_by_16bb((Word32) * (p_h++), (Word32) dec, 0x020L)) >> 6); 348b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *(p_h2++) = (Word16)((amrnb_fxp_mac_16_by_16bb((Word32) * (p_h++), (Word32) dec, 0x020L)) >> 6); 349b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 350b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 351b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* build matrix rr[] */ 352b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 353b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s = 0; 354b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 355b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber p_h2 = h2; 356b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 357b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber rr1 = &rr[L_CODE-1][L_CODE-1]; 358b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 359b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (i = L_CODE >> 1; i != 0 ; i--) 360b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 361b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp1 = *(p_h2++); 362b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s = amrnb_fxp_mac_16_by_16bb((Word32) tmp1, (Word32) tmp1, s); 363b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *rr1 = (Word16)((s + 0x00004000L) >> 15); 364b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber rr1 -= (L_CODE + 1); 365b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp1 = *(p_h2++); 366b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s = amrnb_fxp_mac_16_by_16bb((Word32) tmp1, (Word32) tmp1, s); 367b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *rr1 = (Word16)((s + 0x00004000L) >> 15); 368b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber rr1 -= (L_CODE + 1); 369b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 370b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 371b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 372b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber p_rr_ref1 = rr[L_CODE-1]; 373b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 374b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (dec = 1; dec < L_CODE; dec += 2) 375b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 376b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber rr1 = &p_rr_ref1[L_CODE-1-dec]; 377b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 378b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber rr2 = &rr[L_CODE-1-dec][L_CODE-1]; 379b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber rr3 = &rr[L_CODE-1-(dec+1)][L_CODE-1]; 380b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 381b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s = 0; 382b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s2 = 0; 383b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 384b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber p_sign1 = &sign[L_CODE - 1]; 385b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber p_sign2 = &sign[L_CODE - 1 - dec]; 386b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 387b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber p_h2 = h2; 388b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber p_h = &h2[dec]; 389b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 390b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (i = (L_CODE - dec - 1); i != 0 ; i--) 391b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 392b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s = amrnb_fxp_mac_16_by_16bb((Word32) * (p_h2), (Word32) * (p_h++), s); 393b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s2 = amrnb_fxp_mac_16_by_16bb((Word32) * (p_h2++), (Word32) * (p_h), s2); 394b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 395b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp1 = (Word16)((s + 0x00004000L) >> 15); 396b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp11 = (Word16)((s2 + 0x00004000L) >> 15); 397b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 398b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp2 = ((Word32) * (p_sign1) * *(p_sign2--)) >> 15; 399b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp22 = ((Word32) * (p_sign1--) * *(p_sign2)) >> 15; 400b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 401b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *rr2 = ((Word32) tmp1 * tmp2) >> 15; 402b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *(rr1--) = *rr2; 403b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *rr1 = ((Word32) tmp11 * tmp22) >> 15; 404b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *rr3 = *rr1; 405b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 406b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber rr1 -= (L_CODE); 407b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber rr2 -= (L_CODE + 1); 408b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber rr3 -= (L_CODE + 1); 409b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 410b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 411b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 412b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber s = amrnb_fxp_mac_16_by_16bb((Word32) * (p_h2), (Word32) * (p_h), s); 413b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 414b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp1 = (Word16)((s + 0x00004000L) >> 15); 415b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 416b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp2 = ((Word32) * (p_sign1) * *(p_sign2)) >> 15; 417b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *rr1 = ((Word32) tmp1 * tmp2) >> 15; 418b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 419b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *rr2 = *rr1; 420b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 421b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber rr1 -= (L_CODE + 1); 422b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber rr2 -= (L_CODE + 1); 423b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 424b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 425b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 426b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber return; 427b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 428b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber} 429b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 430