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 Pathname: ./audio/gsm-amr/c/src/lsp_az.c 324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Funtions: Get_lsp_pol 334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Lsp_Az 344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REVISION HISTORY 374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Updated template used to PV coding template. First attempt at 394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber optimizing C code. 404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Deleted all Local stores needed/modified. Optimized Lsp_Az 424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber function by getting rid of call to L_shr_r function. 434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Updated file per comments gathered from Phase 2/3 review. 454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Added setting of Overflow flag in the inlined code. 474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: 1. Optimized Lsp_Az function code. 494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2. Changed Input/Output definitions by adding Word type. 504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Made changes based on review meeting. 524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1. Removed pseudocode. 534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Synchronized file with UMTS version 3.2.0. Updated coding 554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber template. Removed unnecessary include files. 564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Replaced basic_op.h and oper_32b.h with the header files of the 584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber math functions used in the file. 594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Modified to pass overflow flag through to L_add and L_sub. The 614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber flag is passed back to the calling function by pointer reference. 624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Removed the id line since it was removed in the header file by 644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Ken. 654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Added the write-only variable, pOverflow, to the inputs section. 674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: For lsp_az() and Get_lsp_pol() 694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1. Eliminated unused include files. 704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2. Replaced array addressing by pointers 714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3. Eliminated math operations that unnecessary checked for 724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber saturation, in some cases this by shifting before adding and 734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber in other cases by evaluating the operands 744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4. Unrolled loops to speed up processing 754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5. Replaced mpy_32_16 by multpilcations in place 764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 6. Eliminated if-else statements for sign extension when 774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber right-shifting 784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Added casting to eliminate warnings, and eliminated include 804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber files that now are chosen by OSCL definitions 814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Replaced "int" and/or "char" with defined types. 834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Added proper casting (Word32) to some left shifting operations 844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Who: Date: 864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: 874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber MODULE DESCRIPTION 904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber This file contains functions that convert line spectral pairs (LSP) to 924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber linear predictive (LP) coefficients (filter order = 10). The functions 934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber included in this file include Get_lsp_pol, which finds the coefficients of 944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber F1(z) and F2(z), and Lsp_Az, which converts LSP to LPC by multiplying 954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber F1(z) by 1+z^(-1) and F2(z) by 1-z^(-1), then calculating A(z) = (F1(z) + 964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber F2(z))/2. 974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber*/ 1004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 1024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; INCLUDES 1034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 1044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "lsp_az.h" 1054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 1074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; MACROS 1084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Define module specific macros here 1094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 1104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 1124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; DEFINES 1134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Include all pre-processor statements here. Include conditional 1144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; compile variables also. 1154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 1164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 1184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; LOCAL FUNCTION DEFINITIONS 1194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Function Prototype declaration 1204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 1214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 1234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; LOCAL VARIABLE DEFINITIONS 1244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Variable declaration - defined here and used outside this module 1254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 1264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* 1294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber FUNCTION NAME: Get_lsp_pol 1314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber INPUT AND OUTPUT DEFINITIONS 1334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Inputs: 1354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber lsp = pointer to the buffer containing the line spectral pairs (LSP) 1364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber of type Word16 1374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber f = pointer to the polynomial of type Word32 to be generated 1384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber pOverflow = pointer set in case where one of the operations overflows. 1404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber [data type Pointer to Flag] 1414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Outputs: 1434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber buffer pointed to by f contains the polynomial generated 1444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber pOverflow = pointer set in case where one of the operations overflows. 1464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber [data type Pointer to Flag] 1474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Returns: 1494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 1504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Global Variables Used: 1524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 1534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Local Variables Needed: 1554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 1564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber FUNCTION DESCRIPTION 1594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber This function finds the polynomial F1(z) or F2(z) from the LSPs. If the LSP 1614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber vector is passed at address 0, F1(z) is computed and if it is passed at 1624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber address 1, F2(z) is computed. 1634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber This is performed by expanding the product polynomials: 1654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber F1(z) = product ( 1 - 2 lsp[i] z^-1 + z^-2 ) 1674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber i=0,2,4,6,8 1684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber F2(z) = product ( 1 - 2 lsp[i] z^-1 + z^-2 ) 1694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber i=1,3,5,7,9 1704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber where lsp[] is the LSP vector in the cosine domain. 1724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber The expansion is performed using the following recursion: 1744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber f[0] = 1 1764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber b = -2.0 * lsp[0] 1774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber f[1] = b 1784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for i=2 to 5 do 1794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber b = -2.0 * lsp[2*i-2]; 1804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for j=i-1 down to 2 do 1814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber f[j] = f[j] + b*f[j-1] + f[j-2]; 1824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber f[1] = f[1] + b; 1834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REQUIREMENTS 1864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 1884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REFERENCES 1914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber lsp_az.c, UMTS GSM AMR speech codec, R99 - Version 3.2.0, March 2, 2001 1934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber PSEUDO-CODE 1964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huberstatic void Get_lsp_pol (Word16 *lsp, Word32 *f) 1984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 1994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 i, j, hi, lo; 2004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 t0; 2014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // f[0] = 1.0; 2034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *f = L_mult (4096, 2048); 2044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber f++; 2054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *f = L_msu ((Word32) 0, *lsp, 512); // f[1] = -2.0 * lsp[0]; 2064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber f++; 2074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber lsp += 2; // Advance lsp pointer 2084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 2; i <= 5; i++) 2104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 2114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *f = f[-2]; 2124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (j = 1; j < i; j++, f--) 2144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 2154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_Extract (f[-1], &hi, &lo); 2164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber t0 = Mpy_32_16 (hi, lo, *lsp); // t0 = f[-1] * lsp 2174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber t0 = L_shl (t0, 1); 2184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *f = L_add (*f, f[-2]); // *f += f[-2] 2194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *f = L_sub (*f, t0); // *f -= t0 2204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 2214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *f = L_msu (*f, *lsp, 512); // *f -= lsp<<9 2224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber f += i; // Advance f pointer 2234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber lsp += 2; // Advance lsp pointer 2244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 2254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber return; 2274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 2284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 2304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber RESOURCES USED [optional] 2314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber When the code is written for a specific target processor the 2334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber the resources used should be documented below. 2344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber HEAP MEMORY USED: x bytes 2364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber STACK MEMORY USED: x bytes 2384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CLOCK CYCLES: (cycle count equation for this function) + (variable 2404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber used to represent cycle count for each subroutine 2414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber called) 2424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber where: (cycle count variable) = cycle count for [subroutine 2434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber name] 2444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 2464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CAUTION [optional] 2474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber [State any special notes, constraints or cautions for users of this function] 2484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 2504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber*/ 2514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huberstatic void Get_lsp_pol( 2534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *lsp, 2544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 *f, 2554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Flag *pOverflow) 2564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 2574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber register Word16 i; 2584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber register Word16 j; 2594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 hi; 2614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 lo; 2624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 t0; 2634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber OSCL_UNUSED_ARG(pOverflow); 2644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* f[0] = 1.0; */ 2664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *f++ = (Word32) 0x01000000; 2674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *f++ = (Word32) - *(lsp++) << 10; /* f[1] = -2.0 * lsp[0]; */ 2684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber lsp++; /* Advance lsp pointer */ 2694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 2; i <= 5; i++) 2714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 2724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *f = *(f - 2); 2734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (j = 1; j < i; j++) 2754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 2764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber hi = (Word16)(*(f - 1) >> 16); 2774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber lo = (Word16)((*(f - 1) >> 1) - ((Word32) hi << 15)); 2794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber t0 = ((Word32)hi * *lsp); 2814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber t0 += ((Word32)lo * *lsp) >> 15; 2824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *(f) += *(f - 2); /* *f += f[-2] */ 2844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *(f--) -= t0 << 2; /* *f -= t0 */ 2854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 2874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *f -= (Word32)(*lsp++) << 10; 2894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber f += i; 2914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber lsp++; 2924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 2934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber return; 2954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 2964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/****************************************************************************/ 2984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* 3014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 3024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber FUNCTION NAME: Get_lsp_pol_wrapper 3034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 3044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber INPUT AND OUTPUT DEFINITIONS 3054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Inputs: 3074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber lsp = pointer to the buffer containing the line spectral pairs (LSP) 3084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber of type Word16 3094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber f = pointer to the polynomial of type Word32 to be generated 3104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber pOverflow = pointer set in case where one of the operations overflows. 3124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber [data type Pointer to Flag] 3134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Outputs: 3154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber buffer pointed to by f contains the polynomial generated 3164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber pOverflow = pointer set in case where one of the operations overflows. 3184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber [data type Pointer to Flag] 3194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Returns: 3214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 3224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Global Variables Used: 3244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 3254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Local Variables Needed: 3274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 3284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 3304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber FUNCTION DESCRIPTION 3314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber This function provides external access to the static function Get_lsp_pol. 3334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 3354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REQUIREMENTS 3364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 3384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 3404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REFERENCES 3414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 3434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 3454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber PSEUDO-CODE 3464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CALL Get_lsp_pol(lsp = lsp_ptr 3484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber f = f_ptr ) 3494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber MODIFYING(nothing) 3504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber RETURNING(nothing) 3514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 3534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber RESOURCES USED [optional] 3544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber When the code is written for a specific target processor the 3564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber the resources used should be documented below. 3574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber HEAP MEMORY USED: x bytes 3594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber STACK MEMORY USED: x bytes 3614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CLOCK CYCLES: (cycle count equation for this function) + (variable 3634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber used to represent cycle count for each subroutine 3644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber called) 3654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber where: (cycle count variable) = cycle count for [subroutine 3664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber name] 3674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 3694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CAUTION [optional] 3704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber [State any special notes, constraints or cautions for users of this function] 3714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 3734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber*/ 3744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Hubervoid Get_lsp_pol_wrapper( 3764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *lsp, 3774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 *f, 3784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Flag *pOverflow) 3794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 3804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /*---------------------------------------------------------------------------- 3814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CALL Get_lsp_pol(lsp = lsp_ptr 3824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber f = f_ptr ) 3834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ----------------------------------------------------------------------------*/ 3844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Get_lsp_pol(lsp, f, pOverflow); 3854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /*---------------------------------------------------------------------------- 3874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber MODIFYING(nothing) 3884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber RETURNING(nothing) 3894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ----------------------------------------------------------------------------*/ 3904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber return; 3914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 3924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/****************************************************************************/ 3944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* 3974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 3984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber FUNCTION NAME: Lsp_Az 3994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 4004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber INPUT AND OUTPUT DEFINITIONS 4014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Inputs: 4034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber lsp = pointer to the buffer containing the line spectral pairs (LSP) 4044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber of type Word16 4054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber a = pointer to the buffer containing Linear Predictive (LP) 4074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber coefficients of type Word16 to be generated 4084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber pOverflow = pointer set in case where one of the operations overflows. 4104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber [data type Pointer to Flag] 4114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Local Stores/Buffers/Pointers Needed: 4134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 4144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Global Stores/Buffers/Pointers Needed: 4164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 4174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Outputs: 4194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber pOverflow = pointer set in case where one of the operations overflows. 4204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber [data type Pointer to Flag] 4214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Pointers and Buffers Modified: 4234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber a buffer contains the generated Linear Predictive (LP) coefficients 4244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Local Stores Modified: 4264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 4274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Global Stores Modified: 4294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 4304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 4324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber FUNCTION DESCRIPTION 4334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber This function converts from the line spectral pairs (LSP) to LP coefficients 4354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for a 10th order filter. 4364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber This is done by: 4384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber (1) Find the coefficients of F1(z) and F2(z) (see Get_lsp_pol) 4394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber (2) Multiply F1(z) by 1+z^{-1} and F2(z) by 1-z^{-1} 4404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber (3) A(z) = ( F1(z) + F2(z) ) / 2 4414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 4434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REQUIREMENTS 4444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 4464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 4484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REFERENCES 4494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber lsp_az.c, UMTS GSM AMR speech codec, R99 - Version 3.2.0, March 2, 2001 4514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 4534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber PSEUDO-CODE 4544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Hubervoid Lsp_Az ( 4564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 lsp[], // (i) : line spectral frequencies 4574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 a[] // (o) : predictor coefficients (order = 10) 4584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber) 4594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 4604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 i, j; 4614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 f1[6], f2[6]; 4624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 t0; 4634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Get_lsp_pol (&lsp[0], f1); 4654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Get_lsp_pol (&lsp[1], f2); 4664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 5; i > 0; i--) 4684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 4694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber f1[i] = L_add (f1[i], f1[i - 1]); // f1[i] += f1[i-1]; 4704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber f2[i] = L_sub (f2[i], f2[i - 1]); // f2[i] -= f2[i-1]; 4714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 4724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber a[0] = 4096; 4744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 1, j = 10; i <= 5; i++, j--) 4754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 4764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber t0 = L_add (f1[i], f2[i]); // f1[i] + f2[i] 4774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber a[i] = extract_l (L_shr_r (t0, 13)); 4784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber t0 = L_sub (f1[i], f2[i]); // f1[i] - f2[i] 4794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber a[j] = extract_l (L_shr_r (t0, 13)); 4804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 4814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber return; 4834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 4844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 4864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber RESOURCES USED [optional] 4874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber When the code is written for a specific target processor the 4894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber the resources used should be documented below. 4904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber HEAP MEMORY USED: x bytes 4924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber STACK MEMORY USED: x bytes 4944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CLOCK CYCLES: (cycle count equation for this function) + (variable 4964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber used to represent cycle count for each subroutine 4974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber called) 4984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber where: (cycle count variable) = cycle count for [subroutine 4994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber name] 5004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 5024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CAUTION [optional] 5034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber [State any special notes, constraints or cautions for users of this function] 5044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 5064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber*/ 5074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Hubervoid Lsp_Az( 5094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 lsp[], /* (i) : line spectral frequencies */ 5104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 a[], /* (o) : predictor coefficients (order = 10) */ 5114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Flag *pOverflow /* (o) : overflow flag */ 5124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber) 5134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 5144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber register Word16 i; 5154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber register Word16 j; 5164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 f1[6]; 5184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 f2[6]; 5194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 t0; 5204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 t1; 5214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *p_a = &a[0]; 5224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 *p_f1; 5234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 *p_f2; 5244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Get_lsp_pol(&lsp[0], f1, pOverflow); 5264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Get_lsp_pol(&lsp[1], f2, pOverflow); 5284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber p_f1 = &f1[5]; 5304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber p_f2 = &f2[5]; 5314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 5; i > 0; i--) 5334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 5344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *(p_f1--) += f1[i-1]; 5354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *(p_f2--) -= f2[i-1]; 5364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 5374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *(p_a++) = 4096; 5394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber p_f1 = &f1[1]; 5404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber p_f2 = &f2[1]; 5414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 1, j = 10; i <= 5; i++, j--) 5434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 5444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber t0 = *(p_f1) + *(p_f2); /* f1[i] + f2[i] */ 5454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber t1 = *(p_f1++) - *(p_f2++); /* f1[i] - f2[i] */ 5464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber t0 = t0 + ((Word32) 1 << 12); 5484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber t1 = t1 + ((Word32) 1 << 12); 5494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *(p_a++) = (Word16)(t0 >> 13); 5514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber a[j] = (Word16)(t1 >> 13); 5524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 5534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber return; 5554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 556