int main(int argc, char** argv)
{
set_affinity();
printf("FFT/DFT benchmarking tool\n");
printf("Copyright (C) 2016 D Levin https://www.kfrlib.com\n");
printf("Benchmark source code is dual-licensed under MIT and GPL 2 or later\n");
printf("Individual DFT/FFT algorithms have its own licenses. See its source code for details\n");
printf("Usage:\n %s <size> <size> ... <size>\n\n", execfile(argv).c_str());
printf("Data type: %s\n", real_name);
printf("Algorithm: %s\n", fft_benchmark<false>::name().c_str());
printf("Built using %s\n\n", __VERSION__);
printf(">[\n");
for (size_t i = 1; i < argc; i++)
{
const long size = std::atol(argv[i]);
if (size > 1)
test_fft(size);
else
printf("Incorrect size: %s\n", argv[i]);
}
printf(">]\n");
return 0;
}
int main(int argc, char **argv) {
test_fft();
}
int main (void)
{
size_t exponent = 4;
test_fft(exponent);
}
int main(int argc, char *argv[])
{
log_file("Starting arss ...");
test_fft();
FILE *fin;
FILE *fout;
int32_t i;
double **sound, **image, basefreq=0, maxfreq=0, pixpersec=0, bpo=0, brightness=1;
int32_t channels, samplecount=0, samplerate=0, Xsize=0, Ysize=0, format_param=0;
int32_t clockb;
char mode=0, *in_name=NULL, *out_name=NULL;
// initialisation of global using defaults defined in dsp.h
pi=PI_D;
LOGBASE=LOGBASE_D;
LOOP_SIZE_SEC=LOOP_SIZE_SEC_D;
BMSQ_LUT_SIZE=BMSQ_LUT_SIZE_D;
#ifdef QUIET
quiet=1;
#else
quiet=0;
#endif
printf("The Analysis & Resynthesis Sound Spectrograph %s\n", version);
srand(time(NULL));
for (i=1; i<argc; i++)
{
if (strcmp(argv[i], "/?")==0) // DOS friendly help
{
argv[i][0] = '-';
argv[i][1] = 'h';
}
if (argv[i][0] != '-') // if the argument is not a function
{
if (in_name==NULL) // if the input file name hasn't been filled in yet
in_name = argv[i];
else
if (out_name==NULL) // if the input name has been filled but not the output name yet
out_name = argv[i];
else // if both names have already been filled in
{
fprintf(stderr, "You can only have two file names as parameters.\nRemove parameter \"%s\".\nExiting with error.\n", argv[i]);
exit(EXIT_FAILURE);
}
}
else // if the argument is a parameter
{
if (strcmp(argv[i], "--analysis")==0 || strcmp(argv[i], "-a")==0)
mode=1;
if (strcmp(argv[i], "--sine")==0 || strcmp(argv[i], "-s")==0)
mode=2;
if (strcmp(argv[i], "--noise")==0 || strcmp(argv[i], "-n")==0)
mode=3;
if (strcmp(argv[i], "--quiet")==0 || strcmp(argv[i], "-q")==0)
quiet=1;
if (strcmp(argv[i], "--linear")==0 || strcmp(argv[i], "-l")==0)
LOGBASE=1.0;
if (strcmp(argv[i], "--sample-rate")==0 || strcmp(argv[i], "-r")==0)
if (str_isnumber(argv[++i]))
samplerate = atoi(argv[i]);
else
{
fprintf(stderr, MSG_NUMBER_EXPECTED, argv[i-1]);
exit(EXIT_FAILURE);
}
if (strcmp(argv[i], "--min-freq")==0 || strcmp(argv[i], "-min")==0)
if (str_isnumber(argv[++i]))
{
basefreq = atof(argv[i]);
if (basefreq==0)
basefreq = DBL_MIN; // Set it to this extremely close-to-zero number so that it's considered set
}
else
{
fprintf(stderr, MSG_NUMBER_EXPECTED, argv[i-1]);
exit(EXIT_FAILURE);
}
if (strcmp(argv[i], "--max-freq")==0 || strcmp(argv[i], "-max")==0)
if (str_isnumber(argv[++i]))
maxfreq = atof(argv[i]);
else
{
fprintf(stderr, MSG_NUMBER_EXPECTED, argv[i-1]);
exit(EXIT_FAILURE);
}
if (strcmp(argv[i], "--bpo")==0 || strcmp(argv[i], "-b")==0)
if (str_isnumber(argv[++i]))
bpo = atof(argv[i]);
else
{
fprintf(stderr, MSG_NUMBER_EXPECTED, argv[i-1]);
exit(EXIT_FAILURE);
}
if (strcmp(argv[i], "--pps")==0 || strcmp(argv[i], "-p")==0)
if (str_isnumber(argv[++i]))
pixpersec = atof(argv[i]);
else
{
fprintf(stderr, MSG_NUMBER_EXPECTED, argv[i-1]);
exit(EXIT_FAILURE);
}
if (strcmp(argv[i], "--height")==0 || strcmp(argv[i], "-y")==0)
if (str_isnumber(argv[++i]))
Ysize = atoi(argv[i]);
else
{
fprintf(stderr, MSG_NUMBER_EXPECTED, argv[i-1]);
exit(EXIT_FAILURE);
}
if (strcmp(argv[i], "--width")==0 || strcmp(argv[i], "-x")==0)
if (str_isnumber(argv[++i]))
Xsize = atoi(argv[i]);
else
{
fprintf(stderr, MSG_NUMBER_EXPECTED, argv[i-1]);
exit(EXIT_FAILURE);
}
if (strcmp(argv[i], "--loop-size")==0)
if (str_isnumber(argv[++i]))
LOOP_SIZE_SEC = atoi(argv[i]);
else
{
fprintf(stderr, MSG_NUMBER_EXPECTED, argv[i-1]);
exit(EXIT_FAILURE);
}
if (strcmp(argv[i], "--log-base")==0)
if (str_isnumber(argv[++i]))
LOGBASE = atof(argv[i]);
else
{
fprintf(stderr, MSG_NUMBER_EXPECTED, argv[i-1]);
exit(EXIT_FAILURE);
}
if (strcmp(argv[i], "--bmsq-lut-size")==0)
if (str_isnumber(argv[++i]))
BMSQ_LUT_SIZE = atoi(argv[i]);
else
{
fprintf(stderr, MSG_NUMBER_EXPECTED, argv[i-1]);
exit(EXIT_FAILURE);
}
if (strcmp(argv[i], "--pi")==0) // lol
if (str_isnumber(argv[++i]))
pi = atof(argv[i]);
else
{
fprintf(stderr, MSG_NUMBER_EXPECTED, argv[i-1]);
exit(EXIT_FAILURE);
}
if (strcmp(argv[i], "--format-param")==0 || strcmp(argv[i], "-f")==0)
if (str_isnumber(argv[++i]))
format_param = atoi(argv[i]);
else
{
fprintf(stderr, MSG_NUMBER_EXPECTED, argv[i-1]);
exit(EXIT_FAILURE);
}
if (strcmp(argv[i], "--brightness")==0 || strcmp(argv[i], "-g")==0)
if (str_isnumber(argv[++i]))
brightness = atof(argv[i]);
else
{
fprintf(stderr, MSG_NUMBER_EXPECTED, argv[i-1]);
exit(EXIT_FAILURE);
}
// TODO implement --duration, -d
if (strcmp(argv[i], "--version")==0 || strcmp(argv[i], "-v")==0)
{
printf("Copyright (C) 2005-2008 Michel Rouzic\nProgram last modified by its author on %s\n", date);
exit(EXIT_SUCCESS);
}
if (strcmp(argv[i], "--help")==0 || strcmp(argv[i], "-h")==0)
{
print_help();
exit(EXIT_SUCCESS);
}
if (strcmp(argv[i], "--adv-help")==0)
{
print_adv_help();
exit(EXIT_SUCCESS);
}
}
}
if (in_name!=NULL) // if in_name has already been filled in
{
fin=fopen(in_name, "rb"); // try to open it
if (fin==NULL)
{
fprintf(stderr, "The input file %s could not be found\nExiting with error.\n", in_name);
exit(EXIT_FAILURE);
}
printf("Input file : %s\n", in_name);
}
else
{
if (quiet==1)
{
fprintf(stderr, "Please specify an input file.\nExiting with error.\n");
exit(EXIT_FAILURE);
}
printf("Type 'help' to read the manual page\n");
do
{
printf("Input file : ");
in_name=getstring();
if (strcmp(in_name, "help")==0) // if 'help' has been typed
{
fin=NULL;
print_help(); // print the help
}
else
fin=fopen(in_name, "rb");
}
while (fin==NULL);
}
if (out_name!=NULL) // if out_name has already been filled in
{
fout=fopen(out_name, "wb");
if (fout==NULL)
{
fprintf(stderr, "The output file %s could not be opened.\nPlease make sure it isn't opened by any other program and press Return.\nExiting with error.\n", out_name);
exit(EXIT_FAILURE);
}
printf("Output file : %s\n", out_name);
}
else
{
if (quiet==1)
{
fprintf(stderr, "Please specify an output file.\nExiting with error.\n");
exit(EXIT_FAILURE);
}
printf("Output file : ");
out_name=getstring();
fout=fopen(out_name, "wb");
while (fout==NULL)
{
fprintf(stderr, "The output file %s could not be opened.\nPlease make sure it isn't opened by any other program and press Return.\n", out_name);
getchar();
fout=fopen(out_name, "wb");
}
}
for (i=0; i<strlen(in_name); i++) if (in_name[i]>='A' && in_name[i]<='Z') in_name[i] -= 'A' - 'a'; // make the string lowercase
if (mode==0 && strstr(in_name, ".wav")!=NULL && strstr(in_name, ".wav")[4]==0)
mode=1; // Automatic switch to the Analysis mode
if (mode==0)
do
{
if (quiet==1)
{
fprintf(stderr, "Please specify an operation mode.\nUse either --analysis (-a), --sine (-s) or --noise (-n).\nExiting with error.\n");
exit(EXIT_FAILURE);
}
printf("Choose the mode (Press 1, 2 or 3) :\n\t1. Analysis\n\t2. Sine synthesis\n\t3. Noise synthesis\n> ");
mode=getfloat();
}
while (mode!=1 && mode!=2 && mode!=3);
if (mode==1)
{
sound=wav_in(fin, &channels, &samplecount, &samplerate); // Sound input
log_file("Read sound...");
#ifdef DEBUG
printf("samplecount : %i\nchannels : %i\n", samplecount, channels);
#endif
settingsinput(&Ysize, samplecount, &samplerate, &basefreq, maxfreq, &pixpersec, &bpo, Xsize, 0); // User settings input
image = anal(sound[0], samplecount, samplerate, &Xsize, Ysize, bpo, pixpersec, basefreq); // Analysis
if (brightness!=1.0)
brightness_control(image, Ysize, Xsize, 1.0/brightness);
bmp_out(fout, image, Ysize, Xsize); // Image output
}
if (mode==2 || mode==3)
{
sound = calloc (1, sizeof(double *));
image = bmp_in(fin, &Ysize, &Xsize); // Image input
if (format_param==0) // if the output format parameter is undefined
if (quiet==0) // if prompting is allowed
format_param = wav_out_param();
else
format_param = 32; // default is 32
settingsinput(&Ysize, samplecount, &samplerate, &basefreq, maxfreq, &pixpersec, &bpo, Xsize, 1); // User settings input
if (brightness!=1.0) brightness_control(image, Ysize, Xsize, brightness);
if (mode==2)
sound[0] = synt_sine(image, Xsize, Ysize, &samplecount, samplerate, basefreq, pixpersec, bpo); // Sine synthesis
else
sound[0] = synt_noise(image, Xsize, Ysize, &samplecount, samplerate, basefreq, pixpersec, bpo); // Noise synthesis
wav_out(fout, sound, 1, samplecount, samplerate, format_param);
}
clockb=gettime();
printf("Processing time : %.3f s\n", (double) (clockb-clocka)/1000.0);
win_return();
return 0;
}
