e1856096ba2bc520bcd803002f49906eef4a8edc
[Faustine.git] / interpretor / preprocessor / faust-0.9.47mr3 / tools / faust2pd / examples / seqdemo / subtractive.dsp
1
2 declare name "subtractive -- saw wave filtered with resonant lowpass";
3 declare author "Albert Graef";
4 declare version "1.0";
5
6 import("music.lib");
7
8 // control variables
9
10 // master volume and pan
11 vol = hslider("vol", 0.3, 0, 10, 0.01); // %
12 pan = hslider("pan", 0.5, 0, 1, 0.01); // %
13
14 // ADSR envelop
15 attack = hslider("attack", 0.01, 0, 1, 0.001); // sec
16 decay = hslider("decay", 0.3, 0, 1, 0.001); // sec
17 sustain = hslider("sustain", 0.5, 0, 1, 0.01); // %
18 release = hslider("release", 0.2, 0, 1, 0.001); // sec
19
20 // filter parameters
21 res = hslider("resonance (dB)", 3, 0, 20, 0.1);
22 cutoff = hslider("cutoff (harmonic)", 6, 1, 20, 0.1);
23
24 // voice parameters
25 freq = nentry("freq", 440, 20, 20000, 1); // Hz
26 gain = nentry("gain", 1, 0, 10, 0.01); // %
27 gate = button("gate"); // 0/1
28
29 // generic table-driven oscillator with phase modulation
30
31 // n = the size of the table, must be a power of 2
32 // f = the wave function, must be defined on the range [0,2*PI]
33 // freq = the desired frequency in Hz
34 // mod = the phase modulation signal, in radians
35
36 tblosc(n,f,freq,mod) = (1-d)*rdtable(n,waveform,i&(n-1)) +
37 d*rdtable(n,waveform,(i+1)&(n-1))
38 with {
39 waveform = time*(2.0*PI)/n : f;
40 phase = freq/SR : (+ : decimal) ~ _;
41 modphase = decimal(phase+mod/(2*PI))*n;
42 i = int(floor(modphase));
43 d = decimal(modphase);
44 };
45
46 // resonant lowpass
47
48 // This is a tweaked Butterworth filter by David Werner and Patrice Tarrabia,
49 // see http://www.musicdsp.org and http://www.experimentalscene.com for
50 // details.
51
52 // res = resonance in dB above DC gain
53 // freq = cutoff frequency
54
55 lowpass(res,freq) = f : (+ ~ g) : *(a)
56 with {
57 f(x) = a0*x+a1*x'+a2*x'';
58 g(y) = 0-b1*y-b2*y';
59 a = 1/db2linear(0.5*res);
60
61 c = 1.0/tan(PI*(freq/SR));
62 c2 = c*c;
63 r = 1/db2linear(2.0*res);
64 q = sqrt(2.0)*r;
65 a0 = 1.0/(1.0+(q*c)+(c2));
66 a1 = 2.0*a0;
67 a2 = a0;
68 b1 = 2.0*a0*(1.0-c2);
69 b2 = a0*(1.0-q*c+c2);
70 };
71
72 // subtractive synth (saw wave passed through resonant lowpass)
73
74 saw(x) = x/PI-1;
75
76 smooth(c) = *(1-c) : +~*(c);
77
78 process = tblosc(1<<16, saw, freq, 0) : ((env,freq,_) : filter) :
79 *(env * (gain/*:smooth(0.999)*/))
80 : vgroup("3-master", *(vol) : panner(pan))
81 with {
82 env = gate : vgroup("1-adsr", adsr(attack, decay, sustain, release));
83 filter(env,freq)
84 = vgroup("2-filter", lowpass(env*res, fmax(1/cutoff, env)*freq*cutoff));
85 };