blob: a56c2fd58696a8081003aff47db5e4ef13ea2db9 (
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
|
(ARESONVV-ALG
(NAME "aresonvv")
(ARGUMENTS ("sound_type" "s1") ("sound_type" "hz1") ("sound_type" "bw")
("int" "normalization"))
(INLINE-INTERPOLATION T)
(ALWAYS-SCALE hz1 bw)
(START (MAX s1 hz1 bw))
(TERMINATE (MIN s1 hz1 bw))
(LOGICAL-STOP (MIN s1))
(SAMPLE-RATE (MAX s1))
(STEP-FUNCTION hz1 bw)
(LINEAR s1)
(STATE ("double" "scale1" "s1->scale")
("double" "c3co" "0.0")
("double" "c3p1" "0.0")
("double" "c3t4" "0.0")
("double" "omc3" "0.0")
("double" "coshz" "0.0")
("double" "c2" "0.0")
("double" "c1" "0.0")
("boolean" "recompute" "false")
("int" "normalization" "normalization")
("double" "y1" "0.0")
("double" "y2" "0.0;
hz1->scale = (sample_type) (hz1->scale * (PI2 / s1->sr));
bw->scale = (sample_type) (bw->scale * (-PI2 / s1->sr));"))
(DEPENDS ("coshz" "hz1" "cos(hz1)")
("recompute" "hz1" "true")
("c3co" "bw" "exp(bw)")
("c3p1" "bw" "c3co + 1.0")
("c3t4" "bw" "c3co * 4.0")
("omc3" "bw" "1.0 - c3co")
("recompute" "bw" "true"))
(JOINT-DEPENDENCY (("hz1" "bw")
"if (recompute) {"
" recompute = false;"
" c2 = c3t4 * coshz / c3p1;"
" c1 = (normalization == 0 ? 0.0 :"
" (normalization == 1 ? 1.0 - omc3 * sqrt(1.0 - c2 * c2 / c3t4) :"
" 1.0 - sqrt(c3p1 * c3p1 - c2 * c2) * omc3 / c3p1));"
"}"))
(CONSTANT "c1" "c2" "c3co" "coshz" "c3p1" "c3t4" "omc3"
"normalization" "scale1")
(FORCE-INTO-REGISTER recompute) ;c3t4 c3p1 normalization omc3 scale1
(INNER-LOOP-LOCALS " register double y0, current;")
(INNER-LOOP "current = s1;
y0 = c1 * current + c2 * y1 - c3co * y2;
output = (sample_type) y0;
y2 = y1; y1 = y0 - current")
)
|
