| 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 |