/*** Copyright (c) 2002-2016, Erik de Castro Lopo <erikd@mega-nerd.com>** All rights reserved.**** This code is released under 2-clause BSD license. Please see the** file at : https://github.com/libsndfile/libsamplerate/blob/master/COPYING*/ #include <stdio.h>#include <stdlib.h>#include <string.h>#include <ctype.h>#include <math.h> #include	"util.h" #ifndef	M_PI#define	M_PI			3.14159265358979323846264338#endif voidgen_windowed_sines (int freq_count, const double *freqs, double max, float *output, int output_len){	int 	k, freq ;	double	amplitude, phase ; 	amplitude = max / freq_count ; 	for (k = 0 ; k < output_len ; k++)		output [k] = 0.0 ; 	for (freq = 0 ; freq < freq_count ; freq++)	{	phase = 0.9 * M_PI / freq_count ; 		if (freqs [freq] <= 0.0 || freqs [freq] >= 0.5)		{	printf ("\n%s : Error : freq [%d] == %g is out of range. Should be < 0.5.\n", __FILE__, freq, freqs [freq]) ;			exit (1) ;			} ; 		for (k = 0 ; k < output_len ; k++)			output [k] = (float) (output [k] + (amplitude * sin (freqs [freq] * (2 * k) * M_PI + phase))) ;		} ; 	/* Apply Hanning Window. */	for (k = 0 ; k < output_len ; k++)		output [k] = (float) (output [k] * (0.5 - 0.5 * cos ((2 * k) * M_PI / (output_len - 1)))) ; 	/*	data [k] *= 0.3635819 - 0.4891775 * cos ((2 * k) * M_PI / (output_len - 1))					+ 0.1365995 * cos ((4 * k) * M_PI / (output_len - 1))					- 0.0106411 * cos ((6 * k) * M_PI / (output_len - 1)) ;		*/ 	return ;} /* gen_windowed_sines */ voidsave_oct_float (char *filename, float *input, int in_len, float *output, int out_len){	FILE 	*file ;	int		k ; 	printf ("Dumping input and output data to file : %s.\n\n", filename) ; 	if (! (file = fopen (filename, "w")))		return ; 	fprintf (file, "# Not created by Octave\n") ; 	fprintf (file, "# name: input\n") ;	fprintf (file, "# type: matrix\n") ;	fprintf (file, "# rows: %d\n", in_len) ;	fprintf (file, "# columns: 1\n") ; 	for (k = 0 ; k < in_len ; k++)		fprintf (file, "% g\n", input [k]) ; 	fprintf (file, "# name: output\n") ;	fprintf (file, "# type: matrix\n") ;	fprintf (file, "# rows: %d\n", out_len) ;	fprintf (file, "# columns: 1\n") ; 	for (k = 0 ; k < out_len ; k++)		fprintf (file, "% g\n", output [k]) ; 	fclose (file) ;	return ;} /* save_oct_float */ voidsave_oct_double (char *filename, double *input, int in_len, double *output, int out_len){	FILE 	*file ;	int		k ; 	printf ("Dumping input and output data to file : %s.\n\n", filename) ; 	if (! (file = fopen (filename, "w")))		return ; 	fprintf (file, "# Not created by Octave\n") ; 	fprintf (file, "# name: input\n") ;	fprintf (file, "# type: matrix\n") ;	fprintf (file, "# rows: %d\n", in_len) ;	fprintf (file, "# columns: 1\n") ; 	for (k = 0 ; k < in_len ; k++)		fprintf (file, "% g\n", input [k]) ; 	fprintf (file, "# name: output\n") ;	fprintf (file, "# type: matrix\n") ;	fprintf (file, "# rows: %d\n", out_len) ;	fprintf (file, "# columns: 1\n") ; 	for (k = 0 ; k < out_len ; k++)		fprintf (file, "% g\n", output [k]) ; 	fclose (file) ;	return ;} /* save_oct_double */ voidinterleave_data (const float *in, float *out, int frames, int channels){	int fr, ch ; 	for (fr = 0 ; fr < frames ; fr++)		for (ch = 0 ; ch < channels ; ch++)			out [ch + channels * fr] = in [fr + frames * ch] ; 	return ;} /* interleave_data */ voiddeinterleave_data (const float *in, float *out, int frames, int channels){	int fr, ch ; 	for (ch = 0 ; ch < channels ; ch++)		for (fr = 0 ; fr < frames ; fr++)			out [fr + frames * ch] = in [ch + channels * fr] ; 	return ;} /* deinterleave_data */ voidreverse_data (float *data, int datalen){	int left, right ;	float temp ; 	left = 0 ;	right = datalen - 1 ; 	while (left < right)	{	temp = data [left] ;		data [left] = data [right] ;		data [right] = temp ;		left ++ ;		right -- ;		} ; } /* reverse_data */ const char *get_cpu_name (void){	const char *name = "Unknown", *search = NULL ;	static char buffer [512] ;	FILE * file = NULL ;	int is_pipe = 0 ; #if defined (__linux__)	file = fopen ("/proc/cpuinfo", "r") ;	search = "model name" ;#elif defined (__APPLE__)	file = popen ("/usr/sbin/system_profiler -detailLevel full SPHardwareDataType", "r") ;	search = "Processor Name" ;	is_pipe = 1 ;#elif defined (__FreeBSD__)	file = popen ("sysctl -a", "r") ;	search = "hw.model" ;	is_pipe = 1 ;#else	return name;#endif #if defined (__linux__) || defined (__APPLE__) || defined (__FreeBSD__)	if (search == NULL)	{	printf ("Error : search is NULL in function %s.\n", __func__) ;		return name ;		} ; 	while (fgets (buffer, sizeof (buffer), file) != NULL)		if (strstr (buffer, search))		{	char *src, *dest ; 			if ((src = strchr (buffer, ':')) != NULL)			{	src ++ ;				while (isspace (src [0]))					src ++ ;				name = src ; 				/* Remove consecutive spaces. */				src ++ ;				for (dest = src ; src [0] ; src ++)				{	if (isspace (src [0]) && isspace (dest [-1]))						continue ;					dest [0] = src [0] ;					dest ++ ;					} ;				dest [0] = 0 ;				break ;				} ;			} ; 	if (is_pipe)		pclose (file) ;	else		fclose (file) ; 	return name ;#endif} /* get_cpu_name */