1 /* ------------------------------------------------------------------ 2 * Copyright (C) 1998-2009 PacketVideo 3 * 4 * Licensed under the Apache License, Version 2.0 (the "License"); 5 * you may not use this file except in compliance with the License. 6 * You may obtain a copy of the License at 7 * 8 * http://www.apache.org/licenses/LICENSE-2.0 9 * 10 * Unless required by applicable law or agreed to in writing, software 11 * distributed under the License is distributed on an "AS IS" BASIS, 12 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either 13 * express or implied. 14 * See the License for the specific language governing permissions 15 * and limitations under the License. 16 * ------------------------------------------------------------------- 17 */ 18 /**************************************************************************************** 19 Portions of this file are derived from the following 3GPP standard: 20 21 3GPP TS 26.173 22 ANSI-C code for the Adaptive Multi-Rate - Wideband (AMR-WB) speech codec 23 Available from http://www.3gpp.org 24 25 (C) 2007, 3GPP Organizational Partners (ARIB, ATIS, CCSA, ETSI, TTA, TTC) 26 Permission to distribute, modify and use this file under the standard license 27 terms listed above has been obtained from the copyright holder. 28 ****************************************************************************************/ 29 /* 30 ------------------------------------------------------------------------------ 31 32 33 34 Filename: pred_lt4.cpp 35 36 Date: 05/08/2004 37 38 ------------------------------------------------------------------------------ 39 REVISION HISTORY 40 41 42 Description: 43 44 ------------------------------------------------------------------------------ 45 INPUT AND OUTPUT DEFINITIONS 46 47 int16 signal[], input signal / output is divided by 16 48 int16 lg, lenght of signal 49 int16 mem[] in/out: memory (size=30) 50 int16 x[] scratch mem ( size= 60) 51 52 ------------------------------------------------------------------------------ 53 FUNCTION DESCRIPTION 54 55 Compute the result of long term prediction with fractionnal 56 interpolation of resolution 1/4. 57 58 On return exc[0..L_subfr-1] contains the interpolated signal 59 (adaptive codebook excitation) 60 61 62 ------------------------------------------------------------------------------ 63 REQUIREMENTS 64 65 66 ------------------------------------------------------------------------------ 67 REFERENCES 68 69 ------------------------------------------------------------------------------ 70 PSEUDO-CODE 71 72 ------------------------------------------------------------------------------ 73 */ 74 75 76 /*---------------------------------------------------------------------------- 77 ; INCLUDES 78 ----------------------------------------------------------------------------*/ 79 80 #include "pv_amr_wb_type_defs.h" 81 #include "pvamrwbdecoder_basic_op.h" 82 #include "pvamrwbdecoder_acelp.h" 83 84 /*---------------------------------------------------------------------------- 85 ; MACROS 86 ; Define module specific macros here 87 ----------------------------------------------------------------------------*/ 88 89 90 /*---------------------------------------------------------------------------- 91 ; DEFINES 92 ; Include all pre-processor statements here. Include conditional 93 ; compile variables also. 94 ----------------------------------------------------------------------------*/ 95 96 #define UP_SAMP 4 97 #define L_INTERPOL2 16 98 99 /*---------------------------------------------------------------------------- 100 ; LOCAL FUNCTION DEFINITIONS 101 ; Function Prototype declaration 102 ----------------------------------------------------------------------------*/ 103 104 /*---------------------------------------------------------------------------- 105 ; LOCAL STORE/BUFFER/POINTER DEFINITIONS 106 ; Variable declaration - defined here and used outside this module 107 ----------------------------------------------------------------------------*/ 108 109 /* 1/4 resolution interpolation filter (-3 dB at 0.856*fs/2) in Q14 */ 110 111 112 const int16 inter4_2[UP_SAMP][ 2*L_INTERPOL2] = 113 { 114 { 115 0, -2, 4, -2, -10, 38, 116 -88, 165, -275, 424, -619, 871, 117 -1207, 1699, -2598, 5531, 14031, -2147, 118 780, -249, -16, 153, -213, 226, 119 -209, 175, -133, 91, -55, 28, 120 -10, 2 121 }, 122 { 123 1, -7, 19, -33, 47, -52, 124 43, -9, -60, 175, -355, 626, 125 -1044, 1749, -3267, 10359, 10359, -3267, 126 1749, -1044, 626, -355, 175, -60, 127 -9, 43, -52, 47, -33, 19, 128 -7, 1 129 }, 130 { 131 2, -10, 28, -55, 91, -133, 132 175, -209, 226, -213, 153, -16, 133 -249, 780, -2147, 14031, 5531, -2598, 134 1699, -1207, 871, -619, 424, -275, 135 165, -88, 38, -10, -2, 4, 136 -2, 0 137 }, 138 { 139 1, -7, 22, -49, 92, -153, 140 231, -325, 431, -544, 656, -762, 141 853, -923, 968, 15401, 968, -923, 142 853, -762, 656, -544, 431, -325, 143 231, -153, 92, -49, 22, -7, 144 1, 0 145 } 146 }; 147 148 /*---------------------------------------------------------------------------- 149 ; EXTERNAL FUNCTION REFERENCES 150 ; Declare functions defined elsewhere and referenced in this module 151 ----------------------------------------------------------------------------*/ 152 153 /*---------------------------------------------------------------------------- 154 ; EXTERNAL GLOBAL STORE/BUFFER/POINTER REFERENCES 155 ; Declare variables used in this module but defined elsewhere 156 ----------------------------------------------------------------------------*/ 157 158 /*---------------------------------------------------------------------------- 159 ; FUNCTION CODE 160 ----------------------------------------------------------------------------*/ 161 162 void Pred_lt4( 163 int16 exc[], /* in/out: excitation buffer */ 164 int16 T0, /* input : integer pitch lag */ 165 int16 frac, /* input : fraction of lag */ 166 int16 L_subfr /* input : subframe size */ 167 ) 168 { 169 int16 i, j, *pt_exc; 170 int32 L_sum1; 171 int32 L_sum2; 172 int32 L_sum3; 173 int32 L_sum4; 174 pt_exc = &exc[-T0]; 175 176 const int16 *pt_inter4_2; 177 178 frac = -frac; 179 180 if (frac < 0) 181 { 182 frac += UP_SAMP; 183 pt_exc--; 184 185 } 186 pt_exc -= (L_INTERPOL2 - 1); 187 188 pt_inter4_2 = inter4_2[UP_SAMP-1 - frac]; 189 190 for (j = 0; j < (L_subfr >> 2); j++) 191 { 192 193 L_sum1 = 0x00002000; /* pre-roundig */ 194 L_sum2 = 0x00002000; 195 L_sum3 = 0x00002000; 196 L_sum4 = 0x00002000; 197 198 for (i = 0; i < L_INTERPOL2 << 1; i += 4) 199 { 200 int16 tmp1 = pt_exc[i ]; 201 int16 tmp2 = pt_exc[i+1]; 202 int16 tmp3 = pt_exc[i+2]; 203 204 205 L_sum1 = fxp_mac_16by16(tmp1, pt_inter4_2[i ], L_sum1); 206 L_sum2 = fxp_mac_16by16(tmp2, pt_inter4_2[i ], L_sum2); 207 L_sum1 = fxp_mac_16by16(tmp2, pt_inter4_2[i+1], L_sum1); 208 L_sum2 = fxp_mac_16by16(tmp3, pt_inter4_2[i+1], L_sum2); 209 L_sum3 = fxp_mac_16by16(tmp3, pt_inter4_2[i ], L_sum3); 210 L_sum1 = fxp_mac_16by16(tmp3, pt_inter4_2[i+2], L_sum1); 211 212 tmp1 = pt_exc[i+3]; 213 tmp2 = pt_exc[i+4]; 214 215 L_sum4 = fxp_mac_16by16(tmp1, pt_inter4_2[i ], L_sum4); 216 L_sum3 = fxp_mac_16by16(tmp1, pt_inter4_2[i+1], L_sum3); 217 L_sum2 = fxp_mac_16by16(tmp1, pt_inter4_2[i+2], L_sum2); 218 L_sum1 = fxp_mac_16by16(tmp1, pt_inter4_2[i+3], L_sum1); 219 L_sum4 = fxp_mac_16by16(tmp2, pt_inter4_2[i+1], L_sum4); 220 L_sum2 = fxp_mac_16by16(tmp2, pt_inter4_2[i+3], L_sum2); 221 L_sum3 = fxp_mac_16by16(tmp2, pt_inter4_2[i+2], L_sum3); 222 223 tmp1 = pt_exc[i+5]; 224 tmp2 = pt_exc[i+6]; 225 226 L_sum4 = fxp_mac_16by16(tmp1, pt_inter4_2[i+2], L_sum4); 227 L_sum3 = fxp_mac_16by16(tmp1, pt_inter4_2[i+3], L_sum3); 228 L_sum4 = fxp_mac_16by16(tmp2, pt_inter4_2[i+3], L_sum4); 229 230 } 231 232 233 234 exc[(j<<2)] = (int16)(L_sum1 >> 14); 235 exc[(j<<2)+1] = (int16)(L_sum2 >> 14); 236 exc[(j<<2)+2] = (int16)(L_sum3 >> 14); 237 exc[(j<<2)+3] = (int16)(L_sum4 >> 14); 238 239 pt_exc += 4; 240 241 } 242 243 if (L_subfr&1) 244 { 245 L_sum1 = 0x00002000; 246 247 for (i = 0; i < 2*L_INTERPOL2; i += 4) 248 { 249 int16 tmp1 = pt_exc[i ]; 250 int16 tmp2 = pt_exc[i+1]; 251 L_sum1 = fxp_mac_16by16(tmp1, pt_inter4_2[i ], L_sum1); 252 L_sum1 = fxp_mac_16by16(tmp2, pt_inter4_2[i+1], L_sum1); 253 tmp1 = pt_exc[i+2]; 254 tmp2 = pt_exc[i+3]; 255 L_sum1 = fxp_mac_16by16(tmp1, pt_inter4_2[i+2], L_sum1); 256 L_sum1 = fxp_mac_16by16(tmp2, pt_inter4_2[i+3], L_sum1); 257 258 } 259 260 exc[(j<<2)] = (int16)((L_sum1) >> 14); 261 262 } 263 264 265 return; 266 } 267 268