OLD | NEW |
1 /* | 1 /* |
2 * Copyright (c) 2012 The WebRTC project authors. All Rights Reserved. | 2 * Copyright (c) 2012 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 #include "webrtc/common_audio/signal_processing/include/signal_processing_librar
y.h" | 11 #include "webrtc/common_audio/signal_processing/dot_product_with_scale.h" |
| 12 |
| 13 #include "webrtc/base/safe_conversions.h" |
12 | 14 |
13 int32_t WebRtcSpl_DotProductWithScale(const int16_t* vector1, | 15 int32_t WebRtcSpl_DotProductWithScale(const int16_t* vector1, |
14 const int16_t* vector2, | 16 const int16_t* vector2, |
15 size_t length, | 17 size_t length, |
16 int scaling) { | 18 int scaling) { |
17 int32_t sum = 0; | 19 int64_t sum = 0; |
18 size_t i = 0; | 20 size_t i = 0; |
19 | 21 |
20 /* Unroll the loop to improve performance. */ | 22 /* Unroll the loop to improve performance. */ |
21 for (i = 0; i + 3 < length; i += 4) { | 23 for (i = 0; i + 3 < length; i += 4) { |
22 sum += (vector1[i + 0] * vector2[i + 0]) >> scaling; | 24 sum += (vector1[i + 0] * vector2[i + 0]) >> scaling; |
23 sum += (vector1[i + 1] * vector2[i + 1]) >> scaling; | 25 sum += (vector1[i + 1] * vector2[i + 1]) >> scaling; |
24 sum += (vector1[i + 2] * vector2[i + 2]) >> scaling; | 26 sum += (vector1[i + 2] * vector2[i + 2]) >> scaling; |
25 sum += (vector1[i + 3] * vector2[i + 3]) >> scaling; | 27 sum += (vector1[i + 3] * vector2[i + 3]) >> scaling; |
26 } | 28 } |
27 for (; i < length; i++) { | 29 for (; i < length; i++) { |
28 sum += (vector1[i] * vector2[i]) >> scaling; | 30 sum += (vector1[i] * vector2[i]) >> scaling; |
29 } | 31 } |
30 | 32 |
31 return sum; | 33 return rtc::saturated_cast<int32_t>(sum); |
32 } | 34 } |
OLD | NEW |