14f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* ------------------------------------------------------------------ 24f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * Copyright (C) 1998-2009 PacketVideo 34f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * 44f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * Licensed under the Apache License, Version 2.0 (the "License"); 54f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * you may not use this file except in compliance with the License. 64f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * You may obtain a copy of the License at 74f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * 84f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * http://www.apache.org/licenses/LICENSE-2.0 94f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * 104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * Unless required by applicable law or agreed to in writing, software 114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * distributed under the License is distributed on an "AS IS" BASIS, 124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either 134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * express or implied. 144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * See the License for the specific language governing permissions 154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * and limitations under the License. 164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * ------------------------------------------------------------------- 174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber */ 184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/**************************************************************************************** 194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas HuberPortions of this file are derived from the following 3GPP standard: 204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3GPP TS 26.073 224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ANSI-C code for the Adaptive Multi-Rate (AMR) speech codec 234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Available from http://www.3gpp.org 244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber(C) 2004, 3GPP Organizational Partners (ARIB, ATIS, CCSA, ETSI, TTA, TTC) 264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas HuberPermission to distribute, modify and use this file under the standard license 274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huberterms listed above has been obtained from the copyright holder. 284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber****************************************************************************************/ 294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* 304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Pathname: ./audio/gsm-amr/c/src/gc_pred.c 324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Functions: 334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber gc_pred_reset 344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber gc_pred 354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber gc_pred_update 364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber gc_pred_average_limited 374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber MODULE DESCRIPTION 404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber This file contains the functions that perform codebook gain MA prediction. 424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber*/ 454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; INCLUDES 484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "gc_pred.h" 504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "basicop_malloc.h" 514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "basic_op.h" 524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "cnst.h" 534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#include "log2.h" 544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; MACROS 574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Define module specific macros here 584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; DEFINES 624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Include all pre-processor statements here. Include conditional 634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; compile variables also. 644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#define NPRED 4 /* number of prediction taps */ 664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* average innovation energy. */ 684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* MEAN_ENER = 36.0/constant, constant = 20*Log10(2) */ 694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#define MEAN_ENER_MR122 783741L /* 36/(20*log10(2)) (Q17) */ 704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* minimum quantized energy: -14 dB */ 724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#define MIN_ENERGY -14336 /* 14 Q10 */ 734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber#define MIN_ENERGY_MR122 -2381 /* 14 / (20*log10(2)) Q10 */ 744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; LOCAL FUNCTION DEFINITIONS 774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Function Prototype declaration 784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/*---------------------------------------------------------------------------- 814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; LOCAL VARIABLE DEFINITIONS 824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber; Variable declaration - defined here and used outside this module 834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber----------------------------------------------------------------------------*/ 844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* MA prediction coefficients (Q13) */ 864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huberstatic const Word16 pred[NPRED] = {5571, 4751, 2785, 1556}; 874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* MA prediction coefficients (Q6) */ 894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huberstatic const Word16 pred_MR122[NPRED] = {44, 37, 22, 12}; 904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* 924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber FUNCTION NAME: gc_pred_reset 944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber INPUT AND OUTPUT DEFINITIONS 964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Inputs: 984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber state = pointer to a structure of type gc_predState 994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Outputs: 1014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber past_qua_en field in the structure pointed to by state is initialized 1024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber to MIN_ENERGY 1034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber past_qua_en_MR122 field in the structure pointed to by state is 1044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber initialized to MIN_ENERGY_MR122 1054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Returns: 1074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber return_value = 0, if reset was successful; -1, otherwise (int) 1084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Global Variables Used: 1104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 1114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Local Variables Needed: 1134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 1144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber FUNCTION DESCRIPTION 1174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber This function initializes the state memory used by gc_pred to zero. 1194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REQUIREMENTS 1224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 1244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REFERENCES 1274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber gc_pred.c, UMTS GSM AMR speech codec, R99 - Version 3.2.0, March 2, 2001 1294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber PSEUDO-CODE 1324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huberint gc_pred_reset (gc_predState *state) 1344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 1354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 i; 1364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (state == (gc_predState *) NULL){ 1384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber fprintf(stderr, "gc_pred_reset: invalid parameter\n"); 1394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber return -1; 1404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 1414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for(i = 0; i < NPRED; i++) 1434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 1444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber state->past_qua_en[i] = MIN_ENERGY; 1454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber state->past_qua_en_MR122[i] = MIN_ENERGY_MR122; 1464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 1474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber return 0; 1484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 1494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber RESOURCES USED [optional] 1524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber When the code is written for a specific target processor the 1544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber the resources used should be documented below. 1554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber HEAP MEMORY USED: x bytes 1574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber STACK MEMORY USED: x bytes 1594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CLOCK CYCLES: (cycle count equation for this function) + (variable 1614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber used to represent cycle count for each subroutine 1624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber called) 1634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber where: (cycle count variable) = cycle count for [subroutine 1644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber name] 1654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CAUTION [optional] 1684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber [State any special notes, constraints or cautions for users of this function] 1694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber*/ 1724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas HuberWord16 gc_pred_reset(gc_predState *state) 1744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 1754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 i; 1764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (state == (gc_predState *) NULL) 1784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 1794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* fprintf(stderr, "gc_pred_reset: invalid parameter\n"); */ 1804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber return -1; 1814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 1824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 0; i < NPRED; i++) 1844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 1854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber state->past_qua_en[i] = MIN_ENERGY; 1864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber state->past_qua_en_MR122[i] = MIN_ENERGY_MR122; 1874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 1884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber return(0); 1904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 1914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/****************************************************************************/ 1934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 1944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* 1954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber FUNCTION NAME: gc_pred 1974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 1984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber INPUT AND OUTPUT DEFINITIONS 1994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Inputs: 2014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber st = pointer to a structure of type gc_predState 2024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber mode = AMR mode (enum Mode) 2034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber code = pointer to the innovative codebook vector; Q12 in MR122 mode, 2044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber otherwise, Q13 (Word16) 2054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber exp_gcode0 = pointer to the exponent part of predicted gain factor 2064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber (Q0) (Word16) 2074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber frac_gcode0 = pointer to the fractional part of predicted gain factor 2084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber (Q15) (Word16) 2094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber exp_en = pointer to the exponent part of the innovation energy; this 2104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber is calculated for MR795 mode, Q0 (Word16) 2114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber frac_en = pointer to the fractional part of the innovation energy; 2124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber this is calculated for MR795 mode, Q15 (Word16) 2134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber pOverflow = pointer to overflow (Flag) 2144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Outputs: 2164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber store pointed to by exp_gcode0 contains the exponent part of the 2174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber recently calculated predicted gain factor 2184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber store pointed to by frac_gcode0 contains the fractional part of the 2194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber recently calculated predicted gain factor 2204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber store pointed to by exp_en contains the exponent part of the 2214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber recently calculated innovation energy 2224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber store pointed to by frac_en contains the fractional part of the 2234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber recently calculated innovation energy 2244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber pOverflow = 1 if the math functions called by gc_pred 2254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber results in overflow else zero. 2264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Returns: 2284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 2294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Global Variables Used: 2314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 2324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Local Variables Needed: 2344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber pred = table of MA prediction coefficients (Q13) (Word16) 2354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber pred_MR122 = table of MA prediction coefficients (Q6) (Word16) 2364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 2384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber FUNCTION DESCRIPTION 2394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber This function performs the MA prediction of the innovation energy (in 2414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber dB/(20*log10(2))), with the mean removed. 2424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 2444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REQUIREMENTS 2454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 2474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 2494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REFERENCES 2504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber gc_pred.c, UMTS GSM AMR speech codec, R99 - Version 3.2.0, March 2, 2001 2524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 2544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber PSEUDO-CODE 2554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas HuberThe original etsi reference code uses a global flag Overflow. However, in the 2574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huberactual implementation a pointer to a the overflow flag is passed in. 2584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Hubervoid 2604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Hubergc_pred( 2614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber gc_predState *st, // i/o: State struct 2624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber enum Mode mode, // i : AMR mode 2634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *code, // i : innovative codebook vector (L_SUBFR) 2644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // MR122: Q12, other modes: Q13 2654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *exp_gcode0, // o : exponent of predicted gain factor, Q0 2664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *frac_gcode0,// o : fraction of predicted gain factor Q15 2674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *exp_en, // o : exponent of innovation energy, Q0 2684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // (only calculated for MR795) 2694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *frac_en // o : fraction of innovation energy, Q15 2704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // (only calculated for MR795) 2714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber) 2724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 2734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 i; 2744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 ener_code; 2754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 exp, frac; 2764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *-------------------------------------------------------------------* 2784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * energy of code: * 2794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * ~~~~~~~~~~~~~~~ * 2804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * ener_code = sum(code[i]^2) * 2814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *-------------------------------------------------------------------* 2824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ener_code = L_mac((Word32) 0, code[0], code[0]); 2834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // MR122: Q12*Q12 -> Q25 2844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // others: Q13*Q13 -> Q27 2854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 1; i < L_SUBFR; i++) 2864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ener_code = L_mac(ener_code, code[i], code[i]); 2874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (sub (mode, MR122) == 0) 2894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 2904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 ener; 2914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // ener_code = ener_code / lcode; lcode = 40; 1/40 = 26214 Q20 2934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ener_code = L_mult (pv_round (ener_code), 26214); // Q9 * Q20 -> Q30 2944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 2954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *-------------------------------------------------------------------* 2964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * energy of code: * 2974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * ~~~~~~~~~~~~~~~ * 2984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * ener_code(Q17) = 10 * Log10(energy) / constant * 2994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * = 1/2 * Log2(energy) * 3004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * constant = 20*Log10(2) * 3014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *-------------------------------------------------------------------* 3024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // ener_code = 1/2 * Log2(ener_code); Note: Log2=log2+30 3034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Log2(ener_code, &exp, &frac); 3044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ener_code = L_Comp (sub (exp, 30), frac); // Q16 for log() 3054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // ->Q17 for 1/2 log() 3064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *-------------------------------------------------------------------* 3084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * predicted energy: * 3094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * ~~~~~~~~~~~~~~~~~ * 3104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * ener(Q24) = (Emean + sum{pred[i]*past_en[i]})/constant * 3114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * = MEAN_ENER + sum(pred[i]*past_qua_en[i]) * 3124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * constant = 20*Log10(2) * 3134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *-------------------------------------------------------------------* 3144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ener = MEAN_ENER_MR122; // Q24 (Q17) 3164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 0; i < NPRED; i++) 3174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ener = L_mac (ener, st->past_qua_en_MR122[i], pred_MR122[i]); 3194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // Q10 * Q13 -> Q24 3204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // Q10 * Q6 -> Q17 3214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *-------------------------------------------------------------------* 3244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * predicted codebook gain * 3254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * ~~~~~~~~~~~~~~~~~~~~~~~ * 3264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * gc0 = Pow10( (ener*constant - ener_code*constant) / 20 ) * 3274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * = Pow2(ener-ener_code) * 3284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * = Pow2(int(d)+frac(d)) * 3294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * * 3304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * (store exp and frac for pow2()) * 3314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *-------------------------------------------------------------------* 3324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ener = L_shr (L_sub (ener, ener_code), 1); // Q16 3344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_Extract(ener, exp_gcode0, frac_gcode0); 3354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else // all modes except 12.2 3374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 L_tmp; 3394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 exp_code, gcode0; 3404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *-----------------------------------------------------------------* 3424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * Compute: means_ener - 10log10(ener_code/ L_sufr) * 3434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *-----------------------------------------------------------------* 3444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber exp_code = norm_l (ener_code); 3464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ener_code = L_shl (ener_code, exp_code); 3474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // Log2 = log2 + 27 3494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Log2_norm (ener_code, exp_code, &exp, &frac); 3504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // fact = 10/log2(10) = 3.01 = 24660 Q13 3524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = Mpy_32_16(exp, frac, -24660); // Q0.Q15 * Q13 -> Q14 3534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * L_tmp = means_ener - 10log10(ener_code/L_SUBFR) 3554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * = means_ener - 10log10(ener_code) + 10log10(L_SUBFR) 3564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * = K - fact * Log2(ener_code) 3574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * = K - fact * log2(ener_code) - fact*27 3584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * 3594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * ==> K = means_ener + fact*27 + 10log10(L_SUBFR) 3604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * 3614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * means_ener = 33 = 540672 Q14 (MR475, MR515, MR59) 3624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * means_ener = 28.75 = 471040 Q14 (MR67) 3634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * means_ener = 30 = 491520 Q14 (MR74) 3644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * means_ener = 36 = 589824 Q14 (MR795) 3654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * means_ener = 33 = 540672 Q14 (MR102) 3664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * 10log10(L_SUBFR) = 16.02 = 262481.51 Q14 3674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * fact * 27 = 1331640 Q14 3684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * ----------------------------------------- 3694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * (MR475, MR515, MR59) K = 2134793.51 Q14 ~= 16678 * 64 * 2 3704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * (MR67) K = 2065161.51 Q14 ~= 32268 * 32 * 2 3714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * (MR74) K = 2085641.51 Q14 ~= 32588 * 32 * 2 3724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * (MR795) K = 2183945.51 Q14 ~= 17062 * 64 * 2 3734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * (MR102) K = 2134793.51 Q14 ~= 16678 * 64 * 2 3744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (sub (mode, MR102) == 0) 3774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // mean = 33 dB 3794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_mac(L_tmp, 16678, 64); // Q14 3804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else if (sub (mode, MR795) == 0) 3824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // ener_code = <xn xn> * 2^27*2^exp_code 3844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // frac_en = ener_code / 2^16 3854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // = <xn xn> * 2^11*2^exp_code 3864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // <xn xn> = <xn xn>*2^11*2^exp * 2^exp_en 3874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // := frac_en * 2^exp_en 3884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // ==> exp_en = -11-exp_code; 3904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *frac_en = extract_h (ener_code); 3924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *exp_en = sub (-11, exp_code); 3934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 3944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // mean = 36 dB 3954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_mac(L_tmp, 17062, 64); // Q14 3964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 3974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else if (sub (mode, MR74) == 0) 3984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 3994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // mean = 30 dB 4004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_mac(L_tmp, 32588, 32); // Q14 4014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 4024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else if (sub (mode, MR67) == 0) 4034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 4044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // mean = 28.75 dB 4054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_mac(L_tmp, 32268, 32); // Q14 4064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 4074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else // MR59, MR515, MR475 4084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 4094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // mean = 33 dB 4104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_mac(L_tmp, 16678, 64); // Q14 4114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 4124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *-----------------------------------------------------------------* 4144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * Compute gcode0. * 4154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * = Sum(i=0,3) pred[i]*past_qua_en[i] - ener_code + mean_ener * 4164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *-----------------------------------------------------------------* 4174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_shl(L_tmp, 10); // Q24 4194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 0; i < 4; i++) 4204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_mac(L_tmp, pred[i], st->past_qua_en[i]); 4214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // Q13 * Q10 -> Q24 4224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber gcode0 = extract_h(L_tmp); // Q8 4244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *-----------------------------------------------------------------* 4264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * gcode0 = pow(10.0, gcode0/20) * 4274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * = pow(2, 3.3219*gcode0/20) * 4284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * = pow(2, 0.166*gcode0) * 4294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *-----------------------------------------------------------------* 4304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // 5439 Q15 = 0.165985 4324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // (correct: 1/(20*log10(2)) 0.166096 = 5443 Q15) 4334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (sub (mode, MR74) == 0) // For IS641 bitexactness 4344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_mult(gcode0, 5439); // Q8 * Q15 -> Q24 4354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else 4364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_mult(gcode0, 5443); // Q8 * Q15 -> Q24 4374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_shr(L_tmp, 8); // -> Q16 4394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_Extract(L_tmp, exp_gcode0, frac_gcode0); // -> Q0.Q15 4404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 4414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 4424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 4444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber RESOURCES USED [optional] 4454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber When the code is written for a specific target processor the 4474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber the resources used should be documented below. 4484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber HEAP MEMORY USED: x bytes 4504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber STACK MEMORY USED: x bytes 4524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CLOCK CYCLES: (cycle count equation for this function) + (variable 4544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber used to represent cycle count for each subroutine 4554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber called) 4564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber where: (cycle count variable) = cycle count for [subroutine 4574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber name] 4584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 4604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CAUTION [optional] 4614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber [State any special notes, constraints or cautions for users of this function] 4624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 4644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber*/ 4654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Hubervoid gc_pred( 4674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber gc_predState *st, /* i/o: State struct */ 4684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber enum Mode mode, /* i : AMR mode */ 4694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *code, /* i : innovative codebook vector (L_SUBFR) */ 4704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* MR122: Q12, other modes: Q13 */ 4714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *exp_gcode0, /* o : exponent of predicted gain factor, Q0 */ 4724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *frac_gcode0,/* o : fraction of predicted gain factor Q15 */ 4734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *exp_en, /* o : exponent of innovation energy, Q0 */ 4744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* (only calculated for MR795) */ 4754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *frac_en, /* o : fraction of innovation energy, Q15 */ 4764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* (only calculated for MR795) */ 4774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Flag *pOverflow 4784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber) 4794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 4804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber register Word16 i; 4814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber register Word32 L_temp1, L_temp2; 4824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber register Word32 L_tmp; 4834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 ener_code; 4844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word32 ener; 4854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 exp, frac; 4864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 exp_code, gcode0; 4874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 tmp; 4884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *p_code = &code[0]; 4894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /*-------------------------------------------------------------------* 4914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * energy of code: * 4924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * ~~~~~~~~~~~~~~~ * 4934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * ener_code = sum(code[i]^2) * 4944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *-------------------------------------------------------------------*/ 4954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ener_code = 0; 4964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 4974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* MR122: Q12*Q12 -> Q25 */ 4984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* others: Q13*Q13 -> Q27 */ 4994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = L_SUBFR >> 2; i != 0; i--) 5014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 5024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp = *(p_code++); 5034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ener_code += ((Word32) tmp * tmp) >> 3; 5044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp = *(p_code++); 5054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ener_code += ((Word32) tmp * tmp) >> 3; 5064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp = *(p_code++); 5074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ener_code += ((Word32) tmp * tmp) >> 3; 5084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber tmp = *(p_code++); 5094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ener_code += ((Word32) tmp * tmp) >> 3; 5104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 5114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ener_code <<= 4; 5134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (ener_code < 0) /* Check for saturation */ 5154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 5164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ener_code = MAX_32; 5174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 5184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (mode == MR122) 5204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 5214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* ener_code = ener_code / lcode; lcode = 40; 1/40 = 26214 Q20 */ 5224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* Q9 * Q20 -> Q30 */ 5234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ener_code = ((Word32)(pv_round(ener_code, pOverflow) * 26214)) << 1; 5254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /*-------------------------------------------------------------* 5274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * energy of code: * 5284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * ~~~~~~~~~~~~~~~ * 5294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * ener_code(Q17) = 10 * Log10(energy) / constant * 5304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * = 1/2 * Log2(energy) * 5314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * constant = 20*Log10(2) * 5324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *-------------------------------------------------------------*/ 5334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* ener_code = 1/2 * Log2(ener_code); Note: Log2=log2+30 */ 5344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Log2(ener_code, &exp, &frac, pOverflow); 5354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* Q16 for log() */ 5374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* ->Q17 for 1/2 log()*/ 5384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_temp1 = (Word32)(exp - 30) << 16; 5404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ener_code = L_temp1 + ((Word32)frac << 1); 5414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /*-------------------------------------------------------------* 5434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * predicted energy: * 5444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * ~~~~~~~~~~~~~~~~~ * 5454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * ener(Q24) = (Emean + sum{pred[i]*past_en[i]})/constant * 5464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * = MEAN_ENER + sum(pred[i]*past_qua_en[i]) * 5474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * constant = 20*Log10(2) * 5484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *-------------------------------------------------------------*/ 5494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ener = MEAN_ENER_MR122; /* Q24 (Q17) */ 5514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 0; i < NPRED; i++) 5524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 5534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_temp1 = (((Word32) st->past_qua_en_MR122[i]) * 5544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber pred_MR122[i]) << 1; 5554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ener = L_add(ener, L_temp1, pOverflow); 5564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* Q10 * Q13 -> Q24 */ 5584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* Q10 * Q6 -> Q17 */ 5594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 5604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /*---------------------------------------------------------------* 5624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * predicted codebook gain * 5634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * ~~~~~~~~~~~~~~~~~~~~~~~ * 5644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * gc0 = Pow10( (ener*constant - ener_code*constant) / 20 ) * 5654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * = Pow2(ener-ener_code) * 5664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * = Pow2(int(d)+frac(d)) * 5674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * * 5684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * (store exp and frac for pow2()) * 5694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *---------------------------------------------------------------*/ 5704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* Q16 */ 5714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_temp1 = L_sub(ener, ener_code, pOverflow); 5734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *exp_gcode0 = (Word16)(L_temp1 >> 17); 5764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_temp2 = (Word32) * exp_gcode0 << 15; 5784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_temp1 >>= 2; 5794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *frac_gcode0 = (Word16)(L_temp1 - L_temp2); 5814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 5834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else /* all modes except 12.2 */ 5844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 5854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /*-----------------------------------------------------------------* 5864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * Compute: means_ener - 10log10(ener_code/ L_sufr) * 5874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *-----------------------------------------------------------------*/ 5884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber exp_code = norm_l(ener_code); 5904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ener_code = L_shl(ener_code, exp_code, pOverflow); 5914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* Log2 = log2 + 27 */ 5934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Log2_norm(ener_code, exp_code, &exp, &frac); 5944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* fact = 10/log2(10) = 3.01 = 24660 Q13 */ 5964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* Q0.Q15 * Q13 -> Q14 */ 5974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 5984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_temp2 = (((Word32) exp) * -24660) << 1; 5994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = (((Word32) frac) * -24660) >> 15; 6004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 6014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* Sign-extend resulting product */ 6024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (L_tmp & (Word32) 0x00010000L) 6034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 6044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_tmp | (Word32) 0xffff0000L; 6054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 6064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 6074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_tmp << 1; 6084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_add(L_tmp, L_temp2, pOverflow); 6094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 6104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 6114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* L_tmp = means_ener - 10log10(ener_code/L_SUBFR) 6124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * = means_ener - 10log10(ener_code) + 10log10(L_SUBFR) 6134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * = K - fact * Log2(ener_code) 6144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * = K - fact * log2(ener_code) - fact*27 6154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * 6164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * ==> K = means_ener + fact*27 + 10log10(L_SUBFR) 6174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * 6184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * means_ener = 33 = 540672 Q14 (MR475, MR515, MR59) 6194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * means_ener = 28.75 = 471040 Q14 (MR67) 6204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * means_ener = 30 = 491520 Q14 (MR74) 6214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * means_ener = 36 = 589824 Q14 (MR795) 6224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * means_ener = 33 = 540672 Q14 (MR102) 6234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * 10log10(L_SUBFR) = 16.02 = 262481.51 Q14 6244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * fact * 27 = 1331640 Q14 6254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * ----------------------------------------- 6264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * (MR475, MR515, MR59) K = 2134793.51 Q14 ~= 16678 * 64 * 2 6274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * (MR67) K = 2065161.51 Q14 ~= 32268 * 32 * 2 6284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * (MR74) K = 2085641.51 Q14 ~= 32588 * 32 * 2 6294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * (MR795) K = 2183945.51 Q14 ~= 17062 * 64 * 2 6304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * (MR102) K = 2134793.51 Q14 ~= 16678 * 64 * 2 6314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber */ 6324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 6334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (mode == MR102) 6344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 6354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* mean = 33 dB */ 6364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_temp2 = (Word32) 16678 << 7; 6374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_add(L_tmp, L_temp2, pOverflow); /* Q14 */ 6384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 6394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else if (mode == MR795) 6404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 6414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* ener_code = <xn xn> * 2^27*2^exp_code 6424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber frac_en = ener_code / 2^16 6434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber = <xn xn> * 2^11*2^exp_code 6444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber <xn xn> = <xn xn>*2^11*2^exp * 2^exp_en 6454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber : = frac_en * 2^exp_en 6464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ==> exp_en = -11-exp_code; */ 6474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *frac_en = (Word16)(ener_code >> 16); 6484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *exp_en = sub(-11, exp_code, pOverflow); 6494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 6504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* mean = 36 dB */ 6514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_temp2 = (Word32) 17062 << 7; 6524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_add(L_tmp, L_temp2, pOverflow); /* Q14 */ 6534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 6544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else if (mode == MR74) 6554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 6564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* mean = 30 dB */ 6574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_temp2 = (Word32) 32588 << 6; 6584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_add(L_tmp, L_temp2, pOverflow); /* Q14 */ 6594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 6604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else if (mode == MR67) 6614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 6624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* mean = 28.75 dB */ 6634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_temp2 = (Word32) 32268 << 6; 6644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_add(L_tmp, L_temp2, pOverflow); /* Q14 */ 6654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 6664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else /* MR59, MR515, MR475 */ 6674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 6684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* mean = 33 dB */ 6694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_temp2 = (Word32) 16678 << 7; 6704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_add(L_tmp, L_temp2, pOverflow); /* Q14 */ 6714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 6724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 6734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /*-------------------------------------------------------------* 6744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * Compute gcode0. * 6754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * = Sum(i=0,3) pred[i]*past_qua_en[i] - ener_code + mean_ener * 6764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *--------------------------------------------------------------*/ 6774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* Q24 */ 6784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (L_tmp > (Word32) 0X001fffffL) 6794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 6804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *pOverflow = 1; 6814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = MAX_32; 6824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 6834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else if (L_tmp < (Word32) 0xffe00000L) 6844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 6854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *pOverflow = 1; 6864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = MIN_32; 6874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 6884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else 6894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 6904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_tmp << 10; 6914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 6924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 6934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 0; i < 4; i++) 6944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 6954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_temp2 = ((((Word32) pred[i]) * st->past_qua_en[i]) << 1); 6964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_add(L_tmp, L_temp2, pOverflow); /* Q13 * Q10 -> Q24 */ 6974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 6984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 6994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber gcode0 = (Word16)(L_tmp >> 16); /* Q8 */ 7004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 7014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /*-----------------------------------------------------------* 7024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * gcode0 = pow(10.0, gcode0/20) * 7034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * = pow(2, 3.3219*gcode0/20) * 7044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber * = pow(2, 0.166*gcode0) * 7054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *-----------------------------------------------------------*/ 7064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 7074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* 5439 Q15 = 0.165985 */ 7084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* (correct: 1/(20*log10(2)) 0.166096 = 5443 Q15) */ 7094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 7104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (mode == MR74) /* For IS641 bitexactness */ 7114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 7124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = (((Word32) gcode0) * 5439) << 1; /* Q8 * Q15 -> Q24 */ 7134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 7144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else 7154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 7164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = (((Word32) gcode0) * 5443) << 1; /* Q8 * Q15 -> Q24 */ 7174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 7184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 7194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (L_tmp < 0) 7204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 7214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = ~((~L_tmp) >> 8); 7224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 7234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else 7244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 7254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_tmp = L_tmp >> 8; /* -> Q16 */ 7264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 7274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 7284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *exp_gcode0 = (Word16)(L_tmp >> 16); 7294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (L_tmp < 0) 7304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 7314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_temp1 = ~((~L_tmp) >> 1); 7324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 7334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else 7344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 7354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_temp1 = L_tmp >> 1; 7364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 7374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber L_temp2 = (Word32) * exp_gcode0 << 15; 7384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *frac_gcode0 = (Word16)(L_sub(L_temp1, L_temp2, pOverflow)); 7394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* -> Q0.Q15 */ 7404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 7414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 7424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber return; 7434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 7444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 7454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/****************************************************************************/ 7464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 7474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* 7484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 7494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber FUNCTION NAME: gc_pred_update 7504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 7514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber INPUT AND OUTPUT DEFINITIONS 7524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 7534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Inputs: 7544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber st = pointer to a structure of type gc_predState 7554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber qua_ener_MR122 = quantized energy for update (Q10); calculated as 7564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber (log2(qua_err)) (Word16) 7574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber qua_ener = quantized energy for update (Q10); calculated as 7584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber (20*log10(qua_err)) (Word16) 7594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 7604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Outputs: 7614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber structure pointed to by st contains the calculated quantized energy 7624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for update 7634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 7644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Returns: 7654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 7664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 7674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Global Variables Used: 7684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 7694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 7704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Local Variables Needed: 7714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 7724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 7734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 7744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber FUNCTION DESCRIPTION 7754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 7764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber This function updates the MA predictor with the last quantized energy. 7774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 7784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 7794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REQUIREMENTS 7804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 7814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 7824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 7834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 7844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REFERENCES 7854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 7864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber gc_pred.c, UMTS GSM AMR speech codec, R99 - Version 3.2.0, March 2, 2001 7874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 7884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 7894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber PSEUDO-CODE 7904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 7914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Hubervoid gc_pred_update( 7924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber gc_predState *st, // i/o: State struct 7934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 qua_ener_MR122, // i : quantized energy for update, Q10 7944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // (log2(qua_err)) 7954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 qua_ener // i : quantized energy for update, Q10 7964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // (20*log10(qua_err)) 7974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber) 7984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 7994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 i; 8004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 3; i > 0; i--) 8024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 8034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber st->past_qua_en[i] = st->past_qua_en[i - 1]; 8044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber st->past_qua_en_MR122[i] = st->past_qua_en_MR122[i - 1]; 8054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 8064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber st->past_qua_en_MR122[0] = qua_ener_MR122; // log2 (qua_err), Q10 8084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber st->past_qua_en[0] = qua_ener; // 20*log10(qua_err), Q10 8104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 8124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 8144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber RESOURCES USED [optional] 8154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber When the code is written for a specific target processor the 8174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber the resources used should be documented below. 8184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber HEAP MEMORY USED: x bytes 8204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber STACK MEMORY USED: x bytes 8224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CLOCK CYCLES: (cycle count equation for this function) + (variable 8244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber used to represent cycle count for each subroutine 8254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber called) 8264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber where: (cycle count variable) = cycle count for [subroutine 8274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber name] 8284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 8304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CAUTION [optional] 8314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber [State any special notes, constraints or cautions for users of this function] 8324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 8344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber*/ 8354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Hubervoid gc_pred_update( 8374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber gc_predState *st, /* i/o: State struct */ 8384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 qua_ener_MR122, /* i : quantized energy for update, Q10 */ 8394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* (log2(qua_err)) */ 8404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 qua_ener /* i : quantized energy for update, Q10 */ 8414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* (20*log10(qua_err)) */ 8424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber) 8434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 8444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber st->past_qua_en[3] = st->past_qua_en[2]; 8454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber st->past_qua_en_MR122[3] = st->past_qua_en_MR122[2]; 8464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber st->past_qua_en[2] = st->past_qua_en[1]; 8484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber st->past_qua_en_MR122[2] = st->past_qua_en_MR122[1]; 8494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber st->past_qua_en[1] = st->past_qua_en[0]; 8514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber st->past_qua_en_MR122[1] = st->past_qua_en_MR122[0]; 8524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber st->past_qua_en_MR122[0] = qua_ener_MR122; /* log2 (qua_err), Q10 */ 8544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber st->past_qua_en[0] = qua_ener; /* 20*log10(qua_err), Q10 */ 8564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber return; 8584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 8594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/****************************************************************************/ 8614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber/* 8634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 8644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber FUNCTION NAME: gc_pred_average_limited 8654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 8664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber INPUT AND OUTPUT DEFINITIONS 8674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Inputs: 8694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber st = pointer to a structure of type gc_predState 8704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ener_avg_MR122 = pointer to the averaged quantized energy (Q10); 8714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber calculated as (log2(qua_err)) (Word16) 8724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber ener_avg = pointer to the averaged quantized energy (Q10); calculated 8734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber as (20*log10(qua_err)) (Word16) 8744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber pOverflow = pointer to overflow (Flag) 8754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Outputs: 8774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber store pointed to by ener_avg_MR122 contains the new averaged quantized 8784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber energy 8794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber store pointed to by ener_avg contains the new averaged quantized 8804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber energy 8814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber pOverflow = 1 if the math functions called by gc_pred_average_limited 8824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber results in overflow else zero. 8834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Returns: 8854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 8864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Global Variables Used: 8884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 8894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Local Variables Needed: 8914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 8924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 8944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber FUNCTION DESCRIPTION 8954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber This function calculates the average of MA predictor state values (with a 8974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber lower limit) used in error concealment. 8984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 8994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 9004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REQUIREMENTS 9014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 9024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber None 9034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 9044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 9054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber REFERENCES 9064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 9074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber gc_pred.c, UMTS GSM AMR speech codec, R99 - Version 3.2.0, March 2, 2001 9084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 9094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 9104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber PSEUDO-CODE 9114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 9124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas HuberThe original etsi reference code uses a global flag Overflow. However, in the 9134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huberactual implementation a pointer to a the overflow flag is passed in. 9144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 9154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Hubervoid gc_pred_average_limited( 9164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber gc_predState *st, // i: State struct 9174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *ener_avg_MR122, // o: everaged quantized energy, Q10 9184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // (log2(qua_err)) 9194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *ener_avg // o: averaged quantized energy, Q10 9204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // (20*log10(qua_err)) 9214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber) 9224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 9234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 av_pred_en; 9244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 i; 9254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 9264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // do average in MR122 mode (log2() domain) 9274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber av_pred_en = 0; 9284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 0; i < NPRED; i++) 9294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 9304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber av_pred_en = add (av_pred_en, st->past_qua_en_MR122[i]); 9314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 9324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 9334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // av_pred_en = 0.25*av_pred_en 9344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber av_pred_en = mult (av_pred_en, 8192); 9354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 9364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // if (av_pred_en < -14/(20Log10(2))) av_pred_en = .. 9374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 9384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (sub (av_pred_en, MIN_ENERGY_MR122) < 0) 9394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 9404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber av_pred_en = MIN_ENERGY_MR122; 9414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 9424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *ener_avg_MR122 = av_pred_en; 9434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 9444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // do average for other modes (20*log10() domain) 9454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber av_pred_en = 0; 9464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 0; i < NPRED; i++) 9474f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 9484f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber av_pred_en = add (av_pred_en, st->past_qua_en[i]); 9494f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 9504f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 9514f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // av_pred_en = 0.25*av_pred_en 9524f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber av_pred_en = mult (av_pred_en, 8192); 9534f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 9544f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber // if (av_pred_en < -14) av_pred_en = .. 9554f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 9564f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (sub (av_pred_en, MIN_ENERGY) < 0) 9574f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 9584f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber av_pred_en = MIN_ENERGY; 9594f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 9604f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *ener_avg = av_pred_en; 9614f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 9624f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 9634f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 9644f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber RESOURCES USED [optional] 9654f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 9664f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber When the code is written for a specific target processor the 9674f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber the resources used should be documented below. 9684f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 9694f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber HEAP MEMORY USED: x bytes 9704f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 9714f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber STACK MEMORY USED: x bytes 9724f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 9734f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CLOCK CYCLES: (cycle count equation for this function) + (variable 9744f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber used to represent cycle count for each subroutine 9754f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber called) 9764f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber where: (cycle count variable) = cycle count for [subroutine 9774f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber name] 9784f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 9794f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 9804f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber CAUTION [optional] 9814f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber [State any special notes, constraints or cautions for users of this function] 9824f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 9834f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber------------------------------------------------------------------------------ 9844f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber*/ 9854f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 9864f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Hubervoid gc_pred_average_limited( 9874f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber gc_predState *st, /* i: State struct */ 9884f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *ener_avg_MR122, /* o: everaged quantized energy, Q10 */ 9894f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* (log2(qua_err)) */ 9904f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 *ener_avg, /* o: averaged quantized energy, Q10 */ 9914f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* (20*log10(qua_err)) */ 9924f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Flag *pOverflow 9934f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber) 9944f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber{ 9954f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber Word16 av_pred_en; 9964f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber register Word16 i; 9974f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 9984f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* do average in MR122 mode (log2() domain) */ 9994f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber av_pred_en = 0; 10004f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 0; i < NPRED; i++) 10014f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 10024f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber av_pred_en = 10034f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber add(av_pred_en, st->past_qua_en_MR122[i], pOverflow); 10044f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 10054f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 10064f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* av_pred_en = 0.25*av_pred_en (with sign-extension)*/ 10074f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (av_pred_en < 0) 10084f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 10094f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber av_pred_en = (av_pred_en >> 2) | 0xc000; 10104f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 10114f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else 10124f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 10134f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber av_pred_en >>= 2; 10144f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 10154f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 10164f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* if (av_pred_en < -14/(20Log10(2))) av_pred_en = .. */ 10174f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (av_pred_en < MIN_ENERGY_MR122) 10184f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 10194f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber av_pred_en = MIN_ENERGY_MR122; 10204f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 10214f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *ener_avg_MR122 = av_pred_en; 10224f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 10234f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* do average for other modes (20*log10() domain) */ 10244f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber av_pred_en = 0; 10254f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber for (i = 0; i < NPRED; i++) 10264f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 10274f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber av_pred_en = add(av_pred_en, st->past_qua_en[i], pOverflow); 10284f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 10294f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 10304f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* av_pred_en = 0.25*av_pred_en (with sign-extension)*/ 10314f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (av_pred_en < 0) 10324f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 10334f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber av_pred_en = (av_pred_en >> 2) | 0xc000; 10344f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 10354f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber else 10364f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 10374f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber av_pred_en >>= 2; 10384f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 10394f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber 10404f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber /* if (av_pred_en < -14) av_pred_en = .. */ 10414f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber if (av_pred_en < MIN_ENERGY) 10424f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber { 10434f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber av_pred_en = MIN_ENERGY; 10444f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber } 10454f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber *ener_avg = av_pred_en; 10464f1efc098cb5791c3e9f483f2af84aef70d2d0a0Andreas Huber} 1047