| OLD | NEW |
| 1 /* | 1 /* |
| 2 * Copyright (c) 2011 The WebRTC project authors. All Rights Reserved. | 2 * Copyright (c) 2011 The WebRTC project authors. All Rights Reserved. |
| 3 * | 3 * |
| 4 * Use of this source code is governed by a BSD-style license | 4 * Use of this source code is governed by a BSD-style license |
| 5 * that can be found in the LICENSE file in the root of the source | 5 * that can be found in the LICENSE file in the root of the source |
| 6 * tree. An additional intellectual property rights grant can be found | 6 * tree. An additional intellectual property rights grant can be found |
| 7 * in the file PATENTS. All contributing project authors may | 7 * in the file PATENTS. All contributing project authors may |
| 8 * be found in the AUTHORS file in the root of the source tree. | 8 * be found in the AUTHORS file in the root of the source tree. |
| 9 */ | 9 */ |
| 10 | 10 |
| (...skipping 256 matching lines...) Expand 10 before | Expand all | Expand 10 after Loading... |
| 267 sh = WebRtcSpl_NormW32(tmp32); // tmp32 is the gain | 267 sh = WebRtcSpl_NormW32(tmp32); // tmp32 is the gain |
| 268 den16 = (int16_t) WEBRTC_SPL_SHIFT_W32(tmp32, sh-16); //Q(27+sh-16) = Q(sh+1
1) (all 16 bits are value bits) | 268 den16 = (int16_t) WEBRTC_SPL_SHIFT_W32(tmp32, sh-16); //Q(27+sh-16) = Q(sh+1
1) (all 16 bits are value bits) |
| 269 inv_gain32 = WebRtcSpl_DivW32W16((int32_t)2147483647, den16); // 1/gain in Q
31/Q(sh+11) = Q(20-sh) | 269 inv_gain32 = WebRtcSpl_DivW32W16((int32_t)2147483647, den16); // 1/gain in Q
31/Q(sh+11) = Q(20-sh) |
| 270 | 270 |
| 271 //initial conditions | 271 //initial conditions |
| 272 inv_gain16 = (int16_t)(inv_gain32 >> 2); // 1/gain in Q(20-sh-2) = Q(18-sh) | 272 inv_gain16 = (int16_t)(inv_gain32 >> 2); // 1/gain in Q(20-sh-2) = Q(18-sh) |
| 273 | 273 |
| 274 for (i=0;i<HALF_SUBFRAMELEN;i++) | 274 for (i=0;i<HALF_SUBFRAMELEN;i++) |
| 275 { | 275 { |
| 276 | 276 |
| 277 tmp32 = lat_inQ25[i + temp1] << 1; // Q25->Q26 | 277 tmp32 = lat_inQ25[i + temp1] * (1 << 1); // Q25->Q26 |
| 278 tmp32 = WEBRTC_SPL_MUL_16_32_RSFT16(inv_gain16, tmp32); //lat_in[]*inv_gai
n in (Q(18-sh)*Q26)>>16 = Q(28-sh) | 278 tmp32 = WEBRTC_SPL_MUL_16_32_RSFT16(inv_gain16, tmp32); //lat_in[]*inv_gai
n in (Q(18-sh)*Q26)>>16 = Q(28-sh) |
| 279 tmp32 = WEBRTC_SPL_SHIFT_W32(tmp32, -(28-sh)); // lat_in[]*inv_gain in Q0 | 279 tmp32 = WEBRTC_SPL_SHIFT_W32(tmp32, -(28-sh)); // lat_in[]*inv_gain in Q0 |
| 280 | 280 |
| 281 ARfQ0vec[i] = (int16_t)WebRtcSpl_SatW32ToW16(tmp32); // Q0 | 281 ARfQ0vec[i] = (int16_t)WebRtcSpl_SatW32ToW16(tmp32); // Q0 |
| 282 } | 282 } |
| 283 | 283 |
| 284 // Get the state of f & g for the first input, for all orders. | 284 // Get the state of f & g for the first input, for all orders. |
| 285 for (i = orderCoef; i > 0; i--) | 285 for (i = orderCoef; i > 0; i--) |
| 286 { | 286 { |
| 287 tmp32 = (cthQ15[i - 1] * ARfQ0vec[0] - sthQ15[i - 1] * stateGQ0[i - 1] + | 287 tmp32 = (cthQ15[i - 1] * ARfQ0vec[0] - sthQ15[i - 1] * stateGQ0[i - 1] + |
| (...skipping 19 matching lines...) Expand all Loading... |
| 307 /* cannot use memcpy in the following */ | 307 /* cannot use memcpy in the following */ |
| 308 | 308 |
| 309 for (i=0;i<ord_1;i++) | 309 for (i=0;i<ord_1;i++) |
| 310 { | 310 { |
| 311 stateGQ0[i] = ARgQ0vec[i]; | 311 stateGQ0[i] = ARgQ0vec[i]; |
| 312 } | 312 } |
| 313 } | 313 } |
| 314 | 314 |
| 315 return; | 315 return; |
| 316 } | 316 } |
| OLD | NEW |