2 Copyright (C) 2007-2009 Paul Brossier <piem@aubio.org>
4 This file is part of aubio.
6 aubio is free software: you can redistribute it and/or modify
7 it under the terms of the GNU General Public License as published by
8 the Free Software Foundation, either version 3 of the License, or
9 (at your option) any later version.
11 aubio is distributed in the hope that it will be useful,
12 but WITHOUT ANY WARRANTY; without even the implied warranty of
13 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14 GNU General Public License for more details.
16 You should have received a copy of the GNU General Public License
17 along with aubio. If not, see <http://www.gnu.org/licenses/>.
29 uint_t n_filters = 40;
32 unsigned int pos = 0; /*frames%dspblocksize*/
34 static int aubio_process(smpl_t **input, smpl_t **output, int nframes) {
35 unsigned int j; /*frames*/
37 for (j=0;j<(unsigned)nframes;j++) {
39 /* write input to datanew */
40 fvec_write_sample(ibuf, input[0][j], pos);
41 /* put synthnew in output */
42 output[0][j] = fvec_read_sample(obuf, pos);
45 if (pos == overlap_size-1) {
48 //compute mag spectrum
49 aubio_pvoc_do (pv, ibuf, fftgrain);
52 aubio_mfcc_do(mfcc, fftgrain, mfcc_out);
54 /* end of block loop */
55 pos = -1; /* so it will be zero next j loop */
62 static void process_print (void) {
63 /* output times in seconds
68 if (sink_uri == NULL) {
69 outmsg("%f\t",frames*overlap_size/(float)samplerate);
70 for (coef_cnt = 0; coef_cnt < n_coefs; coef_cnt++) {
71 outmsg("%f ", fvec_read_sample (mfcc_out, coef_cnt) );
77 int main(int argc, char **argv) {
82 examples_common_init(argc,argv);
85 pv = new_aubio_pvoc (buffer_size, overlap_size);
87 fftgrain = new_cvec (buffer_size);
89 //populating the filter
90 mfcc = new_aubio_mfcc(buffer_size, n_filters, n_coefs, samplerate);
92 mfcc_out = new_fvec(n_coefs);
95 examples_common_process(aubio_process,process_print);
100 del_aubio_mfcc(mfcc);
103 examples_common_del();
104 debug("End of program.\n");