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  )

extern

Definition at line 99 of file mload.c.

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

References order.

Referenced by analys_().