lpc10/placev.c

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