96ca1c7a82b660cbf578f7215106b10a98443d9a
[Faustine.git] / interpretor / faust-0.9.47mr3 / examples / faust-stk / flute.dsp
1 declare name "Flute";
2 declare description "Nonlinear WaveGuide Flute";
3 declare author "Romain Michon (rmichon@ccrma.stanford.edu)";
4 declare copyright "Romain Michon";
5 declare version "1.0";
6 declare licence "STK-4.3"; // Synthesis Tool Kit 4.3 (MIT style license);
7 declare description "A simple flute based on Smith algorythm: https://ccrma.stanford.edu/~jos/pasp/Flutes_Recorders_Pipe_Organs.html";
8
9 import("music.lib");
10 import("instrument.lib");
11 import("effect.lib");
12
13 //==================== GUI SPECIFICATION ================
14
15 freq = nentry("h:Basic_Parameters/freq [1][unit:Hz] [tooltip:Tone frequency]",440,20,20000,1);
16 gain = nentry("h:Basic_Parameters/gain [1][tooltip:Gain (value between 0 and 1)]",1,0,1,0.01);
17 gate = button("h:Basic_Parameters/gate [1][tooltip:noteOn = 1, noteOff = 0]") : int;
18
19 pressure = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Pressure
20 [2][tooltip:Breath pressure (value bewteen 0 and 1)]",0.9,0,1.5,0.01) : smooth(0.999);
21 breathAmp = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Noise Gain
22 [2][tooltip:Breath noise gain (value between 0 and 1)]",0.1,0,1,0.01)/10;
23
24 typeModulation = nentry("h:Physical_and_Nonlinearity/v:Nonlinear_Filter_Parameters/Modulation_Type
25 [3][tooltip: 0=theta is modulated by the incoming signal; 1=theta is modulated by the averaged incoming signal;
26 2=theta is modulated by the squared incoming signal; 3=theta is modulated by a sine wave of frequency freqMod;
27 4=theta is modulated by a sine wave of frequency freq;]",0,0,4,1);
28 nonLinearity = hslider("h:Physical_and_Nonlinearity/v:Nonlinear_Filter_Parameters/Nonlinearity
29 [3][tooltip:Nonlinearity factor (value between 0 and 1)]",0,0,1,0.01);
30 frequencyMod = hslider("h:Physical_and_Nonlinearity/v:Nonlinear_Filter_Parameters/Modulation_Frequency
31 [3][unit:Hz][tooltip:Frequency of the sine wave for the modulation of theta (works if Modulation Type=3)]",220,20,1000,0.1);
32 nonLinAttack = hslider("h:Physical_and_Nonlinearity/v:Nonlinear_Filter_Parameters/Nonlinearity Attack
33 [3][unit:s][Attack duration of the nonlinearity]",0.1,0,2,0.01);
34
35 vibratoFreq = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Freq
36 [4][unit:Hz]",5,1,15,0.1);
37 vibratoGain = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Gain
38 [4][tooltip:A value between 0 and 1]",0.1,0,1,0.01);
39 vibratoBegin = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Begin
40 [4][unit:s][tooltip:Vibrato silence duration before attack]",0.1,0,2,0.01);
41 vibratoAttack = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Attack
42 [4][unit:s][tooltip:Vibrato attack duration]",0.5,0,2,0.01);
43 vibratoRelease = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Release
44 [4][unit:s][tooltip:Vibrato release duration]",0.2,0,2,0.01);
45
46 pressureEnvelope = checkbox("h:Envelopes_and_Vibrato/v:Pressure_Envelope_Parameters/Pressure_Env
47 [5][unit:s][tooltip:Activate Pressure envelope]") : int;
48 env1Attack = hslider("h:Envelopes_and_Vibrato/v:Pressure_Envelope_Parameters/Press_Env_Attack
49 [5][unit:s][tooltip:Pressure envelope attack duration]",0.05,0,2,0.01);
50 env1Decay = hslider("h:Envelopes_and_Vibrato/v:Pressure_Envelope_Parameters/Press_Env_Decay
51 [5][unit:s][tooltip:Pressure envelope decay duration]",0.2,0,2,0.01);
52 env1Release = hslider("h:Envelopes_and_Vibrato/v:Pressure_Envelope_Parameters/Press_Env_Release
53 [5][unit:s][tooltip:Pressure envelope release duration]",1,0,2,0.01);
54
55 env2Attack = hslider("h:Envelopes_and_Vibrato/v:Global_Envelope_Parameters/Glob_Env_Attack
56 [6][unit:s][tooltip:Global envelope attack duration]",0.1,0,2,0.01);
57 env2Release = hslider("h:Envelopes_and_Vibrato/v:Global_Envelope_Parameters/Glob_Env_Release
58 [6][unit:s][tooltip:Global envelope release duration]",0.1,0,2,0.01);
59
60 //==================== SIGNAL PROCESSING ================
61
62 //----------------------- Nonlinear filter ----------------------------
63 //nonlinearities are created by the nonlinear passive allpass ladder filter declared in filter.lib
64
65 //nonlinear filter order
66 nlfOrder = 6;
67
68 //attack - sustain - release envelope for nonlinearity (declared in instrument.lib)
69 envelopeMod = asr(nonLinAttack,100,0.1,gate);
70
71 //nonLinearModultor is declared in instrument.lib, it adapts allpassnn from filter.lib
72 //for using it with waveguide instruments
73 NLFM = nonLinearModulator((nonLinearity : smooth(0.999)),envelopeMod,freq,
74 typeModulation,(frequencyMod : smooth(0.999)),nlfOrder);
75
76 //----------------------- Synthesis parameters computing and functions declaration ----------------------------
77
78 //Loops feedbacks gains
79 feedBack1 = 0.4;
80 feedBack2 = 0.4;
81
82 //Delay Lines
83 embouchureDelayLength = (SR/freq)/2-2;
84 boreDelayLength = SR/freq-2;
85 embouchureDelay = fdelay(4096,embouchureDelayLength);
86 boreDelay = fdelay(4096,boreDelayLength);
87
88 //Polinomial
89 poly = _ <: _ - _*_*_;
90
91 //jet filter is a lowwpass filter (declared in filter.lib)
92 reflexionFilter = lowpass(1,2000);
93
94 //stereoizer is declared in instrument.lib and implement a stereo spacialisation in function of
95 //the frequency period in number of samples
96 stereo = stereoizer(SR/freq);
97
98 //----------------------- Algorithm implementation ----------------------------
99
100 //Pressure envelope
101 env1 = adsr(env1Attack,env1Decay,90,env1Release,(gate | pressureEnvelope))*pressure*1.1;
102
103 //Global envelope
104 env2 = asr(env2Attack,100,env2Release,gate)*0.5;
105
106 //Vibrato Envelope
107 vibratoEnvelope = envVibrato(vibratoBegin,vibratoAttack,100,vibratoRelease,gate)*vibratoGain;
108
109 vibrato = osc(vibratoFreq)*vibratoEnvelope;
110
111 breath = noise*env1;
112
113 flow = env1 + breath*breathAmp + vibrato;
114
115 //instrReverb is declared in instrument.lib
116 process = (_ <: (flow + *(feedBack1) : embouchureDelay : poly) + *(feedBack2) : reflexionFilter)~(boreDelay : NLFM) : *(env2)*gain :
117 stereo : instrReverb;