Chris@41: /*
Chris@41: ** Copyright (c) 2004-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 <stdio.h>
Chris@41: #include <stdlib.h>
Chris@41: #include <string.h>
Chris@41: #include <time.h>
Chris@41: #include <unistd.h>
Chris@41: 
Chris@41: #include <samplerate.h>
Chris@41: 
Chris@41: #include "config.h"
Chris@41: 
Chris@41: #include "util.h"
Chris@41: #include "float_cast.h"
Chris@41: 
Chris@41: #define BUFFER_LEN	(1<<16)
Chris@41: 
Chris@41: static float input [BUFFER_LEN] ;
Chris@41: static float output [BUFFER_LEN] ;
Chris@41: 
Chris@41: static long
Chris@41: throughput_test (int converter, long best_throughput)
Chris@41: {	SRC_DATA src_data ;
Chris@41: 	clock_t start_time, clock_time ;
Chris@41: 	double duration ;
Chris@41: 	long total_frames = 0, throughput ;
Chris@41: 	int error ;
Chris@41: 
Chris@41: 	printf ("    %-30s    ", src_get_name (converter)) ;
Chris@41: 	fflush (stdout) ;
Chris@41: 
Chris@41: 	src_data.data_in = input ;
Chris@41: 	src_data.input_frames = ARRAY_LEN (input) ;
Chris@41: 
Chris@41: 	src_data.data_out = output ;
Chris@41: 	src_data.output_frames = ARRAY_LEN (output) ;
Chris@41: 
Chris@41: 	src_data.src_ratio = 0.99 ;
Chris@41: 
Chris@41: 	sleep (2) ;
Chris@41: 
Chris@41: 	start_time = clock () ;
Chris@41: 
Chris@41: 	do
Chris@41: 	{
Chris@41: 		if ((error = src_simple (&src_data, converter, 1)) != 0)
Chris@41: 		{	puts (src_strerror (error)) ;
Chris@41: 			exit (1) ;
Chris@41: 			} ;
Chris@41: 
Chris@41: 		total_frames += src_data.output_frames_gen ;
Chris@41: 
Chris@41: 		clock_time = clock () - start_time ;
Chris@41: 		duration = (1.0 * clock_time) / CLOCKS_PER_SEC ;
Chris@41: 	}
Chris@41: 	while (duration < 3.0) ;
Chris@41: 
Chris@41: 	if (src_data.input_frames_used != ARRAY_LEN (input))
Chris@41: 	{	printf ("\n\nLine %d : input frames used %ld should be %d\n", __LINE__, src_data.input_frames_used, ARRAY_LEN (input)) ;
Chris@41: 		exit (1) ;
Chris@41: 		} ;
Chris@41: 
Chris@41: 	if (fabs (src_data.src_ratio * src_data.input_frames_used - src_data.output_frames_gen) > 2)
Chris@41: 	{	printf ("\n\nLine %d : input / output length mismatch.\n\n", __LINE__) ;
Chris@41: 		printf ("    input len  : %d\n", ARRAY_LEN (input)) ;
Chris@41: 		printf ("    output len : %ld (should be %g +/- 2)\n\n", src_data.output_frames_gen,
Chris@41: 				floor (0.5 + src_data.src_ratio * src_data.input_frames_used)) ;
Chris@41: 		exit (1) ;
Chris@41: 		} ;
Chris@41: 
Chris@41: 	throughput = lrint (floor (total_frames / duration)) ;
Chris@41: 
Chris@41: 	if (best_throughput == 0)
Chris@41: 	{	best_throughput = MAX (throughput, best_throughput) ;
Chris@41: 		printf ("%5.2f          %10ld\n", duration, throughput) ;
Chris@41: 		}
Chris@41: 	else
Chris@41: 	{	best_throughput = MAX (throughput, best_throughput) ;
Chris@41: 		printf ("%5.2f          %10ld       %10ld\n", duration, throughput, best_throughput) ;
Chris@41: 		}
Chris@41: 
Chris@41: 
Chris@41: 	return best_throughput ;
Chris@41: } /* throughput_test */
Chris@41: 
Chris@41: static void
Chris@41: single_run (void)
Chris@41: {
Chris@41: 
Chris@41: 	printf ("\n    CPU name : %s\n", get_cpu_name ()) ;
Chris@41: 
Chris@41: 	puts (
Chris@41: 		"\n"
Chris@41: 		"    Converter                        Duration        Throughput\n"
Chris@41: 		"    -----------------------------------------------------------"
Chris@41: 		) ;
Chris@41: 
Chris@41: 	throughput_test (SRC_ZERO_ORDER_HOLD, 0) ;
Chris@41: 	throughput_test (SRC_LINEAR, 0) ;
Chris@41: 	throughput_test (SRC_SINC_FASTEST, 0) ;
Chris@41: 	throughput_test (SRC_SINC_MEDIUM_QUALITY, 0) ;
Chris@41: 	throughput_test (SRC_SINC_BEST_QUALITY, 0) ;
Chris@41: 
Chris@41: 	puts ("") ;
Chris@41: 	return ;
Chris@41: } /* single_run */
Chris@41: 
Chris@41: static void
Chris@41: multi_run (int run_count)
Chris@41: {	long zero_order_hold = 0, linear = 0 ;
Chris@41: 	long sinc_fastest = 0, sinc_medium = 0, sinc_best = 0 ;
Chris@41: 	int k ;
Chris@41: 
Chris@41: 	puts (
Chris@41: 		"\n"
Chris@41: 		"    Converter                        Duration        Throughput      Best Throughput\n"
Chris@41: 		"    --------------------------------------------------------------------------------"
Chris@41: 		) ;
Chris@41: 
Chris@41: 	for (k = 0 ; k < run_count ; k++)
Chris@41: 	{	zero_order_hold =	throughput_test (SRC_ZERO_ORDER_HOLD, zero_order_hold) ;
Chris@41: 		linear =			throughput_test (SRC_LINEAR, linear) ;
Chris@41: 		sinc_fastest =		throughput_test (SRC_SINC_FASTEST, sinc_fastest) ;
Chris@41: 		sinc_medium =		throughput_test (SRC_SINC_MEDIUM_QUALITY, sinc_medium) ;
Chris@41: 		sinc_best =			throughput_test (SRC_SINC_BEST_QUALITY, sinc_best) ;
Chris@41: 
Chris@41: 		puts ("") ;
Chris@41: 
Chris@41: 		/* Let the CPU cool down. We might be running on a laptop. */
Chris@41: 		sleep (10) ;
Chris@41: 		} ;
Chris@41: 
Chris@41: 	printf ("\n    CPU name : %s\n", get_cpu_name ()) ;
Chris@41: 
Chris@41: 	puts (
Chris@41: 		"\n"
Chris@41: 		"    Converter                        Best Throughput\n"
Chris@41: 		"    ------------------------------------------------"
Chris@41: 		) ;
Chris@41: 	printf ("    %-30s    %10ld\n", src_get_name (SRC_ZERO_ORDER_HOLD), zero_order_hold) ;
Chris@41: 	printf ("    %-30s    %10ld\n", src_get_name (SRC_LINEAR), linear) ;
Chris@41: 	printf ("    %-30s    %10ld\n", src_get_name (SRC_SINC_FASTEST), sinc_fastest) ;
Chris@41: 	printf ("    %-30s    %10ld\n", src_get_name (SRC_SINC_MEDIUM_QUALITY), sinc_medium) ;
Chris@41: 	printf ("    %-30s    %10ld\n", src_get_name (SRC_SINC_BEST_QUALITY), sinc_best) ;
Chris@41: 
Chris@41: 	puts ("") ;
Chris@41: } /* multi_run */
Chris@41: 
Chris@41: static void
Chris@41: usage_exit (const char * argv0)
Chris@41: {	const char * cptr ;
Chris@41: 
Chris@41: 	if ((cptr = strrchr (argv0, '/')) != NULL)
Chris@41: 		argv0 = cptr ;
Chris@41: 
Chris@41: 	printf (
Chris@41: 		"Usage :\n"
Chris@41: 	 	"    %s                 - Single run of the throughput test.\n"
Chris@41: 		"    %s --best-of N     - Do N runs of test a print bext result.\n"
Chris@41: 		"\n",
Chris@41: 		argv0, argv0) ;
Chris@41: 
Chris@41: 	exit (0) ;
Chris@41: } /* usage_exit */
Chris@41: 
Chris@41: int
Chris@41: main (int argc, char ** argv)
Chris@41: {	double freq ;
Chris@41: 
Chris@41: 	memset (input, 0, sizeof (input)) ;
Chris@41: 	freq = 0.01 ;
Chris@41: 	gen_windowed_sines (1, &freq, 1.0, input, BUFFER_LEN) ;
Chris@41: 
Chris@41: 	if (argc == 1)
Chris@41: 		single_run () ;
Chris@41: 	else if (argc == 3 && strcmp (argv [1], "--best-of") == 0)
Chris@41: 	{	int run_count = atoi (argv [2]) ;
Chris@41: 
Chris@41: 		if (run_count < 1 || run_count > 20)
Chris@41: 		{	printf ("Please be sensible. Run count should be in range (1, 10].\n") ;
Chris@41: 			exit (1) ;
Chris@41: 			} ;
Chris@41: 
Chris@41: 		multi_run (run_count) ;
Chris@41: 		}
Chris@41: 	else
Chris@41: 		usage_exit (argv [0]) ;
Chris@41: 
Chris@41: 	puts (
Chris@41: 		"            Duration is in seconds.\n"
Chris@41: 		"            Throughput is in samples/sec (more is better).\n"
Chris@41: 		) ;
Chris@41: 
Chris@41: 	return 0 ;
Chris@41: } /* main */
Chris@41: