713b26a1b818501af2588344cc5b65092ee1d685
[asterisk/asterisk.git] / codecs / lpc10 / lpfilt.c
1 /*
2
3 $Log$
4 Revision 1.10  2000/01/05 08:20:39  markster
5 Version 0.1.11 from FTP
6
7 Revision 1.2  2000/01/05 08:20:39  markster
8 Some OSS fixes and a few lpc changes to make it actually work
9
10  * Revision 1.1  1996/08/19  22:31:35  jaf
11  * Initial revision
12  *
13
14 */
15
16 #ifdef P_R_O_T_O_T_Y_P_E_S
17 extern int lpfilt_(real *inbuf, real *lpbuf, integer *len, integer *nsamp);
18 #endif
19
20 /*  -- translated by f2c (version 19951025).
21    You must link the resulting object file with the libraries:
22         -lf2c -lm   (in that order)
23 */
24
25 #include "f2c.h"
26
27 /* *********************************************************************** */
28
29 /*      LPFILT Version 55 */
30
31 /* $Log$
32  * Revision 1.10  2000/01/05 08:20:39  markster
33  * Version 0.1.11 from FTP
34  *
35 /* Revision 1.2  2000/01/05 08:20:39  markster
36 /* Some OSS fixes and a few lpc changes to make it actually work
37 /*
38  * Revision 1.1  1996/08/19  22:31:35  jaf
39  * Initial revision
40  * */
41 /* Revision 1.3  1996/03/15  16:53:49  jaf */
42 /* Just put comment header in standard form. */
43
44 /* Revision 1.2  1996/03/12  23:58:06  jaf */
45 /* Comments added explaining that none of the local variables of this */
46 /* subroutine need to be saved from one invocation to the next. */
47
48 /* Revision 1.1  1996/02/07 14:47:44  jaf */
49 /* Initial revision */
50
51
52 /* *********************************************************************** */
53
54 /*   31 Point Equiripple FIR Low-Pass Filter */
55 /*     Linear phase, delay = 15 samples */
56
57 /*      Passband:  ripple = 0.25 dB, cutoff =  800 Hz */
58 /*      Stopband:  atten. =  40. dB, cutoff = 1240 Hz */
59
60 /* Inputs: */
61 /*  LEN    - Length of speech buffers */
62 /*  NSAMP  - Number of samples to filter */
63 /*  INBUF  - Input speech buffer */
64 /*           Indices len-nsamp-29 through len are read. */
65 /* Output: */
66 /*  LPBUF  - Low passed speech buffer (must be different array than INBUF) */
67 /*           Indices len+1-nsamp through len are written. */
68
69 /* This subroutine has no local state. */
70
71 /* Subroutine */ int lpfilt_(real *inbuf, real *lpbuf, integer *len, integer *
72         nsamp)
73 {
74     /* System generated locals */
75     integer i__1;
76
77     /* Local variables */
78     integer j;
79     real t;
80
81 /*      Arguments */
82 /*      Parameters/constants */
83 /*       Local variables that need not be saved */
84 /*       Local state */
85 /*       None */
86     /* Parameter adjustments */
87     --lpbuf;
88     --inbuf;
89
90     /* Function Body */
91     i__1 = *len;
92     for (j = *len + 1 - *nsamp; j <= i__1; ++j) {
93         t = (inbuf[j] + inbuf[j - 30]) * -.0097201988f;
94         t += (inbuf[j - 1] + inbuf[j - 29]) * -.0105179986f;
95         t += (inbuf[j - 2] + inbuf[j - 28]) * -.0083479648f;
96         t += (inbuf[j - 3] + inbuf[j - 27]) * 5.860774e-4f;
97         t += (inbuf[j - 4] + inbuf[j - 26]) * .0130892089f;
98         t += (inbuf[j - 5] + inbuf[j - 25]) * .0217052232f;
99         t += (inbuf[j - 6] + inbuf[j - 24]) * .0184161253f;
100         t += (inbuf[j - 7] + inbuf[j - 23]) * 3.39723e-4f;
101         t += (inbuf[j - 8] + inbuf[j - 22]) * -.0260797087f;
102         t += (inbuf[j - 9] + inbuf[j - 21]) * -.0455563702f;
103         t += (inbuf[j - 10] + inbuf[j - 20]) * -.040306855f;
104         t += (inbuf[j - 11] + inbuf[j - 19]) * 5.029835e-4f;
105         t += (inbuf[j - 12] + inbuf[j - 18]) * .0729262903f;
106         t += (inbuf[j - 13] + inbuf[j - 17]) * .1572008878f;
107         t += (inbuf[j - 14] + inbuf[j - 16]) * .2247288674f;
108         t += inbuf[j - 15] * .250535965f;
109         lpbuf[j] = t;
110     }
111     return 0;
112 } /* lpfilt_ */
113