12 #include <math.h> /* for isfinite */
15 /* not supported yet */
17 #include <ladcca/ladcca.h>
18 cca_client_t * aubio_cca_client;
19 #endif /* LADCCA_SUPPORT */
22 const char * output_filename = NULL;
23 const char * input_filename = NULL;
24 const char * onset_filename = "/usr/share/sounds/aubio/woodblock.aiff";
31 /* energy,specdiff,hfc,complexdomain,phase */
32 aubio_onsetdetection_type type_onset = hfc;
33 aubio_onsetdetection_type type_onset2 = complexdomain;
34 smpl_t threshold = 0.3;
35 smpl_t threshold2 = -90.;
36 uint_t buffer_size = 1024;
37 uint_t overlap_size = 512;
39 uint_t samplerate = 44100;
42 aubio_file_t * file = NULL;
43 aubio_file_t * fileout = NULL;
50 aubio_onsetdetection_t *o;
51 aubio_onsetdetection_t *o2;
55 aubio_pickpeak_t * parms;
60 aubio_pitchdetection_t * pitchdet;
61 aubio_pitchdetection_type mode = aubio_yin; // aubio_mcomb
64 fvec_t * note_buffer = NULL;
65 fvec_t * note_buffer2 = NULL;
70 aubio_midi_player_t * mplay;
71 aubio_midi_driver_t * mdriver;
72 aubio_midi_event_t * event;
80 /* badly redeclare some things */
81 aubio_onsetdetection_type type_onset;
84 const char * prog_name;
86 void usage (FILE * stream, int exit_code)
88 fprintf(stream, "usage: %s [ options ] \n", prog_name);
90 " -j --jack Use Jack.\n"
91 " -o --output Output type.\n"
92 " -i --input Input type.\n"
93 " -h --help Display this message.\n"
94 " -v --verbose Print verbose message.\n"
99 int parse_args (int argc, char **argv) {
100 const char *options = "hvjo:i:O:t:s:H:a";
102 struct option long_options[] =
104 {"help" , 0, NULL, 'h'},
105 {"verbose" , 0, NULL, 'v'},
106 {"jack" , 0, NULL, 'j'},
107 {"output" , 0, NULL, 'o'},
108 {"input" , 0, NULL, 'i'},
109 {"onset" , 0, NULL, 'O'},
110 {"threshold", 0, NULL, 't'},
111 {"silence" , 0, NULL, 's'},
112 {"averaging", 0, NULL, 'a'},
113 {"hopsize", 0, NULL, 'H'},
122 next_option = getopt_long (argc, argv, options,
124 switch (next_option) {
126 output_filename = optarg;
129 input_filename = optarg;
134 case 'v': /* verbose */
137 case 'j': /* verbose */
140 case 'O': /*onset type*/
141 if (strcmp(optarg,"energy") == 0)
143 else if (strcmp(optarg,"specdiff") == 0)
144 type_onset = specdiff;
145 else if (strcmp(optarg,"hfc") == 0)
147 else if (strcmp(optarg,"complexdomain") == 0)
148 type_onset = complexdomain;
149 else if (strcmp(optarg,"phase") == 0)
152 debug("could not get onset type.\n");
157 case 's': /* threshold value for onset */
158 threshold2 = (smpl_t)atof(optarg);
160 case 't': /* threshold value for onset */
161 threshold = (smpl_t)atof(optarg);
163 if (!isfinite(threshold)) {
164 debug("could not get threshold.\n");
173 overlap_size = atoi(optarg);
175 case '?': /* unknown options */
178 case -1: /* done with options */
180 default: /*something else unexpected */
184 while (next_option != -1);
186 if (input_filename != NULL) {
187 debug ("Input file : %s\n", input_filename );
188 } else if (input_filename != NULL && output_filename != NULL) {
189 debug ("Input file : %s\n", input_filename );
190 debug ("Output file : %s\n", output_filename );
194 debug ("Jack input output\n");
197 debug ("Error: Could not switch to jack mode\n aubio was compiled without jack support\n");
204 void examples_common_init(int argc,char ** argv) {
207 aubio_file_t * onsetfile;
208 /* parse command line arguments */
209 parse_args(argc, argv);
211 woodblock = new_fvec(buffer_size,1);
212 if (output_filename || usejack) {
213 (onsetfile = new_file_ro(onset_filename)) ||
214 (onsetfile = new_file_ro("sounds/woodblock.aiff")) ||
215 (onsetfile = new_file_ro("../sounds/woodblock.aiff"));
216 /* read the output sound once */
217 file_read(onsetfile, overlap_size, woodblock);
222 debug("Opening files ...\n");
223 file = new_file_ro (input_filename);
224 if (verbose) file_info(file);
225 channels = aubio_file_channels(file);
226 samplerate = aubio_file_samplerate(file);
227 if (output_filename != NULL)
228 fileout = new_file_wo(file, output_filename);
231 ibuf = new_fvec(overlap_size, channels);
232 obuf = new_fvec(overlap_size, channels);
233 fftgrain = new_cvec(buffer_size, channels);
236 pitchdet = new_aubio_pitchdetection(buffer_size*4,
237 overlap_size, channels, samplerate, mode, aubio_freq);
240 note_buffer = new_fvec(median, 1);
241 note_buffer2= new_fvec(median, 1);
245 pv = new_aubio_pvoc(buffer_size, overlap_size, channels);
247 parms = new_aubio_peakpicker(threshold);
248 o = new_aubio_onsetdetection(type_onset,buffer_size,channels);
249 onset = new_fvec(1, channels);
251 o2 = new_aubio_onsetdetection(type_onset2,buffer_size,channels);
252 onset2 = new_fvec(1 , channels);
258 void examples_common_del(void){
260 send_noteon(curnote,0);
261 del_aubio_pitchdetection(pitchdet);
263 del_fvec(note_buffer);
264 del_fvec(note_buffer2);
274 void examples_common_process(aubio_process_func_t process_func, aubio_print_func_t print ){
277 aubio_jack_t * jack_setup;
278 debug("Jack init ...\n");
279 jack_setup = new_aubio_jack(channels, channels,
280 (aubio_process_func_t)process_func);
282 debug("Midi init ...\n");
283 mplay = new_aubio_midi_player();
284 mdriver = new_aubio_midi_driver("alsa_seq",
285 (handle_midi_event_func_t)aubio_midi_send_event, mplay);
286 event = new_aubio_midi_event();
288 debug("Jack activation ...\n");
289 aubio_jack_activate(jack_setup);
290 debug("Processing (Ctrl+C to quit) ...\n");
292 aubio_jack_close(jack_setup);
294 send_noteon(curnote,0);
295 del_aubio_midi_driver(mdriver);
299 outmsg("Compiled without jack output, exiting.\n");
304 debug("Processing ...\n");
308 while (overlap_size == file_read(file, overlap_size, ibuf))
311 process_func(ibuf->data, obuf->data, overlap_size);
313 if (output_filename != NULL) {
314 file_write(fileout,overlap_size,obuf);
319 debug("Processed %d frames of %d samples.\n", frames, buffer_size);
322 if (output_filename != NULL)
330 void send_noteon(int pitch, int velo)
332 smpl_t mpitch = (FLOOR)(freqtomidi(pitch)+.5);
333 /* we should check if we use midi here, not jack */
337 aubio_midi_event_set_type(event,NOTE_OFF);
339 aubio_midi_event_set_type(event,NOTE_ON);
341 aubio_midi_event_set_channel(event,0);
342 aubio_midi_event_set_pitch(event,mpitch);
343 aubio_midi_event_set_velocity(event,velo);
344 aubio_midi_direct_output(mdriver,event);
350 outmsg("%f\n",frames*overlap_size/(float)samplerate);
352 outmsg("%f\t%f\t", mpitch,
353 frames*overlap_size/(float)samplerate);
359 void note_append(fvec_t * note_buffer, smpl_t curnote) {
361 for (i = 0; i < note_buffer->length - 1; i++) {
362 note_buffer->data[0][i] = note_buffer->data[0][i+1];
364 note_buffer->data[0][note_buffer->length - 1] = curnote;
368 uint_t get_note(fvec_t *note_buffer, fvec_t *note_buffer2){
370 for (i = 0; i < note_buffer->length; i++) {
371 note_buffer2->data[0][i] = note_buffer->data[0][i];
373 return vec_median(note_buffer2);