blob: ca03055a0d6699cb8e4710305819314bf407264b (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
|
(ALLPOLES-ALG
(NAME "allpoles")
(ARGUMENTS ("sound_type" "x_snd")("LVAL" "ak_array")("double" "gain"))
(START (MIN x_snd))
(NOT-IN-INNER-LOOP "ak_array")
(ALWAYS-SCALE x_snd)
(TERMINATE (MIN x_snd))
(LOGICAL-STOP (MIN x_snd))
(STATE
("long" "ak_len" "0") ; length of coefs ak array
("LVAL" "ak_array" "ak_array")
("double" "gain" "gain")
("double *" "ak_coefs" "NULL") ; coefs array
("double *" "zk_buf" "NULL") ; last values of output
("long" "index" "0")
)
(OUTER-LOOP "
if (susp->ak_array == NULL) {
out: togo = 0; /* indicate termination */
break; /* we're done */
}
else if (!vectorp(susp->ak_array))
xlerror(\"array expected\", susp->ak_array);
else if (susp->ak_coefs == NULL)
{
long i;
susp->ak_len = getsize(susp->ak_array);
if (susp->ak_len < 1) xlerror(\"array has not elements\", susp->ak_array);
susp->ak_coefs = (double *) calloc(susp->ak_len, sizeof(double));
susp->zk_buf = (double *) calloc(susp->ak_len, sizeof(double));
/* at this point we have a new array and a place to put ak coefs */
for(i=0; i < susp->ak_len; i++) {
LVAL elem = getelement(susp->ak_array,i);
if (ntype(elem) != FLONUM) {
xlerror(\"flonum expected\", elem);
}
susp->ak_coefs[i] = getflonum(elem);
}
}
")
(CONSTANT "ak_array" "ak_coefs" "ak_len" "gain")
(SAMPLE-RATE "x_snd->sr")
(INNER-LOOP-LOCALS "double z0; long xi; long xj;")
(INNER-LOOP "
z0 = x_snd*gain;
for (xi=0; xi < ak_len ; xi++)
{
xj = index + xi; if (xj >= ak_len) xj -= ak_len;
z0 += ak_coefs[xi] * zk_buf[xj];
}
zk_buf[index] = z0;
index++; if (index == ak_len) index = 0;
output = (sample_type) z0;
")
(FINALIZATION "
free(susp->zk_buf);
free(susp->ak_coefs);
susp->ak_array = NULL; /* free array */
")
)
|
