OLD | NEW |
1 // Copyright (c) 2012 The Chromium Authors. All rights reserved. | 1 // Copyright (c) 2012 The Chromium Authors. All rights reserved. |
2 // Use of this source code is governed by a BSD-style license that can be | 2 // Use of this source code is governed by a BSD-style license that can be |
3 // found in the LICENSE file. | 3 // found in the LICENSE file. |
4 | 4 |
5 // MSVC++ requires this to be set before any other includes to get M_PI. | 5 // MSVC++ requires this to be set before any other includes to get M_PI. |
6 #define _USE_MATH_DEFINES | 6 #define _USE_MATH_DEFINES |
7 | 7 |
8 #include <cmath> | 8 #include <cmath> |
9 | 9 |
10 #include "base/bind.h" | 10 #include "base/bind.h" |
(...skipping 15 matching lines...) Expand all Loading... |
26 | 26 |
27 static const double kSampleRateRatio = 192000.0 / 44100.0; | 27 static const double kSampleRateRatio = 192000.0 / 44100.0; |
28 static const double kKernelInterpolationFactor = 0.5; | 28 static const double kKernelInterpolationFactor = 0.5; |
29 | 29 |
30 // Command line switch for runtime adjustment of ConvolveBenchmark iterations. | 30 // Command line switch for runtime adjustment of ConvolveBenchmark iterations. |
31 static const char kConvolveIterations[] = "convolve-iterations"; | 31 static const char kConvolveIterations[] = "convolve-iterations"; |
32 | 32 |
33 // Helper class to ensure ChunkedResample() functions properly. | 33 // Helper class to ensure ChunkedResample() functions properly. |
34 class MockSource { | 34 class MockSource { |
35 public: | 35 public: |
36 MOCK_METHOD2(ProvideInput, void(float* destination, int frames)); | 36 MOCK_METHOD2(ProvideInput, void(int frames, float* destination)); |
37 }; | 37 }; |
38 | 38 |
39 ACTION(ClearBuffer) { | 39 ACTION(ClearBuffer) { |
40 memset(arg0, 0, arg1 * sizeof(float)); | 40 memset(arg1, 0, arg0 * sizeof(float)); |
41 } | 41 } |
42 | 42 |
43 ACTION(FillBuffer) { | 43 ACTION(FillBuffer) { |
44 // Value chosen arbitrarily such that SincResampler resamples it to something | 44 // Value chosen arbitrarily such that SincResampler resamples it to something |
45 // easily representable on all platforms; e.g., using kSampleRateRatio this | 45 // easily representable on all platforms; e.g., using kSampleRateRatio this |
46 // becomes 1.81219. | 46 // becomes 1.81219. |
47 memset(arg0, 64, arg1 * sizeof(float)); | 47 memset(arg1, 64, arg0 * sizeof(float)); |
48 } | 48 } |
49 | 49 |
50 // Test requesting multiples of ChunkSize() frames results in the proper number | 50 // Test requesting multiples of ChunkSize() frames results in the proper number |
51 // of callbacks. | 51 // of callbacks. |
52 TEST(SincResamplerTest, ChunkedResample) { | 52 TEST(SincResamplerTest, ChunkedResample) { |
53 MockSource mock_source; | 53 MockSource mock_source; |
54 | 54 |
55 // Choose a high ratio of input to output samples which will result in quick | 55 // Choose a high ratio of input to output samples which will result in quick |
56 // exhaustion of SincResampler's internal buffers. | 56 // exhaustion of SincResampler's internal buffers. |
57 SincResampler resampler( | 57 SincResampler resampler( |
58 kSampleRateRatio, | 58 kSampleRateRatio, SincResampler::kDefaultRequestSize, |
59 base::Bind(&MockSource::ProvideInput, base::Unretained(&mock_source))); | 59 base::Bind(&MockSource::ProvideInput, base::Unretained(&mock_source))); |
60 | 60 |
61 static const int kChunks = 2; | 61 static const int kChunks = 2; |
62 int max_chunk_size = resampler.ChunkSize() * kChunks; | 62 int max_chunk_size = resampler.ChunkSize() * kChunks; |
63 scoped_ptr<float[]> resampled_destination(new float[max_chunk_size]); | 63 scoped_ptr<float[]> resampled_destination(new float[max_chunk_size]); |
64 | 64 |
65 // Verify requesting ChunkSize() frames causes a single callback. | 65 // Verify requesting ChunkSize() frames causes a single callback. |
66 EXPECT_CALL(mock_source, ProvideInput(_, _)) | 66 EXPECT_CALL(mock_source, ProvideInput(_, _)) |
67 .Times(1).WillOnce(ClearBuffer()); | 67 .Times(1).WillOnce(ClearBuffer()); |
68 resampler.Resample(resampled_destination.get(), resampler.ChunkSize()); | 68 resampler.Resample(resampler.ChunkSize(), resampled_destination.get()); |
69 | 69 |
70 // Verify requesting kChunks * ChunkSize() frames causes kChunks callbacks. | 70 // Verify requesting kChunks * ChunkSize() frames causes kChunks callbacks. |
71 testing::Mock::VerifyAndClear(&mock_source); | 71 testing::Mock::VerifyAndClear(&mock_source); |
72 EXPECT_CALL(mock_source, ProvideInput(_, _)) | 72 EXPECT_CALL(mock_source, ProvideInput(_, _)) |
73 .Times(kChunks).WillRepeatedly(ClearBuffer()); | 73 .Times(kChunks).WillRepeatedly(ClearBuffer()); |
74 resampler.Resample(resampled_destination.get(), max_chunk_size); | 74 resampler.Resample(max_chunk_size, resampled_destination.get()); |
75 } | 75 } |
76 | 76 |
77 // Test flush resets the internal state properly. | 77 // Test flush resets the internal state properly. |
78 TEST(SincResamplerTest, Flush) { | 78 TEST(SincResamplerTest, Flush) { |
79 MockSource mock_source; | 79 MockSource mock_source; |
80 SincResampler resampler( | 80 SincResampler resampler( |
81 kSampleRateRatio, | 81 kSampleRateRatio, SincResampler::kDefaultRequestSize, |
82 base::Bind(&MockSource::ProvideInput, base::Unretained(&mock_source))); | 82 base::Bind(&MockSource::ProvideInput, base::Unretained(&mock_source))); |
83 scoped_ptr<float[]> resampled_destination(new float[resampler.ChunkSize()]); | 83 scoped_ptr<float[]> resampled_destination(new float[resampler.ChunkSize()]); |
84 | 84 |
85 // Fill the resampler with junk data. | 85 // Fill the resampler with junk data. |
86 EXPECT_CALL(mock_source, ProvideInput(_, _)) | 86 EXPECT_CALL(mock_source, ProvideInput(_, _)) |
87 .Times(1).WillOnce(FillBuffer()); | 87 .Times(1).WillOnce(FillBuffer()); |
88 resampler.Resample(resampled_destination.get(), resampler.ChunkSize() / 2); | 88 resampler.Resample(resampler.ChunkSize() / 2, resampled_destination.get()); |
89 ASSERT_NE(resampled_destination[0], 0); | 89 ASSERT_NE(resampled_destination[0], 0); |
90 | 90 |
91 // Flush and request more data, which should all be zeros now. | 91 // Flush and request more data, which should all be zeros now. |
92 resampler.Flush(); | 92 resampler.Flush(); |
93 testing::Mock::VerifyAndClear(&mock_source); | 93 testing::Mock::VerifyAndClear(&mock_source); |
94 EXPECT_CALL(mock_source, ProvideInput(_, _)) | 94 EXPECT_CALL(mock_source, ProvideInput(_, _)) |
95 .Times(1).WillOnce(ClearBuffer()); | 95 .Times(1).WillOnce(ClearBuffer()); |
96 resampler.Resample(resampled_destination.get(), resampler.ChunkSize() / 2); | 96 resampler.Resample(resampler.ChunkSize() / 2, resampled_destination.get()); |
97 for (int i = 0; i < resampler.ChunkSize() / 2; ++i) | 97 for (int i = 0; i < resampler.ChunkSize() / 2; ++i) |
98 ASSERT_FLOAT_EQ(resampled_destination[i], 0); | 98 ASSERT_FLOAT_EQ(resampled_destination[i], 0); |
99 } | 99 } |
100 | 100 |
101 // Test flush resets the internal state properly. | 101 // Test flush resets the internal state properly. |
102 TEST(SincResamplerTest, DISABLED_SetRatioBench) { | 102 TEST(SincResamplerTest, DISABLED_SetRatioBench) { |
103 MockSource mock_source; | 103 MockSource mock_source; |
104 SincResampler resampler( | 104 SincResampler resampler( |
105 kSampleRateRatio, | 105 kSampleRateRatio, SincResampler::kDefaultRequestSize, |
106 base::Bind(&MockSource::ProvideInput, base::Unretained(&mock_source))); | 106 base::Bind(&MockSource::ProvideInput, base::Unretained(&mock_source))); |
107 | 107 |
108 base::TimeTicks start = base::TimeTicks::HighResNow(); | 108 base::TimeTicks start = base::TimeTicks::HighResNow(); |
109 for (int i = 1; i < 10000; ++i) | 109 for (int i = 1; i < 10000; ++i) |
110 resampler.SetRatio(1.0 / i); | 110 resampler.SetRatio(1.0 / i); |
111 double total_time_c_ms = | 111 double total_time_c_ms = |
112 (base::TimeTicks::HighResNow() - start).InMillisecondsF(); | 112 (base::TimeTicks::HighResNow() - start).InMillisecondsF(); |
113 printf("SetRatio() took %.2fms.\n", total_time_c_ms); | 113 printf("SetRatio() took %.2fms.\n", total_time_c_ms); |
114 } | 114 } |
115 | 115 |
(...skipping 10 matching lines...) Expand all Loading... |
126 // will be tested by the parameterized SincResampler tests below. | 126 // will be tested by the parameterized SincResampler tests below. |
127 #if defined(CONVOLVE_FUNC) | 127 #if defined(CONVOLVE_FUNC) |
128 TEST(SincResamplerTest, Convolve) { | 128 TEST(SincResamplerTest, Convolve) { |
129 #if defined(ARCH_CPU_X86_FAMILY) | 129 #if defined(ARCH_CPU_X86_FAMILY) |
130 ASSERT_TRUE(base::CPU().has_sse()); | 130 ASSERT_TRUE(base::CPU().has_sse()); |
131 #endif | 131 #endif |
132 | 132 |
133 // Initialize a dummy resampler. | 133 // Initialize a dummy resampler. |
134 MockSource mock_source; | 134 MockSource mock_source; |
135 SincResampler resampler( | 135 SincResampler resampler( |
136 kSampleRateRatio, | 136 kSampleRateRatio, SincResampler::kDefaultRequestSize, |
137 base::Bind(&MockSource::ProvideInput, base::Unretained(&mock_source))); | 137 base::Bind(&MockSource::ProvideInput, base::Unretained(&mock_source))); |
138 | 138 |
139 // The optimized Convolve methods are slightly more precise than Convolve_C(), | 139 // The optimized Convolve methods are slightly more precise than Convolve_C(), |
140 // so comparison must be done using an epsilon. | 140 // so comparison must be done using an epsilon. |
141 static const double kEpsilon = 0.00000005; | 141 static const double kEpsilon = 0.00000005; |
142 | 142 |
143 // Use a kernel from SincResampler as input and kernel data, this has the | 143 // Use a kernel from SincResampler as input and kernel data, this has the |
144 // benefit of already being properly sized and aligned for Convolve_SSE(). | 144 // benefit of already being properly sized and aligned for Convolve_SSE(). |
145 double result = resampler.Convolve_C( | 145 double result = resampler.Convolve_C( |
146 resampler.kernel_storage_.get(), resampler.kernel_storage_.get(), | 146 resampler.kernel_storage_.get(), resampler.kernel_storage_.get(), |
(...skipping 14 matching lines...) Expand all Loading... |
161 } | 161 } |
162 #endif | 162 #endif |
163 | 163 |
164 // Benchmark for the various Convolve() methods. Make sure to build with | 164 // Benchmark for the various Convolve() methods. Make sure to build with |
165 // branding=Chrome so that DCHECKs are compiled out when benchmarking. Original | 165 // branding=Chrome so that DCHECKs are compiled out when benchmarking. Original |
166 // benchmarks were run with --convolve-iterations=50000000. | 166 // benchmarks were run with --convolve-iterations=50000000. |
167 TEST(SincResamplerTest, ConvolveBenchmark) { | 167 TEST(SincResamplerTest, ConvolveBenchmark) { |
168 // Initialize a dummy resampler. | 168 // Initialize a dummy resampler. |
169 MockSource mock_source; | 169 MockSource mock_source; |
170 SincResampler resampler( | 170 SincResampler resampler( |
171 kSampleRateRatio, | 171 kSampleRateRatio, SincResampler::kDefaultRequestSize, |
172 base::Bind(&MockSource::ProvideInput, base::Unretained(&mock_source))); | 172 base::Bind(&MockSource::ProvideInput, base::Unretained(&mock_source))); |
173 | 173 |
174 // Retrieve benchmark iterations from command line. | 174 // Retrieve benchmark iterations from command line. |
175 int convolve_iterations = 10; | 175 int convolve_iterations = 10; |
176 std::string iterations(CommandLine::ForCurrentProcess()->GetSwitchValueASCII( | 176 std::string iterations(CommandLine::ForCurrentProcess()->GetSwitchValueASCII( |
177 kConvolveIterations)); | 177 kConvolveIterations)); |
178 if (!iterations.empty()) | 178 if (!iterations.empty()) |
179 base::StringToInt(iterations, &convolve_iterations); | 179 base::StringToInt(iterations, &convolve_iterations); |
180 | 180 |
181 printf("Benchmarking %d iterations:\n", convolve_iterations); | 181 printf("Benchmarking %d iterations:\n", convolve_iterations); |
(...skipping 45 matching lines...) Expand 10 before | Expand all | Expand 10 after Loading... |
227 #endif | 227 #endif |
228 } | 228 } |
229 | 229 |
230 #undef CONVOLVE_FUNC | 230 #undef CONVOLVE_FUNC |
231 | 231 |
232 // Fake audio source for testing the resampler. Generates a sinusoidal linear | 232 // Fake audio source for testing the resampler. Generates a sinusoidal linear |
233 // chirp (http://en.wikipedia.org/wiki/Chirp) which can be tuned to stress the | 233 // chirp (http://en.wikipedia.org/wiki/Chirp) which can be tuned to stress the |
234 // resampler for the specific sample rate conversion being used. | 234 // resampler for the specific sample rate conversion being used. |
235 class SinusoidalLinearChirpSource { | 235 class SinusoidalLinearChirpSource { |
236 public: | 236 public: |
237 SinusoidalLinearChirpSource(int sample_rate, int samples, | 237 SinusoidalLinearChirpSource(int sample_rate, |
| 238 int samples, |
238 double max_frequency) | 239 double max_frequency) |
239 : sample_rate_(sample_rate), | 240 : sample_rate_(sample_rate), |
240 total_samples_(samples), | 241 total_samples_(samples), |
241 max_frequency_(max_frequency), | 242 max_frequency_(max_frequency), |
242 current_index_(0) { | 243 current_index_(0) { |
243 // Chirp rate. | 244 // Chirp rate. |
244 double duration = static_cast<double>(total_samples_) / sample_rate_; | 245 double duration = static_cast<double>(total_samples_) / sample_rate_; |
245 k_ = (max_frequency_ - kMinFrequency) / duration; | 246 k_ = (max_frequency_ - kMinFrequency) / duration; |
246 } | 247 } |
247 | 248 |
248 virtual ~SinusoidalLinearChirpSource() {} | 249 virtual ~SinusoidalLinearChirpSource() {} |
249 | 250 |
250 void ProvideInput(float* destination, int frames) { | 251 void ProvideInput(int frames, float* destination) { |
251 for (int i = 0; i < frames; ++i, ++current_index_) { | 252 for (int i = 0; i < frames; ++i, ++current_index_) { |
252 // Filter out frequencies higher than Nyquist. | 253 // Filter out frequencies higher than Nyquist. |
253 if (Frequency(current_index_) > 0.5 * sample_rate_) { | 254 if (Frequency(current_index_) > 0.5 * sample_rate_) { |
254 destination[i] = 0; | 255 destination[i] = 0; |
255 } else { | 256 } else { |
256 // Calculate time in seconds. | 257 // Calculate time in seconds. |
257 double t = static_cast<double>(current_index_) / sample_rate_; | 258 double t = static_cast<double>(current_index_) / sample_rate_; |
258 | 259 |
259 // Sinusoidal linear chirp. | 260 // Sinusoidal linear chirp. |
260 destination[i] = sin(2 * M_PI * (kMinFrequency * t + (k_ / 2) * t * t)); | 261 destination[i] = sin(2 * M_PI * (kMinFrequency * t + (k_ / 2) * t * t)); |
(...skipping 49 matching lines...) Expand 10 before | Expand all | Expand 10 after Loading... |
310 | 311 |
311 // Nyquist frequency for the input sampling rate. | 312 // Nyquist frequency for the input sampling rate. |
312 double input_nyquist_freq = 0.5 * input_rate_; | 313 double input_nyquist_freq = 0.5 * input_rate_; |
313 | 314 |
314 // Source for data to be resampled. | 315 // Source for data to be resampled. |
315 SinusoidalLinearChirpSource resampler_source( | 316 SinusoidalLinearChirpSource resampler_source( |
316 input_rate_, input_samples, input_nyquist_freq); | 317 input_rate_, input_samples, input_nyquist_freq); |
317 | 318 |
318 const double io_ratio = input_rate_ / static_cast<double>(output_rate_); | 319 const double io_ratio = input_rate_ / static_cast<double>(output_rate_); |
319 SincResampler resampler( | 320 SincResampler resampler( |
320 io_ratio, | 321 io_ratio, SincResampler::kDefaultRequestSize, |
321 base::Bind(&SinusoidalLinearChirpSource::ProvideInput, | 322 base::Bind(&SinusoidalLinearChirpSource::ProvideInput, |
322 base::Unretained(&resampler_source))); | 323 base::Unretained(&resampler_source))); |
323 | 324 |
324 // Force an update to the sample rate ratio to ensure dyanmic sample rate | 325 // Force an update to the sample rate ratio to ensure dyanmic sample rate |
325 // changes are working correctly. | 326 // changes are working correctly. |
326 scoped_ptr<float[]> kernel(new float[SincResampler::kKernelStorageSize]); | 327 scoped_ptr<float[]> kernel(new float[SincResampler::kKernelStorageSize]); |
327 memcpy(kernel.get(), resampler.get_kernel_for_testing(), | 328 memcpy(kernel.get(), resampler.get_kernel_for_testing(), |
328 SincResampler::kKernelStorageSize); | 329 SincResampler::kKernelStorageSize); |
329 resampler.SetRatio(M_PI); | 330 resampler.SetRatio(M_PI); |
330 ASSERT_NE(0, memcmp(kernel.get(), resampler.get_kernel_for_testing(), | 331 ASSERT_NE(0, memcmp(kernel.get(), resampler.get_kernel_for_testing(), |
331 SincResampler::kKernelStorageSize)); | 332 SincResampler::kKernelStorageSize)); |
332 resampler.SetRatio(io_ratio); | 333 resampler.SetRatio(io_ratio); |
333 ASSERT_EQ(0, memcmp(kernel.get(), resampler.get_kernel_for_testing(), | 334 ASSERT_EQ(0, memcmp(kernel.get(), resampler.get_kernel_for_testing(), |
334 SincResampler::kKernelStorageSize)); | 335 SincResampler::kKernelStorageSize)); |
335 | 336 |
336 // TODO(dalecurtis): If we switch to AVX/SSE optimization, we'll need to | 337 // TODO(dalecurtis): If we switch to AVX/SSE optimization, we'll need to |
337 // allocate these on 32-byte boundaries and ensure they're sized % 32 bytes. | 338 // allocate these on 32-byte boundaries and ensure they're sized % 32 bytes. |
338 scoped_ptr<float[]> resampled_destination(new float[output_samples]); | 339 scoped_ptr<float[]> resampled_destination(new float[output_samples]); |
339 scoped_ptr<float[]> pure_destination(new float[output_samples]); | 340 scoped_ptr<float[]> pure_destination(new float[output_samples]); |
340 | 341 |
341 // Generate resampled signal. | 342 // Generate resampled signal. |
342 resampler.Resample(resampled_destination.get(), output_samples); | 343 resampler.Resample(output_samples, resampled_destination.get()); |
343 | 344 |
344 // Generate pure signal. | 345 // Generate pure signal. |
345 SinusoidalLinearChirpSource pure_source( | 346 SinusoidalLinearChirpSource pure_source( |
346 output_rate_, output_samples, input_nyquist_freq); | 347 output_rate_, output_samples, input_nyquist_freq); |
347 pure_source.ProvideInput(pure_destination.get(), output_samples); | 348 pure_source.ProvideInput(output_samples, pure_destination.get()); |
348 | 349 |
349 // Range of the Nyquist frequency (0.5 * min(input rate, output_rate)) which | 350 // Range of the Nyquist frequency (0.5 * min(input rate, output_rate)) which |
350 // we refer to as low and high. | 351 // we refer to as low and high. |
351 static const double kLowFrequencyNyquistRange = 0.7; | 352 static const double kLowFrequencyNyquistRange = 0.7; |
352 static const double kHighFrequencyNyquistRange = 0.9; | 353 static const double kHighFrequencyNyquistRange = 0.9; |
353 | 354 |
354 // Calculate Root-Mean-Square-Error and maximum error for the resampling. | 355 // Calculate Root-Mean-Square-Error and maximum error for the resampling. |
355 double sum_of_squares = 0; | 356 double sum_of_squares = 0; |
356 double low_freq_max_error = 0; | 357 double low_freq_max_error = 0; |
357 double high_freq_max_error = 0; | 358 double high_freq_max_error = 0; |
(...skipping 76 matching lines...) Expand 10 before | Expand all | Expand 10 after Loading... |
434 std::tr1::make_tuple(11025, 192000, kResamplingRMSError, -62.61), | 435 std::tr1::make_tuple(11025, 192000, kResamplingRMSError, -62.61), |
435 std::tr1::make_tuple(16000, 192000, kResamplingRMSError, -63.14), | 436 std::tr1::make_tuple(16000, 192000, kResamplingRMSError, -63.14), |
436 std::tr1::make_tuple(22050, 192000, kResamplingRMSError, -62.42), | 437 std::tr1::make_tuple(22050, 192000, kResamplingRMSError, -62.42), |
437 std::tr1::make_tuple(32000, 192000, kResamplingRMSError, -63.38), | 438 std::tr1::make_tuple(32000, 192000, kResamplingRMSError, -63.38), |
438 std::tr1::make_tuple(44100, 192000, kResamplingRMSError, -62.63), | 439 std::tr1::make_tuple(44100, 192000, kResamplingRMSError, -62.63), |
439 std::tr1::make_tuple(48000, 192000, kResamplingRMSError, -73.44), | 440 std::tr1::make_tuple(48000, 192000, kResamplingRMSError, -73.44), |
440 std::tr1::make_tuple(96000, 192000, kResamplingRMSError, -73.52), | 441 std::tr1::make_tuple(96000, 192000, kResamplingRMSError, -73.52), |
441 std::tr1::make_tuple(192000, 192000, kResamplingRMSError, -73.52))); | 442 std::tr1::make_tuple(192000, 192000, kResamplingRMSError, -73.52))); |
442 | 443 |
443 } // namespace media | 444 } // namespace media |
OLD | NEW |