cyfaust 0.1.0__cp311-cp311-macosx_15_0_arm64.whl

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

@@ -0,0 +1,116 @@
+declare name "flute";
+declare description "Nonlinear WaveGuide Flute";
+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 flute based on Smith algorithm: https://ccrma.stanford.edu/~jos/pasp/Flutes_Recorders_Pipe_Organs.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]") : int;
+
+pressure = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Pressure 
+[2][tooltip:Breath pressure (value bewteen 0 and 1)]",0.9,0,1.5,0.01) : si.smoo;
+breathAmp = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Noise Gain 
+[2][tooltip:Breath noise gain (value between 0 and 1)]",0.1,0,1,0.01)/10;
+
+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.1,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.2,0,2,0.01);
+
+pressureEnvelope = checkbox("h:Envelopes_and_Vibrato/v:Pressure_Envelope_Parameters/Pressure_Env 
+[5][unit:s][tooltip:Activate Pressure envelope]") : int;
+env1Attack = hslider("h:Envelopes_and_Vibrato/v:Pressure_Envelope_Parameters/Press_Env_Attack 
+[5][unit:s][tooltip:Pressure envelope attack duration]",0.05,0,2,0.01);
+env1Decay = hslider("h:Envelopes_and_Vibrato/v:Pressure_Envelope_Parameters/Press_Env_Decay 
+[5][unit:s][tooltip:Pressure envelope decay duration]",0.2,0,2,0.01);
+env1Release = hslider("h:Envelopes_and_Vibrato/v:Pressure_Envelope_Parameters/Press_Env_Release 
+[5][unit:s][tooltip:Pressure envelope release duration]",1,0,2,0.01);
+
+env2Attack = hslider("h:Envelopes_and_Vibrato/v:Global_Envelope_Parameters/Glob_Env_Attack 
+[6][unit:s][tooltip:Global envelope attack duration]",0.1,0,2,0.01);
+env2Release = hslider("h:Envelopes_and_Vibrato/v:Global_Envelope_Parameters/Glob_Env_Release 
+[6][unit:s][tooltip:Global 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,0.1,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 ----------------------------
+
+//Loops feedbacks gains
+feedBack1 = 0.4;
+feedBack2 = 0.4;
+
+//Delay Lines
+embouchureDelayLength = (ma.SR/freq)/2-2;
+boreDelayLength = ma.SR/freq-2;
+embouchureDelay = de.fdelay(4096,embouchureDelayLength);
+boreDelay = de.fdelay(4096,boreDelayLength);
+
+//Polynomial
+poly = _ <: _ - _*_*_;
+
+//jet filter is a lowpass filter (declared in miscfilter.lib)
+reflexionFilter = fi.lowpass(1,2000);
+
+//stereoizer is declared in instruments.lib and implement a stereo spatialisation in function of 
+//the frequency period in number of samples 
+stereo = stereoizer(ma.SR/freq);
+
+//----------------------- Algorithm implementation ----------------------------
+
+//Pressure envelope
+env1 = en.adsr(env1Attack,env1Decay,0.9,env1Release,(gate | pressureEnvelope))*pressure*1.1; 
+
+//Global envelope
+env2 = en.asr(env2Attack,1,env2Release,gate)*0.5;
+
+//Vibrato Envelope
+vibratoEnvelope = envVibrato(vibratoBegin,vibratoAttack,100,vibratoRelease,gate)*vibratoGain; 
+
+vibrato = os.osc(vibratoFreq)*vibratoEnvelope;
+
+breath = no.noise*env1;
+
+flow = env1 + breath*breathAmp + vibrato;
+
+//instrReverb is declared in instruments.lib
+process = (_ <: (flow + *(feedBack1) : embouchureDelay : poly) + *(feedBack2) : reflexionFilter)~(boreDelay : NLFM) : *(env2)*gain : 
+stereo : instrReverb;
+

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

@@ -0,0 +1,121 @@
+declare name "fluteStk";
+declare description "Nonlinear WaveGuide Flute from STK";
+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 flute physical model, as discussed by Karjalainen, Smith, Waryznyk, etc.  The jet model uses a polynomial, a la Cook.";
+declare reference "https://ccrma.stanford.edu/~jos/pasp/Flutes_Recorders_Pipe_Organs.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]");
+
+embouchureAjust = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Embouchure_Ajust
+[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
+[2][tooltip:A value between 0 and 1]",0.03,0,1,0.01);
+pressure = hslider("h:Physical_and_Nonlinearity/v:Physical_Parameters/Pressure
+[2][tooltip:Breath pressure (value between 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]",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.03,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.3,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 ----------------------------
+
+jetReflexion = 0.5;
+//jetRatio = 0.08 + (0.48*embouchureAjust); //original stk function
+jetRatio = 1+(0.5-embouchureAjust); //corrected function
+endReflexion = 0.5;
+
+//Delay lines lengths in number of samples
+//jetDelayLength = (SR/freq-2)*jetRatio; //original stk function for jet delay length
+jetDelayLength = (ma.SR/(freq*2)-2)*jetRatio; //corrected function for jet delay length
+boreDelayLength = ma.SR/(freq*2)-2; //original function for bore delay length
+//boreDelayLength = SR/(freq)-2; //corrected function for bore delay length
+filterPole = 0.7 - (0.1*22050/ma.SR);
+
+//One Pole Filter (declared in instruments.lib)
+onePoleFilter = _*gain : onePole(b0,a1)
+	with{
+		gain = -1;
+		pole = 0.7 - (0.1*22050/ma.SR);
+		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 ----------------------------
+
+//the vibrato amplitude is controlled by an envelope generator (declared in instruments.lib)
+vibrato = vibratoGain*envVibrato(vibratoBegin,vibratoAttack,100,vibratoRelease,gate)*os.osc(vibratoFreq);
+
+//Breath pressure is controlled by an Attack / Decay / Sustain / Release envelope
+envelopeBreath = pressure*en.adsr(pressure*envelopeAttack,envelopeDecay,0.8,envelopeRelease,gate);
+breathPressure = envelopeBreath + envelopeBreath*(noiseGain*no.noise + vibrato) + 10.0^(-15.0);
+
+//delay lines
+jetDelay = de.fdelay(4096,jetDelayLength);
+boreDelay = de.fdelay(4096,boreDelayLength);
+
+//reflexion filter is a one pole and a dcblocker
+reflexionFilters = onePoleFilter : fi.dcblocker;
+
+process =
+	(reflexionFilters <: 
+	//Differential Pressure
+	((breathPressure - _*jetReflexion) : 
+	jetDelay : jetTable) + (_*endReflexion)) ~ (boreDelay : NLFM) : 
+	//output scaling and stereo signal
+	*(0.3*gain) : stereo : instrReverb; 
+
+

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

@@ -0,0 +1,131 @@
+declare name "glassHarmonica";
+declare description "Nonlinear Banded Waveguide Modeled Glass Harmonica";
+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 = 3;
+
+nMode(3) = 6;
+
+modes(3,0) = 1.0;
+basegains(3,0) = pow(0.999,1);
+excitation(3,0) = 1*gain*gate/(nMode(3) - 1);
+
+modes(3,1) = 2.32;
+basegains(3,1) = pow(0.999,2);
+excitation(3,1) = 1*gain*gate/(nMode(3) - 1);
+
+modes(3,2) = 4.25;
+basegains(3,2) = pow(0.999,3);
+excitation(3,2) = 1*gain*gate/(nMode(3) - 1);
+
+modes(3,3) = 6.63;
+basegains(3,3) = pow(0.999,4);
+excitation(3,3) = 1*gain*gate/(nMode(3) - 1);
+
+modes(3,4) = 9.38;
+basegains(3,4) = pow(0.999,5);
+excitation(3,4) = 1*gain*gate/(nMode(3) - 1);
+
+modes(3,5) = 9 : float;
+basegains(3,5) = 0 : float;
+excitation(3,5) = 0 : float;
+
+//==================== 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
+		*(4) : NLFM : stereo : instrReverb;
+

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

@@ -0,0 +1,90 @@
+declare name "harpsi";
+declare description "Nonlinear WaveGuide Commuted Harpsichord";
+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 commuted WaveGuide Harpsichord."; 
+
+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]");
+
+typeModulation = nentry("h:Physical_and_Nonlinearity/v:Nonlinear_Filter_Parameters/Modulation_Type 
+[2][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 
+[2][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 
+[2][unit:Hz][tooltip:Frequency of the sine wave for the modulation of theta (works if Modulation Type=3)]",220,20,1000,0.1);
+
+//==================== 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 ----------------------------
+
+//convert a frequency in a midi note number
+freqToNoteNumber = (log-log(440))/log(2)*12+69+0.5 : int;
+freqn = freq : freqToNoteNumber;
+
+//string excitation
+soundBoard = dryTapAmp*no.noise
+	with{
+		dryTapAmpT60 = ffunction(float getValueDryTapAmpT60(float), <harpsichord.h>,"");
+		noteCutOffTime = freqn : dryTapAmpT60*gain;
+		dryTapAmp = asympT60(0.15,0,noteCutOffTime,gate);
+	};	
+
+//loopfilter is a biquad filter whose coefficients are extracted from a C++ file using the foreign function mechanism
+loopFilter = fi.TF2(b0,b1,b2,a1,a2)
+	   with{
+		//functions are imported from the C++ file
+		loopFilterb0 = ffunction(float getValueLoopFilterb0(float), <harpsichord.h>,"");
+		loopFilterb1 = ffunction(float getValueLoopFilterb1(float), <harpsichord.h>,"");
+		loopFilterb2 = ffunction(float getValueLoopFilterb2(float), <harpsichord.h>,"");
+		loopFiltera1 = ffunction(float getValueLoopFiltera1(float), <harpsichord.h>,"");
+		loopFiltera2 = ffunction(float getValueLoopFiltera2(float), <harpsichord.h>,"");
+		//coefficients are extracted from the functions
+		b0 = loopFilterb0(freqn);
+		b1 = loopFilterb1(freqn);
+		b2 = loopFilterb2(freqn);
+		a1 = loopFiltera1(freqn);
+		a2 = loopFiltera2(freqn);
+	   };
+
+//delay length as a number of samples
+delayLength = 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);
+
+//----------------------- Algorithm implementation ----------------------------
+
+//envelope for string loop resonance time
+stringLoopGainT = gate*0.9996 + (gate<1)*releaseLoopGain(freqn)*0.9 : si.smoo
+		with{
+			releaseLoopGain = ffunction(float getValueReleaseLoopGain(float), <harpsichord.h>,"");
+		};
+
+//one string
+string = (*(stringLoopGainT)+_ : de.delay(4096,delayLength) : loopFilter)~NLFM;
+
+process = soundBoard : string : stereo : instrReverb;
+
+

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

@@ -0,0 +1,185 @@
+//Parameters for harpsichord.dsp
+//©Romain Michon (rmichon@ccrma.stanford.edu), 2011
+//licence: FAUST-STK
+
+#include "instrument.h"
+
+#define _LOOKUP_TABLE_H_
+
+double DryTapAmpT60_points[16*2] = {
+//  NoteNumber(A440=69), T60(sec)
+	21.001,0.491,
+	26.587,0.498,
+	34.249,0.470,
+	40.794,0.441,
+	47.977,0.392,
+	55.000,0.370,
+	60.000,0.370,
+	66.000,0.370,
+	71.934,0.370,
+	78.000,0.370,
+	83.936,0.390,
+	88.557,0.387, 
+	92.858,0.400,
+	97.319,0.469,
+	102.400,0.500,
+	107.198,0.494	};
+extern LookupTable DryTapAmpT60;
+LookupTable DryTapAmpT60(&DryTapAmpT60_points[0], 16);
+
+float getValueDryTapAmpT60 (float index){
+	return DryTapAmpT60.getValue(index);
+}
+
+double releaseLoopGain_points[10*2] = {
+	21.000,0.865,
+	24.000,0.880,
+	29.000,0.896,
+	36.000,0.910,
+	48.000,0.920,
+	60.000,0.950,
+	72.000,0.965,
+	84.000,0.988,
+	88.000,0.997,
+	99.000,0.988	};
+extern LookupTable releaseLoopGain;
+LookupTable releaseLoopGain(&releaseLoopGain_points[0], 10);
+
+float getValueReleaseLoopGain(float index){
+	return releaseLoopGain.getValue(index);
+}
+
+double loopFilterb0_points[18*2] = {
+	35.000,0.94373,
+	36.000,0.94731,
+	46.000,0.94731,
+	47.000,0.96202,
+	52.000,0.96202,
+	53.000,0.97477,
+	58.000,0.97477,
+	59.000,0.97733,
+	64.000,0.97733,
+	65.000,0.97971,
+	70.000,0.97971,
+	71.000,0.97971,
+	76.000,0.97971,
+	77.000,0.98698,
+	82.000,0.98698,
+	83.000,0.98462,
+	86.000,0.98462,
+	87.000,0.98611,
+};
+extern LookupTable loopFilterb0;
+LookupTable loopFilterb0(&loopFilterb0_points[0], 18);
+
+float getValueLoopFilterb0(float index){
+	return loopFilterb0.getValue(index);
+}
+
+double loopFilterb1_points[18*2] = {
+	35.000,0.60010,
+	36.000,-0.59124,
+	46.000,-0.59124,
+	47.000,-0.21243,
+	52.000,-0.21243,
+	53.000,-0.39280,
+	58.000,-0.39280,
+	59.000,-0.48307,
+	64.000,-0.48307,
+	65.000,0.51965,
+	70.000,0.51965,
+	71.000,0.51965,
+	76.000,0.51965,
+	77.000,-0.42463,
+	82.000,-0.42463,
+	83.000,0.85655,
+	86.000,0.85655,
+	87.000,0.68851,
+};
+extern LookupTable loopFilterb1;
+LookupTable loopFilterb1(&loopFilterb1_points[0], 18);
+
+float getValueLoopFilterb1(float index){
+	return loopFilterb1.getValue(index);
+}
+
+double loopFilterb2_points[18*2] = {
+	35.000,-0.00360,
+	36.000,-0.12249,
+	46.000,-0.12249,
+	47.000,-0.16044,
+	52.000,-0.16044,
+	53.000,-0.21680,
+	58.000,-0.21680,
+	59.000,-0.16346,
+	64.000,-0.16346,
+	65.000,0.22162,
+	70.000,0.22162,
+	71.000,0.22162,
+	76.000,0.22162,
+	77.000,-0.14973,
+	82.000,-0.14973,
+	83.000,0.24937,
+	86.000,0.24937,
+	87.000,0.14838,
+};
+extern LookupTable loopFilterb2;
+LookupTable loopFilterb2(&loopFilterb2_points[0], 18);
+
+float getValueLoopFilterb2(float index){
+	return loopFilterb2.getValue(index);
+}
+
+double loopFiltera1_points[18*2] = {
+	35.000,0.5941,
+	36.000,-0.65928,
+	46.000,-0.65928,
+	47.000,-0.24222,
+	52.000,-0.24222,
+	53.000,-0.41402,
+	58.000,-0.41402,
+	59.000,-0.50837,
+	64.000,-0.50837,
+	65.000,0.51263,
+	70.000,0.51263,
+	71.000,0.51263,
+	76.000,0.51263,
+	77.000,-0.43976,
+	82.000,-0.43976,
+	83.000,0.85396,
+	86.000,0.85396,
+	87.000,0.68332,
+};
+extern LookupTable loopFiltera1;
+LookupTable loopFiltera1(&loopFiltera1_points[0], 18);
+
+float getValueLoopFiltera1(float index){
+	return loopFiltera1.getValue(index);
+}
+
+double loopFiltera2_points[18*2] = {
+	35.000,-0.02641,
+	36.000,-0.10275,
+	46.000,-0.10275,
+	47.000,-0.15842,
+	52.000,-0.15842,
+	53.000,-0.21653,
+	58.000,-0.21653,
+	59.000,-0.15833,
+	64.000,-0.15833,
+	65.000,0.22025,
+	70.000,0.22025,
+	71.000,0.22025,
+	76.000,0.22025,
+	77.000,-0.14583,
+	82.000,-0.14583,
+	83.000,0.24405,
+	86.000,0.24405,
+	87.000,0.14370,
+};
+extern LookupTable loopFiltera2;
+LookupTable loopFiltera2(&loopFiltera2_points[0], 18);
+
+float getValueLoopFiltera2(float index){
+	return loopFiltera2.getValue(index);
+}