Bug Summary

File:modules/hts21_engine/vocoder.cc
Location:line 205, column 4
Description:Value stored to 'aa' is never read

Annotated Source Code

1/* --------------------------------------------------------------- */
2/* The HMM-Based Speech Synthesis System (HTS): version 1.1b */
3/* HTS Working Group */
4/* */
5/* Department of Computer Science */
6/* Nagoya Institute of Technology */
7/* and */
8/* Interdisciplinary Graduate School of Science and Engineering */
9/* Tokyo Institute of Technology */
10/* Copyright (c) 2001-2003 */
11/* All Rights Reserved. */
12/* */
13/* Permission is hereby granted, free of charge, to use and */
14/* distribute this software and its documentation without */
15/* restriction, including without limitation the rights to use, */
16/* copy, modify, merge, publish, distribute, sublicense, and/or */
17/* sell copies of this work, and to permit persons to whom this */
18/* work is furnished to do so, subject to the following conditions: */
19/* */
20/* 1. The code must retain the above copyright notice, this list */
21/* of conditions and the following disclaimer. */
22/* */
23/* 2. Any modifications must be clearly marked as such. */
24/* */
25/* NAGOYA INSTITUTE OF TECHNOLOGY, TOKYO INSITITUTE OF TECHNOLOGY, */
26/* HTS WORKING GROUP, AND THE CONTRIBUTORS TO THIS WORK DISCLAIM */
27/* ALL WARRANTIES WITH REGARD TO THIS SOFTWARE, INCLUDING ALL */
28/* IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS, IN NO EVENT */
29/* SHALL NAGOYA INSTITUTE OF TECHNOLOGY, TOKYO INSITITUTE OF */
30/* TECHNOLOGY, HTS WORKING GROUP, NOR THE CONTRIBUTORS BE LIABLE */
31/* FOR ANY SPECIAL, INDIRECT OR CONSEQUENTIAL DAMAGES OR ANY */
32/* DAMAGES WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, */
33/* WHETHER IN AN ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTUOUS */
34/* ACTION, ARISING OUT OF OR IN CONNECTION WITH THE USE OR */
35/* PERFORMANCE OF THIS SOFTWARE. */
36/* */
37/* --------------------------------------------------------------- */
38/* vocoder.c : mel-cepstral vocoder */
39/* (pulse/noise excitation & MLSA filter) */
40/* */
41/* 2003/06/11 by Heiga Zen */
42/* --------------------------------------------------------------- */
43
44#include <cstdio>
45#include <cstdlib>
46#include <cmath>
47#include "EST_walloc.h"
48
49#include "misc.h"
50#include "model.h"
51#include "defaults.h"
52#include "global.h"
53#include "vocoder.h"
54
55void init_vocoder(int m, VocoderSetup *vs)
56{
57 vs->fprd = FPERIOD80;
58 vs->iprd = IPERIOD1;
59 vs->seed = SEED1;
60 vs->pd = PADEORDER4;
61
62 vs->next = SEED1;
63 vs->gauss = GAUSS1;
64
65 vs->pade[ 0]=1.0;
66 vs->pade[ 1]=1.0; vs->pade[ 2]=0.0;
67 vs->pade[ 3]=1.0; vs->pade[ 4]=0.0; vs->pade[ 5]=0.0;
68 vs->pade[ 6]=1.0; vs->pade[ 7]=0.0; vs->pade[ 8]=0.0; vs->pade[ 9]=0.0;
69 vs->pade[10]=1.0; vs->pade[11]=0.4999273; vs->pade[12]=0.1067005; vs->pade[13]=0.01170221; vs->pade[14]=0.0005656279;
70 vs->pade[15]=1.0; vs->pade[16]=0.4999391; vs->pade[17]=0.1107098; vs->pade[18]=0.01369984; vs->pade[19]=0.0009564853;
71 vs->pade[20]=0.00003041721;
72
73 vs->rate=RATE16000;
74
75 vs->c = walloc(double,3*(m+1)+3*(vs->pd+1)+vs->pd*(m+2))((double *)safe_walloc(sizeof(double)*(3*(m+1)+3*(vs->pd+1
)+vs->pd*(m+2))));
76
77 vs->p1 = -1;
78 vs->sw = 0;
79 vs->x = 0x55555555;
80}
81
82void vocoder (double p, float *mc, int m, FILE *rawfp, globalP *gp, VocoderSetup *vs)
83{
84 double inc, x;
85 int i, j, k;
86 short xs;
87 double a = gp->ALPHA;
88
89 if (p!=0.0)
90 p = vs->rate / p; /* f0 -> pitch */
91
92 if (vs->p1 < 0) {
93 if (vs->gauss & (vs->seed != 1)) vs->next = srnd ((unsigned)vs->seed);
94
95 vs->p1 = p;
96 vs->pc = vs->p1;
97 vs->cc = vs->c + m + 1;
98 vs->cinc = vs->cc + m + 1;
99 vs->d1 = vs->cinc + m + 1;
100
101 mc2b(mc, vs->c, m, a);
102
103 return;
104 }
105
106 mc2b(mc, vs->cc, m, a);
107
108 for (k=0; k<=m; k++)
109 vs->cinc[k] = (vs->cc[k]-vs->c[k])*(double)vs->iprd/(double)vs->fprd;
110
111 if (vs->p1!=0.0 && p!=0.0) {
112 inc = (p-vs->p1)*(double)vs->iprd/(double)vs->fprd;
113 }
114 else {
115 inc = 0.0;
116 vs->pc = p;
117 vs->p1 = 0.0;
118 }
119
120 for (j=vs->fprd, i=(vs->iprd+1)/2; j--;) {
121 if (vs->p1 == 0.0) {
122 if (vs->gauss)
123 x = (double) nrandom(vs);
124 else
125 x = mseq(vs);
126 }
127 else {
128 if ((vs->pc += 1.0)>=vs->p1) {
129 x = sqrt (vs->p1);
130 vs->pc = vs->pc - vs->p1;
131 }
132 else
133 x = 0.0;
134 }
135
136 x *= exp(vs->c[0]);
137
138 x = mlsadf(x, vs->c, m, a, vs->pd, vs->d1, vs);
139 xs = (short) x;
140
141 fwrite(&xs, sizeof(short), 1, rawfp);
142
143 fflush(stdoutstdout);
144
145 if (!--i) {
146 vs->p1 += inc;
147 for (k=0;k<=m;k++) vs->c[k] += vs->cinc[k];
148 i = vs->iprd;
149 }
150 }
151
152 vs->p1 = p;
153 movem(vs->cc,vs->c,m+1);
154}
155
156double mlsafir (double x, double *b, int m, double a, double *d)
157{
158 double y = 0.0;
159 double aa;
160 register int i;
161
162 aa = 1 - a*a;
163
164 d[0] = x;
165 d[1] = aa*d[0] + a*d[1];
166
167 for (i=2; i<=m; i++) {
168 d[i] = d[i] + a*(d[i+1]-d[i-1]);
169 y += d[i]*b[i];
170 }
171
172 for (i=m+1; i>1; i--) d[i] = d[i-1];
173
174 return (y);
175}
176
177double mlsadf1(double x, double *b, int m, double a, int pd, double *d, VocoderSetup *vs)
178{
179 double v, out = 0.0, *pt, aa;
180 register int i;
181
182 aa = 1 - a*a;
183 pt = &d[pd+1];
184
185 for (i=pd; i>=1; i--) {
186 d[i] = aa*pt[i-1] + a*d[i];
187 pt[i] = d[i] * b[1];
188 v = pt[i] * vs->ppade[i];
189
190 x += (1 & i) ? v : -v;
191 out += v;
192 }
193
194 pt[0] = x;
195 out += x;
196
197 return(out);
198}
199
200double mlsadf2 (double x, double *b, int m, double a, int pd, double *d, VocoderSetup *vs)
201{
202 double v, out = 0.0, *pt, aa;
203 register int i;
204
205 aa = 1 - a*a;
Value stored to 'aa' is never read
206 pt = &d[pd * (m+2)];
207
208 for (i=pd; i>=1; i--) {
209 pt[i] = mlsafir (pt[i-1], b, m, a, &d[(i-1)*(m+2)]);
210 v = pt[i] * vs->ppade[i];
211
212 x += (1&i) ? v : -v;
213 out += v;
214 }
215
216 pt[0] = x;
217 out += x;
218
219 return(out);
220}
221
222double mlsadf(double x, double *b, int m, double a, int pd, double *d, VocoderSetup *vs)
223{
224
225 vs->ppade = &(vs->pade[pd*(pd+1)/2]);
226
227 x = mlsadf1 (x, b, m, a, pd, d, vs);
228 x = mlsadf2 (x, b, m, a, pd, &d[2*(pd+1)], vs);
229
230 return (x);
231}
232
233double nrandom (VocoderSetup *vs)
234{
235 unsigned long rr;
236 if (vs->sw == 0) {
237 vs->sw = 1;
238 do {
239 rr = vs->next;
240 vs->r1 = 2 * rnd(&rr) - 1;
241 vs->r2 = 2 * rnd(&rr) - 1;
242 vs->next = rr;
243 vs->s = vs->r1 * vs->r1 + vs->r2 * vs->r2;
244 } while (vs->s > 1 || vs->s == 0);
245
246 vs->s = sqrt (-2 * log(vs->s) / vs->s);
247 return ( vs->r1 * vs->s );
248 }
249 else {
250 vs->sw = 0;
251 return ( vs->r2 * vs->s );
252 }
253}
254
255double rnd (unsigned long *next)
256{
257 double r;
258
259 *next = *next * 1103515245L + 12345;
260 r = (*next / 65536L) % 32768L;
261
262 return ( r/RANDMAX32767 );
263}
264
265unsigned long srnd ( unsigned long seed )
266{
267 return (seed);
268}
269
270
271int mseq (VocoderSetup *vs)
272{
273 register int x0, x28;
274
275 vs->x >>= 1;
276
277 if (vs->x & B00x00000001)
278 x0 = 1;
279 else
280 x0 = -1;
281
282 if (vs->x & B280x10000000)
283 x28 = 1;
284 else
285 x28 = -1;
286
287 if (x0 + x28)
288 vs->x &= B31_0x7fffffff;
289 else
290 vs->x |= B310x80000000;
291
292 return (x0);
293}
294
295void mc2b( float *mc, double *b, int m, double a)
296{
297 b[m] = mc[m];
298
299 for (m--; m>=0; m--)
300 b[m] = mc[m] - a * b[m+1];
301}
302
303/* -------------------- End of "vocoder.c" -------------------- */