eae5fcdec1d2355992fd407aa53d2f27f25547d9
[asterisk/asterisk.git] / codecs / lpc10 / placev.c
1 /*
2
3 $Log$
4 Revision 1.14  2003/02/12 13:59:15  matteo
5 mer feb 12 14:56:57 CET 2003
6
7 Revision 1.1.1.1  2003/02/12 13:59:15  matteo
8 mer feb 12 14:56:57 CET 2003
9
10 Revision 1.3  2001/04/12 21:27:53  markh
11 app_record now supports wildcards of sort so your output file is not overwritten every time it's run.  File.h got a documentation update on the ast_fileexists to include the return call.  Watch out for the placea.c placev.c code, it's updates have not been tested yet.  Just a few parenthesis to make it compile nicer on newer gcc versions with all the -W flags set.
12
13 Revision 1.2  2000/01/05 08:20:39  markster
14 Some OSS fixes and a few lpc changes to make it actually work
15
16  * Revision 1.1  1996/08/19  22:31:02  jaf
17  * Initial revision
18  *
19
20 */
21
22 #ifdef P_R_O_T_O_T_Y_P_E_S
23 extern int placev_(integer *osbuf, integer *osptr, integer *oslen, integer *obound, integer *vwin, integer *af, integer *lframe, integer *minwin, integer *maxwin, integer *dvwinl, integer *dvwinh);
24 #endif
25
26 /*  -- translated by f2c (version 19951025).
27    You must link the resulting object file with the libraries:
28         -lf2c -lm   (in that order)
29 */
30
31 #include "f2c.h"
32
33 /* ****************************************************************** */
34
35 /*      PLACEV Version 48 */
36
37 /* $Log$
38  * Revision 1.14  2003/02/12 13:59:15  matteo
39  * mer feb 12 14:56:57 CET 2003
40  *
41 /* Revision 1.1.1.1  2003/02/12 13:59:15  matteo
42 /* mer feb 12 14:56:57 CET 2003
43 /*
44 /* Revision 1.3  2001/04/12 21:27:53  markh
45 /* app_record now supports wildcards of sort so your output file is not overwritten every time it's run.  File.h got a documentation update on the ast_fileexists to include the return call.  Watch out for the placea.c placev.c code, it's updates have not been tested yet.  Just a few parenthesis to make it compile nicer on newer gcc versions with all the -W flags set.
46 /*
47 /* Revision 1.2  2000/01/05 08:20:39  markster
48 /* Some OSS fixes and a few lpc changes to make it actually work
49 /*
50  * Revision 1.1  1996/08/19  22:31:02  jaf
51  * Initial revision
52  * */
53 /* Revision 1.6  1996/03/19  20:42:19  jaf */
54 /* Added some conditions satisfied by the output values in VWIN. */
55
56 /* Revision 1.5  1996/03/19  18:37:56  jaf */
57 /* Strengthened the specification of which indices of VWIN are read and */
58 /* written. */
59
60 /* Revision 1.4  1996/03/15  16:38:33  jaf */
61 /* One tiny comment added. */
62
63 /* Revision 1.3  1996/03/15  16:36:13  jaf */
64 /* Added comments giving In/Out status of arguments. */
65
66 /* Revision 1.2  1996/03/12  23:56:01  jaf */
67 /* Comments added explaining that none of the local variables of this */
68 /* subroutine need to be saved from one invocation to the next. */
69
70 /* Revision 1.1  1996/02/07 14:48:39  jaf */
71 /* Initial revision */
72
73
74 /* ****************************************************************** */
75
76 /* Input: */
77 /*  OSBUF       Buffer which holds sorted indexes of onsets */
78 /*              I believe that only indices 1 through OSPTR-1 can be read. */
79 /*  OSLEN */
80 /*  OSPTR       Free pointer into OSBUF */
81 /*  AF */
82 /*  LFRAME */
83 /*  MINWIN */
84 /*  MAXWIN */
85 /*  DVWINL */
86 /*  DVWINH      (This argument is never used.  Should it be?) */
87 /* Input/Output: */
88 /*  VWIN                Buffer of Voicing Window Positions (Modified) */
89 /*              Index (2,AF-1) is read. */
90 /*              Indices (1,AF) and (2,AF) are written, */
91 /*              and then possibly read. */
92 /*              All other indices are unused. */
93 /*              In all cases, the final values will satsify the condition:*/
94 /*               VWIN(2,AF)-VWIN(1,AF)+1 .LE. MAXWIN */
95 /*               I'm not certain yet, but they may also satisfy: */
96 /*               MINWIN .LE. VWIN(2,AF)-VWIN(1,AF)+1 */
97 /* Output: */
98 /*  OBOUND      This variable is set by this procedure and used */
99 /*              in placing analysis windows (PLACEA).  Bit 1 */
100 /*              indicates whether an onset bounds the left side */
101 /*              of the voicing window, and bit 2 indicates whether */
102 /*              an onset bounds the right side of the voicing window. */
103
104 /* This subroutine has no local state. */
105
106 /* Subroutine */ int placev_(integer *osbuf, integer *osptr, integer *oslen, 
107         integer *obound, integer *vwin, integer *af, integer *lframe, integer 
108         *minwin, integer *maxwin, integer *dvwinl, integer *dvwinh)
109 {
110     /* System generated locals */
111     integer i__1, i__2;
112
113     /* Local variables */
114     logical crit;
115     integer i__, q, osptr1, hrange, lrange;
116
117 /*       Arguments */
118 /*       Local variables that need not be saved */
119 /*   Variables */
120 /*    LRANGE, HRANGE  Range in which window is placed */
121 /*    OSPTR1     OSPTR excluding samples in 3F */
122 /*       Local state */
123 /*       None */
124 /*   Voicing Window Placement */
125
126 /*         __________________ __________________ ______________ */
127 /*        |                  |                  | */
128 /*        |        1F        |        2F        |        3F ... */
129 /*        |__________________|__________________|______________ */
130
131 /*    Previous | */
132 /*      Window | */
133 /*  ...________| */
134
135 /*             |                                | */
136 /*      ------>| This window's placement range  |<------ */
137 /*             |                                | */
138
139 /*   There are three cases.  Note that these are different from those */
140 /*   given in the LPC-10e phase 1 report. */
141
142 /*   1.  If there are no onsets in this range, then the voicing window */
143 /*   is centered in the pitch window.  If such a placement is not within 
144 */
145 /*   the window's placement range, then the window is placed in the left- 
146 */
147 /*   most portion of the placement range.  Its length is always MAXWIN. */
148
149 /*   2.  If the first onset is in 2F and there is sufficient room to place
150  */
151 /*   the window immediately before this onset, then the window is placed 
152 */
153 /*   there, and its length is set to the maximum possible under these */
154 /*   constraints. */
155
156 /*      "Critical Region Exception":  If there is another onset in 2F */
157 /*      such that a window can be placed between the two onsets, the */
158 /*      window is placed there (ie, as in case 3). */
159
160 /*   3.  Otherwise, the window is placed immediately after the onset.  The
161  */
162 /*   window's length */
163 /*  is the longest length that can fit in the range under these constraint
164 s,*/
165 /*  except that the window may be shortened even further to avoid overlapp
166 ing*/
167 /*  other onsets in the placement range.  In any case, the window's length
168 */
169 /*   is at least MINWIN. */
170
171 /*   Note that the values of MINWIN and LFRAME must be chosen such */
172 /*   that case 2 = false implies case 3 = true.   This means that */
173 /*   MINWIN <= LFRAME/2.  If this were not the case, then a fourth case */
174 /*   would have to be added for when the window cannot fit either before 
175 */
176 /*   or after the onset. */
177
178 /*   Note also that onsets which weren't in 2F last time may be in 1F this
179  */
180 /*  time, due to the filter delays in computing onsets.  The result is tha
181 t*/
182 /*   occasionally a voicing window will overlap that onset.  The only way 
183 */
184 /*   to circumvent this problem is to add more delay in processing input 
185 */
186 /*   speech.  In the trade-off between delay and window-placement, window 
187 */
188 /*   placement lost. */
189 /* Compute the placement range */
190     /* Parameter adjustments */
191     --osbuf;
192     vwin -= 3;
193
194     /* Function Body */
195 /* Computing MAX */
196     i__1 = vwin[((*af - 1) << 1) + 2] + 1, i__2 = (*af - 2) * *lframe + 1;
197     lrange = max(i__1,i__2);
198     hrange = *af * *lframe;
199 /* Compute OSPTR1, so the following code only looks at relevant onsets. */
200     for (osptr1 = *osptr - 1; osptr1 >= 1; --osptr1) {
201         if (osbuf[osptr1] <= hrange) {
202             goto L90;
203         }
204     }
205 L90:
206     ++osptr1;
207 /* Check for case 1 first (fast case): */
208     if (osptr1 <= 1 || osbuf[osptr1 - 1] < lrange) {
209 /* Computing MAX */
210         i__1 = vwin[((*af - 1) << 1) + 2] + 1;
211         vwin[(*af << 1) + 1] = max(i__1,*dvwinl);
212         vwin[(*af << 1) + 2] = vwin[(*af << 1) + 1] + *maxwin - 1;
213         *obound = 0;
214     } else {
215 /* Search backward in OSBUF for first onset in range. */
216 /* This code relies on the above check being performed first. */
217         for (q = osptr1 - 1; q >= 1; --q) {
218             if (osbuf[q] < lrange) {
219                 goto L100;
220             }
221         }
222 L100:
223         ++q;
224 /* Check for case 2 (placement before onset): */
225 /* Check for critical region exception: */
226         i__1 = osptr1 - 1;
227         for (i__ = q + 1; i__ <= i__1; ++i__) {
228             if (osbuf[i__] - osbuf[q] >= *minwin) {
229                 crit = TRUE_;
230                 goto L105;
231             }
232         }
233         crit = FALSE_;
234 L105:
235 /* Computing MAX */
236         i__1 = (*af - 1) * *lframe, i__2 = lrange + *minwin - 1;
237         if (! crit && osbuf[q] > max(i__1,i__2)) {
238             vwin[(*af << 1) + 2] = osbuf[q] - 1;
239 /* Computing MAX */
240             i__1 = lrange, i__2 = vwin[(*af << 1) + 2] - *maxwin + 1;
241             vwin[(*af << 1) + 1] = max(i__1,i__2);
242             *obound = 2;
243 /* Case 3 (placement after onset) */
244         } else {
245             vwin[(*af << 1) + 1] = osbuf[q];
246 L110:
247             ++q;
248             if (q >= osptr1) {
249                 goto L120;
250             }
251             if (osbuf[q] > vwin[(*af << 1) + 1] + *maxwin) {
252                 goto L120;
253             }
254             if (osbuf[q] < vwin[(*af << 1) + 1] + *minwin) {
255                 goto L110;
256             }
257             vwin[(*af << 1) + 2] = osbuf[q] - 1;
258             *obound = 3;
259             return 0;
260 L120:
261 /* Computing MIN */
262             i__1 = vwin[(*af << 1) + 1] + *maxwin - 1;
263             vwin[(*af << 1) + 2] = min(i__1,hrange);
264             *obound = 1;
265         }
266     }
267     return 0;
268 } /* placev_ */
269