OLD | NEW |
1 /* | 1 /* |
2 * Copyright (c) 2015 The WebRTC project authors. All Rights Reserved. | 2 * Copyright (c) 2015 The WebRTC project authors. All Rights Reserved. |
3 * | 3 * |
4 * Use of this source code is governed by a BSD-style license | 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 | 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 | 6 * tree. An additional intellectual property rights grant can be found |
7 * in the file PATENTS. All contributing project authors may | 7 * in the file PATENTS. All contributing project authors may |
8 * be found in the AUTHORS file in the root of the source tree. | 8 * be found in the AUTHORS file in the root of the source tree. |
9 */ | 9 */ |
10 | 10 |
(...skipping 24 matching lines...) Expand all Loading... |
35 } // namespace | 35 } // namespace |
36 | 36 |
37 RealFourierOoura::RealFourierOoura(int fft_order) | 37 RealFourierOoura::RealFourierOoura(int fft_order) |
38 : order_(fft_order), | 38 : order_(fft_order), |
39 length_(FftLength(order_)), | 39 length_(FftLength(order_)), |
40 complex_length_(ComplexLength(order_)), | 40 complex_length_(ComplexLength(order_)), |
41 // Zero-initializing work_ip_ will cause rdft to initialize these work | 41 // Zero-initializing work_ip_ will cause rdft to initialize these work |
42 // arrays on the first call. | 42 // arrays on the first call. |
43 work_ip_(new size_t[ComputeWorkIpSize(length_)]()), | 43 work_ip_(new size_t[ComputeWorkIpSize(length_)]()), |
44 work_w_(new float[complex_length_]()) { | 44 work_w_(new float[complex_length_]()) { |
45 CHECK_GE(fft_order, 1); | 45 RTC_CHECK_GE(fft_order, 1); |
46 } | 46 } |
47 | 47 |
48 void RealFourierOoura::Forward(const float* src, complex<float>* dest) const { | 48 void RealFourierOoura::Forward(const float* src, complex<float>* dest) const { |
49 { | 49 { |
50 // This cast is well-defined since C++11. See "Non-static data members" at: | 50 // This cast is well-defined since C++11. See "Non-static data members" at: |
51 // http://en.cppreference.com/w/cpp/numeric/complex | 51 // http://en.cppreference.com/w/cpp/numeric/complex |
52 auto dest_float = reinterpret_cast<float*>(dest); | 52 auto dest_float = reinterpret_cast<float*>(dest); |
53 std::copy(src, src + length_, dest_float); | 53 std::copy(src, src + length_, dest_float); |
54 WebRtc_rdft(length_, 1, dest_float, work_ip_.get(), work_w_.get()); | 54 WebRtc_rdft(length_, 1, dest_float, work_ip_.get(), work_w_.get()); |
55 } | 55 } |
(...skipping 20 matching lines...) Expand all Loading... |
76 } | 76 } |
77 | 77 |
78 WebRtc_rdft(length_, -1, dest, work_ip_.get(), work_w_.get()); | 78 WebRtc_rdft(length_, -1, dest, work_ip_.get(), work_w_.get()); |
79 | 79 |
80 // Ooura returns a scaled version. | 80 // Ooura returns a scaled version. |
81 const float scale = 2.0f / length_; | 81 const float scale = 2.0f / length_; |
82 std::for_each(dest, dest + length_, [scale](float& v) { v *= scale; }); | 82 std::for_each(dest, dest + length_, [scale](float& v) { v *= scale; }); |
83 } | 83 } |
84 | 84 |
85 } // namespace webrtc | 85 } // namespace webrtc |
OLD | NEW |