| 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  | 
|   11  |   11  | 
|   12 /* |   12 /* | 
|   13  * This file contains the function WebRtcSpl_LevinsonDurbin(). |   13  * This file contains the function WebRtcSpl_LevinsonDurbin(). | 
|   14  * The description header can be found in signal_processing_library.h |   14  * The description header can be found in signal_processing_library.h | 
|   15  * |   15  * | 
|   16  */ |   16  */ | 
|   17  |   17  | 
|   18 #include "webrtc/common_audio/signal_processing/include/signal_processing_librar
     y.h" |   18 #include "webrtc/common_audio/signal_processing/include/signal_processing_librar
     y.h" | 
 |   19 #include "webrtc/rtc_base/sanitizer.h" | 
|   19  |   20  | 
|   20 #define SPL_LEVINSON_MAXORDER 20 |   21 #define SPL_LEVINSON_MAXORDER 20 | 
|   21  |   22  | 
|   22 int16_t WebRtcSpl_LevinsonDurbin(const int32_t* R, int16_t* A, int16_t* K, |   23 int16_t RTC_NO_SANITIZE("signed-integer-overflow")  // bugs.webrtc.org/5486 | 
|   23                                  size_t order) |   24 WebRtcSpl_LevinsonDurbin(const int32_t* R, int16_t* A, int16_t* K, | 
 |   25                          size_t order) | 
|   24 { |   26 { | 
|   25     size_t i, j; |   27     size_t i, j; | 
|   26     // Auto-correlation coefficients in high precision |   28     // Auto-correlation coefficients in high precision | 
|   27     int16_t R_hi[SPL_LEVINSON_MAXORDER + 1], R_low[SPL_LEVINSON_MAXORDER + 1]; |   29     int16_t R_hi[SPL_LEVINSON_MAXORDER + 1], R_low[SPL_LEVINSON_MAXORDER + 1]; | 
|   28     // LPC coefficients in high precision |   30     // LPC coefficients in high precision | 
|   29     int16_t A_hi[SPL_LEVINSON_MAXORDER + 1], A_low[SPL_LEVINSON_MAXORDER + 1]; |   31     int16_t A_hi[SPL_LEVINSON_MAXORDER + 1], A_low[SPL_LEVINSON_MAXORDER + 1]; | 
|   30     // LPC coefficients for next iteration |   32     // LPC coefficients for next iteration | 
|   31     int16_t A_upd_hi[SPL_LEVINSON_MAXORDER + 1], A_upd_low[SPL_LEVINSON_MAXORDER
      + 1]; |   33     int16_t A_upd_hi[SPL_LEVINSON_MAXORDER + 1], A_upd_low[SPL_LEVINSON_MAXORDER
      + 1]; | 
|   32     // Reflection coefficient in high precision |   34     // Reflection coefficient in high precision | 
|   33     int16_t K_hi, K_low; |   35     int16_t K_hi, K_low; | 
|   34     // Prediction gain Alpha in high precision and with scale factor |   36     // Prediction gain Alpha in high precision and with scale factor | 
|   35     int16_t Alpha_hi, Alpha_low, Alpha_exp; |   37     int16_t Alpha_hi, Alpha_low, Alpha_exp; | 
|   36     int16_t tmp_hi, tmp_low; |   38     int16_t tmp_hi, tmp_low; | 
|   37     int32_t temp1W32, temp2W32, temp3W32; |   39     int32_t temp1W32, temp2W32, temp3W32; | 
|   38     int16_t norm; |   40     int16_t norm; | 
|   39  |   41  | 
|   40     // Normalize the autocorrelation R[0]...R[order+1] |   42     // Normalize the autocorrelation R[0]...R[order+1] | 
|   41  |   43  | 
|   42     norm = WebRtcSpl_NormW32(R[0]); |   44     norm = WebRtcSpl_NormW32(R[0]); | 
|   43  |   45  | 
|   44     for (i = 0; i <= order; ++i) |   46     for (i = 0; i <= order; ++i) | 
|   45     { |   47     { | 
|   46         temp1W32 = R[i] * (1 << norm); |   48         temp1W32 = R[i] * (1 << norm); | 
 |   49         // UBSan: 12 * 268435456 cannot be represented in type 'int' | 
 |   50  | 
|   47         // Put R in hi and low format |   51         // Put R in hi and low format | 
|   48         R_hi[i] = (int16_t)(temp1W32 >> 16); |   52         R_hi[i] = (int16_t)(temp1W32 >> 16); | 
|   49         R_low[i] = (int16_t)((temp1W32 - ((int32_t)R_hi[i] * 65536)) >> 1); |   53         R_low[i] = (int16_t)((temp1W32 - ((int32_t)R_hi[i] * 65536)) >> 1); | 
|   50     } |   54     } | 
|   51  |   55  | 
|   52     // K = A[1] = -R[1] / R[0] |   56     // K = A[1] = -R[1] / R[0] | 
|   53  |   57  | 
|   54     temp2W32 = R[1] * (1 << norm); // R[1] in Q31 |   58     temp2W32 = R[1] * (1 << norm); // R[1] in Q31 | 
|   55     temp3W32 = WEBRTC_SPL_ABS_W32(temp2W32); // abs R[1] |   59     temp3W32 = WEBRTC_SPL_ABS_W32(temp2W32); // abs R[1] | 
|   56     temp1W32 = WebRtcSpl_DivW32HiLow(temp3W32, R_hi[0], R_low[0]); // abs(R[1])/
     R[0] in Q31 |   60     temp1W32 = WebRtcSpl_DivW32HiLow(temp3W32, R_hi[0], R_low[0]); // abs(R[1])/
     R[0] in Q31 | 
| (...skipping 179 matching lines...) Expand 10 before | Expand all | Expand 10 after  Loading... | 
|  236     for (i = 1; i <= order; i++) |  240     for (i = 1; i <= order; i++) | 
|  237     { |  241     { | 
|  238         // temp1W32 in Q27 |  242         // temp1W32 in Q27 | 
|  239         temp1W32 = (int32_t)A_hi[i] * 65536 |  243         temp1W32 = (int32_t)A_hi[i] * 65536 | 
|  240                 + WEBRTC_SPL_LSHIFT_W32((int32_t)A_low[i], 1); |  244                 + WEBRTC_SPL_LSHIFT_W32((int32_t)A_low[i], 1); | 
|  241         // Round and store upper word |  245         // Round and store upper word | 
|  242         A[i] = (int16_t)(((temp1W32 * 2) + 32768) >> 16); |  246         A[i] = (int16_t)(((temp1W32 * 2) + 32768) >> 16); | 
|  243     } |  247     } | 
|  244     return 1; // Stable filters |  248     return 1; // Stable filters | 
|  245 } |  249 } | 
| OLD | NEW |