1 /*
2
3 * Revision 1.1 1996/08/19 22:31:25 jaf
4 * Initial revision
5 *
6
7 */
8
9 /* -- translated by f2c (version 19951025).
10 You must link the resulting object file with the libraries:
11 -lf2c -lm (in that order)
12 */
13
14 #include "f2c.h"
15
16 extern int mload_(integer *order, integer *awins, integer *awinf, real *speech, real *phi, real *psi);
17
18 /* ***************************************************************** */
19
20 /* MLOAD Version 48 */
21
22 /*
23 * Revision 1.1 1996/08/19 22:31:25 jaf
24 * Initial revision
25 * */
26 /* Revision 1.5 1996/03/27 23:59:51 jaf */
27 /* Added some more accurate comments about which indices of the argument */
28 /* array SPEECH are read. I thought that this might be the cause of a */
29 /* problem I've been having, but it isn't. */
30
31 /* Revision 1.4 1996/03/26 19:16:53 jaf */
32 /* Commented out the code at the end that copied the lower triangular */
33 /* half of PHI into the upper triangular half (making the resulting */
34 /* matrix symmetric). The upper triangular half was never used by later */
35 /* code in subroutine ANALYS. */
36
37 /* Revision 1.3 1996/03/18 21:16:00 jaf */
38 /* Just added a few comments about which array indices of the arguments */
39 /* are used, and mentioning that this subroutine has no local state. */
40
41 /* Revision 1.2 1996/03/13 16:47:41 jaf */
42 /* Comments added explaining that none of the local variables of this */
43 /* subroutine need to be saved from one invocation to the next. */
44
45 /* Revision 1.1 1996/02/07 14:48:01 jaf */
46 /* Initial revision */
47
48
49 /* ***************************************************************** */
50
51 /* Load a covariance matrix. */
52
53 /* Input: */
54 /* ORDER - Analysis order */
55 /* AWINS - Analysis window start */
56 /* AWINF - Analysis window finish */
57 /* SPEECH(AWINF) - Speech buffer */
58 /* Indices MIN(AWINS, AWINF-(ORDER-1)) through */
59 /* MAX(AWINF, AWINS+(ORDER-1)) read. */
60 /* As long as (AWINF-AWINS) .GE. (ORDER-1), */
61 /* this is just indices AWINS through AWINF. */
62 /* Output: */
63 /* PHI(ORDER,ORDER) - Covariance matrix */
64 /* Lower triangular half and diagonal written, and read.*/
65 /* Upper triangular half untouched. */
66 /* PSI(ORDER) - Prediction vector */
67 /* Indices 1 through ORDER written, */
68 /* and most are read after that. */
69
70 /* This subroutine has no local state. */
71
mload_(integer * order,integer * awins,integer * awinf,real * speech,real * phi,real * psi)72 /* Subroutine */ int mload_(integer *order, integer *awins, integer *awinf,
73 real *speech, real *phi, real *psi)
74 {
75 /* System generated locals */
76 integer phi_dim1, phi_offset, i__1, i__2;
77
78 /* Local variables */
79 integer c__, i__, r__, start;
80
81 /* Arguments */
82 /* Local variables that need not be saved */
83 /* Load first column of triangular covariance matrix PHI */
84 /* Parameter adjustments */
85 --psi;
86 phi_dim1 = *order;
87 phi_offset = phi_dim1 + 1;
88 phi -= phi_offset;
89 --speech;
90
91 /* Function Body */
92 start = *awins + *order;
93 i__1 = *order;
94 for (r__ = 1; r__ <= i__1; ++r__) {
95 phi[r__ + phi_dim1] = 0.f;
96 i__2 = *awinf;
97 for (i__ = start; i__ <= i__2; ++i__) {
98 phi[r__ + phi_dim1] += speech[i__ - 1] * speech[i__ - r__];
99 }
100 }
101 /* Load last element of vector PSI */
102 psi[*order] = 0.f;
103 i__1 = *awinf;
104 for (i__ = start; i__ <= i__1; ++i__) {
105 psi[*order] += speech[i__] * speech[i__ - *order];
106 }
107 /* End correct to get additional columns of PHI */
108 i__1 = *order;
109 for (r__ = 2; r__ <= i__1; ++r__) {
110 i__2 = r__;
111 for (c__ = 2; c__ <= i__2; ++c__) {
112 phi[r__ + c__ * phi_dim1] = phi[r__ - 1 + (c__ - 1) * phi_dim1] -
113 speech[*awinf + 1 - r__] * speech[*awinf + 1 - c__] +
114 speech[start - r__] * speech[start - c__];
115 }
116 }
117 /* End correct to get additional elements of PSI */
118 i__1 = *order - 1;
119 for (c__ = 1; c__ <= i__1; ++c__) {
120 psi[c__] = phi[c__ + 1 + phi_dim1] - speech[start - 1] * speech[start
121 - 1 - c__] + speech[*awinf] * speech[*awinf - c__];
122 }
123 /* Copy lower triangular section into upper (why bother?) */
124 /* I'm commenting this out, since the upper triangular half of PHI
125 */
126 /* is never used by later code, unless a sufficiently high level of
127 */
128 /* tracing is turned on. */
129 /* DO R = 1,ORDER */
130 /* DO C = 1,R-1 */
131 /* PHI(C,R) = PHI(R,C) */
132 /* END DO */
133 /* END DO */
134 return 0;
135 } /* mload_ */
136
137