0cbb58e046
This is a reland of 87a7b82520b83a6cf42da27cdc46142c2eb6248c Original change's description: > Refactoring of the noise suppressor and adding true multichannel support > > This CL adds proper multichannel support to the noise suppressor. > To accomplish that in a safe way, a full refactoring of the noise > suppressor code has been done. > > Due to floating point precision, the changes made are not entirely > bitexact. They are, however, very close to being bitexact. > > As a safety measure, the former noise suppressor code is preserved > and a kill-switch is added to allow revering to the legacy noise > suppressor in case issues arise. > > Bug: webrtc:10895, b/143344262 > Change-Id: I0b071011b23265ac12e6d4b3956499d122286657 > Reviewed-on: https://webrtc-review.googlesource.com/c/src/+/158407 > Commit-Queue: Per Åhgren <peah@webrtc.org> > Reviewed-by: Gustaf Ullberg <gustaf@webrtc.org> > Cr-Commit-Position: refs/heads/master@{#29646} Bug: webrtc:10895, b/143344262 Change-Id: I236f1e67bb0baa4e30908a4cf7a8a7bb55fbced3 Reviewed-on: https://webrtc-review.googlesource.com/c/src/+/158747 Reviewed-by: Gustaf Ullberg <gustaf@webrtc.org> Commit-Queue: Per Åhgren <peah@webrtc.org> Cr-Commit-Position: refs/heads/master@{#29663}
65 lines
2.0 KiB
C++
65 lines
2.0 KiB
C++
/*
|
|
* Copyright (c) 2019 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 "modules/audio_processing/ns/ns_fft.h"
|
|
|
|
#include "common_audio/third_party/fft4g/fft4g.h"
|
|
|
|
namespace webrtc {
|
|
|
|
NrFft::NrFft() : bit_reversal_state_(kFftSize / 2), tables_(kFftSize / 2) {
|
|
// Initialize WebRtc_rdt (setting (bit_reversal_state_[0] to 0 triggers
|
|
// initialization)
|
|
bit_reversal_state_[0] = 0.f;
|
|
std::array<float, kFftSize> tmp_buffer;
|
|
tmp_buffer.fill(0.f);
|
|
WebRtc_rdft(kFftSize, 1, tmp_buffer.data(), bit_reversal_state_.data(),
|
|
tables_.data());
|
|
}
|
|
|
|
void NrFft::Fft(rtc::ArrayView<float, kFftSize> time_data,
|
|
rtc::ArrayView<float, kFftSize> real,
|
|
rtc::ArrayView<float, kFftSize> imag) {
|
|
WebRtc_rdft(kFftSize, 1, time_data.data(), bit_reversal_state_.data(),
|
|
tables_.data());
|
|
|
|
imag[0] = 0;
|
|
real[0] = time_data[0];
|
|
|
|
imag[kFftSizeBy2Plus1 - 1] = 0;
|
|
real[kFftSizeBy2Plus1 - 1] = time_data[1];
|
|
|
|
for (size_t i = 1; i < kFftSizeBy2Plus1 - 1; ++i) {
|
|
real[i] = time_data[2 * i];
|
|
imag[i] = time_data[2 * i + 1];
|
|
}
|
|
}
|
|
|
|
void NrFft::Ifft(rtc::ArrayView<const float> real,
|
|
rtc::ArrayView<const float> imag,
|
|
rtc::ArrayView<float> time_data) {
|
|
time_data[0] = real[0];
|
|
time_data[1] = real[kFftSizeBy2Plus1 - 1];
|
|
for (size_t i = 1; i < kFftSizeBy2Plus1 - 1; ++i) {
|
|
time_data[2 * i] = real[i];
|
|
time_data[2 * i + 1] = imag[i];
|
|
}
|
|
WebRtc_rdft(kFftSize, -1, time_data.data(), bit_reversal_state_.data(),
|
|
tables_.data());
|
|
|
|
// Scale the output
|
|
constexpr float kScaling = 2.f / kFftSize;
|
|
for (float& d : time_data) {
|
|
d *= kScaling;
|
|
}
|
|
}
|
|
|
|
} // namespace webrtc
|