cyfaust 0.1.0__cp311-cp311-macosx_15_0_arm64.whl

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

@@ -0,0 +1,114 @@
+declare name "saxophony";
+declare description "Nonlinear WaveGuide Saxophone";
+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 "This class implements a hybrid digital waveguide instrument that can generate a variety of wind-like sounds.  It has also been referred to as the blowed string model. The waveguide section is essentially that of a string, with one rigid and one lossy termination.  The non-linear function is a reed table. The string can be blown at any point between the terminations, though just as with strings, it is impossible to excite the system at either end. If the excitation is placed at the string mid-point, the sound is that of a clarinet.  At points closer to the bridge, the sound is closer to that of a saxophone.  See Scavone (2002) for more details.";
+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]");
+
+pressure = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Pressure
+[2][tooltip:Breath pressure (a value between 0 and 1)]",1,0,1,0.01);
+reedStiffness = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Reed_Stiffness
+[2][tooltip:A value between 0 and 1]",0.3,0,1,0.01);
+blowPosition = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Blow_Position
+[2][tooltip:A value between 0 and 1]",0.5,0,1,0.01);
+noiseGain = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Noise_Gain",0.05,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.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.3,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.05,0,2,0.01);
+envelopeRelease = hslider("h:Envelopes_and_Vibrato/v:Envelope_Parameters/Envelope_Release 
+[5][unit:s][tooltip:Envelope release duration]",0.01,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 ----------------------------
+
+//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);
+
+//reed table parameters
+reedTableOffset = 0.7;
+reedTableSlope = 0.1 + (0.4*reedStiffness);
+
+//the reed function is declared in instruments.lib
+reedTable = reed(reedTableOffset,reedTableSlope);
+
+//Delay lines length in number of samples
+fdel1 = (1-blowPosition) * (ma.SR/freq - 3);
+fdel2 = (ma.SR/freq - 3)*blowPosition +1 ;
+
+//Delay lines
+delay1 = de.fdelay(4096,fdel1);
+delay2 = de.fdelay(4096,fdel2);
+
+//Breath pressure is controlled by an attack / sustain / release envelope (asr is declared in instruments.lib)
+envelope = (0.55+pressure*0.3)*en.asr(pressure*envelopeAttack,1,pressure*envelopeRelease,gate);
+breath = envelope + envelope*noiseGain*no.noise;
+
+//envVibrato is decalred in instruments.lib
+vibrato = vibratoGain*envVibrato(vibratoBegin,vibratoAttack,100,vibratoRelease,gate)*osc(vibratoFreq);
+breathPressure = breath + breath*vibratoGain*os.osc(vibratoFreq);
+
+//Body filter is a one zero filter (declared in instruments.lib)
+bodyFilter = *(gain) : oneZero1(b0,b1)
+	with{
+		gain = -0.95;
+		b0 = 0.5;
+		b1 = 0.5;	
+	};
+
+instrumentBody(delay1FeedBack,breathP) = delay1FeedBack <: -(delay2) <: 
+	((breathP - _ <: breathP - _*reedTable) - delay1FeedBack),_;
+
+process =
+	(bodyFilter,breathPressure : instrumentBody) ~ 
+	(delay1 : NLFM) : !,
+	//Scaling Output and stereo
+	*(gain) : stereo : instrReverb;
+

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

@@ -0,0 +1,48 @@
+declare name "Sitar";
+declare description "WaveGuide Sitar";
+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 instrument implements a sitar plucked string physical model based on the Karplus-Strong algorithm using a randomly modulated delay line.";
+
+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]");
+
+resonance = hslider("v:Physical_Parameters/Resonance
+[2][tooltip:A value between 0 and 1]",0.7,0,1,0.01)*0.1;
+
+//==================== SIGNAL PROCESSING ================
+
+//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);
+
+//excitation envelope (adsr)
+envelope = en.adsr(0.001,0.04,1,0.5,gate);
+
+//the delay length is randomly modulated
+targetDelay = ma.SR/freq;
+delayLength = targetDelay*((1+(0.5*no.noise)) : si.smooth(0.9992));
+delayLine = de.delay(4096,delayLength);
+
+//the loop gain control the resonance duration
+loopGain = 0.895 + resonance + (freq*0.0000005);
+amGain = 0.1*gain;
+
+//feedback filter is a one zero (declared in instruments.lib)
+filter = oneZero1(b0,b1)
+	with{
+		zero = 0.01;
+		b0 = 1/(1 + zero);
+		b1 = -zero*b0;
+	};
+
+process = (*(loopGain) : filter + (envelope*no.noise*amGain))~delayLine : *(8) : 
+stereo : instrReverb;
+

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

@@ -0,0 +1,155 @@
+declare name "tibetanBowl";
+declare description "Banded Waveguide Modeld Tibetan Bowl";
+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 "This instrument uses banded waveguide. For more information, see Essl, G. and Cook, P. Banded Waveguides: Towards Physical Modelling of Bar Percussion Instruments, Proceedings of the 1999 International Computer Music Conference.";
+
+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)]",0.8,0,1,0.01); 
+gate = button("h:Basic_Parameters/gate [1][tooltip:noteOn = 1, noteOff = 0]");
+
+select = nentry("h:Physical_and_Nonlinearity/v:Physical_Parameters/Excitation_Selector
+[2][tooltip:0=Bow; 1=Strike]",0,0,1,1);
+integrationConstant = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Integration_Constant
+[2][tooltip:A value between 0 and 1]",0,0,1,0.01);
+baseGain = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Base_Gain
+[2][tooltip:A value between 0 and 1]",1,0,1,0.01);
+bowPressure = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Bow_Pressure
+[2][tooltip:Bow pressure on the instrument (Value between 0 and 1)]",0.2,0,1,0.01);
+bowPosition = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Bow_Position
+[2][tooltip:Bow position on the instrument (Value between 0 and 1)]",0,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);
+
+//==================== MODAL PARAMETERS ================
+
+preset = 0;
+
+nMode(0) = 12;
+
+modes(0,0) = 0.996108344;
+basegains(0,0) = 0.999925960128219;
+excitation(0,0) = 11.900357 / 10;
+    
+modes(0,1) = 1.0038916562;
+basegains(0,1) = 0.999925960128219;
+excitation(0,1) = 11.900357 / 10;
+
+modes(0,2) = 2.979178;
+basegains(0,2) = 0.999982774366897;
+excitation(0,2) = 10.914886 / 10;
+
+modes(0,3) = 2.99329767;
+basegains(0,3) = 0.999982774366897;
+excitation(0,3) = 10.914886 / 10;
+    
+modes(0,4) = 5.704452;
+basegains(0,4) = 1.0;
+excitation(0,4) = 42.995041 / 10;
+    
+modes(0,5) = 5.704452;
+basegains(0,5) = 1.0;
+excitation(0,5) = 42.995041 / 10;
+    
+modes(0,6) = 8.9982;
+basegains(0,6) = 1.0;
+excitation(0,6) = 40.063034 / 10;
+    
+modes(0,7) = 9.01549726;
+basegains(0,7) = 1.0;
+excitation(0,7) = 40.063034 / 10;
+    
+modes(0,8) = 12.83303;
+basegains(0,8) = 0.999965497558225;
+excitation(0,8) = 7.063034 / 10;
+   
+modes(0,9) = 12.807382;
+basegains(0,9) = 0.999965497558225;
+excitation(0,9) = 7.063034 / 10;
+    
+modes(0,10) = 17.2808219;
+basegains(0,10) = 0.9999999999999999999965497558225;
+excitation(0,10) = 57.063034 / 10;
+    
+modes(0,11) = 21.97602739726;
+basegains(0,11) = 0.999999999999999965497558225;
+excitation(0,11) = 57.063034 / 10;
+
+//==================== 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 ----------------------------
+
+//the number of modes depends on the preset being used
+nModes = nMode(preset);
+
+//bow table parameters
+tableOffset = 0;
+tableSlope = 10 - (9*bowPressure);
+
+delayLengthBase = ma.SR/freq;
+
+//delay lengths in number of samples
+delayLength(x) = delayLengthBase/modes(preset,x);
+
+//delay lines
+delayLine(x) = de.delay(4096,delayLength(x));
+
+//Filter bank: bandpass filters (declared in instruments.lib)
+radius = 1 - ma.PI*32/ma.SR;
+bandPassFilter(x) = bandPass(freq*modes(preset,x),radius);
+
+//Delay lines feedback for bow table lookup control
+baseGainApp = 0.8999999999999999 + (0.1*baseGain);
+velocityInputApp = integrationConstant;
+velocityInput = velocityInputApp + _*baseGainApp,par(i,(nModes-1),(_*baseGainApp)) :> +;
+
+//Bow velocity is controlled by an ADSR envelope
+maxVelocity = 0.03 + 0.1*gain;
+bowVelocity = maxVelocity*en.adsr(0.02,0.005,1,0.01,gate);
+
+//stereoizer is declared in instruments.lib and implement a stereo spacialisation in function of 
+//the frequency period in number of samples 
+stereo = stereoizer(delayLengthBase);
+
+//----------------------- Algorithm implementation ----------------------------
+
+//Bow table lookup (bow is decalred in instruments.lib)
+bowing = bowVelocity - velocityInput <: *(bow(tableOffset,tableSlope)) : /(nModes);
+
+//One resonance
+resonance(x) = + : + (excitation(preset,x)*select) : delayLine(x) : *(basegains(preset,x)) : bandPassFilter(x);
+
+process =
+		//Bowed Excitation
+		(bowing*((select-1)*-1) <:
+		//nModes resonances with nModes feedbacks for bow table look-up 
+		par(i,nModes,(resonance(i)~_)))~par(i,nModes,_) :> + : 
+		//Signal Scaling and stereo
+		NLFM : stereo : instrReverb;
+

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

@@ -0,0 +1,123 @@
+declare name "tunedBar";
+declare description "Nonlinear Banded Waveguide Models";
+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 "This instrument uses banded waveguide. For more information, see Essl, G. and Cook, P. Banded Waveguides: Towards Physical Modelling of Bar Percussion Instruments, Proceedings of the 1999 International Computer Music Conference.";
+
+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)]",0.8,0,1,0.01); 
+gate = button("h:Basic_Parameters/gate [1][tooltip:noteOn = 1, noteOff = 0]");
+
+select = nentry("h:Physical_and_Nonlinearity/v:Physical_Parameters/Excitation_Selector
+[2][tooltip:0=Bow; 1=Strike]",0,0,1,1);
+integrationConstant = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Integration_Constant
+[2][tooltip:A value between 0 and 1]",0,0,1,0.01);
+baseGain = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Base_Gain
+[2][tooltip:A value between 0 and 1]",1,0,1,0.01);
+bowPressure = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Bow_Pressure
+[2][tooltip:Bow pressure on the instrument (Value between 0 and 1)]",0.2,0,1,0.01);
+bowPosition = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Bow_Position
+[2][tooltip:Bow position on the instrument (Value between 0 and 1)]",0,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);
+
+//==================== MODAL PARAMETERS ================
+
+preset = 2;
+
+nMode(2) = 4;
+
+modes(2,0) = 1;
+basegains(2,0) = pow(0.999,1);
+excitation(2,0) = 1*gain*gate/nMode(2);
+
+modes(2,1) = 4.0198391420;
+basegains(2,1) = pow(0.999,2);
+excitation(2,1) = 1*gain*gate/nMode(2);
+
+modes(2,2) = 10.7184986595;
+basegains(2,2) = pow(0.999,3);
+excitation(2,2) = 1*gain*gate/nMode(2);
+
+modes(2,3) = 18.0697050938;
+basegains(2,3) = pow(0.999,4);
+excitation(2,3) = 1*gain*gate/nMode(2);
+
+//==================== 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 ----------------------------
+
+//the number of modes depends on the preset being used
+nModes = nMode(preset);
+
+//bow table parameters
+tableOffset = 0;
+tableSlope = 10 - (9*bowPressure);
+
+delayLengthBase = ma.SR/freq;
+
+//delay lengths in number of samples
+delayLength(x) = delayLengthBase/modes(preset,x);
+
+//delay lines
+delayLine(x) = de.delay(4096,delayLength(x));
+
+//Filter bank: bandpass filters (declared in instruments.lib)
+radius = 1 - ma.PI*32/ma.SR;
+bandPassFilter(x) = bandPass(freq*modes(preset,x),radius);
+
+//Delay lines feedback for bow table lookup control
+baseGainApp = 0.8999999999999999 + (0.1*baseGain);
+velocityInputApp = integrationConstant;
+velocityInput = velocityInputApp + _*baseGainApp,par(i,(nModes-1),(_*baseGainApp)) :> +;
+
+//Bow velocity is controlled by an ADSR envelope
+maxVelocity = 0.03 + 0.1*gain;
+bowVelocity = maxVelocity*en.adsr(0.02,0.005,0.9,0.01,gate);
+
+//stereoizer is declared in instruments.lib and implement a stereo spacialisation in function of 
+//the frequency period in number of samples 
+stereo = stereoizer(delayLengthBase);
+
+//----------------------- Algorithm implementation ----------------------------
+
+//Bow table lookup (bow is decalred in instruments.lib)
+bowing = bowVelocity - velocityInput <: *(bow(tableOffset,tableSlope)) : /(nModes);
+
+//One resonance
+resonance(x) = + : + (excitation(preset,x)*select) : delayLine(x) : *(basegains(preset,x)) : bandPassFilter(x);
+
+process =
+		//Bowed Excitation
+		(bowing*((select-1)*-1) <:
+		//nModes resonances with nModes feedbacks for bow table look-up 
+		par(i,nModes,(resonance(i)~_)))~par(i,nModes,_) :> + : 
+		//Signal Scaling and stereo
+		*(4) : NLFM : stereo : instrReverb;
+

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

@@ -0,0 +1,100 @@
+declare name "UniBar";
+declare description "Nonlinear Banded Waveguide Models";
+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 "This instrument uses banded waveguide. For more information, see Essl, G. and Cook, P. Banded Waveguides: Towards Physical Modelling of Bar Percussion Instruments, Proceedings of the 1999 International Computer Music Conference.";
+
+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)]",0.8,0,1,0.01); 
+gate = button("h:Basic_Parameters/gate [1][tooltip:noteOn = 1, noteOff = 0]");
+
+select = nentry("h:Physical_and_Nonlinearity/v:Physical_Parameters/Excitation_Selector
+[2][tooltip:0=Bow; 1=Strike]",0,0,1,1);
+integrationConstant = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Integration_Constant
+[2][tooltip:A value between 0 and 1]",0,0,1,0.01);
+baseGain = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Base_Gain
+[2][tooltip:A value between 0 and 1]",1,0,1,0.01);
+bowPressure = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Bow_Pressure
+[2][tooltip:Bow pressure on the instrument (Value between 0 and 1)]",0.2,0,1,0.01);
+bowPosition = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Bow_Position
+[2][tooltip:Bow position on the instrument (Value between 0 and 1)]",0,0,1,0.01);
+
+//==================== MODAL PARAMETERS ================
+
+preset = 1;
+
+nMode(1) = 4;
+
+modes(1,0) = 1;
+basegains(1,0) = pow(0.9,1);
+excitation(1,0) = 1*gain/nMode(1);
+
+modes(1,1) = 2.756;
+basegains(1,1) = pow(0.9,2);
+excitation(1,1) = 1*gain/nMode(1);
+
+modes(1,2) = 5.404;
+basegains(1,2) = pow(0.9,3);
+excitation(1,2) = 1*gain/nMode(1);
+
+modes(1,3) = 8.933;
+basegains(1,3) = pow(0.9,4);
+excitation(1,3) = 1*gain/nMode(1);
+
+//==================== SIGNAL PROCESSING ================
+
+//----------------------- Synthesis parameters computing and functions declaration ----------------------------
+
+//the number of modes depends on the preset being used
+nModes = nMode(preset);
+
+//bow table parameters
+tableOffset = 0;
+tableSlope = 10 - (9*bowPressure);
+
+delayLengthBase = ma.SR/freq;
+
+//delay lengths in number of samples
+delayLength(x) = delayLengthBase/modes(preset,x);
+
+//delay lines
+delayLine(x) = de.delay(4096,delayLength(x));
+
+//Filter bank: bandpass filters (declared in instruments.lib)
+radius = 1 - ma.PI*32/ma.SR;
+bandPassFilter(x) = bandPass(freq*modes(preset,x),radius);
+
+//Delay lines feedback for bow table lookup control
+baseGainApp = 0.8999999999999999 + (0.1*baseGain);
+velocityInputApp = integrationConstant;
+velocityInput = velocityInputApp + _*baseGainApp,par(i,(nModes-1),(_*baseGainApp)) :> +;
+
+//Bow velocity is controlled by an ADSR envelope
+maxVelocity = 0.03 + 0.1*gain;
+bowVelocity = maxVelocity*en.adsr(0.02,0.005,0.9,0.01,gate);
+
+//stereoizer is declared in instruments.lib and implement a stereo spacialisation in function of 
+//the frequency period in number of samples 
+stereo = stereoizer(delayLengthBase);
+
+//----------------------- Algorithm implementation ----------------------------
+
+//Bow table lookup (bow is decalred in instruments.lib)
+bowing = bowVelocity - velocityInput <: *(bow(tableOffset,tableSlope)) : /(nModes);
+
+//One resonance
+resonance(x) = + : + (excitation(preset,x)*select) : delayLine(x) : *(basegains(preset,x)) : bandPassFilter(x);
+
+process =
+		//Bowed Excitation
+		(bowing*((select-1)*-1) <:
+		//nModes resonances with nModes feedbacks for bow table look-up 
+		par(i,nModes,(resonance(i)~_)))~par(i,nModes,_) :> + : 
+		//Signal Scaling and stereo
+		*(14) : stereo : instrReverb;

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

@@ -0,0 +1,125 @@
+declare name "voiceForm";
+declare description "Voice Formant 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 instrument contains an excitation singing wavetable (looping wave with random and periodic vibrato, smoothing on frequency, etc.), excitation noise, and four sweepable complex resonances. Phoneme preset numbers: 0->eee (beet), 1->ihh (bit), 2->ehh (bet), 3->aaa (bat), 4->ahh (father), 5->aww (bought), 6->ohh (bone), 7->uhh (but), 8->uuu (foot), 9->ooo (boot), 10->rrr (bird), 11->lll (lull), 12->mmm (mom), 13->nnn (nun), 14->nng (sang), 15->ngg (bong), 16->fff, 17->sss, 18->thh, 19->shh, 20->xxx, 21->hee (beet), 22->hoo (boot), 23->hah (father), 24->bbb, 25->ddd, 26->jjj, 27->ggg, 28->vvv, 29->zzz, 30->thz, 31->zhh";
+
+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]");
+
+phoneme = hslider("v:Physical_Parameters/Phoneme
+[2][tooltip:0->eee, 1->ihh, 2->ehh, 3->aaa, 4->ahh, 5->aww, 6->ohh, 7->uhh, 8->uuu, 9->ooo, 10->rrr, 11->lll, 12->mmm, 13->nnn, 14->nng, 15->ngg, 16->fff, 17->sss, 18->thh, 19->shh, 20->xxx, 21->hee, 22->hoo, 23->hah, 24->bbb, 25->ddd, 26->jjj, 27->ggg, 28->vvv, 29->zzz, 30->thz, 31->zhh]",4,0,31,1);
+
+vibratoFreq = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Freq 
+[3][unit:Hz]",6,1,15,0.1);
+vibratoGain = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Gain
+[3][tooltip:A value between 0 and 1]",0.05,0,1,0.01);
+vibratoBegin = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Begin
+[3][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 
+[3][unit:s][tooltip:Vibrato attack duration]",0.5,0,2,0.01);
+vibratoRelease = hslider("h:Envelopes_and_Vibrato/v:Vibrato_Parameters/Vibrato_Release 
+[3][unit:s][tooltip:Vibrato release duration]",0.1,0,2,0.01);
+
+voicedEnvelopeAttack = hslider("h:Envelopes_and_Vibrato/v:Envelope_Parameters/Voiced_Attack
+[4][unit:s][tooltip:Voiced sounds attack duration]",0.01,0,2,0.01);
+voicedEnvelopeRelease = hslider("h:Envelopes_and_Vibrato/v:Envelope_Parameters/Voiced_Release
+[4][unit:s][tooltip:Voiced sounds release duration]",0.01,0,2,0.01);
+
+noiseEnvelopeAttack = hslider("h:Envelopes_and_Vibrato/v:Envelope_Parameters/Noised_Attack
+[4][unit:s][tooltip:Noised sounds attack duration]",0.001,0,2,0.001);
+noiseEnvelopeRelease = hslider("h:Envelopes_and_Vibrato/v:Envelope_Parameters/Noised_Release
+[4][unit:s][tooltip:Noised sounds release duration]",0.001,0,2,0.001);
+
+//==================== SIGNAL PROCESSING ================
+
+//----------------------- Synthesis parameters computing and functions declaration ----------------------------
+
+//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);
+
+//exitation filters (declared in instruments.lib)
+onePoleFilter = onePole(b0,a1)
+	with{
+		pole = 0.97 - (gain*0.2);
+		b0 = 1 - pole;
+		a1 = -pole;	
+	};
+oneZeroFilter = oneZero1(b0,b1)
+	with{
+		zero = -0.9;
+		b0 = 1/(1 - zero);
+		b1 = -zero*b0;
+	};
+
+//implements a formant (resonance) which can be "swept" over time from one frequency setting to another
+formSwep(frequency,radius,filterGain) = *(gain_) : bandPass(frequency_,radius)
+	with{
+		//filter's radius, gain and frequency are interpolated
+		radius_ = radius : si.smoo;
+		frequency_ = frequency : si.smoo;
+		gain_ = filterGain : si.smoo;
+	};
+
+//formants parameters are countained in a C++ file
+phonemeGains = ffunction(float loadPhonemeGains(int,int), <phonemes.h>,"");
+phonemeParameters = ffunction(float loadPhonemeParameters(int,int,int), <phonemes.h>,"");
+
+//formants frequencies
+ffreq0 = phonemeParameters(phoneme,0,0);
+ffreq1 = phonemeParameters(phoneme,1,0);
+ffreq2 = phonemeParameters(phoneme,2,0);
+ffreq3 = phonemeParameters(phoneme,3,0);
+
+//formants radius
+frad0 = phonemeParameters(phoneme,0,1);
+frad1 = phonemeParameters(phoneme,1,1);
+frad2 = phonemeParameters(phoneme,2,1);
+frad3 = phonemeParameters(phoneme,3,1);
+
+//formants gains
+fgain0 = phonemeParameters(phoneme,0,2) : pow(10,(_/20));
+fgain1 = phonemeParameters(phoneme,1,2) : pow(10,(_/20));
+fgain2 = phonemeParameters(phoneme,2,2) : pow(10,(_/20));
+fgain3 = phonemeParameters(phoneme,3,2) : pow(10,(_/20));
+
+//gain of the voiced part od the sound
+voiceGain = phonemeGains(phoneme,0) : si.smoo;
+
+//gain of the fricative part of the sound 
+noiseGain = phonemeGains(phoneme,1) : si.smoo;
+
+//formants filters
+filter0 = formSwep(ffreq0,frad0,fgain0);
+filter1 = formSwep(ffreq1,frad1,fgain1);
+filter2 = formSwep(ffreq2,frad2,fgain2);
+filter3 = formSwep(ffreq3,frad3,fgain3);
+
+//----------------------- Algorithm implementation ----------------------------
+
+//envelopes (declared in instruments.lib) and vibrato
+vibratoEnvelope = envVibrato(vibratoBegin,vibratoAttack,100,vibratoRelease,gate);
+voicedEnvelope = en.asr(voicedEnvelopeAttack,1,voicedEnvelopeRelease,gate);
+noiseEnvelope = en.asr(noiseEnvelopeAttack,1,noiseEnvelopeRelease,gate);
+vibrato = os.osc(vibratoFreq)*vibratoGain*100*vibratoEnvelope;
+
+//the voice source is generated by an impulse train
+//(imptrain defined in oscillator.lib) that is lowpass filtered
+voiced = os.imptrain(freq+vibrato) : fi.lowpass3e(3300) : *(voiceGain*voicedEnvelope);
+
+//ficative sounds are produced by a noise generator
+frica = no.noise*noiseEnvelope*noiseGain;
+
+process = voiced : oneZeroFilter : onePoleFilter : 
+		 +(frica) <: filter0,filter1,filter2,filter3 :> + : stereo : 
+		 instrReverb;
+
+