1 /*
2
3 * Revision 1.2 1996/08/20 20:35:41 jaf
4 * Added functions for allocating and initializing lpc10_encoder_state
5 * and lpc10_decoder_state structures.
6 *
7 * Revision 1.1 1996/08/19 22:31:40 jaf
8 * Initial revision
9 *
10
11 */
12
13 /* -- translated by f2c (version 19951025).
14 You must link the resulting object file with the libraries:
15 -lf2c -lm (in that order)
16 */
17
18 #include <stdlib.h>
19
20 #include "f2c.h"
21
22 extern int lpcini_(void);
23
24 /* Common Block Declarations */
25
26 struct {
27 integer order, lframe;
28 logical corrp;
29 } contrl_;
30
31 #define contrl_1 contrl_
32
33 /* ***************************************************************** */
34
35 /*
36 * Revision 1.2 1996/08/20 20:35:41 jaf
37 * Added functions for allocating and initializing lpc10_encoder_state
38 * and lpc10_decoder_state structures.
39 *
40 * Revision 1.1 1996/08/19 22:31:40 jaf
41 * Initial revision
42 * */
43 /* Revision 1.1 1996/03/28 00:04:05 jaf */
44 /* Initial revision */
45
46
47 /* ***************************************************************** */
48
49 /* Initialize COMMON block variables used by LPC-10 encoder and decoder, */
50 /* and call initialization routines for both of them. */
51
lpcini_(void)52 /* Subroutine */ int lpcini_(void)
53 {
54
55 /* LPC Processing control variables: */
56
57 /* *** Read-only: initialized in setup */
58
59 /* Files for Speech, Parameter, and Bitstream Input & Output, */
60 /* and message and debug outputs. */
61
62 /* Here are the only files which use these variables: */
63
64 /* lpcsim.f setup.f trans.f error.f vqsetup.f */
65
66 /* Many files which use fdebug are not listed, since it is only used in */
67 /* those other files conditionally, to print trace statements. */
68 /* integer fsi, fso, fpi, fpo, fbi, fbo, pbin, fmsg, fdebug */
69 /* LPC order, Frame size, Quantization rate, Bits per frame, */
70 /* Error correction */
71 /* Subroutine SETUP is the only place where order is assigned a value, */
72 /* and that value is 10. It could increase efficiency 1% or so to */
73 /* declare order as a constant (i.e., a Fortran PARAMETER) instead of as
74 */
75 /* a variable in a COMMON block, since it is used in many places in the */
76 /* core of the coding and decoding routines. Actually, I take that back.
77 */
78 /* At least when compiling with f2c, the upper bound of DO loops is */
79 /* stored in a local variable before the DO loop begins, and then that is
80 */
81 /* compared against on each iteration. */
82 /* Similarly for lframe, which is given a value of MAXFRM in SETUP. */
83 /* Similarly for quant, which is given a value of 2400 in SETUP. quant */
84 /* is used in only a few places, and never in the core coding and */
85 /* decoding routines, so it could be eliminated entirely. */
86 /* nbits is similar to quant, and is given a value of 54 in SETUP. */
87 /* corrp is given a value of .TRUE. in SETUP, and is only used in the */
88 /* subroutines ENCODE and DECODE. It doesn't affect the speed of the */
89 /* coder significantly whether it is .TRUE. or .FALSE., or whether it is
90 */
91 /* a constant or a variable, since it is only examined once per frame. */
92 /* Leaving it as a variable that is set to .TRUE. seems like a good */
93 /* idea, since it does enable some error-correction capability for */
94 /* unvoiced frames, with no change in the coding rate, and no noticeable
95 */
96 /* quality difference in the decoded speech. */
97 /* integer quant, nbits */
98 /* *** Read/write: variables for debugging, not needed for LPC algorithm
99 */
100
101 /* Current frame, Unstable frames, Output clip count, Max onset buffer,
102 */
103 /* Debug listing detail level, Line count on listing page */
104
105 /* nframe is not needed for an embedded LPC10 at all. */
106 /* nunsfm is initialized to 0 in SETUP, and incremented in subroutine */
107 /* ERROR, which is only called from RCCHK. When LPC10 is embedded into */
108 /* an application, I would recommend removing the call to ERROR in RCCHK,
109 */
110 /* and remove ERROR and nunsfm completely. */
111 /* iclip is initialized to 0 in SETUP, and incremented in entry SWRITE in
112 */
113 /* sread.f. When LPC10 is embedded into an application, one might want */
114 /* to cause it to be incremented in a routine that takes the output of */
115 /* SYNTHS and sends it to an audio device. It could be optionally */
116 /* displayed, for those that might want to know what it is. */
117 /* maxosp is never initialized to 0 in SETUP, although it probably should
118 */
119 /* be, and it is updated in subroutine ANALYS. I doubt that its value */
120 /* would be of much interest to an application in which LPC10 is */
121 /* embedded. */
122 /* listl and lincnt are not needed for an embedded LPC10 at all. */
123 /* integer nframe, nunsfm, iclip, maxosp, listl, lincnt */
124 /* common /contrl/ fsi, fso, fpi, fpo, fbi, fbo, pbin, fmsg, fdebug */
125 /* common /contrl/ quant, nbits */
126 /* common /contrl/ nframe, nunsfm, iclip, maxosp, listl, lincnt */
127 contrl_1.order = 10;
128 contrl_1.lframe = 180;
129 contrl_1.corrp = TRUE_;
130 return 0;
131 } /* lpcini_ */
132
133
134
135 /* Allocate memory for, and initialize, the state that needs to be
136 kept from encoding one frame to the next for a single
137 LPC-10-compressed audio stream. Return 0 if malloc fails,
138 otherwise return pointer to new structure. */
139
140 struct lpc10_encoder_state *
create_lpc10_encoder_state(void)141 create_lpc10_encoder_state(void)
142 {
143 struct lpc10_encoder_state *st;
144
145 st = (struct lpc10_encoder_state *)
146 malloc((unsigned) sizeof (struct lpc10_encoder_state));
147 if (st != 0) {
148 init_lpc10_encoder_state(st);
149 }
150 return (st);
151 }
152
153
154
init_lpc10_encoder_state(struct lpc10_encoder_state * st)155 void init_lpc10_encoder_state(struct lpc10_encoder_state *st)
156 {
157 int i;
158
159 lpcini_();
160
161 /* State used only by function hp100 */
162 st->z11 = 0.0f;
163 st->z21 = 0.0f;
164 st->z12 = 0.0f;
165 st->z22 = 0.0f;
166
167 /* State used by function analys */
168 for (i = 0; i < 540; i++) {
169 st->inbuf[i] = 0.0f;
170 st->pebuf[i] = 0.0f;
171 }
172 for (i = 0; i < 696; i++) {
173 st->lpbuf[i] = 0.0f;
174 }
175 for (i = 0; i < 312; i++) {
176 st->ivbuf[i] = 0.0f;
177 }
178 st->bias = 0.0f;
179 /* integer osbuf[10]; no initial value necessary */
180 st->osptr = 1;
181 for (i = 0; i < 3; i++) {
182 st->obound[i] = 0;
183 }
184 st->vwin[4] = 307;
185 st->vwin[5] = 462;
186 st->awin[4] = 307;
187 st->awin[5] = 462;
188 for (i = 0; i < 8; i++) {
189 st->voibuf[i] = 0;
190 }
191 for (i = 0; i < 3; i++) {
192 st->rmsbuf[i] = 0.0f;
193 }
194 for (i = 0; i < 30; i++) {
195 st->rcbuf[i] = 0.0f;
196 }
197 st->zpre = 0.0f;
198
199
200 /* State used by function onset */
201 st->n = 0.0f;
202 st->d__ = 1.0f;
203 /* real fpc; no initial value necessary */
204 for (i = 0; i < 16; i++) {
205 st->l2buf[i] = 0.0f;
206 }
207 st->l2sum1 = 0.0f;
208 st->l2ptr1 = 1;
209 st->l2ptr2 = 9;
210 /* integer lasti; no initial value necessary */
211 st->hyst = FALSE_;
212
213 /* State used by function voicin */
214 st->dither = 20.0f;
215 st->maxmin = 0.0f;
216 for (i = 0; i < 6; i++) {
217 st->voice[i] = 0.0f;
218 }
219 st->lbve = 3000;
220 st->fbve = 3000;
221 st->fbue = 187;
222 st->ofbue = 187;
223 st->sfbue = 187;
224 st->lbue = 93;
225 st->olbue = 93;
226 st->slbue = 93;
227 st->snr = (real) (st->fbve / st->fbue << 6);
228
229 /* State used by function dyptrk */
230 for (i = 0; i < 60; i++) {
231 st->s[i] = 0.0f;
232 }
233 for (i = 0; i < 120; i++) {
234 st->p[i] = 0;
235 }
236 st->ipoint = 0;
237 st->alphax = 0.0f;
238
239 /* State used by function chanwr */
240 st->isync = 0;
241
242 }
243
244
245
246 /* Allocate memory for, and initialize, the state that needs to be
247 kept from decoding one frame to the next for a single
248 LPC-10-compressed audio stream. Return 0 if malloc fails,
249 otherwise return pointer to new structure. */
250
251 struct lpc10_decoder_state *
create_lpc10_decoder_state(void)252 create_lpc10_decoder_state(void)
253 {
254 struct lpc10_decoder_state *st;
255
256 st = (struct lpc10_decoder_state *)
257 malloc((unsigned) sizeof (struct lpc10_decoder_state));
258 if (st != 0) {
259 init_lpc10_decoder_state(st);
260 }
261 return (st);
262 }
263
264
265
init_lpc10_decoder_state(struct lpc10_decoder_state * st)266 void init_lpc10_decoder_state(struct lpc10_decoder_state *st)
267 {
268 int i;
269
270 lpcini_();
271
272 /* State used by function decode */
273 st->iptold = 60;
274 st->first = TRUE_;
275 st->ivp2h = 0;
276 st->iovoic = 0;
277 st->iavgp = 60;
278 st->erate = 0;
279 for (i = 0; i < 30; i++) {
280 st->drc[i] = 0;
281 }
282 for (i = 0; i < 3; i++) {
283 st->dpit[i] = 0;
284 st->drms[i] = 0;
285 }
286
287 /* State used by function synths */
288 for (i = 0; i < 360; i++) {
289 st->buf[i] = 0.0f;
290 }
291 st->buflen = 180;
292
293 /* State used by function pitsyn */
294 /* ivoico; no initial value necessary as long as first_pitsyn is initially TRUE_ */
295 /* ipito; no initial value necessary as long as first_pitsyn is initially TRUE_ */
296 st->rmso = 1.0f;
297 /* rco[10]; no initial value necessary as long as first_pitsyn is initially TRUE_ */
298 /* integer jsamp; no initial value necessary as long as first_pitsyn is initially TRUE_ */
299 st->first_pitsyn = TRUE_;
300
301 /* State used by function bsynz */
302 st->ipo = 0;
303 for (i = 0; i < 166; i++) {
304 st->exc[i] = 0.0f;
305 st->exc2[i] = 0.0f;
306 }
307 st->lpi1 = 0.0f;
308 st->lpi2 = 0.0f;
309 st->lpi3 = 0.0f;
310 st->hpi1 = 0.0f;
311 st->hpi2 = 0.0f;
312 st->hpi3 = 0.0f;
313 st->rmso_bsynz = 0.0f;
314
315 /* State used by function random */
316 st->j = 2;
317 st->k = 5;
318 st->y[0] = (shortint) -21161;
319 st->y[1] = (shortint) -8478;
320 st->y[2] = (shortint) 30892;
321 st->y[3] = (shortint) -10216;
322 st->y[4] = (shortint) 16950;
323
324 /* State used by function deemp */
325 st->dei1 = 0.0f;
326 st->dei2 = 0.0f;
327 st->deo1 = 0.0f;
328 st->deo2 = 0.0f;
329 st->deo3 = 0.0f;
330 }
331