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