src/temporal/filter.c: reset filter memory in filtfilt
[aubio.git] / examples / aubiotrack.c
1 /*
2    Copyright (C) 2003 Paul Brossier
3
4    This program is free software; you can redistribute it and/or modify
5    it under the terms of the GNU General Public License as published by
6    the Free Software Foundation; either version 2 of the License, or
7    (at your option) any later version.
8
9    This program is distributed in the hope that it will be useful,
10    but WITHOUT ANY WARRANTY; without even the implied warranty of
11    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
12    GNU General Public License for more details.
13
14    You should have received a copy of the GNU General Public License
15    along with this program; if not, write to the Free Software
16    Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
17 */
18
19 #include <aubio.h>
20 #include "utils.h"
21
22 uint_t pos = 0;    /* frames%dspblocksize */
23 fvec_t * tempo_out = NULL;
24 aubio_tempo_t * bt = NULL;
25 smpl_t istactus = 0;
26 smpl_t isonset = 0;
27
28 static int aubio_process(smpl_t **input, smpl_t **output, int nframes) {
29   unsigned int i;       /*channels*/
30   unsigned int j;       /*frames*/
31   for (j=0;j<(unsigned)nframes;j++) {
32     if(usejack) {
33       for (i=0;i<channels;i++) {
34         /* write input to datanew */
35         fvec_write_sample(ibuf, input[i][j], i, pos);
36         /* put synthnew in output */
37         output[i][j] = fvec_read_sample(obuf, i, pos);
38       }
39     }
40     /*time for fft*/
41     if (pos == overlap_size-1) {         
42       /* block loop */
43       aubio_tempo_do (bt,ibuf,tempo_out);
44       if (tempo_out->data[0][0]>0) 
45         istactus = tempo_out->data[0][0];
46       else 
47         istactus = 0;
48       if (tempo_out->data[0][1]>0) 
49         isonset = tempo_out->data[0][0];
50       else 
51         isonset = 0;
52       if (istactus) {
53               for (pos = 0; pos < overlap_size; pos++)
54                       obuf->data[0][pos] = woodblock->data[0][pos];
55       } else {
56               for (pos = 0; pos < overlap_size; pos++)
57                       obuf->data[0][pos] = 0.;
58       }
59       /* end of block loop */
60       pos = -1; /* so it will be zero next j loop */
61     }
62     pos++;
63   }
64   return 1;
65 }
66
67 static void process_print (void) {
68         if (output_filename == NULL) {
69                 if (istactus) {
70                         outmsg("%f\n",((smpl_t)(frames*overlap_size)+(istactus-1.)*overlap_size)/(smpl_t)samplerate); 
71                 }
72                 if (isonset && verbose)
73                         outmsg(" \t \t%f\n",(frames)*overlap_size/(float)samplerate);
74         }
75 }
76
77 int main(int argc, char **argv) {
78   
79   buffer_size = 1024;
80   overlap_size = 512;
81   /* override default settings */
82   examples_common_init(argc,argv);
83
84   tempo_out = new_fvec(2,channels);
85   bt = new_aubio_tempo(onset_mode,buffer_size,overlap_size,channels, samplerate);
86
87   examples_common_process(aubio_process,process_print);
88
89   del_aubio_tempo(bt);
90   del_fvec(tempo_out);
91
92   examples_common_del();
93
94   debug("End of program.\n");
95
96   fflush(stderr);
97
98   return 0;
99 }
100