1a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber/* ------------------------------------------------------------------ 2a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber * Copyright (C) 1998-2009 PacketVideo 3a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber * 4a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber * Licensed under the Apache License, Version 2.0 (the "License"); 5a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber * you may not use this file except in compliance with the License. 6a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber * You may obtain a copy of the License at 7a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber * 8a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber * http://www.apache.org/licenses/LICENSE-2.0 9a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber * 10a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber * Unless required by applicable law or agreed to in writing, software 11a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber * distributed under the License is distributed on an "AS IS" BASIS, 12a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either 13a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber * express or implied. 14a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber * See the License for the specific language governing permissions 15a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber * and limitations under the License. 16a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber * ------------------------------------------------------------------- 17a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber */ 18a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber/**************************************************************************************** 19a30d40083856cb4edd225faf8b488fab156e5976Andreas HuberPortions of this file are derived from the following 3GPP standard: 20a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 21a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 3GPP TS 26.073 22a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber ANSI-C code for the Adaptive Multi-Rate (AMR) speech codec 23a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Available from http://www.3gpp.org 24a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 25a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber(C) 2004, 3GPP Organizational Partners (ARIB, ATIS, CCSA, ETSI, TTA, TTC) 26a30d40083856cb4edd225faf8b488fab156e5976Andreas HuberPermission to distribute, modify and use this file under the standard license 27a30d40083856cb4edd225faf8b488fab156e5976Andreas Huberterms listed above has been obtained from the copyright holder. 28a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber****************************************************************************************/ 29a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber/* 30a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber------------------------------------------------------------------------------ 31a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 32a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 33a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 34a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Pathname: ./audio/gsm-amr/c/src/c_g_aver.c 35a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Functions: 36a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Cb_gain_average_reset 37a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Cb_gain_average 38a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 39a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Date: 03/28/2000 40a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 41a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber------------------------------------------------------------------------------ 42a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber REVISION HISTORY 43a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 44a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Description: Made some changes to the comments to match the comments from 45a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber other modules. 46a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 47a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Description: Made changes based on comments from the review meeting. 48a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 49a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Description: Synchronized file with UMTS version 3.2.0. Updated coding 50a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber template. 51a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 52a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Description: Made the following changes per comments from Phase 2/3 review: 53a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 1. Defined one local variable per line. 54a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 55a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Description: Removed the functions Cb_gain_average_init and 56a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Cb_gain_average_exit. The Cb_gain_average related structure is no longer 57a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber dynamically allocated. 58a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 59a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Description: Passing in pOverflow to comply with changes needed for EPOC 60a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Updated the include files for the module. 61a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 62a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Description: Changed round function name to pv_round to avoid conflict with 63a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber round function in C standard library. 64a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 65a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 66a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Description: Replaced OSCL mem type functions and eliminated include 67a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber files that now are chosen by OSCL definitions 68a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 69a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Description: 70a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber------------------------------------------------------------------------------ 71a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber MODULE DESCRIPTION 72a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 73a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber This file contains functions that reset and perform 74a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber codebook gain calculations. 75a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 76a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber------------------------------------------------------------------------------ 77a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber*/ 78a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 79a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 80a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber/*---------------------------------------------------------------------------- 81a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber; INCLUDES 82a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber----------------------------------------------------------------------------*/ 83a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber#include <string.h> 84a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 85a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber#include "c_g_aver.h" 86a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber#include "typedef.h" 87a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber#include "mode.h" 88a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber#include "cnst.h" 89a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 90a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber#include "basic_op.h" 91a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 92a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber/*---------------------------------------------------------------------------- 93a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber; MACROS 94a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber; Define module specific macros here 95a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber----------------------------------------------------------------------------*/ 96a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 97a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 98a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber/*---------------------------------------------------------------------------- 99a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber; DEFINES 100a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber; Include all pre-processor statements here. Include conditional 101a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber; compile variables also. 102a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber----------------------------------------------------------------------------*/ 103a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 104a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber/*---------------------------------------------------------------------------- 105a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber; LOCAL FUNCTION DEFINITIONS 106a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber; Function Prototype declaration 107a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber----------------------------------------------------------------------------*/ 108a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 109a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber/*---------------------------------------------------------------------------- 110a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber; LOCAL VARIABLE DEFINITIONS 111a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber; Variable declaration - defined here and used outside this module 112a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber----------------------------------------------------------------------------*/ 113a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 114a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber/* 115a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber------------------------------------------------------------------------------ 116a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber FUNCTION NAME: Cb_gain_average_reset 117a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber------------------------------------------------------------------------------ 118a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber INPUT AND OUTPUT DEFINITIONS 119a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 120a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Inputs: 121a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber state = pointer to a structure of type Cb_gain_averageState 122a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 123a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Outputs: 124a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Structure pointed to by state is initialized to zeros 125a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 126a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Returns: 127a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Returns 0 if memory was successfully initialized, 128a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber otherwise returns -1. 129a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 130a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Global Variables Used: 131a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber None. 132a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 133a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Local Variables Needed: 134a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber None. 135a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 136a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber------------------------------------------------------------------------------ 137a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber FUNCTION DESCRIPTION 138a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 139a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Resets state memory 140a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 141a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber------------------------------------------------------------------------------ 142a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber REQUIREMENTS 143a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 144a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber None. 145a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 146a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber------------------------------------------------------------------------------ 147a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber REFERENCES 148a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 149a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber c_g_aver.c, UMTS GSM AMR speech codec, R99 - Version 3.2.0, March 2, 2001 150a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 151a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber------------------------------------------------------------------------------ 152a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber PSEUDO-CODE 153a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 154a30d40083856cb4edd225faf8b488fab156e5976Andreas HuberWord16 Cb_gain_average_reset (Cb_gain_averageState *state) 155a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber{ 156a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber if (state == (Cb_gain_averageState *) NULL){ 157a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber fprintf(stderr, "Cb_gain_average_reset: invalid parameter\n"); 158a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber return -1; 159a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 160a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 161a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber // Static vectors to zero 162a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Set_zero (state->cbGainHistory, L_CBGAINHIST); 163a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 164a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber // Initialize hangover handling 165a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber state->hangVar = 0; 166a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber state->hangCount= 0; 167a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 168a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber return 0; 169a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber} 170a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 171a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber------------------------------------------------------------------------------ 172a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber RESOURCES USED [optional] 173a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 174a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber When the code is written for a specific target processor the 175a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber the resources used should be documented below. 176a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 177a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber HEAP MEMORY USED: x bytes 178a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 179a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber STACK MEMORY USED: x bytes 180a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 181a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber CLOCK CYCLES: (cycle count equation for this function) + (variable 182a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber used to represent cycle count for each subroutine 183a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber called) 184a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber where: (cycle count variable) = cycle count for [subroutine 185a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber name] 186a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 187a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber------------------------------------------------------------------------------ 188a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber CAUTION [optional] 189a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber [State any special notes, constraints or cautions for users of this function] 190a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 191a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber------------------------------------------------------------------------------ 192a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber*/ 193a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 194a30d40083856cb4edd225faf8b488fab156e5976Andreas HuberWord16 Cb_gain_average_reset(Cb_gain_averageState *state) 195a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber{ 196a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber if (state == (Cb_gain_averageState *) NULL) 197a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 198a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* fprint(stderr, "Cb_gain_average_reset: invalid parameter\n"); */ 199a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber return(-1); 200a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 201a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 202a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* Static vectors to zero */ 203a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber memset(state->cbGainHistory, 0, L_CBGAINHIST*sizeof(Word16)); 204a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 205a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* Initialize hangover handling */ 206a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber state->hangVar = 0; 207a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber state->hangCount = 0; 208a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 209a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber return(0); 210a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber} 211a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 212a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber/****************************************************************************/ 213a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 214a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber/* 215a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber------------------------------------------------------------------------------ 216a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber FUNCTION NAME: Cb_gain_average 217a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber------------------------------------------------------------------------------ 218a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber INPUT AND OUTPUT DEFINITIONS 219a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 220a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Inputs: 221a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber st = pointer to structure of type Cb_gain_averageState 222a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber mode = AMR mode (enum Mode) 223a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber gain_code = CB gain (Word16) 224a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber lsp = the LSP for the current frame (Word16) 225a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber lspAver = the average of LSP for 8 frames (Word16) 226a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber bfi = bad frame indication flag (Word16) 227a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber prev_bf = previous bad frame indication flag (Word16) 228a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber pdfi = potential degraded bad frame ind flag (Word16) 229a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber prev_pdf = prev pot. degraded bad frame ind flag (Word16) 230a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber inBackgroundNoise = background noise decision (Word16) 231a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber voicedHangover = # of frames after last voiced frame (Word16) 232a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber pOverflow = address of overflow (Flag) 233a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 234a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Returns: 235a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber cbGainMix = codebook gain (Word16) 236a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 237a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Outputs: 238a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber None. 239a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 240a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Global Variables Used: 241a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber None. 242a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 243a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Local Variables Needed: 244a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber None. 245a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 246a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber------------------------------------------------------------------------------ 247a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber FUNCTION DESCRIPTION 248a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 249a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber The mix cb gains for MR475, MR515, MR59, MR67, MR102; gain_code other modes 250a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 251a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber------------------------------------------------------------------------------ 252a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber REQUIREMENTS 253a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 254a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber None. 255a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 256a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber------------------------------------------------------------------------------ 257a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber REFERENCES 258a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 259a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber c_g_aver.c, UMTS GSM AMR speech codec, R99 - Version 3.2.0, March 2, 2001 260a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 261a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber------------------------------------------------------------------------------ 262a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber PSEUDO-CODE 263a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 264a30d40083856cb4edd225faf8b488fab156e5976Andreas HuberWord16 Cb_gain_average ( 265a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Cb_gain_averageState *st, // i/o : State variables for CB gain avergeing 266a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber enum Mode mode, // i : AMR mode 267a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 gain_code, // i : CB gain Q1 268a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 lsp[], // i : The LSP for the current frame Q15 269a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 lspAver[], // i : The average of LSP for 8 frames Q15 270a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 bfi, // i : bad frame indication flag 271a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 prev_bf, // i : previous bad frame indication flag 272a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 pdfi, // i : potential degraded bad frame ind flag 273a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 prev_pdf, // i : prev pot. degraded bad frame ind flag 274a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 inBackgroundNoise, // i : background noise decision 275a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 voicedHangover // i : # of frames after last voiced frame 276a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber ) 277a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber{ 278a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber //---------------------------------------------------------* 279a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber * Compute mixed cb gain, used to make cb gain more * 280a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber * smooth in background noise for modes 5.15, 5.9 and 6.7 * 281a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber * states that needs to be updated by all * 282a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber *--------------------------------------------------------- 283a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 i; 284a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 cbGainMix, diff, tmp_diff, bgMix, cbGainMean; 285a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word32 L_sum; 286a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 tmp[M], tmp1, tmp2, shift1, shift2, shift; 287a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 288a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber // set correct cbGainMix for MR74, MR795, MR122 289a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber cbGainMix = gain_code; 290a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 291a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber *-------------------------------------------------------* 292a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber * Store list of CB gain needed in the CB gain * 293a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber * averaging * 294a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber *-------------------------------------------------------* 295a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber for (i = 0; i < (L_CBGAINHIST-1); i++) 296a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 297a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber st->cbGainHistory[i] = st->cbGainHistory[i+1]; 298a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 299a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber st->cbGainHistory[L_CBGAINHIST-1] = gain_code; 300a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 301a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber // compute lsp difference 302a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber for (i = 0; i < M; i++) { 303a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber tmp1 = abs_s(sub(lspAver[i], lsp[i])); // Q15 304a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber shift1 = sub(norm_s(tmp1), 1); // Qn 305a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber tmp1 = shl(tmp1, shift1); // Q15+Qn 306a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber shift2 = norm_s(lspAver[i]); // Qm 307a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber tmp2 = shl(lspAver[i], shift2); // Q15+Qm 308a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber tmp[i] = div_s(tmp1, tmp2); // Q15+(Q15+Qn)-(Q15+Qm) 309a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber shift = sub(add(2, shift1), shift2); 310a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber if (shift >= 0) 311a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 312a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber tmp[i] = shr(tmp[i], shift); // Q15+Qn-Qm-Qx=Q13 313a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 314a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber else 315a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 316a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber tmp[i] = shl(tmp[i], negate(shift)); // Q15+Qn-Qm-Qx=Q13 317a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 318a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 319a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 320a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber diff = tmp[0]; 321a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber for (i = 1; i < M; i++) { 322a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber diff = add(diff, tmp[i]); // Q13 323a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 324a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 325a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber // Compute hangover 326a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber if (sub(diff, 5325) > 0) // 0.65 in Q11 327a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 328a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber st->hangVar = add(st->hangVar, 1); 329a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 330a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber else 331a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 332a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber st->hangVar = 0; 333a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 334a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 335a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber if (sub(st->hangVar, 10) > 0) 336a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 337a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber st->hangCount = 0; // Speech period, reset hangover variable 338a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 339a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 340a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber // Compute mix constant (bgMix) 341a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber bgMix = 8192; // 1 in Q13 342a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber if ((sub(mode, MR67) <= 0) || (sub(mode, MR102) == 0)) 343a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber // MR475, MR515, MR59, MR67, MR102 344a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 345a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber // if errors and presumed noise make smoothing probability stronger 346a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber if (((((pdfi != 0) && (prev_pdf != 0)) || (bfi != 0) || (prev_bf != 0)) && 347a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber (sub(voicedHangover, 1) > 0) && (inBackgroundNoise != 0) && 348a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber ((sub(mode, MR475) == 0) || 349a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber (sub(mode, MR515) == 0) || 350a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber (sub(mode, MR59) == 0)) )) 351a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 352a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber // bgMix = min(0.25, max(0.0, diff-0.55)) / 0.25; 353a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber tmp_diff = sub(diff, 4506); // 0.55 in Q13 354a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 355a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber // max(0.0, diff-0.55) 356a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber if (tmp_diff > 0) 357a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 358a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber tmp1 = tmp_diff; 359a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 360a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber else 361a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 362a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber tmp1 = 0; 363a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 364a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 365a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber // min(0.25, tmp1) 366a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber if (sub(2048, tmp1) < 0) 367a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 368a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber bgMix = 8192; 369a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 370a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber else 371a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 372a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber bgMix = shl(tmp1, 2); 373a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 374a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 375a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber else 376a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 377a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber // bgMix = min(0.25, max(0.0, diff-0.40)) / 0.25; 378a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber tmp_diff = sub(diff, 3277); // 0.4 in Q13 379a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 380a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber // max(0.0, diff-0.40) 381a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber if (tmp_diff > 0) 382a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 383a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber tmp1 = tmp_diff; 384a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 385a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber else 386a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 387a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber tmp1 = 0; 388a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 389a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 390a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber // min(0.25, tmp1) 391a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber if (sub(2048, tmp1) < 0) 392a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 393a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber bgMix = 8192; 394a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 395a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber else 396a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 397a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber bgMix = shl(tmp1, 2); 398a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 399a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 400a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 401a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber if ((sub(st->hangCount, 40) < 0) || (sub(diff, 5325) > 0)) // 0.65 in Q13 402a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 403a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber bgMix = 8192; // disable mix if too short time since 404a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 405a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 406a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber // Smoothen the cb gain trajectory 407a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber // smoothing depends on mix constant bgMix 408a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber L_sum = L_mult(6554, st->cbGainHistory[2]); // 0.2 in Q15; L_sum in Q17 409a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber for (i = 3; i < L_CBGAINHIST; i++) 410a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 411a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber L_sum = L_mac(L_sum, 6554, st->cbGainHistory[i]); 412a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 413a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber cbGainMean = pv_round(L_sum); // Q1 414a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 415a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber // more smoothing in error and bg noise (NB no DFI used here) 416a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber if (((bfi != 0) || (prev_bf != 0)) && (inBackgroundNoise != 0) && 417a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber ((sub(mode, MR475) == 0) || 418a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber (sub(mode, MR515) == 0) || 419a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber (sub(mode, MR59) == 0)) ) 420a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 421a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber L_sum = L_mult(4681, st->cbGainHistory[0]); // 0.143 in Q15; L_sum in Q17 422a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber for (i = 1; i < L_CBGAINHIST; i++) 423a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 424a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber L_sum = L_mac(L_sum, 4681, st->cbGainHistory[i]); 425a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 426a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber cbGainMean = pv_round(L_sum); // Q1 427a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 428a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 429a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber // cbGainMix = bgMix*cbGainMix + (1-bgMix)*cbGainMean; 430a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber L_sum = L_mult(bgMix, cbGainMix); // L_sum in Q15 431a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber L_sum = L_mac(L_sum, 8192, cbGainMean); 432a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber L_sum = L_msu(L_sum, bgMix, cbGainMean); 433a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber cbGainMix = pv_round(L_shl(L_sum, 2)); // Q1 434a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 435a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 436a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber st->hangCount = add(st->hangCount, 1); 437a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber return cbGainMix; 438a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber} 439a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 440a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber------------------------------------------------------------------------------ 441a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber RESOURCES USED [optional] 442a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 443a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber When the code is written for a specific target processor the 444a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber the resources used should be documented below. 445a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 446a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber HEAP MEMORY USED: x bytes 447a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 448a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber STACK MEMORY USED: x bytes 449a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 450a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber CLOCK CYCLES: (cycle count equation for this function) + (variable 451a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber used to represent cycle count for each subroutine 452a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber called) 453a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber where: (cycle count variable) = cycle count for [subroutine 454a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber name] 455a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 456a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber------------------------------------------------------------------------------ 457a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber CAUTION [optional] 458a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber [State any special notes, constraints or cautions for users of this function] 459a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 460a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber------------------------------------------------------------------------------ 461a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber*/ 462a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 463a30d40083856cb4edd225faf8b488fab156e5976Andreas HuberWord16 Cb_gain_average( 464a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Cb_gain_averageState *st, /* i/o : State variables for CB gain averaging */ 465a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber enum Mode mode, /* i : AMR mode */ 466a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 gain_code, /* i : CB gain Q1 */ 467a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 lsp[], /* i : The LSP for the current frame Q15 */ 468a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 lspAver[], /* i : The average of LSP for 8 frames Q15 */ 469a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 bfi, /* i : bad frame indication flag */ 470a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 prev_bf, /* i : previous bad frame indication flag */ 471a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 pdfi, /* i : potential degraded bad frame ind flag */ 472a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 prev_pdf, /* i : prev pot. degraded bad frame ind flag */ 473a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 inBackgroundNoise, /* i : background noise decision */ 474a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 voicedHangover, /* i : # of frames after last voiced frame */ 475a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Flag *pOverflow 476a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber) 477a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber{ 478a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 i; 479a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 cbGainMix; 480a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 diff; 481a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 tmp_diff; 482a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 bgMix; 483a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 cbGainMean; 484a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word32 L_sum; 485a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 tmp[M]; 486a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 tmp1; 487a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 tmp2; 488a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 shift1; 489a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 shift2; 490a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber Word16 shift; 491a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 492a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /*---------------------------------------------------------* 493a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber * Compute mixed cb gain, used to make cb gain more * 494a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber * smooth in background noise for modes 5.15, 5.9 and 6.7 * 495a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber * states that needs to be updated by all * 496a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber *---------------------------------------------------------*/ 497a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 498a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* set correct cbGainMix for MR74, MR795, MR122 */ 499a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber cbGainMix = gain_code; 500a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 501a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /*-------------------------------------------------------* 502a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber * Store list of CB gain needed in the CB gain * 503a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber * averaging * 504a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber *-------------------------------------------------------*/ 505a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber for (i = 0; i < (L_CBGAINHIST - 1); i++) 506a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 507a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber st->cbGainHistory[i] = st->cbGainHistory[i+1]; 508a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 509a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber st->cbGainHistory[L_CBGAINHIST-1] = gain_code; 510a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 511a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber diff = 0; 512a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 513a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* compute lsp difference */ 514a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber for (i = 0; i < M; i++) 515a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 516a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber tmp1 = abs_s(sub(*(lspAver + i), *(lsp + i), pOverflow)); 517a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* Q15 */ 518a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber shift1 = sub(norm_s(tmp1), 1, pOverflow); /* Qn */ 519a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber tmp1 = shl(tmp1, shift1, pOverflow); /* Q15+Qn */ 520a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber shift2 = norm_s(*(lspAver + i)); /* Qm */ 521a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber tmp2 = shl(*(lspAver + i), shift2, pOverflow); /* Q15+Qm */ 522a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber tmp[i] = div_s(tmp1, tmp2); /* Q15+(Q15+Qn)-(Q15+Qm) */ 523a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 524a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber shift = 2 + shift1 - shift2; 525a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 526a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber if (shift >= 0) 527a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 528a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber *(tmp + i) = shr(*(tmp + i), shift, pOverflow); 529a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* Q15+Qn-Qm-Qx=Q13 */ 530a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 531a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber else 532a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 533a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber *(tmp + i) = shl(*(tmp + i), negate(shift), pOverflow); 534a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* Q15+Qn-Qm-Qx=Q13 */ 535a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 536a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 537a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber diff = add(diff, *(tmp + i), pOverflow); /* Q13 */ 538a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 539a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 540a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* Compute hangover */ 541a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 542a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber if (diff > 5325) /* 0.65 in Q11 */ 543a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 544a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber st->hangVar += 1; 545a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 546a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber else 547a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 548a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber st->hangVar = 0; 549a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 550a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 551a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 552a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber if (st->hangVar > 10) 553a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 554a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* Speech period, reset hangover variable */ 555a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber st->hangCount = 0; 556a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 557a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 558a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* Compute mix constant (bgMix) */ 559a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber bgMix = 8192; /* 1 in Q13 */ 560a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 561a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber if ((mode <= MR67) || (mode == MR102)) 562a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* MR475, MR515, MR59, MR67, MR102 */ 563a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 564a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* if errors and presumed noise make smoothing probability stronger */ 565a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 566a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber if (((((pdfi != 0) && (prev_pdf != 0)) || (bfi != 0) || 567a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber (prev_bf != 0)) 568a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber && (voicedHangover > 1) 569a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber && (inBackgroundNoise != 0) 570a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber && ((mode == MR475) || (mode == MR515) || 571a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber (mode == MR59)))) 572a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 573a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* bgMix = min(0.25, max(0.0, diff-0.55)) / 0.25; */ 574a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber tmp_diff = sub(diff, 4506, pOverflow); /* 0.55 in Q13 */ 575a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 576a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber else 577a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 578a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* bgMix = min(0.25, max(0.0, diff-0.40)) / 0.25; */ 579a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber tmp_diff = sub(diff, 3277, pOverflow); /* 0.4 in Q13 */ 580a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 581a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 582a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* max(0.0, diff-0.55) or */ 583a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* max(0.0, diff-0.40) */ 584a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber if (tmp_diff > 0) 585a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 586a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber tmp1 = tmp_diff; 587a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 588a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber else 589a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 590a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber tmp1 = 0; 591a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 592a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 593a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* min(0.25, tmp1) */ 594a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber if (2048 < tmp1) 595a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 596a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber bgMix = 8192; 597a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 598a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber else 599a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 600a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber bgMix = shl(tmp1, 2, pOverflow); 601a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 602a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 603a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber if ((st->hangCount < 40) || (diff > 5325)) /* 0.65 in Q13 */ 604a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 605a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* disable mix if too short time since */ 606a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber bgMix = 8192; 607a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 608a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 609a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* Smoothen the cb gain trajectory */ 610a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* smoothing depends on mix constant bgMix */ 611a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber L_sum = L_mult(6554, st->cbGainHistory[2], pOverflow); 612a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* 0.2 in Q15; L_sum in Q17 */ 613a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 614a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber for (i = 3; i < L_CBGAINHIST; i++) 615a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 616a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber L_sum = L_mac(L_sum, 6554, st->cbGainHistory[i], pOverflow); 617a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 618a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber cbGainMean = pv_round(L_sum, pOverflow); /* Q1 */ 619a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 620a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* more smoothing in error and bg noise (NB no DFI used here) */ 621a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 622a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber if (((bfi != 0) || (prev_bf != 0)) && (inBackgroundNoise != 0) 623a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber && ((mode == MR475) || (mode == MR515) 624a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber || (mode == MR59))) 625a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 626a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* 0.143 in Q15; L_sum in Q17 */ 627a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber L_sum = L_mult(4681, st->cbGainHistory[0], pOverflow); 628a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber for (i = 1; i < L_CBGAINHIST; i++) 629a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber { 630a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber L_sum = 631a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber L_mac(L_sum, 4681, st->cbGainHistory[i], pOverflow); 632a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 633a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber cbGainMean = pv_round(L_sum, pOverflow); /* Q1 */ 634a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 635a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 636a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* cbGainMix = bgMix*cbGainMix + (1-bgMix)*cbGainMean; */ 637a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber /* L_sum in Q15 */ 638a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber L_sum = L_mult(bgMix, cbGainMix, pOverflow); 639a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber L_sum = L_mac(L_sum, 8192, cbGainMean, pOverflow); 640a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber L_sum = L_msu(L_sum, bgMix, cbGainMean, pOverflow); 641a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber cbGainMix = pv_round(L_shl(L_sum, 2, pOverflow), pOverflow); /* Q1 */ 642a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber } 643a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 644a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber st->hangCount += 1; 645a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 646a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber return (cbGainMix); 647a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber} 648a30d40083856cb4edd225faf8b488fab156e5976Andreas Huber 649