VariMu

3 files changed, 764 insertions, 0 deletions

diff --git a/plugins/LinuxVST/src/VariMu/VariMu.cpp b/plugins/LinuxVST/src/VariMu/VariMu.cpp
new file mode 100755
index 0000000..c4d63b2
--- /dev/null
+++ b/plugins/LinuxVST/src/VariMu/VariMu.cpp
@@ -0,0 +1,159 @@
+/* ========================================
+ *  VariMu - VariMu.h
+ *  Copyright (c) 2016 airwindows, All rights reserved
+ * ======================================== */
+
+#ifndef __VariMu_H
+#include "VariMu.h"
+#endif
+
+AudioEffect* createEffectInstance(audioMasterCallback audioMaster) {return new VariMu(audioMaster);}
+
+VariMu::VariMu(audioMasterCallback audioMaster) :
+    AudioEffectX(audioMaster, kNumPrograms, kNumParameters)
+{
+	
+	muSpeedAL = 10000;
+	muSpeedBL = 10000;
+	muCoefficientAL = 1;
+	muCoefficientBL = 1;
+	muVaryL = 1;
+	previousL = 0.0;
+	
+	muSpeedAR = 10000;
+	muSpeedBR = 10000;
+	muCoefficientAR = 1;
+	muCoefficientBR = 1;
+	muVaryR = 1;
+	previousR = 0.0;
+	flip = false;
+	
+	A = 0.0;
+	B = 0.5;
+	C = 1.0;
+	D = 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
+}
+
+VariMu::~VariMu() {}
+VstInt32 VariMu::getVendorVersion () {return 1000;}
+void VariMu::setProgramName(char *name) {vst_strncpy (_programName, name, kVstMaxProgNameLen);}
+void VariMu::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 VariMu::getChunk (void** data, bool isPreset)
+{
+	float *chunkData = (float *)calloc(kNumParameters, sizeof(float));
+	chunkData[0] = A;
+	chunkData[1] = B;
+	chunkData[2] = C;
+	chunkData[3] = D;
+	/* 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 VariMu::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]);
+	/* 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 VariMu::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;
+        default: throw; // unknown parameter, shouldn't happen!
+    }
+}
+
+float VariMu::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;
+        default: break; // unknown parameter, shouldn't happen!
+    } return 0.0; //we only need to update the relevant name, this is simple to manage
+}
+
+void VariMu::getParameterName(VstInt32 index, char *text) {
+    switch (index) {
+        case kParamA: vst_strncpy (text, "Intensty", kVstMaxParamStrLen); break;
+		case kParamB: vst_strncpy (text, "Speed", kVstMaxParamStrLen); break;
+		case kParamC: vst_strncpy (text, "Output", kVstMaxParamStrLen); break;
+		case kParamD: 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 VariMu::getParameterDisplay(VstInt32 index, char *text) {
+    switch (index) {
+        case kParamA: float2string (A, text, kVstMaxParamStrLen); break;
+        case kParamB: float2string (B, text, kVstMaxParamStrLen); break;
+        case kParamC: float2string (C, text, kVstMaxParamStrLen); break;
+        case kParamD: float2string (D, text, kVstMaxParamStrLen); break;
+        default: break; // unknown parameter, shouldn't happen!
+	} //this displays the values and handles 'popups' where it's discrete choices
+}
+
+void VariMu::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;
+		default: break; // unknown parameter, shouldn't happen!
+    }
+}
+
+VstInt32 VariMu::canDo(char *text) 
+{ return (_canDo.find(text) == _canDo.end()) ? -1: 1; } // 1 = yes, -1 = no, 0 = don't know
+
+bool VariMu::getEffectName(char* name) {
+    vst_strncpy(name, "VariMu", kVstMaxProductStrLen); return true;
+}
+
+VstPlugCategory VariMu::getPlugCategory() {return kPlugCategEffect;}
+
+bool VariMu::getProductString(char* text) {
+  	vst_strncpy (text, "airwindows VariMu", kVstMaxProductStrLen); return true;
+}
+
+bool VariMu::getVendorString(char* text) {
+  	vst_strncpy (text, "airwindows", kVstMaxVendorStrLen); return true;
+}
diff --git a/plugins/LinuxVST/src/VariMu/VariMu.h b/plugins/LinuxVST/src/VariMu/VariMu.h
new file mode 100755
index 0000000..d79dff1
--- /dev/null
+++ b/plugins/LinuxVST/src/VariMu/VariMu.h
@@ -0,0 +1,87 @@
+/* ========================================
+ *  VariMu - VariMu.h
+ *  Created 8/12/11 by SPIAdmin 
+ *  Copyright (c) 2011 __MyCompanyName__, All rights reserved
+ * ======================================== */
+
+#ifndef __VariMu_H
+#define __VariMu_H
+
+#ifndef __audioeffect__
+#include "audioeffectx.h"
+#endif
+
+#include <set>
+#include <string>
+#include <math.h>
+
+enum {
+	kParamA = 0,
+	kParamB = 1,
+	kParamC = 2,
+	kParamD = 3,
+  kNumParameters = 4
+}; //
+
+const int kNumPrograms = 0;
+const int kNumInputs = 2;
+const int kNumOutputs = 2;
+const unsigned long kUniqueId = 'vari';    //Change this to what the AU identity is!
+
+class VariMu : 
+    public AudioEffectX 
+{
+public:
+    VariMu(audioMasterCallback audioMaster);
+    ~VariMu();
+    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 muVaryL;
+	double muAttackL;
+	double muNewSpeedL;
+	double muSpeedAL;
+	double muSpeedBL;
+	double muCoefficientAL;
+	double muCoefficientBL;
+	double previousL;
+
+	double muVaryR;
+	double muAttackR;
+	double muNewSpeedR;
+	double muSpeedAR;
+	double muSpeedBR;
+	double muCoefficientAR;
+	double muCoefficientBR;
+	double previousR;
+	bool flip;
+	
+	long double fpNShapeL;
+	long double fpNShapeR;
+	//default stuff
+
+    float A;
+    float B;
+    float C;
+    float D;
+};
+
+#endif
diff --git a/plugins/LinuxVST/src/VariMu/VariMuProc.cpp b/plugins/LinuxVST/src/VariMu/VariMuProc.cpp
new file mode 100755
index 0000000..b966cd8
--- /dev/null
+++ b/plugins/LinuxVST/src/VariMu/VariMuProc.cpp
@@ -0,0 +1,518 @@
+/* ========================================
+ *  VariMu - VariMu.h
+ *  Copyright (c) 2016 airwindows, All rights reserved
+ * ======================================== */
+
+#ifndef __VariMu_H
+#include "VariMu.h"
+#endif
+
+void VariMu::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 = 2.0;
+	overallscale /= 44100.0;
+	overallscale *= getSampleRate();
+
+	double threshold = 1.001 - (1.0-pow(1.0-A,3));
+	double muMakeupGain = sqrt(1.0 / threshold);
+	muMakeupGain = (muMakeupGain + sqrt(muMakeupGain))/2.0;
+	muMakeupGain = sqrt(muMakeupGain);
+	double outGain = sqrt(muMakeupGain);
+	//gain settings around threshold
+	double release = pow((1.15-B),5)*32768.0;
+	release /= overallscale;
+	double fastest = sqrt(release);
+	//speed settings around release
+	double coefficient;
+	double output = outGain * C;
+	double wet = D;
+	long double squaredSampleL;
+	long double squaredSampleR;
+	
+	// µ µ µ µ µ µ µ µ µ µ µ µ is the kitten song o/~
+    
+    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 aVariMu. We want a 'air' hiss
+		long double drySampleL = inputSampleL;
+		long double drySampleR = inputSampleR;
+		
+		if (fabs(inputSampleL) > fabs(previousL)) squaredSampleL = previousL * previousL;
+		else squaredSampleL = inputSampleL * inputSampleL;
+		previousL = inputSampleL;
+		inputSampleL *= muMakeupGain;
+
+		if (fabs(inputSampleR) > fabs(previousR)) squaredSampleR = previousR * previousR;
+		else squaredSampleR = inputSampleR * inputSampleR;
+		previousR = inputSampleR;
+		inputSampleR *= muMakeupGain;
+		
+		//adjust coefficients for L
+		if (flip)
+		{
+			if (fabs(squaredSampleL) > threshold)
+			{
+				muVaryL = threshold / fabs(squaredSampleL);
+				muAttackL = sqrt(fabs(muSpeedAL));
+				muCoefficientAL = muCoefficientAL * (muAttackL-1.0);
+				if (muVaryL < threshold)
+				{
+					muCoefficientAL = muCoefficientAL + threshold;
+				}
+				else
+				{
+					muCoefficientAL = muCoefficientAL + muVaryL;
+				}
+				muCoefficientAL = muCoefficientAL / muAttackL;
+			}
+			else
+			{
+				muCoefficientAL = muCoefficientAL * ((muSpeedAL * muSpeedAL)-1.0);
+				muCoefficientAL = muCoefficientAL + 1.0;
+				muCoefficientAL = muCoefficientAL / (muSpeedAL * muSpeedAL);
+			}
+			muNewSpeedL = muSpeedAL * (muSpeedAL-1);
+			muNewSpeedL = muNewSpeedL + fabs(squaredSampleL*release)+fastest;
+			muSpeedAL = muNewSpeedL / muSpeedAL;
+		}
+		else
+		{
+			if (fabs(squaredSampleL) > threshold)
+			{
+				muVaryL = threshold / fabs(squaredSampleL);
+				muAttackL = sqrt(fabs(muSpeedBL));
+				muCoefficientBL = muCoefficientBL * (muAttackL-1);
+				if (muVaryL < threshold)
+				{
+					muCoefficientBL = muCoefficientBL + threshold;
+				}
+				else
+				{
+					muCoefficientBL = muCoefficientBL + muVaryL;
+				}
+				muCoefficientBL = muCoefficientBL / muAttackL;
+			}
+			else
+			{
+				muCoefficientBL = muCoefficientBL * ((muSpeedBL * muSpeedBL)-1.0);
+				muCoefficientBL = muCoefficientBL + 1.0;
+				muCoefficientBL = muCoefficientBL / (muSpeedBL * muSpeedBL);
+			}
+			muNewSpeedL = muSpeedBL * (muSpeedBL-1);
+			muNewSpeedL = muNewSpeedL + fabs(squaredSampleL*release)+fastest;
+			muSpeedBL = muNewSpeedL / muSpeedBL;
+		}
+		//got coefficients, adjusted speeds for L
+		
+		//adjust coefficients for R
+		if (flip)
+		{
+			if (fabs(squaredSampleR) > threshold)
+			{
+				muVaryR = threshold / fabs(squaredSampleR);
+				muAttackR = sqrt(fabs(muSpeedAR));
+				muCoefficientAR = muCoefficientAR * (muAttackR-1.0);
+				if (muVaryR < threshold)
+				{
+					muCoefficientAR = muCoefficientAR + threshold;
+				}
+				else
+				{
+					muCoefficientAR = muCoefficientAR + muVaryR;
+				}
+				muCoefficientAR = muCoefficientAR / muAttackR;
+			}
+			else
+			{
+				muCoefficientAR = muCoefficientAR * ((muSpeedAR * muSpeedAR)-1.0);
+				muCoefficientAR = muCoefficientAR + 1.0;
+				muCoefficientAR = muCoefficientAR / (muSpeedAR * muSpeedAR);
+			}
+			muNewSpeedR = muSpeedAR * (muSpeedAR-1);
+			muNewSpeedR = muNewSpeedR + fabs(squaredSampleR*release)+fastest;
+			muSpeedAR = muNewSpeedR / muSpeedAR;
+		}
+		else
+		{
+			if (fabs(squaredSampleR) > threshold)
+			{
+				muVaryR = threshold / fabs(squaredSampleR);
+				muAttackR = sqrt(fabs(muSpeedBR));
+				muCoefficientBR = muCoefficientBR * (muAttackR-1);
+				if (muVaryR < threshold)
+				{
+					muCoefficientBR = muCoefficientBR + threshold;
+				}
+				else
+				{
+					muCoefficientBR = muCoefficientBR + muVaryR;
+				}
+				muCoefficientBR = muCoefficientBR / muAttackR;
+			}
+			else
+			{
+				muCoefficientBR = muCoefficientBR * ((muSpeedBR * muSpeedBR)-1.0);
+				muCoefficientBR = muCoefficientBR + 1.0;
+				muCoefficientBR = muCoefficientBR / (muSpeedBR * muSpeedBR);
+			}
+			muNewSpeedR = muSpeedBR * (muSpeedBR-1);
+			muNewSpeedR = muNewSpeedR + fabs(squaredSampleR*release)+fastest;
+			muSpeedBR = muNewSpeedR / muSpeedBR;
+		}
+		//got coefficients, adjusted speeds for R
+				
+		if (flip)
+		{
+			coefficient = (muCoefficientAL + pow(muCoefficientAL,2))/2.0;
+			inputSampleL *= coefficient;
+			coefficient = (muCoefficientAR + pow(muCoefficientAR,2))/2.0;
+			inputSampleR *= coefficient;
+		}
+		else
+		{
+			coefficient = (muCoefficientBL + pow(muCoefficientBL,2))/2.0;
+			inputSampleL *= coefficient;
+			coefficient = (muCoefficientBR + pow(muCoefficientBR,2))/2.0;
+			inputSampleR *= coefficient;
+		}
+		//applied compression with vari-vari-µ-µ-µ-µ-µ-µ-is-the-kitten-song o/~
+		//applied gain correction to control output level- tends to constrain sound rather than inflate it
+		flip = !flip;		
+
+		if (output < 1.0) {
+			inputSampleL *= output;
+			inputSampleR *= output;
+		}
+		if (wet < 1.0) {
+			inputSampleL = (drySampleL * (1.0-wet)) + (inputSampleL * wet);
+			inputSampleR = (drySampleR * (1.0-wet)) + (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
+		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 VariMu::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 = 2.0;
+	overallscale /= 44100.0;
+	overallscale *= getSampleRate();
+	
+	double threshold = 1.001 - (1.0-pow(1.0-A,3));
+	double muMakeupGain = sqrt(1.0 / threshold);
+	muMakeupGain = (muMakeupGain + sqrt(muMakeupGain))/2.0;
+	muMakeupGain = sqrt(muMakeupGain);
+	double outGain = sqrt(muMakeupGain);
+	//gain settings around threshold
+	double release = pow((1.15-B),5)*32768.0;
+	release /= overallscale;
+	double fastest = sqrt(release);
+	//speed settings around release
+	double coefficient;
+	double output = outGain * C;
+	double wet = D;
+	long double squaredSampleL;
+	long double squaredSampleR;
+	
+	// µ µ µ µ µ µ µ µ µ µ µ µ is the kitten song o/~
+    
+    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 aVariMu. We want a 'air' hiss
+		long double drySampleL = inputSampleL;
+		long double drySampleR = inputSampleR;
+		
+		if (fabs(inputSampleL) > fabs(previousL)) squaredSampleL = previousL * previousL;
+		else squaredSampleL = inputSampleL * inputSampleL;
+		previousL = inputSampleL;
+		inputSampleL *= muMakeupGain;
+		
+		if (fabs(inputSampleR) > fabs(previousR)) squaredSampleR = previousR * previousR;
+		else squaredSampleR = inputSampleR * inputSampleR;
+		previousR = inputSampleR;
+		inputSampleR *= muMakeupGain;
+		
+		//adjust coefficients for L
+		if (flip)
+		{
+			if (fabs(squaredSampleL) > threshold)
+			{
+				muVaryL = threshold / fabs(squaredSampleL);
+				muAttackL = sqrt(fabs(muSpeedAL));
+				muCoefficientAL = muCoefficientAL * (muAttackL-1.0);
+				if (muVaryL < threshold)
+				{
+					muCoefficientAL = muCoefficientAL + threshold;
+				}
+				else
+				{
+					muCoefficientAL = muCoefficientAL + muVaryL;
+				}
+				muCoefficientAL = muCoefficientAL / muAttackL;
+			}
+			else
+			{
+				muCoefficientAL = muCoefficientAL * ((muSpeedAL * muSpeedAL)-1.0);
+				muCoefficientAL = muCoefficientAL + 1.0;
+				muCoefficientAL = muCoefficientAL / (muSpeedAL * muSpeedAL);
+			}
+			muNewSpeedL = muSpeedAL * (muSpeedAL-1);
+			muNewSpeedL = muNewSpeedL + fabs(squaredSampleL*release)+fastest;
+			muSpeedAL = muNewSpeedL / muSpeedAL;
+		}
+		else
+		{
+			if (fabs(squaredSampleL) > threshold)
+			{
+				muVaryL = threshold / fabs(squaredSampleL);
+				muAttackL = sqrt(fabs(muSpeedBL));
+				muCoefficientBL = muCoefficientBL * (muAttackL-1);
+				if (muVaryL < threshold)
+				{
+					muCoefficientBL = muCoefficientBL + threshold;
+				}
+				else
+				{
+					muCoefficientBL = muCoefficientBL + muVaryL;
+				}
+				muCoefficientBL = muCoefficientBL / muAttackL;
+			}
+			else
+			{
+				muCoefficientBL = muCoefficientBL * ((muSpeedBL * muSpeedBL)-1.0);
+				muCoefficientBL = muCoefficientBL + 1.0;
+				muCoefficientBL = muCoefficientBL / (muSpeedBL * muSpeedBL);
+			}
+			muNewSpeedL = muSpeedBL * (muSpeedBL-1);
+			muNewSpeedL = muNewSpeedL + fabs(squaredSampleL*release)+fastest;
+			muSpeedBL = muNewSpeedL / muSpeedBL;
+		}
+		//got coefficients, adjusted speeds for L
+		
+		//adjust coefficients for R
+		if (flip)
+		{
+			if (fabs(squaredSampleR) > threshold)
+			{
+				muVaryR = threshold / fabs(squaredSampleR);
+				muAttackR = sqrt(fabs(muSpeedAR));
+				muCoefficientAR = muCoefficientAR * (muAttackR-1.0);
+				if (muVaryR < threshold)
+				{
+					muCoefficientAR = muCoefficientAR + threshold;
+				}
+				else
+				{
+					muCoefficientAR = muCoefficientAR + muVaryR;
+				}
+				muCoefficientAR = muCoefficientAR / muAttackR;
+			}
+			else
+			{
+				muCoefficientAR = muCoefficientAR * ((muSpeedAR * muSpeedAR)-1.0);
+				muCoefficientAR = muCoefficientAR + 1.0;
+				muCoefficientAR = muCoefficientAR / (muSpeedAR * muSpeedAR);
+			}
+			muNewSpeedR = muSpeedAR * (muSpeedAR-1);
+			muNewSpeedR = muNewSpeedR + fabs(squaredSampleR*release)+fastest;
+			muSpeedAR = muNewSpeedR / muSpeedAR;
+		}
+		else
+		{
+			if (fabs(squaredSampleR) > threshold)
+			{
+				muVaryR = threshold / fabs(squaredSampleR);
+				muAttackR = sqrt(fabs(muSpeedBR));
+				muCoefficientBR = muCoefficientBR * (muAttackR-1);
+				if (muVaryR < threshold)
+				{
+					muCoefficientBR = muCoefficientBR + threshold;
+				}
+				else
+				{
+					muCoefficientBR = muCoefficientBR + muVaryR;
+				}
+				muCoefficientBR = muCoefficientBR / muAttackR;
+			}
+			else
+			{
+				muCoefficientBR = muCoefficientBR * ((muSpeedBR * muSpeedBR)-1.0);
+				muCoefficientBR = muCoefficientBR + 1.0;
+				muCoefficientBR = muCoefficientBR / (muSpeedBR * muSpeedBR);
+			}
+			muNewSpeedR = muSpeedBR * (muSpeedBR-1);
+			muNewSpeedR = muNewSpeedR + fabs(squaredSampleR*release)+fastest;
+			muSpeedBR = muNewSpeedR / muSpeedBR;
+		}
+		//got coefficients, adjusted speeds for R
+		
+		if (flip)
+		{
+			coefficient = (muCoefficientAL + pow(muCoefficientAL,2))/2.0;
+			inputSampleL *= coefficient;
+			coefficient = (muCoefficientAR + pow(muCoefficientAR,2))/2.0;
+			inputSampleR *= coefficient;
+		}
+		else
+		{
+			coefficient = (muCoefficientBL + pow(muCoefficientBL,2))/2.0;
+			inputSampleL *= coefficient;
+			coefficient = (muCoefficientBR + pow(muCoefficientBR,2))/2.0;
+			inputSampleR *= coefficient;
+		}
+		//applied compression with vari-vari-µ-µ-µ-µ-µ-µ-is-the-kitten-song o/~
+		//applied gain correction to control output level- tends to constrain sound rather than inflate it
+		flip = !flip;		
+		
+		if (output < 1.0) {
+			inputSampleL *= output;
+			inputSampleR *= output;
+		}
+		if (wet < 1.0) {
+			inputSampleL = (drySampleL * (1.0-wet)) + (inputSampleL * wet);
+			inputSampleR = (drySampleR * (1.0-wet)) + (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
+		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.	
+}