ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 1 | /* |
| 2 | * Copyright (c) 2014 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 | |
| 11 | #ifndef WEBRTC_MODULES_AUDIO_PROCESSING_INTELLIGIBILITY_INTELLIGIBILITY_UTILS_H_ |
| 12 | #define WEBRTC_MODULES_AUDIO_PROCESSING_INTELLIGIBILITY_INTELLIGIBILITY_UTILS_H_ |
| 13 | |
| 14 | #include <complex> |
kwiberg | 88788ad | 2016-02-19 07:04:49 -0800 | [diff] [blame] | 15 | #include <memory> |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 16 | |
| 17 | namespace webrtc { |
| 18 | |
| 19 | namespace intelligibility { |
| 20 | |
Alejandro Luebs | 3234819 | 2016-02-17 20:04:19 -0800 | [diff] [blame] | 21 | // Internal helper for computing the power of a stream of arrays. |
| 22 | // The result is an array of power per position: the i-th power is the power of |
| 23 | // the stream of data on the i-th positions in the input arrays. |
| 24 | class PowerEstimator { |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 25 | public: |
Alejandro Luebs | 3234819 | 2016-02-17 20:04:19 -0800 | [diff] [blame] | 26 | // Construct an instance for the given input array length (|freqs|), with the |
| 27 | // appropriate parameters. |decay| is the forgetting factor. |
| 28 | PowerEstimator(size_t freqs, float decay); |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 29 | |
Alejandro Luebs | 3234819 | 2016-02-17 20:04:19 -0800 | [diff] [blame] | 30 | // Add a new data point to the series. |
| 31 | void Step(const std::complex<float>* data); |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 32 | |
Alejandro Luebs | 3234819 | 2016-02-17 20:04:19 -0800 | [diff] [blame] | 33 | // The current power array. |
| 34 | const float* Power(); |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 35 | |
| 36 | private: |
ekm | db4fecf | 2015-06-22 17:49:08 -0700 | [diff] [blame] | 37 | // TODO(ekmeyerson): Switch the following running means |
kwiberg | 88788ad | 2016-02-19 07:04:49 -0800 | [diff] [blame] | 38 | // and histories from std::unique_ptr to std::vector. |
| 39 | std::unique_ptr<std::complex<float>[]> running_mean_sq_; |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 40 | |
Alejandro Luebs | 3234819 | 2016-02-17 20:04:19 -0800 | [diff] [blame] | 41 | // The current magnitude array. |
kwiberg | 88788ad | 2016-02-19 07:04:49 -0800 | [diff] [blame] | 42 | std::unique_ptr<float[]> magnitude_; |
Alejandro Luebs | 3234819 | 2016-02-17 20:04:19 -0800 | [diff] [blame] | 43 | // The current power array. |
kwiberg | 88788ad | 2016-02-19 07:04:49 -0800 | [diff] [blame] | 44 | std::unique_ptr<float[]> power_; |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 45 | |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 46 | const size_t num_freqs_; |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 47 | const float decay_; |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 48 | }; |
| 49 | |
Alejandro Luebs | 3234819 | 2016-02-17 20:04:19 -0800 | [diff] [blame] | 50 | // Helper class for smoothing gain changes. On each application step, the |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 51 | // currently used gains are changed towards a set of settable target gains, |
| 52 | // constrained by a limit on the magnitude of the changes. |
| 53 | class GainApplier { |
| 54 | public: |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 55 | GainApplier(size_t freqs, float change_limit); |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 56 | |
| 57 | // Copy |in_block| to |out_block|, multiplied by the current set of gains, |
| 58 | // and step the current set of gains towards the target set. |
| 59 | void Apply(const std::complex<float>* in_block, |
| 60 | std::complex<float>* out_block); |
| 61 | |
| 62 | // Return the current target gain set. Modify this array to set the targets. |
ekm | db4fecf | 2015-06-22 17:49:08 -0700 | [diff] [blame] | 63 | float* target() const { return target_.get(); } |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 64 | |
| 65 | private: |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 66 | const size_t num_freqs_; |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 67 | const float change_limit_; |
kwiberg | 88788ad | 2016-02-19 07:04:49 -0800 | [diff] [blame] | 68 | std::unique_ptr<float[]> target_; |
| 69 | std::unique_ptr<float[]> current_; |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 70 | }; |
| 71 | |
| 72 | } // namespace intelligibility |
| 73 | |
| 74 | } // namespace webrtc |
| 75 | |
| 76 | #endif // WEBRTC_MODULES_AUDIO_PROCESSING_INTELLIGIBILITY_INTELLIGIBILITY_UTILS_H_ |