Chris@41: /*
Chris@41: ** Copyright (c) 2002-2016, Erik de Castro Lopo <erikd@mega-nerd.com>
Chris@41: ** All rights reserved.
Chris@41: **
Chris@41: ** This code is released under 2-clause BSD license. Please see the
Chris@41: ** file at : https://github.com/erikd/libsamplerate/blob/master/COPYING
Chris@41: */
Chris@41: 
Chris@41: #include "config.h"
Chris@41: 
Chris@41: #include <stdio.h>
Chris@41: #include <stdlib.h>
Chris@41: #include <string.h>
Chris@41: #include <math.h>
Chris@41: #include <time.h>
Chris@41: 
Chris@41: #if (HAVE_FFTW3)
Chris@41: 
Chris@41: #include <fftw3.h>
Chris@41: 
Chris@41: #include <samplerate.h>
Chris@41: 
Chris@41: #include "util.h"
Chris@41: 
Chris@41: #define	BUFFER_LEN		50000
Chris@41: #define	MAX_FREQS		4
Chris@41: #define	MAX_RATIOS		6
Chris@41: #define	MAX_SPEC_LEN	(1<<15)
Chris@41: 
Chris@41: #ifndef	M_PI
Chris@41: #define	M_PI			3.14159265358979323846264338
Chris@41: #endif
Chris@41: 
Chris@41: enum
Chris@41: {	BOOLEAN_FALSE	= 0,
Chris@41: 	BOOLEAN_TRUE	= 1
Chris@41: } ;
Chris@41: 
Chris@41: typedef struct
Chris@41: {	int		freq_count ;
Chris@41: 	double	freqs [MAX_FREQS] ;
Chris@41: 
Chris@41: 	double	src_ratio ;
Chris@41: 	int		pass_band_peaks ;
Chris@41: 
Chris@41: 	double	snr ;
Chris@41: 	double	peak_value ;
Chris@41: } SINGLE_TEST ;
Chris@41: 
Chris@41: typedef struct
Chris@41: {	int			converter ;
Chris@41: 	int			tests ;
Chris@41: 	int			do_bandwidth_test ;
Chris@41: 	SINGLE_TEST	test_data [10] ;
Chris@41: } CONVERTER_TEST ;
Chris@41: 
Chris@41: static double snr_test (SINGLE_TEST *snr_test_data, int number, int converter, int verbose) ;
Chris@41: static double find_peak (float *output, int output_len) ;
Chris@41: static double bandwidth_test (int converter, int verbose) ;
Chris@41: 
Chris@41: int
Chris@41: main (int argc, char *argv [])
Chris@41: {	CONVERTER_TEST snr_test_data [] =
Chris@41: 	{
Chris@41: 		{	SRC_ZERO_ORDER_HOLD,
Chris@41: 			8,
Chris@41: 			BOOLEAN_FALSE,
Chris@41: 			{	{	1,	{ 0.01111111111 },		3.0,		1,	 28.0,	1.0 },
Chris@41: 				{	1,	{ 0.01111111111 },		0.6,		1,	 36.0,	1.0 },
Chris@41: 				{	1,	{ 0.01111111111 },		0.3,		1,	 36.0,	1.0 },
Chris@41: 				{	1,	{ 0.01111111111 },		1.0,		1,	150.0,	1.0 },
Chris@41: 				{	1,	{ 0.01111111111 },		1.001,		1,	 38.0,	1.0 },
Chris@41: 				{	2,	{ 0.011111, 0.324 },	1.9999,		2,	 14.0,	1.0 },
Chris@41: 				{	2,	{ 0.012345, 0.457 },	0.456789,	1,	 12.0,	1.0 },
Chris@41: 				{	1,	{ 0.3511111111 },		1.33,		1,	 10.0,	1.0 }
Chris@41: 				}
Chris@41: 			},
Chris@41: 
Chris@41: 		{	SRC_LINEAR,
Chris@41: 			8,
Chris@41: 			BOOLEAN_FALSE,
Chris@41: 			{	{	1,	{ 0.01111111111 },		3.0,		1,	 73.0,	1.0 },
Chris@41: 				{	1,	{ 0.01111111111 },		0.6,		1,	 73.0,	1.0 },
Chris@41: 				{	1,	{ 0.01111111111 },		0.3,		1,	 73.0,	1.0 },
Chris@41: 				{	1,	{ 0.01111111111 },		1.0,		1,	150.0,	1.0 },
Chris@41: 				{	1,	{ 0.01111111111 },		1.001,		1,	 77.0,	1.0 },
Chris@41: 				{	2,	{ 0.011111, 0.324 },	1.9999,		2,	 15.0,	0.94 },
Chris@41: 				{	2,	{ 0.012345, 0.457 },	0.456789,	1,	 25.0,	0.96 },
Chris@41: 				{	1,	{ 0.3511111111 },		1.33,		1,	 22.0,	0.99 }
Chris@41: 				}
Chris@41: 			},
Chris@41: 
Chris@41: 		{	SRC_SINC_FASTEST,
Chris@41: 			9,
Chris@41: 			BOOLEAN_TRUE,
Chris@41: 			{	{	1,	{ 0.01111111111 },		3.0,		1,	100.0,	1.0 },
Chris@41: 				{	1,	{ 0.01111111111 },		0.6,		1,	 99.0,	1.0 },
Chris@41: 				{	1,	{ 0.01111111111 },		0.3,		1,	100.0,	1.0 },
Chris@41: 				{	1,	{ 0.01111111111 },		1.0,		1,	150.0,	1.0 },
Chris@41: 				{	1,	{ 0.01111111111 },		1.001,		1,	100.0,	1.0 },
Chris@41: 				{	2,	{ 0.011111, 0.324 },	1.9999,		2,	 97.0,	1.0 },
Chris@41: 				{	2,	{ 0.012345, 0.457 },	0.456789,	1,	100.0,	0.5 },
Chris@41: 				{	2,	{ 0.011111, 0.45 },		0.6,		1,	 97.0,	0.5 },
Chris@41: 				{	1,	{ 0.3511111111 },		1.33,		1,	 97.0,	1.0 }
Chris@41: 				}
Chris@41: 			},
Chris@41: 
Chris@41: 		{	SRC_SINC_MEDIUM_QUALITY,
Chris@41: 			9,
Chris@41: 			BOOLEAN_TRUE,
Chris@41: 			{	{	1,	{ 0.01111111111 },		3.0,		1,	145.0,	1.0 },
Chris@41: 				{	1,	{ 0.01111111111 },		0.6,		1,	132.0,	1.0 },
Chris@41: 				{	1,	{ 0.01111111111 },		0.3,		1,	138.0,	1.0 },
Chris@41: 				{	1,	{ 0.01111111111 },		1.0,		1,	157.0,	1.0 },
Chris@41: 				{	1,	{ 0.01111111111 },		1.001,		1,	148.0,	1.0 },
Chris@41: 				{	2,	{ 0.011111, 0.324 },	1.9999,		2,	127.0,	1.0 },
Chris@41: 				{	2,	{ 0.012345, 0.457 },	0.456789,	1,	123.0,	0.5 },
Chris@41: 				{	2,	{ 0.011111, 0.45 },		0.6,		1,	126.0,	0.5 },
Chris@41: 				{	1,	{ 0.43111111111 },		1.33,		1,	121.0,	1.0 }
Chris@41: 				}
Chris@41: 			},
Chris@41: 
Chris@41: 		{	SRC_SINC_BEST_QUALITY,
Chris@41: 			9,
Chris@41: 			BOOLEAN_TRUE,
Chris@41: 			{	{	1,	{ 0.01111111111 },		3.0,		1,	147.0,	1.0 },
Chris@41: 				{	1,	{ 0.01111111111 },		0.6,		1,	147.0,	1.0 },
Chris@41: 				{	1,	{ 0.01111111111 },		0.3,		1,	148.0,	1.0 },
Chris@41: 				{	1,	{ 0.01111111111 },		1.0,		1,	155.0,	1.0 },
Chris@41: 				{	1,	{ 0.01111111111 },		1.001,		1,	148.0,	1.0 },
Chris@41: 				{	2,	{ 0.011111, 0.324 },	1.9999,		2,	146.0,	1.0 },
Chris@41: 				{	2,	{ 0.012345, 0.457 },	0.456789,	1,	147.0,	0.5 },
Chris@41: 				{	2,	{ 0.011111, 0.45 },		0.6,		1,	144.0,	0.5 },
Chris@41: 				{	1,	{ 0.43111111111 },		1.33,		1,	145.0,	1.0 }
Chris@41: 				}
Chris@41: 			},
Chris@41: 		} ; /* snr_test_data */
Chris@41: 
Chris@41: 	double	best_snr, snr, freq3dB ;
Chris@41: 	int 	j, k, converter, verbose = 0 ;
Chris@41: 
Chris@41: 	if (argc == 2 && strcmp (argv [1], "--verbose") == 0)
Chris@41: 		verbose = 1 ;
Chris@41: 
Chris@41: 	puts ("") ;
Chris@41: 
Chris@41: 	for (j = 0 ; j < ARRAY_LEN (snr_test_data) ; j++)
Chris@41: 	{	best_snr = 5000.0 ;
Chris@41: 
Chris@41: 		converter = snr_test_data [j].converter ;
Chris@41: 
Chris@41: 		printf ("    Converter %d : %s\n", converter, src_get_name (converter)) ;
Chris@41: 		printf ("    %s\n", src_get_description (converter)) ;
Chris@41: 
Chris@41: 		for (k = 0 ; k < snr_test_data [j].tests ; k++)
Chris@41: 		{	snr = snr_test (&(snr_test_data [j].test_data [k]), k, converter, verbose) ;
Chris@41: 			if (best_snr > snr)
Chris@41: 				best_snr = snr ;
Chris@41: 			} ;
Chris@41: 
Chris@41: 		printf ("    Worst case Signal-to-Noise Ratio : %.2f dB.\n", best_snr) ;
Chris@41: 
Chris@41: 		if (snr_test_data [j].do_bandwidth_test == BOOLEAN_FALSE)
Chris@41: 		{	puts ("    Bandwith test not performed on this converter.\n") ;
Chris@41: 			continue ;
Chris@41: 			}
Chris@41: 
Chris@41: 		freq3dB = bandwidth_test (converter, verbose) ;
Chris@41: 
Chris@41: 		printf ("    Measured -3dB rolloff point      : %5.2f %%.\n\n", freq3dB) ;
Chris@41: 		} ;
Chris@41: 
Chris@41: 	fftw_cleanup () ;
Chris@41: 
Chris@41: 	return 0 ;
Chris@41: } /* main */
Chris@41: 
Chris@41: /*==============================================================================
Chris@41: */
Chris@41: 
Chris@41: static double
Chris@41: snr_test (SINGLE_TEST *test_data, int number, int converter, int verbose)
Chris@41: {	static float data [BUFFER_LEN + 1] ;
Chris@41: 	static float output [MAX_SPEC_LEN] ;
Chris@41: 
Chris@41: 	SRC_STATE	*src_state ;
Chris@41: 	SRC_DATA	src_data ;
Chris@41: 
Chris@41: 	double		output_peak, snr ;
Chris@41: 	int 		k, output_len, input_len, error ;
Chris@41: 
Chris@41: 	if (verbose != 0)
Chris@41: 	{	printf ("\tSignal-to-Noise Ratio Test %d.\n"
Chris@41: 				"\t=====================================\n", number) ;
Chris@41: 		printf ("\tFrequencies : [ ") ;
Chris@41: 		for (k = 0 ; k < test_data->freq_count ; k++)
Chris@41: 			printf ("%6.4f ", test_data->freqs [k]) ;
Chris@41: 
Chris@41: 		printf ("]\n\tSRC Ratio   : %8.4f\n", test_data->src_ratio) ;
Chris@41: 		}
Chris@41: 	else
Chris@41: 	{	printf ("\tSignal-to-Noise Ratio Test %d : ", number) ;
Chris@41: 		fflush (stdout) ;
Chris@41: 		} ;
Chris@41: 
Chris@41: 	/* Set up the output array. */
Chris@41: 	if (test_data->src_ratio >= 1.0)
Chris@41: 	{	output_len = MAX_SPEC_LEN ;
Chris@41: 		input_len = (int) ceil (MAX_SPEC_LEN / test_data->src_ratio) ;
Chris@41: 		if (input_len > BUFFER_LEN)
Chris@41: 			input_len = BUFFER_LEN ;
Chris@41: 		}
Chris@41: 	else
Chris@41: 	{	input_len = BUFFER_LEN ;
Chris@41: 		output_len = (int) ceil (BUFFER_LEN * test_data->src_ratio) ;
Chris@41: 		output_len &= ((~0u) << 4) ;
Chris@41: 		if (output_len > MAX_SPEC_LEN)
Chris@41: 			output_len = MAX_SPEC_LEN ;
Chris@41: 		input_len = (int) ceil (output_len / test_data->src_ratio) ;
Chris@41: 		} ;
Chris@41: 
Chris@41: 	memset (output, 0, sizeof (output)) ;
Chris@41: 
Chris@41: 	/* Generate input data array. */
Chris@41: 	gen_windowed_sines (test_data->freq_count, test_data->freqs, 1.0, data, input_len) ;
Chris@41: 
Chris@41: 	/* Perform sample rate conversion. */
Chris@41: 	if ((src_state = src_new (converter, 1, &error)) == NULL)
Chris@41: 	{	printf ("\n\nLine %d : src_new() failed : %s.\n\n", __LINE__, src_strerror (error)) ;
Chris@41: 		exit (1) ;
Chris@41: 		} ;
Chris@41: 
Chris@41: 	src_data.end_of_input = 1 ; /* Only one buffer worth of input. */
Chris@41: 
Chris@41: 	src_data.data_in = data ;
Chris@41: 	src_data.input_frames = input_len ;
Chris@41: 
Chris@41: 	src_data.src_ratio = test_data->src_ratio ;
Chris@41: 
Chris@41: 	src_data.data_out = output ;
Chris@41: 	src_data.output_frames = output_len ;
Chris@41: 
Chris@41: 	if ((error = src_process (src_state, &src_data)))
Chris@41: 	{	printf ("\n\nLine %d : %s\n\n", __LINE__, src_strerror (error)) ;
Chris@41: 		exit (1) ;
Chris@41: 		} ;
Chris@41: 
Chris@41: 	src_state = src_delete (src_state) ;
Chris@41: 
Chris@41: 	if (verbose != 0)
Chris@41: 		printf ("\tOutput Len  :   %ld\n", src_data.output_frames_gen) ;
Chris@41: 
Chris@41: 	if (abs (src_data.output_frames_gen - output_len) > 4)
Chris@41: 	{	printf ("\n\nLine %d : output data length should be %d.\n\n", __LINE__, output_len) ;
Chris@41: 		exit (1) ;
Chris@41: 		} ;
Chris@41: 
Chris@41: 	/* Check output peak. */
Chris@41: 	output_peak = find_peak (output, src_data.output_frames_gen) ;
Chris@41: 
Chris@41: 	if (verbose != 0)
Chris@41: 		printf ("\tOutput Peak :   %6.4f\n", output_peak) ;
Chris@41: 
Chris@41: 	if (fabs (output_peak - test_data->peak_value) > 0.01)
Chris@41: 	{	printf ("\n\nLine %d : output peak (%6.4f) should be %6.4f\n\n", __LINE__, output_peak, test_data->peak_value) ;
Chris@41: 		save_oct_float ("snr_test.dat", data, BUFFER_LEN, output, output_len) ;
Chris@41: 		exit (1) ;
Chris@41: 		} ;
Chris@41: 
Chris@41: 	/* Calculate signal-to-noise ratio. */
Chris@41: 	snr = calculate_snr (output, src_data.output_frames_gen, test_data->pass_band_peaks) ;
Chris@41: 
Chris@41: 	if (snr < 0.0)
Chris@41: 	{	/* An error occurred. */
Chris@41: 		save_oct_float ("snr_test.dat", data, BUFFER_LEN, output, src_data.output_frames_gen) ;
Chris@41: 		exit (1) ;
Chris@41: 		} ;
Chris@41: 
Chris@41: 	if (verbose != 0)
Chris@41: 		printf ("\tSNR Ratio   :   %.2f dB\n", snr) ;
Chris@41: 
Chris@41: 	if (snr < test_data->snr)
Chris@41: 	{	printf ("\n\nLine %d : SNR (%5.2f) should be > %6.2f dB\n\n", __LINE__, snr, test_data->snr) ;
Chris@41: 		exit (1) ;
Chris@41: 		} ;
Chris@41: 
Chris@41: 	if (verbose != 0)
Chris@41: 		puts ("\t-------------------------------------\n\tPass\n") ;
Chris@41: 	else
Chris@41: 		puts ("Pass") ;
Chris@41: 
Chris@41: 	return snr ;
Chris@41: } /* snr_test */
Chris@41: 
Chris@41: static double
Chris@41: find_peak (float *data, int len)
Chris@41: {	double 	peak = 0.0 ;
Chris@41: 	int		k = 0 ;
Chris@41: 
Chris@41: 	for (k = 0 ; k < len ; k++)
Chris@41: 		if (fabs (data [k]) > peak)
Chris@41: 			peak = fabs (data [k]) ;
Chris@41: 
Chris@41: 	return peak ;
Chris@41: } /* find_peak */
Chris@41: 
Chris@41: 
Chris@41: static double
Chris@41: find_attenuation (double freq, int converter, int verbose)
Chris@41: {	static float input	[BUFFER_LEN] ;
Chris@41: 	static float output [2 * BUFFER_LEN] ;
Chris@41: 
Chris@41: 	SRC_DATA	src_data ;
Chris@41: 	double 		output_peak ;
Chris@41: 	int			error ;
Chris@41: 
Chris@41: 	gen_windowed_sines (1, &freq, 1.0, input, BUFFER_LEN) ;
Chris@41: 
Chris@41: 	src_data.end_of_input = 1 ; /* Only one buffer worth of input. */
Chris@41: 
Chris@41: 	src_data.data_in = input ;
Chris@41: 	src_data.input_frames = BUFFER_LEN ;
Chris@41: 
Chris@41: 	src_data.src_ratio = 1.999 ;
Chris@41: 
Chris@41: 	src_data.data_out = output ;
Chris@41: 	src_data.output_frames = ARRAY_LEN (output) ;
Chris@41: 
Chris@41: 	if ((error = src_simple (&src_data, converter, 1)))
Chris@41: 	{	printf ("\n\nLine %d : %s\n\n", __LINE__, src_strerror (error)) ;
Chris@41: 		exit (1) ;
Chris@41: 		} ;
Chris@41: 
Chris@41: 	output_peak = find_peak (output, ARRAY_LEN (output)) ;
Chris@41: 
Chris@41: 	if (verbose)
Chris@41: 		printf ("\tFreq : %6f   InPeak : %6f    OutPeak : %6f   Atten : %6.2f dB\n",
Chris@41: 				freq, 1.0, output_peak, 20.0 * log10 (1.0 / output_peak)) ;
Chris@41: 
Chris@41: 	return 20.0 * log10 (1.0 / output_peak) ;
Chris@41: } /* find_attenuation */
Chris@41: 
Chris@41: static double
Chris@41: bandwidth_test (int converter, int verbose)
Chris@41: {	double	f1, f2, a1, a2 ;
Chris@41: 	double	freq, atten ;
Chris@41: 
Chris@41: 	f1 = 0.35 ;
Chris@41: 	a1 = find_attenuation (f1, converter, verbose) ;
Chris@41: 
Chris@41: 	f2 = 0.495 ;
Chris@41: 	a2 = find_attenuation (f2, converter, verbose) ;
Chris@41: 
Chris@41: 	if (a1 > 3.0 || a2 < 3.0)
Chris@41: 	{	printf ("\n\nLine %d : cannot bracket 3dB point.\n\n", __LINE__) ;
Chris@41: 		exit (1) ;
Chris@41: 		} ;
Chris@41: 
Chris@41: 	while (a2 - a1 > 1.0)
Chris@41: 	{	freq = f1 + 0.5 * (f2 - f1) ;
Chris@41: 		atten = find_attenuation (freq, converter, verbose) ;
Chris@41: 
Chris@41: 		if (atten < 3.0)
Chris@41: 		{	f1 = freq ;
Chris@41: 			a1 = atten ;
Chris@41: 			}
Chris@41: 		else
Chris@41: 		{	f2 = freq ;
Chris@41: 			a2 = atten ;
Chris@41: 			} ;
Chris@41: 		} ;
Chris@41: 
Chris@41: 	freq = f1 + (3.0 - a1) * (f2 - f1) / (a2 - a1) ;
Chris@41: 
Chris@41: 	return 200.0 * freq ;
Chris@41: } /* bandwidth_test */
Chris@41: 
Chris@41: #else /* (HAVE_FFTW3) == 0 */
Chris@41: 
Chris@41: /* Alternative main function when librfftw is not available. */
Chris@41: 
Chris@41: int
Chris@41: main (void)
Chris@41: {	puts ("\n"
Chris@41: 		"****************************************************************\n"
Chris@41: 		" This test cannot be run without FFTW (http://www.fftw.org/).\n"
Chris@41: 		" Both the real and the complex versions of the library are\n"
Chris@41: 		" required.") ;
Chris@41: 	puts ("****************************************************************\n") ;
Chris@41: 
Chris@41: 	return 0 ;
Chris@41: } /* main */
Chris@41: 
Chris@41: #endif
Chris@41: