14f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* ------------------------------------------------------------------ 24f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * Copyright (C) 1998-2009 PacketVideo 34f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * 44f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * Licensed under the Apache License, Version 2.0 (the "License"); 54f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * you may not use this file except in compliance with the License. 64f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * You may obtain a copy of the License at 74f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * 84f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * http://www.apache.org/licenses/LICENSE-2.0 94f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * 104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * Unless required by applicable law or agreed to in writing, software 114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * distributed under the License is distributed on an "AS IS" BASIS, 124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either 134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * express or implied. 144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * See the License for the specific language governing permissions 154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * and limitations under the License. 164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * ------------------------------------------------------------------- 174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber */ 184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/**************************************************************************************** 194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas HuberPortions of this file are derived from the following 3GPP standard: 204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3GPP TS 26.073 224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ANSI-C code for the Adaptive Multi-Rate (AMR) speech codec 234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Available from http://www.3gpp.org 244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber(C) 2004, 3GPP Organizational Partners (ARIB, ATIS, CCSA, ETSI, TTA, TTC) 264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas HuberPermission to distribute, modify and use this file under the standard license 274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huberterms listed above has been obtained from the copyright holder. 284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber****************************************************************************************/ 294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* 304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Filename: /audio/gsm-amr/c/src/a_refl.c 354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Functions: a_refl 364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Date: 02/05/2002 384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REVISION HISTORY 414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Removing unneeded include files and the goto statement. 434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Changed function name to pv_round to avoid conflict with 464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber round function in C standard library. 474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Replaced "int" and/or "char" with OSCL defined types. 494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Using inline functions from basic_op.h . 514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Removing unneeded include files. 524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: 544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber*/ 574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; INCLUDES 614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "a_refl.h" 634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "typedef.h" 644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "cnst.h" 654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "basic_op.h" 664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; MACROS [optional] 694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; [Define module specific macros here] 704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; DEFINES [optional] 744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; [Include all pre-processor statements here. Include conditional 754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; compile variables also.] 764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; LOCAL FUNCTION DEFINITIONS 804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; [List function prototypes here] 814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; LOCAL VARIABLE DEFINITIONS 854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; [Variable declaration - defined here and used outside this module] 864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* 894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber FUNCTION NAME: AMREncode 914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber INPUT AND OUTPUT DEFINITIONS 934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Inputs: 954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber a[] = pointer to directform coefficients of type Word16 964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber refl[] = pointer to reflection coefficients of type Word16 974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Outputs: 994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber pOverflow = 1 if overflow exists in the math operations else zero. 1004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Returns: 1024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 1034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Global Variables Used: 1054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 1064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Local Variables Needed: 1084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 1094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber FUNCTION DESCRIPTION 1124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber File : a_refl.c 1144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Purpose : Convert from direct form coefficients to 1154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber reflection coefficients 1164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REQUIREMENTS 1194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 1214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REFERENCES 1244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber [1] a_refl.c , 3GPP TS 26.101 version 4.1.0 Release 4, June 2001 1264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber PSEUDO-CODE 1294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Hubervoid A_Refl( 1324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 a[], // i : Directform coefficients 1334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 refl[] // o : Reflection coefficients 1344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber) 1354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 1364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // local variables 1374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 i,j; 1384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 aState[M]; 1394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 bState[M]; 1404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 normShift; 1414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 normProd; 1424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 L_acc; 1434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 scale; 1444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 L_temp; 1454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 temp; 1464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 mult; 1474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // initialize states 1494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 0; i < M; i++) 1504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 1514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber aState[i] = a[i]; 1524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 1534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // backward Levinson recursion 1554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = M-1; i >= 0; i--) 1564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 1574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (sub(abs_s(aState[i]), 4096) >= 0) 1584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 1594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber goto ExitRefl; 1604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 1614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber refl[i] = shl(aState[i], 3); 1634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_temp = L_mult(refl[i], refl[i]); 1654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_acc = L_sub(MAX_32, L_temp); 1664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber normShift = norm_l(L_acc); 1684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber scale = sub(15, normShift); 1694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_acc = L_shl(L_acc, normShift); 1714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber normProd = pv_round(L_acc); 1724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber mult = div_s(16384, normProd); 1744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (j = 0; j < i; j++) 1764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 1774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_acc = L_deposit_h(aState[j]); 1784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_acc = L_msu(L_acc, refl[i], aState[i-j-1]); 1794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber temp = pv_round(L_acc); 1814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_temp = L_mult(mult, temp); 1824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_temp = L_shr_r(L_temp, scale); 1834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (L_sub(L_abs(L_temp), 32767) > 0) 1854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 1864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber goto ExitRefl; 1874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 1884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber bState[j] = extract_l(L_temp); 1904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 1914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (j = 0; j < i; j++) 1934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 1944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber aState[j] = bState[j]; 1954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 1964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 1974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber return; 1984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas HuberExitRefl: 2004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 0; i < M; i++) 2014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 2024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber refl[i] = 0; 2034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 2044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 2054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 2074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber RESOURCES USED [optional] 2084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber When the code is written for a specific target processor the 2104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber the resources used should be documented below. 2114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber HEAP MEMORY USED: x bytes 2134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber STACK MEMORY USED: x bytes 2154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CLOCK CYCLES: (cycle count equation for this function) + (variable 2174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber used to represent cycle count for each subroutine 2184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber called) 2194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber where: (cycle count variable) = cycle count for [subroutine 2204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber name] 2214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 2234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CAUTION [optional] 2244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber [State any special notes, constraints or cautions for users of this function] 2254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 2274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber*/ 2284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Hubervoid A_Refl( 2304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 a[], /* i : Directform coefficients */ 2314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 refl[], /* o : Reflection coefficients */ 2324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Flag *pOverflow 2334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber) 2344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 2354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* local variables */ 2364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 i; 2374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 j; 2384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 aState[M]; 2394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 bState[M]; 2404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 normShift; 2414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 normProd; 2424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 L_acc; 2434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 scale; 2444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 L_temp; 2454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 temp; 2464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 mult; 2474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* initialize states */ 2494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 0; i < M; i++) 2504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 2514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber aState[i] = a[i]; 2524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 2534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* backward Levinson recursion */ 2554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = M - 1; i >= 0; i--) 2564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 2574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (abs_s(aState[i]) >= 4096) 2584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 2594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 0; i < M; i++) 2604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 2614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber refl[i] = 0; 2624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 2634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber break; 2644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 2654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber refl[i] = shl(aState[i], 3, pOverflow); 2674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_temp = L_mult(refl[i], refl[i], pOverflow); 2694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_acc = L_sub(MAX_32, L_temp, pOverflow); 2704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber normShift = norm_l(L_acc); 2724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber scale = sub(15, normShift, pOverflow); 2734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_acc = L_shl(L_acc, normShift, pOverflow); 2754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber normProd = pv_round(L_acc, pOverflow); 2764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber mult = div_s(16384, normProd); 2784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (j = 0; j < i; j++) 2804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 2814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_acc = L_deposit_h(aState[j]); 2824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_acc = L_msu(L_acc, refl[i], aState[i-j-1], pOverflow); 2834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber temp = pv_round(L_acc, pOverflow); 2854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_temp = L_mult(mult, temp, pOverflow); 2864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_temp = L_shr_r(L_temp, scale, pOverflow); 2874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (L_abs(L_temp) > 32767) 2894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 2904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 0; i < M; i++) 2914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 2924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber refl[i] = 0; 2934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 2944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber break; 2954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 2964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber bState[j] = extract_l(L_temp); 2984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 2994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (j = 0; j < i; j++) 3014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber aState[j] = bState[j]; 3034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber return; 3064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 3074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 314