lpc10/vparms.c

   1 /*
   2
   3  * Revision 1.1  1996/08/19  22:30:04  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 vparms_(integer *vwin, real *inbuf, real *lpbuf, integer *buflim, integer *half, real *dither, integer *mintau, integer *zc, integer *lbe, integer *fbe, real *qs, real *rc1, real *ar_b__, real *ar_f__);
  17
  18 /* Table of constant values */
  19
  20 static real c_b2 = 1.f;
  21
  22 /* ********************************************************************* */
  23
  24 /*      VPARMS Version 50 */
  25
  26 /*
  27  * Revision 1.1  1996/08/19  22:30:04  jaf
  28  * Initial revision
  29  * */
  30 /* Revision 1.6  1996/03/29  18:01:16  jaf */
  31 /* Added some more comments about the range of INBUF and LPBUF that can */
  32 /* be read.  Note that it is possible for index VWIN(2)+1 to be read from */
  33 /* INBUF, which might be outside of its defined range, although that will */
  34 /* require more careful checking. */
  35
  36 /* Revision 1.5  1996/03/19  00:02:02  jaf */
  37 /* I just noticed that the argument DITHER is modified inside of this */
  38 /* subroutine.  Comments were added explaining the possible final values. */
  39
  40 /* Revision 1.4  1996/03/18  22:22:59  jaf */
  41 /* Finishing the job I said I did with the last check-in comments. */
  42
  43 /* Revision 1.3  1996/03/18  22:22:17  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  15:02:58  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:50:42  jaf */
  52 /* Initial revision */
  53
  54
  55 /* ********************************************************************* */
  56
  57 /*  Calculate voicing parameters: */
  58
  59 /* Input: */
  60 /*  VWIN   - Voicing window limits */
  61 /*           Indices 1 through 2 read. */
  62 /*  INBUF  - Input speech buffer */
  63 /*           Indices START-1 through STOP read, */
  64 /*          where START and STOP are defined in the code (only written once).
  65 */
  66 /*           Note that STOP can be as large as VWIN(2)+1 ! */
  67 /*  LPBUF  - Low pass filtered speech */
  68 /*           Indices START-MINTAU through STOP+MINTAU read, */
  69 /*          where START and STOP are defined in the code (only written once).
  70 */
  71 /*  BUFLIM - Array bounds for INBUF and LPBUF */
  72 /*           Indices 1 through 4 read. */
  73 /*  HALF   - Half frame (1 or 2) */
  74 /*  MINTAU - Lag corresponding to minimum AMDF value (pitch estimate) */
  75 /* Input/Output: */
  76 /*  DITHER - Zero crossing threshold */
  77 /*           The resulting value might be the negation of the input */
  78 /*           value.  It might always be the same as the input value, */
  79 /*           if the DO loop below always executes an even number of times. */
  80 /* Output: (all of them are written on every call) */
  81 /*  ZC     - Zero crossing rate */
  82 /*  LBE    - Low band energy (sum of magnitudes - SM) */
  83 /*  FBE    - Full band energy (SM) */
  84 /*  QS     - Ratio of 6 dB/oct preemphasized energy to full band energy */
  85 /*  RC1    - First reflection coefficient */
  86 /*  AR_B   - Product of the causal forward and reverse pitch */
  87 /*           prediction gains */
  88 /*  AR_F   - Product of the noncausal forward and reverse pitch */
  89 /*           prediction gains */
  90 /* Internal: */
  91 /*  OLDSGN - Previous sign of dithered signal */
  92 /*  VLEN   - Length of voicing window */
  93 /*  START  - Lower address of current half of voicing window */
  94 /*  STOP   - Upper address of current half of voicing window */
  95 /*  E_0    - Energy of LPF speech (sum of squares - SS) */
  96 /*  E_B    - Energy of LPF speech backward one pitch period (SS) */
  97 /*  E_F    - Energy of LPF speech forward one pitch period (SS) */
  98 /*  R_B    - Autocovariance of LPF speech backward one pitch period */
  99 /*  R_F    - Autocovariance of LPF speech forward one pitch period */
 100 /*  LP_RMS - Energy of LPF speech (sum of magnitudes - SM) */
 101 /*  AP_RMS - Energy of all-pass speech (SM) */
 102 /*  E_PRE  - Energy of 6dB preemphasized speech (SM) */
 103 /*  E0AP   - Energy of all-pass speech (SS) */
 104
 105 /* This subroutine has no local state. */
 106
 107 /* Subroutine */ int vparms_(integer *vwin, real *inbuf, real *lpbuf, integer
 108         *buflim, integer *half, real *dither, integer *mintau, integer *zc,
 109         integer *lbe, integer *fbe, real *qs, real *rc1, real *ar_b__, real *
 110         ar_f__)
 111 {
 112     /* System generated locals */
 113     integer inbuf_offset, lpbuf_offset, i__1;
 114     real r__1, r__2;
 115
 116     /* Builtin functions */
 117     double r_sign(real *, real *);
 118     integer i_nint(real *);
 119
 120     /* Local variables */
 121     integer vlen, stop, i__;
 122     real e_pre__;
 123     integer start;
 124     real ap_rms__, e_0__, oldsgn, lp_rms__, e_b__, e_f__, r_b__, r_f__, e0ap;
 125
 126 /*       Arguments */
 127 /*       Local variables that need not be saved */
 128 /*   Calculate zero crossings (ZC) and several energy and correlation */
 129 /*   measures on low band and full band speech.  Each measure is taken */
 130 /*   over either the first or the second half of the voicing window, */
 131 /*   depending on the variable HALF. */
 132     /* Parameter adjustments */
 133     --vwin;
 134     --buflim;
 135     lpbuf_offset = buflim[3];
 136     lpbuf -= lpbuf_offset;
 137     inbuf_offset = buflim[1];
 138     inbuf -= inbuf_offset;
 139
 140     /* Function Body */
 141     lp_rms__ = 0.f;
 142     ap_rms__ = 0.f;
 143     e_pre__ = 0.f;
 144     e0ap = 0.f;
 145     *rc1 = 0.f;
 146     e_0__ = 0.f;
 147     e_b__ = 0.f;
 148     e_f__ = 0.f;
 149     r_f__ = 0.f;
 150     r_b__ = 0.f;
 151     *zc = 0;
 152     vlen = vwin[2] - vwin[1] + 1;
 153     start = vwin[1] + (*half - 1) * vlen / 2 + 1;
 154     stop = start + vlen / 2 - 1;
 155
 156 /* I'll use the symbol HVL in the table below to represent the value */
 157 /* VLEN/2.  Note that if VLEN is odd, then HVL should be rounded down, */
 158 /* i.e., HVL = (VLEN-1)/2. */
 159
 160 /* HALF  START          STOP */
 161
 162 /* 1     VWIN(1)+1      VWIN(1)+HVL */
 163 /* 2     VWIN(1)+HVL+1  VWIN(1)+2*HVL */
 164
 165 /* Note that if VLEN is even and HALF is 2, then STOP will be */
 166 /* VWIN(1)+VLEN = VWIN(2)+1.  That could be bad, if that index of INBUF */
 167 /* is undefined. */
 168
 169     r__1 = inbuf[start - 1] - *dither;
 170     oldsgn = r_sign(&c_b2, &r__1);
 171     i__1 = stop;
 172     for (i__ = start; i__ <= i__1; ++i__) {
 173         lp_rms__ += (r__1 = lpbuf[i__], abs(r__1));
 174         ap_rms__ += (r__1 = inbuf[i__], abs(r__1));
 175         e_pre__ += (r__1 = inbuf[i__] - inbuf[i__ - 1], abs(r__1));
 176 /* Computing 2nd power */
 177         r__1 = inbuf[i__];
 178         e0ap += r__1 * r__1;
 179         *rc1 += inbuf[i__] * inbuf[i__ - 1];
 180 /* Computing 2nd power */
 181         r__1 = lpbuf[i__];
 182         e_0__ += r__1 * r__1;
 183 /* Computing 2nd power */
 184         r__1 = lpbuf[i__ - *mintau];
 185         e_b__ += r__1 * r__1;
 186 /* Computing 2nd power */
 187         r__1 = lpbuf[i__ + *mintau];
 188         e_f__ += r__1 * r__1;
 189         r_f__ += lpbuf[i__] * lpbuf[i__ + *mintau];
 190         r_b__ += lpbuf[i__] * lpbuf[i__ - *mintau];
 191         r__1 = inbuf[i__] + *dither;
 192         if (r_sign(&c_b2, &r__1) != oldsgn) {
 193             ++(*zc);
 194             oldsgn = -oldsgn;
 195         }
 196         *dither = -(*dither);
 197     }
 198 /*   Normalized short-term autocovariance coefficient at unit sample delay
 199  */
 200     *rc1 /= max(e0ap,1.f);
 201 /*  Ratio of the energy of the first difference signal (6 dB/oct preemphas
 202 is)*/
 203 /*   to the energy of the full band signal */
 204 /* Computing MAX */
 205     r__1 = ap_rms__ * 2.f;
 206     *qs = e_pre__ / max(r__1,1.f);
 207 /*   aR_b is the product of the forward and reverse prediction gains, */
 208 /*   looking backward in time (the causal case). */
 209     *ar_b__ = r_b__ / max(e_b__,1.f) * (r_b__ / max(e_0__,1.f));
 210 /*  aR_f is the same as aR_b, but looking forward in time (non causal case
 211 ).*/
 212     *ar_f__ = r_f__ / max(e_f__,1.f) * (r_f__ / max(e_0__,1.f));
 213 /*   Normalize ZC, LBE, and FBE to old fixed window length of 180. */
 214 /*   (The fraction 90/VLEN has a range of .58 to 1) */
 215     r__2 = (real) (*zc << 1);
 216     r__1 = r__2 * (90.f / vlen);
 217     *zc = i_nint(&r__1);
 218 /* Computing MIN */
 219     r__1 = lp_rms__ / 4 * (90.f / vlen);
 220     i__1 = i_nint(&r__1);
 221     *lbe = min(i__1,32767);
 222 /* Computing MIN */
 223     r__1 = ap_rms__ / 4 * (90.f / vlen);
 224     i__1 = i_nint(&r__1);
 225     *fbe = min(i__1,32767);
 226     return 0;
 227 } /* vparms_ */
 228