#include <stdarg.h>
#include <getopt.h>
#include <unistd.h>
-#include <math.h> // how do i do a floorf with a mask again ?
+#include <math.h>
#include "aubio.h"
#include "aubioext.h"
#include "utils.h"
int usejack = 0;
int usedoubled = 1;
int usemidi = 1;
+int median = 0;
/* energy,specdiff,hfc,complexdomain,phase */
aubio_onsetdetection_type type_onset = hfc;
-aubio_onsetdetection_type type_onset2 = hfc;
+aubio_onsetdetection_type type_onset2 = complexdomain;
smpl_t threshold = 0.3;
smpl_t threshold2 = -90.;
uint_t buffer_size = 1024;
aubio_pickpeak_t * parms;
/* pitch objects */
-int curnote = 0;
smpl_t pitch = 0.;
aubio_pitchdetection_t * pitchdet;
-int bufpos = 0;
-smpl_t curlevel = 0;
+
+fvec_t * note_buffer = NULL;
+fvec_t * note_buffer2 = NULL;
+uint_t medianfiltlen = 6;
+smpl_t curlevel = 0.;
smpl_t maxonset = 0.;
/* midi objects */
}
}
+
+/** append new note candidate to the note_buffer and return filtered value. we
+ * need to copy the input array as vec_median destroy its input data.*/
+
+void note_append(fvec_t * note_buffer, smpl_t curnote);
+void note_append(fvec_t * note_buffer, smpl_t curnote) {
+ uint_t i = 0;
+ for (i = 0; i < note_buffer->length - 1; i++) {
+ note_buffer->data[0][i] = note_buffer->data[0][i+1];
+ }
+ note_buffer->data[0][note_buffer->length - 1] = curnote;
+ return;
+}
+
+uint_t get_note(fvec_t *note_buffer, fvec_t *note_buffer2);
+uint_t get_note(fvec_t *note_buffer, fvec_t *note_buffer2){
+ uint_t i = 0;
+ for (i = 0; i < note_buffer->length; i++) {
+ note_buffer2->data[0][i] = note_buffer->data[0][i];
+ }
+ return vec_median(note_buffer2);
+}
+
+
+smpl_t curnote = 0.;
+smpl_t newnote = 0.;
+uint_t isready = 0;
+
int aubio_process(float **input, float **output, int nframes);
int aubio_process(float **input, float **output, int nframes) {
unsigned int i; /*channels*/
isonset = aubio_peakpick_pimrt(onset,parms);
pitch = aubio_pitchdetection(pitchdet,ibuf);
+ if(median){
+ newnote = (FLOOR)(bintomidi(pitch,samplerate,buffer_size*4)+.5);
+ note_append(note_buffer, newnote);
+ }
/* curlevel is negatif or 1 if silence */
curlevel = aubio_level_detection(ibuf, threshold2);
/* test for silence */
if (curlevel == 1.) {
isonset=0;
+ if (median) isready = 0;
/* send note off */
send_noteon(mdriver,curnote,0);
} else {
+ if (median) {
+ isready = 1;
+ } else {
/* kill old note */
send_noteon(mdriver,curnote,0);
//curnote = (int)FLOOR(bintomidi(pitch,samplerate,buffer_size*4) + .5);
if (curnote>45){
send_noteon(mdriver,curnote,127+(int)FLOOR(curlevel));
}
+ }
+
for (pos = 0; pos < overlap_size; pos++){
obuf->data[0][pos] = woodblock->data[0][pos];
}
}
} else {
+ if (median) {
+// if (curlevel != 1) {
+// if (bufpos == note_buffer->length)
+// curnote = aubio_getnote(note_buffer);
+// }
+//
+ if (isready > 0)
+ isready++;
+ if (isready == note_buffer->length)
+ {
+ //outmsg("%d, %f\n", isready, curnote);
+ send_noteon(mdriver,curnote,0);
+ /* kill old note */
+ newnote = get_note(note_buffer, note_buffer2);
+ curnote = newnote;
+ /* get and send new one */
+ /*if (curnote<45){
+ send_noteon(mdriver,curnote,0);
+ } else {*/
+ if (curnote>45){
+ send_noteon(mdriver,curnote,127+(int)FLOOR(curlevel));
+ }
+ }
+ } // if median
for (pos = 0; pos < overlap_size; pos++)
obuf->data[0][pos] = 0.;
}
woodblock = new_fvec(buffer_size,1);
fftgrain = new_cvec(buffer_size, channels);
- pitchdet = new_aubio_pitchdetection(buffer_size*4, overlap_size, channels, samplerate, mcomb, freq);
- //pitchdet = new_aubio_pitchdetection(buffer_size*2, overlap_size, channels, samplerate, yin, freq);
+ aubio_pitchdetection_type mode = yin; // mcomb
+ pitchdet = new_aubio_pitchdetection(buffer_size*4, overlap_size, channels, samplerate, mode, freq);
+if (median) {
+ note_buffer = new_fvec(medianfiltlen, 1);
+ note_buffer2= new_fvec(medianfiltlen, 1);
+}
/* read the output sound once */
file_read(onsetfile, overlap_size, woodblock);
/* phase vocoder */
del_cvec(fftgrain);
del_aubio_pitchdetection(pitchdet);
del_fvec(onset);
+ if (median) {
+ del_fvec(note_buffer);
+ del_fvec(note_buffer2);
+ }
debug("End of program.\n");
fflush(stderr);
+++ /dev/null
-/*
- 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 "aubioext.h"
-#include "utils.h"
-
-#include <math.h> //required for FLOORF :(
-
-/* settings */
-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 usedoubled = 1;
-int usemidi = 1;
-
-/* energy,specdiff,hfc,complexdomain,phase */
-aubio_onsetdetection_type type_onset = hfc;
-aubio_onsetdetection_type type_onset2 = complexdomain;
-smpl_t threshold = 0.3;
-smpl_t threshold2 = -90.;
-uint_t buffer_size = 1024;
-uint_t overlap_size = 512;
-uint_t channels = 1;
-uint_t samplerate = 44100;
-
-/* global objects */
-int frames;
-aubio_pvoc_t * pv;
-fvec_t * ibuf;
-fvec_t * obuf;
-cvec_t * fftgrain;
-fvec_t * woodblock;
-aubio_onsetdetection_t *o;
-aubio_onsetdetection_t *o2;
-fvec_t *onset;
-fvec_t *onset2;
-int isonset = 0;
-aubio_pickpeak_t * parms;
-
-/* pitch objects */
-smpl_t pitch = 0.;
-cvec_t * fftpitch = NULL;
-aubio_filter_t * filter = NULL;
-aubio_pvoc_t * pvi = NULL;
-aubio_pitchmcomb_t * p = NULL;
-
-fvec_t * note_buffer = NULL;
-fvec_t * note_buffer2 = NULL;
-uint_t medianfiltlen = 6;
-smpl_t curlevel = 0;
-smpl_t maxonset = 0.;
-
-/* midi objects */
-aubio_midi_player_t * mplay;
-aubio_midi_driver_t * mdriver;
-aubio_midi_event_t * event;
-
-void send_noteon(aubio_midi_driver_t * d, int pitch, int velo);
-void send_noteon(aubio_midi_driver_t * d, int pitch, int velo)
-{
- /* we should check if we use midi here, not jack */
-#if ALSA_SUPPORT
- if (usejack) {
- if (velo==0) {
- aubio_midi_event_set_type(event,NOTE_OFF);
- } else {
- aubio_midi_event_set_type(event,NOTE_ON);
- }
- aubio_midi_event_set_channel(event,0);
- aubio_midi_event_set_pitch(event,pitch);
- aubio_midi_event_set_velocity(event,velo);
- aubio_midi_direct_output(mdriver,event);
- } else
-#endif
- {
- if (velo==0) {
- outmsg("%f\n",frames*overlap_size/(float)samplerate);
- } else {
- outmsg("%d\t%f\t",pitch,frames*overlap_size/(float)samplerate);
- }
- }
-}
-
-
-/** append new note candidate to the note_buffer and return filtered value. we
- * need to copy the input array as vec_median destroy its input data.*/
-
-void note_append(fvec_t * note_buffer, smpl_t curnote);
-void note_append(fvec_t * note_buffer, smpl_t curnote) {
- uint_t i = 0;
- for (i = 0; i < note_buffer->length - 1; i++) {
- note_buffer->data[0][i] = note_buffer->data[0][i+1];
- }
- note_buffer->data[0][note_buffer->length - 1] = curnote;
- return;
-}
-
-uint_t get_note(fvec_t *note_buffer, fvec_t *note_buffer2);
-uint_t get_note(fvec_t *note_buffer, fvec_t *note_buffer2){
- uint_t i = 0;
- for (i = 0; i < note_buffer->length; i++) {
- note_buffer2->data[0][i] = note_buffer->data[0][i];
- }
- return vec_median(note_buffer2);
-}
-
-
-smpl_t curnote = 0.;
-smpl_t newnote = 0.;
-uint_t isready = 0;
-
-int aubio_process(float **input, float **output, int nframes);
-int aubio_process(float **input, float **output, int nframes) {
- unsigned int i; /*channels*/
- unsigned int j; /*frames*/
- unsigned int pos = 0; /*frames%dspblocksize*/
- for (j=0;j<nframes;j++) {
- if(usejack) {
- 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) {
- /* block loop */
- aubio_pvoc_do (pv,ibuf, fftgrain);
- aubio_onsetdetection(o,fftgrain, onset);
- if (usedoubled) {
- aubio_onsetdetection(o2,fftgrain, onset2);
- onset->data[0][0] *= onset2->data[0][0];
- }
- isonset = aubio_peakpick_pimrt(onset,parms);
-
- aubio_filter_do(filter,ibuf);
- aubio_filter_do(filter,ibuf);
- aubio_pvoc_do(pvi,ibuf, fftpitch);
- pitch = aubio_pitchmcomb_detect(p,fftpitch);
- newnote = (FLOOR)(bintomidi(pitch,samplerate,buffer_size*4)+.5);
- note_append(note_buffer, newnote);
-
- //outmsg("%f %f\n", pitch, newnote);
-
- /* curlevel is negatif or 1 if silence */
- curlevel = aubio_level_detection(ibuf, threshold2);
- if (isonset) {
- /* test for silence */
- if (curlevel == 1.) {
- isonset=0;
- isready=0;
- /* send note off */
- send_noteon(mdriver,curnote,0);
- } else {
- isready = 1;
- for (pos = 0; pos < overlap_size; pos++){
- obuf->data[0][pos] = woodblock->data[0][pos];
- }
- }
- } else {
-// if (curlevel != 1) {
-// if (bufpos == note_buffer->length)
-// curnote = aubio_getnote(note_buffer);
-// }
-//
- if (isready > 0)
- isready++;
- if (isready == note_buffer->length)
- {
- //outmsg("%d, %f\n", isready, curnote);
- send_noteon(mdriver,curnote,0);
- /* kill old note */
- newnote = get_note(note_buffer, note_buffer2);
- curnote = newnote;
- /* get and send new one */
- /*if (curnote<45){
- send_noteon(mdriver,curnote,0);
- } else {*/
- if (curnote>45){
- send_noteon(mdriver,curnote,127+(int)FLOOR(curlevel));
- }
- }
- for (pos = 0; pos < overlap_size; pos++)
- obuf->data[0][pos] = 0.;
- }
- /* end of block loop */
- //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) {
-
- aubio_file_t * file = NULL;
- aubio_file_t * fileout = NULL;
-
- aubio_file_t * onsetfile = new_file_ro(onset_filename);
- /* parse command line arguments */
- 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);
-
- pvi = new_aubio_pvoc(buffer_size*4, overlap_size, channels);
- fftpitch = new_cvec(buffer_size*4, channels);
- filter = new_aubio_adsgn_filter((smpl_t)samplerate);
- p = new_aubio_pitchmcomb(buffer_size*4,channels);
-
- note_buffer = new_fvec(medianfiltlen, 1);
- note_buffer2= new_fvec(medianfiltlen, 1);
- /* read the output sound once */
- file_read(onsetfile, overlap_size, woodblock);
-
- /* phase vocoder */
- debug("Phase voc init ... \n");
- pv = new_aubio_pvoc(buffer_size, overlap_size, channels);
-
- /* onsets */
- parms = new_aubio_peakpicker(threshold);
- o = new_aubio_onsetdetection(type_onset,buffer_size,channels);
- onset = new_fvec(1, channels);
- if (usedoubled) {
- o2 = new_aubio_onsetdetection(type_onset2,buffer_size,channels);
- onset2 = new_fvec(1 , channels);
- }
-
- if(usejack) {
-#ifdef JACK_SUPPORT
- aubio_jack_t * jack_setup;
- debug("Midi init ...\n");
- debug("Jack init ...\n");
- jack_setup = new_aubio_jack(channels, channels,
- (aubio_process_func_t)aubio_process);
-
- mplay = new_aubio_midi_player();
-
- mdriver = new_aubio_midi_driver("alsa_seq",
- (handle_midi_event_func_t)aubio_midi_send_event, mplay);
-
- event = new_aubio_midi_event();
-
- debug("Jack activation ...\n");
- aubio_jack_activate(jack_setup);
- debug("Processing (Ctrl+C to quit) ...\n");
- pause();
- send_noteon(mdriver,curnote,0);
- aubio_jack_close(jack_setup);
- del_aubio_midi_driver(mdriver);
-#else
- outmsg("Compiled without jack output, exiting.");
-#endif
-
- } else {
- /* phasevoc */
- debug("Processing ...\n");
-
- frames = 0;
-
- while (overlap_size == file_read(file, overlap_size, ibuf))
- {
- isonset=0;
- aubio_process(ibuf->data, obuf->data, overlap_size);
- if (output_filename != NULL) {
- file_write(fileout,overlap_size,obuf);
- }
- frames++;
- }
- send_noteon(mdriver,curnote,0);
-
- 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_cvec(fftpitch);
- del_aubio_pvoc(pvi);
- del_fvec(onset);
- del_fvec(note_buffer);
- del_fvec(note_buffer2);
-
- debug("End of program.\n");
- fflush(stderr);
- return 0;
-}
-
projects / aubio.git / commitdiff