StudioTan

3 files changed, 789 insertions, 0 deletions

diff --git a/plugins/MacVST/StudioTan/source/StudioTan.cpp b/plugins/MacVST/StudioTan/source/StudioTan.cpp
new file mode 100755
index 0000000..31242f7
--- /dev/null
+++ b/plugins/MacVST/StudioTan/source/StudioTan.cpp
@@ -0,0 +1,156 @@
+/* ========================================
+ *  StudioTan - StudioTan.h
+ *  Copyright (c) 2016 airwindows, All rights reserved
+ * ======================================== */
+
+#ifndef __StudioTan_H
+#include "StudioTan.h"
+#endif
+
+AudioEffect* createEffectInstance(audioMasterCallback audioMaster) {return new StudioTan(audioMaster);}
+
+StudioTan::StudioTan(audioMasterCallback audioMaster) :
+    AudioEffectX(audioMaster, kNumPrograms, kNumParameters)
+{
+	A = 0.0;
+
+	bynL[0] = 1000.0;
+	bynL[1] = 301.0;
+	bynL[2] = 176.0;
+	bynL[3] = 125.0;
+	bynL[4] = 97.0;
+	bynL[5] = 79.0;
+	bynL[6] = 67.0;
+	bynL[7] = 58.0;
+	bynL[8] = 51.0;
+	bynL[9] = 46.0;
+	bynL[10] = 1000.0;
+	noiseShapingL = 0.0;
+	lastSampleL = 0.0;
+	lastSample2L = 0.0;
+
+	bynR[0] = 1000.0;
+	bynR[1] = 301.0;
+	bynR[2] = 176.0;
+	bynR[3] = 125.0;
+	bynR[4] = 97.0;
+	bynR[5] = 79.0;
+	bynR[6] = 67.0;
+	bynR[7] = 58.0;
+	bynR[8] = 51.0;
+	bynR[9] = 46.0;
+	bynR[10] = 1000.0;
+	noiseShapingR = 0.0;
+	lastSampleR = 0.0;
+	lastSample2R = 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
+}
+
+StudioTan::~StudioTan() {}
+VstInt32 StudioTan::getVendorVersion () {return 1000;}
+void StudioTan::setProgramName(char *name) {vst_strncpy (_programName, name, kVstMaxProgNameLen);}
+void StudioTan::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 StudioTan::getChunk (void** data, bool isPreset)
+{
+	float *chunkData = (float *)calloc(kNumParameters, sizeof(float));
+	chunkData[0] = A;
+	/* 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 StudioTan::setChunk (void* data, VstInt32 byteSize, bool isPreset)
+{	
+	float *chunkData = (float *)data;
+	A = pinParameter(chunkData[0]);
+	/* 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 StudioTan::setParameter(VstInt32 index, float value) {
+    switch (index) {
+        case kParamA: A = value; break;
+        default: throw; // unknown parameter, shouldn't happen!
+    }
+}
+
+float StudioTan::getParameter(VstInt32 index) {
+    switch (index) {
+        case kParamA: return A; 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 StudioTan::getParameterName(VstInt32 index, char *text) {
+    switch (index) {
+        case kParamA: vst_strncpy (text, "Quantzr", kVstMaxParamStrLen); break;
+        default: break; // unknown parameter, shouldn't happen!
+    } //this is our labels for displaying in the VST host
+}
+
+void StudioTan::getParameterDisplay(VstInt32 index, char *text) {
+    switch (index) {
+        case kParamA: switch((VstInt32)( A * 5.999 )) //0 to almost edge of # of params
+		{case 0: vst_strncpy (text, "ST 24", kVstMaxParamStrLen); break;
+		 case 1: vst_strncpy (text, "DMT 24", kVstMaxParamStrLen); break;
+		 case 2: vst_strncpy (text, "NJAD 24", kVstMaxParamStrLen); break;
+		 case 3: vst_strncpy (text, "ST 16", kVstMaxParamStrLen); break;
+		 case 4: vst_strncpy (text, "DMT 16", kVstMaxParamStrLen); break;
+		 case 5: vst_strncpy (text, "NJAD 16", kVstMaxParamStrLen); break;
+		 default: break; // unknown parameter, shouldn't happen!
+		} break; //E as example 'popup' parameter with four values  */
+			
+        default: break; // unknown parameter, shouldn't happen!
+	} //this displays the values and handles 'popups' where it's discrete choices
+}
+
+void StudioTan::getParameterLabel(VstInt32 index, char *text) {
+    switch (index) {
+        case kParamA: vst_strncpy (text, "", kVstMaxParamStrLen); break;
+		default: break; // unknown parameter, shouldn't happen!
+    }
+}
+
+VstInt32 StudioTan::canDo(char *text) 
+{ return (_canDo.find(text) == _canDo.end()) ? -1: 1; } // 1 = yes, -1 = no, 0 = don't know
+
+bool StudioTan::getEffectName(char* name) {
+    vst_strncpy(name, "StudioTan", kVstMaxProductStrLen); return true;
+}
+
+VstPlugCategory StudioTan::getPlugCategory() {return kPlugCategEffect;}
+
+bool StudioTan::getProductString(char* text) {
+  	vst_strncpy (text, "airwindows StudioTan", kVstMaxProductStrLen); return true;
+}
+
+bool StudioTan::getVendorString(char* text) {
+  	vst_strncpy (text, "airwindows", kVstMaxVendorStrLen); return true;
+}
diff --git a/plugins/MacVST/StudioTan/source/StudioTan.h b/plugins/MacVST/StudioTan/source/StudioTan.h
new file mode 100755
index 0000000..66caf48
--- /dev/null
+++ b/plugins/MacVST/StudioTan/source/StudioTan.h
@@ -0,0 +1,69 @@
+/* ========================================
+ *  StudioTan - StudioTan.h
+ *  Created 8/12/11 by SPIAdmin 
+ *  Copyright (c) 2011 __MyCompanyName__, All rights reserved
+ * ======================================== */
+
+#ifndef __StudioTan_H
+#define __StudioTan_H
+
+#ifndef __audioeffect__
+#include "audioeffectx.h"
+#endif
+
+#include <set>
+#include <string>
+#include <math.h>
+
+enum {
+	kParamA = 0,
+  kNumParameters = 1
+}; //
+
+const int kNumPrograms = 0;
+const int kNumInputs = 2;
+const int kNumOutputs = 2;
+const unsigned long kUniqueId = 'stan';    //Change this to what the AU identity is!
+
+class StudioTan : 
+    public AudioEffectX 
+{
+public:
+    StudioTan(audioMasterCallback audioMaster);
+    ~StudioTan();
+    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;
+    
+	long double bynL[13];
+	long double noiseShapingL;
+	long double lastSampleL;
+	long double lastSample2L;
+
+	long double bynR[13];
+	long double noiseShapingR;
+	long double lastSampleR;
+	long double lastSample2R;
+	
+    float A;
+
+};
+
+#endif
diff --git a/plugins/MacVST/StudioTan/source/StudioTanProc.cpp b/plugins/MacVST/StudioTan/source/StudioTanProc.cpp
new file mode 100755
index 0000000..3258f7b
--- /dev/null
+++ b/plugins/MacVST/StudioTan/source/StudioTanProc.cpp
@@ -0,0 +1,564 @@
+/* ========================================
+ *  StudioTan - StudioTan.h
+ *  Copyright (c) 2016 airwindows, All rights reserved
+ * ======================================== */
+
+#ifndef __StudioTan_H
+#include "StudioTan.h"
+#endif
+
+void StudioTan::processReplacing(float **inputs, float **outputs, VstInt32 sampleFrames) 
+{
+    float* in1  =  inputs[0];
+    float* in2  =  inputs[1];
+    float* out1 = outputs[0];
+    float* out2 = outputs[1];
+	
+	bool highres = true; //for 24 bit: false for 16 bit
+	bool brightfloor = true; //for Studio Tan: false for Dither Me Timbers
+	bool benford = true; //for Not Just Another Dither: false for newer two
+	bool cutbins = false; //for NJAD: only attenuate bins if one gets very full
+	switch ((VstInt32)( A * 5.999 ))
+	{
+		case 0: benford = false; break; //Studio Tan 24
+		case 1: benford = false; brightfloor = false; break; //Dither Me Timbers 24
+		case 2: break; //Not Just Another Dither 24
+		case 3: benford = false; highres = false; break; //Studio Tan 16
+		case 4: benford = false; brightfloor = false; highres = false; break; //Dither Me Timbers 16
+		case 5: highres = false; break; //Not Just Another Dither 16
+	}
+	
+    while (--sampleFrames >= 0)
+    {
+		long double inputSampleL;
+		long double outputSampleL;
+		long double drySampleL;
+		long double inputSampleR;
+		long double outputSampleR;
+		long double drySampleR;
+		
+		if (highres) {
+			inputSampleL = *in1 * 8388608.0;
+			inputSampleR = *in2 * 8388608.0;
+		} else {
+			inputSampleL = *in1 * 32768.0;
+			inputSampleR = *in2 * 32768.0;
+		}
+		//shared input stage
+		
+		if (benford) {
+			//begin Not Just Another Dither
+			drySampleL = inputSampleL;
+			drySampleR = inputSampleR;
+			inputSampleL -= noiseShapingL;
+			inputSampleR -= noiseShapingR;
+			
+			cutbins = false;
+			long double benfordize; //we get to re-use this for each channel
+			
+			//begin left channel NJAD
+			benfordize = floor(inputSampleL);
+			while (benfordize >= 1.0) {benfordize /= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			int hotbinA = floor(benfordize);
+			//hotbin becomes the Benford bin value for this number floored
+			long double totalA = 0;
+			if ((hotbinA > 0) && (hotbinA < 10))
+			{
+				bynL[hotbinA] += 1;
+				if (bynL[hotbinA] > 982) cutbins = true;
+				totalA += (301-bynL[1]);
+				totalA += (176-bynL[2]);
+				totalA += (125-bynL[3]);
+				totalA += (97-bynL[4]);
+				totalA += (79-bynL[5]);
+				totalA += (67-bynL[6]);
+				totalA += (58-bynL[7]);
+				totalA += (51-bynL[8]);
+				totalA += (46-bynL[9]);
+				bynL[hotbinA] -= 1;
+			} else {hotbinA = 10;}
+			//produce total number- smaller is closer to Benford real
+			
+			benfordize = ceil(inputSampleL);
+			while (benfordize >= 1.0) {benfordize /= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			int hotbinB = floor(benfordize);
+			//hotbin becomes the Benford bin value for this number ceiled
+			long double totalB = 0;
+			if ((hotbinB > 0) && (hotbinB < 10))
+			{
+				bynL[hotbinB] += 1;
+				if (bynL[hotbinB] > 982) cutbins = true;
+				totalB += (301-bynL[1]);
+				totalB += (176-bynL[2]);
+				totalB += (125-bynL[3]);
+				totalB += (97-bynL[4]);
+				totalB += (79-bynL[5]);
+				totalB += (67-bynL[6]);
+				totalB += (58-bynL[7]);
+				totalB += (51-bynL[8]);
+				totalB += (46-bynL[9]);
+				bynL[hotbinB] -= 1;
+			} else {hotbinB = 10;}
+			//produce total number- smaller is closer to Benford real
+			
+			if (totalA < totalB)
+			{
+				bynL[hotbinA] += 1;
+				outputSampleL = floor(inputSampleL);
+			}
+			else
+			{
+				bynL[hotbinB] += 1;
+				outputSampleL = floor(inputSampleL+1);
+			}
+			//assign the relevant one to the delay line
+			//and floor/ceil signal accordingly
+			if (cutbins) {
+				bynL[1] *= 0.99;
+				bynL[2] *= 0.99;
+				bynL[3] *= 0.99;
+				bynL[4] *= 0.99;
+				bynL[5] *= 0.99;
+				bynL[6] *= 0.99;
+				bynL[7] *= 0.99;
+				bynL[8] *= 0.99;
+				bynL[9] *= 0.99;
+				bynL[10] *= 0.99; //catchall for garbage data
+			}
+			noiseShapingL += outputSampleL - drySampleL;
+			//end left channel NJAD
+			
+			//begin right channel NJAD
+			cutbins = false;
+			benfordize = floor(inputSampleR);
+			while (benfordize >= 1.0) {benfordize /= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			hotbinA = floor(benfordize);
+			//hotbin becomes the Benford bin value for this number floored
+			totalA = 0;
+			if ((hotbinA > 0) && (hotbinA < 10))
+			{
+				bynR[hotbinA] += 1;
+				if (bynR[hotbinA] > 982) cutbins = true;
+				totalA += (301-bynR[1]);
+				totalA += (176-bynR[2]);
+				totalA += (125-bynR[3]);
+				totalA += (97-bynR[4]);
+				totalA += (79-bynR[5]);
+				totalA += (67-bynR[6]);
+				totalA += (58-bynR[7]);
+				totalA += (51-bynR[8]);
+				totalA += (46-bynR[9]);
+				bynR[hotbinA] -= 1;
+			} else {hotbinA = 10;}
+			//produce total number- smaller is closer to Benford real
+			
+			benfordize = ceil(inputSampleR);
+			while (benfordize >= 1.0) {benfordize /= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			hotbinB = floor(benfordize);
+			//hotbin becomes the Benford bin value for this number ceiled
+			totalB = 0;
+			if ((hotbinB > 0) && (hotbinB < 10))
+			{
+				bynR[hotbinB] += 1;
+				if (bynR[hotbinB] > 982) cutbins = true;
+				totalB += (301-bynR[1]);
+				totalB += (176-bynR[2]);
+				totalB += (125-bynR[3]);
+				totalB += (97-bynR[4]);
+				totalB += (79-bynR[5]);
+				totalB += (67-bynR[6]);
+				totalB += (58-bynR[7]);
+				totalB += (51-bynR[8]);
+				totalB += (46-bynR[9]);
+				bynR[hotbinB] -= 1;
+			} else {hotbinB = 10;}
+			//produce total number- smaller is closer to Benford real
+			
+			if (totalA < totalB)
+			{
+				bynR[hotbinA] += 1;
+				outputSampleR = floor(inputSampleR);
+			}
+			else
+			{
+				bynR[hotbinB] += 1;
+				outputSampleR = floor(inputSampleR+1);
+			}
+			//assign the relevant one to the delay line
+			//and floor/ceil signal accordingly
+			if (cutbins) {
+				bynR[1] *= 0.99;
+				bynR[2] *= 0.99;
+				bynR[3] *= 0.99;
+				bynR[4] *= 0.99;
+				bynR[5] *= 0.99;
+				bynR[6] *= 0.99;
+				bynR[7] *= 0.99;
+				bynR[8] *= 0.99;
+				bynR[9] *= 0.99;
+				bynR[10] *= 0.99; //catchall for garbage data
+			}
+			noiseShapingR += outputSampleR - drySampleR;
+			//end right channel NJAD			
+			
+			//end Not Just Another Dither
+		} else {		
+			//begin StudioTan or Dither Me Timbers
+			if (brightfloor) {
+				lastSampleL -= (noiseShapingL*0.8);
+				lastSampleR -= (noiseShapingR*0.8);
+				if ((lastSampleL+lastSampleL) <= (inputSampleL+lastSample2L)) outputSampleL = floor(lastSampleL); //StudioTan
+				else outputSampleL = floor(lastSampleL+1.0); //round down or up based on whether it softens treble angles
+				if ((lastSampleR+lastSampleR) <= (inputSampleR+lastSample2R)) outputSampleR = floor(lastSampleR); //StudioTan
+				else outputSampleR = floor(lastSampleR+1.0); //round down or up based on whether it softens treble angles
+			} else {
+				lastSampleL -= (noiseShapingL*0.11);
+				lastSampleR -= (noiseShapingR*0.11);
+				if ((lastSampleL+lastSampleL) >= (inputSampleL+lastSample2L)) outputSampleL = floor(lastSampleL); //DitherMeTimbers
+				else outputSampleL = floor(lastSampleL+1.0); //round down or up based on whether it softens treble angles
+				if ((lastSampleR+lastSampleR) >= (inputSampleR+lastSample2R)) outputSampleR = floor(lastSampleR); //DitherMeTimbers
+				else outputSampleR = floor(lastSampleR+1.0); //round down or up based on whether it softens treble angles
+			}
+			noiseShapingL += outputSampleL;
+			noiseShapingL -= lastSampleL; //apply noise shaping
+			lastSample2L = lastSampleL;
+			lastSampleL = inputSampleL; //we retain three samples in a row
+			
+			noiseShapingR += outputSampleR;
+			noiseShapingR -= lastSampleR; //apply noise shaping
+			lastSample2R = lastSampleR;
+			lastSampleR = inputSampleR; //we retain three samples in a row
+			//end StudioTan or Dither Me Timbers
+		}
+		
+		//shared output stage
+		long double noiseSuppressL = fabs(inputSampleL);
+		if (noiseShapingL > noiseSuppressL) noiseShapingL = noiseSuppressL;
+		if (noiseShapingL < -noiseSuppressL) noiseShapingL = -noiseSuppressL;
+
+		long double noiseSuppressR = fabs(inputSampleR);
+		if (noiseShapingR > noiseSuppressR) noiseShapingR = noiseSuppressR;
+		if (noiseShapingR < -noiseSuppressR) noiseShapingR = -noiseSuppressR;
+		
+		float ironBarL;
+		float ironBarR;
+		if (highres) {
+			ironBarL = outputSampleL / 8388608.0;
+			ironBarR = outputSampleR / 8388608.0;
+		} else {
+			ironBarL = outputSampleL / 32768.0;
+			ironBarR = outputSampleR / 32768.0;
+		}
+		
+		if (ironBarL > 1.0) ironBarL = 1.0;
+		if (ironBarL < -1.0) ironBarL = -1.0;
+		if (ironBarR > 1.0) ironBarR = 1.0;
+		if (ironBarR < -1.0) ironBarR = -1.0;
+		
+		*out1 = ironBarL;
+		*out2 = ironBarR;
+
+		*in1++;
+		*in2++;
+		*out1++;
+		*out2++;
+    }
+}
+
+void StudioTan::processDoubleReplacing(double **inputs, double **outputs, VstInt32 sampleFrames) 
+{
+    double* in1  =  inputs[0];
+    double* in2  =  inputs[1];
+    double* out1 = outputs[0];
+    double* out2 = outputs[1];
+
+	bool highres = true; //for 24 bit: false for 16 bit
+	bool brightfloor = true; //for Studio Tan: false for Dither Me Timbers
+	bool benford = true; //for Not Just Another Dither: false for newer two
+	bool cutbins = false; //for NJAD: only attenuate bins if one gets very full
+	switch ((VstInt32)( A * 5.999 ))
+	{
+		case 0: benford = false; break; //Studio Tan 24
+		case 1: benford = false; brightfloor = false; break; //Dither Me Timbers 24
+		case 2: break; //Not Just Another Dither 24
+		case 3: benford = false; highres = false; break; //Studio Tan 16
+		case 4: benford = false; brightfloor = false; highres = false; break; //Dither Me Timbers 16
+		case 5: highres = false; break; //Not Just Another Dither 16
+	}
+		
+    while (--sampleFrames >= 0)
+    {
+		long double inputSampleL;
+		long double outputSampleL;
+		long double drySampleL;
+		long double inputSampleR;
+		long double outputSampleR;
+		long double drySampleR;
+		
+		if (highres) {
+			inputSampleL = *in1 * 8388608.0;
+			inputSampleR = *in2 * 8388608.0;
+		} else {
+			inputSampleL = *in1 * 32768.0;
+			inputSampleR = *in2 * 32768.0;
+		}
+		//shared input stage
+		
+		if (benford) {
+			//begin Not Just Another Dither
+			drySampleL = inputSampleL;
+			drySampleR = inputSampleR;
+			inputSampleL -= noiseShapingL;
+			inputSampleR -= noiseShapingR;
+			
+			cutbins = false;
+			long double benfordize; //we get to re-use this for each channel
+			
+			//begin left channel NJAD
+			benfordize = floor(inputSampleL);
+			while (benfordize >= 1.0) {benfordize /= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			int hotbinA = floor(benfordize);
+			//hotbin becomes the Benford bin value for this number floored
+			long double totalA = 0;
+			if ((hotbinA > 0) && (hotbinA < 10))
+			{
+				bynL[hotbinA] += 1;
+				if (bynL[hotbinA] > 982) cutbins = true;
+				totalA += (301-bynL[1]);
+				totalA += (176-bynL[2]);
+				totalA += (125-bynL[3]);
+				totalA += (97-bynL[4]);
+				totalA += (79-bynL[5]);
+				totalA += (67-bynL[6]);
+				totalA += (58-bynL[7]);
+				totalA += (51-bynL[8]);
+				totalA += (46-bynL[9]);
+				bynL[hotbinA] -= 1;
+			} else {hotbinA = 10;}
+			//produce total number- smaller is closer to Benford real
+			
+			benfordize = ceil(inputSampleL);
+			while (benfordize >= 1.0) {benfordize /= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			int hotbinB = floor(benfordize);
+			//hotbin becomes the Benford bin value for this number ceiled
+			long double totalB = 0;
+			if ((hotbinB > 0) && (hotbinB < 10))
+			{
+				bynL[hotbinB] += 1;
+				if (bynL[hotbinB] > 982) cutbins = true;
+				totalB += (301-bynL[1]);
+				totalB += (176-bynL[2]);
+				totalB += (125-bynL[3]);
+				totalB += (97-bynL[4]);
+				totalB += (79-bynL[5]);
+				totalB += (67-bynL[6]);
+				totalB += (58-bynL[7]);
+				totalB += (51-bynL[8]);
+				totalB += (46-bynL[9]);
+				bynL[hotbinB] -= 1;
+			} else {hotbinB = 10;}
+			//produce total number- smaller is closer to Benford real
+			
+			if (totalA < totalB)
+			{
+				bynL[hotbinA] += 1;
+				outputSampleL = floor(inputSampleL);
+			}
+			else
+			{
+				bynL[hotbinB] += 1;
+				outputSampleL = floor(inputSampleL+1);
+			}
+			//assign the relevant one to the delay line
+			//and floor/ceil signal accordingly
+			if (cutbins) {
+				bynL[1] *= 0.99;
+				bynL[2] *= 0.99;
+				bynL[3] *= 0.99;
+				bynL[4] *= 0.99;
+				bynL[5] *= 0.99;
+				bynL[6] *= 0.99;
+				bynL[7] *= 0.99;
+				bynL[8] *= 0.99;
+				bynL[9] *= 0.99;
+				bynL[10] *= 0.99; //catchall for garbage data
+			}
+			noiseShapingL += outputSampleL - drySampleL;
+			//end left channel NJAD
+			
+			//begin right channel NJAD
+			cutbins = false;
+			benfordize = floor(inputSampleR);
+			while (benfordize >= 1.0) {benfordize /= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			hotbinA = floor(benfordize);
+			//hotbin becomes the Benford bin value for this number floored
+			totalA = 0;
+			if ((hotbinA > 0) && (hotbinA < 10))
+			{
+				bynR[hotbinA] += 1;
+				if (bynR[hotbinA] > 982) cutbins = true;
+				totalA += (301-bynR[1]);
+				totalA += (176-bynR[2]);
+				totalA += (125-bynR[3]);
+				totalA += (97-bynR[4]);
+				totalA += (79-bynR[5]);
+				totalA += (67-bynR[6]);
+				totalA += (58-bynR[7]);
+				totalA += (51-bynR[8]);
+				totalA += (46-bynR[9]);
+				bynR[hotbinA] -= 1;
+			} else {hotbinA = 10;}
+			//produce total number- smaller is closer to Benford real
+			
+			benfordize = ceil(inputSampleR);
+			while (benfordize >= 1.0) {benfordize /= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			if (benfordize < 1.0) {benfordize *= 10;}
+			hotbinB = floor(benfordize);
+			//hotbin becomes the Benford bin value for this number ceiled
+			totalB = 0;
+			if ((hotbinB > 0) && (hotbinB < 10))
+			{
+				bynR[hotbinB] += 1;
+				if (bynR[hotbinB] > 982) cutbins = true;
+				totalB += (301-bynR[1]);
+				totalB += (176-bynR[2]);
+				totalB += (125-bynR[3]);
+				totalB += (97-bynR[4]);
+				totalB += (79-bynR[5]);
+				totalB += (67-bynR[6]);
+				totalB += (58-bynR[7]);
+				totalB += (51-bynR[8]);
+				totalB += (46-bynR[9]);
+				bynR[hotbinB] -= 1;
+			} else {hotbinB = 10;}
+			//produce total number- smaller is closer to Benford real
+			
+			if (totalA < totalB)
+			{
+				bynR[hotbinA] += 1;
+				outputSampleR = floor(inputSampleR);
+			}
+			else
+			{
+				bynR[hotbinB] += 1;
+				outputSampleR = floor(inputSampleR+1);
+			}
+			//assign the relevant one to the delay line
+			//and floor/ceil signal accordingly
+			if (cutbins) {
+				bynR[1] *= 0.99;
+				bynR[2] *= 0.99;
+				bynR[3] *= 0.99;
+				bynR[4] *= 0.99;
+				bynR[5] *= 0.99;
+				bynR[6] *= 0.99;
+				bynR[7] *= 0.99;
+				bynR[8] *= 0.99;
+				bynR[9] *= 0.99;
+				bynR[10] *= 0.99; //catchall for garbage data
+			}
+			noiseShapingR += outputSampleR - drySampleR;
+			//end right channel NJAD			
+			
+			//end Not Just Another Dither
+		} else {		
+			//begin StudioTan or Dither Me Timbers
+			if (brightfloor) {
+				lastSampleL -= (noiseShapingL*0.8);
+				lastSampleR -= (noiseShapingR*0.8);
+				if ((lastSampleL+lastSampleL) <= (inputSampleL+lastSample2L)) outputSampleL = floor(lastSampleL); //StudioTan
+				else outputSampleL = floor(lastSampleL+1.0); //round down or up based on whether it softens treble angles
+				if ((lastSampleR+lastSampleR) <= (inputSampleR+lastSample2R)) outputSampleR = floor(lastSampleR); //StudioTan
+				else outputSampleR = floor(lastSampleR+1.0); //round down or up based on whether it softens treble angles
+			} else {
+				lastSampleL -= (noiseShapingL*0.11);
+				lastSampleR -= (noiseShapingR*0.11);
+				if ((lastSampleL+lastSampleL) >= (inputSampleL+lastSample2L)) outputSampleL = floor(lastSampleL); //DitherMeTimbers
+				else outputSampleL = floor(lastSampleL+1.0); //round down or up based on whether it softens treble angles
+				if ((lastSampleR+lastSampleR) >= (inputSampleR+lastSample2R)) outputSampleR = floor(lastSampleR); //DitherMeTimbers
+				else outputSampleR = floor(lastSampleR+1.0); //round down or up based on whether it softens treble angles
+			}
+			noiseShapingL += outputSampleL;
+			noiseShapingL -= lastSampleL; //apply noise shaping
+			lastSample2L = lastSampleL;
+			lastSampleL = inputSampleL; //we retain three samples in a row
+			
+			noiseShapingR += outputSampleR;
+			noiseShapingR -= lastSampleR; //apply noise shaping
+			lastSample2R = lastSampleR;
+			lastSampleR = inputSampleR; //we retain three samples in a row
+			//end StudioTan or Dither Me Timbers
+		}
+		
+		//shared output stage
+		long double noiseSuppressL = fabs(inputSampleL);
+		if (noiseShapingL > noiseSuppressL) noiseShapingL = noiseSuppressL;
+		if (noiseShapingL < -noiseSuppressL) noiseShapingL = -noiseSuppressL;
+		
+		long double noiseSuppressR = fabs(inputSampleR);
+		if (noiseShapingR > noiseSuppressR) noiseShapingR = noiseSuppressR;
+		if (noiseShapingR < -noiseSuppressR) noiseShapingR = -noiseSuppressR;
+		
+		double ironBarL;
+		double ironBarR;
+		if (highres) {
+			ironBarL = outputSampleL / 8388608.0;
+			ironBarR = outputSampleR / 8388608.0;
+		} else {
+			ironBarL = outputSampleL / 32768.0;
+			ironBarR = outputSampleR / 32768.0;
+		}
+		
+		if (ironBarL > 1.0) ironBarL = 1.0;
+		if (ironBarL < -1.0) ironBarL = -1.0;
+		if (ironBarR > 1.0) ironBarR = 1.0;
+		if (ironBarR < -1.0) ironBarR = -1.0;
+		
+		*out1 = ironBarL;
+		*out2 = ironBarR;
+		
+		*in1++;
+		*in2++;
+		*out1++;
+		*out2++;
+    }
+}