diff options
Diffstat (limited to '1.4.23-rc4/codecs/lpc10/tbdm.c')
| -rw-r--r-- | 1.4.23-rc4/codecs/lpc10/tbdm.c | 188 |
1 files changed, 0 insertions, 188 deletions
diff --git a/1.4.23-rc4/codecs/lpc10/tbdm.c b/1.4.23-rc4/codecs/lpc10/tbdm.c deleted file mode 100644 index 2f6f3d692..000000000 --- a/1.4.23-rc4/codecs/lpc10/tbdm.c +++ /dev/null @@ -1,188 +0,0 @@ -/* - -$Log$ -Revision 1.15 2004/06/26 03:50:14 markster -Merge source cleanups (bug #1911) - -Revision 1.14 2003/02/12 13:59:15 matteo -mer feb 12 14:56:57 CET 2003 - -Revision 1.1.1.1 2003/02/12 13:59:15 matteo -mer feb 12 14:56:57 CET 2003 - -Revision 1.2 2000/01/05 08:20:40 markster -Some OSS fixes and a few lpc changes to make it actually work - - * Revision 1.1 1996/08/19 22:30:26 jaf - * Initial revision - * - -*/ - -/* -- translated by f2c (version 19951025). - You must link the resulting object file with the libraries: - -lf2c -lm (in that order) -*/ - -#include "f2c.h" - -#ifdef P_R_O_T_O_T_Y_P_E_S -extern int tbdm_(real *speech, integer *lpita, integer *tau, integer *ltau, real *amdf, integer *minptr, integer *maxptr, integer *mintau); -/*:ref: difmag_ 14 8 6 4 4 4 4 6 4 4 */ -#endif - -/* ********************************************************************** */ - -/* TBDM Version 49 */ - -/* $Log$ - * Revision 1.15 2004/06/26 03:50:14 markster - * Merge source cleanups (bug #1911) - * - * Revision 1.14 2003/02/12 13:59:15 matteo - * mer feb 12 14:56:57 CET 2003 - * - * Revision 1.1.1.1 2003/02/12 13:59:15 matteo - * mer feb 12 14:56:57 CET 2003 - * - * Revision 1.2 2000/01/05 08:20:40 markster - * Some OSS fixes and a few lpc changes to make it actually work - * - * Revision 1.1 1996/08/19 22:30:26 jaf - * Initial revision - * */ -/* Revision 1.3 1996/03/18 22:14:00 jaf */ -/* Just added a few comments about which array indices of the arguments */ -/* are used, and mentioning that this subroutine has no local state. */ - -/* Revision 1.2 1996/03/13 14:48:37 jaf */ -/* Comments added explaining that none of the local variables of this */ -/* subroutine need to be saved from one invocation to the next. */ - -/* Revision 1.1 1996/02/07 14:49:54 jaf */ -/* Initial revision */ - - -/* ********************************************************************* */ - -/*TURBO DIFMAG: Compute High Resolution Average Magnitude Difference Function -*/ - -/* Note: There are several constants in here that appear to depend on a */ -/* particular TAU table. That's not a problem for the LPC10 coder, but */ -/* watch out if you change the contents of TAU in the subroutine ANALYS. */ - -/* Input: */ -/* SPEECH - Low pass filtered speech */ -/* Indices 1 through MAX+LPITA-1 are read, where: */ -/* MAX = (TAU(LTAU)-TAU(1))/2+1 */ -/* (If TAU(1) .LT. 39, then larger indices could be read */ -/* by the last call to DIFMAG below.) */ -/* LPITA - Length of speech buffer */ -/* TAU - Table of lags, sorted in increasing order. */ -/* Indices 1 through LTAU read. */ -/* LTAU - Number of lag values to compute */ -/* Output: */ -/* AMDF - Average Magnitude Difference for each lag in TAU */ -/* Indices 1 through LTAU written, and several might then be read.*/ -/* MINPTR - Index of minimum AMDF value */ -/* MAXPTR - Index of maximum AMDF value within +/- 1/2 octave of min */ -/* MINTAU - Lag corresponding to minimum AMDF value */ - -/* This subroutine has no local state. */ - -/* Subroutine */ int tbdm_(real *speech, integer *lpita, integer *tau, - integer *ltau, real *amdf, integer *minptr, integer *maxptr, integer * - mintau) -{ - /* 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, <au2, &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, <au2, &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_ */ - |