move libc aliases from mathutils.h to to aubio_priv.h
[aubio.git] / src / resample.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
20
21 #include <samplerate.h> /* from libsamplerate */
22
23 #include "aubio_priv.h"
24 #include "sample.h"
25 #include "resample.h"
26
27 struct _aubio_resampler_t {
28         SRC_DATA  *proc;
29         SRC_STATE *stat;
30         float ratio;
31         uint_t type;
32 };
33
34 aubio_resampler_t * new_aubio_resampler(float ratio, uint_t type) {
35         aubio_resampler_t * s  = AUBIO_NEW(aubio_resampler_t);
36         int error = 0;
37         s->stat = src_new (type, 1, &error) ; /* only one channel */
38         s->proc = AUBIO_NEW(SRC_DATA);
39         if (error) AUBIO_ERR("%s\n",src_strerror(error));
40         s->ratio = ratio;
41         return s;
42 }
43
44 void del_aubio_resampler(aubio_resampler_t *s) {
45         src_delete(s->stat);
46         AUBIO_FREE(s->proc);
47         AUBIO_FREE(s);
48 }
49
50 uint_t aubio_resampler_process(aubio_resampler_t *s, 
51     fvec_t * input,  fvec_t * output) {
52         uint_t i ;
53         s->proc->input_frames = input->length;
54         s->proc->output_frames = output->length;
55         s->proc->src_ratio = s->ratio;
56         for (i = 0 ; i< input->channels; i++) 
57         {
58                 /* make SRC_PROC data point to input outputs */
59                 s->proc->data_in = input->data[i];
60                 s->proc->data_out= output->data[i];
61                 /* do resampling */
62                 src_process (s->stat, s->proc) ;
63         }
64         return AUBIO_OK;
65 }