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 <stdio.h>
Chris@41: #include <stdlib.h>
Chris@41: #include <string.h>
Chris@41: #include <ctype.h>
Chris@41: #include <math.h>
Chris@41: 
Chris@41: #include	"util.h"
Chris@41: 
Chris@41: #ifndef	M_PI
Chris@41: #define	M_PI			3.14159265358979323846264338
Chris@41: #endif
Chris@41: 
Chris@41: void
Chris@41: gen_windowed_sines (int freq_count, const double *freqs, double max, float *output, int output_len)
Chris@41: {	int 	k, freq ;
Chris@41: 	double	amplitude, phase ;
Chris@41: 
Chris@41: 	amplitude = max / freq_count ;
Chris@41: 
Chris@41: 	for (k = 0 ; k < output_len ; k++)
Chris@41: 		output [k] = 0.0 ;
Chris@41: 
Chris@41: 	for (freq = 0 ; freq < freq_count ; freq++)
Chris@41: 	{	phase = 0.9 * M_PI / freq_count ;
Chris@41: 
Chris@41: 		if (freqs [freq] <= 0.0 || freqs [freq] >= 0.5)
Chris@41: 		{	printf ("\n%s : Error : freq [%d] == %g is out of range. Should be < 0.5.\n", __FILE__, freq, freqs [freq]) ;
Chris@41: 			exit (1) ;
Chris@41: 			} ;
Chris@41: 
Chris@41: 		for (k = 0 ; k < output_len ; k++)
Chris@41: 			output [k] += amplitude * sin (freqs [freq] * (2 * k) * M_PI + phase) ;
Chris@41: 		} ;
Chris@41: 
Chris@41: 	/* Apply Hanning Window. */
Chris@41: 	for (k = 0 ; k < output_len ; k++)
Chris@41: 		output [k] *= 0.5 - 0.5 * cos ((2 * k) * M_PI / (output_len - 1)) ;
Chris@41: 
Chris@41: 	/*	data [k] *= 0.3635819 - 0.4891775 * cos ((2 * k) * M_PI / (output_len - 1))
Chris@41: 					+ 0.1365995 * cos ((4 * k) * M_PI / (output_len - 1))
Chris@41: 					- 0.0106411 * cos ((6 * k) * M_PI / (output_len - 1)) ;
Chris@41: 		*/
Chris@41: 
Chris@41: 	return ;
Chris@41: } /* gen_windowed_sines */
Chris@41: 
Chris@41: void
Chris@41: save_oct_float (char *filename, float *input, int in_len, float *output, int out_len)
Chris@41: {	FILE 	*file ;
Chris@41: 	int		k ;
Chris@41: 
Chris@41: 	printf ("Dumping input and output data to file : %s.\n\n", filename) ;
Chris@41: 
Chris@41: 	if (! (file = fopen (filename, "w")))
Chris@41: 		return ;
Chris@41: 
Chris@41: 	fprintf (file, "# Not created by Octave\n") ;
Chris@41: 
Chris@41: 	fprintf (file, "# name: input\n") ;
Chris@41: 	fprintf (file, "# type: matrix\n") ;
Chris@41: 	fprintf (file, "# rows: %d\n", in_len) ;
Chris@41: 	fprintf (file, "# columns: 1\n") ;
Chris@41: 
Chris@41: 	for (k = 0 ; k < in_len ; k++)
Chris@41: 		fprintf (file, "% g\n", input [k]) ;
Chris@41: 
Chris@41: 	fprintf (file, "# name: output\n") ;
Chris@41: 	fprintf (file, "# type: matrix\n") ;
Chris@41: 	fprintf (file, "# rows: %d\n", out_len) ;
Chris@41: 	fprintf (file, "# columns: 1\n") ;
Chris@41: 
Chris@41: 	for (k = 0 ; k < out_len ; k++)
Chris@41: 		fprintf (file, "% g\n", output [k]) ;
Chris@41: 
Chris@41: 	fclose (file) ;
Chris@41: 	return ;
Chris@41: } /* save_oct_float */
Chris@41: 
Chris@41: void
Chris@41: save_oct_double (char *filename, double *input, int in_len, double *output, int out_len)
Chris@41: {	FILE 	*file ;
Chris@41: 	int		k ;
Chris@41: 
Chris@41: 	printf ("Dumping input and output data to file : %s.\n\n", filename) ;
Chris@41: 
Chris@41: 	if (! (file = fopen (filename, "w")))
Chris@41: 		return ;
Chris@41: 
Chris@41: 	fprintf (file, "# Not created by Octave\n") ;
Chris@41: 
Chris@41: 	fprintf (file, "# name: input\n") ;
Chris@41: 	fprintf (file, "# type: matrix\n") ;
Chris@41: 	fprintf (file, "# rows: %d\n", in_len) ;
Chris@41: 	fprintf (file, "# columns: 1\n") ;
Chris@41: 
Chris@41: 	for (k = 0 ; k < in_len ; k++)
Chris@41: 		fprintf (file, "% g\n", input [k]) ;
Chris@41: 
Chris@41: 	fprintf (file, "# name: output\n") ;
Chris@41: 	fprintf (file, "# type: matrix\n") ;
Chris@41: 	fprintf (file, "# rows: %d\n", out_len) ;
Chris@41: 	fprintf (file, "# columns: 1\n") ;
Chris@41: 
Chris@41: 	for (k = 0 ; k < out_len ; k++)
Chris@41: 		fprintf (file, "% g\n", output [k]) ;
Chris@41: 
Chris@41: 	fclose (file) ;
Chris@41: 	return ;
Chris@41: } /* save_oct_double */
Chris@41: 
Chris@41: void
Chris@41: interleave_data (const float *in, float *out, int frames, int channels)
Chris@41: {	int fr, ch ;
Chris@41: 
Chris@41: 	for (fr = 0 ; fr < frames ; fr++)
Chris@41: 		for (ch = 0 ; ch < channels ; ch++)
Chris@41: 			out [ch + channels * fr] = in [fr + frames * ch] ;
Chris@41: 
Chris@41: 	return ;
Chris@41: } /* interleave_data */
Chris@41: 
Chris@41: void
Chris@41: deinterleave_data (const float *in, float *out, int frames, int channels)
Chris@41: {	int fr, ch ;
Chris@41: 
Chris@41: 	for (ch = 0 ; ch < channels ; ch++)
Chris@41: 		for (fr = 0 ; fr < frames ; fr++)
Chris@41: 			out [fr + frames * ch] = in [ch + channels * fr] ;
Chris@41: 
Chris@41: 	return ;
Chris@41: } /* deinterleave_data */
Chris@41: 
Chris@41: void
Chris@41: reverse_data (float *data, int datalen)
Chris@41: {	int left, right ;
Chris@41: 	float temp ;
Chris@41: 
Chris@41: 	left = 0 ;
Chris@41: 	right = datalen - 1 ;
Chris@41: 
Chris@41: 	while (left < right)
Chris@41: 	{	temp = data [left] ;
Chris@41: 		data [left] = data [right] ;
Chris@41: 		data [right] = temp ;
Chris@41: 		left ++ ;
Chris@41: 		right -- ;
Chris@41: 		} ;
Chris@41: 
Chris@41: } /* reverse_data */
Chris@41: 
Chris@41: const char *
Chris@41: get_cpu_name (void)
Chris@41: {
Chris@41: 	const char *name = "Unknown", *search = NULL ;
Chris@41: 	static char buffer [512] ;
Chris@41: 	FILE * file = NULL ;
Chris@41: 	int is_pipe = 0 ;
Chris@41: 
Chris@41: #if defined (__linux__)
Chris@41: 	file = fopen ("/proc/cpuinfo", "r") ;
Chris@41: 	search = "model name" ;
Chris@41: #elif defined (__APPLE__)
Chris@41: 	file = popen ("/usr/sbin/system_profiler -detailLevel full SPHardwareDataType", "r") ;
Chris@41: 	search = "Processor Name" ;
Chris@41: 	is_pipe = 1 ;
Chris@41: #elif defined (__FreeBSD__)
Chris@41: 	file = popen ("sysctl -a", "r") ;
Chris@41: 	search = "hw.model" ;
Chris@41: 	is_pipe = 1 ;
Chris@41: #else
Chris@41: 	file = NULL ;
Chris@41: #endif
Chris@41: 
Chris@41: 	if (file == NULL)
Chris@41: 		return name ;
Chris@41: 
Chris@41: 	if (search == NULL)
Chris@41: 	{	printf ("Error : search is NULL in function %s.\n", __func__) ;
Chris@41: 		return name ;
Chris@41: 		} ;
Chris@41: 
Chris@41: 	while (fgets (buffer, sizeof (buffer), file) != NULL)
Chris@41: 		if (strstr (buffer, search))
Chris@41: 		{	char *src, *dest ;
Chris@41: 
Chris@41: 			if ((src = strchr (buffer, ':')) != NULL)
Chris@41: 			{	src ++ ;
Chris@41: 				while (isspace (src [0]))
Chris@41: 					src ++ ;
Chris@41: 				name = src ;
Chris@41: 
Chris@41: 				/* Remove consecutive spaces. */
Chris@41: 				src ++ ;
Chris@41: 				for (dest = src ; src [0] ; src ++)
Chris@41: 				{	if (isspace (src [0]) && isspace (dest [-1]))
Chris@41: 						continue ;
Chris@41: 					dest [0] = src [0] ;
Chris@41: 					dest ++ ;
Chris@41: 					} ;
Chris@41: 				dest [0] = 0 ;
Chris@41: 				break ;
Chris@41: 				} ;
Chris@41: 			} ;
Chris@41: 
Chris@41: 	if (is_pipe)
Chris@41: 		pclose (file) ;
Chris@41: 	else
Chris@41: 		fclose (file) ;
Chris@41: 
Chris@41: 	return name ;
Chris@41: } /* get_cpu_name */
Chris@41: