2 Copyright (C) 2003-2014 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/>.
21 #include "aubio_priv.h"
24 #include "mathutils.h"
25 #include "spectral/fft.h"
26 #include "spectral/phasevoc.h"
28 /** phasevocoder internal object */
29 struct _aubio_pvoc_t {
30 uint_t win_s; /** grain length */
31 uint_t hop_s; /** overlap step */
32 aubio_fft_t * fft; /** fft object */
33 fvec_t * data; /** current input grain, [win_s] frames */
34 fvec_t * dataold; /** memory of past grain, [win_s-hop_s] frames */
35 fvec_t * synth; /** current output grain, [win_s] frames */
36 fvec_t * synthold; /** memory of past grain, [win_s-hop_s] frames */
37 fvec_t * w; /** grain window [win_s] */
38 uint_t start; /** where to start additive synthesis */
39 uint_t end; /** where to end it */
40 smpl_t scale; /** scaling factor for synthesis */
41 uint_t end_datasize; /** size of memory to end */
42 uint_t hop_datasize; /** size of memory to hop_s */
46 /** returns data and dataold slided by hop_s */
47 static void aubio_pvoc_swapbuffers(aubio_pvoc_t *pv, const fvec_t *new);
49 /** do additive synthesis from 'old' and 'cur' */
50 static void aubio_pvoc_addsynth(aubio_pvoc_t *pv, fvec_t * synthnew);
52 void aubio_pvoc_do(aubio_pvoc_t *pv, const fvec_t * datanew, cvec_t *fftgrain) {
54 aubio_pvoc_swapbuffers(pv, datanew);
56 fvec_weight(pv->data, pv->w);
60 aubio_fft_do (pv->fft,pv->data,fftgrain);
63 void aubio_pvoc_rdo(aubio_pvoc_t *pv,cvec_t * fftgrain, fvec_t * synthnew) {
65 aubio_fft_rdo(pv->fft,fftgrain,pv->synth);
67 fvec_ishift(pv->synth);
69 // if overlap = 50%, do not apply window (identity)
70 if (pv->hop_s * 2 < pv->win_s) {
71 fvec_weight(pv->synth, pv->w);
73 /* additive synthesis */
74 aubio_pvoc_addsynth(pv, synthnew);
77 aubio_pvoc_t * new_aubio_pvoc (uint_t win_s, uint_t hop_s) {
78 aubio_pvoc_t * pv = AUBIO_NEW(aubio_pvoc_t);
80 /* if (win_s < 2*hop_s) {
81 AUBIO_WRN("Hop size bigger than half the window size!\n");
84 if ((sint_t)hop_s < 1) {
85 AUBIO_ERR("pvoc: got hop_size %d, but can not be < 1\n", hop_s);
87 } else if ((sint_t)win_s < 2) {
88 AUBIO_ERR("pvoc: got buffer_size %d, but can not be < 2\n", win_s);
90 } else if (win_s < hop_s) {
91 AUBIO_ERR("pvoc: hop size (%d) is larger than win size (%d)\n", hop_s, win_s);
95 pv->fft = new_aubio_fft (win_s);
96 if (pv->fft == NULL) {
101 pv->data = new_fvec (win_s);
102 pv->synth = new_fvec (win_s);
104 /* new input output */
106 pv->dataold = new_fvec (win_s-hop_s);
107 pv->synthold = new_fvec (win_s-hop_s);
109 pv->dataold = new_fvec (1);
110 pv->synthold = new_fvec (1);
112 pv->w = new_aubio_window ("hanningz", win_s);
117 /* more than 50% overlap, overlap anyway */
118 if (win_s < 2 * hop_s) pv->start = 0;
119 /* less than 50% overlap, reset latest grain trail */
120 else pv->start = win_s - hop_s - hop_s;
122 if (win_s > hop_s) pv->end = win_s - hop_s;
125 pv->end_datasize = pv->end * sizeof(smpl_t);
126 pv->hop_datasize = pv->hop_s * sizeof(smpl_t);
128 // for reconstruction with 75% overlap
129 if (win_s == hop_s * 4) {
131 } else if (win_s == hop_s * 8) {
133 } else if (win_s == hop_s * 2) {
146 uint_t aubio_pvoc_set_window(aubio_pvoc_t *pv, const char_t *window) {
147 return fvec_set_window(pv->w, (char_t*)window);
150 void del_aubio_pvoc(aubio_pvoc_t *pv) {
153 del_fvec(pv->dataold);
154 del_fvec(pv->synthold);
156 del_aubio_fft(pv->fft);
160 static void aubio_pvoc_swapbuffers(aubio_pvoc_t *pv, const fvec_t *new)
162 /* some convenience pointers */
163 smpl_t * data = pv->data->data;
164 smpl_t * dataold = pv->dataold->data;
165 smpl_t * datanew = new->data;
166 #ifndef HAVE_MEMCPY_HACKS
168 for (i = 0; i < pv->end; i++)
169 data[i] = dataold[i];
170 for (i = 0; i < pv->hop_s; i++)
171 data[pv->end + i] = datanew[i];
172 for (i = 0; i < pv->end; i++)
173 dataold[i] = data[i + pv->hop_s];
175 memcpy(data, dataold, pv->end_datasize);
177 memcpy(data, datanew, pv->hop_datasize);
180 memcpy(dataold, data, pv->end_datasize);
184 static void aubio_pvoc_addsynth(aubio_pvoc_t *pv, fvec_t *synth_new)
187 /* some convenience pointers */
188 smpl_t * synth = pv->synth->data;
189 smpl_t * synthold = pv->synthold->data;
190 smpl_t * synthnew = synth_new->data;
192 /* put new result in synthnew */
193 for (i = 0; i < pv->hop_s; i++)
194 synthnew[i] = synth[i] * pv->scale;
196 /* no overlap, nothing else to do */
197 if (pv->end == 0) return;
199 /* add new synth to old one */
200 for (i = 0; i < pv->hop_s; i++)
201 synthnew[i] += synthold[i];
204 for (i = 0; i < pv->start; i++)
205 synthold[i] = synthold[i + pv->hop_s];
207 /* erase last frame in synthold */
208 for (i = pv->start; i < pv->end; i++)
212 for (i = 0; i < pv->end; i++)
213 synthold[i] += synth[i + pv->hop_s] * pv->scale;