73 lines
2.4 KiB
C++
73 lines
2.4 KiB
C++
|
/*
|
||
|
|
* Copyright (c) 2017 The WebRTC project authors. All Rights Reserved.
|
||
|
|
*
|
||
|
|
* Use of this source code is governed by a BSD-style license
|
||
|
|
* that can be found in the LICENSE file in the root of the source
|
||
|
|
* tree. An additional intellectual property rights grant can be found
|
||
|
|
* in the file PATENTS. All contributing project authors may
|
||
|
|
* be found in the AUTHORS file in the root of the source tree.
|
||
|
|
*/
|
||
|
|
|
||
|
|
#include "webrtc/modules/audio_processing/aec3/fft_buffer.h"
|
||
|
|
|
||
|
|
#include <algorithm>
|
||
|
|
|
||
|
|
#include "webrtc/base/checks.h"
|
||
|
|
#include "webrtc/modules/audio_processing/aec3/aec3_common.h"
|
||
|
|
|
||
|
|
namespace webrtc {
|
||
|
|
|
||
|
|
FftBuffer::FftBuffer(Aec3Optimization optimization,
|
||
|
|
size_t num_partitions,
|
||
|
|
const std::vector<size_t> num_ffts_for_spectral_sums)
|
||
|
|
: optimization_(optimization),
|
||
|
|
fft_buffer_(num_partitions),
|
||
|
|
spectrum_buffer_(num_partitions, std::array<float, kFftLengthBy2Plus1>()),
|
||
|
|
spectral_sums_(num_ffts_for_spectral_sums.size(),
|
||
|
|
std::array<float, kFftLengthBy2Plus1>()) {
|
||
|
|
// Current implementation only allows a maximum of one spectral sum lengths.
|
||
|
|
RTC_DCHECK_EQ(1, num_ffts_for_spectral_sums.size());
|
||
|
|
spectral_sums_length_ = num_ffts_for_spectral_sums[0];
|
||
|
|
RTC_DCHECK_GE(fft_buffer_.size(), spectral_sums_length_);
|
||
|
|
|
||
|
|
for (auto& sum : spectral_sums_) {
|
||
|
|
sum.fill(0.f);
|
||
|
|
}
|
||
|
|
|
||
|
|
for (auto& spectrum : spectrum_buffer_) {
|
||
|
|
spectrum.fill(0.f);
|
||
|
|
}
|
||
|
|
|
||
|
|
for (auto& fft : fft_buffer_) {
|
||
|
|
fft.Clear();
|
||
|
|
}
|
||
|
|
}
|
||
|
|
|
||
|
|
FftBuffer::~FftBuffer() = default;
|
||
|
|
|
||
|
|
void FftBuffer::Insert(const FftData& fft) {
|
||
|
|
// Insert the fft into the buffer.
|
||
|
|
position_ = (position_ - 1 + fft_buffer_.size()) % fft_buffer_.size();
|
||
|
|
fft_buffer_[position_].Assign(fft);
|
||
|
|
|
||
|
|
// Compute and insert the spectrum for the FFT into the spectrum buffer.
|
||
|
|
fft.Spectrum(optimization_, &spectrum_buffer_[position_]);
|
||
|
|
|
||
|
|
// Pre-compute and cachec the spectral sums.
|
||
|
|
std::copy(spectrum_buffer_[position_].begin(),
|
||
|
|
spectrum_buffer_[position_].end(), spectral_sums_[0].begin());
|
||
|
|
size_t position = (position_ + 1) % fft_buffer_.size();
|
||
|
|
for (size_t j = 1; j < spectral_sums_length_; ++j) {
|
||
|
|
const std::array<float, kFftLengthBy2Plus1>& spectrum =
|
||
|
|
spectrum_buffer_[position];
|
||
|
|
|
||
|
|
for (size_t k = 0; k < spectral_sums_[0].size(); ++k) {
|
||
|
|
spectral_sums_[0][k] += spectrum[k];
|
||
|
|
}
|
||
|
|
|
||
|
|
position = position < (fft_buffer_.size() - 1) ? position + 1 : 0;
|
||
|
|
}
|
||
|
|
}
|
||
|
|
|
||
|
|
} // namespace webrtc
|