cyfaust 0.1.0__cp311-cp311-macosx_15_0_arm64.whl

cyfaust/resources/libraries/examples/physicalModeling/faust-stk/bass.dsp

@@ -0,0 +1,84 @@
+declare name "Bass";
+declare description "Nonlinear WaveGuide Acoustic Bass";
+declare author "Romain Michon";
+declare copyright "Romain Michon (rmichon@ccrma.stanford.edu)";
+declare version "1.0";
+declare licence "STK-4.3"; // Synthesis Tool Kit 4.3 (MIT style license);
+
+import("instruments.lib");
+
+//==================== GUI SPECIFICATION ================
+
+freq = nentry("h:Basic_Parameters/freq [1][unit:Hz] [tooltip:Tone frequency]",120,20,20000,1);
+gain = nentry("h:Basic_Parameters/gain [1][tooltip:Gain (value between 0 and 1)]",1,0,1,0.01); 
+gate = button("h:Basic_Parameters/gate [1][tooltip:noteOn = 1, noteOff = 0]");
+
+touchLength = hslider("v:Physical_Parameters/Touch_Length
+[2][tooltip:A value between 0 and 1]",0.15,0,1,0.01)*2;
+
+typeModulation = nentry("v:Nonlinear_Filter_Parameters/Modulation_Type 
+[3][tooltip: 0=theta is modulated by the incoming signal; 1=theta is modulated by the averaged incoming signal;
+2=theta is modulated by the squared incoming signal; 3=theta is modulated by a sine wave of frequency freqMod;
+4=theta is modulated by a sine wave of frequency freq;]",0,0,4,1);
+nonLinearity = hslider("v:Nonlinear_Filter_Parameters/Nonlinearity 
+[3][tooltip:Nonlinearity factor (value between 0 and 1)]",0,0,1,0.01);
+frequencyMod = hslider("v:Nonlinear_Filter_Parameters/Modulation_Frequency 
+[3][unit:Hz][tooltip:Frequency of the sine wave for the modulation of theta (works if Modulation Type=3)]",220,20,1000,0.1);
+
+//==================== SIGNAL PROCESSING ======================
+
+//----------------------- Nonlinear filter ----------------------------
+//nonlinearities are created by the nonlinear passive allpass ladder filter declared in miscfilter.lib
+
+//nonlinear filter order
+nlfOrder = 6; 
+
+//nonLinearModultor is declared in instruments.lib, it adapts allpassnn from miscfilter.lib 
+//for using it with waveguide instruments
+NLFM =  nonLinearModulator((nonLinearity : si.smoo),1,freq,
+     typeModulation,(frequencyMod : si.smoo),nlfOrder);
+
+//----------------------- Synthesis parameters computing and functions declaration ----------------------------
+
+//delay length in number of samples
+delayLength = float(ma.SR)/freq;
+
+//stereoizer is declared in instruments.lib and implement a stereo spacialisation in function of 
+//the frequency period in number of samples 
+stereo = stereoizer(delayLength);
+
+//string excitation
+excitation = asympT60(-0.5,-0.985,0.02,gate),no.noise*asympT60(gain,0,touchLength,gate) : 
+	   onePoleSwep : excitationFilter : excitationFilter
+	   with{
+		//the exitation filter is a one pole filter (declared in instruments.lib)
+		excitationFilter = onePole(0.035,-0.965);
+	   };
+
+//the bodyfilter is a bandpass filter (declared in instruments.lib)
+bodyFilter = bandPass(108,0.997);
+
+//the reflexion filter is pole zero filter (declared in instruments.lib) whose coefficients are 
+//modulated in function of the tone being played
+reflexionFilter = poleZero(b0,b1,a1)
+	   with{
+		//filter coefficients are stored in a C++ function
+		loopFilterb0 = ffunction(float getValueBassLoopFilterb0(float), <bass.h>,"");
+		loopFilterb1 = ffunction(float getValueBassLoopFilterb1(float), <bass.h>,"");
+		loopFiltera1 = ffunction(float getValueBassLoopFiltera1(float), <bass.h>,"");
+		freqToNoteNumber = (log - log(440))/log(2)*12 + 69 + 0.5 : int;
+		freqn = freq : freqToNoteNumber;
+		b0 = loopFilterb0(freqn);
+		b1 = loopFilterb1(freqn);
+		a1 = loopFiltera1(freqn);
+	   };
+
+delayLine = asympT60(0,delayLength,0.01,gate),_ : de.fdelay(4096);
+
+//the resonance duration is different whether a note-on signal is sent or not 
+resonanceGain = gate + (gate < 1 <: *(asympT60(1,0.9,0.05)));
+
+process = excitation : 
+	(+)~(delayLine : NLFM : reflexionFilter*resonanceGain) <: 
+	bodyFilter*1.5 + *(0.5) : *(4) : stereo : instrReverb;
+

cyfaust/resources/libraries/examples/physicalModeling/faust-stk/bass.h

@@ -0,0 +1,91 @@
+//Parameters for bass.dsp
+//©Romain Michon (rmichon@ccrma.stanford.edu), 2011
+//licence: FAUST-STK
+
+#include "instrument.h"
+
+#define _LOOKUP_TABLE_H_
+
+double bassLoopFilterb0_points[19*2] = {
+	24.000,0.54355,
+	26.000,0.54355,
+	27.000,0.55677,
+	29.000,0.55677,
+	32.000,0.55677,
+	33.000,0.83598,
+	36.000,0.83598,
+	43.000,0.83598,
+	44.000,0.88292,
+	48.000,0.88292,
+	51.000,0.88292,
+	52.000,0.77805,
+	54.000,0.77805,
+	57.000,0.77805,
+	58.000,0.91820,
+	60.000,0.91820,
+	61.000,0.91820,
+	63.000,0.94594,
+	65.000,0.91820,
+};
+extern LookupTable bassLoopFilterb0;
+LookupTable bassLoopFilterb0(&bassLoopFilterb0_points[0], 18);
+
+float getValueBassLoopFilterb0(float index){
+	return bassLoopFilterb0.getValue(index);
+}
+
+double bassLoopFilterb1_points[19*2] = {
+	24.000,-0.36586,
+	26.000,-0.36586,
+	27.000,-0.37628,
+	29.000,-0.37628,
+	32.000,-0.37628,
+	33.000,-0.60228,
+	36.000,-0.60228,
+	43.000,-0.60228,
+	44.000,-0.65721,
+	48.000,-0.65721,
+	51.000,-0.65721,
+	52.000,-0.51902,
+	54.000,-0.51902,
+	57.000,-0.51902,
+	58.000,-0.80765,
+	60.000,-0.80765,
+	61.000,-0.80765,
+	63.000,-0.83230,
+	65.000,-0.83230,
+};
+extern LookupTable bassLoopFilterb1;
+LookupTable bassLoopFilterb1(&bassLoopFilterb1_points[0], 18);
+
+float getValueBassLoopFilterb1(float index){
+	return bassLoopFilterb1.getValue(index);
+}
+
+double bassLoopFiltera1_points[19*2] = {
+	24.000,-0.81486,
+	26.000,-0.81486,
+	27.000,-0.81147,
+	29.000,-0.81147,
+	32.000,-0.81147,
+	33.000,-0.76078,
+	36.000,-0.76078,
+	43.000,-0.76078,
+	44.000,-0.77075,
+	48.000,-0.77075,
+	51.000,-0.77075,
+	52.000,-0.73548,
+	54.000,-0.73548,
+	57.000,-0.73548,
+	58.000,-0.88810,
+	60.000,-0.88810,
+	61.000,-0.88810,
+	63.000,-0.88537,
+	65.000,-0.88537,
+};
+extern LookupTable bassLoopFiltera1;
+LookupTable bassLoopFiltera1(&bassLoopFiltera1_points[0], 18);
+
+float getValueBassLoopFiltera1(float index){
+	return bassLoopFiltera1.getValue(index);
+}

cyfaust/resources/libraries/examples/physicalModeling/faust-stk/blowBottle.dsp

@@ -0,0 +1,102 @@
+declare name "blowBottle";
+declare description "Blown Bottle Instrument";
+declare author "Romain Michon (rmichon@ccrma.stanford.edu)";
+declare copyright "Romain Michon";
+declare version "1.0";
+declare licence "STK-4.3"; // Synthesis Tool Kit 4.3 (MIT style license);
+declare description "This object implements a helmholtz resonator (biquad filter) with a polynomial jet excitation (a la Cook).";
+
+import("instruments.lib");
+
+//==================== GUI SPECIFICATION ================
+
+freq = nentry("h:Basic_Parameters/freq [1][unit:Hz] [tooltip:Tone frequency]",440,20,20000,1);
+gain = nentry("h:Basic_Parameters/gain [1][tooltip:Gain (value between 0 and 1)]",1,0,1,0.01); 
+gate = button("h:Basic_Parameters/gate [1][tooltip:noteOn = 1, noteOff = 0]");
+
+noiseGain = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Noise_Gain 
+[2][tooltip:Breath noise gain (value between 0 and 1)]",0.5,0,1,0.01)*2;
+pressure = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Pressure 
+[2][tooltip:Breath pressure (value bewteen 0 and 1)]",1,0,1,0.01);
+
+typeModulation = nentry("h:Physical_and_Nonlinearity/v:Nonlinear_Filter_Parameters/Modulation_Type 
+[3][tooltip: 0=theta is modulated by the incoming signal; 1=theta is modulated by the averaged incoming signal;
+2=theta is modulated by the squared incoming signal; 3=theta is modulated by a sine wave of frequency freqMod;
+4=theta is modulated by a sine wave of frequency freq;]",0,0,4,1);
+nonLinearity = hslider("h:Physical_and_Nonlinearity/v:Nonlinear_Filter_Parameters/Nonlinearity 
+[3][tooltip:Nonlinearity factor (value between 0 and 1)]",0,0,1,0.01);
+frequencyMod = hslider("h:Physical_and_Nonlinearity/v:Nonlinear_Filter_Parameters/Modulation_Frequency 
+[3][unit:Hz][tooltip:Frequency of the sine wave for the modulation of theta (works if Modulation Type=3)]",220,20,1000,0.1);
+nonLinAttack = hslider("h:Physical_and_Nonlinearity/v:Nonlinear_Filter_Parameters/Nonlinearity_Attack
+[3][unit:s][Attack duration of the nonlinearity]",0.1,0,2,0.01);
+
+vibratoFreq = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Freq 
+[4][unit:Hz]",5,1,15,0.1);
+vibratoGain = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Gain
+[4][tooltip:A value between 0 and 1]",0.1,0,1,0.01);
+vibratoBegin = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Begin
+[4][unit:s][tooltip:Vibrato silence duration before attack]",0.05,0,2,0.01);
+vibratoAttack = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Attack 
+[4][unit:s][tooltip:Vibrato attack duration]",0.5,0,2,0.01);
+vibratoRelease = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Release 
+[4][unit:s][tooltip:Vibrato release duration]",0.01,0,2,0.01);
+
+envelopeAttack = hslider("h:Envelopes_and_Vibrato/v:Envelope_Parameters/Envelope_Attack 
+[5][unit:s][tooltip:Envelope attack duration]",0.01,0,2,0.01);
+envelopeDecay = hslider("h:Envelopes_and_Vibrato/v:Envelope_Parameters/Envelope_Decay 
+[5][unit:s][tooltip:Envelope decay duration]",0.01,0,2,0.01);
+envelopeRelease = hslider("h:Envelopes_and_Vibrato/v:Envelope_Parameters/Envelope_Release 
+[5][unit:s][tooltip:Envelope release duration]",0.5,0,2,0.01);
+
+
+//==================== SIGNAL PROCESSING ================
+
+//----------------------- Nonlinear filter ----------------------------
+//nonlinearities are created by the nonlinear passive allpass ladder filter declared in miscfilter.lib
+
+//nonlinear filter order
+nlfOrder = 6; 
+
+//attack - sustain - release envelope for nonlinearity (declared in instruments.lib)
+envelopeMod = en.asr(nonLinAttack,1,envelopeRelease,gate);
+
+//nonLinearModultor is declared in instruments.lib, it adapts allpassnn from miscfilter.lib 
+//for using it with waveguide instruments
+NLFM =  nonLinearModulator((nonLinearity : si.smoo),envelopeMod,freq,
+typeModulation,(frequencyMod : si.smoo),nlfOrder);
+
+//----------------------- Synthesis parameters computing and functions declaration ----------------------------
+
+//botlle radius
+bottleRadius = 0.999;
+
+//stereoizer is declared in instruments.lib and implement a stereo spacialisation in function of 
+//the frequency period in number of samples 
+stereo = stereoizer(ma.SR/freq);
+
+bandPassFilter = bandPass(freq,bottleRadius);
+
+//----------------------- Algorithm implementation ----------------------------
+
+//global envelope is of type attack - decay - sustain - release
+envelopeG =  gain*en.adsr(gain*envelopeAttack,envelopeDecay,1,envelopeRelease,gate);
+
+//pressure envelope is also ADSR
+envelope = pressure*en.adsr(gain*0.02,0.01,1,gain*0.2,gate);
+
+//vibrato
+vibrato = os.osc(vibratoFreq)*vibratoGain*envVibrato(vibratoBegin,vibratoAttack,100,vibratoRelease,gate)*os.osc(vibratoFreq);
+
+//breat pressure
+breathPressure = envelope + vibrato;
+
+//breath noise
+randPressure = noiseGain*no.noise*breathPressure ;
+
+process = 
+	//differential pressure
+	(-(breathPressure) <: 
+	((+(1))*randPressure : +(breathPressure)) - *(jetTable),_ : bandPassFilter,_)~NLFM : !,_ : 
+	//signal scaling
+	fi.dcblocker*envelopeG*0.5 : stereo : instrReverb;
+

cyfaust/resources/libraries/examples/physicalModeling/faust-stk/bowed.dsp

@@ -0,0 +1,114 @@
+declare name "bowed";
+declare description "Nonlinear WaveGuide Bowed Instrument";
+declare author "Romain Michon";
+declare copyright "Romain Michon (rmichon@ccrma.stanford.edu)";
+declare version "1.0";
+declare licence "STK-4.3"; // Synthesis Tool Kit 4.3 (MIT style license);
+declare description "A bowed string model, a la Smith (1986), after McIntyre, Schumacher, Woodhouse (1983).";
+declare reference "https://ccrma.stanford.edu/~jos/pasp/Bowed_Strings.html";
+
+import("instruments.lib");
+
+//==================== GUI SPECIFICATION ================
+
+freq = nentry("h:Basic_Parameters/freq [1][unit:Hz] [tooltip:Tone frequency]",440,20,20000,1);
+gain = nentry("h:Basic_Parameters/gain [1][tooltip:Gain (value between 0 and 1)]",1,0,1,0.01); 
+gate = button("h:Basic_Parameters/gate [1][tooltip:noteOn = 1, noteOff = 0]");
+
+bowPosition = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Bow_Position
+[2][tooltip:Bow position along the string (value between 0 and 1)]",0.7,0.01,1,0.01);
+bowPressure = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Bow_Pressure
+[2][tooltip:Bow pressure on the string (value between 0 and 1)]",0.75,0,1,0.01);
+
+typeModulation = nentry("h:Physical_and_Nonlinearity/v:Nonlinear_Filter_Parameters/Modulation_Type 
+[3][tooltip: 0=theta is modulated by the incoming signal; 1=theta is modulated by the averaged incoming signal;
+2=theta is modulated by the squared incoming signal; 3=theta is modulated by a sine wave of frequency freqMod;
+4=theta is modulated by a sine wave of frequency freq;]",0,0,4,1);
+nonLinearity = hslider("h:Physical_and_Nonlinearity/v:Nonlinear_Filter_Parameters/Nonlinearity 
+[3][tooltip:Nonlinearity factor (value between 0 and 1)]",0,0,1,0.01);
+frequencyMod = hslider("h:Physical_and_Nonlinearity/v:Nonlinear_Filter_Parameters/Modulation_Frequency 
+[3][unit:Hz][tooltip:Frequency of the sine wave for the modulation of theta (works if Modulation Type=3)]",220,20,1000,0.1);
+nonLinAttack = hslider("h:Physical_and_Nonlinearity/v:Nonlinear_Filter_Parameters/Nonlinearity_Attack
+[3][unit:s][Attack duration of the nonlinearity]",0.1,0,2,0.01);
+
+vibratoFreq = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Freq 
+[4][unit:Hz]",6,1,15,0.1);
+vibratoGain = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Gain
+[4][tooltip:A value between 0 and 1]",0.01,0,1,0.01);
+vibratoBegin = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Begin
+[4][unit:s][tooltip:Vibrato silence duration before attack]",0.05,0,2,0.01);
+vibratoAttack = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Attack 
+[4][unit:s][tooltip:Vibrato attack duration]",0.5,0,2,0.01);
+vibratoRelease = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Release 
+[4][unit:s][tooltip:Vibrato release duration]",0.01,0,2,0.01);
+
+envelopeAttack = hslider("h:Envelopes_and_Vibrato/v:Envelope_Parameters/Envelope_Attack 
+[5][unit:s][tooltip:Envelope attack duration]",0.01,0,2,0.01);
+envelopeDecay = hslider("h:Envelopes_and_Vibrato/v:Envelope_Parameters/Envelope_Decay 
+[5][unit:s][tooltip:Envelope decay duration]",0.05,0,2,0.01);
+envelopeRelease = hslider("h:Envelopes_and_Vibrato/v:Envelope_Parameters/Envelope_Release 
+[5][unit:s][tooltip:Envelope release duration]",0.1,0,2,0.01);
+
+//==================== SIGNAL PROCESSING ================
+
+//----------------------- Nonlinear filter ----------------------------
+//nonlinearities are created by the nonlinear passive allpass ladder filter declared in miscfilter.lib
+
+//nonlinear filter order
+nlfOrder = 6; 
+
+//attack - sustain - release envelope for nonlinearity (declared in instruments.lib)
+envelopeMod = en.asr(nonLinAttack,1,envelopeRelease,gate);
+
+//nonLinearModultor is declared in instruments.lib, it adapts allpassnn from miscfilter.lib 
+//for using it with waveguide instruments
+NLFM = nonLinearModulator((nonLinearity : si.smoo),envelopeMod,freq,
+     typeModulation,(frequencyMod : si.smoo),nlfOrder);
+
+//----------------------- Synthesis parameters computing and functions declaration ----------------------------
+
+//Parameters for the bow table
+tableOffset =  0;
+tableSlope = 5 - (4*bowPressure);
+
+//the bow table is declared in instruments.lib
+bowTable = bow(tableOffset,tableSlope);
+
+//a attack - decay - sustain - release envelope is used 
+envelope = en.adsr(gain*envelopeAttack,envelopeDecay,1, (1-gain)*envelopeRelease,gate);
+maxVelocity = 0.03 + (0.2 * gain);
+
+//Delay lines declaration and vibrato, the length of the two delay lines are evolving propotionally
+betaRatio = 0.027236 + (0.2*bowPosition);
+fdelneck = (ma.SR/freq-4)*(1 - betaRatio);
+vibratoEnvelope = envVibrato(vibratoBegin,vibratoAttack,100,vibratoRelease,gate);
+vibrato = fdelneck + ((ma.SR/freq - 4)*vibratoGain*vibratoEnvelope*os.osc(vibratoFreq));
+neckDelay = de.fdelay(4096,vibrato);
+fdelbridge = (ma.SR/freq - 4)*betaRatio;
+bridgeDelay = de.delay(4096,fdelbridge);
+
+//Body Filter: a biquad filter with a normalized pick gain (declared in instruments.lib)
+bodyFilter = bandPass(500,0.85);
+
+//String Filter: a lowpass filter (declared in instruments.lib)
+stringFilter = *(0.95) : -onePole(b0,a1)
+	with{
+		pole = 0.6 - (0.1*22050/ma.SR);
+		gain = 0.95;
+		b0 = 1-pole;
+		a1 = -pole;	
+	};
+
+//stereoizer is declared in instruments.lib and implement a stereo spacialisation in function of 
+//the frequency period in number of samples 
+stereo = stereoizer(ma.SR/freq);
+
+//----------------------- Algorithm implementation ----------------------------
+
+bowVelocity = envelope*maxVelocity;
+instrumentBody(feedBckBridge) = (*(-1) <: +(feedBckBridge),_ : (bowVelocity-_ <: *(bowTable) <: _,_),_ : 
+	_, + : +(feedBckBridge),_) ~ (neckDelay) : !,_;
+
+process = (stringFilter : instrumentBody) ~ (bridgeDelay : NLFM) : bodyFilter(*(0.2)) : 
+	_*gain*8 : stereo : instrReverb;
+

cyfaust/resources/libraries/examples/physicalModeling/faust-stk/brass.dsp

@@ -0,0 +1,103 @@
+declare name "brass";
+declare description "WaveGuide Brass instrument from STK";
+declare author "Romain Michon (rmichon@ccrma.stanford.edu)";
+declare copyright "Romain Michon";
+declare version "1.0";
+declare licence "STK-4.3"; // Synthesis Tool Kit 4.3 (MIT style license);
+declare description "A simple brass instrument waveguide model, a la Cook (TBone, HosePlayer).";
+declare reference "https://ccrma.stanford.edu/~jos/pasp/Brasses.html"; 
+
+import("instruments.lib");
+
+//==================== GUI SPECIFICATION ================
+
+freq = nentry("h:Basic_Parameters/freq [1][unit:Hz] [tooltip:Tone frequency]",440,20,20000,1);
+gain = nentry("h:Basic_Parameters/gain [1][tooltip:Gain (value between 0 and 1)]",1,0,1,0.01); 
+gate = button("h:Basic_Parameters/gate [1][tooltip:noteOn = 1, noteOff = 0]");
+
+pressure = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Pressure
+[2][tooltip:A value between 0 and 1]",1,0.01,1,0.01);
+lipTension = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Lip_Tension
+[2][tooltip:A value between 0 and 1]",0.780,0.01,1,0.001);
+slideLength = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Slide_Length
+[2][tooltip:A value between 0 and 1]",0.041,0.01,1,0.001);
+
+typeModulation = nentry("h:Physical_and_Nonlinearity/v:Nonlinear_Filter_Parameters/Modulation_Type 
+[3][tooltip: 0=theta is modulated by the incoming signal; 1=theta is modulated by the averaged incoming signal;
+2=theta is modulated by the squared incoming signal; 3=theta is modulated by a sine wave of frequency freqMod;
+4=theta is modulated by a sine wave of frequency freq;]",0,0,4,1);
+nonLinearity = hslider("h:Physical_and_Nonlinearity/v:Nonlinear_Filter_Parameters/Nonlinearity 
+[3][tooltip:Nonlinearity factor (value between 0 and 1)]",0,0,1,0.01);
+frequencyMod = hslider("h:Physical_and_Nonlinearity/v:Nonlinear_Filter_Parameters/Modulation_Frequency 
+[3][unit:Hz][tooltip:Frequency of the sine wave for the modulation of theta (works if Modulation Type=3)]",220,20,1000,0.1);
+nonLinAttack = hslider("h:Physical_and_Nonlinearity/v:Nonlinear_Filter_Parameters/Nonlinearity_Attack
+[3][unit:s][Attack duration of the nonlinearity]",0.1,0,2,0.01);
+
+vibratoFreq = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Freq 
+[4][unit:Hz]",6,1,15,0.1);
+vibratoGain = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Gain
+[4][tooltip:A value between 0 and 1]",0.05,0,1,0.01);
+vibratoBegin = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Begin
+[4][unit:s][tooltip:Vibrato silence duration before attack]",0.05,0,2,0.01);
+vibratoAttack = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Attack 
+[4][unit:s][tooltip:Vibrato attack duration]",0.5,0,2,0.01);
+vibratoRelease = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Release 
+[4][unit:s][tooltip:Vibrato release duration]",0.1,0,2,0.01);
+
+envelopeAttack = hslider("h:Envelopes_and_Vibrato/v:Envelope_Parameters/Envelope_Attack 
+[5][unit:s][tooltip:Envelope attack duration]",0.005,0,2,0.01);
+envelopeDecay = hslider("h:Envelopes_and_Vibrato/v:Envelope_Parameters/Envelope_Decay 
+[5][unit:s][tooltip:Envelope decay duration]",0.001,0,2,0.01);
+envelopeRelease = hslider("h:Envelopes_and_Vibrato/v:Envelope_Parameters/Envelope_Release 
+[5][unit:s][tooltip:Envelope release duration]",0.07,0,2,0.01);
+
+//==================== SIGNAL PROCESSING ================
+
+//----------------------- Nonlinear filter ----------------------------
+//nonlinearities are created by the nonlinear passive allpass ladder filter declared in miscfilter.lib
+
+//nonlinear filter order
+nlfOrder = 6; 
+
+//attack - sustain - release envelope for nonlinearity (declared in instruments.lib)
+envelopeMod = en.asr(nonLinAttack,1,envelopeRelease,gate);
+
+//nonLinearModultor is declared in instruments.lib, it adapts allpassnn from miscfilter.lib 
+//for using it with waveguide instruments
+NLFM =  nonLinearModulator((nonLinearity : si.smoo),envelopeMod,freq,
+     typeModulation,(frequencyMod : si.smoo),nlfOrder);
+
+//----------------------- Synthesis parameters computing and functions declaration ----------------------------
+
+//lips are simulated by a biquad filter whose output is squared and hard-clipped, bandPassH and saturationPos are declared in instruments.lib
+lipFilterFrequency = freq*pow(4,(2*lipTension)-1);
+lipFilter = *(0.03) : bandPassH(lipFilterFrequency,0.997) <: * : saturationPos;
+
+//stereoizer is declared in instruments.lib and implement a stereo spacialisation in function of 
+//the frequency period in number of samples 
+stereo = stereoizer(ma.SR/freq);
+
+//delay times in number of samples
+slideTarget = ((ma.SR/freq)*2 + 3)*(0.5 + slideLength);
+boreDelay = de.fdelay(4096,slideTarget);
+
+//----------------------- Algorithm implementation ----------------------------
+
+//vibrato
+vibrato = vibratoGain*os.osc(vibratoFreq)*envVibrato(vibratoBegin,vibratoAttack,100,vibratoRelease,gate);
+
+//envelope (Attack / Decay / Sustain / Release), breath pressure and vibrato
+breathPressure = pressure*en.adsr(envelopeAttack,envelopeDecay,1,envelopeRelease,gate) + vibrato;
+mouthPressure = 0.3*breathPressure;
+
+//scale the delay feedback
+borePressure = *(0.85);
+
+//differencial presure
+deltaPressure = mouthPressure - _;
+
+process = (borePressure <: deltaPressure,_ : 
+	  (lipFilter <: *(mouthPressure),(1-_)),_ : _, * :> + :
+	  fi.dcblocker) ~ (boreDelay : NLFM) :
+	  *(gain)*4 : stereo : instrReverb;
+

cyfaust/resources/libraries/examples/physicalModeling/faust-stk/clarinet.dsp

@@ -0,0 +1,110 @@
+declare name "clarinet";
+declare description "Nonlinear WaveGuide Clarinet";
+declare author "Romain Michon";
+declare copyright "Romain Michon (rmichon@ccrma.stanford.edu)";
+declare version "1.0";
+declare licence "STK-4.3"; // Synthesis Tool Kit 4.3 (MIT style license);
+declare description "A simple clarinet physical model, as discussed by Smith (1986), McIntyre, Schumacher, Woodhouse (1983), and others.";
+declare reference "https://ccrma.stanford.edu/~jos/pasp/Woodwinds.html";
+
+import("instruments.lib");
+
+//==================== GUI SPECIFICATION ================
+
+freq = nentry("h:Basic_Parameters/freq [1][unit:Hz] [tooltip:Tone frequency]",440,20,20000,1);
+gain = nentry("h:Basic_Parameters/gain [1][tooltip:Gain (value between 0 and 1)]",1,0,1,0.01); 
+gate = button("h:Basic_Parameters/gate [1][tooltip:noteOn = 1, noteOff = 0]");
+
+reedStiffness = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Reed_Stiffness 
+[2][tooltip:Reed stiffness (value between 0 and 1)]",0.5,0,1,0.01);
+noiseGain = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Noise_Gain 
+[2][tooltip:Breath noise gain (value between 0 and 1)]",0,0,1,0.01);
+pressure = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Pressure 
+[2][tooltip:Breath pressure (value bewteen 0 and 1)]",1,0,1,0.01);
+
+typeModulation = nentry("h:Physical_and_Nonlinearity/v:Nonlinear_Filter_Parameters/Modulation_Type 
+[3][tooltip: 0=theta is modulated by the incoming signal; 1=theta is modulated by the averaged incoming signal;
+2=theta is modulated by the squared incoming signal; 3=theta is modulated by a sine wave of frequency freqMod;
+4=theta is modulated by a sine wave of frequency freq;]",0,0,4,1);
+nonLinearity = hslider("h:Physical_and_Nonlinearity/v:Nonlinear_Filter_Parameters/Nonlinearity 
+[3][tooltip:Nonlinearity factor (value between 0 and 1)]",0,0,1,0.01);
+frequencyMod = hslider("h:Physical_and_Nonlinearity/v:Nonlinear_Filter_Parameters/Modulation_Frequency 
+[3][unit:Hz][tooltip:Frequency of the sine wave for the modulation of theta (works if Modulation Type=3)]",220,20,1000,0.1);
+nonLinAttack = hslider("h:Physical_and_Nonlinearity/v:Nonlinear_Filter_Parameters/Nonlinearity_Attack
+[3][unit:s][Attack duration of the nonlinearity]",0.1,0,2,0.01);
+
+vibratoFreq = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Freq 
+[4][unit:Hz]",5,1,15,0.1);
+vibratoGain = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Gain
+[4][tooltip:A value between 0 and 1]",0.1,0,1,0.01);
+vibratoAttack = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Attack 
+[4][unit:s][tooltip:Vibrato attack duration]",0.5,0,2,0.01);
+vibratoRelease = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Release 
+[4][unit:s][tooltip:Vibrato release duration]",0.01,0,2,0.01);
+
+envelopeAttack = hslider("h:Envelopes_and_Vibrato/v:Envelope_Parameters/Envelope_Attack 
+[5][unit:s][tooltip:Envelope attack duration]",0.01,0,2,0.01);
+envelopeDecay = hslider("h:Envelopes_and_Vibrato/v:Envelope_Parameters/Envelope_Decay 
+[5][unit:s][tooltip:Envelope decay duration]",0.05,0,2,0.01);
+envelopeRelease = hslider("h:Envelopes_and_Vibrato/v:Envelope_Parameters/Envelope_Release 
+[5][unit:s][tooltip:Envelope release duration]",0.1,0,2,0.01);
+
+//==================== SIGNAL PROCESSING ======================
+
+//----------------------- Nonlinear filter ----------------------------
+//nonlinearities are created by the nonlinear passive allpass ladder filter declared in miscfilter.lib
+
+//nonlinear filter order
+nlfOrder = 6; 
+
+//attack - sustain - release envelope for nonlinearity (declared in instruments.lib)
+envelopeMod = en.asr(nonLinAttack,1,envelopeRelease,gate);
+
+//nonLinearModultor is declared in instruments.lib, it adapts allpassnn from miscfilter.lib 
+//for using it with waveguide instruments
+NLFM =  nonLinearModulator((nonLinearity : si.smoo),envelopeMod,freq,
+     typeModulation,(frequencyMod : si.smoo),nlfOrder);
+
+//----------------------- Synthesis parameters computing and functions declaration ----------------------------
+
+//reed table parameters
+reedTableOffset = 0.7;
+reedTableSlope = -0.44 + (0.26*reedStiffness);
+
+//the reed function is declared in instruments.lib
+reedTable = reed(reedTableOffset,reedTableSlope);
+
+//delay line with a length adapted in function of the order of nonlinear filter
+delayLength = ma.SR/freq*0.5 - 1.5 - (nlfOrder*nonLinearity)*(typeModulation < 2);
+delayLine = de.fdelay(4096,delayLength);
+
+//one zero filter used as a allpass: pole is set to -1
+filter = oneZero0(0.5,0.5);
+
+//stereoizer is declared in instruments.lib and implement a stereo spacialisation in function of 
+//the frequency period in number of samples 
+stereo = stereoizer(ma.SR/freq);
+
+//----------------------- Algorithm implementation ----------------------------
+
+//Breath pressure + vibrato + breath noise + envelope (Attack / Decay / Sustain / Release)
+envelope = en.adsr(envelopeAttack,envelopeDecay,1,envelopeRelease,gate)*pressure*0.9;
+
+vibrato = os.osc(vibratoFreq)*vibratoGain*
+	envVibrato(0.1*2*vibratoAttack,0.9*2*vibratoAttack,100,vibratoRelease,gate);
+breath = envelope + envelope*no.noise*noiseGain;
+breathPressure = breath + breath*vibrato;
+
+process =
+	//Commuted Loss Filtering
+	(_,(breathPressure <: _,_) : (filter*-0.95 - _ <: 
+	
+	//Non-Linear Scattering
+	*(reedTable)) + _) ~ 
+	
+	//Delay with Feedback
+	(delayLine : NLFM) : 
+	
+	//scaling and stereo
+	*(gain)*1.5 : stereo : instrReverb; 
+