Finalization of the first version of EchoCanceller 3

webrtc/modules/audio_processing/aec3/erle_estimator.cc

Index: webrtc/modules/audio_processing/aec3/erle_estimator.cc
diff --git a/webrtc/modules/audio_processing/aec3/erle_estimator.cc b/webrtc/modules/audio_processing/aec3/erle_estimator.cc
new file mode 100644
index 0000000000000000000000000000000000000000..58d3b123d4403bf3f40781f2bff54f2141878ed6
--- /dev/null
+++ b/webrtc/modules/audio_processing/aec3/erle_estimator.cc
@@ -0,0 +1,63 @@
+/*
+ *  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/erle_estimator.h"
+
+#include <algorithm>
+
+namespace webrtc {
+
+namespace {
+constexpr float kMinErle = 1.f;
+constexpr float kMaxErle = 8.f;
+
+}  // namespace
+
+void ErleEstimator::Update(
+    const std::array<float, kFftLengthBy2Plus1>& render_spectrum,
+    const std::array<float, kFftLengthBy2Plus1>& capture_spectrum,
+    const std::array<float, kFftLengthBy2Plus1>& linear_aec_output_spectrum) {
+  const std::array<float, kFftLengthBy2Plus1>& X2 = render_spectrum;

[ivoc, 2017/02/10 13:52:35]

I would prefer const auto&.

[peah-webrtc, 2017/02/20 07:37:16]

Done.

+  const std::array<float, kFftLengthBy2Plus1>& Y2 = capture_spectrum;
+  const std::array<float, kFftLengthBy2Plus1>& E2 = linear_aec_output_spectrum;
+
+  const float kX2Min = 44015068.0f;

[ivoc, 2017/02/10 13:52:34]

Where does this value come from? A comment would be helpful. Also constexpr.

[peah-webrtc, 2017/02/20 07:37:17]

Done.

+
+  // Update the estimates in a clamped minimum statistics manner.
+  for (size_t k = 1; k < kFftLengthBy2; ++k) {
+    if (X2[k] > kX2Min && E2[k] > 0.f) {
+      const float new_erle = Y2[k] / E2[k];
+      if (new_erle > erle_[k]) {
+        hold_counters_[k - 1] = 150;
+        erle_[k] += 0.1f * (new_erle - erle_[k]);
+        erle_[k] = std::max(kMinErle, std::min(erle_[k], kMaxErle));
+      }
+    }
+  }
+
+  std::for_each(hold_counters_.begin(), hold_counters_.end(),
+                [](int& a) { --a; });
+  std::transform(hold_counters_.begin(), hold_counters_.end(),
+                 erle_.begin() + 1, erle_.begin() + 1, [](int a, float b) {
+                   return a > 0 ? b : std::max(kMinErle, 0.99f * b);
+                 });
+
+  erle_[0] = erle_[1];
+  erle_[kFftLengthBy2] = erle_[kFftLengthBy2 - 1];
+}
+
+ErleEstimator::ErleEstimator() {

[ivoc, 2017/02/10 13:52:34]

Move the constructor and destructor to the top of the file please.

[peah-webrtc, 2017/02/20 07:37:16]

Done.

+  erle_.fill(kMinErle);
+  hold_counters_.fill(0);
+}
+
+ErleEstimator::~ErleEstimator() = default;
+
+}  // namespace webrtc