1 /*
2
3 $Log: placea.c,v $
4 Revision 1.2 2007/04/18 14:00:38 rrt
5 Remove $Log tokens and associated log messages (in many files, several
6 copies of every log message were being written) and lots of warnings.
7
8 Revision 1.1 2007/04/16 21:57:47 rrt
9 LPC-10 support, documentation still to come; I wanted to land the code
10 before 14.0.0 went into test, and I'll be busy tomorrow.
11
12 Not highly tested either, but it's just a format, doesn't interfere
13 with anything else, and I'll get on that case before we go stable.
14
15 * Revision 1.1 1996/08/19 22:31:07 jaf
16 * Initial revision
17 *
18
19 */
20
21 /* -- translated by f2c (version 19951025).
22 You must link the resulting object file with the libraries:
23 -lf2c -lm (in that order)
24 */
25
26 #include "f2c.h"
27
28 extern int placea_(integer *ipitch, integer *voibuf, integer *obound, integer *af, integer *vwin, integer *awin, integer *ewin, integer *lframe, integer *maxwin);
29
30 /* *********************************************************************** */
31
32 /* PLACEA Version 48 */
33
34 /* Revision 1.1 1996/08/19 22:31:07 jaf
35 * Initial revision
36 * */
37 /* Revision 1.5 1996/03/19 20:41:55 jaf */
38 /* Added some conditions satisfied by the output values in EWIN. */
39
40 /* Revision 1.4 1996/03/19 20:24:17 jaf */
41 /* Added some conditions satisfied by the output values in AWIN. */
42
43 /* Revision 1.3 1996/03/18 21:40:04 jaf */
44 /* Just added a few comments about which array indices of the arguments */
45 /* are used, and mentioning that this subroutine has no local state. */
46
47 /* Revision 1.2 1996/03/13 16:43:09 jaf */
48 /* Comments added explaining that none of the local variables of this */
49 /* subroutine need to be saved from one invocation to the next. */
50
51 /* Revision 1.1 1996/02/07 14:48:31 jaf */
52 /* Initial revision */
53
54
55 /* *********************************************************************** */
56 /* Input: */
57 /* IPITCH */
58 /* VOIBUF */
59 /* Indices (2,AF-2), (1,AF-1), (2,AF-1), (1,AF), and (2,AF) read.*/
60 /* All other indices untouched. */
61 /* OBOUND */
62 /* AF */
63 /* VWIN */
64 /* Indices (1,AF) and (2,AF) read. */
65 /* All other indices untouched. */
66 /* LFRAME */
67 /* MAXWIN */
68 /* Input/Output: */
69 /* AWIN */
70 /* Index (1,AF-1) read. */
71 /* Indices (1,AF) and (2,AF) written, and then read. */
72 /* All other indices untouched. */
73 /* In all cases (except possibly one), the final values will */
74 /* satisfy the condition: AWIN(2,AF)-AWIN(1,AF)+1 = MAXWIN. */
75 /* In that other case, */
76 /* AWIN(1,AF)=VWIN(1,AF) and AWIN(2,AF)=VWIN(2,AF). */
77 /* Output: */
78 /* EWIN */
79 /* Indices (1,AF) and (2,AF) written. */
80 /* All other indices untouched. */
81 /* In all cases, the final values will satisfy the condition: */
82 /* AWIN(1,AF) .LE. EWIN(1,AF) .LE. EWIN(2,AF) .LE. AWIN(2,AF) */
83 /* In other words, the energy window is a sub-window of */
84 /* the analysis window. */
85
86 /* This subroutine has no local state. */
87
placea_(integer * ipitch,integer * voibuf,integer * obound,integer * af,integer * vwin,integer * awin,integer * ewin,integer * lframe,integer * maxwin)88 /* Subroutine */ int placea_(integer *ipitch, integer *voibuf, integer *
89 obound, integer *af, integer *vwin, integer *awin, integer *ewin,
90 integer *lframe, integer *maxwin)
91 {
92 /* System generated locals */
93 real r__1;
94
95 /* Builtin functions */
96 integer i_nint(real *);
97
98 /* Local variables */
99 logical allv, winv;
100 integer i__, j, k, l, hrange;
101 logical ephase;
102 integer lrange;
103
104 /* Arguments */
105 /* Local variables that need not be saved */
106 /* Parameter adjustments */
107 ewin -= 3;
108 awin -= 3;
109 vwin -= 3;
110 --voibuf;
111
112 /* Function Body */
113 lrange = (*af - 2) * *lframe + 1;
114 hrange = *af * *lframe;
115 /* Place the Analysis window based on the voicing window */
116 /* placement, onsets, tentative voicing decision, and pitch. */
117
118 /* Case 1: Sustained Voiced Speech */
119 /* If the five most recent voicing decisions are */
120 /* voiced, then the window is placed phase-synchronously with the */
121 /* previous window, as close to the present voicing window if possible.
122 */
123 /* If onsets bound the voicing window, then preference is given to */
124 /* a phase-synchronous placement which does not overlap these onsets. */
125
126 /* Case 2: Voiced Transition */
127 /* If at least one voicing decision in AF is voicied, and there are no
128 */
129 /* onsets, then the window is placed as in case 1. */
130
131 /* Case 3: Unvoiced Speech or Onsets */
132 /* If both voicing decisions in AF are unvoiced, or there are onsets, */
133 /* then the window is placed coincident with the voicing window. */
134
135 /* Note: During phase-synchronous placement of windows, the length */
136 /* is not altered from MAXWIN, since this would defeat the purpose */
137 /* of phase-synchronous placement. */
138 /* Check for case 1 and case 2 */
139 allv = voibuf[((*af - 2) << 1) + 2] == 1;
140 allv = allv && voibuf[((*af - 1) << 1) + 1] == 1;
141 allv = allv && voibuf[((*af - 1) << 1) + 2] == 1;
142 allv = allv && voibuf[(*af << 1) + 1] == 1;
143 allv = allv && voibuf[(*af << 1) + 2] == 1;
144 winv = voibuf[(*af << 1) + 1] == 1 || voibuf[(*af << 1) + 2] == 1;
145 if (allv || (winv && *obound == 0)) {
146 /* APHASE: Phase synchronous window placement. */
147 /* Get minimum lower index of the window. */
148 i__ = (lrange + *ipitch - 1 - awin[((*af - 1) << 1) + 1]) / *ipitch;
149 i__ *= *ipitch;
150 i__ += awin[((*af - 1) << 1) + 1];
151 /* L = the actual length of this frame's analysis window. */
152 l = *maxwin;
153 /* Calculate the location where a perfectly centered window would star
154 t. */
155 k = (vwin[(*af << 1) + 1] + vwin[(*af << 1) + 2] + 1 - l) / 2;
156 /* Choose the actual location to be the pitch multiple closest to this
157 . */
158 r__1 = (real) (k - i__) / *ipitch;
159 awin[(*af << 1) + 1] = i__ + i_nint(&r__1) * *ipitch;
160 awin[(*af << 1) + 2] = awin[(*af << 1) + 1] + l - 1;
161 /* If there is an onset bounding the right of the voicing window and t
162 he */
163 /* analysis window overlaps that, then move the analysis window backwa
164 rd */
165 /* to avoid this onset. */
166 if (*obound >= 2 && awin[(*af << 1) + 2] > vwin[(*af << 1) + 2]) {
167 awin[(*af << 1) + 1] -= *ipitch;
168 awin[(*af << 1) + 2] -= *ipitch;
169 }
170 /* Similarly for the left of the voicing window. */
171 if ((*obound == 1 || *obound == 3) && awin[(*af << 1) + 1] < vwin[(*
172 af << 1) + 1]) {
173 awin[(*af << 1) + 1] += *ipitch;
174 awin[(*af << 1) + 2] += *ipitch;
175 }
176 /* If this placement puts the analysis window above HRANGE, then */
177 /* move it backward an integer number of pitch periods. */
178 while(awin[(*af << 1) + 2] > hrange) {
179 awin[(*af << 1) + 1] -= *ipitch;
180 awin[(*af << 1) + 2] -= *ipitch;
181 }
182 /* Similarly if the placement puts the analysis window below LRANGE.
183 */
184 while(awin[(*af << 1) + 1] < lrange) {
185 awin[(*af << 1) + 1] += *ipitch;
186 awin[(*af << 1) + 2] += *ipitch;
187 }
188 /* Make Energy window be phase-synchronous. */
189 ephase = TRUE_;
190 /* Case 3 */
191 } else {
192 awin[(*af << 1) + 1] = vwin[(*af << 1) + 1];
193 awin[(*af << 1) + 2] = vwin[(*af << 1) + 2];
194 ephase = FALSE_;
195 }
196 /* RMS is computed over an integer number of pitch periods in the analysis
197 */
198 /*window. When it is not placed phase-synchronously, it is placed as clos
199 e*/
200 /* as possible to onsets. */
201 j = (awin[(*af << 1) + 2] - awin[(*af << 1) + 1] + 1) / *ipitch * *ipitch;
202 if (j == 0 || ! winv) {
203 ewin[(*af << 1) + 1] = vwin[(*af << 1) + 1];
204 ewin[(*af << 1) + 2] = vwin[(*af << 1) + 2];
205 } else if (! ephase && *obound == 2) {
206 ewin[(*af << 1) + 1] = awin[(*af << 1) + 2] - j + 1;
207 ewin[(*af << 1) + 2] = awin[(*af << 1) + 2];
208 } else {
209 ewin[(*af << 1) + 1] = awin[(*af << 1) + 1];
210 ewin[(*af << 1) + 2] = awin[(*af << 1) + 1] + j - 1;
211 }
212 return 0;
213 } /* placea_ */
214
215