tbdm.c File Reference

#include "f2c.h"

Include dependency graph for tbdm.c:

Go to the source code of this file.

Functions
int	tbdm_ (real *speech, integer *lpita, integer *tau, integer *ltau, real *amdf, integer *minptr, integer *maxptr, integer *mintau)

Function Documentation

tbdm_()

int tbdm_	(	real *	speech,
		integer *	lpita,
		integer *	tau,
		integer *	ltau,
		real *	amdf,
		integer *	minptr,
		integer *	maxptr,
		integer *	mintau
	)

Definition at line 94 of file tbdm.c.

References difmag_(), max, and min.

Referenced by analys_().

{
    /* System generated locals */
    integer i__1, i__2, i__3, i__4;

    /* Local variables */
    real amdf2[6];
    integer minp2, ltau2, maxp2, i__;
    extern /* Subroutine */ int difmag_(real *, integer *, integer *, integer
       *, integer *, real *, integer *, integer *);
    integer minamd, ptr, tau2[6];

/*    Arguments */
/*    REAL SPEECH(LPITA+TAU(LTAU)), AMDF(LTAU) */
/*   Stupid TOAST doesn't understand expressions */
/*       Local variables that need not be saved */
/*       Local state */
/*       None */
/*   Compute full AMDF using log spaced lags, find coarse minimum */
    /* Parameter adjustments */
    --speech;
    --amdf;
    --tau;

    /* Function Body */
    difmag_(&speech[1], lpita, &tau[1], ltau, &tau[*ltau], &amdf[1], minptr,
       maxptr);
    *mintau = tau[*minptr];
    minamd = (integer)amdf[*minptr];
/*   Build table containing all lags within +/- 3 of the AMDF minimum */
/*    excluding all that have already been computed */
    ltau2 = 0;
    ptr = *minptr - 2;
/* Computing MAX */
    i__1 = *mintau - 3;
/* Computing MIN */
    i__3 = *mintau + 3, i__4 = tau[*ltau] - 1;
    i__2 = min(i__3,i__4);
    for (i__ = max(i__1,41); i__ <= i__2; ++i__) {
   while(tau[ptr] < i__) {
       ++ptr;
   }
   if (tau[ptr] != i__) {
       ++ltau2;
       tau2[ltau2 - 1] = i__;
   }
    }
/*   Compute AMDF of the new lags, if there are any, and choose one */
/*    if it is better than the coarse minimum */
    if (ltau2 > 0) {
   difmag_(&speech[1], lpita, tau2, &ltau2, &tau[*ltau], amdf2, &minp2, &
      maxp2);
   if (amdf2[minp2 - 1] < (real) minamd) {
       *mintau = tau2[minp2 - 1];
       minamd = (integer)amdf2[minp2 - 1];
   }
    }
/*   Check one octave up, if there are any lags not yet computed */
    if (*mintau >= 80) {
   i__ = *mintau / 2;
   if ((i__ & 1) == 0) {
       ltau2 = 2;
       tau2[0] = i__ - 1;
       tau2[1] = i__ + 1;
   } else {
       ltau2 = 1;
       tau2[0] = i__;
   }
   difmag_(&speech[1], lpita, tau2, &ltau2, &tau[*ltau], amdf2, &minp2, &
      maxp2);
   if (amdf2[minp2 - 1] < (real) minamd) {
       *mintau = tau2[minp2 - 1];
       minamd = (integer)amdf2[minp2 - 1];
       *minptr += -20;
   }
    }
/*   Force minimum of the AMDF array to the high resolution minimum */
    amdf[*minptr] = (real) minamd;
/*   Find maximum of AMDF within 1/2 octave of minimum */
/* Computing MAX */
    i__2 = *minptr - 5;
    *maxptr = max(i__2,1);
/* Computing MIN */
    i__1 = *minptr + 5;
    i__2 = min(i__1,*ltau);
    for (i__ = *maxptr + 1; i__ <= i__2; ++i__) {
   if (amdf[i__] > amdf[*maxptr]) {
       *maxptr = i__;
   }
    }
    return 0;
} /* tbdm_ */