cannam@126: /*
cannam@126: ** Copyright (c) 2008-2016, Erik de Castro Lopo <erikd@mega-nerd.com>
cannam@126: ** All rights reserved.
cannam@126: **
cannam@126: ** This code is released under 2-clause BSD license. Please see the
cannam@126: ** file at : https://github.com/erikd/libsamplerate/blob/master/COPYING
cannam@126: */
cannam@126: 
cannam@126: #include <stdio.h>
cannam@126: #include <stdlib.h>
cannam@126: #include <string.h>
cannam@126: #include <time.h>
cannam@126: #include <unistd.h>
cannam@126: 
cannam@126: #include <samplerate.h>
cannam@126: 
cannam@126: #include "config.h"
cannam@126: 
cannam@126: #include "util.h"
cannam@126: #include "float_cast.h"
cannam@126: 
cannam@126: #define BUFFER_LEN	(1<<17)
cannam@126: 
cannam@126: static float input [BUFFER_LEN] ;
cannam@126: static float output [BUFFER_LEN] ;
cannam@126: 
cannam@126: static long
cannam@126: throughput_test (int converter, int channels, long best_throughput)
cannam@126: {	SRC_DATA src_data ;
cannam@126: 	clock_t start_time, clock_time ;
cannam@126: 	double duration ;
cannam@126: 	long total_frames = 0, throughput ;
cannam@126: 	int error ;
cannam@126: 
cannam@126: 	printf ("    %-30s     %2d         ", src_get_name (converter), channels) ;
cannam@126: 	fflush (stdout) ;
cannam@126: 
cannam@126: 	src_data.data_in = input ;
cannam@126: 	src_data.input_frames = ARRAY_LEN (input) / channels ;
cannam@126: 
cannam@126: 	src_data.data_out = output ;
cannam@126: 	src_data.output_frames = ARRAY_LEN (output) / channels ;
cannam@126: 
cannam@126: 	src_data.src_ratio = 0.99 ;
cannam@126: 
cannam@126: 	sleep (2) ;
cannam@126: 
cannam@126: 	start_time = clock () ;
cannam@126: 
cannam@126: 	do
cannam@126: 	{
cannam@126: 		if ((error = src_simple (&src_data, converter, channels)) != 0)
cannam@126: 		{	puts (src_strerror (error)) ;
cannam@126: 			exit (1) ;
cannam@126: 			} ;
cannam@126: 
cannam@126: 		total_frames += src_data.output_frames_gen ;
cannam@126: 
cannam@126: 		clock_time = clock () - start_time ;
cannam@126: 		duration = (1.0 * clock_time) / CLOCKS_PER_SEC ;
cannam@126: 	}
cannam@126: 	while (duration < 5.0) ;
cannam@126: 
cannam@126: 	if (src_data.input_frames_used != src_data.input_frames)
cannam@126: 	{	printf ("\n\nLine %d : input frames used %ld should be %ld\n", __LINE__, src_data.input_frames_used, src_data.input_frames) ;
cannam@126: 		exit (1) ;
cannam@126: 		} ;
cannam@126: 
cannam@126: 	if (fabs (src_data.src_ratio * src_data.input_frames_used - src_data.output_frames_gen) > 2)
cannam@126: 	{	printf ("\n\nLine %d : input / output length mismatch.\n\n", __LINE__) ;
cannam@126: 		printf ("    input len  : %d\n", ARRAY_LEN (input) / channels) ;
cannam@126: 		printf ("    output len : %ld (should be %g +/- 2)\n\n", src_data.output_frames_gen,
cannam@126: 				floor (0.5 + src_data.src_ratio * src_data.input_frames_used)) ;
cannam@126: 		exit (1) ;
cannam@126: 		} ;
cannam@126: 
cannam@126: 	throughput = lrint (floor (total_frames / duration)) ;
cannam@126: 
cannam@126: 	if (best_throughput == 0)
cannam@126: 	{	best_throughput = MAX (throughput, best_throughput) ;
cannam@126: 		printf ("%5.2f      %10ld\n", duration, throughput) ;
cannam@126: 		}
cannam@126: 	else
cannam@126: 	{	best_throughput = MAX (throughput, best_throughput) ;
cannam@126: 		printf ("%5.2f      %10ld       %10ld\n", duration, throughput, best_throughput) ;
cannam@126: 		}
cannam@126: 
cannam@126: 	return best_throughput ;
cannam@126: } /* throughput_test */
cannam@126: 
cannam@126: static void
cannam@126: single_run (void)
cannam@126: {	const int max_channels = 10 ;
cannam@126: 	int k ;
cannam@126: 
cannam@126: 	printf ("\n    CPU name : %s\n", get_cpu_name ()) ;
cannam@126: 
cannam@126: 	puts (
cannam@126: 		"\n"
cannam@126: 		"    Converter                        Channels    Duration      Throughput\n"
cannam@126: 		"    ---------------------------------------------------------------------"
cannam@126: 		) ;
cannam@126: 
cannam@126: 	for (k = 1 ; k <= max_channels / 2 ; k++)
cannam@126: 		throughput_test (SRC_SINC_FASTEST, k, 0) ;
cannam@126: 
cannam@126: 	puts ("") ;
cannam@126: 	for (k = 1 ; k <= max_channels / 2 ; k++)
cannam@126: 		throughput_test (SRC_SINC_MEDIUM_QUALITY, k, 0) ;
cannam@126: 
cannam@126: 	puts ("") ;
cannam@126: 	for (k = 1 ; k <= max_channels ; k++)
cannam@126: 		throughput_test (SRC_SINC_BEST_QUALITY, k, 0) ;
cannam@126: 
cannam@126: 	puts ("") ;
cannam@126: 	return ;
cannam@126: } /* single_run */
cannam@126: 
cannam@126: static void
cannam@126: multi_run (int run_count)
cannam@126: {	int k, ch ;
cannam@126: 
cannam@126: 	printf ("\n    CPU name : %s\n", get_cpu_name ()) ;
cannam@126: 
cannam@126: 	puts (
cannam@126: 		"\n"
cannam@126: 		"    Converter                        Channels    Duration      Throughput    Best Throughput\n"
cannam@126: 		"    ----------------------------------------------------------------------------------------"
cannam@126: 		) ;
cannam@126: 
cannam@126: 	for (ch = 1 ; ch <= 5 ; ch++)
cannam@126: 	{	long sinc_fastest = 0, sinc_medium = 0, sinc_best = 0 ;
cannam@126: 
cannam@126: 		for (k = 0 ; k < run_count ; k++)
cannam@126: 		{	sinc_fastest =		throughput_test (SRC_SINC_FASTEST, ch, sinc_fastest) ;
cannam@126: 			sinc_medium =		throughput_test (SRC_SINC_MEDIUM_QUALITY, ch, sinc_medium) ;
cannam@126: 			sinc_best =			throughput_test (SRC_SINC_BEST_QUALITY, ch, sinc_best) ;
cannam@126: 
cannam@126: 			puts ("") ;
cannam@126: 
cannam@126: 			/* Let the CPU cool down. We might be running on a laptop. */
cannam@126: 			sleep (10) ;
cannam@126: 			} ;
cannam@126: 
cannam@126: 		puts (
cannam@126: 			"\n"
cannam@126: 			"    Converter                        Best Throughput\n"
cannam@126: 			"    ------------------------------------------------"
cannam@126: 			) ;
cannam@126: 
cannam@126: 		printf ("    %-30s    %10ld\n", src_get_name (SRC_SINC_FASTEST), sinc_fastest) ;
cannam@126: 		printf ("    %-30s    %10ld\n", src_get_name (SRC_SINC_MEDIUM_QUALITY), sinc_medium) ;
cannam@126: 		printf ("    %-30s    %10ld\n", src_get_name (SRC_SINC_BEST_QUALITY), sinc_best) ;
cannam@126: 		} ;
cannam@126: 
cannam@126: 	puts ("") ;
cannam@126: } /* multi_run */
cannam@126: 
cannam@126: static void
cannam@126: usage_exit (const char * argv0)
cannam@126: {	const char * cptr ;
cannam@126: 
cannam@126: 	if ((cptr = strrchr (argv0, '/')) != NULL)
cannam@126: 		argv0 = cptr ;
cannam@126: 
cannam@126: 	printf (
cannam@126: 		"Usage :\n"
cannam@126: 	 	"    %s                 - Single run of the throughput test.\n"
cannam@126: 		"    %s --best-of N     - Do N runs of test a print bext result.\n"
cannam@126: 		"\n",
cannam@126: 		argv0, argv0) ;
cannam@126: 
cannam@126: 	exit (0) ;
cannam@126: } /* usage_exit */
cannam@126: 
cannam@126: int
cannam@126: main (int argc, char ** argv)
cannam@126: {	double freq ;
cannam@126: 
cannam@126: 	memset (input, 0, sizeof (input)) ;
cannam@126: 	freq = 0.01 ;
cannam@126: 	gen_windowed_sines (1, &freq, 1.0, input, BUFFER_LEN) ;
cannam@126: 
cannam@126: 	if (argc == 1)
cannam@126: 		single_run () ;
cannam@126: 	else if (argc == 3 && strcmp (argv [1], "--best-of") == 0)
cannam@126: 	{	int run_count = atoi (argv [2]) ;
cannam@126: 
cannam@126: 		if (run_count < 1 || run_count > 20)
cannam@126: 		{	printf ("Please be sensible. Run count should be in range (1, 10].\n") ;
cannam@126: 			exit (1) ;
cannam@126: 			} ;
cannam@126: 
cannam@126: 		multi_run (run_count) ;
cannam@126: 		}
cannam@126: 	else
cannam@126: 		usage_exit (argv [0]) ;
cannam@126: 
cannam@126: 	puts (
cannam@126: 		"            Duration is in seconds.\n"
cannam@126: 		"            Throughput is in frames/sec (more is better).\n"
cannam@126: 		) ;
cannam@126: 
cannam@126: 	return 0 ;
cannam@126: } /* main */
cannam@126: