2 Copyright (C) 2007 Amaury Hazan
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.
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.
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.
21 unsigned int pos = 0; /*frames%dspblocksize*/
24 int aubio_process(float **input, float **output, int nframes);
25 int aubio_process(float **input, float **output, int nframes) {
26 unsigned int i; /*channels*/
27 unsigned int j; /*frames*/
28 for (j=0;j<(unsigned)nframes;j++) {
30 for (i=0;i<channels;i++) {
31 /* write input to datanew */
32 fvec_write_sample(ibuf, input[i][j], i, pos);
33 /* put synthnew in output */
34 output[i][j] = fvec_read_sample(obuf, i, pos);
38 if (pos == overlap_size-1) {
41 //compute mag spectrum
42 aubio_pvoc_do (pv,ibuf, fftgrain);
46 aubio_mffc_do(fftgrain->norm, nframes, filterbank, outbuf);
50 /* end of block loop */
51 pos = -1; /* so it will be zero next j loop */
58 void process_print (void);
59 void process_print (void) {
60 /* output times in seconds
64 if (output_filename == NULL) {
66 outmsg("%f\n",(frames-4)*overlap_size/(float)samplerate);
67 } else if (frames < 4) {
73 int main(int argc, char **argv) {
74 examples_common_init(argc,argv);
76 //allocate and initialize mel filter bank
78 uint_t nyquist= samplerate / 2.;
80 uint_t banksize = (uint) ( sizeof(aubio_mel_filter));
81 aubio_mel_filter * mf = (aubio_mel_filter *)getbytes(banksize);
83 mfilterbank->n_filters = 20;
84 mfilterbank->filters = (smpl_t **)getbytes(mf->n_filters * sizeof(smpl_t *));
85 for(n = 0; n < mf->n_filters; n++)
86 mf->filters[n] = (smpl_t *)getbytes((buffer_size/2+1) * sizeof(smpl_t));
88 //populating the filter
89 new_aubio_mfcc(buffer_size, nyquist, XTRACT_EQUAL_GAIN, 80.0f, 18000.0f, mf->n_filters, mf->filters);
92 examples_common_process(aubio_process,process_print);
93 examples_common_del();
94 debug("End of program.\n");
97 //destroying filterbank