mload.c File Reference

#include "f2c.h"

Include dependency graph for mload.c:

Go to the source code of this file.

Functions
int	mload_ (integer *order, integer *awins, integer *awinf, real *speech, real *phi, real *psi)

Function Documentation

mload_()

int mload_	(	integer *	order,
		integer *	awins,
		integer *	awinf,
		real *	speech,
		real *	phi,
		real *	psi
	)

Definition at line 99 of file mload.c.

References order.

Referenced by analys_().

{
 /* 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_ */