Updates (in case my plane crashes)

3 files changed, 898 insertions, 0 deletions

diff --git a/plugins/LinuxVST/src/Air/Air.cpp b/plugins/LinuxVST/src/Air/Air.cpp
new file mode 100755
index 0000000..ebe8adc
--- /dev/null
+++ b/plugins/LinuxVST/src/Air/Air.cpp
@@ -0,0 +1,205 @@
+/* ========================================
+ *  Air - Air.h
+ *  Copyright (c) 2016 airwindows, All rights reserved
+ * ======================================== */
+
+#ifndef __Air_H
+#include "Air.h"
+#endif
+
+AudioEffect* createEffectInstance(audioMasterCallback audioMaster) {return new Air(audioMaster);}
+
+Air::Air(audioMasterCallback audioMaster) :
+    AudioEffectX(audioMaster, kNumPrograms, kNumParameters)
+{
+	airPrevAL = 0.0;
+	airEvenAL = 0.0;
+	airOddAL = 0.0;
+	airFactorAL = 0.0;
+	airPrevBL = 0.0;
+	airEvenBL = 0.0;
+	airOddBL = 0.0;
+	airFactorBL = 0.0;
+	airPrevCL = 0.0;
+	airEvenCL = 0.0;
+	airOddCL = 0.0;
+	airFactorCL = 0.0;
+	tripletPrevL = 0.0;
+	tripletMidL = 0.0;
+	tripletAL = 0.0;
+	tripletBL = 0.0;
+	tripletCL = 0.0;
+	tripletFactorL = 0.0;
+
+	airPrevAR = 0.0;
+	airEvenAR = 0.0;
+	airOddAR = 0.0;
+	airFactorAR = 0.0;
+	airPrevBR = 0.0;
+	airEvenBR = 0.0;
+	airOddBR = 0.0;
+	airFactorBR = 0.0;
+	airPrevCR = 0.0;
+	airEvenCR = 0.0;
+	airOddCR = 0.0;
+	airFactorCR = 0.0;
+	tripletPrevR = 0.0;
+	tripletMidR = 0.0;
+	tripletAR = 0.0;
+	tripletBR = 0.0;
+	tripletCR = 0.0;
+	tripletFactorR = 0.0;
+	
+	flipA = false;
+	flipB = false;
+	flop = false;
+	count = 1;
+
+	A = 0.5;
+	B = 0.5;
+	C = 0.5;
+	D = 0.0;
+	E = 1.0;
+	F = 1.0;
+	fpNShapeLA = 0.0;
+	fpNShapeLB = 0.0;
+	fpNShapeRA = 0.0;
+	fpNShapeRB = 0.0;
+	fpFlip = true;
+	//this is reset: values being initialized only once. Startup values, whatever they are.
+	
+    _canDo.insert("plugAsChannelInsert"); // plug-in can be used as a channel insert effect.
+    _canDo.insert("plugAsSend"); // plug-in can be used as a send effect.
+    _canDo.insert("x2in2out"); 
+    setNumInputs(kNumInputs);
+    setNumOutputs(kNumOutputs);
+    setUniqueID(kUniqueId);
+    canProcessReplacing();     // supports output replacing
+    canDoubleReplacing();      // supports double precision processing
+	programsAreChunks(true);
+    vst_strncpy (_programName, "Default", kVstMaxProgNameLen); // default program name
+}
+
+Air::~Air() {}
+VstInt32 Air::getVendorVersion () {return 1000;}
+void Air::setProgramName(char *name) {vst_strncpy (_programName, name, kVstMaxProgNameLen);}
+void Air::getProgramName(char *name) {vst_strncpy (name, _programName, kVstMaxProgNameLen);}
+//airwindows likes to ignore this stuff. Make your own programs, and make a different plugin rather than
+//trying to do versioning and preventing people from using older versions. Maybe they like the old one!
+
+static float pinParameter(float data)
+{
+	if (data < 0.0f) return 0.0f;
+	if (data > 1.0f) return 1.0f;
+	return data;
+}
+
+VstInt32 Air::getChunk (void** data, bool isPreset)
+{
+	float *chunkData = (float *)calloc(kNumParameters, sizeof(float));
+	chunkData[0] = A;
+	chunkData[1] = B;
+	chunkData[2] = C;
+	chunkData[3] = D;
+	chunkData[4] = E;
+	chunkData[5] = F;
+	/* Note: The way this is set up, it will break if you manage to save settings on an Intel
+	 machine and load them on a PPC Mac. However, it's fine if you stick to the machine you 
+	 started with. */
+	
+	*data = chunkData;
+	return kNumParameters * sizeof(float);
+}
+
+VstInt32 Air::setChunk (void* data, VstInt32 byteSize, bool isPreset)
+{	
+	float *chunkData = (float *)data;
+	A = pinParameter(chunkData[0]);
+	B = pinParameter(chunkData[1]);
+	C = pinParameter(chunkData[2]);
+	D = pinParameter(chunkData[3]);
+	E = pinParameter(chunkData[4]);
+	F = pinParameter(chunkData[5]);
+	/* We're ignoring byteSize as we found it to be a filthy liar */
+	
+	/* calculate any other fields you need here - you could copy in 
+	 code from setParameter() here. */
+	return 0;
+}
+
+void Air::setParameter(VstInt32 index, float value) {
+    switch (index) {
+        case kParamA: A = value; break;
+        case kParamB: B = value; break; //percent. Using this value, it'll be 0-100 everywhere
+        case kParamC: C = value; break;
+        case kParamD: D = value; break;
+        case kParamE: E = value; break;
+        case kParamF: F = value; break;
+        default: throw; // unknown parameter, shouldn't happen!
+    }
+}
+
+float Air::getParameter(VstInt32 index) {
+    switch (index) {
+        case kParamA: return A; break;
+        case kParamB: return B; break;
+        case kParamC: return C; break;
+        case kParamD: return D; break;
+        case kParamE: return E; break;
+        case kParamF: return F; break;
+        default: break; // unknown parameter, shouldn't happen!
+    } return 0.0; //we only need to update the relevant name, this is simple to manage
+}
+
+void Air::getParameterName(VstInt32 index, char *text) {
+    switch (index) {
+        case kParamA: vst_strncpy (text, "22K tap", kVstMaxParamStrLen); break;
+		case kParamB: vst_strncpy (text, "15K tap", kVstMaxParamStrLen); break;
+		case kParamC: vst_strncpy (text, "11K tap", kVstMaxParamStrLen); break;
+		case kParamD: vst_strncpy (text, "filters Q", kVstMaxParamStrLen); break;
+		case kParamE: vst_strncpy (text, "Output Level", kVstMaxParamStrLen); break;
+		case kParamF: vst_strncpy (text, "Dry/Wet", kVstMaxParamStrLen); break;
+        default: break; // unknown parameter, shouldn't happen!
+    } //this is our labels for displaying in the VST host
+}
+
+void Air::getParameterDisplay(VstInt32 index, char *text) {
+    switch (index) {
+        case kParamA: float2string ((A*2.0)-1.0, text, kVstMaxParamStrLen); break;
+        case kParamB: float2string ((B*2.0)-1.0, text, kVstMaxParamStrLen); break;
+        case kParamC: float2string ((C*2.0)-1.0, text, kVstMaxParamStrLen); break;
+        case kParamD: float2string (D, text, kVstMaxParamStrLen); break;		
+        case kParamE: float2string (E, text, kVstMaxParamStrLen); break;		
+        case kParamF: float2string (F, text, kVstMaxParamStrLen); break;		
+        default: break; // unknown parameter, shouldn't happen!
+	} //this displays the values and handles 'popups' where it's discrete choices
+}
+
+void Air::getParameterLabel(VstInt32 index, char *text) {
+    switch (index) {
+        case kParamA: vst_strncpy (text, " ", kVstMaxParamStrLen); break;
+        case kParamB: vst_strncpy (text, " ", kVstMaxParamStrLen); break;
+        case kParamC: vst_strncpy (text, " ", kVstMaxParamStrLen); break;
+        case kParamD: vst_strncpy (text, " ", kVstMaxParamStrLen); break;
+        case kParamE: vst_strncpy (text, " ", kVstMaxParamStrLen); break;
+        case kParamF: vst_strncpy (text, " ", kVstMaxParamStrLen); break;
+        default: break; // unknown parameter, shouldn't happen!
+    }
+}
+
+VstInt32 Air::canDo(char *text) 
+{ return (_canDo.find(text) == _canDo.end()) ? -1: 1; } // 1 = yes, -1 = no, 0 = don't know
+
+bool Air::getEffectName(char* name) {
+    vst_strncpy(name, "Air", kVstMaxProductStrLen); return true;
+}
+
+VstPlugCategory Air::getPlugCategory() {return kPlugCategEffect;}
+
+bool Air::getProductString(char* text) {
+  	vst_strncpy (text, "airwindows Air", kVstMaxProductStrLen); return true;
+}
+
+bool Air::getVendorString(char* text) {
+  	vst_strncpy (text, "airwindows", kVstMaxVendorStrLen); return true;
+}
diff --git a/plugins/LinuxVST/src/Air/Air.h b/plugins/LinuxVST/src/Air/Air.h
new file mode 100755
index 0000000..00c0b76
--- /dev/null
+++ b/plugins/LinuxVST/src/Air/Air.h
@@ -0,0 +1,119 @@
+/* ========================================
+ *  Air - Air.h
+ *  Created 8/12/11 by SPIAdmin 
+ *  Copyright (c) 2011 __MyCompanyName__, All rights reserved
+ * ======================================== */
+
+#ifndef __Air_H
+#define __Air_H
+
+#ifndef __audioeffect__
+#include "audioeffectx.h"
+#endif
+
+#include <set>
+#include <string>
+#include <math.h>
+
+enum {
+	kParamA = 0,
+	kParamB = 1,
+	kParamC = 2,
+	kParamD = 3,
+	kParamE = 4,
+	kParamF = 5,
+  kNumParameters = 6
+}; //
+
+const int kNumPrograms = 0;
+const int kNumInputs = 2;
+const int kNumOutputs = 2;
+const unsigned long kUniqueId = 'aiyr';    //Change this to what the AU identity is!
+
+class Air : 
+    public AudioEffectX 
+{
+public:
+    Air(audioMasterCallback audioMaster);
+    ~Air();
+    virtual bool getEffectName(char* name);                       // The plug-in name
+    virtual VstPlugCategory getPlugCategory();                    // The general category for the plug-in
+    virtual bool getProductString(char* text);                    // This is a unique plug-in string provided by Steinberg
+    virtual bool getVendorString(char* text);                     // Vendor info
+    virtual VstInt32 getVendorVersion();                          // Version number
+    virtual void processReplacing (float** inputs, float** outputs, VstInt32 sampleFrames);
+    virtual void processDoubleReplacing (double** inputs, double** outputs, VstInt32 sampleFrames);
+    virtual void getProgramName(char *name);                      // read the name from the host
+    virtual void setProgramName(char *name);                      // changes the name of the preset displayed in the host
+	virtual VstInt32 getChunk (void** data, bool isPreset);
+	virtual VstInt32 setChunk (void* data, VstInt32 byteSize, bool isPreset);
+    virtual float getParameter(VstInt32 index);                   // get the parameter value at the specified index
+    virtual void setParameter(VstInt32 index, float value);       // set the parameter at index to value
+    virtual void getParameterLabel(VstInt32 index, char *text);  // label for the parameter (eg dB)
+    virtual void getParameterName(VstInt32 index, char *text);    // name of the parameter
+    virtual void getParameterDisplay(VstInt32 index, char *text); // text description of the current value    
+    virtual VstInt32 canDo(char *text);
+private:
+    char _programName[kVstMaxProgNameLen + 1];
+    std::set< std::string > _canDo;
+	
+	double airPrevAL;
+	double airEvenAL;
+	double airOddAL;
+	double airFactorAL;
+	double airPrevBL;
+	double airEvenBL;
+	double airOddBL;
+	double airFactorBL;
+	double airPrevCL;
+	double airEvenCL;
+	double airOddCL;
+	double airFactorCL;
+	double tripletPrevL;
+	double tripletMidL;
+	double tripletAL;
+	double tripletBL;
+	double tripletCL;
+	double tripletFactorL;
+
+	double airPrevAR;
+	double airEvenAR;
+	double airOddAR;
+	double airFactorAR;
+	double airPrevBR;
+	double airEvenBR;
+	double airOddBR;
+	double airFactorBR;
+	double airPrevCR;
+	double airEvenCR;
+	double airOddCR;
+	double airFactorCR;
+	double tripletPrevR;
+	double tripletMidR;
+	double tripletAR;
+	double tripletBR;
+	double tripletCR;
+	double tripletFactorR;
+	
+	bool flipA;
+	bool flipB;
+	bool flop;
+	int count;
+    
+	long double fpNShapeLA;
+	long double fpNShapeLB;
+	long double fpNShapeRA;
+	long double fpNShapeRB;
+	bool fpFlip;
+	//default stuff
+
+    float A;
+    float B;
+    float C;
+    float D; //parameters. Always 0-1, and we scale/alter them elsewhere.
+	float E;
+	float F;
+
+};
+
+#endif
diff --git a/plugins/LinuxVST/src/Air/AirProc.cpp b/plugins/LinuxVST/src/Air/AirProc.cpp
new file mode 100755
index 0000000..f73e7d8
--- /dev/null
+++ b/plugins/LinuxVST/src/Air/AirProc.cpp
@@ -0,0 +1,574 @@
+/* ========================================
+ *  Air - Air.h
+ *  Copyright (c) 2016 airwindows, All rights reserved
+ * ======================================== */
+
+#ifndef __Air_H
+#include "Air.h"
+#endif
+
+void Air::processReplacing(float **inputs, float **outputs, VstInt32 sampleFrames) 
+{
+    float* in1  =  inputs[0];
+    float* in2  =  inputs[1];
+    float* out1 = outputs[0];
+    float* out2 = outputs[1];
+
+	double hiIntensity = -pow(((A*2.0)-1.0),3)*2;
+	double tripletintensity = -pow(((B*2.0)-1.0),3);
+	double airIntensity = -pow(((C*2.0)-1.0),3)/2;
+	double filterQ = 2.1-D;
+	double output = E;
+	double wet = F;
+	double dry = 1.0-wet;
+	
+	float fpTemp;
+	long double fpOld = 0.618033988749894848204586; //golden ratio!
+	long double fpNew = 1.0 - fpOld;	
+
+	long double inputSampleL;
+	long double inputSampleR;
+	double drySampleL;
+	double drySampleR;
+	double correctionL;
+	double correctionR;
+    
+    while (--sampleFrames >= 0)
+    {
+		inputSampleL = *in1;
+		inputSampleR = *in2;
+		if (inputSampleL<1.2e-38 && -inputSampleL<1.2e-38) {
+			static int noisesource = 0;
+			//this declares a variable before anything else is compiled. It won't keep assigning
+			//it to 0 for every sample, it's as if the declaration doesn't exist in this context,
+			//but it lets me add this denormalization fix in a single place rather than updating
+			//it in three different locations. The variable isn't thread-safe but this is only
+			//a random seed and we can share it with whatever.
+			noisesource = noisesource % 1700021; noisesource++;
+			int residue = noisesource * noisesource;
+			residue = residue % 170003; residue *= residue;
+			residue = residue % 17011; residue *= residue;
+			residue = residue % 1709; residue *= residue;
+			residue = residue % 173; residue *= residue;
+			residue = residue % 17;
+			double applyresidue = residue;
+			applyresidue *= 0.00000001;
+			applyresidue *= 0.00000001;
+			inputSampleL = applyresidue;
+		}
+		if (inputSampleR<1.2e-38 && -inputSampleR<1.2e-38) {
+			static int noisesource = 0;
+			noisesource = noisesource % 1700021; noisesource++;
+			int residue = noisesource * noisesource;
+			residue = residue % 170003; residue *= residue;
+			residue = residue % 17011; residue *= residue;
+			residue = residue % 1709; residue *= residue;
+			residue = residue % 173; residue *= residue;
+			residue = residue % 17;
+			double applyresidue = residue;
+			applyresidue *= 0.00000001;
+			applyresidue *= 0.00000001;
+			inputSampleR = applyresidue;
+			//this denormalization routine produces a white noise at -300 dB which the noise
+			//shaping will interact with to produce a bipolar output, but the noise is actually
+			//all positive. That should stop any variables from going denormal, and the routine
+			//only kicks in if digital black is input. As a final touch, if you save to 24-bit
+			//the silence will return to being digital black again.
+		}
+		drySampleL = inputSampleL;
+		drySampleR = inputSampleR;
+		
+		correctionL = 0.0;
+		correctionR = 0.0;  //from here on down, please add L and R to the code
+		
+		if (count < 1 || count > 3) count = 1;
+		tripletFactorL = tripletPrevL - inputSampleL;
+		tripletFactorR = tripletPrevR - inputSampleR;
+		switch (count)
+		{
+			case 1:
+				tripletAL += tripletFactorL;
+				tripletCL -= tripletFactorL;
+				tripletFactorL = tripletAL * tripletintensity;
+				tripletPrevL = tripletMidL;
+				tripletMidL = inputSampleL;
+				
+				tripletAR += tripletFactorR;
+				tripletCR -= tripletFactorR;
+				tripletFactorR = tripletAR * tripletintensity;
+				tripletPrevR = tripletMidR;
+				tripletMidR = inputSampleR;
+				break;
+			case 2:
+				tripletBL += tripletFactorL;
+				tripletAL -= tripletFactorL;
+				tripletFactorL = tripletBL * tripletintensity;
+				tripletPrevL = tripletMidL;
+				tripletMidL = inputSampleL;
+				
+				tripletBR += tripletFactorR;
+				tripletAR -= tripletFactorR;
+				tripletFactorR = tripletBR * tripletintensity;
+				tripletPrevR = tripletMidR;
+				tripletMidR = inputSampleR;
+				break;
+			case 3:
+				tripletCL += tripletFactorL;
+				tripletBL -= tripletFactorL;
+				tripletFactorL = tripletCL * tripletintensity;
+				tripletPrevL = tripletMidL;
+				tripletMidL = inputSampleL;
+				
+				tripletCR += tripletFactorR;
+				tripletBR -= tripletFactorR;
+				tripletFactorR = tripletCR * tripletintensity;
+				tripletPrevR = tripletMidR;
+				tripletMidR = inputSampleR;
+				break;
+		}
+		tripletAL /= filterQ;
+		tripletBL /= filterQ;
+		tripletCL /= filterQ;
+		correctionL = correctionL + tripletFactorL;
+
+		tripletAR /= filterQ;
+		tripletBR /= filterQ;
+		tripletCR /= filterQ;
+		correctionR = correctionR + tripletFactorR;
+		
+		count++;
+		//finished Triplet section- 15K
+		
+		if (flop)
+		{
+			airFactorAL = airPrevAL - inputSampleL;
+			airFactorAR = airPrevAR - inputSampleR;
+			if (flipA)
+			{
+				airEvenAL += airFactorAL;
+				airOddAL -= airFactorAL;
+				airFactorAL = airEvenAL * airIntensity;
+
+				airEvenAR += airFactorAR;
+				airOddAR -= airFactorAR;
+				airFactorAR = airEvenAR * airIntensity;
+			}
+			else
+			{
+				airOddAL += airFactorAL;
+				airEvenAL -= airFactorAL;
+				airFactorAL = airOddAL * airIntensity;
+
+				airOddAR += airFactorAR;
+				airEvenAR -= airFactorAR;
+				airFactorAR = airOddAR * airIntensity;
+			}
+			airOddAL = (airOddAL - ((airOddAL - airEvenAL)/256.0)) / filterQ;
+			airEvenAL = (airEvenAL - ((airEvenAL - airOddAL)/256.0)) / filterQ;
+			airPrevAL = inputSampleL;
+			correctionL = correctionL + airFactorAL;
+
+			airOddAR = (airOddAR - ((airOddAR - airEvenAR)/256.0)) / filterQ;
+			airEvenAR = (airEvenAR - ((airEvenAR - airOddAR)/256.0)) / filterQ;
+			airPrevAR = inputSampleR;
+			correctionR = correctionR + airFactorAR;
+
+			flipA = !flipA;
+		}
+		else
+		{
+			airFactorBL = airPrevBL - inputSampleL;
+			airFactorBR = airPrevBR - inputSampleR;
+			if (flipB)
+			{
+				airEvenBL += airFactorBL;
+				airOddBL -= airFactorBL;
+				airFactorBL = airEvenBL * airIntensity;
+
+				airEvenBR += airFactorBR;
+				airOddBR -= airFactorBR;
+				airFactorBR = airEvenBR * airIntensity;
+			}
+			else
+			{
+				airOddBL += airFactorBL;
+				airEvenBL -= airFactorBL;
+				airFactorBL = airOddBL * airIntensity;
+
+				airOddBR += airFactorBR;
+				airEvenBR -= airFactorBR;
+				airFactorBR = airOddBR * airIntensity;
+			}
+			airOddBL = (airOddBL - ((airOddBL - airEvenBL)/256.0)) / filterQ;
+			airEvenBL = (airEvenBL - ((airEvenBL - airOddBL)/256.0)) / filterQ;
+			airPrevBL = inputSampleL;
+			correctionL = correctionL + airFactorBL;
+
+			airOddBR = (airOddBR - ((airOddBR - airEvenBR)/256.0)) / filterQ;
+			airEvenBR = (airEvenBR - ((airEvenBR - airOddBR)/256.0)) / filterQ;
+			airPrevBR = inputSampleR;
+			correctionR = correctionR + airFactorBR;
+
+			flipB = !flipB;
+		}
+		//11K one
+		airFactorCL = airPrevCL - inputSampleL;
+		airFactorCR = airPrevCR - inputSampleR;
+		if (flop)
+		{
+			airEvenCL += airFactorCL;
+			airOddCL -= airFactorCL;
+			airFactorCL = airEvenCL * hiIntensity;
+
+			airEvenCR += airFactorCR;
+			airOddCR -= airFactorCR;
+			airFactorCR = airEvenCR * hiIntensity;
+		}
+		else
+		{
+			airOddCL += airFactorCL;
+			airEvenCL -= airFactorCL;
+			airFactorCL = airOddCL * hiIntensity;
+
+			airOddCR += airFactorCR;
+			airEvenCR -= airFactorCR;
+			airFactorCR = airOddCR * hiIntensity;
+		}
+		airOddCL = (airOddCL - ((airOddCL - airEvenCL)/256.0)) / filterQ;
+		airEvenCL = (airEvenCL - ((airEvenCL - airOddCL)/256.0)) / filterQ;
+		airPrevCL = inputSampleL;
+		correctionL = correctionL + airFactorCL;
+
+		airOddCR = (airOddCR - ((airOddCR - airEvenCR)/256.0)) / filterQ;
+		airEvenCR = (airEvenCR - ((airEvenCR - airOddCR)/256.0)) / filterQ;
+		airPrevCR = inputSampleR;
+		correctionR = correctionR + airFactorCR;
+
+		flop = !flop;
+		
+		inputSampleL += correctionL;
+		inputSampleR += correctionR;
+		
+		if (output < 1.0) {
+			inputSampleL *= output;
+			inputSampleR *= output;
+		}
+		if (wet < 1.0) {
+			inputSampleL = (drySampleL*dry)+(inputSampleL*wet);
+			inputSampleR = (drySampleR*dry)+(inputSampleR*wet);
+		}
+		//nice little output stage template: if we have another scale of floating point
+		//number, we really don't want to meaninglessly multiply that by 1.0.
+		
+		//noise shaping to 32-bit floating point
+		if (fpFlip) {
+			fpTemp = inputSampleL;
+			fpNShapeLA = (fpNShapeLA*fpOld)+((inputSampleL-fpTemp)*fpNew);
+			inputSampleL += fpNShapeLA;
+			fpTemp = inputSampleR;
+			fpNShapeRA = (fpNShapeRA*fpOld)+((inputSampleR-fpTemp)*fpNew);
+			inputSampleR += fpNShapeRA;
+		}
+		else {
+			fpTemp = inputSampleL;
+			fpNShapeLB = (fpNShapeLB*fpOld)+((inputSampleL-fpTemp)*fpNew);
+			inputSampleL += fpNShapeLB;
+			fpTemp = inputSampleR;
+			fpNShapeRB = (fpNShapeRB*fpOld)+((inputSampleR-fpTemp)*fpNew);
+			inputSampleR += fpNShapeRB;
+		}
+		fpFlip = !fpFlip;
+		//end noise shaping on 32 bit output
+
+		*out1 = inputSampleL;
+		*out2 = inputSampleR;
+
+		*in1++;
+		*in2++;
+		*out1++;
+		*out2++;
+    }
+}
+
+void Air::processDoubleReplacing(double **inputs, double **outputs, VstInt32 sampleFrames) 
+{
+    double* in1  =  inputs[0];
+    double* in2  =  inputs[1];
+    double* out1 = outputs[0];
+    double* out2 = outputs[1];
+
+	double hiIntensity = -pow(((A*2.0)-1.0),3)*2;
+	double tripletintensity = -pow(((B*2.0)-1.0),3);
+	double airIntensity = -pow(((C*2.0)-1.0),3)/2;
+	double filterQ = 2.1-D;
+	double output = E;
+	double wet = F;
+	double dry = 1.0-wet;
+	
+	double fpTemp; //this is different from singlereplacing
+	long double fpOld = 0.618033988749894848204586; //golden ratio!
+	long double fpNew = 1.0 - fpOld;	
+
+	long double inputSampleL;
+	long double inputSampleR;
+	double drySampleL;
+	double drySampleR;
+	double correctionL;
+	double correctionR;
+
+    while (--sampleFrames >= 0)
+    {
+		inputSampleL = *in1;
+		inputSampleR = *in2;
+		if (inputSampleL<1.2e-38 && -inputSampleL<1.2e-38) {
+			static int noisesource = 0;
+			//this declares a variable before anything else is compiled. It won't keep assigning
+			//it to 0 for every sample, it's as if the declaration doesn't exist in this context,
+			//but it lets me add this denormalization fix in a single place rather than updating
+			//it in three different locations. The variable isn't thread-safe but this is only
+			//a random seed and we can share it with whatever.
+			noisesource = noisesource % 1700021; noisesource++;
+			int residue = noisesource * noisesource;
+			residue = residue % 170003; residue *= residue;
+			residue = residue % 17011; residue *= residue;
+			residue = residue % 1709; residue *= residue;
+			residue = residue % 173; residue *= residue;
+			residue = residue % 17;
+			double applyresidue = residue;
+			applyresidue *= 0.00000001;
+			applyresidue *= 0.00000001;
+			inputSampleL = applyresidue;
+		}
+		if (inputSampleR<1.2e-38 && -inputSampleR<1.2e-38) {
+			static int noisesource = 0;
+			noisesource = noisesource % 1700021; noisesource++;
+			int residue = noisesource * noisesource;
+			residue = residue % 170003; residue *= residue;
+			residue = residue % 17011; residue *= residue;
+			residue = residue % 1709; residue *= residue;
+			residue = residue % 173; residue *= residue;
+			residue = residue % 17;
+			double applyresidue = residue;
+			applyresidue *= 0.00000001;
+			applyresidue *= 0.00000001;
+			inputSampleR = applyresidue;
+			//this denormalization routine produces a white noise at -300 dB which the noise
+			//shaping will interact with to produce a bipolar output, but the noise is actually
+			//all positive. That should stop any variables from going denormal, and the routine
+			//only kicks in if digital black is input. As a final touch, if you save to 24-bit
+			//the silence will return to being digital black again.
+		}
+		drySampleL = inputSampleL;
+		drySampleR = inputSampleR;
+		
+		correctionL = 0.0;
+		correctionR = 0.0;  //from here on down, please add L and R to the code
+		
+		if (count < 1 || count > 3) count = 1;
+		tripletFactorL = tripletPrevL - inputSampleL;
+		tripletFactorR = tripletPrevR - inputSampleR;
+		switch (count)
+		{
+			case 1:
+				tripletAL += tripletFactorL;
+				tripletCL -= tripletFactorL;
+				tripletFactorL = tripletAL * tripletintensity;
+				tripletPrevL = tripletMidL;
+				tripletMidL = inputSampleL;
+				
+				tripletAR += tripletFactorR;
+				tripletCR -= tripletFactorR;
+				tripletFactorR = tripletAR * tripletintensity;
+				tripletPrevR = tripletMidR;
+				tripletMidR = inputSampleR;
+				break;
+			case 2:
+				tripletBL += tripletFactorL;
+				tripletAL -= tripletFactorL;
+				tripletFactorL = tripletBL * tripletintensity;
+				tripletPrevL = tripletMidL;
+				tripletMidL = inputSampleL;
+				
+				tripletBR += tripletFactorR;
+				tripletAR -= tripletFactorR;
+				tripletFactorR = tripletBR * tripletintensity;
+				tripletPrevR = tripletMidR;
+				tripletMidR = inputSampleR;
+				break;
+			case 3:
+				tripletCL += tripletFactorL;
+				tripletBL -= tripletFactorL;
+				tripletFactorL = tripletCL * tripletintensity;
+				tripletPrevL = tripletMidL;
+				tripletMidL = inputSampleL;
+				
+				tripletCR += tripletFactorR;
+				tripletBR -= tripletFactorR;
+				tripletFactorR = tripletCR * tripletintensity;
+				tripletPrevR = tripletMidR;
+				tripletMidR = inputSampleR;
+				break;
+		}
+		tripletAL /= filterQ;
+		tripletBL /= filterQ;
+		tripletCL /= filterQ;
+		correctionL = correctionL + tripletFactorL;
+		
+		tripletAR /= filterQ;
+		tripletBR /= filterQ;
+		tripletCR /= filterQ;
+		correctionR = correctionR + tripletFactorR;
+		
+		count++;
+		//finished Triplet section- 15K
+		
+		if (flop)
+		{
+			airFactorAL = airPrevAL - inputSampleL;
+			airFactorAR = airPrevAR - inputSampleR;
+			if (flipA)
+			{
+				airEvenAL += airFactorAL;
+				airOddAL -= airFactorAL;
+				airFactorAL = airEvenAL * airIntensity;
+				
+				airEvenAR += airFactorAR;
+				airOddAR -= airFactorAR;
+				airFactorAR = airEvenAR * airIntensity;
+			}
+			else
+			{
+				airOddAL += airFactorAL;
+				airEvenAL -= airFactorAL;
+				airFactorAL = airOddAL * airIntensity;
+				
+				airOddAR += airFactorAR;
+				airEvenAR -= airFactorAR;
+				airFactorAR = airOddAR * airIntensity;
+			}
+			airOddAL = (airOddAL - ((airOddAL - airEvenAL)/256.0)) / filterQ;
+			airEvenAL = (airEvenAL - ((airEvenAL - airOddAL)/256.0)) / filterQ;
+			airPrevAL = inputSampleL;
+			correctionL = correctionL + airFactorAL;
+			
+			airOddAR = (airOddAR - ((airOddAR - airEvenAR)/256.0)) / filterQ;
+			airEvenAR = (airEvenAR - ((airEvenAR - airOddAR)/256.0)) / filterQ;
+			airPrevAR = inputSampleR;
+			correctionR = correctionR + airFactorAR;
+			
+			flipA = !flipA;
+		}
+		else
+		{
+			airFactorBL = airPrevBL - inputSampleL;
+			airFactorBR = airPrevBR - inputSampleR;
+			if (flipB)
+			{
+				airEvenBL += airFactorBL;
+				airOddBL -= airFactorBL;
+				airFactorBL = airEvenBL * airIntensity;
+				
+				airEvenBR += airFactorBR;
+				airOddBR -= airFactorBR;
+				airFactorBR = airEvenBR * airIntensity;
+			}
+			else
+			{
+				airOddBL += airFactorBL;
+				airEvenBL -= airFactorBL;
+				airFactorBL = airOddBL * airIntensity;
+				
+				airOddBR += airFactorBR;
+				airEvenBR -= airFactorBR;
+				airFactorBR = airOddBR * airIntensity;
+			}
+			airOddBL = (airOddBL - ((airOddBL - airEvenBL)/256.0)) / filterQ;
+			airEvenBL = (airEvenBL - ((airEvenBL - airOddBL)/256.0)) / filterQ;
+			airPrevBL = inputSampleL;
+			correctionL = correctionL + airFactorBL;
+			
+			airOddBR = (airOddBR - ((airOddBR - airEvenBR)/256.0)) / filterQ;
+			airEvenBR = (airEvenBR - ((airEvenBR - airOddBR)/256.0)) / filterQ;
+			airPrevBR = inputSampleR;
+			correctionR = correctionR + airFactorBR;
+			
+			flipB = !flipB;
+		}
+		//11K one
+		airFactorCL = airPrevCL - inputSampleL;
+		airFactorCR = airPrevCR - inputSampleR;
+		if (flop)
+		{
+			airEvenCL += airFactorCL;
+			airOddCL -= airFactorCL;
+			airFactorCL = airEvenCL * hiIntensity;
+			
+			airEvenCR += airFactorCR;
+			airOddCR -= airFactorCR;
+			airFactorCR = airEvenCR * hiIntensity;
+		}
+		else
+		{
+			airOddCL += airFactorCL;
+			airEvenCL -= airFactorCL;
+			airFactorCL = airOddCL * hiIntensity;
+			
+			airOddCR += airFactorCR;
+			airEvenCR -= airFactorCR;
+			airFactorCR = airOddCR * hiIntensity;
+		}
+		airOddCL = (airOddCL - ((airOddCL - airEvenCL)/256.0)) / filterQ;
+		airEvenCL = (airEvenCL - ((airEvenCL - airOddCL)/256.0)) / filterQ;
+		airPrevCL = inputSampleL;
+		correctionL = correctionL + airFactorCL;
+		
+		airOddCR = (airOddCR - ((airOddCR - airEvenCR)/256.0)) / filterQ;
+		airEvenCR = (airEvenCR - ((airEvenCR - airOddCR)/256.0)) / filterQ;
+		airPrevCR = inputSampleR;
+		correctionR = correctionR + airFactorCR;
+		
+		flop = !flop;
+		
+		inputSampleL += correctionL;
+		inputSampleR += correctionR;
+		
+		if (output < 1.0) {
+			inputSampleL *= output;
+			inputSampleR *= output;
+		}
+		if (wet < 1.0) {
+			inputSampleL = (drySampleL*dry)+(inputSampleL*wet);
+			inputSampleR = (drySampleR*dry)+(inputSampleR*wet);
+		}
+		//nice little output stage template: if we have another scale of floating point
+		//number, we really don't want to meaninglessly multiply that by 1.0.
+
+		//noise shaping to 64-bit floating point
+		if (fpFlip) {
+			fpTemp = inputSampleL;
+			fpNShapeLA = (fpNShapeLA*fpOld)+((inputSampleL-fpTemp)*fpNew);
+			inputSampleL += fpNShapeLA;
+			fpTemp = inputSampleR;
+			fpNShapeRA = (fpNShapeRA*fpOld)+((inputSampleR-fpTemp)*fpNew);
+			inputSampleR += fpNShapeRA;
+		}
+		else {
+			fpTemp = inputSampleL;
+			fpNShapeLB = (fpNShapeLB*fpOld)+((inputSampleL-fpTemp)*fpNew);
+			inputSampleL += fpNShapeLB;
+			fpTemp = inputSampleR;
+			fpNShapeRB = (fpNShapeRB*fpOld)+((inputSampleR-fpTemp)*fpNew);
+			inputSampleR += fpNShapeRB;
+		}
+		fpFlip = !fpFlip;
+		//end noise shaping on 64 bit output
+
+		*out1 = inputSampleL;
+		*out2 = inputSampleR;
+
+		*in1++;
+		*in2++;
+		*out1++;
+		*out2++;
+    }
+}
\ No newline at end of file