OLD | NEW |
1 /* | 1 /* |
2 * Copyright (c) 2014 The WebRTC project authors. All Rights Reserved. | 2 * Copyright (c) 2014 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 #define _USE_MATH_DEFINES | 11 #define _USE_MATH_DEFINES |
12 | 12 |
13 #include "webrtc/common_audio/window_generator.h" | 13 #include "webrtc/common_audio/window_generator.h" |
14 | 14 |
15 #include <cmath> | 15 #include <cmath> |
16 #include <complex> | 16 #include <complex> |
17 | 17 |
18 #include "webrtc/base/checks.h" | 18 #include "webrtc/rtc_base/checks.h" |
19 | 19 |
20 using std::complex; | 20 using std::complex; |
21 | 21 |
22 namespace { | 22 namespace { |
23 | 23 |
24 // Modified Bessel function of order 0 for complex inputs. | 24 // Modified Bessel function of order 0 for complex inputs. |
25 complex<float> I0(complex<float> x) { | 25 complex<float> I0(complex<float> x) { |
26 complex<float> y = x / 3.75f; | 26 complex<float> y = x / 3.75f; |
27 y *= y; | 27 y *= y; |
28 return 1.0f + y * ( | 28 return 1.0f + y * ( |
(...skipping 34 matching lines...) Expand 10 before | Expand all | Expand 10 after Loading... |
63 window[length - i - 1] = sqrtf(window[length - i - 1] / sum); | 63 window[length - i - 1] = sqrtf(window[length - i - 1] / sum); |
64 window[i] = window[length - i - 1]; | 64 window[i] = window[length - i - 1]; |
65 } | 65 } |
66 if (length % 2 == 1) { | 66 if (length % 2 == 1) { |
67 window[half - 1] = sqrtf(window[half - 1] / sum); | 67 window[half - 1] = sqrtf(window[half - 1] / sum); |
68 } | 68 } |
69 } | 69 } |
70 | 70 |
71 } // namespace webrtc | 71 } // namespace webrtc |
72 | 72 |
OLD | NEW |