peah | 522d71b | 2017-02-23 05:16:26 -0800 | [diff] [blame] | 1 | /* |
| 2 | * Copyright (c) 2017 The WebRTC project authors. All Rights Reserved. |
| 3 | * |
| 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 |
| 6 | * tree. An additional intellectual property rights grant can be found |
| 7 | * in the file PATENTS. All contributing project authors may |
| 8 | * be found in the AUTHORS file in the root of the source tree. |
| 9 | */ |
| 10 | |
Mirko Bonadei | 92ea95e | 2017-09-15 06:47:31 +0200 | [diff] [blame] | 11 | #ifndef MODULES_AUDIO_PROCESSING_AEC3_AEC3_FFT_H_ |
| 12 | #define MODULES_AUDIO_PROCESSING_AEC3_AEC3_FFT_H_ |
peah | 522d71b | 2017-02-23 05:16:26 -0800 | [diff] [blame] | 13 | |
| 14 | #include <array> |
| 15 | |
Mirko Bonadei | 92ea95e | 2017-09-15 06:47:31 +0200 | [diff] [blame] | 16 | #include "api/array_view.h" |
Mirko Bonadei | f0d64a5 | 2020-04-17 12:25:19 +0200 | [diff] [blame] | 17 | #include "common_audio/third_party/ooura/fft_size_128/ooura_fft.h" |
Mirko Bonadei | 92ea95e | 2017-09-15 06:47:31 +0200 | [diff] [blame] | 18 | #include "modules/audio_processing/aec3/aec3_common.h" |
| 19 | #include "modules/audio_processing/aec3/fft_data.h" |
Yves Gerey | 988cc08 | 2018-10-23 12:03:01 +0200 | [diff] [blame] | 20 | #include "rtc_base/checks.h" |
Steve Anton | 10542f2 | 2019-01-11 09:11:00 -0800 | [diff] [blame] | 21 | #include "rtc_base/constructor_magic.h" |
peah | 522d71b | 2017-02-23 05:16:26 -0800 | [diff] [blame] | 22 | |
| 23 | namespace webrtc { |
| 24 | |
| 25 | // Wrapper class that provides 128 point real valued FFT functionality with the |
| 26 | // FftData type. |
| 27 | class Aec3Fft { |
| 28 | public: |
Per Åhgren | 47d7fbd | 2018-04-24 12:44:29 +0200 | [diff] [blame] | 29 | enum class Window { kRectangular, kHanning, kSqrtHanning }; |
Per Åhgren | d20639f | 2018-01-11 10:29:49 +0100 | [diff] [blame] | 30 | |
Per Åhgren | 80b3012 | 2020-06-02 14:28:25 +0200 | [diff] [blame] | 31 | Aec3Fft(); |
| 32 | |
peah | 522d71b | 2017-02-23 05:16:26 -0800 | [diff] [blame] | 33 | // Computes the FFT. Note that both the input and output are modified. |
| 34 | void Fft(std::array<float, kFftLength>* x, FftData* X) const { |
| 35 | RTC_DCHECK(x); |
| 36 | RTC_DCHECK(X); |
| 37 | ooura_fft_.Fft(x->data()); |
| 38 | X->CopyFromPackedArray(*x); |
| 39 | } |
| 40 | // Computes the inverse Fft. |
| 41 | void Ifft(const FftData& X, std::array<float, kFftLength>* x) const { |
| 42 | RTC_DCHECK(x); |
| 43 | X.CopyToPackedArray(x); |
| 44 | ooura_fft_.InverseFft(x->data()); |
| 45 | } |
| 46 | |
Per Åhgren | d20639f | 2018-01-11 10:29:49 +0100 | [diff] [blame] | 47 | // Windows the input using a Hanning window, and then adds padding of |
| 48 | // kFftLengthBy2 initial zeros before computing the Fft. |
| 49 | void ZeroPaddedFft(rtc::ArrayView<const float> x, |
| 50 | Window window, |
| 51 | FftData* X) const; |
peah | 522d71b | 2017-02-23 05:16:26 -0800 | [diff] [blame] | 52 | |
| 53 | // Concatenates the kFftLengthBy2 values long x and x_old before computing the |
| 54 | // Fft. After that, x is copied to x_old. |
| 55 | void PaddedFft(rtc::ArrayView<const float> x, |
Per Åhgren | 658ad88 | 2018-04-27 14:22:37 +0200 | [diff] [blame] | 56 | rtc::ArrayView<const float> x_old, |
Per Åhgren | 47d7fbd | 2018-04-24 12:44:29 +0200 | [diff] [blame] | 57 | FftData* X) const { |
| 58 | PaddedFft(x, x_old, Window::kRectangular, X); |
| 59 | } |
| 60 | |
| 61 | // Padded Fft using a time-domain window. |
| 62 | void PaddedFft(rtc::ArrayView<const float> x, |
Per Åhgren | 658ad88 | 2018-04-27 14:22:37 +0200 | [diff] [blame] | 63 | rtc::ArrayView<const float> x_old, |
Per Åhgren | 47d7fbd | 2018-04-24 12:44:29 +0200 | [diff] [blame] | 64 | Window window, |
peah | 522d71b | 2017-02-23 05:16:26 -0800 | [diff] [blame] | 65 | FftData* X) const; |
| 66 | |
| 67 | private: |
| 68 | const OouraFft ooura_fft_; |
| 69 | |
| 70 | RTC_DISALLOW_COPY_AND_ASSIGN(Aec3Fft); |
| 71 | }; |
| 72 | |
| 73 | } // namespace webrtc |
| 74 | |
Mirko Bonadei | 92ea95e | 2017-09-15 06:47:31 +0200 | [diff] [blame] | 75 | #endif // MODULES_AUDIO_PROCESSING_AEC3_AEC3_FFT_H_ |