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 Pathname: ./audio/gsm-amr/c/src/syn_filt.c 314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REVISION HISTORY 344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Making changes based on comments from the review meeting. 364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Added typedef to Input/Output Definition section. 384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Synchronized file with UMTS version 3.2.0. Updated coding 404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber template. 414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Fixed typecasting issue with the TI C compiler. 434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Modified FOR loops to count down. 454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Modified FOR loop to count up again so that the correct values 474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber are stored in the tmp buffer. Updated copyright year. 484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: 504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber - Modified for loop and introduced pointers to avoid adding 514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber offsets 524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber - Eliminated check for saturation given that the max values of input 534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber data and coefficients will not saturate the multiply and 544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber accumulation 554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber - eliminated memcpy to update history buffer in every pass. This is 564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber done now just updating the pointers. 574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: 594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1. Eliminated unused include files. 604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2. Unrolled loops to process twice as many samples as before, 614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber this saves on memory accesses to the vector coeff. a[] and 624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber elements in the history buffer of this recursive filter 634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: 654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1. Added overflow check inside both loops. (this is needed just 664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber to satisfy bit exactness on the decoder, a faster 674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber implementation will add an extra shift, do the same, 684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber but will not be bit exact, and it may have better audio 694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber quality because will avoid clipping) 704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2. Added include file for constant definitions 714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Replaced OSCL mem type functions and eliminated include 734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber files that now are chosen by OSCL definitions 744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Replaced "int" and/or "char" with OSCL defined types. 764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Changed round function name to pv_round to avoid conflict with 784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber round function in C standard library. 794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Using fxp_arithmetic.h that includes inline assembly functions 814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for ARM and linux-arm. 824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Replacing fxp_arithmetic.h with basic_op.h. 844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Who: Date: 864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: 874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber*/ 904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; INCLUDES 934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include <string.h> 964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "syn_filt.h" 984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "cnst.h" 994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "basic_op.h" 1004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "basic_op.h" 1024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 1044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; MACROS 1054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Define module specific macros here 1064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 1074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 1104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; DEFINES 1114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Include all pre-processor statements here. Include conditional 1124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; compile variables also. 1134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 1144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 1164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; LOCAL FUNCTION DEFINITIONS 1174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Function Prototype declaration 1184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 1194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 1214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; LOCAL STORE/BUFFER/POINTER DEFINITIONS 1224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Variable declaration - defined here and used outside this module 1234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 1244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* 1274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber FUNCTION NAME: Syn_filt 1294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber INPUT AND OUTPUT DEFINITIONS 1314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Inputs: 1334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber a = buffer containing the prediction coefficients (Word16) max 2^12 1344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber x = input signal buffer (Word16) max 2^15 1354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber y = output signal buffer (Word16) 1364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber lg = size of filtering (Word16) 1374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber mem = memory buffer associated with this filtering (Word16) 1384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber update = flag to indicate memory update; 0=no update, 1=update memory 1394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber (Word16) 1404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Outputs: 1424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber mem buffer is changed to be the last M data points of the output signal 1434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if update was set to 1 1444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber y buffer contains the newly calculated filter output 1454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Returns: 1474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 1484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Global Variables Used: 1504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 1514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Local Variables Needed: 1534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 1544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber FUNCTION DESCRIPTION 1574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Perform synthesis filtering through 1/A(z) 1594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REQUIREMENTS 1624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 1644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REFERENCES 1674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber syn_filt.c, UMTS GSM AMR speech codec, R99 - Version 3.2.0, March 2, 2001 1694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber PSEUDO-CODE 1724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Hubervoid Syn_filt ( 1744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 a[], // (i) : a[M+1] prediction coefficients (M=10) 1754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 x[], // (i) : input signal 1764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 y[], // (o) : output signal 1774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 lg, // (i) : size of filtering 1784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 mem[], // (i/o) : memory associated with this filtering. 1794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 update // (i) : 0=no update, 1=update of memory. 1804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber) 1814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 1824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 i, j; 1834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 s; 1844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 tmp[80]; // This is usually done by memory allocation (lg+M) 1854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *yy; 1864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // Copy mem[] to yy[] 1884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber yy = tmp; 1904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 0; i < M; i++) 1924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 1934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *yy++ = mem[i]; 1944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 1954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // Do the filtering. 1974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 0; i < lg; i++) 1994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 2004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber s = L_mult (x[i], a[0]); 2014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (j = 1; j <= M; j++) 2024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 2034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber s = L_msu (s, a[j], yy[-j]); 2044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 2054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber s = L_shl (s, 3); 2064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *yy++ = pv_round (s); 2074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 2084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 0; i < lg; i++) 2104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 2114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber y[i] = tmp[i + M]; 2124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 2134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // Update of memory if update==1 2154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (update != 0) 2174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 2184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 0; i < M; i++) 2194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 2204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber mem[i] = y[lg - M + i]; 2214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 2224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 2234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber return; 2244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 2254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 2274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber RESOURCES USED [optional] 2284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber When the code is written for a specific target processor the 2304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber the resources used should be documented below. 2314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber HEAP MEMORY USED: x bytes 2334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber STACK MEMORY USED: x bytes 2354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CLOCK CYCLES: (cycle count equation for this function) + (variable 2374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber used to represent cycle count for each subroutine 2384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber called) 2394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber where: (cycle count variable) = cycle count for [subroutine 2404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber name] 2414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 2434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CAUTION [optional] 2444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber [State any special notes, constraints or cautions for users of this function] 2454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 2474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber*/ 248b875f69d8e867cb64bd101e66d85a880537c2b72Dan Austin#ifdef __clang__ 249b875f69d8e867cb64bd101e66d85a880537c2b72Dan Austin__attribute__((no_sanitize("integer"))) 250b875f69d8e867cb64bd101e66d85a880537c2b72Dan Austin#endif 2514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Hubervoid Syn_filt( 2524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 a[], /* (i) : a[M+1] prediction coefficients (M=10) */ 2534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 x[], /* (i) : input signal */ 2544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 y[], /* (o) : output signal */ 2554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 lg, /* (i) : size of filtering (40) */ 2564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 mem[], /* (i/o) : memory associated with this filtering. */ 2574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 update /* (i) : 0=no update, 1=update of memory. */ 2584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber) 2594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 2604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 i, j; 2614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 s1; 2624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 s2; 2634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 tmp[2*M]; /* This is usually done by memory allocation (lg+M) */ 2644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *yy; 2654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *p_a; 2674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *p_yy1; 2684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *p_y; 2694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *p_x; 2704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 temp; 2714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* Copy mem[] to yy[] */ 2724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber yy = tmp; 2744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber memcpy(yy, mem, M*sizeof(Word16)); 2764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber yy = yy + M; 2784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* Do the filtering. */ 2804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber p_y = y; 2824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber p_x = x; 2834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber p_yy1 = &yy[-1]; 2844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = M >> 1; i != 0; i--) 2864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 2874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber p_a = a; 2884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber s1 = amrnb_fxp_mac_16_by_16bb((Word32) * (p_x++), (Word32) * (p_a), 0x00000800L); 2904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber s2 = amrnb_fxp_mac_16_by_16bb((Word32) * (p_x++), (Word32) * (p_a++), 0x00000800L); 2914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber s1 = amrnb_fxp_msu_16_by_16bb((Word32) * (p_a++), (Word32) * (p_yy1), s1); 2924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (j = (M >> 1) - 2; j != 0; j--) 2944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 2954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber s2 = amrnb_fxp_msu_16_by_16bb((Word32) * (p_a), (Word32) * (p_yy1--), s2); 2964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber s1 = amrnb_fxp_msu_16_by_16bb((Word32) * (p_a++), (Word32) * (p_yy1), s1); 2974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber s2 = amrnb_fxp_msu_16_by_16bb((Word32) * (p_a), (Word32) * (p_yy1--), s2); 2984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber s1 = amrnb_fxp_msu_16_by_16bb((Word32) * (p_a++), (Word32) * (p_yy1), s1); 2994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber s2 = amrnb_fxp_msu_16_by_16bb((Word32) * (p_a), (Word32) * (p_yy1--), s2); 3004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber s1 = amrnb_fxp_msu_16_by_16bb((Word32) * (p_a++), (Word32) * (p_yy1), s1); 3014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* check for overflow on s1 */ 3044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if ((UWord32)(s1 - 0xf8000000L) < 0x0fffffffL) 3054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber temp = (Word16)(s1 >> 12); 3074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else if (s1 > 0x07ffffffL) 3094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber temp = MAX_16; 3114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else 3134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber temp = MIN_16; 3154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber s2 = amrnb_fxp_msu_16_by_16bb((Word32)a[1], (Word32)temp, s2); 3184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *(yy++) = temp; 3204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *(p_y++) = temp; 3214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber p_yy1 = yy; 3234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* check for overflow on s2 */ 3254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if ((UWord32)(s2 - 0xf8000000L) < 0x0fffffffL) 3264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber temp = (Word16)(s2 >> 12); 3284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else if (s2 > 0x07ffffffL) 3304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber temp = MAX_16; 3324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else 3344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber temp = MIN_16; 3364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *(yy++) = temp; 3394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *(p_y++) = temp; 3404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber p_yy1 = &y[M-1]; 3434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = (lg - M) >> 1; i != 0; i--) 3454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber p_a = a; 3474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber s1 = amrnb_fxp_mac_16_by_16bb((Word32) * (p_x++), (Word32) * (p_a), 0x00000800L); 3494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber s2 = amrnb_fxp_mac_16_by_16bb((Word32) * (p_x++), (Word32) * (p_a++), 0x00000800L); 3504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber s1 = amrnb_fxp_msu_16_by_16bb((Word32) * (p_a++), (Word32) * (p_yy1), s1); 3514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (j = (M >> 1) - 2; j != 0; j--) 3534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber s2 = amrnb_fxp_msu_16_by_16bb((Word32) * (p_a), (Word32) * (p_yy1--), s2); 3554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber s1 = amrnb_fxp_msu_16_by_16bb((Word32) * (p_a++), (Word32) * (p_yy1), s1); 3564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber s2 = amrnb_fxp_msu_16_by_16bb((Word32) * (p_a), (Word32) * (p_yy1--), s2); 3574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber s1 = amrnb_fxp_msu_16_by_16bb((Word32) * (p_a++), (Word32) * (p_yy1), s1); 3584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber s2 = amrnb_fxp_msu_16_by_16bb((Word32) * (p_a), (Word32) * (p_yy1--), s2); 3594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber s1 = amrnb_fxp_msu_16_by_16bb((Word32) * (p_a++), (Word32) * (p_yy1), s1); 3604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if ((UWord32)(s1 - 0xf8000000L) < 0x0fffffffL) 3634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber temp = (Word16)(s1 >> 12); 3654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else if (s1 > 0x07ffffffL) 3674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber temp = MAX_16; 3694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else 3714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber temp = MIN_16; 3734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber s2 = amrnb_fxp_msu_16_by_16bb((Word32)a[1], (Word32)temp, s2); 3764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *(p_y++) = temp; 3784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber p_yy1 = p_y; 3794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if ((UWord32)(s2 - 0xf8000000L) < 0x0fffffffL) 3814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *(p_y++) = (Word16)(s2 >> 12); 3834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else if (s2 > 0x07ffffffL) 3854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *(p_y++) = MAX_16; 3874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else 3894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *(p_y++) = MIN_16; 3914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* Update of memory if update==1 */ 3954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (update != 0) 3964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber memcpy(mem, &y[lg-M], M*sizeof(Word16)); 3984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber return; 4014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 402