1b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/* ------------------------------------------------------------------ 2b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * Copyright (C) 1998-2009 PacketVideo 3b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * 4b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * Licensed under the Apache License, Version 2.0 (the "License"); 5b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * you may not use this file except in compliance with the License. 6b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * You may obtain a copy of the License at 7b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * 8b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * http://www.apache.org/licenses/LICENSE-2.0 9b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * 10b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * Unless required by applicable law or agreed to in writing, software 11b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * distributed under the License is distributed on an "AS IS" BASIS, 12b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either 13b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * express or implied. 14b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * See the License for the specific language governing permissions 15b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * and limitations under the License. 16b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * ------------------------------------------------------------------- 17b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber */ 18b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/**************************************************************************************** 19b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas HuberPortions of this file are derived from the following 3GPP standard: 20b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 21b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 3GPP TS 26.073 22b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber ANSI-C code for the Adaptive Multi-Rate (AMR) speech codec 23b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Available from http://www.3gpp.org 24b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 25b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber(C) 2004, 3GPP Organizational Partners (ARIB, ATIS, CCSA, ETSI, TTA, TTC) 26b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas HuberPermission to distribute, modify and use this file under the standard license 27b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huberterms listed above has been obtained from the copyright holder. 28b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber****************************************************************************************/ 29b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/* 30b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 31b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 32b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 33b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 34b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Pathname: ./audio/gsm-amr/c/src/qgain475.c 35b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Funtions: MR475_quant_store_results 36b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber MR475_update_unq_pred 37b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber MR475_gain_quant 38b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 39b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 40b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber MODULE DESCRIPTION 41b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 42b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber These modules handle the quantization of pitch and codebook gains for MR475. 43b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 44b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 45b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber*/ 46b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 47b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 48b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*---------------------------------------------------------------------------- 49b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; INCLUDES 50b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/ 51b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "qgain475.h" 52b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "typedef.h" 53b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "basic_op.h" 54b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "mode.h" 55b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "cnst.h" 56b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "pow2.h" 57b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#include "log2.h" 58b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 59b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*---------------------------------------------------------------------------- 60b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; MACROS 61b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; Define module specific macros here 62b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/ 63b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 64b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 65b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*---------------------------------------------------------------------------- 66b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; DEFINES 67b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; Include all pre-processor statements here. Include conditional 68b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; compile variables also. 69b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/ 70b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber#define MR475_VQ_SIZE 256 71b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 72b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*---------------------------------------------------------------------------- 73b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; LOCAL FUNCTION DEFINITIONS 74b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; Function Prototype declaration 75b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/ 76b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 77b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/*---------------------------------------------------------------------------- 78b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; LOCAL VARIABLE DEFINITIONS 79b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber; Variable declaration - defined here and used outside this module 80b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber----------------------------------------------------------------------------*/ 81b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 82b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/* The table contains the following data: 83b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * 84b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * g_pitch(0) (Q14) // for sub- 85b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * g_fac(0) (Q12) // frame 0 and 2 86b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * g_pitch(1) (Q14) // for sub- 87b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * g_fac(2) (Q12) // frame 1 and 3 88b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * 89b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber */ 90b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huberstatic const Word16 table_gain_MR475[MR475_VQ_SIZE*4] = 91b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber{ 92b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /*g_pit(0), g_fac(0), g_pit(1), g_fac(1) */ 93b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 812, 128, 542, 140, 94b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 2873, 1135, 2266, 3402, 95b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 2067, 563, 12677, 647, 96b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 4132, 1798, 5601, 5285, 97b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 7689, 374, 3735, 441, 98b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 10912, 2638, 11807, 2494, 99b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 20490, 797, 5218, 675, 100b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 6724, 8354, 5282, 1696, 101b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1488, 428, 5882, 452, 102b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 5332, 4072, 3583, 1268, 103b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 2469, 901, 15894, 1005, 104b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 14982, 3271, 10331, 4858, 105b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 3635, 2021, 2596, 835, 106b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 12360, 4892, 12206, 1704, 107b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 13432, 1604, 9118, 2341, 108b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 3968, 1538, 5479, 9936, 109b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 3795, 417, 1359, 414, 110b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 3640, 1569, 7995, 3541, 111b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 11405, 645, 8552, 635, 112b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 4056, 1377, 16608, 6124, 113b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 11420, 700, 2007, 607, 114b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 12415, 1578, 11119, 4654, 115b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 13680, 1708, 11990, 1229, 116b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 7996, 7297, 13231, 5715, 117b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 2428, 1159, 2073, 1941, 118b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 6218, 6121, 3546, 1804, 119b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 8925, 1802, 8679, 1580, 120b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 13935, 3576, 13313, 6237, 121b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 6142, 1130, 5994, 1734, 122b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 14141, 4662, 11271, 3321, 123b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 12226, 1551, 13931, 3015, 124b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 5081, 10464, 9444, 6706, 125b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1689, 683, 1436, 1306, 126b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 7212, 3933, 4082, 2713, 127b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 7793, 704, 15070, 802, 128b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 6299, 5212, 4337, 5357, 129b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 6676, 541, 6062, 626, 130b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 13651, 3700, 11498, 2408, 131b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 16156, 716, 12177, 751, 132b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 8065, 11489, 6314, 2256, 133b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 4466, 496, 7293, 523, 134b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 10213, 3833, 8394, 3037, 135b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 8403, 966, 14228, 1880, 136b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 8703, 5409, 16395, 4863, 137b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 7420, 1979, 6089, 1230, 138b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 9371, 4398, 14558, 3363, 139b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 13559, 2873, 13163, 1465, 140b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 5534, 1678, 13138, 14771, 141b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 7338, 600, 1318, 548, 142b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 4252, 3539, 10044, 2364, 143b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 10587, 622, 13088, 669, 144b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 14126, 3526, 5039, 9784, 145b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 15338, 619, 3115, 590, 146b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 16442, 3013, 15542, 4168, 147b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 15537, 1611, 15405, 1228, 148b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 16023, 9299, 7534, 4976, 149b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1990, 1213, 11447, 1157, 150b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 12512, 5519, 9475, 2644, 151b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 7716, 2034, 13280, 2239, 152b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 16011, 5093, 8066, 6761, 153b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 10083, 1413, 5002, 2347, 154b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 12523, 5975, 15126, 2899, 155b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 18264, 2289, 15827, 2527, 156b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 16265, 10254, 14651, 11319, 157b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1797, 337, 3115, 397, 158b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 3510, 2928, 4592, 2670, 159b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 7519, 628, 11415, 656, 160b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 5946, 2435, 6544, 7367, 161b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 8238, 829, 4000, 863, 162b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 10032, 2492, 16057, 3551, 163b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 18204, 1054, 6103, 1454, 164b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 5884, 7900, 18752, 3468, 165b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1864, 544, 9198, 683, 166b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 11623, 4160, 4594, 1644, 167b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 3158, 1157, 15953, 2560, 168b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 12349, 3733, 17420, 5260, 169b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 6106, 2004, 2917, 1742, 170b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 16467, 5257, 16787, 1680, 171b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 17205, 1759, 4773, 3231, 172b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 7386, 6035, 14342, 10012, 173b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 4035, 442, 4194, 458, 174b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 9214, 2242, 7427, 4217, 175b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 12860, 801, 11186, 825, 176b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 12648, 2084, 12956, 6554, 177b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 9505, 996, 6629, 985, 178b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 10537, 2502, 15289, 5006, 179b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 12602, 2055, 15484, 1653, 180b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 16194, 6921, 14231, 5790, 181b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 2626, 828, 5615, 1686, 182b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 13663, 5778, 3668, 1554, 183b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 11313, 2633, 9770, 1459, 184b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 14003, 4733, 15897, 6291, 185b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 6278, 1870, 7910, 2285, 186b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 16978, 4571, 16576, 3849, 187b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 15248, 2311, 16023, 3244, 188b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 14459, 17808, 11847, 2763, 189b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1981, 1407, 1400, 876, 190b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 4335, 3547, 4391, 4210, 191b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 5405, 680, 17461, 781, 192b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 6501, 5118, 8091, 7677, 193b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 7355, 794, 8333, 1182, 194b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 15041, 3160, 14928, 3039, 195b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 20421, 880, 14545, 852, 196b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 12337, 14708, 6904, 1920, 197b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 4225, 933, 8218, 1087, 198b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 10659, 4084, 10082, 4533, 199b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 2735, 840, 20657, 1081, 200b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 16711, 5966, 15873, 4578, 201b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 10871, 2574, 3773, 1166, 202b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 14519, 4044, 20699, 2627, 203b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 15219, 2734, 15274, 2186, 204b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 6257, 3226, 13125, 19480, 205b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 7196, 930, 2462, 1618, 206b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 4515, 3092, 13852, 4277, 207b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 10460, 833, 17339, 810, 208b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 16891, 2289, 15546, 8217, 209b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 13603, 1684, 3197, 1834, 210b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 15948, 2820, 15812, 5327, 211b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 17006, 2438, 16788, 1326, 212b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 15671, 8156, 11726, 8556, 213b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 3762, 2053, 9563, 1317, 214b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 13561, 6790, 12227, 1936, 215b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 8180, 3550, 13287, 1778, 216b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 16299, 6599, 16291, 7758, 217b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 8521, 2551, 7225, 2645, 218b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 18269, 7489, 16885, 2248, 219b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 17882, 2884, 17265, 3328, 220b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 9417, 20162, 11042, 8320, 221b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1286, 620, 1431, 583, 222b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 5993, 2289, 3978, 3626, 223b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 5144, 752, 13409, 830, 224b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 5553, 2860, 11764, 5908, 225b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 10737, 560, 5446, 564, 226b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 13321, 3008, 11946, 3683, 227b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 19887, 798, 9825, 728, 228b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 13663, 8748, 7391, 3053, 229b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 2515, 778, 6050, 833, 230b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 6469, 5074, 8305, 2463, 231b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 6141, 1865, 15308, 1262, 232b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 14408, 4547, 13663, 4515, 233b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 3137, 2983, 2479, 1259, 234b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 15088, 4647, 15382, 2607, 235b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 14492, 2392, 12462, 2537, 236b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 7539, 2949, 12909, 12060, 237b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 5468, 684, 3141, 722, 238b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 5081, 1274, 12732, 4200, 239b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 15302, 681, 7819, 592, 240b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 6534, 2021, 16478, 8737, 241b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 13364, 882, 5397, 899, 242b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 14656, 2178, 14741, 4227, 243b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 14270, 1298, 13929, 2029, 244b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 15477, 7482, 15815, 4572, 245b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 2521, 2013, 5062, 1804, 246b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 5159, 6582, 7130, 3597, 247b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 10920, 1611, 11729, 1708, 248b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 16903, 3455, 16268, 6640, 249b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 9306, 1007, 9369, 2106, 250b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 19182, 5037, 12441, 4269, 251b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 15919, 1332, 15357, 3512, 252b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 11898, 14141, 16101, 6854, 253b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 2010, 737, 3779, 861, 254b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 11454, 2880, 3564, 3540, 255b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 9057, 1241, 12391, 896, 256b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 8546, 4629, 11561, 5776, 257b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 8129, 589, 8218, 588, 258b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 18728, 3755, 12973, 3149, 259b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 15729, 758, 16634, 754, 260b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 15222, 11138, 15871, 2208, 261b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 4673, 610, 10218, 678, 262b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 15257, 4146, 5729, 3327, 263b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 8377, 1670, 19862, 2321, 264b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 15450, 5511, 14054, 5481, 265b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 5728, 2888, 7580, 1346, 266b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 14384, 5325, 16236, 3950, 267b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 15118, 3744, 15306, 1435, 268b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 14597, 4070, 12301, 15696, 269b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 7617, 1699, 2170, 884, 270b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 4459, 4567, 18094, 3306, 271b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 12742, 815, 14926, 907, 272b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 15016, 4281, 15518, 8368, 273b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 17994, 1087, 2358, 865, 274b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 16281, 3787, 15679, 4596, 275b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 16356, 1534, 16584, 2210, 276b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 16833, 9697, 15929, 4513, 277b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 3277, 1085, 9643, 2187, 278b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 11973, 6068, 9199, 4462, 279b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 8955, 1629, 10289, 3062, 280b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 16481, 5155, 15466, 7066, 281b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 13678, 2543, 5273, 2277, 282b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 16746, 6213, 16655, 3408, 283b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 20304, 3363, 18688, 1985, 284b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 14172, 12867, 15154, 15703, 285b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 4473, 1020, 1681, 886, 286b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 4311, 4301, 8952, 3657, 287b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 5893, 1147, 11647, 1452, 288b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 15886, 2227, 4582, 6644, 289b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 6929, 1205, 6220, 799, 290b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 12415, 3409, 15968, 3877, 291b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 19859, 2109, 9689, 2141, 292b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 14742, 8830, 14480, 2599, 293b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1817, 1238, 7771, 813, 294b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 19079, 4410, 5554, 2064, 295b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 3687, 2844, 17435, 2256, 296b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 16697, 4486, 16199, 5388, 297b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 8028, 2763, 3405, 2119, 298b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 17426, 5477, 13698, 2786, 299b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 19879, 2720, 9098, 3880, 300b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 18172, 4833, 17336, 12207, 301b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 5116, 996, 4935, 988, 302b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 9888, 3081, 6014, 5371, 303b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 15881, 1667, 8405, 1183, 304b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 15087, 2366, 19777, 7002, 305b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 11963, 1562, 7279, 1128, 306b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 16859, 1532, 15762, 5381, 307b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 14708, 2065, 20105, 2155, 308b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 17158, 8245, 17911, 6318, 309b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 5467, 1504, 4100, 2574, 310b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 17421, 6810, 5673, 2888, 311b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 16636, 3382, 8975, 1831, 312b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 20159, 4737, 19550, 7294, 313b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 6658, 2781, 11472, 3321, 314b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 19397, 5054, 18878, 4722, 315b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 16439, 2373, 20430, 4386, 316b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 11353, 26526, 11593, 3068, 317b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 2866, 1566, 5108, 1070, 318b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 9614, 4915, 4939, 3536, 319b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 7541, 878, 20717, 851, 320b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 6938, 4395, 16799, 7733, 321b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 10137, 1019, 9845, 964, 322b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 15494, 3955, 15459, 3430, 323b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 18863, 982, 20120, 963, 324b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 16876, 12887, 14334, 4200, 325b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 6599, 1220, 9222, 814, 326b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 16942, 5134, 5661, 4898, 327b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 5488, 1798, 20258, 3962, 328b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 17005, 6178, 17929, 5929, 329b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 9365, 3420, 7474, 1971, 330b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 19537, 5177, 19003, 3006, 331b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 16454, 3788, 16070, 2367, 332b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 8664, 2743, 9445, 26358, 333b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 10856, 1287, 3555, 1009, 334b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 5606, 3622, 19453, 5512, 335b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 12453, 797, 20634, 911, 336b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 15427, 3066, 17037, 10275, 337b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 18883, 2633, 3913, 1268, 338b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 19519, 3371, 18052, 5230, 339b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 19291, 1678, 19508, 3172, 340b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 18072, 10754, 16625, 6845, 341b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 3134, 2298, 10869, 2437, 342b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 15580, 6913, 12597, 3381, 343b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 11116, 3297, 16762, 2424, 344b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 18853, 6715, 17171, 9887, 345b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 12743, 2605, 8937, 3140, 346b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 19033, 7764, 18347, 3880, 347b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 20475, 3682, 19602, 3380, 348b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 13044, 19373, 10526, 23124 349b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber}; 350b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 351b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/* 352b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 353b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber FUNCTION NAME: MR475_quant_store_results 354b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 355b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber INPUT AND OUTPUT DEFINITIONS 356b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 357b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Inputs: 358b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber pred_st = pointer to structure of type gc_predState 359b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber p = pointer to selected quantizer table entry (const Word16) 360b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber gcode0 = predicted CB gain (Word16) 361b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp_gcode0 = exponent of predicted CB gain (Word16) 362b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber gain_pit = pointer to Pitch gain (Word16) 363b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber gain_cod = pointer to Code gain (Word16) 364b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 365b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Outputs: 366b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber pred_st points to the updated structure of type gc_predState 367b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber gain_pit points to Pitch gain 368b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber gain_cod points to Code gain 369b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber pOverflow points to overflow indicator (Flag) 370b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 371b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Returns: 372b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber None. 373b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 374b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Global Variables Used: 375b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber None. 376b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 377b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Local Variables Needed: 378b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber None. 379b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 380b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 381b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber FUNCTION DESCRIPTION 382b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 383b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber This function calculates the final fixed codebook gain and the predictor 384b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber update values, and updates the gain predictor. 385b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 386b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 387b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber REQUIREMENTS 388b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 389b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber None. 390b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 391b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 392b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber REFERENCES 393b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 394b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber qgain475.c, UMTS GSM AMR speech codec, R99 - Version 3.2.0, March 2, 2001 395b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 396b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 397b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber PSEUDO-CODE 398b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 399b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huberstatic void MR475_quant_store_results( 400b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 401b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber gc_predState *pred_st, // i/o: gain predictor state struct 402b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber const Word16 *p, // i : pointer to selected quantizer table entry 403b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 gcode0, // i : predicted CB gain, Q(14 - exp_gcode0) 404b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 exp_gcode0, // i : exponent of predicted CB gain, Q0 405b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *gain_pit, // o : Pitch gain, Q14 406b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *gain_cod // o : Code gain, Q1 407b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber) 408b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber{ 409b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 410b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 g_code, exp, frac, tmp; 411b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word32 L_tmp; 412b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 413b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 qua_ener_MR122; // o : quantized energy error, MR122 version Q10 414b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 qua_ener; // o : quantized energy error, Q10 415b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 416b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // Read the quantized gains 417b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *gain_pit = *p++; 418b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g_code = *p++; 419b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 420b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber //------------------------------------------------------------------* 421b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * calculate final fixed codebook gain: * 422b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ * 423b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * * 424b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * gc = gc0 * g * 425b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *------------------------------------------------------------------ 426b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 427b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = L_mult(g_code, gcode0); 428b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = L_shr(L_tmp, sub(10, exp_gcode0)); 429b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *gain_cod = extract_h(L_tmp); 430b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 431b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber //------------------------------------------------------------------* 432b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * calculate predictor update values and update gain predictor: * 433b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ * 434b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * * 435b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * qua_ener = log2(g) * 436b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * qua_ener_MR122 = 20*log10(g) * 437b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *------------------------------------------------------------------ 438b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 439b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Log2 (L_deposit_l (g_code), &exp, &frac); // Log2(x Q12) = log2(x) + 12 440b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp = sub(exp, 12); 441b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 442b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp = shr_r (frac, 5); 443b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber qua_ener_MR122 = add (tmp, shl (exp, 10)); 444b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 445b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = Mpy_32_16(exp, frac, 24660); // 24660 Q12 ~= 6.0206 = 20*log10(2) 446b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber qua_ener = pv_round (L_shl (L_tmp, 13)); // Q12 * Q0 = Q13 -> Q10 447b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 448b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber gc_pred_update(pred_st, qua_ener_MR122, qua_ener); 449b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber} 450b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 451b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 452b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber RESOURCES USED [optional] 453b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 454b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber When the code is written for a specific target processor the 455b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber the resources used should be documented below. 456b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 457b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber HEAP MEMORY USED: x bytes 458b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 459b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber STACK MEMORY USED: x bytes 460b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 461b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber CLOCK CYCLES: (cycle count equation for this function) + (variable 462b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber used to represent cycle count for each subroutine 463b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber called) 464b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber where: (cycle count variable) = cycle count for [subroutine 465b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber name] 466b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 467b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 468b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber CAUTION [optional] 469b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber [State any special notes, constraints or cautions for users of this function] 470b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 471b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 472b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber*/ 473b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 474b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huberstatic void MR475_quant_store_results( 475b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber gc_predState *pred_st, /* i/o: gain predictor state struct */ 476b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber const Word16 *p, /* i : pointer to selected quantizer table entry */ 477b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 gcode0, /* i : predicted CB gain, Q(14 - exp_gcode0) */ 478b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 exp_gcode0, /* i : exponent of predicted CB gain, Q0 */ 479b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *gain_pit, /* o : Pitch gain, Q14 */ 480b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *gain_cod, /* o : Code gain, Q1 */ 481b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Flag *pOverflow /* o : overflow indicator */ 482b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber) 483b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber{ 484b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 g_code; 485b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 exp; 486b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 frac; 487b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 tmp; 488b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word32 L_tmp; 489b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 490b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 qua_ener_MR122; /* o : quantized energy error, MR122 version Q10 */ 491b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 qua_ener; /* o : quantized energy error, Q10 */ 492b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 493b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 494b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* Read the quantized gains */ 495b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *gain_pit = *p++; 496b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g_code = *p++; 497b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 498b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /*------------------------------------------------------------------* 499b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * calculate final fixed codebook gain: * 500b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ * 501b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * * 502b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * gc = gc0 * g * 503b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *------------------------------------------------------------------*/ 504b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 505b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = ((Word32) g_code * gcode0) << 1; 506b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp = 10 - exp_gcode0; 507b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = L_shr(L_tmp, tmp, pOverflow); 508b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *gain_cod = (Word16)(L_tmp >> 16); 509b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 510b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /*------------------------------------------------------------------* 511b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * calculate predictor update values and update gain predictor: * 512b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ * 513b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * * 514b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * qua_ener = log2(g) * 515b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * qua_ener_MR122 = 20*log10(g) * 516b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *------------------------------------------------------------------*/ 517b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 518b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* Log2(x Q12) = log2(x) + 12 */ 519b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Log2((Word32) g_code, &exp, &frac, pOverflow); 520b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp -= 12; 521b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 522b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp = shr_r(frac, 5, pOverflow); 523b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber qua_ener_MR122 = exp << 10; 524b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber qua_ener_MR122 = tmp + qua_ener_MR122; 525b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 526b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* 24660 Q12 ~= 6.0206 = 20*log10(2) */ 527b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = Mpy_32_16(exp, frac, 24660, pOverflow); 528b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = L_tmp << 13; 529b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 530b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* Q12 * Q0 = Q13 -> Q10 */ 531b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber qua_ener = (Word16)((L_tmp + (Word32) 0x00008000L) >> 16); 532b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 533b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber gc_pred_update(pred_st, qua_ener_MR122, qua_ener); 534b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 535b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber return; 536b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber} 537b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 538b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/****************************************************************************/ 539b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 540b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 541b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/* 542b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 543b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber FUNCTION NAME: MR475_update_unq_pred 544b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 545b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber INPUT AND OUTPUT DEFINITIONS 546b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 547b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Inputs: 548b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber pred_st = pointer to structure of type gc_predState 549b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp_gcode0 = predicted CB gain (exponent MSW) (Word16) 550b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber frac_gcode0 = predicted CB gain (exponent LSW) (Word16) 551b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber cod_gain_exp = optimum codebook gain (exponent)(Word16) 552b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber cod_gain_frac = optimum codebook gain (fraction) (Word16) 553b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 554b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Outputs: 555b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber pred_st points to the updated structure of type gc_predState 556b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber pOverflow points to overflow indicator (Flag) 557b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 558b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Returns: 559b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber None. 560b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 561b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Global Variables Used: 562b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber None. 563b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 564b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Local Variables Needed: 565b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber None. 566b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 567b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 568b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber FUNCTION DESCRIPTION 569b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 570b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber This module uses the optimum codebook gain and updates the "unquantized" 571b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber gain predictor with the (bounded) prediction error. 572b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 573b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 574b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber REQUIREMENTS 575b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 576b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber None. 577b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 578b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 579b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber REFERENCES 580b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 581b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber qgain475.c, UMTS GSM AMR speech codec, R99 - Version 3.2.0, March 2, 2001 582b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 583b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 584b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber PSEUDO-CODE 585b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 586b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Hubervoid 587b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas HuberMR475_update_unq_pred( 588b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber gc_predState *pred_st, // i/o: gain predictor state struct 589b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 exp_gcode0, // i : predicted CB gain (exponent MSW), Q0 590b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 frac_gcode0, // i : predicted CB gain (exponent LSW), Q15 591b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 cod_gain_exp, // i : optimum codebook gain (exponent), Q0 592b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 cod_gain_frac // i : optimum codebook gain (fraction), Q15 593b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber) 594b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber{ 595b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 tmp, exp, frac; 596b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 qua_ener, qua_ener_MR122; 597b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word32 L_tmp; 598b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 599b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // calculate prediction error factor (given optimum CB gain gcu): 600b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // predErrFact = gcu / gcode0 601b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // (limit to MIN_PRED_ERR_FACT <= predErrFact <= MAX_PRED_ERR_FACT 602b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // -> limit qua_ener*) 603b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // 604b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // calculate prediction error (log): 605b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // 606b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // qua_ener_MR122 = log2(predErrFact) 607b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // qua_ener = 20*log10(predErrFact) 608b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 609b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber if (cod_gain_frac <= 0) 610b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 611b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // if gcu <= 0 -> predErrFact = 0 < MIN_PRED_ERR_FACT 612b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // -> set qua_ener(_MR122) directly 613b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber qua_ener = MIN_QUA_ENER; 614b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber qua_ener_MR122 = MIN_QUA_ENER_MR122; 615b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 616b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber else 617b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 618b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // convert gcode0 from DPF to standard fraction/exponent format 619b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // with normalized frac, i.e. 16384 <= frac <= 32767 620b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // Note: exponent correction (exp=exp-14) is done after div_s 621b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber frac_gcode0 = extract_l (Pow2 (14, frac_gcode0)); 622b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 623b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // make sure cod_gain_frac < frac_gcode0 for div_s 624b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber if (sub(cod_gain_frac, frac_gcode0) >= 0) 625b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 626b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber cod_gain_frac = shr (cod_gain_frac, 1); 627b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber cod_gain_exp = add (cod_gain_exp, 1); 628b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 629b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 630b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // predErrFact 631b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // = gcu / gcode0 632b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // = cod_gain_frac/frac_gcode0 * 2^(cod_gain_exp-(exp_gcode0-14)) 633b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // = div_s (c_g_f, frac_gcode0)*2^-15 * 2^(c_g_e-exp_gcode0+14) 634b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // = div_s * 2^(cod_gain_exp-exp_gcode0 - 1) 635b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 636b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber frac = div_s (cod_gain_frac, frac_gcode0); 637b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp = sub (sub (cod_gain_exp, exp_gcode0), 1); 638b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 639b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Log2 (L_deposit_l (frac), &exp, &frac); 640b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp = add (exp, tmp); 641b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 642b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // calculate prediction error (log2, Q10) 643b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber qua_ener_MR122 = shr_r (frac, 5); 644b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber qua_ener_MR122 = add (qua_ener_MR122, shl (exp, 10)); 645b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 646b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber if (sub(qua_ener_MR122, MIN_QUA_ENER_MR122) < 0) 647b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 648b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber qua_ener = MIN_QUA_ENER; 649b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber qua_ener_MR122 = MIN_QUA_ENER_MR122; 650b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 651b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber else if (sub(qua_ener_MR122, MAX_QUA_ENER_MR122) > 0) 652b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 653b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber qua_ener = MAX_QUA_ENER; 654b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber qua_ener_MR122 = MAX_QUA_ENER_MR122; 655b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 656b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber else 657b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 658b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // calculate prediction error (20*log10, Q10) 659b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = Mpy_32_16(exp, frac, 24660); 660b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // 24660 Q12 ~= 6.0206 = 20*log10(2) 661b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber qua_ener = pv_round (L_shl (L_tmp, 13)); 662b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // Q12 * Q0 = Q13 -> Q26 -> Q10 663b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 664b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 665b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 666b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // update MA predictor memory 667b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber gc_pred_update(pred_st, qua_ener_MR122, qua_ener); 668b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber} 669b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 670b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 671b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber RESOURCES USED [optional] 672b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 673b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber When the code is written for a specific target processor the 674b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber the resources used should be documented below. 675b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 676b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber HEAP MEMORY USED: x bytes 677b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 678b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber STACK MEMORY USED: x bytes 679b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 680b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber CLOCK CYCLES: (cycle count equation for this function) + (variable 681b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber used to represent cycle count for each subroutine 682b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber called) 683b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber where: (cycle count variable) = cycle count for [subroutine 684b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber name] 685b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 686b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 687b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber CAUTION [optional] 688b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber [State any special notes, constraints or cautions for users of this function] 689b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 690b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 691b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber*/ 692b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 693b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Hubervoid MR475_update_unq_pred( 694b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber gc_predState *pred_st, /* i/o: gain predictor state struct */ 695b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 exp_gcode0, /* i : predicted CB gain (exponent MSW), Q0 */ 696b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 frac_gcode0, /* i : predicted CB gain (exponent LSW), Q15 */ 697b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 cod_gain_exp, /* i : optimum codebook gain (exponent), Q0 */ 698b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 cod_gain_frac, /* i : optimum codebook gain (fraction), Q15 */ 699b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Flag *pOverflow /* o : overflow indicator */ 700b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber) 701b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber{ 702b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 tmp; 703b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 exp; 704b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 frac; 705b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 qua_ener; 706b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 qua_ener_MR122; 707b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word32 L_tmp; 708b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 709b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* calculate prediction error factor (given optimum CB gain gcu): 710b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * 711b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * predErrFact = gcu / gcode0 712b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * (limit to MIN_PRED_ERR_FACT <= predErrFact <= MAX_PRED_ERR_FACT 713b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * -> limit qua_ener*) 714b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * 715b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * calculate prediction error (log): 716b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * 717b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * qua_ener_MR122 = log2(predErrFact) 718b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * qua_ener = 20*log10(predErrFact) 719b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * 720b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber */ 721b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 722b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber if (cod_gain_frac <= 0) 723b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 724b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* if gcu <= 0 -> predErrFact = 0 < MIN_PRED_ERR_FACT */ 725b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* -> set qua_ener(_MR122) directly */ 726b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber qua_ener = MIN_QUA_ENER; 727b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber qua_ener_MR122 = MIN_QUA_ENER_MR122; 728b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 729b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber else 730b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 731b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* convert gcode0 from DPF to standard fraction/exponent format */ 732b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* with normalized frac, i.e. 16384 <= frac <= 32767 */ 733b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* Note: exponent correction (exp=exp-14) is done after div_s */ 734b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber frac_gcode0 = (Word16)(Pow2(14, frac_gcode0, pOverflow)); 735b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 736b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* make sure cod_gain_frac < frac_gcode0 for div_s */ 737b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber if (cod_gain_frac >= frac_gcode0) 738b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 739b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber cod_gain_frac >>= 1; 740b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber cod_gain_exp += 1; 741b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 742b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 743b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* 744b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber predErrFact 745b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber = gcu / gcode0 746b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber = cod_gain_frac/frac_gcode0 * 2^(cod_gain_exp-(exp_gcode0-14)) 747b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber = div_s (c_g_f, frac_gcode0)*2^-15 * 2^(c_g_e-exp_gcode0+14) 748b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber = div_s * 2^(cod_gain_exp-exp_gcode0 - 1) 749b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber */ 750b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber frac = div_s(cod_gain_frac, frac_gcode0); 751b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp = cod_gain_exp - exp_gcode0; 752b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp -= 1; 753b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 754b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Log2((Word32) frac, &exp, &frac, pOverflow); 755b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp += tmp; 756b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 757b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* calculate prediction error (log2, Q10) */ 758b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber qua_ener_MR122 = shr_r(frac, 5, pOverflow); 759b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp = exp << 10; 760b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber qua_ener_MR122 += tmp; 761b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 762b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber if (qua_ener_MR122 > MAX_QUA_ENER_MR122) 763b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 764b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber qua_ener = MAX_QUA_ENER; 765b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber qua_ener_MR122 = MAX_QUA_ENER_MR122; 766b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 767b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber else 768b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 769b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* calculate prediction error (20*log10, Q10) */ 770b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = Mpy_32_16(exp, frac, 24660, pOverflow); 771b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* 24660 Q12 ~= 6.0206 = 20*log10(2) */ 772b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = L_shl(L_tmp, 13, pOverflow); 773b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber qua_ener = pv_round(L_tmp, pOverflow); 774b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 775b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* Q12 * Q0 = Q13 -> Q26 -> Q10 */ 776b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 777b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 778b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 779b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* update MA predictor memory */ 780b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber gc_pred_update(pred_st, qua_ener_MR122, qua_ener); 781b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 782b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 783b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber return; 784b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber} 785b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 786b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/****************************************************************************/ 787b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 788b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 789b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber/* 790b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 791b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber FUNCTION NAME: MR475_gain_quant 792b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 793b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber INPUT AND OUTPUT DEFINITIONS 794b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 795b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Inputs: 796b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber pred_st = pointer to structure of type gc_predState 797b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf0_exp_gcode0 = predicted CB gain (exponent) (Word16) 798b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber f0_frac_gcode0 = predicted CB gain (fraction) (Word16) 799b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf0_exp_coeff = energy coeff. (exponent part) (Word16) 800b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf0_frac_coeff = energy coeff. ((fraction part) (Word16) 801b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf0_exp_target_en = exponent of target energy (Word16) 802b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf0_frac_target_en = fraction of target energy (Word16) 803b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf1_code_nosharp = innovative codebook vector (Word16) 804b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf1_exp_gcode0 = predicted CB gain (exponent) (Word16) 805b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf1_frac_gcode0 = predicted CB gain (fraction) (Word16) 806b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf1_exp_coeff = energy coeff. (exponent part) (Word16) 807b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf1_frac_coeff = energy coeff. (fraction part) (Word16) 808b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf1_exp_target_en = exponent of target energy (Word16) 809b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf1_frac_target_en = fraction of target energy (Word16) 810b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber gp_limit = pitch gain limit (Word16) 811b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf0_gain_pit = pointer to Pitch gain (Word16) 812b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf0_gain_cod = pointer to Code gain (Word16) 813b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf1_gain_pit = pointer to Pitch gain (Word16) 814b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf1_gain_cod = pointer to Code gain (Word16) 815b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 816b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Outputs: 817b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber pred_st points to the updated structure of type gc_predState 818b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf0_gain_pit points to Pitch gain 819b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf0_gain_cod points to Code gain 820b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf1_gain_pit points to Pitch gain 821b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf1_gain_cod points to Code gain 822b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 823b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Returns: 824b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber index = index of quantization 825b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 826b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Global Variables Used: 827b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber None. 828b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 829b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Local Variables Needed: 830b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber None. 831b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 832b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 833b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber FUNCTION DESCRIPTION 834b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 835b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber This module provides quantization of pitch and codebook gains for two 836b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber subframes using the predicted codebook gain. 837b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 838b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 839b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber REQUIREMENTS 840b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 841b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber None. 842b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 843b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 844b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber REFERENCES 845b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 846b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber qgain475.c, UMTS GSM AMR speech codec, R99 - Version 3.2.0, March 2, 2001 847b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 848b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 849b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber PSEUDO-CODE 850b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 851b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas HuberWord16 852b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas HuberMR475_gain_quant( // o : index of quantization. 853b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber gc_predState *pred_st, // i/o: gain predictor state struct 854b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 855b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // data from subframe 0 (or 2) 856b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf0_exp_gcode0, // i : predicted CB gain (exponent), Q0 857b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf0_frac_gcode0, // i : predicted CB gain (fraction), Q15 858b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf0_exp_coeff[], // i : energy coeff. (5), exponent part, Q0 859b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf0_frac_coeff[], // i : energy coeff. (5), fraction part, Q15 860b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // (frac_coeff and exp_coeff computed in 861b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // calc_filt_energies()) 862b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf0_exp_target_en, // i : exponent of target energy, Q0 863b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf0_frac_target_en, // i : fraction of target energy, Q15 864b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 865b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // data from subframe 1 (or 3) 866b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf1_code_nosharp[], // i : innovative codebook vector (L_SUBFR) 867b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // (whithout pitch sharpening) 868b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf1_exp_gcode0, // i : predicted CB gain (exponent), Q0 869b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf1_frac_gcode0, // i : predicted CB gain (fraction), Q15 870b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf1_exp_coeff[], // i : energy coeff. (5), exponent part, Q0 871b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf1_frac_coeff[], // i : energy coeff. (5), fraction part, Q15 872b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // (frac_coeff and exp_coeff computed in 873b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // calc_filt_energies()) 874b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf1_exp_target_en, // i : exponent of target energy, Q0 875b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf1_frac_target_en, // i : fraction of target energy, Q15 876b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 877b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 gp_limit, // i : pitch gain limit 878b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 879b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *sf0_gain_pit, // o : Pitch gain, Q14 880b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *sf0_gain_cod, // o : Code gain, Q1 881b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 882b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *sf1_gain_pit, // o : Pitch gain, Q14 883b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *sf1_gain_cod // o : Code gain, Q1 884b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber) 885b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber{ 886b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber const Word16 *p; 887b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 i, index = 0; 888b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 tmp; 889b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 exp; 890b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf0_gcode0, sf1_gcode0; 891b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 g_pitch, g2_pitch, g_code, g2_code, g_pit_cod; 892b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 coeff[10], coeff_lo[10], exp_max[10]; // 0..4: sf0; 5..9: sf1 893b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word32 L_tmp, dist_min; 894b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 895b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *-------------------------------------------------------------------* 896b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * predicted codebook gain * 897b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * ~~~~~~~~~~~~~~~~~~~~~~~ * 898b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * gc0 = 2^exp_gcode0 + 2^frac_gcode0 * 899b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * * 900b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * gcode0 (Q14) = 2^14*2^frac_gcode0 = gc0 * 2^(14-exp_gcode0) * 901b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *-------------------------------------------------------------------* 902b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 903b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf0_gcode0 = extract_l(Pow2(14, sf0_frac_gcode0)); 904b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf1_gcode0 = extract_l(Pow2(14, sf1_frac_gcode0)); 905b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 906b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * For each subframe, the error energy (sum) to be minimized consists 907b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * of five terms, t[0..4]. 908b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * 909b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * t[0] = gp^2 * <y1 y1> 910b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * t[1] = -2*gp * <xn y1> 911b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * t[2] = gc^2 * <y2 y2> 912b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * t[3] = -2*gc * <xn y2> 913b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * t[4] = 2*gp*gc * <y1 y2> 914b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * 915b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 916b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // sf 0 917b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // determine the scaling exponent for g_code: ec = ec0 - 11 918b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp = sub(sf0_exp_gcode0, 11); 919b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 920b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // calculate exp_max[i] = s[i]-1 921b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp_max[0] = sub(sf0_exp_coeff[0], 13); 922b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp_max[1] = sub(sf0_exp_coeff[1], 14); 923b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp_max[2] = add(sf0_exp_coeff[2], add(15, shl(exp, 1))); 924b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp_max[3] = add(sf0_exp_coeff[3], exp); 925b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp_max[4] = add(sf0_exp_coeff[4], add(1, exp)); 926b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 927b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // sf 1 928b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // determine the scaling exponent for g_code: ec = ec0 - 11 929b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp = sub(sf1_exp_gcode0, 11); 930b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 931b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // calculate exp_max[i] = s[i]-1 932b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp_max[5] = sub(sf1_exp_coeff[0], 13); 933b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp_max[6] = sub(sf1_exp_coeff[1], 14); 934b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp_max[7] = add(sf1_exp_coeff[2], add(15, shl(exp, 1))); 935b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp_max[8] = add(sf1_exp_coeff[3], exp); 936b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp_max[9] = add(sf1_exp_coeff[4], add(1, exp)); 937b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 938b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *-------------------------------------------------------------------* 939b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * Gain search equalisation: * 940b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * ~~~~~~~~~~~~~~~~~~~~~~~~~ * 941b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * The MSE for the two subframes is weighted differently if there * 942b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * is a big difference in the corresponding target energies * 943b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *-------------------------------------------------------------------* 944b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 945b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // make the target energy exponents the same by de-normalizing the 946b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // fraction of the smaller one. This is necessary to be able to compare 947b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // them 948b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 949b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp = sf0_exp_target_en - sf1_exp_target_en; 950b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber if (exp > 0) 951b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 952b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf1_frac_target_en = shr (sf1_frac_target_en, exp); 953b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 954b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber else 955b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 956b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf0_frac_target_en = shl (sf0_frac_target_en, exp); 957b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 958b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 959b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // assume no change of exponents 960b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp = 0; 961b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 962b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // test for target energy difference; set exp to +1 or -1 to scale 963b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // up/down coefficients for sf 1 964b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 965b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp = shr_r (sf1_frac_target_en, 1); // tmp = ceil(0.5*en(sf1)) 966b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber if (sub (tmp, sf0_frac_target_en) > 0) // tmp > en(sf0)? 967b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 968b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // target_energy(sf1) > 2*target_energy(sf0) 969b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // -> scale up MSE(sf0) by 2 by adding 1 to exponents 0..4 970b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp = 1; 971b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 972b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber else 973b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 974b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp = shr (add (sf0_frac_target_en, 3), 2); // tmp=ceil(0.25*en(sf0)) 975b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber if (sub (tmp, sf1_frac_target_en) > 0) // tmp > en(sf1)? 976b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 977b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // target_energy(sf1) < 0.25*target_energy(sf0) 978b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // -> scale down MSE(sf0) by 0.5 by subtracting 1 from 979b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // coefficients 0..4 980b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp = -1; 981b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 982b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 983b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 984b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (i = 0; i < 5; i++) 985b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 986b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp_max[i] = add (exp_max[i], exp); 987b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 988b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 989b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *-------------------------------------------------------------------* 990b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * Find maximum exponent: * 991b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * ~~~~~~~~~~~~~~~~~~~~~~ * 992b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * * 993b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * For the sum operation, all terms must have the same scaling; * 994b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * that scaling should be low enough to prevent overflow. There- * 995b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * fore, the maximum scale is determined and all coefficients are * 996b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * re-scaled: * 997b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * * 998b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * exp = max(exp_max[i]) + 1; * 999b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * e = exp_max[i]-exp; e <= 0! * 1000b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * c[i] = c[i]*2^e * 1001b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *-------------------------------------------------------------------* 1002b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1003b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp = exp_max[0]; 1004b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (i = 1; i < 10; i++) 1005b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 1006b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber if (sub(exp_max[i], exp) > 0) 1007b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 1008b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp = exp_max[i]; 1009b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 1010b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 1011b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp = add(exp, 1); // To avoid overflow 1012b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1013b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber p = &sf0_frac_coeff[0]; 1014b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (i = 0; i < 5; i++) { 1015b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp = sub(exp, exp_max[i]); 1016b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = L_deposit_h(*p++); 1017b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = L_shr(L_tmp, tmp); 1018b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_Extract(L_tmp, &coeff[i], &coeff_lo[i]); 1019b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 1020b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber p = &sf1_frac_coeff[0]; 1021b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (; i < 10; i++) { 1022b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp = sub(exp, exp_max[i]); 1023b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = L_deposit_h(*p++); 1024b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = L_shr(L_tmp, tmp); 1025b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_Extract(L_tmp, &coeff[i], &coeff_lo[i]); 1026b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 1027b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1028b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber //-------------------------------------------------------------------* 1029b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * Codebook search: * 1030b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * ~~~~~~~~~~~~~~~~ * 1031b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * * 1032b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * For each pair (g_pitch, g_fac) in the table calculate the * 1033b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * terms t[0..4] and sum them up; the result is the mean squared * 1034b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * error for the quantized gains from the table. The index for the * 1035b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * minimum MSE is stored and finally used to retrieve the quantized * 1036b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * gains * 1037b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *------------------------------------------------------------------- 1038b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1039b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // start with "infinite" MSE 1040b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber dist_min = MAX_32; 1041b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1042b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber p = &table_gain_MR475[0]; 1043b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1044b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (i = 0; i < MR475_VQ_SIZE; i++) 1045b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 1046b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // subframe 0 (and 2) calculations 1047b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g_pitch = *p++; 1048b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g_code = *p++; 1049b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1050b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g_code = mult(g_code, sf0_gcode0); 1051b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g2_pitch = mult(g_pitch, g_pitch); 1052b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g2_code = mult(g_code, g_code); 1053b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g_pit_cod = mult(g_code, g_pitch); 1054b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1055b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = Mpy_32_16( coeff[0], coeff_lo[0], g2_pitch); 1056b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = Mac_32_16(L_tmp, coeff[1], coeff_lo[1], g_pitch); 1057b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = Mac_32_16(L_tmp, coeff[2], coeff_lo[2], g2_code); 1058b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = Mac_32_16(L_tmp, coeff[3], coeff_lo[3], g_code); 1059b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = Mac_32_16(L_tmp, coeff[4], coeff_lo[4], g_pit_cod); 1060b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1061b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp = sub (g_pitch, gp_limit); 1062b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1063b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // subframe 1 (and 3) calculations 1064b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g_pitch = *p++; 1065b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g_code = *p++; 1066b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1067b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber if (tmp <= 0 && sub(g_pitch, gp_limit) <= 0) 1068b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 1069b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g_code = mult(g_code, sf1_gcode0); 1070b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g2_pitch = mult(g_pitch, g_pitch); 1071b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g2_code = mult(g_code, g_code); 1072b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g_pit_cod = mult(g_code, g_pitch); 1073b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1074b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = Mac_32_16(L_tmp, coeff[5], coeff_lo[5], g2_pitch); 1075b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = Mac_32_16(L_tmp, coeff[6], coeff_lo[6], g_pitch); 1076b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = Mac_32_16(L_tmp, coeff[7], coeff_lo[7], g2_code); 1077b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = Mac_32_16(L_tmp, coeff[8], coeff_lo[8], g_code); 1078b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = Mac_32_16(L_tmp, coeff[9], coeff_lo[9], g_pit_cod); 1079b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1080b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // store table index if MSE for this index is lower 1081b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber than the minimum MSE seen so far 1082b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber if (L_sub(L_tmp, dist_min) < (Word32) 0) 1083b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 1084b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber dist_min = L_tmp; 1085b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber index = i; 1086b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 1087b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 1088b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 1089b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1090b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *------------------------------------------------------------------* 1091b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * read quantized gains and update MA predictor memories * 1092b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ * 1093b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *------------------------------------------------------------------* 1094b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1095b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // for subframe 0, the pre-calculated gcode0/exp_gcode0 are the same 1096b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // as those calculated from the "real" predictor using quantized gains 1097b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp = shl(index, 2); 1098b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber MR475_quant_store_results(pred_st, 1099b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber &table_gain_MR475[tmp], 1100b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf0_gcode0, 1101b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf0_exp_gcode0, 1102b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf0_gain_pit, 1103b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf0_gain_cod); 1104b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1105b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // calculate new predicted gain for subframe 1 (this time using 1106b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber // the real, quantized gains) 1107b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber gc_pred(pred_st, MR475, sf1_code_nosharp, 1108b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber &sf1_exp_gcode0, &sf1_frac_gcode0, 1109b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber &sf0_exp_gcode0, &sf0_gcode0); // last two args are dummy 1110b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf1_gcode0 = extract_l(Pow2(14, sf1_frac_gcode0)); 1111b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1112b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp = add (tmp, 2); 1113b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber MR475_quant_store_results(pred_st, 1114b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber &table_gain_MR475[tmp], 1115b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf1_gcode0, 1116b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf1_exp_gcode0, 1117b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf1_gain_pit, 1118b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf1_gain_cod); 1119b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1120b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber return index; 1121b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber} 1122b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1123b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 1124b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber RESOURCES USED [optional] 1125b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1126b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber When the code is written for a specific target processor the 1127b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber the resources used should be documented below. 1128b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1129b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber HEAP MEMORY USED: x bytes 1130b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1131b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber STACK MEMORY USED: x bytes 1132b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1133b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber CLOCK CYCLES: (cycle count equation for this function) + (variable 1134b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber used to represent cycle count for each subroutine 1135b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber called) 1136b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber where: (cycle count variable) = cycle count for [subroutine 1137b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber name] 1138b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1139b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 1140b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber CAUTION [optional] 1141b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber [State any special notes, constraints or cautions for users of this function] 1142b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1143b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber------------------------------------------------------------------------------ 1144b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber*/ 1145b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1146b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas HuberWord16 MR475_gain_quant( /* o : index of quantization. */ 1147b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber gc_predState *pred_st, /* i/o: gain predictor state struct */ 1148b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1149b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* data from subframe 0 (or 2) */ 1150b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf0_exp_gcode0, /* i : predicted CB gain (exponent), Q0 */ 1151b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf0_frac_gcode0, /* i : predicted CB gain (fraction), Q15 */ 1152b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf0_exp_coeff[], /* i : energy coeff. (5), exponent part, Q0 */ 1153b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf0_frac_coeff[], /* i : energy coeff. (5), fraction part, Q15 */ 1154b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* (frac_coeff and exp_coeff computed in */ 1155b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* calc_filt_energies()) */ 1156b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf0_exp_target_en, /* i : exponent of target energy, Q0 */ 1157b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf0_frac_target_en, /* i : fraction of target energy, Q15 */ 1158b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1159b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* data from subframe 1 (or 3) */ 1160b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf1_code_nosharp[], /* i : innovative codebook vector (L_SUBFR) */ 1161b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* (whithout pitch sharpening) */ 1162b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf1_exp_gcode0, /* i : predicted CB gain (exponent), Q0 */ 1163b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf1_frac_gcode0, /* i : predicted CB gain (fraction), Q15 */ 1164b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf1_exp_coeff[], /* i : energy coeff. (5), exponent part, Q0 */ 1165b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf1_frac_coeff[], /* i : energy coeff. (5), fraction part, Q15 */ 1166b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* (frac_coeff and exp_coeff computed in */ 1167b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* calc_filt_energies()) */ 1168b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf1_exp_target_en, /* i : exponent of target energy, Q0 */ 1169b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf1_frac_target_en, /* i : fraction of target energy, Q15 */ 1170b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1171b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 gp_limit, /* i : pitch gain limit */ 1172b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1173b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *sf0_gain_pit, /* o : Pitch gain, Q14 */ 1174b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *sf0_gain_cod, /* o : Code gain, Q1 */ 1175b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1176b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *sf1_gain_pit, /* o : Pitch gain, Q14 */ 1177b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 *sf1_gain_cod, /* o : Code gain, Q1 */ 1178b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Flag *pOverflow /* o : overflow indicator */ 1179b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber) 1180b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber{ 1181b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber const Word16 *p; 1182b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 i; 1183b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 index = 0; 1184b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 tmp; 1185b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 exp; 1186b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf0_gcode0; 1187b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 sf1_gcode0; 1188b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 g_pitch; 1189b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 g2_pitch; 1190b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 g_code; 1191b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 g2_code; 1192b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 g_pit_cod; 1193b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 coeff[10]; 1194b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 coeff_lo[10]; 1195b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word16 exp_max[10]; /* 0..4: sf0; 5..9: sf1 */ 1196b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word32 L_tmp; 1197b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Word32 dist_min; 1198b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1199b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /*-------------------------------------------------------------------* 1200b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * predicted codebook gain * 1201b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * ~~~~~~~~~~~~~~~~~~~~~~~ * 1202b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * gc0 = 2^exp_gcode0 + 2^frac_gcode0 * 1203b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * * 1204b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * gcode0 (Q14) = 2^14*2^frac_gcode0 = gc0 * 2^(14-exp_gcode0) * 1205b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *-------------------------------------------------------------------*/ 1206b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1207b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf0_gcode0 = (Word16)(Pow2(14, sf0_frac_gcode0, pOverflow)); 1208b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf1_gcode0 = (Word16)(Pow2(14, sf1_frac_gcode0, pOverflow)); 1209b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1210b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* 1211b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * For each subframe, the error energy (sum) to be minimized consists 1212b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * of five terms, t[0..4]. 1213b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * 1214b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * t[0] = gp^2 * <y1 y1> 1215b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * t[1] = -2*gp * <xn y1> 1216b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * t[2] = gc^2 * <y2 y2> 1217b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * t[3] = -2*gc * <xn y2> 1218b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * t[4] = 2*gp*gc * <y1 y2> 1219b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * 1220b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber */ 1221b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1222b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* sf 0 */ 1223b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* determine the scaling exponent for g_code: ec = ec0 - 11 */ 1224b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp = sf0_exp_gcode0 - 11; 1225b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1226b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* calculate exp_max[i] = s[i]-1 */ 1227b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp_max[0] = (sf0_exp_coeff[0] - 13); 1228b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp_max[1] = (sf0_exp_coeff[1] - 14); 1229b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp_max[2] = (sf0_exp_coeff[2] + (15 + (exp << 1))); 1230b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp_max[3] = (sf0_exp_coeff[3] + exp); 1231b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp_max[4] = (sf0_exp_coeff[4] + (1 + exp)); 1232b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1233b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* sf 1 */ 1234b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* determine the scaling exponent for g_code: ec = ec0 - 11 */ 1235b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp = sf1_exp_gcode0 - 11; 1236b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1237b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* calculate exp_max[i] = s[i]-1 */ 1238b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp_max[5] = (sf1_exp_coeff[0] - 13); 1239b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp_max[6] = (sf1_exp_coeff[1] - 14); 1240b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp_max[7] = (sf1_exp_coeff[2] + (15 + (exp << 1))); 1241b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp_max[8] = (sf1_exp_coeff[3] + exp); 1242b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp_max[9] = (sf1_exp_coeff[4] + (1 + exp)); 1243b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1244b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /*-------------------------------------------------------------------* 1245b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * Gain search equalisation: * 1246b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * ~~~~~~~~~~~~~~~~~~~~~~~~~ * 1247b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * The MSE for the two subframes is weighted differently if there * 1248b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * is a big difference in the corresponding target energies * 1249b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *-------------------------------------------------------------------*/ 1250b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1251b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* make the target energy exponents the same by de-normalizing the 1252b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber fraction of the smaller one. This is necessary to be able to compare 1253b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber them 1254b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber */ 1255b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp = sf0_exp_target_en - sf1_exp_target_en; 1256b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber if (exp > 0) 1257b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 1258b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf1_frac_target_en >>= exp; 1259b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 1260b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber else 1261b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 1262b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf0_frac_target_en >>= (-exp); 1263b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 1264b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1265b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* assume no change of exponents */ 1266b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp = 0; 1267b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1268b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* test for target energy difference; set exp to +1 or -1 to scale 1269b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * up/down coefficients for sf 1 1270b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber */ 1271b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp = shr_r(sf1_frac_target_en, 1, pOverflow); /* tmp = ceil(0.5*en(sf1)) */ 1272b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1273b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber if (tmp > sf0_frac_target_en) /* tmp > en(sf0)? */ 1274b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 1275b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* 1276b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * target_energy(sf1) > 2*target_energy(sf0) 1277b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * -> scale up MSE(sf0) by 2 by adding 1 to exponents 0..4 1278b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber */ 1279b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp = 1; 1280b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 1281b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber else 1282b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 1283b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp = ((sf0_frac_target_en + 3) >> 2); /* tmp=ceil(0.25*en(sf0)) */ 1284b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1285b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber if (tmp > sf1_frac_target_en) /* tmp > en(sf1)? */ 1286b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 1287b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* 1288b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * target_energy(sf1) < 0.25*target_energy(sf0) 1289b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * -> scale down MSE(sf0) by 0.5 by subtracting 1 from 1290b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * coefficients 0..4 1291b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber */ 1292b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp = -1; 1293b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 1294b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 1295b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1296b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (i = 0; i < 5; i++) 1297b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 1298b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp_max[i] += exp; 1299b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 1300b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1301b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /*-------------------------------------------------------------------* 1302b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * Find maximum exponent: * 1303b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * ~~~~~~~~~~~~~~~~~~~~~~ * 1304b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * * 1305b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * For the sum operation, all terms must have the same scaling; * 1306b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * that scaling should be low enough to prevent overflow. There- * 1307b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * fore, the maximum scale is determined and all coefficients are * 1308b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * re-scaled: * 1309b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * * 1310b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * exp = max(exp_max[i]) + 1; * 1311b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * e = exp_max[i]-exp; e <= 0! * 1312b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * c[i] = c[i]*2^e * 1313b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *-------------------------------------------------------------------*/ 1314b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1315b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp = exp_max[0]; 1316b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (i = 9; i > 0; i--) 1317b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 1318b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber if (exp_max[i] > exp) 1319b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 1320b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp = exp_max[i]; 1321b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 1322b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 1323b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber exp++; /* To avoid overflow */ 1324b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1325b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber p = &sf0_frac_coeff[0]; 1326b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (i = 0; i < 5; i++) 1327b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 1328b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp = (exp - exp_max[i]); 1329b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = ((Word32)(*p++) << 16); 1330b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = L_shr(L_tmp, tmp, pOverflow); 1331b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber coeff[i] = (Word16)(L_tmp >> 16); 1332b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber coeff_lo[i] = (Word16)((L_tmp >> 1) - ((L_tmp >> 16) << 15)); 1333b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 1334b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber p = &sf1_frac_coeff[0]; 1335b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (; i < 10; i++) 1336b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 1337b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp = exp - exp_max[i]; 1338b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = ((Word32)(*p++) << 16); 1339b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = L_shr(L_tmp, tmp, pOverflow); 1340b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber coeff[i] = (Word16)(L_tmp >> 16); 1341b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber coeff_lo[i] = (Word16)((L_tmp >> 1) - ((L_tmp >> 16) << 15)); 1342b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 1343b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1344b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1345b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /*-------------------------------------------------------------------* 1346b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * Codebook search: * 1347b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * ~~~~~~~~~~~~~~~~ * 1348b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * * 1349b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * For each pair (g_pitch, g_fac) in the table calculate the * 1350b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * terms t[0..4] and sum them up; the result is the mean squared * 1351b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * error for the quantized gains from the table. The index for the * 1352b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * minimum MSE is stored and finally used to retrieve the quantized * 1353b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * gains * 1354b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *-------------------------------------------------------------------*/ 1355b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1356b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* start with "infinite" MSE */ 1357b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber dist_min = MAX_32; 1358b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1359b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber p = &table_gain_MR475[0]; 1360b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1361b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber for (i = 0; i < MR475_VQ_SIZE; i++) 1362b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 1363b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* subframe 0 (and 2) calculations */ 1364b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g_pitch = *p++; 1365b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g_code = *p++; 1366b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1367b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* Need to be there OKA */ 1368b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g_code = (Word16)(((Word32) g_code * sf0_gcode0) >> 15); 1369b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g2_pitch = (Word16)(((Word32) g_pitch * g_pitch) >> 15); 1370b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g2_code = (Word16)(((Word32) g_code * g_code) >> 15); 1371b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g_pit_cod = (Word16)(((Word32) g_code * g_pitch) >> 15); 1372b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1373b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1374b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp = Mpy_32_16(coeff[0], coeff_lo[0], g2_pitch, pOverflow) + 1375b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Mpy_32_16(coeff[1], coeff_lo[1], g_pitch, pOverflow) + 1376b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Mpy_32_16(coeff[2], coeff_lo[2], g2_code, pOverflow) + 1377b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Mpy_32_16(coeff[3], coeff_lo[3], g_code, pOverflow) + 1378b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Mpy_32_16(coeff[4], coeff_lo[4], g_pit_cod, pOverflow); 1379b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1380b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp = (g_pitch - gp_limit); 1381b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1382b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* subframe 1 (and 3) calculations */ 1383b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g_pitch = *p++; 1384b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g_code = *p++; 1385b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1386b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber if ((tmp <= 0) && (g_pitch <= gp_limit)) 1387b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 1388b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g_code = (Word16)(((Word32) g_code * sf1_gcode0) >> 15); 1389b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g2_pitch = (Word16)(((Word32) g_pitch * g_pitch) >> 15); 1390b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g2_code = (Word16)(((Word32) g_code * g_code) >> 15); 1391b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber g_pit_cod = (Word16)(((Word32) g_code * g_pitch) >> 15); 1392b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1393b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber L_tmp += (Mpy_32_16(coeff[5], coeff_lo[5], g2_pitch, pOverflow) + 1394b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Mpy_32_16(coeff[6], coeff_lo[6], g_pitch, pOverflow) + 1395b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Mpy_32_16(coeff[7], coeff_lo[7], g2_code, pOverflow) + 1396b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Mpy_32_16(coeff[8], coeff_lo[8], g_code, pOverflow) + 1397b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber Mpy_32_16(coeff[9], coeff_lo[9], g_pit_cod, pOverflow)); 1398b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1399b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* store table index if MSE for this index is lower 1400b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber than the minimum MSE seen so far */ 1401b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber if (L_tmp < dist_min) 1402b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber { 1403b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber dist_min = L_tmp; 1404b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber index = i; 1405b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 1406b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 1407b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber } 1408b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1409b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /*------------------------------------------------------------------* 1410b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * read quantized gains and update MA predictor memories * 1411b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber * ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ * 1412b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber *------------------------------------------------------------------*/ 1413b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1414b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* for subframe 0, the pre-calculated gcode0/exp_gcode0 are the same 1415b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber as those calculated from the "real" predictor using quantized gains */ 1416b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp = index << 2; 1417b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber MR475_quant_store_results(pred_st, 1418b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber &table_gain_MR475[tmp], 1419b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf0_gcode0, 1420b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf0_exp_gcode0, 1421b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf0_gain_pit, 1422b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf0_gain_cod, 1423b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber pOverflow); 1424b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1425b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber /* calculate new predicted gain for subframe 1 (this time using 1426b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber the real, quantized gains) */ 1427b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber gc_pred(pred_st, MR475, sf1_code_nosharp, 1428b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber &sf1_exp_gcode0, &sf1_frac_gcode0, 1429b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber &sf0_exp_gcode0, &sf0_gcode0, /* dummy args */ 1430b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber pOverflow); 1431b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1432b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf1_gcode0 = (Word16)(Pow2(14, sf1_frac_gcode0, pOverflow)); 1433b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1434b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber tmp += 2; 1435b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber MR475_quant_store_results( 1436b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber pred_st, 1437b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber &table_gain_MR475[tmp], 1438b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf1_gcode0, 1439b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf1_exp_gcode0, 1440b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf1_gain_pit, 1441b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber sf1_gain_cod, 1442b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber pOverflow); 1443b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber 1444b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber return(index); 1445b841f14f8e51f2365945281fbfa54ef6a1b1b5a6Andreas Huber} 1446