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 |
(...skipping 20 matching lines...) Expand all Loading... |
31 1.2067492f + y * ( | 31 1.2067492f + y * ( |
32 0.2659732f + y * ( | 32 0.2659732f + y * ( |
33 0.360768e-1f + y * 0.45813e-2f))))); | 33 0.360768e-1f + y * 0.45813e-2f))))); |
34 } | 34 } |
35 | 35 |
36 } // namespace | 36 } // namespace |
37 | 37 |
38 namespace webrtc { | 38 namespace webrtc { |
39 | 39 |
40 void WindowGenerator::Hanning(int length, float* window) { | 40 void WindowGenerator::Hanning(int length, float* window) { |
41 CHECK_GT(length, 1); | 41 RTC_CHECK_GT(length, 1); |
42 CHECK(window != nullptr); | 42 RTC_CHECK(window != nullptr); |
43 for (int i = 0; i < length; ++i) { | 43 for (int i = 0; i < length; ++i) { |
44 window[i] = 0.5f * (1 - cosf(2 * static_cast<float>(M_PI) * i / | 44 window[i] = 0.5f * (1 - cosf(2 * static_cast<float>(M_PI) * i / |
45 (length - 1))); | 45 (length - 1))); |
46 } | 46 } |
47 } | 47 } |
48 | 48 |
49 void WindowGenerator::KaiserBesselDerived(float alpha, size_t length, | 49 void WindowGenerator::KaiserBesselDerived(float alpha, size_t length, |
50 float* window) { | 50 float* window) { |
51 CHECK_GT(length, 1U); | 51 RTC_CHECK_GT(length, 1U); |
52 CHECK(window != nullptr); | 52 RTC_CHECK(window != nullptr); |
53 | 53 |
54 const size_t half = (length + 1) / 2; | 54 const size_t half = (length + 1) / 2; |
55 float sum = 0.0f; | 55 float sum = 0.0f; |
56 | 56 |
57 for (size_t i = 0; i <= half; ++i) { | 57 for (size_t i = 0; i <= half; ++i) { |
58 complex<float> r = (4.0f * i) / length - 1.0f; | 58 complex<float> r = (4.0f * i) / length - 1.0f; |
59 sum += I0(static_cast<float>(M_PI) * alpha * sqrt(1.0f - r * r)).real(); | 59 sum += I0(static_cast<float>(M_PI) * alpha * sqrt(1.0f - r * r)).real(); |
60 window[i] = sum; | 60 window[i] = sum; |
61 } | 61 } |
62 for (size_t i = length - 1; i >= half; --i) { | 62 for (size_t i = length - 1; i >= half; --i) { |
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 |