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