ref: 183610fcad977f4fead56178cf731a2e1c581e65
dir: /tests/util.c/
/* ** Copyright (c) 2002-2016, Erik de Castro Lopo <[email protected]> ** 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 void gen_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 */ void save_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 */ void save_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 */ void interleave_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 */ void deinterleave_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 */ void reverse_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 */