diff options
Diffstat (limited to '1.2-netsec/codecs/lpc10/mload.c')
-rw-r--r-- | 1.2-netsec/codecs/lpc10/mload.c | 163 |
1 files changed, 0 insertions, 163 deletions
diff --git a/1.2-netsec/codecs/lpc10/mload.c b/1.2-netsec/codecs/lpc10/mload.c deleted file mode 100644 index 749646097..000000000 --- a/1.2-netsec/codecs/lpc10/mload.c +++ /dev/null @@ -1,163 +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: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_ */ - |