Updates (in case my plane crashes)

3 files changed, 1215 insertions, 0 deletions

diff --git a/plugins/LinuxVST/src/IronOxide5/IronOxide5.cpp b/plugins/LinuxVST/src/IronOxide5/IronOxide5.cpp
new file mode 100755
index 0000000..2d7c84f
--- /dev/null
+++ b/plugins/LinuxVST/src/IronOxide5/IronOxide5.cpp
@@ -0,0 +1,196 @@
+/* ========================================
+ *  IronOxide5 - IronOxide5.h
+ *  Copyright (c) 2016 airwindows, All rights reserved
+ * ======================================== */
+
+#ifndef __IronOxide5_H
+#include "IronOxide5.h"
+#endif
+
+AudioEffect* createEffectInstance(audioMasterCallback audioMaster) {return new IronOxide5(audioMaster);}
+
+IronOxide5::IronOxide5(audioMasterCallback audioMaster) :
+    AudioEffectX(audioMaster, kNumPrograms, kNumParameters)
+{
+	A = 0.5; //0.0 input trim in dB -18 to +18, default 0 ((A*36.0)-18.0)
+	B = 0.562341325190349; //15.0 ips 1.5 to 150.0 logarithmic. B*B
+	C = 0.562341325190349; // (C*C) * (C*C) * 150 gives 15 ips (clamp to 1.5 minimum)
+	//15.0  (0.316227766016838)squared * 150 gives 15 ips
+	D = 0.5; //0.5 flutter 0 to 1
+	E = 0.5; //0.5 noise 0 to 1
+	F = 0.5; //0.0 output trim in dB -18 to +18, default 0 ((E*36.0)-18.0)
+	G = 1.0; //1.0 inv/dry/wet -1 0 1 ((F*2.0)-1.0)
+	//needs very fussy defaults to comply with unusual defaults
+	
+	for (int temp = 0; temp < 263; temp++) {dL[temp] = 0.0; dR[temp] = 0.0;}
+	gcount = 0;
+	
+	fastIIRAL = fastIIRBL = slowIIRAL = slowIIRBL = 0.0;
+	fastIIHAL = fastIIHBL = slowIIHAL = slowIIHBL = 0.0;
+	iirSamplehAL = iirSamplehBL = 0.0;
+	iirSampleAL = iirSampleBL = 0.0;
+	prevInputSampleL = 0.0;
+	
+	fastIIRAR = fastIIRBR = slowIIRAR = slowIIRBR = 0.0;
+	fastIIHAR = fastIIHBR = slowIIHAR = slowIIHBR = 0.0;
+	iirSamplehAR = iirSamplehBR = 0.0;
+	iirSampleAR = iirSampleBR = 0.0;
+	prevInputSampleR = 0.0;
+		
+	flip = false;
+	for (int temp = 0; temp < 99; temp++) {flL[temp] = 0.0; flR[temp] = 0.0;}
+	
+	fstoredcount = 0;	
+	sweep = 0.0;
+	rateof = 0.5;
+	nextmax = 0.5;
+	
+	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
+}
+
+IronOxide5::~IronOxide5() {}
+VstInt32 IronOxide5::getVendorVersion () {return 1000;}
+void IronOxide5::setProgramName(char *name) {vst_strncpy (_programName, name, kVstMaxProgNameLen);}
+void IronOxide5::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 IronOxide5::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;
+	/* 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 IronOxide5::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]);
+	/* 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 IronOxide5::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;
+        default: throw; // unknown parameter, shouldn't happen!
+    }
+}
+
+float IronOxide5::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;
+        default: break; // unknown parameter, shouldn't happen!
+    } return 0.0; //we only need to update the relevant name, this is simple to manage
+}
+
+void IronOxide5::getParameterName(VstInt32 index, char *text) {
+    switch (index) {
+        case kParamA: vst_strncpy (text, "Input Trim", kVstMaxParamStrLen); break;
+		case kParamB: vst_strncpy (text, "Tape High", kVstMaxParamStrLen); break;
+		case kParamC: vst_strncpy (text, "Tape Low", kVstMaxParamStrLen); break;
+		case kParamD: vst_strncpy (text, "Flutter", kVstMaxParamStrLen); break;
+		case kParamE: vst_strncpy (text, "Noise", kVstMaxParamStrLen); break;
+		case kParamF: vst_strncpy (text, "Output Trim", kVstMaxParamStrLen); break;
+		case kParamG: vst_strncpy (text, "Inv/Dry/Wet", kVstMaxParamStrLen); break;
+        default: break; // unknown parameter, shouldn't happen!
+    } //this is our labels for displaying in the VST host
+}
+
+void IronOxide5::getParameterDisplay(VstInt32 index, char *text) {
+    switch (index) {
+        case kParamA: float2string (((A*36.0)-18.0), text, kVstMaxParamStrLen); break;
+        case kParamB: float2string (((B*B)*(B*B)*148.5)+1.5, text, kVstMaxParamStrLen); break;
+        case kParamC: float2string (((C*C)*(C*C)*148.5)+1.5, text, kVstMaxParamStrLen); break;
+        case kParamD: float2string (D, text, kVstMaxParamStrLen); break;
+        case kParamE: float2string (E, text, kVstMaxParamStrLen); break;
+        case kParamF: float2string (((F*36.0)-18.0), text, kVstMaxParamStrLen); break;
+        case kParamG: float2string (((G*2.0)-1.0), text, kVstMaxParamStrLen); break;
+        default: break; // unknown parameter, shouldn't happen!
+	} //this displays the values and handles 'popups' where it's discrete choices
+}
+
+void IronOxide5::getParameterLabel(VstInt32 index, char *text) {
+    switch (index) {
+        case kParamA: vst_strncpy (text, "dB", kVstMaxParamStrLen); break;
+        case kParamB: vst_strncpy (text, "ips", kVstMaxParamStrLen); break;
+        case kParamC: vst_strncpy (text, "ips", kVstMaxParamStrLen); break;
+        case kParamD: vst_strncpy (text, "", kVstMaxParamStrLen); break;
+        case kParamE: vst_strncpy (text, "", kVstMaxParamStrLen); break;
+        case kParamF: vst_strncpy (text, "dB", kVstMaxParamStrLen); break;
+        case kParamG: vst_strncpy (text, "", kVstMaxParamStrLen); break;
+        default: break; // unknown parameter, shouldn't happen!
+    }
+}
+
+VstInt32 IronOxide5::canDo(char *text) 
+{ return (_canDo.find(text) == _canDo.end()) ? -1: 1; } // 1 = yes, -1 = no, 0 = don't know
+
+bool IronOxide5::getEffectName(char* name) {
+    vst_strncpy(name, "IronOxide5", kVstMaxProductStrLen); return true;
+}
+
+VstPlugCategory IronOxide5::getPlugCategory() {return kPlugCategEffect;}
+
+bool IronOxide5::getProductString(char* text) {
+  	vst_strncpy (text, "airwindows IronOxide5", kVstMaxProductStrLen); return true;
+}
+
+bool IronOxide5::getVendorString(char* text) {
+  	vst_strncpy (text, "airwindows", kVstMaxVendorStrLen); return true;
+}
diff --git a/plugins/LinuxVST/src/IronOxide5/IronOxide5.h b/plugins/LinuxVST/src/IronOxide5/IronOxide5.h
new file mode 100755
index 0000000..a7cad5d
--- /dev/null
+++ b/plugins/LinuxVST/src/IronOxide5/IronOxide5.h
@@ -0,0 +1,121 @@
+/* ========================================
+ *  IronOxide5 - IronOxide5.h
+ *  Created 8/12/11 by SPIAdmin 
+ *  Copyright (c) 2011 __MyCompanyName__, All rights reserved
+ * ======================================== */
+
+#ifndef __IronOxide5_H
+#define __IronOxide5_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,
+	kNumParameters = 7
+	
+}; //
+
+const int kNumPrograms = 0;
+const int kNumInputs = 2;
+const int kNumOutputs = 2;
+const unsigned long kUniqueId = 'feob';    //Change this to what the AU identity is!
+
+class IronOxide5 : 
+    public AudioEffectX 
+{
+public:
+    IronOxide5(audioMasterCallback audioMaster);
+    ~IronOxide5();
+    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 iirSamplehAL;
+	double iirSamplehBL;
+	double iirSampleAL;
+	double iirSampleBL;
+	double dL[264];
+	double fastIIRAL;
+	double fastIIRBL;
+	double slowIIRAL;
+	double slowIIRBL;
+	double fastIIHAL;
+	double fastIIHBL;
+	double slowIIHAL;
+	double slowIIHBL;
+	double prevInputSampleL;
+
+	double iirSamplehAR;
+	double iirSamplehBR;
+	double iirSampleAR;
+	double iirSampleBR;
+	double dR[264];
+	double fastIIRAR;
+	double fastIIRBR;
+	double slowIIRAR;
+	double slowIIRBR;
+	double fastIIHAR;
+	double fastIIHBR;
+	double slowIIHAR;
+	double slowIIHBR;
+	double prevInputSampleR;
+	
+	int gcount;
+	bool flip;
+	
+	double flL[100];
+	double flR[100];
+	
+	int fstoredcount;
+	double rateof;
+	double sweep;
+	double nextmax;
+	
+    
+	long double fpNShapeLA;
+	long double fpNShapeLB;
+	long double fpNShapeRA;
+	long double fpNShapeRB;
+	bool fpFlip;
+	//default stuff
+
+    float A;
+    float B;
+    float C;
+    float D;
+    float E;
+    float F;
+    float G;
+
+};
+
+#endif
diff --git a/plugins/LinuxVST/src/IronOxide5/IronOxide5Proc.cpp b/plugins/LinuxVST/src/IronOxide5/IronOxide5Proc.cpp
new file mode 100755
index 0000000..c0fbb35
--- /dev/null
+++ b/plugins/LinuxVST/src/IronOxide5/IronOxide5Proc.cpp
@@ -0,0 +1,898 @@
+/* ========================================
+ *  IronOxide5 - IronOxide5.h
+ *  Copyright (c) 2016 airwindows, All rights reserved
+ * ======================================== */
+
+#ifndef __IronOxide5_H
+#include "IronOxide5.h"
+#endif
+
+void IronOxide5::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();
+	float fpTemp;
+	long double fpOld = 0.618033988749894848204586; //golden ratio!
+	long double fpNew = 1.0 - fpOld;	
+
+	double inputgain = pow(10.0,((A*36.0)-18.0)/20.0);
+	double outputgain = pow(10.0,((F*36.0)-18.0)/20.0);
+	double ips = (((B*B)*(B*B)*148.5)+1.5) * 1.1;
+	//slight correction to dial in convincing ips settings
+	if (ips < 1 || ips > 200){ips=33.0;}
+	//sanity checks are always key
+	double tempRandy = 0.04+(0.11/sqrt(ips));
+	double randy;
+	double lps = (((C*C)*(C*C)*148.5)+1.5) * 1.1;
+	//slight correction to dial in convincing ips settings
+	if (lps < 1 || lps > 200){lps=33.0;}
+	//sanity checks are always key
+	double iirAmount = lps/430.0; //for low leaning
+	double bridgerectifierL;
+	double bridgerectifierR;
+	double fastTaper = ips/15.0;
+	double slowTaper = 2.0/(lps*lps);
+	double lowspeedscale = (5.0/ips);
+	int count;
+	int flutcount;
+	double flutterrandy;
+	double temp;
+	double depth = pow(D,2)*overallscale;
+	double fluttertrim = 0.00581/overallscale;
+	double sweeptrim = (0.0005*depth)/overallscale;
+	double offset;	
+	double flooroffset;	
+	double tupi = 3.141592653589793238 * 2.0;
+	double newrate = 0.006/overallscale;
+	double oldrate = 1.0-newrate;	
+	if (overallscale == 0) {fastTaper += 1.0; slowTaper += 1.0;}
+	else
+	{
+		iirAmount /= overallscale;
+		lowspeedscale *= overallscale;
+		fastTaper = 1.0 + (fastTaper / overallscale);
+		slowTaper = 1.0 + (slowTaper / overallscale);
+	}
+	
+	double noise = E * 0.5;
+	double invdrywet = ((G*2.0)-1.0);
+	double dry = 1.0;
+	if (invdrywet > 0.0) dry -= invdrywet;
+	
+	long double inputSampleL;
+	long double inputSampleR;
+	double drySampleL;
+	double drySampleR;
+	    
+    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;
+
+		flutterrandy = (rand()/(double)RAND_MAX);
+		//part of flutter section
+		//now we've got a random flutter, so we're messing with the pitch before tape effects go on
+		if (fstoredcount < 0 || fstoredcount > 30) {fstoredcount = 30;}
+		flutcount = fstoredcount;
+		flL[flutcount+31] = flL[flutcount] = inputSampleL;
+		flR[flutcount+31] = flR[flutcount] = inputSampleR;
+		offset = (1.0 + sin(sweep)) * depth;
+		flooroffset = offset-floor(offset);
+		flutcount += (int)floor(offset);
+		bridgerectifierL = (flL[flutcount] * (1-flooroffset));
+		bridgerectifierR = (flR[flutcount] * (1-flooroffset));
+		bridgerectifierL += (flL[flutcount+1] * flooroffset);
+		bridgerectifierR += (flR[flutcount+1] * flooroffset);
+		rateof = (nextmax * newrate) + (rateof * oldrate);
+		sweep += rateof * fluttertrim;
+		sweep += sweep * sweeptrim;
+		if (sweep >= tupi){sweep = 0.0; nextmax = 0.02 + (flutterrandy*0.98);}
+		fstoredcount--;
+		inputSampleL = bridgerectifierL;
+		inputSampleR = bridgerectifierR;
+		//apply to input signal, interpolate samples
+		//all the funky renaming is just trying to fix how I never reassigned the control numbers
+		
+		if (flip)
+		{
+			iirSampleAL = (iirSampleAL * (1 - iirAmount)) + (inputSampleL * iirAmount);
+			iirSampleAR = (iirSampleAR * (1 - iirAmount)) + (inputSampleR * iirAmount);
+			inputSampleL -= iirSampleAL;
+			inputSampleR -= iirSampleAR;
+		}
+		else
+		{
+			iirSampleBL = (iirSampleBL * (1 - iirAmount)) + (inputSampleL * iirAmount);
+			iirSampleBR = (iirSampleBR * (1 - iirAmount)) + (inputSampleR * iirAmount);
+			inputSampleL -= iirSampleBL;
+			inputSampleR -= iirSampleBR;
+		}
+		//do IIR highpass for leaning out
+		
+		inputSampleL *= inputgain;
+		inputSampleR *= inputgain;
+		
+		bridgerectifierL = fabs(inputSampleL);
+		if (bridgerectifierL > 1.57079633) bridgerectifierL = 1.57079633;
+		bridgerectifierL = sin(bridgerectifierL);
+		if (inputSampleL > 0.0) inputSampleL = bridgerectifierL;
+		else inputSampleL = -bridgerectifierL;
+		//preliminary gain stage using antialiasing
+		
+		bridgerectifierR = fabs(inputSampleR);
+		if (bridgerectifierR > 1.57079633) bridgerectifierR = 1.57079633;
+		bridgerectifierR = sin(bridgerectifierR);
+		if (inputSampleR > 0.0) inputSampleR = bridgerectifierR;
+		else inputSampleR = -bridgerectifierR;
+		//preliminary gain stage using antialiasing
+		
+		//over to the Iron Oxide shaping code using inputsample
+		if (gcount < 0 || gcount > 131) {gcount = 131;}
+		count = gcount;
+		//increment the counter
+		
+		dL[count+131] = dL[count] = inputSampleL;
+		dR[count+131] = dR[count] = inputSampleR;
+		
+		if (flip)
+		{
+			fastIIRAL = fastIIRAL/fastTaper;
+			slowIIRAL = slowIIRAL/slowTaper;
+			fastIIRAL += dL[count];
+			//scale stuff down
+			
+			fastIIRAR = fastIIRAR/fastTaper;
+			slowIIRAR = slowIIRAR/slowTaper;
+			fastIIRAR += dR[count];
+			//scale stuff down
+			count += 3;
+			
+			temp = dL[count+127];
+			temp += dL[count+113];
+			temp += dL[count+109];
+			temp += dL[count+107];
+			temp += dL[count+103];
+			temp += dL[count+101];
+			temp += dL[count+97];
+			temp += dL[count+89];
+			temp += dL[count+83];
+			temp /= 2;
+			temp += dL[count+79];
+			temp += dL[count+73];
+			temp += dL[count+71];
+			temp += dL[count+67];
+			temp += dL[count+61];
+			temp += dL[count+59];
+			temp += dL[count+53];
+			temp += dL[count+47];
+			temp += dL[count+43];
+			temp += dL[count+41];
+			temp += dL[count+37];
+			temp += dL[count+31];
+			temp += dL[count+29];
+			temp /= 2;
+			temp += dL[count+23];
+			temp += dL[count+19];
+			temp += dL[count+17];
+			temp += dL[count+13];
+			temp += dL[count+11];
+			temp /= 2;
+			temp += dL[count+7];
+			temp += dL[count+5];
+			temp += dL[count+3];
+			temp /= 2;
+			temp += dL[count+2];
+			temp += dL[count+1]; //end L
+			slowIIRAL += (temp/128);
+		
+			temp = dR[count+127];
+			temp += dR[count+113];
+			temp += dR[count+109];
+			temp += dR[count+107];
+			temp += dR[count+103];
+			temp += dR[count+101];
+			temp += dR[count+97];
+			temp += dR[count+89];
+			temp += dR[count+83];
+			temp /= 2;
+			temp += dR[count+79];
+			temp += dR[count+73];
+			temp += dR[count+71];
+			temp += dR[count+67];
+			temp += dR[count+61];
+			temp += dR[count+59];
+			temp += dR[count+53];
+			temp += dR[count+47];
+			temp += dR[count+43];
+			temp += dR[count+41];
+			temp += dR[count+37];
+			temp += dR[count+31];
+			temp += dR[count+29];
+			temp /= 2;
+			temp += dR[count+23];
+			temp += dR[count+19];
+			temp += dR[count+17];
+			temp += dR[count+13];
+			temp += dR[count+11];
+			temp /= 2;
+			temp += dR[count+7];
+			temp += dR[count+5];
+			temp += dR[count+3];
+			temp /= 2;
+			temp += dR[count+2];
+			temp += dR[count+1]; //end R
+			slowIIRAR += (temp/128);
+			
+			inputSampleL = fastIIRAL - (slowIIRAL / slowTaper);
+			inputSampleR = fastIIRAR - (slowIIRAR / slowTaper);
+		}
+		else
+		{
+			fastIIRBL = fastIIRBL/fastTaper;
+			slowIIRBL = slowIIRBL/slowTaper;
+			fastIIRBL += dL[count];
+			//scale stuff down
+			
+			fastIIRBR = fastIIRBR/fastTaper;
+			slowIIRBR = slowIIRBR/slowTaper;
+			fastIIRBR += dR[count];
+			//scale stuff down
+			count += 3;
+			
+			temp = dL[count+127];
+			temp += dL[count+113];
+			temp += dL[count+109];
+			temp += dL[count+107];
+			temp += dL[count+103];
+			temp += dL[count+101];
+			temp += dL[count+97];
+			temp += dL[count+89];
+			temp += dL[count+83];
+			temp /= 2;
+			temp += dL[count+79];
+			temp += dL[count+73];
+			temp += dL[count+71];
+			temp += dL[count+67];
+			temp += dL[count+61];
+			temp += dL[count+59];
+			temp += dL[count+53];
+			temp += dL[count+47];
+			temp += dL[count+43];
+			temp += dL[count+41];
+			temp += dL[count+37];
+			temp += dL[count+31];
+			temp += dL[count+29];
+			temp /= 2;
+			temp += dL[count+23];
+			temp += dL[count+19];
+			temp += dL[count+17];
+			temp += dL[count+13];
+			temp += dL[count+11];
+			temp /= 2;
+			temp += dL[count+7];
+			temp += dL[count+5];
+			temp += dL[count+3];
+			temp /= 2;
+			temp += dL[count+2];
+			temp += dL[count+1];
+			slowIIRBL += (temp/128);
+		
+			temp = dR[count+127];
+			temp += dR[count+113];
+			temp += dR[count+109];
+			temp += dR[count+107];
+			temp += dR[count+103];
+			temp += dR[count+101];
+			temp += dR[count+97];
+			temp += dR[count+89];
+			temp += dR[count+83];
+			temp /= 2;
+			temp += dR[count+79];
+			temp += dR[count+73];
+			temp += dR[count+71];
+			temp += dR[count+67];
+			temp += dR[count+61];
+			temp += dR[count+59];
+			temp += dR[count+53];
+			temp += dR[count+47];
+			temp += dR[count+43];
+			temp += dR[count+41];
+			temp += dR[count+37];
+			temp += dR[count+31];
+			temp += dR[count+29];
+			temp /= 2;
+			temp += dR[count+23];
+			temp += dR[count+19];
+			temp += dR[count+17];
+			temp += dR[count+13];
+			temp += dR[count+11];
+			temp /= 2;
+			temp += dR[count+7];
+			temp += dR[count+5];
+			temp += dR[count+3];
+			temp /= 2;
+			temp += dR[count+2];
+			temp += dR[count+1];
+			slowIIRBR += (temp/128);
+			
+			inputSampleL = fastIIRBL - (slowIIRBL / slowTaper);
+			inputSampleR = fastIIRBR - (slowIIRBR / slowTaper);
+		}
+		
+		inputSampleL /= fastTaper;
+		inputSampleR /= fastTaper;
+		inputSampleL /= lowspeedscale;
+		inputSampleR /= lowspeedscale;
+		
+		bridgerectifierL = fabs(inputSampleL);
+		if (bridgerectifierL > 1.57079633) bridgerectifierL = 1.57079633;
+		bridgerectifierL = sin(bridgerectifierL);
+		//can use as an output limiter
+		if (inputSampleL > 0.0) inputSampleL = bridgerectifierL;
+		else inputSampleL = -bridgerectifierL;
+		//second stage of overdrive to prevent overs and allow bloody loud extremeness		
+
+		bridgerectifierR = fabs(inputSampleR);
+		if (bridgerectifierR > 1.57079633) bridgerectifierR = 1.57079633;
+		bridgerectifierR = sin(bridgerectifierR);
+		//can use as an output limiter
+		if (inputSampleR > 0.0) inputSampleR = bridgerectifierR;
+		else inputSampleR = -bridgerectifierR;
+		//second stage of overdrive to prevent overs and allow bloody loud extremeness		
+		
+		randy = (0.55 + tempRandy + ((rand()/(double)RAND_MAX)*tempRandy))*noise; //0 to 2
+		
+		inputSampleL *= (1.0 - randy);
+		inputSampleL += (prevInputSampleL*randy);
+		prevInputSampleL = drySampleL;
+		
+		inputSampleR *= (1.0 - randy);
+		inputSampleR += (prevInputSampleR*randy);
+		prevInputSampleR = drySampleR;
+		//a slew-based noise generator: noise is only relative to how much change between samples there is.
+		//This will cause high sample rates to be a little 'smoother' and clearer. I see this mechanic as something that
+		//interacts with the sample rate. The 'dust' is scaled to the size of the samples.
+		
+		flip = !flip;
+		
+		//begin invdrywet block with outputgain
+		if (outputgain != 1.0) {
+			inputSampleL *= outputgain;
+			inputSampleR *= outputgain;
+		}
+		if (invdrywet != 1.0) {
+			inputSampleL *= invdrywet;
+			inputSampleR *= invdrywet;
+		}
+		if (dry != 1.0) {
+			drySampleL *= dry;
+			drySampleR *= dry;
+		}
+		if (fabs(drySampleL) > 0.0) {
+			inputSampleL += drySampleL;
+		}
+		if (fabs(drySampleR) > 0.0) {
+			inputSampleR += drySampleR;
+		}
+		//end invdrywet block with outputgain
+		
+		//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 IronOxide5::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 fpTemp;
+	long double fpOld = 0.618033988749894848204586; //golden ratio!
+	long double fpNew = 1.0 - fpOld;	
+	
+	double inputgain = pow(10.0,((A*36.0)-18.0)/20.0);
+	double outputgain = pow(10.0,((F*36.0)-18.0)/20.0);
+	double ips = (((B*B)*(B*B)*148.5)+1.5) * 1.1;
+	//slight correction to dial in convincing ips settings
+	if (ips < 1 || ips > 200){ips=33.0;}
+	//sanity checks are always key
+	double tempRandy = 0.04+(0.11/sqrt(ips));
+	double randy;
+	double lps = (((C*C)*(C*C)*148.5)+1.5) * 1.1;
+	//slight correction to dial in convincing ips settings
+	if (lps < 1 || lps > 200){lps=33.0;}
+	//sanity checks are always key
+	double iirAmount = lps/430.0; //for low leaning
+	double bridgerectifierL;
+	double bridgerectifierR;
+	double fastTaper = ips/15.0;
+	double slowTaper = 2.0/(lps*lps);
+	double lowspeedscale = (5.0/ips);
+	int count;
+	int flutcount;
+	double flutterrandy;
+	double temp;
+	double depth = pow(D,2)*overallscale;
+	double fluttertrim = 0.00581/overallscale;
+	double sweeptrim = (0.0005*depth)/overallscale;
+	double offset;	
+	double flooroffset;	
+	double tupi = 3.141592653589793238 * 2.0;
+	double newrate = 0.006/overallscale;
+	double oldrate = 1.0-newrate;	
+	if (overallscale == 0) {fastTaper += 1.0; slowTaper += 1.0;}
+	else
+	{
+		iirAmount /= overallscale;
+		lowspeedscale *= overallscale;
+		fastTaper = 1.0 + (fastTaper / overallscale);
+		slowTaper = 1.0 + (slowTaper / overallscale);
+	}
+	
+	double noise = E * 0.5;
+	double invdrywet = ((G*2.0)-1.0);
+	double dry = 1.0;
+	if (invdrywet > 0.0) dry -= invdrywet;
+	
+	long double inputSampleL;
+	long double inputSampleR;
+	double drySampleL;
+	double drySampleR;
+	
+    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;
+		
+		flutterrandy = (rand()/(double)RAND_MAX);
+		//part of flutter section
+		//now we've got a random flutter, so we're messing with the pitch before tape effects go on
+		if (fstoredcount < 0 || fstoredcount > 30) {fstoredcount = 30;}
+		flutcount = fstoredcount;
+		flL[flutcount+31] = flL[flutcount] = inputSampleL;
+		flR[flutcount+31] = flR[flutcount] = inputSampleR;
+		offset = (1.0 + sin(sweep)) * depth;
+		flooroffset = offset-floor(offset);
+		flutcount += (int)floor(offset);
+		bridgerectifierL = (flL[flutcount] * (1-flooroffset));
+		bridgerectifierR = (flR[flutcount] * (1-flooroffset));
+		bridgerectifierL += (flL[flutcount+1] * flooroffset);
+		bridgerectifierR += (flR[flutcount+1] * flooroffset);
+		rateof = (nextmax * newrate) + (rateof * oldrate);
+		sweep += rateof * fluttertrim;
+		sweep += sweep * sweeptrim;
+		if (sweep >= tupi){sweep = 0.0; nextmax = 0.02 + (flutterrandy*0.98);}
+		fstoredcount--;
+		inputSampleL = bridgerectifierL;
+		inputSampleR = bridgerectifierR;
+		//apply to input signal, interpolate samples
+		//all the funky renaming is just trying to fix how I never reassigned the control numbers
+		
+		if (flip)
+		{
+			iirSampleAL = (iirSampleAL * (1 - iirAmount)) + (inputSampleL * iirAmount);
+			iirSampleAR = (iirSampleAR * (1 - iirAmount)) + (inputSampleR * iirAmount);
+			inputSampleL -= iirSampleAL;
+			inputSampleR -= iirSampleAR;
+		}
+		else
+		{
+			iirSampleBL = (iirSampleBL * (1 - iirAmount)) + (inputSampleL * iirAmount);
+			iirSampleBR = (iirSampleBR * (1 - iirAmount)) + (inputSampleR * iirAmount);
+			inputSampleL -= iirSampleBL;
+			inputSampleR -= iirSampleBR;
+		}
+		//do IIR highpass for leaning out
+		
+		inputSampleL *= inputgain;
+		inputSampleR *= inputgain;
+		
+		bridgerectifierL = fabs(inputSampleL);
+		if (bridgerectifierL > 1.57079633) bridgerectifierL = 1.57079633;
+		bridgerectifierL = sin(bridgerectifierL);
+		if (inputSampleL > 0.0) inputSampleL = bridgerectifierL;
+		else inputSampleL = -bridgerectifierL;
+		//preliminary gain stage using antialiasing
+		
+		bridgerectifierR = fabs(inputSampleR);
+		if (bridgerectifierR > 1.57079633) bridgerectifierR = 1.57079633;
+		bridgerectifierR = sin(bridgerectifierR);
+		if (inputSampleR > 0.0) inputSampleR = bridgerectifierR;
+		else inputSampleR = -bridgerectifierR;
+		//preliminary gain stage using antialiasing
+		
+		//over to the Iron Oxide shaping code using inputsample
+		if (gcount < 0 || gcount > 131) {gcount = 131;}
+		count = gcount;
+		//increment the counter
+		
+		dL[count+131] = dL[count] = inputSampleL;
+		dR[count+131] = dR[count] = inputSampleR;
+		
+		if (flip)
+		{
+			fastIIRAL = fastIIRAL/fastTaper;
+			slowIIRAL = slowIIRAL/slowTaper;
+			fastIIRAL += dL[count];
+			//scale stuff down
+			
+			fastIIRAR = fastIIRAR/fastTaper;
+			slowIIRAR = slowIIRAR/slowTaper;
+			fastIIRAR += dR[count];
+			//scale stuff down
+			count += 3;
+			
+			temp = dL[count+127];
+			temp += dL[count+113];
+			temp += dL[count+109];
+			temp += dL[count+107];
+			temp += dL[count+103];
+			temp += dL[count+101];
+			temp += dL[count+97];
+			temp += dL[count+89];
+			temp += dL[count+83];
+			temp /= 2;
+			temp += dL[count+79];
+			temp += dL[count+73];
+			temp += dL[count+71];
+			temp += dL[count+67];
+			temp += dL[count+61];
+			temp += dL[count+59];
+			temp += dL[count+53];
+			temp += dL[count+47];
+			temp += dL[count+43];
+			temp += dL[count+41];
+			temp += dL[count+37];
+			temp += dL[count+31];
+			temp += dL[count+29];
+			temp /= 2;
+			temp += dL[count+23];
+			temp += dL[count+19];
+			temp += dL[count+17];
+			temp += dL[count+13];
+			temp += dL[count+11];
+			temp /= 2;
+			temp += dL[count+7];
+			temp += dL[count+5];
+			temp += dL[count+3];
+			temp /= 2;
+			temp += dL[count+2];
+			temp += dL[count+1]; //end L
+			slowIIRAL += (temp/128);
+			
+			temp = dR[count+127];
+			temp += dR[count+113];
+			temp += dR[count+109];
+			temp += dR[count+107];
+			temp += dR[count+103];
+			temp += dR[count+101];
+			temp += dR[count+97];
+			temp += dR[count+89];
+			temp += dR[count+83];
+			temp /= 2;
+			temp += dR[count+79];
+			temp += dR[count+73];
+			temp += dR[count+71];
+			temp += dR[count+67];
+			temp += dR[count+61];
+			temp += dR[count+59];
+			temp += dR[count+53];
+			temp += dR[count+47];
+			temp += dR[count+43];
+			temp += dR[count+41];
+			temp += dR[count+37];
+			temp += dR[count+31];
+			temp += dR[count+29];
+			temp /= 2;
+			temp += dR[count+23];
+			temp += dR[count+19];
+			temp += dR[count+17];
+			temp += dR[count+13];
+			temp += dR[count+11];
+			temp /= 2;
+			temp += dR[count+7];
+			temp += dR[count+5];
+			temp += dR[count+3];
+			temp /= 2;
+			temp += dR[count+2];
+			temp += dR[count+1]; //end R
+			slowIIRAR += (temp/128);
+			
+			inputSampleL = fastIIRAL - (slowIIRAL / slowTaper);
+			inputSampleR = fastIIRAR - (slowIIRAR / slowTaper);
+		}
+		else
+		{
+			fastIIRBL = fastIIRBL/fastTaper;
+			slowIIRBL = slowIIRBL/slowTaper;
+			fastIIRBL += dL[count];
+			//scale stuff down
+			
+			fastIIRBR = fastIIRBR/fastTaper;
+			slowIIRBR = slowIIRBR/slowTaper;
+			fastIIRBR += dR[count];
+			//scale stuff down
+			count += 3;
+			
+			temp = dL[count+127];
+			temp += dL[count+113];
+			temp += dL[count+109];
+			temp += dL[count+107];
+			temp += dL[count+103];
+			temp += dL[count+101];
+			temp += dL[count+97];
+			temp += dL[count+89];
+			temp += dL[count+83];
+			temp /= 2;
+			temp += dL[count+79];
+			temp += dL[count+73];
+			temp += dL[count+71];
+			temp += dL[count+67];
+			temp += dL[count+61];
+			temp += dL[count+59];
+			temp += dL[count+53];
+			temp += dL[count+47];
+			temp += dL[count+43];
+			temp += dL[count+41];
+			temp += dL[count+37];
+			temp += dL[count+31];
+			temp += dL[count+29];
+			temp /= 2;
+			temp += dL[count+23];
+			temp += dL[count+19];
+			temp += dL[count+17];
+			temp += dL[count+13];
+			temp += dL[count+11];
+			temp /= 2;
+			temp += dL[count+7];
+			temp += dL[count+5];
+			temp += dL[count+3];
+			temp /= 2;
+			temp += dL[count+2];
+			temp += dL[count+1];
+			slowIIRBL += (temp/128);
+			
+			temp = dR[count+127];
+			temp += dR[count+113];
+			temp += dR[count+109];
+			temp += dR[count+107];
+			temp += dR[count+103];
+			temp += dR[count+101];
+			temp += dR[count+97];
+			temp += dR[count+89];
+			temp += dR[count+83];
+			temp /= 2;
+			temp += dR[count+79];
+			temp += dR[count+73];
+			temp += dR[count+71];
+			temp += dR[count+67];
+			temp += dR[count+61];
+			temp += dR[count+59];
+			temp += dR[count+53];
+			temp += dR[count+47];
+			temp += dR[count+43];
+			temp += dR[count+41];
+			temp += dR[count+37];
+			temp += dR[count+31];
+			temp += dR[count+29];
+			temp /= 2;
+			temp += dR[count+23];
+			temp += dR[count+19];
+			temp += dR[count+17];
+			temp += dR[count+13];
+			temp += dR[count+11];
+			temp /= 2;
+			temp += dR[count+7];
+			temp += dR[count+5];
+			temp += dR[count+3];
+			temp /= 2;
+			temp += dR[count+2];
+			temp += dR[count+1];
+			slowIIRBR += (temp/128);
+			
+			inputSampleL = fastIIRBL - (slowIIRBL / slowTaper);
+			inputSampleR = fastIIRBR - (slowIIRBR / slowTaper);
+		}
+		
+		inputSampleL /= fastTaper;
+		inputSampleR /= fastTaper;
+		inputSampleL /= lowspeedscale;
+		inputSampleR /= lowspeedscale;
+		
+		bridgerectifierL = fabs(inputSampleL);
+		if (bridgerectifierL > 1.57079633) bridgerectifierL = 1.57079633;
+		bridgerectifierL = sin(bridgerectifierL);
+		//can use as an output limiter
+		if (inputSampleL > 0.0) inputSampleL = bridgerectifierL;
+		else inputSampleL = -bridgerectifierL;
+		//second stage of overdrive to prevent overs and allow bloody loud extremeness		
+		
+		bridgerectifierR = fabs(inputSampleR);
+		if (bridgerectifierR > 1.57079633) bridgerectifierR = 1.57079633;
+		bridgerectifierR = sin(bridgerectifierR);
+		//can use as an output limiter
+		if (inputSampleR > 0.0) inputSampleR = bridgerectifierR;
+		else inputSampleR = -bridgerectifierR;
+		//second stage of overdrive to prevent overs and allow bloody loud extremeness		
+		
+		randy = (0.55 + tempRandy + ((rand()/(double)RAND_MAX)*tempRandy))*noise; //0 to 2
+		
+		inputSampleL *= (1.0 - randy);
+		inputSampleL += (prevInputSampleL*randy);
+		prevInputSampleL = drySampleL;
+		
+		inputSampleR *= (1.0 - randy);
+		inputSampleR += (prevInputSampleR*randy);
+		prevInputSampleR = drySampleR;
+		//a slew-based noise generator: noise is only relative to how much change between samples there is.
+		//This will cause high sample rates to be a little 'smoother' and clearer. I see this mechanic as something that
+		//interacts with the sample rate. The 'dust' is scaled to the size of the samples.
+		
+		flip = !flip;
+		
+		//begin invdrywet block with outputgain
+		if (outputgain != 1.0) {
+			inputSampleL *= outputgain;
+			inputSampleR *= outputgain;
+		}
+		if (invdrywet != 1.0) {
+			inputSampleL *= invdrywet;
+			inputSampleR *= invdrywet;
+		}
+		if (dry != 1.0) {
+			drySampleL *= dry;
+			drySampleR *= dry;
+		}
+		if (fabs(drySampleL) > 0.0) {
+			inputSampleL += drySampleL;
+		}
+		if (fabs(drySampleR) > 0.0) {
+			inputSampleR += drySampleR;
+		}
+		//end invdrywet block with outputgain
+		
+		//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