DubCenter

4 files changed, 1050 insertions, 0 deletions

diff --git a/plugins/LinuxVST/CMakeLists.txt b/plugins/LinuxVST/CMakeLists.txt
index 1a450b7..9708555 100755
--- a/plugins/LinuxVST/CMakeLists.txt
+++ b/plugins/LinuxVST/CMakeLists.txt
@@ -27,6 +27,7 @@ add_airwindows_plugin(Distance2)
 add_airwindows_plugin(Ditherbox)
 add_airwindows_plugin(DoublePaul)
 add_airwindows_plugin(DrumSlam)
+add_airwindows_plugin(DubCenter)
 add_airwindows_plugin(DubSub)
 add_airwindows_plugin(ElectroHat)
 add_airwindows_plugin(Energy)
diff --git a/plugins/LinuxVST/src/DubCenter/DubCenter.cpp b/plugins/LinuxVST/src/DubCenter/DubCenter.cpp
new file mode 100755
index 0000000..098e723
--- /dev/null
+++ b/plugins/LinuxVST/src/DubCenter/DubCenter.cpp
@@ -0,0 +1,249 @@
+/* ========================================
+ *  DubCenter - DubCenter.h
+ *  Copyright (c) 2016 airwindows, All rights reserved
+ * ======================================== */
+
+#ifndef __DubCenter_H
+#include "DubCenter.h"
+#endif
+
+AudioEffect* createEffectInstance(audioMasterCallback audioMaster) {return new DubCenter(audioMaster);}
+
+DubCenter::DubCenter(audioMasterCallback audioMaster) :
+    AudioEffectX(audioMaster, kNumPrograms, kNumParameters)
+{
+	A = 0.9;
+	B = 0.5;
+	C = 0.74;
+	D = 1.0;
+	E = 0.95;
+	F = 0.5;
+	G = 0.2;
+	H = 0.2;
+	I = 0.5;
+	J = 1.0;
+
+	WasNegative = false;
+	SubOctave = false;
+	flip = false;
+	bflip = 0;
+	iirDriveSampleAL = 0.0;
+	iirDriveSampleBL = 0.0;
+	iirDriveSampleCL = 0.0;
+	iirDriveSampleDL = 0.0;
+	iirDriveSampleEL = 0.0;
+	iirDriveSampleFL = 0.0;
+	iirDriveSampleAR = 0.0;
+	iirDriveSampleBR = 0.0;
+	iirDriveSampleCR = 0.0;
+	iirDriveSampleDR = 0.0;
+	iirDriveSampleER = 0.0;
+	iirDriveSampleFR = 0.0;
+	
+	iirHeadBumpA = 0.0;
+	iirHeadBumpB = 0.0;
+	iirHeadBumpC = 0.0;
+	
+	iirSubBumpA = 0.0;
+	iirSubBumpB = 0.0;
+	iirSubBumpC = 0.0;
+	
+	lastHeadBump = 0.0;
+	lastSubBump = 0.0;
+	
+	iirSampleA = 0.0;
+	iirSampleB = 0.0;
+	iirSampleC = 0.0;
+	iirSampleD = 0.0;
+	iirSampleE = 0.0;
+	iirSampleF = 0.0;
+	iirSampleG = 0.0;
+	iirSampleH = 0.0;
+	iirSampleI = 0.0;
+	iirSampleJ = 0.0;
+	iirSampleK = 0.0;
+	iirSampleL = 0.0;
+	iirSampleM = 0.0;
+	iirSampleN = 0.0;
+	iirSampleO = 0.0;
+	iirSampleP = 0.0;
+	iirSampleQ = 0.0;
+	iirSampleR = 0.0;
+	iirSampleS = 0.0;
+	iirSampleT = 0.0;
+	iirSampleU = 0.0;
+	iirSampleV = 0.0;
+	iirSampleW = 0.0;
+	iirSampleX = 0.0;
+	iirSampleY = 0.0;
+	iirSampleZ = 0.0;	
+	
+	oscGate = 1.0;
+	
+	fpNShapeL = 0.0;
+	fpNShapeR = 0.0;
+	//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
+}
+
+DubCenter::~DubCenter() {}
+VstInt32 DubCenter::getVendorVersion () {return 1000;}
+void DubCenter::setProgramName(char *name) {vst_strncpy (_programName, name, kVstMaxProgNameLen);}
+void DubCenter::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 DubCenter::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;
+	chunkData[6] = G;
+	chunkData[7] = H;
+	chunkData[8] = I;
+	chunkData[9] = J;
+	/* 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 DubCenter::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]);
+	G = pinParameter(chunkData[6]);
+	H = pinParameter(chunkData[7]);
+	I = pinParameter(chunkData[8]);
+	J = pinParameter(chunkData[9]);
+	/* 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 DubCenter::setParameter(VstInt32 index, float value) {
+    switch (index) {
+        case kParamA: A = value; break;
+        case kParamB: B = value; break;
+        case kParamC: C = value; break;
+        case kParamD: D = value; break;
+        case kParamE: E = value; break;
+        case kParamF: F = value; break;
+        case kParamG: G = value; break;
+        case kParamH: H = value; break;
+        case kParamI: I = value; break;
+        case kParamJ: J = value; break;
+        default: throw; // unknown parameter, shouldn't happen!
+    }
+}
+
+float DubCenter::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;
+        case kParamG: return G; break;
+        case kParamH: return H; break;
+        case kParamI: return I; break;
+        case kParamJ: return J; 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 DubCenter::getParameterName(VstInt32 index, char *text) {
+    switch (index) {
+        case kParamA: vst_strncpy (text, "TGrind", kVstMaxParamStrLen); break;
+		case kParamB: vst_strncpy (text, "Grd/Out", kVstMaxParamStrLen); break;
+		case kParamC: vst_strncpy (text, "XOver", kVstMaxParamStrLen); break;
+		case kParamD: vst_strncpy (text, "BsDrive", kVstMaxParamStrLen); break;
+		case kParamE: vst_strncpy (text, "BsVoice", kVstMaxParamStrLen); break;
+		case kParamF: vst_strncpy (text, "BassOut", kVstMaxParamStrLen); break;
+		case kParamG: vst_strncpy (text, "SbDrive", kVstMaxParamStrLen); break;
+		case kParamH: vst_strncpy (text, "SbVoice", kVstMaxParamStrLen); break;
+		case kParamI: vst_strncpy (text, "SubOut", kVstMaxParamStrLen); break;
+		case kParamJ: 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 DubCenter::getParameterDisplay(VstInt32 index, char *text) {
+    switch (index) {
+        case kParamA: float2string (A, text, kVstMaxParamStrLen); break;
+        case kParamB: float2string ((B*2.0)-1.0, text, kVstMaxParamStrLen); break;
+        case kParamC: float2string (C, text, kVstMaxParamStrLen); break;
+        case kParamD: float2string (D, text, kVstMaxParamStrLen); break;
+        case kParamE: float2string (E, text, kVstMaxParamStrLen); break;
+        case kParamF: float2string ((F*2.0)-1.0, text, kVstMaxParamStrLen); break;
+        case kParamG: float2string (G, text, kVstMaxParamStrLen); break;
+        case kParamH: float2string (H, text, kVstMaxParamStrLen); break;
+        case kParamI: float2string ((I*2.0)-1.0, text, kVstMaxParamStrLen); break;
+        case kParamJ: float2string (J, text, kVstMaxParamStrLen); break;
+        default: break; // unknown parameter, shouldn't happen!
+	} //this displays the values and handles 'popups' where it's discrete choices
+}
+
+void DubCenter::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;
+        case kParamG: vst_strncpy (text, "", kVstMaxParamStrLen); break;
+        case kParamH: vst_strncpy (text, "", kVstMaxParamStrLen); break;
+        case kParamI: vst_strncpy (text, "", kVstMaxParamStrLen); break;
+        case kParamJ: vst_strncpy (text, "", kVstMaxParamStrLen); break;
+		default: break; // unknown parameter, shouldn't happen!
+    }
+}
+
+VstInt32 DubCenter::canDo(char *text) 
+{ return (_canDo.find(text) == _canDo.end()) ? -1: 1; } // 1 = yes, -1 = no, 0 = don't know
+
+bool DubCenter::getEffectName(char* name) {
+    vst_strncpy(name, "DubCenter", kVstMaxProductStrLen); return true;
+}
+
+VstPlugCategory DubCenter::getPlugCategory() {return kPlugCategEffect;}
+
+bool DubCenter::getProductString(char* text) {
+  	vst_strncpy (text, "airwindows DubCenter", kVstMaxProductStrLen); return true;
+}
+
+bool DubCenter::getVendorString(char* text) {
+  	vst_strncpy (text, "airwindows", kVstMaxVendorStrLen); return true;
+}
diff --git a/plugins/LinuxVST/src/DubCenter/DubCenter.h b/plugins/LinuxVST/src/DubCenter/DubCenter.h
new file mode 100755
index 0000000..9d7c7e6
--- /dev/null
+++ b/plugins/LinuxVST/src/DubCenter/DubCenter.h
@@ -0,0 +1,137 @@
+/* ========================================
+ *  DubCenter - DubCenter.h
+ *  Created 8/12/11 by SPIAdmin 
+ *  Copyright (c) 2011 __MyCompanyName__, All rights reserved
+ * ======================================== */
+
+#ifndef __DubCenter_H
+#define __DubCenter_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,
+	kParamG = 6,
+	kParamH = 7,
+	kParamI = 8,
+	kParamJ = 9,
+	kNumParameters = 10
+}; //
+
+const int kNumPrograms = 0;
+const int kNumInputs = 2;
+const int kNumOutputs = 2;
+const unsigned long kUniqueId = 'dbct';    //Change this to what the AU identity is!
+
+class DubCenter : 
+    public AudioEffectX 
+{
+public:
+    DubCenter(audioMasterCallback audioMaster);
+    ~DubCenter();
+    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 iirDriveSampleAL;
+	double iirDriveSampleBL;
+	double iirDriveSampleCL;
+	double iirDriveSampleDL;
+	double iirDriveSampleEL;
+	double iirDriveSampleFL;
+	double iirDriveSampleAR;
+	double iirDriveSampleBR;
+	double iirDriveSampleCR;
+	double iirDriveSampleDR;
+	double iirDriveSampleER;
+	double iirDriveSampleFR;
+	bool flip; //drive things
+	
+	int bflip;
+	bool WasNegative;
+	bool SubOctave;
+	double iirHeadBumpA;
+	double iirHeadBumpB;
+	double iirHeadBumpC;
+	
+	double iirSubBumpA;
+	double iirSubBumpB;
+	double iirSubBumpC;
+	
+	double lastHeadBump;
+	double lastSubBump;
+	
+	double iirSampleA;
+	double iirSampleB;
+	double iirSampleC;
+	double iirSampleD;
+	double iirSampleE;
+	double iirSampleF;
+	double iirSampleG;
+	double iirSampleH;
+	double iirSampleI;
+	double iirSampleJ;
+	double iirSampleK;
+	double iirSampleL;
+	double iirSampleM;
+	double iirSampleN;
+	double iirSampleO;
+	double iirSampleP;
+	double iirSampleQ;
+	double iirSampleR;
+	double iirSampleS;
+	double iirSampleT;
+	double iirSampleU;
+	double iirSampleV;
+	double iirSampleW;
+	double iirSampleX;
+	double iirSampleY;
+	double iirSampleZ;		
+	double oscGate;
+	
+	long double fpNShapeL;
+	long double fpNShapeR;
+	//default stuff
+	
+    float A;
+    float B;
+    float C;
+    float D;
+    float E;
+    float F;
+    float G;
+    float H;
+    float I;
+    float J;
+	
+};
+
+#endif
diff --git a/plugins/LinuxVST/src/DubCenter/DubCenterProc.cpp b/plugins/LinuxVST/src/DubCenter/DubCenterProc.cpp
new file mode 100755
index 0000000..75c162b
--- /dev/null
+++ b/plugins/LinuxVST/src/DubCenter/DubCenterProc.cpp
@@ -0,0 +1,663 @@
+/* ========================================
+ *  DubCenter - DubCenter.h
+ *  Copyright (c) 2016 airwindows, All rights reserved
+ * ======================================== */
+
+#ifndef __DubCenter_H
+#include "DubCenter.h"
+#endif
+
+void DubCenter::processReplacing(float **inputs, float **outputs, VstInt32 sampleFrames) 
+{
+    float* in1  =  inputs[0];
+    float* in2  =  inputs[1];
+    float* out1 = outputs[0];
+    float* out2 = outputs[1];
+
+	double overallscale = 1.0;
+	overallscale /= 44100.0;
+	overallscale *= getSampleRate();
+
+	double driveone = pow(A*3.0,2);
+	double driveoutput = (B*2.0)-1.0;
+	double iirAmount = ((C*0.33)+0.1)/overallscale;
+	double ataLowpassL;
+	double ataLowpassR;
+	double ataLowpass;
+	double randy;
+	double invrandy;
+	double HeadBump = 0.0;
+	double BassGain = D * 0.1;
+	double HeadBumpFreq = ((E*0.1)+0.0001)/overallscale;
+	double iirBmount = HeadBumpFreq/44.1;
+	double altBmount = 1.0 - iirBmount;
+	double BassOutGain = (F*2.0)-1.0;
+	double SubBump = 0.0;
+	double SubGain = G * 0.1;
+	double SubBumpFreq = ((H*0.1)+0.0001)/overallscale;
+	double iirCmount = SubBumpFreq/44.1;
+	double altCmount = 1.0 - iirCmount;
+	double SubOutGain = (I*2.0)-1.0;
+	double clamp = 0.0;
+	double out;
+	double fuzz = 0.111;
+	double wet = J;
+	double dry = 1.0-wet;
+	double glitch = 0.60;
+	double tempSample;
+	
+    
+    while (--sampleFrames >= 0)
+    {
+		long double inputSampleL = *in1;
+		long double inputSampleR = *in2;
+
+		static int noisesourceL = 0;
+		static int noisesourceR = 850010;
+		int residue;
+		double applyresidue;
+		
+		noisesourceL = noisesourceL % 1700021; noisesourceL++;
+		residue = noisesourceL * noisesourceL;
+		residue = residue % 170003; residue *= residue;
+		residue = residue % 17011; residue *= residue;
+		residue = residue % 1709; residue *= residue;
+		residue = residue % 173; residue *= residue;
+		residue = residue % 17;
+		applyresidue = residue;
+		applyresidue *= 0.00000001;
+		applyresidue *= 0.00000001;
+		inputSampleL += applyresidue;
+		if (inputSampleL<1.2e-38 && -inputSampleL<1.2e-38) {
+			inputSampleL -= applyresidue;
+		}
+		
+		noisesourceR = noisesourceR % 1700021; noisesourceR++;
+		residue = noisesourceR * noisesourceR;
+		residue = residue % 170003; residue *= residue;
+		residue = residue % 17011; residue *= residue;
+		residue = residue % 1709; residue *= residue;
+		residue = residue % 173; residue *= residue;
+		residue = residue % 17;
+		applyresidue = residue;
+		applyresidue *= 0.00000001;
+		applyresidue *= 0.00000001;
+		inputSampleR += applyresidue;
+		if (inputSampleR<1.2e-38 && -inputSampleR<1.2e-38) {
+			inputSampleR -= applyresidue;
+		}
+		//for live air, we always apply the dither noise. Then, if our result is 
+		//effectively digital black, we'll subtract it again. We want a 'air' hiss
+		long double drySampleL = inputSampleL;
+		long double drySampleR = inputSampleR;
+		
+		// here's the plan.
+		// Grind Boost
+		// Grind Output Level
+		// Bass Split Freq
+		// Bass Drive
+		// Bass Voicing
+		// Bass Output Level
+		// Sub Oct Drive
+		// Sub Voicing
+		// Sub Output Level
+		// Dry/Wet
+		
+		oscGate += fabs((inputSampleL + inputSampleR) * 5.0);
+		oscGate -= 0.001;
+		if (oscGate > 1.0) oscGate = 1.0;
+		if (oscGate < 0) oscGate = 0;
+		//got a value that only goes down low when there's silence or near silence on input
+		clamp = 1.0-oscGate;
+		clamp *= 0.00001;
+		//set up the thing to choke off oscillations- belt and suspenders affair
+		
+		
+		if (flip)
+		{
+			tempSample = inputSampleL;
+			iirDriveSampleAL = (iirDriveSampleAL * (1 - iirAmount)) + (inputSampleL * iirAmount);
+			inputSampleL -= iirDriveSampleAL;
+			iirDriveSampleCL = (iirDriveSampleCL * (1 - iirAmount)) + (inputSampleL * iirAmount);
+			inputSampleL -= iirDriveSampleCL;
+			iirDriveSampleEL = (iirDriveSampleEL * (1 - iirAmount)) + (inputSampleL * iirAmount);
+			inputSampleL -= iirDriveSampleEL;
+			ataLowpassL = tempSample - inputSampleL;
+			
+			tempSample = inputSampleR;
+			iirDriveSampleAR = (iirDriveSampleAR * (1 - iirAmount)) + (inputSampleR * iirAmount);
+			inputSampleR -= iirDriveSampleAR;
+			iirDriveSampleCR = (iirDriveSampleCR * (1 - iirAmount)) + (inputSampleR * iirAmount);
+			inputSampleR -= iirDriveSampleCR;
+			iirDriveSampleER = (iirDriveSampleER * (1 - iirAmount)) + (inputSampleR * iirAmount);
+			inputSampleR -= iirDriveSampleER;
+			ataLowpassR = tempSample - inputSampleR;
+		}
+		else
+		{
+			tempSample = inputSampleL;
+			iirDriveSampleBL = (iirDriveSampleBL * (1 - iirAmount)) + (inputSampleL * iirAmount);
+			inputSampleL -= iirDriveSampleBL;
+			iirDriveSampleDL = (iirDriveSampleDL * (1 - iirAmount)) + (inputSampleL * iirAmount);
+			inputSampleL -= iirDriveSampleDL;
+			iirDriveSampleFL = (iirDriveSampleFL * (1 - iirAmount)) + (inputSampleL * iirAmount);
+			inputSampleL -= iirDriveSampleFL;
+			ataLowpassL = tempSample - inputSampleL;
+			
+			tempSample = inputSampleR;
+			iirDriveSampleBR = (iirDriveSampleBR * (1 - iirAmount)) + (inputSampleR * iirAmount);
+			inputSampleR -= iirDriveSampleBR;
+			iirDriveSampleDR = (iirDriveSampleDR * (1 - iirAmount)) + (inputSampleR * iirAmount);
+			inputSampleR -= iirDriveSampleDR;
+			iirDriveSampleFR = (iirDriveSampleFR * (1 - iirAmount)) + (inputSampleR * iirAmount);
+			inputSampleR -= iirDriveSampleFR;
+			ataLowpassR = tempSample - inputSampleR;
+		}
+		//highpass section
+		
+		if (inputSampleL > 1.0) {inputSampleL = 1.0;}
+		if (inputSampleL < -1.0) {inputSampleL = -1.0;}
+		if (inputSampleR > 1.0) {inputSampleR = 1.0;}
+		if (inputSampleR < -1.0) {inputSampleR = -1.0;}
+		
+		out = driveone;
+		while (out > glitch)
+		{
+			out -= glitch;
+			inputSampleL -= (inputSampleL * (fabs(inputSampleL) * glitch) * (fabs(inputSampleL) * glitch) );
+			inputSampleL *= (1.0+glitch);
+		}
+		//that's taken care of the really high gain stuff
+		inputSampleL -= (inputSampleL * (fabs(inputSampleL) * out) * (fabs(inputSampleL) * out) );
+		inputSampleL *= (1.0+out);
+		
+		
+		out = driveone;
+		while (out > glitch)
+		{
+			out -= glitch;
+			inputSampleR -= (inputSampleR * (fabs(inputSampleR) * glitch) * (fabs(inputSampleR) * glitch) );
+			inputSampleR *= (1.0+glitch);
+		}
+		//that's taken care of the really high gain stuff
+		inputSampleR -= (inputSampleR * (fabs(inputSampleR) * out) * (fabs(inputSampleR) * out) );
+		inputSampleR *= (1.0+out);
+		
+		ataLowpass = (ataLowpassL + ataLowpassR) / 2.0;
+		
+		if (ataLowpass > 0)
+		{if (WasNegative){SubOctave = !SubOctave;} WasNegative = false;}
+		else {WasNegative = true;}
+		//set up polarities for sub-bass version
+		randy = (rand()/(double)RAND_MAX)*fuzz; //0 to 1 the noise, may not be needed
+		invrandy = (1.0-randy);
+		randy /= 2.0;
+		//set up the noise
+		
+		iirSampleA = (iirSampleA * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleA;
+		iirSampleB = (iirSampleB * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleB;
+		iirSampleC = (iirSampleC * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleC;
+		iirSampleD = (iirSampleD * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleD;
+		iirSampleE = (iirSampleE * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleE;
+		iirSampleF = (iirSampleF * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleF;
+		iirSampleG = (iirSampleG * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleG;
+		iirSampleH = (iirSampleH * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleH;
+		iirSampleI = (iirSampleI * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleI;
+		iirSampleJ = (iirSampleJ * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleJ;
+		iirSampleK = (iirSampleK * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleK;
+		iirSampleL = (iirSampleL * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleL;
+		iirSampleM = (iirSampleM * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleM;
+		iirSampleN = (iirSampleN * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleN;
+		iirSampleO = (iirSampleO * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleO;
+		iirSampleP = (iirSampleP * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleP;
+		iirSampleQ = (iirSampleQ * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleQ;
+		iirSampleR = (iirSampleR * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleR;
+		iirSampleS = (iirSampleS * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleS;
+		iirSampleT = (iirSampleT * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleT;
+		iirSampleU = (iirSampleU * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleU;
+		iirSampleV = (iirSampleV * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleV;
+		
+		switch (bflip)
+		{
+			case 1:				
+				iirHeadBumpA += (ataLowpass * BassGain);
+				iirHeadBumpA -= (iirHeadBumpA * iirHeadBumpA * iirHeadBumpA * HeadBumpFreq);
+				iirHeadBumpA = (invrandy * iirHeadBumpA) + (randy * iirHeadBumpB) + (randy * iirHeadBumpC);
+				if (iirHeadBumpA > 0) iirHeadBumpA -= clamp;
+				if (iirHeadBumpA < 0) iirHeadBumpA += clamp;
+				HeadBump = iirHeadBumpA;
+				break;
+			case 2:
+				iirHeadBumpB += (ataLowpass * BassGain);
+				iirHeadBumpB -= (iirHeadBumpB * iirHeadBumpB * iirHeadBumpB * HeadBumpFreq);
+				iirHeadBumpB = (randy * iirHeadBumpA) + (invrandy * iirHeadBumpB) + (randy * iirHeadBumpC);
+				if (iirHeadBumpB > 0) iirHeadBumpB -= clamp;
+				if (iirHeadBumpB < 0) iirHeadBumpB += clamp;
+				HeadBump = iirHeadBumpB;
+				break;
+			case 3:
+				iirHeadBumpC += (ataLowpass * BassGain);
+				iirHeadBumpC -= (iirHeadBumpC * iirHeadBumpC * iirHeadBumpC * HeadBumpFreq);
+				iirHeadBumpC = (randy * iirHeadBumpA) + (randy * iirHeadBumpB) + (invrandy * iirHeadBumpC);
+				if (iirHeadBumpC > 0) iirHeadBumpC -= clamp;
+				if (iirHeadBumpC < 0) iirHeadBumpC += clamp;
+				HeadBump = iirHeadBumpC;
+				break;
+		}
+		
+		iirSampleW = (iirSampleW * altBmount) + (HeadBump * iirBmount); HeadBump -= iirSampleW;
+		iirSampleX = (iirSampleX * altBmount) + (HeadBump * iirBmount); HeadBump -= iirSampleX;
+		
+		SubBump = HeadBump;
+		iirSampleY = (iirSampleY * altCmount) + (SubBump * iirCmount); SubBump -= iirSampleY;
+		
+		SubBump = fabs(SubBump);
+		if (SubOctave == false) {SubBump = -SubBump;}
+		
+		switch (bflip)
+		{
+			case 1:				
+				iirSubBumpA += (SubBump * SubGain);
+				iirSubBumpA -= (iirSubBumpA * iirSubBumpA * iirSubBumpA * SubBumpFreq);
+				iirSubBumpA = (invrandy * iirSubBumpA) + (randy * iirSubBumpB) + (randy * iirSubBumpC);
+				if (iirSubBumpA > 0) iirSubBumpA -= clamp;
+				if (iirSubBumpA < 0) iirSubBumpA += clamp;
+				SubBump = iirSubBumpA;
+				break;
+			case 2:
+				iirSubBumpB += (SubBump * SubGain);
+				iirSubBumpB -= (iirSubBumpB * iirSubBumpB * iirSubBumpB * SubBumpFreq);
+				iirSubBumpB = (randy * iirSubBumpA) + (invrandy * iirSubBumpB) + (randy * iirSubBumpC);
+				if (iirSubBumpB > 0) iirSubBumpB -= clamp;
+				if (iirSubBumpB < 0) iirSubBumpB += clamp;
+				SubBump = iirSubBumpB;
+				break;
+			case 3:
+				iirSubBumpC += (SubBump * SubGain);
+				iirSubBumpC -= (iirSubBumpC * iirSubBumpC * iirSubBumpC * SubBumpFreq);
+				iirSubBumpC = (randy * iirSubBumpA) + (randy * iirSubBumpB) + (invrandy * iirSubBumpC);
+				if (iirSubBumpC > 0) iirSubBumpC -= clamp;
+				if (iirSubBumpC < 0) iirSubBumpC += clamp;
+				SubBump = iirSubBumpC;
+				break;
+		}
+		
+		iirSampleZ = (iirSampleZ * altCmount) + (SubBump * iirCmount); SubBump -= iirSampleZ;
+		
+		inputSampleL *= driveoutput; //start with the drive section then add lows and subs
+		inputSampleR *= driveoutput; //start with the drive section then add lows and subs
+		
+		inputSampleL += ((HeadBump + lastHeadBump) * BassOutGain);
+		inputSampleL += ((SubBump + lastSubBump) * SubOutGain);
+		
+		inputSampleR += ((HeadBump + lastHeadBump) * BassOutGain);
+		inputSampleR += ((SubBump + lastSubBump) * SubOutGain);
+		
+		lastHeadBump = HeadBump;
+		lastSubBump = SubBump;
+		
+		if (wet != 1.0) {
+			inputSampleL = (inputSampleL * wet) + (drySampleL * dry);
+			inputSampleR = (inputSampleR * wet) + (drySampleR * dry);
+		}
+		//Dry/Wet control, defaults to the last slider
+		
+		flip = !flip;
+		bflip++;
+		if (bflip < 1 || bflip > 3) bflip = 1;
+		
+		//noise shaping to 32-bit floating point
+		float fpTemp = inputSampleL;
+		fpNShapeL += (inputSampleL-fpTemp);
+		inputSampleL += fpNShapeL;
+		//if this confuses you look at the wordlength for fpTemp :)
+		fpTemp = inputSampleR;
+		fpNShapeR += (inputSampleR-fpTemp);
+		inputSampleR += fpNShapeR;
+		//for deeper space and warmth, we try a non-oscillating noise shaping
+		//that is kind of ruthless: it will forever retain the rounding errors
+		//except we'll dial it back a hair at the end of every buffer processed
+		//end noise shaping on 32 bit output
+		
+		*out1 = inputSampleL;
+		*out2 = inputSampleR;
+
+		*in1++;
+		*in2++;
+		*out1++;
+		*out2++;
+    }
+	fpNShapeL *= 0.999999;
+	fpNShapeR *= 0.999999;
+	//we will just delicately dial back the FP noise shaping, not even every sample
+	//this is a good place to put subtle 'no runaway' calculations, though bear in mind
+	//that it will be called more often when you use shorter sample buffers in the DAW.
+	//So, very low latency operation will call these calculations more often.	
+}
+
+void DubCenter::processDoubleReplacing(double **inputs, double **outputs, VstInt32 sampleFrames) 
+{
+    double* in1  =  inputs[0];
+    double* in2  =  inputs[1];
+    double* out1 = outputs[0];
+    double* out2 = outputs[1];
+
+	double overallscale = 1.0;
+	overallscale /= 44100.0;
+	overallscale *= getSampleRate();
+	
+	double driveone = pow(A*3.0,2);
+	double driveoutput = (B*2.0)-1.0;
+	double iirAmount = ((C*0.33)+0.1)/overallscale;
+	double ataLowpassL;
+	double ataLowpassR;
+	double ataLowpass;
+	double randy;
+	double invrandy;
+	double HeadBump = 0.0;
+	double BassGain = D * 0.1;
+	double HeadBumpFreq = ((E*0.1)+0.0001)/overallscale;
+	double iirBmount = HeadBumpFreq/44.1;
+	double altBmount = 1.0 - iirBmount;
+	double BassOutGain = (F*2.0)-1.0;
+	double SubBump = 0.0;
+	double SubGain = G * 0.1;
+	double SubBumpFreq = ((H*0.1)+0.0001)/overallscale;
+	double iirCmount = SubBumpFreq/44.1;
+	double altCmount = 1.0 - iirCmount;
+	double SubOutGain = (I*2.0)-1.0;
+	double clamp = 0.0;
+	double out;
+	double fuzz = 0.111;
+	double wet = J;
+	double dry = 1.0-wet;
+	double glitch = 0.60;
+	double tempSample;
+	
+	while (--sampleFrames >= 0)
+    {
+		long double inputSampleL = *in1;
+		long double inputSampleR = *in2;
+
+		static int noisesourceL = 0;
+		static int noisesourceR = 850010;
+		int residue;
+		double applyresidue;
+		
+		noisesourceL = noisesourceL % 1700021; noisesourceL++;
+		residue = noisesourceL * noisesourceL;
+		residue = residue % 170003; residue *= residue;
+		residue = residue % 17011; residue *= residue;
+		residue = residue % 1709; residue *= residue;
+		residue = residue % 173; residue *= residue;
+		residue = residue % 17;
+		applyresidue = residue;
+		applyresidue *= 0.00000001;
+		applyresidue *= 0.00000001;
+		inputSampleL += applyresidue;
+		if (inputSampleL<1.2e-38 && -inputSampleL<1.2e-38) {
+			inputSampleL -= applyresidue;
+		}
+		
+		noisesourceR = noisesourceR % 1700021; noisesourceR++;
+		residue = noisesourceR * noisesourceR;
+		residue = residue % 170003; residue *= residue;
+		residue = residue % 17011; residue *= residue;
+		residue = residue % 1709; residue *= residue;
+		residue = residue % 173; residue *= residue;
+		residue = residue % 17;
+		applyresidue = residue;
+		applyresidue *= 0.00000001;
+		applyresidue *= 0.00000001;
+		inputSampleR += applyresidue;
+		if (inputSampleR<1.2e-38 && -inputSampleR<1.2e-38) {
+			inputSampleR -= applyresidue;
+		}
+		//for live air, we always apply the dither noise. Then, if our result is 
+		//effectively digital black, we'll subtract it again. We want a 'air' hiss
+		long double drySampleL = inputSampleL;
+		long double drySampleR = inputSampleR;
+		
+		// here's the plan.
+		// Grind Boost
+		// Grind Output Level
+		// Bass Split Freq
+		// Bass Drive
+		// Bass Voicing
+		// Bass Output Level
+		// Sub Oct Drive
+		// Sub Voicing
+		// Sub Output Level
+		// Dry/Wet
+		
+		oscGate += fabs((inputSampleL + inputSampleR) * 5.0);
+		oscGate -= 0.001;
+		if (oscGate > 1.0) oscGate = 1.0;
+		if (oscGate < 0) oscGate = 0;
+		//got a value that only goes down low when there's silence or near silence on input
+		clamp = 1.0-oscGate;
+		clamp *= 0.00001;
+		//set up the thing to choke off oscillations- belt and suspenders affair
+		
+		
+		if (flip)
+		{
+			tempSample = inputSampleL;
+			iirDriveSampleAL = (iirDriveSampleAL * (1 - iirAmount)) + (inputSampleL * iirAmount);
+			inputSampleL -= iirDriveSampleAL;
+			iirDriveSampleCL = (iirDriveSampleCL * (1 - iirAmount)) + (inputSampleL * iirAmount);
+			inputSampleL -= iirDriveSampleCL;
+			iirDriveSampleEL = (iirDriveSampleEL * (1 - iirAmount)) + (inputSampleL * iirAmount);
+			inputSampleL -= iirDriveSampleEL;
+			ataLowpassL = tempSample - inputSampleL;
+			
+			tempSample = inputSampleR;
+			iirDriveSampleAR = (iirDriveSampleAR * (1 - iirAmount)) + (inputSampleR * iirAmount);
+			inputSampleR -= iirDriveSampleAR;
+			iirDriveSampleCR = (iirDriveSampleCR * (1 - iirAmount)) + (inputSampleR * iirAmount);
+			inputSampleR -= iirDriveSampleCR;
+			iirDriveSampleER = (iirDriveSampleER * (1 - iirAmount)) + (inputSampleR * iirAmount);
+			inputSampleR -= iirDriveSampleER;
+			ataLowpassR = tempSample - inputSampleR;
+		}
+		else
+		{
+			tempSample = inputSampleL;
+			iirDriveSampleBL = (iirDriveSampleBL * (1 - iirAmount)) + (inputSampleL * iirAmount);
+			inputSampleL -= iirDriveSampleBL;
+			iirDriveSampleDL = (iirDriveSampleDL * (1 - iirAmount)) + (inputSampleL * iirAmount);
+			inputSampleL -= iirDriveSampleDL;
+			iirDriveSampleFL = (iirDriveSampleFL * (1 - iirAmount)) + (inputSampleL * iirAmount);
+			inputSampleL -= iirDriveSampleFL;
+			ataLowpassL = tempSample - inputSampleL;
+			
+			tempSample = inputSampleR;
+			iirDriveSampleBR = (iirDriveSampleBR * (1 - iirAmount)) + (inputSampleR * iirAmount);
+			inputSampleR -= iirDriveSampleBR;
+			iirDriveSampleDR = (iirDriveSampleDR * (1 - iirAmount)) + (inputSampleR * iirAmount);
+			inputSampleR -= iirDriveSampleDR;
+			iirDriveSampleFR = (iirDriveSampleFR * (1 - iirAmount)) + (inputSampleR * iirAmount);
+			inputSampleR -= iirDriveSampleFR;
+			ataLowpassR = tempSample - inputSampleR;
+		}
+		//highpass section
+		
+		if (inputSampleL > 1.0) {inputSampleL = 1.0;}
+		if (inputSampleL < -1.0) {inputSampleL = -1.0;}
+		if (inputSampleR > 1.0) {inputSampleR = 1.0;}
+		if (inputSampleR < -1.0) {inputSampleR = -1.0;}
+		
+		out = driveone;
+		while (out > glitch)
+		{
+			out -= glitch;
+			inputSampleL -= (inputSampleL * (fabs(inputSampleL) * glitch) * (fabs(inputSampleL) * glitch) );
+			inputSampleL *= (1.0+glitch);
+		}
+		//that's taken care of the really high gain stuff
+		inputSampleL -= (inputSampleL * (fabs(inputSampleL) * out) * (fabs(inputSampleL) * out) );
+		inputSampleL *= (1.0+out);
+		
+		
+		out = driveone;
+		while (out > glitch)
+		{
+			out -= glitch;
+			inputSampleR -= (inputSampleR * (fabs(inputSampleR) * glitch) * (fabs(inputSampleR) * glitch) );
+			inputSampleR *= (1.0+glitch);
+		}
+		//that's taken care of the really high gain stuff
+		inputSampleR -= (inputSampleR * (fabs(inputSampleR) * out) * (fabs(inputSampleR) * out) );
+		inputSampleR *= (1.0+out);
+		
+		ataLowpass = (ataLowpassL + ataLowpassR) / 2.0;
+		
+		if (ataLowpass > 0)
+		{if (WasNegative){SubOctave = !SubOctave;} WasNegative = false;}
+		else {WasNegative = true;}
+		//set up polarities for sub-bass version
+		randy = (rand()/(double)RAND_MAX)*fuzz; //0 to 1 the noise, may not be needed
+		invrandy = (1.0-randy);
+		randy /= 2.0;
+		//set up the noise
+		
+		iirSampleA = (iirSampleA * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleA;
+		iirSampleB = (iirSampleB * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleB;
+		iirSampleC = (iirSampleC * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleC;
+		iirSampleD = (iirSampleD * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleD;
+		iirSampleE = (iirSampleE * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleE;
+		iirSampleF = (iirSampleF * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleF;
+		iirSampleG = (iirSampleG * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleG;
+		iirSampleH = (iirSampleH * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleH;
+		iirSampleI = (iirSampleI * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleI;
+		iirSampleJ = (iirSampleJ * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleJ;
+		iirSampleK = (iirSampleK * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleK;
+		iirSampleL = (iirSampleL * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleL;
+		iirSampleM = (iirSampleM * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleM;
+		iirSampleN = (iirSampleN * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleN;
+		iirSampleO = (iirSampleO * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleO;
+		iirSampleP = (iirSampleP * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleP;
+		iirSampleQ = (iirSampleQ * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleQ;
+		iirSampleR = (iirSampleR * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleR;
+		iirSampleS = (iirSampleS * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleS;
+		iirSampleT = (iirSampleT * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleT;
+		iirSampleU = (iirSampleU * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleU;
+		iirSampleV = (iirSampleV * altBmount) + (ataLowpass * iirBmount); ataLowpass -= iirSampleV;
+		
+		switch (bflip)
+		{
+			case 1:				
+				iirHeadBumpA += (ataLowpass * BassGain);
+				iirHeadBumpA -= (iirHeadBumpA * iirHeadBumpA * iirHeadBumpA * HeadBumpFreq);
+				iirHeadBumpA = (invrandy * iirHeadBumpA) + (randy * iirHeadBumpB) + (randy * iirHeadBumpC);
+				if (iirHeadBumpA > 0) iirHeadBumpA -= clamp;
+				if (iirHeadBumpA < 0) iirHeadBumpA += clamp;
+				HeadBump = iirHeadBumpA;
+				break;
+			case 2:
+				iirHeadBumpB += (ataLowpass * BassGain);
+				iirHeadBumpB -= (iirHeadBumpB * iirHeadBumpB * iirHeadBumpB * HeadBumpFreq);
+				iirHeadBumpB = (randy * iirHeadBumpA) + (invrandy * iirHeadBumpB) + (randy * iirHeadBumpC);
+				if (iirHeadBumpB > 0) iirHeadBumpB -= clamp;
+				if (iirHeadBumpB < 0) iirHeadBumpB += clamp;
+				HeadBump = iirHeadBumpB;
+				break;
+			case 3:
+				iirHeadBumpC += (ataLowpass * BassGain);
+				iirHeadBumpC -= (iirHeadBumpC * iirHeadBumpC * iirHeadBumpC * HeadBumpFreq);
+				iirHeadBumpC = (randy * iirHeadBumpA) + (randy * iirHeadBumpB) + (invrandy * iirHeadBumpC);
+				if (iirHeadBumpC > 0) iirHeadBumpC -= clamp;
+				if (iirHeadBumpC < 0) iirHeadBumpC += clamp;
+				HeadBump = iirHeadBumpC;
+				break;
+		}
+		
+		iirSampleW = (iirSampleW * altBmount) + (HeadBump * iirBmount); HeadBump -= iirSampleW;
+		iirSampleX = (iirSampleX * altBmount) + (HeadBump * iirBmount); HeadBump -= iirSampleX;
+		
+		SubBump = HeadBump;
+		iirSampleY = (iirSampleY * altCmount) + (SubBump * iirCmount); SubBump -= iirSampleY;
+		
+		SubBump = fabs(SubBump);
+		if (SubOctave == false) {SubBump = -SubBump;}
+		
+		switch (bflip)
+		{
+			case 1:				
+				iirSubBumpA += (SubBump * SubGain);
+				iirSubBumpA -= (iirSubBumpA * iirSubBumpA * iirSubBumpA * SubBumpFreq);
+				iirSubBumpA = (invrandy * iirSubBumpA) + (randy * iirSubBumpB) + (randy * iirSubBumpC);
+				if (iirSubBumpA > 0) iirSubBumpA -= clamp;
+				if (iirSubBumpA < 0) iirSubBumpA += clamp;
+				SubBump = iirSubBumpA;
+				break;
+			case 2:
+				iirSubBumpB += (SubBump * SubGain);
+				iirSubBumpB -= (iirSubBumpB * iirSubBumpB * iirSubBumpB * SubBumpFreq);
+				iirSubBumpB = (randy * iirSubBumpA) + (invrandy * iirSubBumpB) + (randy * iirSubBumpC);
+				if (iirSubBumpB > 0) iirSubBumpB -= clamp;
+				if (iirSubBumpB < 0) iirSubBumpB += clamp;
+				SubBump = iirSubBumpB;
+				break;
+			case 3:
+				iirSubBumpC += (SubBump * SubGain);
+				iirSubBumpC -= (iirSubBumpC * iirSubBumpC * iirSubBumpC * SubBumpFreq);
+				iirSubBumpC = (randy * iirSubBumpA) + (randy * iirSubBumpB) + (invrandy * iirSubBumpC);
+				if (iirSubBumpC > 0) iirSubBumpC -= clamp;
+				if (iirSubBumpC < 0) iirSubBumpC += clamp;
+				SubBump = iirSubBumpC;
+				break;
+		}
+		
+		iirSampleZ = (iirSampleZ * altCmount) + (SubBump * iirCmount); SubBump -= iirSampleZ;
+		
+		inputSampleL *= driveoutput; //start with the drive section then add lows and subs
+		inputSampleR *= driveoutput; //start with the drive section then add lows and subs
+		
+		inputSampleL += ((HeadBump + lastHeadBump) * BassOutGain);
+		inputSampleL += ((SubBump + lastSubBump) * SubOutGain);
+		
+		inputSampleR += ((HeadBump + lastHeadBump) * BassOutGain);
+		inputSampleR += ((SubBump + lastSubBump) * SubOutGain);
+		
+		lastHeadBump = HeadBump;
+		lastSubBump = SubBump;
+		
+		if (wet != 1.0) {
+			inputSampleL = (inputSampleL * wet) + (drySampleL * dry);
+			inputSampleR = (inputSampleR * wet) + (drySampleR * dry);
+		}
+		//Dry/Wet control, defaults to the last slider
+		
+		flip = !flip;
+		bflip++;
+		if (bflip < 1 || bflip > 3) bflip = 1;
+		
+		//noise shaping to 64-bit floating point
+		double fpTemp = inputSampleL;
+		fpNShapeL += (inputSampleL-fpTemp);
+		inputSampleL += fpNShapeL;
+		//if this confuses you look at the wordlength for fpTemp :)
+		fpTemp = inputSampleR;
+		fpNShapeR += (inputSampleR-fpTemp);
+		inputSampleR += fpNShapeR;
+		//for deeper space and warmth, we try a non-oscillating noise shaping
+		//that is kind of ruthless: it will forever retain the rounding errors
+		//except we'll dial it back a hair at the end of every buffer processed
+		//end noise shaping on 64 bit output
+		
+		*out1 = inputSampleL;
+		*out2 = inputSampleR;
+
+		*in1++;
+		*in2++;
+		*out1++;
+		*out2++;
+    }
+	fpNShapeL *= 0.999999;
+	fpNShapeR *= 0.999999;
+	//we will just delicately dial back the FP noise shaping, not even every sample
+	//this is a good place to put subtle 'no runaway' calculations, though bear in mind
+	//that it will be called more often when you use shorter sample buffers in the DAW.
+	//So, very low latency operation will call these calculations more often.	
+}