cannam@85: /*
cannam@85: ** Copyright (C) 2008-2011 Erik de Castro Lopo <erikd@mega-nerd.com>
cannam@85: **
cannam@85: ** This program is free software; you can redistribute it and/or modify
cannam@85: ** it under the terms of the GNU General Public License as published by
cannam@85: ** the Free Software Foundation; either version 2 of the License, or
cannam@85: ** (at your option) any later version.
cannam@85: **
cannam@85: ** This program is distributed in the hope that it will be useful,
cannam@85: ** but WITHOUT ANY WARRANTY; without even the implied warranty of
cannam@85: ** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
cannam@85: ** GNU General Public License for more details.
cannam@85: **
cannam@85: ** You should have received a copy of the GNU General Public License
cannam@85: ** along with this program; if not, write to the Free Software
cannam@85: ** Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307, USA.
cannam@85: */
cannam@85: 
cannam@85: #include <stdio.h>
cannam@85: #include <stdlib.h>
cannam@85: #include <string.h>
cannam@85: #include <time.h>
cannam@85: #include <unistd.h>
cannam@85: 
cannam@85: #include <samplerate.h>
cannam@85: 
cannam@85: #include "config.h"
cannam@85: 
cannam@85: #include "util.h"
cannam@85: #include "float_cast.h"
cannam@85: 
cannam@85: #define BUFFER_LEN	(1<<17)
cannam@85: 
cannam@85: static float input [BUFFER_LEN] ;
cannam@85: static float output [BUFFER_LEN] ;
cannam@85: 
cannam@85: static long
cannam@85: throughput_test (int converter, int channels, long best_throughput)
cannam@85: {	SRC_DATA src_data ;
cannam@85: 	clock_t start_time, clock_time ;
cannam@85: 	double duration ;
cannam@85: 	long total_frames = 0, throughput ;
cannam@85: 	int error ;
cannam@85: 
cannam@85: 	printf ("    %-30s     %2d         ", src_get_name (converter), channels) ;
cannam@85: 	fflush (stdout) ;
cannam@85: 
cannam@85: 	src_data.data_in = input ;
cannam@85: 	src_data.input_frames = ARRAY_LEN (input) / channels ;
cannam@85: 
cannam@85: 	src_data.data_out = output ;
cannam@85: 	src_data.output_frames = ARRAY_LEN (output) / channels ;
cannam@85: 
cannam@85: 	src_data.src_ratio = 0.99 ;
cannam@85: 
cannam@85: 	sleep (2) ;
cannam@85: 
cannam@85: 	start_time = clock () ;
cannam@85: 
cannam@85: 	do
cannam@85: 	{
cannam@85: 		if ((error = src_simple (&src_data, converter, channels)) != 0)
cannam@85: 		{	puts (src_strerror (error)) ;
cannam@85: 			exit (1) ;
cannam@85: 			} ;
cannam@85: 
cannam@85: 		total_frames += src_data.output_frames_gen ;
cannam@85: 
cannam@85: 		clock_time = clock () - start_time ;
cannam@85: 		duration = (1.0 * clock_time) / CLOCKS_PER_SEC ;
cannam@85: 	}
cannam@85: 	while (duration < 5.0) ;
cannam@85: 
cannam@85: 	if (src_data.input_frames_used != src_data.input_frames)
cannam@85: 	{	printf ("\n\nLine %d : input frames used %ld should be %ld\n", __LINE__, src_data.input_frames_used, src_data.input_frames) ;
cannam@85: 		exit (1) ;
cannam@85: 		} ;
cannam@85: 
cannam@85: 	if (fabs (src_data.src_ratio * src_data.input_frames_used - src_data.output_frames_gen) > 2)
cannam@85: 	{	printf ("\n\nLine %d : input / output length mismatch.\n\n", __LINE__) ;
cannam@85: 		printf ("    input len  : %d\n", ARRAY_LEN (input) / channels) ;
cannam@85: 		printf ("    output len : %ld (should be %g +/- 2)\n\n", src_data.output_frames_gen,
cannam@85: 				floor (0.5 + src_data.src_ratio * src_data.input_frames_used)) ;
cannam@85: 		exit (1) ;
cannam@85: 		} ;
cannam@85: 
cannam@85: 	throughput = lrint (floor (total_frames / duration)) ;
cannam@85: 
cannam@85: 	if (best_throughput == 0)
cannam@85: 	{	best_throughput = MAX (throughput, best_throughput) ;
cannam@85: 		printf ("%5.2f      %10ld\n", duration, throughput) ;
cannam@85: 		}
cannam@85: 	else
cannam@85: 	{	best_throughput = MAX (throughput, best_throughput) ;
cannam@85: 		printf ("%5.2f      %10ld       %10ld\n", duration, throughput, best_throughput) ;
cannam@85: 		}
cannam@85: 
cannam@85: 	return best_throughput ;
cannam@85: } /* throughput_test */
cannam@85: 
cannam@85: static void
cannam@85: single_run (void)
cannam@85: {	const int max_channels = 10 ;
cannam@85: 	int k ;
cannam@85: 
cannam@85: 	printf ("\n    CPU name : %s\n", get_cpu_name ()) ;
cannam@85: 
cannam@85: 	puts (
cannam@85: 		"\n"
cannam@85: 		"    Converter                        Channels    Duration      Throughput\n"
cannam@85: 		"    ---------------------------------------------------------------------"
cannam@85: 		) ;
cannam@85: 
cannam@85: 	for (k = 1 ; k <= max_channels / 2 ; k++)
cannam@85: 		throughput_test (SRC_SINC_FASTEST, k, 0) ;
cannam@85: 
cannam@85: 	puts ("") ;
cannam@85: 	for (k = 1 ; k <= max_channels / 2 ; k++)
cannam@85: 		throughput_test (SRC_SINC_MEDIUM_QUALITY, k, 0) ;
cannam@85: 
cannam@85: 	puts ("") ;
cannam@85: 	for (k = 1 ; k <= max_channels ; k++)
cannam@85: 		throughput_test (SRC_SINC_BEST_QUALITY, k, 0) ;
cannam@85: 
cannam@85: 	puts ("") ;
cannam@85: 	return ;
cannam@85: } /* single_run */
cannam@85: 
cannam@85: static void
cannam@85: multi_run (int run_count)
cannam@85: {	int k, ch ;
cannam@85: 
cannam@85: 	printf ("\n    CPU name : %s\n", get_cpu_name ()) ;
cannam@85: 
cannam@85: 	puts (
cannam@85: 		"\n"
cannam@85: 		"    Converter                        Channels    Duration      Throughput    Best Throughput\n"
cannam@85: 		"    ----------------------------------------------------------------------------------------"
cannam@85: 		) ;
cannam@85: 
cannam@85: 	for (ch = 1 ; ch <= 5 ; ch++)
cannam@85: 	{	long sinc_fastest = 0, sinc_medium = 0, sinc_best = 0 ;
cannam@85: 
cannam@85: 		for (k = 0 ; k < run_count ; k++)
cannam@85: 		{	sinc_fastest =		throughput_test (SRC_SINC_FASTEST, ch, sinc_fastest) ;
cannam@85: 			sinc_medium =		throughput_test (SRC_SINC_MEDIUM_QUALITY, ch, sinc_medium) ;
cannam@85: 			sinc_best =			throughput_test (SRC_SINC_BEST_QUALITY, ch, sinc_best) ;
cannam@85: 
cannam@85: 			puts ("") ;
cannam@85: 
cannam@85: 			/* Let the CPU cool down. We might be running on a laptop. */
cannam@85: 			sleep (10) ;
cannam@85: 			} ;
cannam@85: 
cannam@85: 		puts (
cannam@85: 			"\n"
cannam@85: 			"    Converter                        Best Throughput\n"
cannam@85: 			"    ------------------------------------------------"
cannam@85: 			) ;
cannam@85: 
cannam@85: 		printf ("    %-30s    %10ld\n", src_get_name (SRC_SINC_FASTEST), sinc_fastest) ;
cannam@85: 		printf ("    %-30s    %10ld\n", src_get_name (SRC_SINC_MEDIUM_QUALITY), sinc_medium) ;
cannam@85: 		printf ("    %-30s    %10ld\n", src_get_name (SRC_SINC_BEST_QUALITY), sinc_best) ;
cannam@85: 		} ;
cannam@85: 
cannam@85: 	puts ("") ;
cannam@85: } /* multi_run */
cannam@85: 
cannam@85: static void
cannam@85: usage_exit (const char * argv0)
cannam@85: {	const char * cptr ;
cannam@85: 
cannam@85: 	if ((cptr = strrchr (argv0, '/')) != NULL)
cannam@85: 		argv0 = cptr ;
cannam@85: 
cannam@85: 	printf (
cannam@85: 		"Usage :\n"
cannam@85: 	 	"    %s                 - Single run of the throughput test.\n"
cannam@85: 		"    %s --best-of N     - Do N runs of test a print bext result.\n"
cannam@85: 		"\n",
cannam@85: 		argv0, argv0) ;
cannam@85: 
cannam@85: 	exit (0) ;
cannam@85: } /* usage_exit */
cannam@85: 
cannam@85: int
cannam@85: main (int argc, char ** argv)
cannam@85: {	double freq ;
cannam@85: 
cannam@85: 	memset (input, 0, sizeof (input)) ;
cannam@85: 	freq = 0.01 ;
cannam@85: 	gen_windowed_sines (1, &freq, 1.0, input, BUFFER_LEN) ;
cannam@85: 
cannam@85: 	if (argc == 1)
cannam@85: 		single_run () ;
cannam@85: 	else if (argc == 3 && strcmp (argv [1], "--best-of") == 0)
cannam@85: 	{	int run_count = atoi (argv [2]) ;
cannam@85: 
cannam@85: 		if (run_count < 1 || run_count > 20)
cannam@85: 		{	printf ("Please be sensible. Run count should be in range (1, 10].\n") ;
cannam@85: 			exit (1) ;
cannam@85: 			} ;
cannam@85: 
cannam@85: 		multi_run (run_count) ;
cannam@85: 		}
cannam@85: 	else
cannam@85: 		usage_exit (argv [0]) ;
cannam@85: 
cannam@85: 	puts (
cannam@85: 		"            Duration is in seconds.\n"
cannam@85: 		"            Throughput is in frames/sec (more is better).\n"
cannam@85: 		) ;
cannam@85: 
cannam@85: 	return 0 ;
cannam@85: } /* main */
cannam@85: