/* $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:39 markster Some OSS fixes and a few lpc changes to make it actually work * Revision 1.1 1996/08/19 22:31:25 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 mload_(integer *order, integer *awins, integer *awinf, real *speech, real *phi, real *psi); #endif /* ***************************************************************** */ /* MLOAD Version 48 */ /* $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:39 markster * Some OSS fixes and a few lpc changes to make it actually work * * Revision 1.1 1996/08/19 22:31:25 jaf * Initial revision * */ /* Revision 1.5 1996/03/27 23:59:51 jaf */ /* Added some more accurate comments about which indices of the argument */ /* array SPEECH are read. I thought that this might be the cause of a */ /* problem I've been having, but it isn't. */ /* Revision 1.4 1996/03/26 19:16:53 jaf */ /* Commented out the code at the end that copied the lower triangular */ /* half of PHI into the upper triangular half (making the resulting */ /* matrix symmetric). The upper triangular half was never used by later */ /* code in subroutine ANALYS. */ /* Revision 1.3 1996/03/18 21:16: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 16:47:41 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:48:01 jaf */ /* Initial revision */ /* ***************************************************************** */ /* Load a covariance matrix. */ /* Input: */ /* ORDER - Analysis order */ /* AWINS - Analysis window start */ /* AWINF - Analysis window finish */ /* SPEECH(AWINF) - Speech buffer */ /* Indices MIN(AWINS, AWINF-(ORDER-1)) through */ /* MAX(AWINF, AWINS+(ORDER-1)) read. */ /* As long as (AWINF-AWINS) .GE. (ORDER-1), */ /* this is just indices AWINS through AWINF. */ /* Output: */ /* PHI(ORDER,ORDER) - Covariance matrix */ /* Lower triangular half and diagonal written, and read.*/ /* Upper triangular half untouched. */ /* PSI(ORDER) - Prediction vector */ /* Indices 1 through ORDER written, */ /* and most are read after that. */ /* This subroutine has no local state. */ /* Subroutine */ int mload_(integer *order, integer *awins, integer *awinf, real *speech, real *phi, real *psi) { /* System generated locals */ integer phi_dim1, phi_offset, i__1, i__2; /* Local variables */ integer c__, i__, r__, start; /* Arguments */ /* Local variables that need not be saved */ /* Load first column of triangular covariance matrix PHI */ /* Parameter adjustments */ --psi; phi_dim1 = *order; phi_offset = phi_dim1 + 1; phi -= phi_offset; --speech; /* Function Body */ start = *awins + *order; i__1 = *order; for (r__ = 1; r__ <= i__1; ++r__) { phi[r__ + phi_dim1] = 0.f; i__2 = *awinf; for (i__ = start; i__ <= i__2; ++i__) { phi[r__ + phi_dim1] += speech[i__ - 1] * speech[i__ - r__]; } } /* Load last element of vector PSI */ psi[*order] = 0.f; i__1 = *awinf; for (i__ = start; i__ <= i__1; ++i__) { psi[*order] += speech[i__] * speech[i__ - *order]; } /* End correct to get additional columns of PHI */ i__1 = *order; for (r__ = 2; r__ <= i__1; ++r__) { i__2 = r__; for (c__ = 2; c__ <= i__2; ++c__) { phi[r__ + c__ * phi_dim1] = phi[r__ - 1 + (c__ - 1) * phi_dim1] - speech[*awinf + 1 - r__] * speech[*awinf + 1 - c__] + speech[start - r__] * speech[start - c__]; } } /* End correct to get additional elements of PSI */ i__1 = *order - 1; for (c__ = 1; c__ <= i__1; ++c__) { psi[c__] = phi[c__ + 1 + phi_dim1] - speech[start - 1] * speech[start - 1 - c__] + speech[*awinf] * speech[*awinf - c__]; } /* Copy lower triangular section into upper (why bother?) */ /* I'm commenting this out, since the upper triangular half of PHI */ /* is never used by later code, unless a sufficiently high level of */ /* tracing is turned on. */ /* DO R = 1,ORDER */ /* DO C = 1,R-1 */ /* PHI(C,R) = PHI(R,C) */ /* END DO */ /* END DO */ return 0; } /* mload_ */