15821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)// Copyright (c) 2012 The Chromium Authors. All rights reserved. 25821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)// Use of this source code is governed by a BSD-style license that can be 35821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)// found in the LICENSE file. 45821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 55821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)// MSVC++ requires this to be set before any other includes to get M_PI. 65821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)#define _USE_MATH_DEFINES 75821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 85821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)#include <cmath> 95821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 105821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)#include "base/bind.h" 115821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)#include "base/bind_helpers.h" 122a99a7e74a7f215066514fe81d2bfa6639d9edddTorne (Richard Coles)#include "base/cpu.h" 13868fa2fe829687343ffae624259930155e16dbd8Torne (Richard Coles)#include "base/strings/string_number_conversions.h" 14eb525c5499e34cc9c4b825d6d9e75bb07cc06aceBen Murdoch#include "base/time/time.h" 155821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)#include "build/build_config.h" 165821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)#include "media/base/sinc_resampler.h" 175821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)#include "testing/gmock/include/gmock/gmock.h" 185821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)#include "testing/gtest/include/gtest/gtest.h" 195821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 205821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)using testing::_; 215821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 225821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)namespace media { 235821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 245821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)static const double kSampleRateRatio = 192000.0 / 44100.0; 255821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 265821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)// Helper class to ensure ChunkedResample() functions properly. 275821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)class MockSource { 285821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) public: 29b2df76ea8fec9e32f6f3718986dba0d95315b29cTorne (Richard Coles) MOCK_METHOD2(ProvideInput, void(int frames, float* destination)); 305821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)}; 315821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 325821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)ACTION(ClearBuffer) { 33b2df76ea8fec9e32f6f3718986dba0d95315b29cTorne (Richard Coles) memset(arg1, 0, arg0 * sizeof(float)); 345821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)} 355821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 365821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)ACTION(FillBuffer) { 375821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // Value chosen arbitrarily such that SincResampler resamples it to something 385821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // easily representable on all platforms; e.g., using kSampleRateRatio this 395821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // becomes 1.81219. 40b2df76ea8fec9e32f6f3718986dba0d95315b29cTorne (Richard Coles) memset(arg1, 64, arg0 * sizeof(float)); 415821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)} 425821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 435821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)// Test requesting multiples of ChunkSize() frames results in the proper number 445821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)// of callbacks. 455821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)TEST(SincResamplerTest, ChunkedResample) { 465821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) MockSource mock_source; 475821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 485821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // Choose a high ratio of input to output samples which will result in quick 495821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // exhaustion of SincResampler's internal buffers. 505821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) SincResampler resampler( 51b2df76ea8fec9e32f6f3718986dba0d95315b29cTorne (Richard Coles) kSampleRateRatio, SincResampler::kDefaultRequestSize, 525821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) base::Bind(&MockSource::ProvideInput, base::Unretained(&mock_source))); 535821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 545821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) static const int kChunks = 2; 555821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) int max_chunk_size = resampler.ChunkSize() * kChunks; 56c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) scoped_ptr<float[]> resampled_destination(new float[max_chunk_size]); 575821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 585821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // Verify requesting ChunkSize() frames causes a single callback. 595821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) EXPECT_CALL(mock_source, ProvideInput(_, _)) 605821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) .Times(1).WillOnce(ClearBuffer()); 61b2df76ea8fec9e32f6f3718986dba0d95315b29cTorne (Richard Coles) resampler.Resample(resampler.ChunkSize(), resampled_destination.get()); 625821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 635821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // Verify requesting kChunks * ChunkSize() frames causes kChunks callbacks. 645821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) testing::Mock::VerifyAndClear(&mock_source); 655821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) EXPECT_CALL(mock_source, ProvideInput(_, _)) 665821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) .Times(kChunks).WillRepeatedly(ClearBuffer()); 67b2df76ea8fec9e32f6f3718986dba0d95315b29cTorne (Richard Coles) resampler.Resample(max_chunk_size, resampled_destination.get()); 685821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)} 695821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 705821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)// Test flush resets the internal state properly. 715821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)TEST(SincResamplerTest, Flush) { 725821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) MockSource mock_source; 735821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) SincResampler resampler( 74b2df76ea8fec9e32f6f3718986dba0d95315b29cTorne (Richard Coles) kSampleRateRatio, SincResampler::kDefaultRequestSize, 755821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) base::Bind(&MockSource::ProvideInput, base::Unretained(&mock_source))); 76c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) scoped_ptr<float[]> resampled_destination(new float[resampler.ChunkSize()]); 775821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 785821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // Fill the resampler with junk data. 795821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) EXPECT_CALL(mock_source, ProvideInput(_, _)) 805821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) .Times(1).WillOnce(FillBuffer()); 81b2df76ea8fec9e32f6f3718986dba0d95315b29cTorne (Richard Coles) resampler.Resample(resampler.ChunkSize() / 2, resampled_destination.get()); 825821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) ASSERT_NE(resampled_destination[0], 0); 835821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 845821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // Flush and request more data, which should all be zeros now. 855821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) resampler.Flush(); 865821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) testing::Mock::VerifyAndClear(&mock_source); 875821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) EXPECT_CALL(mock_source, ProvideInput(_, _)) 885821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) .Times(1).WillOnce(ClearBuffer()); 89b2df76ea8fec9e32f6f3718986dba0d95315b29cTorne (Richard Coles) resampler.Resample(resampler.ChunkSize() / 2, resampled_destination.get()); 905821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) for (int i = 0; i < resampler.ChunkSize() / 2; ++i) 915821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) ASSERT_FLOAT_EQ(resampled_destination[i], 0); 925821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)} 935821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 94c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles)// Test flush resets the internal state properly. 95c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles)TEST(SincResamplerTest, DISABLED_SetRatioBench) { 96c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) MockSource mock_source; 97c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) SincResampler resampler( 98b2df76ea8fec9e32f6f3718986dba0d95315b29cTorne (Richard Coles) kSampleRateRatio, SincResampler::kDefaultRequestSize, 99c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) base::Bind(&MockSource::ProvideInput, base::Unretained(&mock_source))); 100c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) 101c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) base::TimeTicks start = base::TimeTicks::HighResNow(); 102c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) for (int i = 1; i < 10000; ++i) 103c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) resampler.SetRatio(1.0 / i); 104c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) double total_time_c_ms = 105c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) (base::TimeTicks::HighResNow() - start).InMillisecondsF(); 106c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) printf("SetRatio() took %.2fms.\n", total_time_c_ms); 107c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles)} 108c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) 109c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) 1105821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)// Define platform independent function name for Convolve* tests. 1112a99a7e74a7f215066514fe81d2bfa6639d9edddTorne (Richard Coles)#if defined(ARCH_CPU_X86_FAMILY) 1125821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)#define CONVOLVE_FUNC Convolve_SSE 1135821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)#elif defined(ARCH_CPU_ARM_FAMILY) && defined(USE_NEON) 1145821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)#define CONVOLVE_FUNC Convolve_NEON 1155821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)#endif 1165821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 1175821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)// Ensure various optimized Convolve() methods return the same value. Only run 1185821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)// this test if other optimized methods exist, otherwise the default Convolve() 1195821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)// will be tested by the parameterized SincResampler tests below. 1205821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)#if defined(CONVOLVE_FUNC) 1211e9bf3e0803691d0a228da41fc608347b6db4340Torne (Richard Coles)static const double kKernelInterpolationFactor = 0.5; 1221e9bf3e0803691d0a228da41fc608347b6db4340Torne (Richard Coles) 1235821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)TEST(SincResamplerTest, Convolve) { 1242a99a7e74a7f215066514fe81d2bfa6639d9edddTorne (Richard Coles)#if defined(ARCH_CPU_X86_FAMILY) 1252a99a7e74a7f215066514fe81d2bfa6639d9edddTorne (Richard Coles) ASSERT_TRUE(base::CPU().has_sse()); 1262a99a7e74a7f215066514fe81d2bfa6639d9edddTorne (Richard Coles)#endif 1272a99a7e74a7f215066514fe81d2bfa6639d9edddTorne (Richard Coles) 1285821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // Initialize a dummy resampler. 1295821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) MockSource mock_source; 1305821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) SincResampler resampler( 131b2df76ea8fec9e32f6f3718986dba0d95315b29cTorne (Richard Coles) kSampleRateRatio, SincResampler::kDefaultRequestSize, 1325821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) base::Bind(&MockSource::ProvideInput, base::Unretained(&mock_source))); 1335821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 1345821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // The optimized Convolve methods are slightly more precise than Convolve_C(), 1355821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // so comparison must be done using an epsilon. 1365821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) static const double kEpsilon = 0.00000005; 1375821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 1385821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // Use a kernel from SincResampler as input and kernel data, this has the 1395821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // benefit of already being properly sized and aligned for Convolve_SSE(). 1405821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) double result = resampler.Convolve_C( 1415821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) resampler.kernel_storage_.get(), resampler.kernel_storage_.get(), 1425821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) resampler.kernel_storage_.get(), kKernelInterpolationFactor); 1435821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) double result2 = resampler.CONVOLVE_FUNC( 1445821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) resampler.kernel_storage_.get(), resampler.kernel_storage_.get(), 1455821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) resampler.kernel_storage_.get(), kKernelInterpolationFactor); 1465821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) EXPECT_NEAR(result2, result, kEpsilon); 1475821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 1485821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // Test Convolve() w/ unaligned input pointer. 1495821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) result = resampler.Convolve_C( 1505821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) resampler.kernel_storage_.get() + 1, resampler.kernel_storage_.get(), 1515821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) resampler.kernel_storage_.get(), kKernelInterpolationFactor); 1525821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) result2 = resampler.CONVOLVE_FUNC( 1535821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) resampler.kernel_storage_.get() + 1, resampler.kernel_storage_.get(), 1545821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) resampler.kernel_storage_.get(), kKernelInterpolationFactor); 1555821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) EXPECT_NEAR(result2, result, kEpsilon); 1565821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)} 1575821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)#endif 1585821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 1595821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)// Fake audio source for testing the resampler. Generates a sinusoidal linear 1605821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)// chirp (http://en.wikipedia.org/wiki/Chirp) which can be tuned to stress the 1615821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)// resampler for the specific sample rate conversion being used. 1625821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)class SinusoidalLinearChirpSource { 1635821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) public: 164b2df76ea8fec9e32f6f3718986dba0d95315b29cTorne (Richard Coles) SinusoidalLinearChirpSource(int sample_rate, 165b2df76ea8fec9e32f6f3718986dba0d95315b29cTorne (Richard Coles) int samples, 1665821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) double max_frequency) 1675821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) : sample_rate_(sample_rate), 1685821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) total_samples_(samples), 1695821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) max_frequency_(max_frequency), 1705821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) current_index_(0) { 1715821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // Chirp rate. 1725821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) double duration = static_cast<double>(total_samples_) / sample_rate_; 1735821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) k_ = (max_frequency_ - kMinFrequency) / duration; 1745821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) } 1755821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 1765821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) virtual ~SinusoidalLinearChirpSource() {} 1775821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 178b2df76ea8fec9e32f6f3718986dba0d95315b29cTorne (Richard Coles) void ProvideInput(int frames, float* destination) { 1795821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) for (int i = 0; i < frames; ++i, ++current_index_) { 1805821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // Filter out frequencies higher than Nyquist. 1815821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) if (Frequency(current_index_) > 0.5 * sample_rate_) { 1825821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) destination[i] = 0; 1835821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) } else { 1845821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // Calculate time in seconds. 1855821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) double t = static_cast<double>(current_index_) / sample_rate_; 1865821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 1875821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // Sinusoidal linear chirp. 1885821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) destination[i] = sin(2 * M_PI * (kMinFrequency * t + (k_ / 2) * t * t)); 1895821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) } 1905821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) } 1915821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) } 1925821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 1935821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) double Frequency(int position) { 1945821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) return kMinFrequency + position * (max_frequency_ - kMinFrequency) 1955821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) / total_samples_; 1965821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) } 1975821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 1985821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) private: 1995821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) enum { 2005821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) kMinFrequency = 5 2015821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) }; 2025821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 2035821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) double sample_rate_; 2045821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) int total_samples_; 2055821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) double max_frequency_; 2065821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) double k_; 2075821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) int current_index_; 2085821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 2095821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) DISALLOW_COPY_AND_ASSIGN(SinusoidalLinearChirpSource); 2105821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)}; 2115821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 2125821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)typedef std::tr1::tuple<int, int, double, double> SincResamplerTestData; 2135821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)class SincResamplerTest 2145821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) : public testing::TestWithParam<SincResamplerTestData> { 2155821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) public: 2165821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) SincResamplerTest() 2175821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) : input_rate_(std::tr1::get<0>(GetParam())), 2185821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) output_rate_(std::tr1::get<1>(GetParam())), 2195821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) rms_error_(std::tr1::get<2>(GetParam())), 2205821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) low_freq_error_(std::tr1::get<3>(GetParam())) { 2215821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) } 2225821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 2235821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) virtual ~SincResamplerTest() {} 2245821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 2255821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) protected: 2265821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) int input_rate_; 2275821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) int output_rate_; 2285821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) double rms_error_; 2295821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) double low_freq_error_; 2305821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)}; 2315821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 2325821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)// Tests resampling using a given input and output sample rate. 2335821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)TEST_P(SincResamplerTest, Resample) { 2345821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // Make comparisons using one second of data. 2355821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) static const double kTestDurationSecs = 1; 2365821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) int input_samples = kTestDurationSecs * input_rate_; 2375821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) int output_samples = kTestDurationSecs * output_rate_; 2385821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 2395821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // Nyquist frequency for the input sampling rate. 2405821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) double input_nyquist_freq = 0.5 * input_rate_; 2415821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 2425821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // Source for data to be resampled. 2435821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) SinusoidalLinearChirpSource resampler_source( 2445821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) input_rate_, input_samples, input_nyquist_freq); 2455821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 246c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) const double io_ratio = input_rate_ / static_cast<double>(output_rate_); 2475821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) SincResampler resampler( 248b2df76ea8fec9e32f6f3718986dba0d95315b29cTorne (Richard Coles) io_ratio, SincResampler::kDefaultRequestSize, 2495821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) base::Bind(&SinusoidalLinearChirpSource::ProvideInput, 2505821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) base::Unretained(&resampler_source))); 2515821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 252c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) // Force an update to the sample rate ratio to ensure dyanmic sample rate 253c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) // changes are working correctly. 254c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) scoped_ptr<float[]> kernel(new float[SincResampler::kKernelStorageSize]); 255c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) memcpy(kernel.get(), resampler.get_kernel_for_testing(), 256c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) SincResampler::kKernelStorageSize); 257c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) resampler.SetRatio(M_PI); 258c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) ASSERT_NE(0, memcmp(kernel.get(), resampler.get_kernel_for_testing(), 259c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) SincResampler::kKernelStorageSize)); 260c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) resampler.SetRatio(io_ratio); 261c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) ASSERT_EQ(0, memcmp(kernel.get(), resampler.get_kernel_for_testing(), 262c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) SincResampler::kKernelStorageSize)); 263c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) 2645821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // TODO(dalecurtis): If we switch to AVX/SSE optimization, we'll need to 2655821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // allocate these on 32-byte boundaries and ensure they're sized % 32 bytes. 266c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) scoped_ptr<float[]> resampled_destination(new float[output_samples]); 267c2e0dbddbe15c98d52c4786dac06cb8952a8ae6dTorne (Richard Coles) scoped_ptr<float[]> pure_destination(new float[output_samples]); 2685821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 2695821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // Generate resampled signal. 270b2df76ea8fec9e32f6f3718986dba0d95315b29cTorne (Richard Coles) resampler.Resample(output_samples, resampled_destination.get()); 2715821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 2725821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // Generate pure signal. 2735821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) SinusoidalLinearChirpSource pure_source( 2745821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) output_rate_, output_samples, input_nyquist_freq); 275b2df76ea8fec9e32f6f3718986dba0d95315b29cTorne (Richard Coles) pure_source.ProvideInput(output_samples, pure_destination.get()); 2765821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 2775821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // Range of the Nyquist frequency (0.5 * min(input rate, output_rate)) which 2785821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // we refer to as low and high. 2795821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) static const double kLowFrequencyNyquistRange = 0.7; 2805821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) static const double kHighFrequencyNyquistRange = 0.9; 2815821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 2825821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // Calculate Root-Mean-Square-Error and maximum error for the resampling. 2835821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) double sum_of_squares = 0; 2845821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) double low_freq_max_error = 0; 2855821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) double high_freq_max_error = 0; 2865821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) int minimum_rate = std::min(input_rate_, output_rate_); 2875821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) double low_frequency_range = kLowFrequencyNyquistRange * 0.5 * minimum_rate; 2885821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) double high_frequency_range = kHighFrequencyNyquistRange * 0.5 * minimum_rate; 2895821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) for (int i = 0; i < output_samples; ++i) { 2905821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) double error = fabs(resampled_destination[i] - pure_destination[i]); 2915821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 2925821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) if (pure_source.Frequency(i) < low_frequency_range) { 2935821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) if (error > low_freq_max_error) 2945821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) low_freq_max_error = error; 2955821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) } else if (pure_source.Frequency(i) < high_frequency_range) { 2965821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) if (error > high_freq_max_error) 2975821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) high_freq_max_error = error; 2985821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) } 2995821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // TODO(dalecurtis): Sanity check frequencies > kHighFrequencyNyquistRange. 3005821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 3015821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) sum_of_squares += error * error; 3025821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) } 3035821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 3045821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) double rms_error = sqrt(sum_of_squares / output_samples); 3055821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 3065821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // Convert each error to dbFS. 3075821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) #define DBFS(x) 20 * log10(x) 3085821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) rms_error = DBFS(rms_error); 3095821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) low_freq_max_error = DBFS(low_freq_max_error); 3105821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) high_freq_max_error = DBFS(high_freq_max_error); 3115821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 3125821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) EXPECT_LE(rms_error, rms_error_); 3135821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) EXPECT_LE(low_freq_max_error, low_freq_error_); 3145821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 3155821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // All conversions currently have a high frequency error around -6 dbFS. 3165821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) static const double kHighFrequencyMaxError = -6.02; 3175821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) EXPECT_LE(high_freq_max_error, kHighFrequencyMaxError); 3185821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)} 3195821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 3205821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)// Almost all conversions have an RMS error of around -14 dbFS. 3215821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)static const double kResamplingRMSError = -14.58; 3225821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 3235821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)// Thresholds chosen arbitrarily based on what each resampling reported during 3245821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)// testing. All thresholds are in dbFS, http://en.wikipedia.org/wiki/DBFS. 3255821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)INSTANTIATE_TEST_CASE_P( 3265821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) SincResamplerTest, SincResamplerTest, testing::Values( 3275821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // To 44.1kHz 3285821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(8000, 44100, kResamplingRMSError, -62.73), 3295821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(11025, 44100, kResamplingRMSError, -72.19), 3305821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(16000, 44100, kResamplingRMSError, -62.54), 3315821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(22050, 44100, kResamplingRMSError, -73.53), 3325821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(32000, 44100, kResamplingRMSError, -63.32), 3335821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(44100, 44100, kResamplingRMSError, -73.53), 3345821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(48000, 44100, -15.01, -64.04), 3355821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(96000, 44100, -18.49, -25.51), 3365821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(192000, 44100, -20.50, -13.31), 3375821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 3385821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // To 48kHz 3395821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(8000, 48000, kResamplingRMSError, -63.43), 3405821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(11025, 48000, kResamplingRMSError, -62.61), 3415821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(16000, 48000, kResamplingRMSError, -63.96), 3425821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(22050, 48000, kResamplingRMSError, -62.42), 3435821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(32000, 48000, kResamplingRMSError, -64.04), 3445821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(44100, 48000, kResamplingRMSError, -62.63), 3455821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(48000, 48000, kResamplingRMSError, -73.52), 3465821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(96000, 48000, -18.40, -28.44), 3475821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(192000, 48000, -20.43, -14.11), 3485821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 3495821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // To 96kHz 3505821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(8000, 96000, kResamplingRMSError, -63.19), 3515821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(11025, 96000, kResamplingRMSError, -62.61), 3525821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(16000, 96000, kResamplingRMSError, -63.39), 3535821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(22050, 96000, kResamplingRMSError, -62.42), 3545821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(32000, 96000, kResamplingRMSError, -63.95), 3555821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(44100, 96000, kResamplingRMSError, -62.63), 3565821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(48000, 96000, kResamplingRMSError, -73.52), 3575821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(96000, 96000, kResamplingRMSError, -73.52), 3585821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(192000, 96000, kResamplingRMSError, -28.41), 3595821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 3605821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) // To 192kHz 3615821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(8000, 192000, kResamplingRMSError, -63.10), 3625821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(11025, 192000, kResamplingRMSError, -62.61), 3635821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(16000, 192000, kResamplingRMSError, -63.14), 3645821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(22050, 192000, kResamplingRMSError, -62.42), 3655821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(32000, 192000, kResamplingRMSError, -63.38), 3665821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(44100, 192000, kResamplingRMSError, -62.63), 3675821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(48000, 192000, kResamplingRMSError, -73.44), 3685821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(96000, 192000, kResamplingRMSError, -73.52), 3695821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) std::tr1::make_tuple(192000, 192000, kResamplingRMSError, -73.52))); 3705821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles) 3715821806d5e7f356e8fa4b058a389a808ea183019Torne (Richard Coles)} // namespace media 372