-/*
- Copyright (C) 2003 Paul Brossier
-
- This program is free software; you can redistribute it and/or modify
- it under the terms of the GNU General Public License as published by
- the Free Software Foundation; either version 2 of the License, or
- (at your option) any later version.
-
- This program is distributed in the hope that it will be useful,
- but WITHOUT ANY WARRANTY; without even the implied warranty of
- MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
- GNU General Public License for more details.
-
- You should have received a copy of the GNU General Public License
- along with this program; if not, write to the Free Software
- Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
-
-*/
-
-#include <stdio.h>
-#include <stdlib.h>
-#include <stdarg.h>
-#include <getopt.h>
-#include <unistd.h>
-
-#include "aubio.h"
-
-#include "utils.h"
-
-
-const char * output_filename = NULL;
-const char * input_filename = NULL;
-
-const char * onset_filename = "/usr/share/sounds/aubio/woodblock.aiff";
-
-int verbose = 0;
-int usejack = 0;
-
-int aubio_process(float **input, float **output, int nframes);
-
-/* settings */
-
-/* energy,specdiff,hfc,complexdomain,phase */
-aubio_onsetdetection_type type_onset;
-smpl_t threshold = 0.1;
-
-smpl_t threshold2 = -80.;
-
-int buffer_size = 1024;
-int overlap_size = 512;
-int channels = 1;
-int samplerate = 44100;
-
-/* global objects */
-aubio_pvoc_t * pv;
-
-fvec_t * ibuf;
-fvec_t * obuf;
-cvec_t * fftgrain;
-
-fvec_t * woodblock;
-
-aubio_onsetdetection_t *o;
-fvec_t *onset;
-//fvec_t *onset_keep;
-//fvec_t *onset_proc;
-int post = 5;
-int pre = 1;
-
-int isonset = 0;
-
-aubio_pickpeak_t * parms;
-
-int aubio_process(float **input, float **output, int nframes) {
- unsigned int i; /*channels*/
- unsigned int j; /*frames*/
- unsigned int pos = 0; /*frames%dspblocksize*/
- //unsigned int ons = 0; /*onset counter*/
- for (j=0;j<nframes;j++) {
-
- for (i=0;i<channels;i++) {
- /* write input to datanew */
- fvec_write_sample(ibuf, input[i][j], i, pos);
- /* put synthnew in output */
- output[i][j] = fvec_read_sample(obuf, i, pos);
- }
- /*time for fft*/
- if (pos == overlap_size-1) {
- aubio_pvoc_do (pv,ibuf, fftgrain);
- aubio_onsetdetection(o,fftgrain, onset);
- fprintf(stderr,"onsetfunc: %f\n",onset->data[0][0]);
- isonset = aubio_peakpick_pimrt(onset,parms);
- if (isonset) {
- /* test for silence */
- if (aubio_silence_detection(ibuf, threshold2)==1)
- isonset=0;
- else
- for (pos = 0; pos < overlap_size; pos++)
- obuf->data[0][pos] =
- woodblock->data[0][pos];
- } else {
- for (pos = 0; pos < overlap_size; pos++)
- obuf->data[0][pos] = 0.;
- }
- //aubio_pvoc_rdo(pv,fftgrain, obuf);
- pos = -1; /* so it will be zero next j loop */
- }
- pos++;
- }
- return 1;
-}
-
-int main (int argc, char **argv) {
- int frames;
- //int pos;
-
- aubio_file_t * file = NULL;
- aubio_file_t * fileout = NULL;
-
- aubio_file_t * onsetfile = new_file_ro(onset_filename);
- parse_args(argc, argv);
-
- if(!usejack)
- {
- debug("Opening files ...\n");
- file = new_file_ro (input_filename);
- file_info(file);
- channels = aubio_file_channels(file);
- if (output_filename != NULL)
- fileout = new_file_wo(file, output_filename);
- }
-
- ibuf = new_fvec(overlap_size, channels);
- obuf = new_fvec(overlap_size, channels);
- woodblock = new_fvec(buffer_size,1);
- fftgrain = new_cvec(buffer_size, channels);
-
- file_read(onsetfile, overlap_size, woodblock);
-
- /* phase vocoder */
- debug("Phase voc init ... \n");
- pv = new_aubio_pvoc(buffer_size, overlap_size, channels);
-
- /* onset init */
- debug("type_onset: %d\n", (int)type_onset);
- debug("threshold : %f\n", threshold);
- parms = new_aubio_peakpicker(threshold);
-
- o = new_aubio_onsetdetection(type_onset,buffer_size,channels);
- /*onset = new_fvec(overlap_size, channels);
- for (frames=0;frames<overlap_size;frames++){
- onset->data[0][frames] = 0.;
- }*/
- onset = new_fvec(0, channels);
-
- //onset_keep = new_fvec(post+pre+1, channels);
- //onset_proc = new_fvec(post+pre+1, channels);
-
- // command line argument parsing
- if(usejack) {
-#ifdef JACK_SUPPORT
- aubio_jack_t * jack_setup;
- debug("Jack init ...\n");
- jack_setup = new_aubio_jack(channels, channels,
- (aubio_process_func_t)aubio_process);
-
- /* initialise fftgrain */
- aubio_pvoc_do (pv,ibuf, fftgrain);
-
- debug("Jack activation ...\n");
- aubio_jack_activate(jack_setup);
-
- debug("Processing (Ctrl+C to quit) ...\n");
- pause();
-
- aubio_jack_close(jack_setup);
-#endif
-
- } else {
- fvec_t * inputbuf ;
- fvec_t * outputbuf;
- /* file objects */
-
- inputbuf = new_fvec(overlap_size, channels);
- outputbuf = new_fvec(overlap_size, channels);
- /* phasevoc init */
- debug("Processing ...\n");
-
- /** \bug phasevoc not initialized yet */
- /* initialise fftgrain and input bug */
- for (frames=0;frames<2;frames++) {
- aubio_pvoc_do (pv,inputbuf, fftgrain);
- }
-
- /** start filling the phasevoc */
- //for (frames=0;frames<1;frames++) {
- //file_read(file, overlap_size, inputbuf);
- //aubio_pvoc_do (pv,inputbuf, fftgrain);
- /** \bug onset is still badly initialized */
- //aubio_onsetdetection(o,fftgrain, onset);
- //}
-
- //debug("coucou\n");
- frames = 0;
- while ((overlap_size == file_read(file, overlap_size, inputbuf)) )
- {
- isonset=0;
- aubio_process(inputbuf->data, outputbuf->data, overlap_size);
- if (isonset) {
- /*output times in seconds*/
- //errmsg("%f\n",(frames-1)*overlap_size/(float)samplerate);
- errmsg("%f\n",(frames-3)*overlap_size/(float)samplerate);
- }
- if (output_filename != NULL)
- file_write(fileout,overlap_size,outputbuf);
- frames++;
- }
- debug("Processed %d frames of %d samples.\n", frames, buffer_size);
- del_file(file);
- if (output_filename != NULL)
- del_file(fileout);
- }
-
- del_aubio_pvoc(pv);
-
- del_fvec(obuf);
- del_fvec(ibuf);
- del_cvec(fftgrain);
- //del_fvec(onset_keep);
- //del_fvec(onset_proc);
- del_fvec(onset);
-
- debug("End of program.\n");
- fflush(stderr);
- return 0;
-}
-
-