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 Pathname: ./audio/gsm-amr/c/src/pstfilt.c 354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Functions: 364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Post_Filter_reset 374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Post_Filter 384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Date: 04/14/2000 404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REVISION HISTORY 434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Changed template used to PV coding template. First attempt at 454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber optimizing C code. 464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Updated file per comments gathered from Phase 2/3 review. 484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Added setting of Overflow flag in inlined code. 504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Synchronized file with UMTS version 3.2.0. Updated coding 524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber template. Removed unnecessary include files. 534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Replaced basic_op.h with the header file of the math functions 554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber used in the file. 564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Made the following changes per comments from Phase 2/3 review: 584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1. Updated copyright year. 594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2. Modified FOR loops to count down. 604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3. Fixed typecasting issue with TI C compiler. 614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4. Added "break" statement after overflow condition occurs. 624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Removed the functions pstfilt_init and pstfilt_exit. 644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber The pst_filt related structure is no longer dynamically allocated. 654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Modified code for EPOC changes where pOverflow is passed in 674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber rather than allowing overflow to be a global variable. 684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Replaced OSCL mem type functions and eliminated include 704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber files that now are chosen by OSCL definitions 714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: Replaced "int" and/or "char" with defined types. 734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Added proper casting (Word32) to some left shifting operations 744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Description: 764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber MODULE DESCRIPTION 794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber This file contains the function that performs adaptive post-filtering on the 814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber synthesized speech. It also contains the functions that initialize, reset, 824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber and exit the post-filtering function. 834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber*/ 864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; INCLUDES 894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include <string.h> 914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "pstfilt.h" 934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "typedef.h" 944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "mode.h" 954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "basicop_malloc.h" 964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "basic_op.h" 974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "weight_a.h" 984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "residu.h" 994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "copy.h" 1004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "syn_filt.h" 1014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "preemph.h" 1024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "cnst.h" 1034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 1054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; MACROS 1064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Define module specific macros here 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#define L_H 22 /* size of truncated impulse response of A(z/g1)/A(z/g2) */ 1154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 1174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; LOCAL FUNCTION DEFINITIONS 1184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Function Prototype declaration 1194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 1204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 1224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; LOCAL VARIABLE DEFINITIONS 1234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Variable declaration - defined here and used outside this module 1244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 1254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* Spectral expansion factors */ 1274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huberstatic const Word16 gamma3_MR122[M] = 1284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 1294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 22938, 16057, 11240, 7868, 5508, 1304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3856, 2699, 1889, 1322, 925 1314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber}; 1324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huberstatic const Word16 gamma3[M] = 1344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 1354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 18022, 9912, 5451, 2998, 1649, 907, 499, 274, 151, 83 1364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber}; 1374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huberstatic const Word16 gamma4_MR122[M] = 1394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 1404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 24576, 18432, 13824, 10368, 7776, 1414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5832, 4374, 3281, 2461, 1846 1424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber}; 1434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huberstatic const Word16 gamma4[M] = 1454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 1464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 22938, 16057, 11240, 7868, 5508, 3856, 2699, 1889, 1322, 925 1474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber}; 1484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 1504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; EXTERNAL FUNCTION REFERENCES 1514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Declare functions defined elsewhere and referenced in this module 1524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 1534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 1554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; EXTERNAL VARIABLES REFERENCES 1564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Declare variables used in this module but defined elsewhere 1574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 1584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* 1604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber FUNCTION NAME: Post_Filter_reset 1624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber INPUT AND OUTPUT DEFINITIONS 1644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Inputs: 1664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber state = pointer to structure of type Post_FilterState 1674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Outputs: 1694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber fields of the structure pointed to by state is initialized to zero 1704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Returns: 1724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber return_value = 0, if reset was successful; -1, otherwise (int) 1734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Global Variables Used: 1754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 1764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Local Variables Needed: 1784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 1794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber FUNCTION DESCRIPTION 1824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber This function initializes the state memory used by the Post_Filter function 1844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber to zero. 1854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REQUIREMENTS 1884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 1904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REFERENCES 1934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber pstfilt.c, UMTS GSM AMR speech codec, R99 - Version 3.2.0, March 2, 2001 1954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber PSEUDO-CODE 1984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huberint Post_Filter_reset (Post_FilterState *state) 2004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 2014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (state == (Post_FilterState *) NULL){ 2024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber fprintf(stderr, "Post_Filter_reset: invalid parameter\n"); 2034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber return -1; 2044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 2054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Set_zero (state->mem_syn_pst, M); 2074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Set_zero (state->res2, L_SUBFR); 2084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Set_zero (state->synth_buf, L_FRAME + M); 2094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber agc_reset(state->agc_state); 2104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber preemphasis_reset(state->preemph_state); 2114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber return 0; 2134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 2144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 2164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber RESOURCES USED [optional] 2174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber When the code is written for a specific target processor the 2194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber the resources used should be documented below. 2204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber HEAP MEMORY USED: x bytes 2224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber STACK MEMORY USED: x bytes 2244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CLOCK CYCLES: (cycle count equation for this function) + (variable 2264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber used to represent cycle count for each subroutine 2274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber called) 2284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber where: (cycle count variable) = cycle count for [subroutine 2294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber name] 2304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 2324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CAUTION [optional] 2334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber [State any special notes, constraints or cautions for users of this function] 2344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 2364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber*/ 2374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas HuberWord16 Post_Filter_reset(Post_FilterState *state) 2394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 2404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (state == (Post_FilterState *) NULL) 2414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 2424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /*fprintf(stderr, "Post_Filter_reset: invalid parameter\n"); */ 2434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber return(-1); 2444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 2454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber memset(state->mem_syn_pst, 0, sizeof(Word16)*M); 2474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber memset(state->res2, 0, sizeof(Word16)*L_SUBFR); 2484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber memset(state->synth_buf, 0, sizeof(Word16)*(L_FRAME + M)); 2494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber agc_reset(&(state->agc_state)); 2504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber preemphasis_reset(&(state->preemph_state)); 2514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber return(0); 2534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 2544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/****************************************************************************/ 2564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* 2584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 2594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber FUNCTION NAME: Post_Filter 2604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 2614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber INPUT AND OUTPUT DEFINITIONS 2624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Inputs: 2644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber st = pointer to a structure of type Post_FilterState 2654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber mode = AMR mode 2664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber syn = pointer to a buffer containing synthesized speech; upon 2674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber exiting this function, it will contain the post-filtered 2684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber synthesized speech 2694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Az_4 = pointer to the interpolated LPC parameters for all subframes 2704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber pOverflow = pointer to overflow indicator of type Flag 2714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Outputs: 2734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber fields of the structure pointed to by st contains the updated field 2744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber values 2754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber syn buffer contains the post-filtered synthesized speech 2764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber pOverflow = 1 if overflow occurrs in the math functions called else 2774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber it is zero. 2784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Returns: 2804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber return_value = 0 (int) 2814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Global Variables Used: 2834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 2844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Local Variables Needed: 2864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 2874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 2894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber FUNCTION DESCRIPTION 2904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber This function performs the post-filtering on the synthesized speech. The 2924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber post-filtering process is described as follows: 2934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber (1) inverse filtering of syn[] through A(z/0.7) to get res2[] 2944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber (2) tilt compensation filtering; 1 - MU*k*z^-1 2954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber (3) synthesis filtering through 1/A(z/0.75) 2964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber (4) adaptive gain control 2974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 2994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REQUIREMENTS 3004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 3024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 3044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REFERENCES 3054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber pstfilt.c, UMTS GSM AMR speech codec, R99 - Version 3.2.0, March 2, 2001 3074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 3094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber PSEUDO-CODE 3104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huberint Post_Filter ( 3124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Post_FilterState *st, // i/o : post filter states 3134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber enum Mode mode, // i : AMR mode 3144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *syn, // i/o : synthesis speech (postfiltered is output) 3154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *Az_4 // i : interpolated LPC parameters in all subfr. 3164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber) 3174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 3184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *-------------------------------------------------------------------* 3194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * Declaration of parameters * 3204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *-------------------------------------------------------------------* 3214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 Ap3[MP1], Ap4[MP1]; // bandwidth expanded LP parameters 3234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *Az; // pointer to Az_4: 3244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // LPC parameters in each subframe 3254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 i_subfr; // index for beginning of subframe 3264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 h[L_H]; 3274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 i; 3294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 temp1, temp2; 3304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 L_tmp; 3314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *syn_work = &st->synth_buf[M]; 3324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *-----------------------------------------------------* 3354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * Post filtering * 3364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *-----------------------------------------------------* 3374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Copy (syn, syn_work , L_FRAME); 3394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Az = Az_4; 3414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i_subfr = 0; i_subfr < L_FRAME; i_subfr += L_SUBFR) 3434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // Find weighted filter coefficients Ap3[] and ap[4] 3454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (sub(mode, MR122) == 0 || sub(mode, MR102) == 0) 3474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Weight_Ai (Az, gamma3_MR122, Ap3); 3494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Weight_Ai (Az, gamma4_MR122, Ap4); 3504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else 3524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Weight_Ai (Az, gamma3, Ap3); 3544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Weight_Ai (Az, gamma4, Ap4); 3554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // filtering of synthesis speech by A(z/0.7) to find res2[] 3584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Residu (Ap3, &syn_work[i_subfr], st->res2, L_SUBFR); 3604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // tilt compensation filter 3624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // impulse response of A(z/0.7)/A(z/0.75) 3644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Copy (Ap3, h, M + 1); 3664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Set_zero (&h[M + 1], L_H - M - 1); 3674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Syn_filt (Ap4, h, h, L_H, &h[M + 1], 0); 3684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // 1st correlation of h[] 3704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_mult (h[0], h[0]); 3724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 1; i < L_H; i++) 3734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_mac (L_tmp, h[i], h[i]); 3754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber temp1 = extract_h (L_tmp); 3774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_mult (h[0], h[1]); 3794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 1; i < L_H - 1; i++) 3804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_mac (L_tmp, h[i], h[i + 1]); 3824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber temp2 = extract_h (L_tmp); 3844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (temp2 <= 0) 3864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber temp2 = 0; 3884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else 3904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber temp2 = mult (temp2, MU); 3924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber temp2 = div_s (temp2, temp1); 3934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber preemphasis (st->preemph_state, st->res2, temp2, L_SUBFR); 3964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // filtering through 1/A(z/0.75) 3984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Syn_filt (Ap4, st->res2, &syn[i_subfr], L_SUBFR, st->mem_syn_pst, 1); 4004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // scale output to input 4024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber agc (st->agc_state, &syn_work[i_subfr], &syn[i_subfr], 4044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber AGC_FAC, L_SUBFR); 4054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Az += MP1; 4074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 4084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // update syn_work[] buffer 4104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Copy (&syn_work[L_FRAME - M], &syn_work[-M], M); 4124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber return 0; 4144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 4154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 4174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber RESOURCES USED [optional] 4184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber When the code is written for a specific target processor the 4204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber the resources used should be documented below. 4214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber HEAP MEMORY USED: x bytes 4234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber STACK MEMORY USED: x bytes 4254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CLOCK CYCLES: (cycle count equation for this function) + (variable 4274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber used to represent cycle count for each subroutine 4284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber called) 4294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber where: (cycle count variable) = cycle count for [subroutine 4304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber name] 4314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 4334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CAUTION [optional] 4344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber [State any special notes, constraints or cautions for users of this function] 4354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 4374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber*/ 4384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Hubervoid Post_Filter( 4404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Post_FilterState *st, /* i/o : post filter states */ 4414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber enum Mode mode, /* i : AMR mode */ 4424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *syn, /* i/o : synthesis speech (postfiltered is output) */ 4434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *Az_4, /* i : interpolated LPC parameters in all subfr. */ 4444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Flag *pOverflow 4454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber) 4464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 4474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 Ap3[MP1]; 4484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 Ap4[MP1]; /* bandwidth expanded LP parameters */ 4494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *Az; /* pointer to Az_4: */ 4504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* LPC parameters in each subframe */ 4514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber register Word16 i_subfr; /* index for beginning of subframe */ 4524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 h[L_H]; 4534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber register Word16 i; 4554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 temp1; 4564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 temp2; 4574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 L_tmp; 4584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 L_tmp2; 4594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *syn_work = &st->synth_buf[M]; 4604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /*-----------------------------------------------------* 4634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * Post filtering * 4644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *-----------------------------------------------------*/ 4654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Copy(syn, syn_work , L_FRAME); 4674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Az = Az_4; 4694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i_subfr = 0; i_subfr < L_FRAME; i_subfr += L_SUBFR) 4714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 4724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* Find weighted filter coefficients Ap3[] and ap[4] */ 4734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (mode == MR122 || mode == MR102) 4754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 4764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Weight_Ai(Az, gamma3_MR122, Ap3); 4774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Weight_Ai(Az, gamma4_MR122, Ap4); 4784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 4794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else 4804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 4814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Weight_Ai(Az, gamma3, Ap3); 4824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Weight_Ai(Az, gamma4, Ap4); 4834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 4844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* filtering of synthesis speech by A(z/0.7) to find res2[] */ 4864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Residu(Ap3, &syn_work[i_subfr], st->res2, L_SUBFR); 4884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* tilt compensation filter */ 4904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* impulse response of A(z/0.7)/A(z/0.75) */ 4924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Copy(Ap3, h, M + 1); 4944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber memset(&h[M + 1], 0, sizeof(Word16)*(L_H - M - 1)); 4954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Syn_filt(Ap4, h, h, L_H, &h[M + 1], 0); 4964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* 1st correlation of h[] */ 4984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = 0; 5004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = L_H - 1; i >= 0; i--) 5024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 5034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp2 = ((Word32) h[i]) * h[i]; 5044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (L_tmp2 != (Word32) 0x40000000L) 5064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 5074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp2 = L_tmp2 << 1; 5084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 5094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else 5104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 5114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *pOverflow = 1; 5124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp2 = MAX_32; 5134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber break; 5144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 5154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_add(L_tmp, L_tmp2, pOverflow); 5174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 5184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber temp1 = (Word16)(L_tmp >> 16); 5194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = 0; 5214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = L_H - 2; i >= 0; i--) 5234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 5244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp2 = ((Word32) h[i]) * h[i + 1]; 5254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (L_tmp2 != (Word32) 0x40000000L) 5274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 5284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp2 = L_tmp2 << 1; 5294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 5304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else 5314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 5324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *pOverflow = 1; 5334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp2 = MAX_32; 5344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber break; 5354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 5364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_add(L_tmp, L_tmp2, pOverflow); 5384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 5394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber temp2 = (Word16)(L_tmp >> 16); 5404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (temp2 <= 0) 5424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 5434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber temp2 = 0; 5444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 5454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else 5464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 5474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = (((Word32) temp2) * MU) >> 15; 5484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* Sign-extend product */ 5504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (L_tmp & (Word32) 0x00010000L) 5514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 5524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_tmp | (Word32) 0xffff0000L; 5534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 5544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber temp2 = (Word16) L_tmp; 5554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber temp2 = div_s(temp2, temp1); 5574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 5584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber preemphasis(&(st->preemph_state), st->res2, temp2, L_SUBFR, pOverflow); 5604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* filtering through 1/A(z/0.75) */ 5624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Syn_filt(Ap4, st->res2, &syn[i_subfr], L_SUBFR, st->mem_syn_pst, 1); 5644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* scale output to input */ 5664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber agc(&(st->agc_state), &syn_work[i_subfr], &syn[i_subfr], 5684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber AGC_FAC, L_SUBFR, pOverflow); 5694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Az += MP1; 5714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 5724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* update syn_work[] buffer */ 5744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Copy(&syn_work[L_FRAME - M], &syn_work[-M], M); 5764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber return; 5784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 579