Ditherbox… and individual dithers

3 files changed, 669 insertions, 0 deletions

diff --git a/plugins/MacVST/NotJustAnotherDither/source/NotJustAnotherDither.cpp b/plugins/MacVST/NotJustAnotherDither/source/NotJustAnotherDither.cpp
new file mode 100755
index 0000000..22a3cc7
--- /dev/null
+++ b/plugins/MacVST/NotJustAnotherDither/source/NotJustAnotherDither.cpp
@@ -0,0 +1,103 @@
+/* ========================================
+ *  NotJustAnotherDither - NotJustAnotherDither.h
+ *  Copyright (c) 2016 airwindows, All rights reserved
+ * ======================================== */
+
+#ifndef __NotJustAnotherDither_H
+#include "NotJustAnotherDither.h"
+#endif
+
+AudioEffect* createEffectInstance(audioMasterCallback audioMaster) {return new NotJustAnotherDither(audioMaster);}
+
+NotJustAnotherDither::NotJustAnotherDither(audioMasterCallback audioMaster) :
+    AudioEffectX(audioMaster, kNumPrograms, kNumParameters)
+{
+	bynL[0] = 1000;
+	bynL[1] = 301;
+	bynL[2] = 176;
+	bynL[3] = 125;
+	bynL[4] = 97;
+	bynL[5] = 79;
+	bynL[6] = 67;
+	bynL[7] = 58;
+	bynL[8] = 51;
+	bynL[9] = 46;
+	bynL[10] = 1000;
+	noiseShapingL = 0.0;
+	
+	bynR[0] = 1000;
+	bynR[1] = 301;
+	bynR[2] = 176;
+	bynR[3] = 125;
+	bynR[4] = 97;
+	bynR[5] = 79;
+	bynR[6] = 67;
+	bynR[7] = 58;
+	bynR[8] = 51;
+	bynR[9] = 46;
+	bynR[10] = 1000;
+	noiseShapingR = 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
+}
+
+NotJustAnotherDither::~NotJustAnotherDither() {}
+VstInt32 NotJustAnotherDither::getVendorVersion () {return 1000;}
+void NotJustAnotherDither::setProgramName(char *name) {vst_strncpy (_programName, name, kVstMaxProgNameLen);}
+void NotJustAnotherDither::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!
+
+
+VstInt32 NotJustAnotherDither::getChunk (void** data, bool isPreset)
+{
+	return kNumParameters * sizeof(float);
+}
+
+VstInt32 NotJustAnotherDither::setChunk (void* data, VstInt32 byteSize, bool isPreset)
+{	
+	return 0;
+}
+
+void NotJustAnotherDither::setParameter(VstInt32 index, float value) {
+ }
+
+float NotJustAnotherDither::getParameter(VstInt32 index) {
+    return 0.0; //we only need to update the relevant name, this is simple to manage
+}
+
+void NotJustAnotherDither::getParameterName(VstInt32 index, char *text) {
+}
+
+void NotJustAnotherDither::getParameterDisplay(VstInt32 index, char *text) {
+}
+
+void NotJustAnotherDither::getParameterLabel(VstInt32 index, char *text) {
+}
+
+VstInt32 NotJustAnotherDither::canDo(char *text) 
+{ return (_canDo.find(text) == _canDo.end()) ? -1: 1; } // 1 = yes, -1 = no, 0 = don't know
+
+bool NotJustAnotherDither::getEffectName(char* name) {
+    vst_strncpy(name, "NotJustAnotherDither", kVstMaxProductStrLen); return true;
+}
+
+VstPlugCategory NotJustAnotherDither::getPlugCategory() {return kPlugCategEffect;}
+
+bool NotJustAnotherDither::getProductString(char* text) {
+  	vst_strncpy (text, "airwindows NotJustAnotherDither", kVstMaxProductStrLen); return true;
+}
+
+bool NotJustAnotherDither::getVendorString(char* text) {
+  	vst_strncpy (text, "airwindows", kVstMaxVendorStrLen); return true;
+}
diff --git a/plugins/MacVST/NotJustAnotherDither/source/NotJustAnotherDither.h b/plugins/MacVST/NotJustAnotherDither/source/NotJustAnotherDither.h
new file mode 100755
index 0000000..37abf4c
--- /dev/null
+++ b/plugins/MacVST/NotJustAnotherDither/source/NotJustAnotherDither.h
@@ -0,0 +1,61 @@
+/* ========================================
+ *  NotJustAnotherDither - NotJustAnotherDither.h
+ *  Created 8/12/11 by SPIAdmin 
+ *  Copyright (c) 2011 __MyCompanyName__, All rights reserved
+ * ======================================== */
+
+#ifndef __NotJustAnotherDither_H
+#define __NotJustAnotherDither_H
+
+#ifndef __audioeffect__
+#include "audioeffectx.h"
+#endif
+
+#include <set>
+#include <string>
+#include <math.h>
+
+enum {
+  kNumParameters = 0
+}; //
+
+const int kNumPrograms = 0;
+const int kNumInputs = 2;
+const int kNumOutputs = 2;
+const unsigned long kUniqueId = 'njad';    //Change this to what the AU identity is!
+
+class NotJustAnotherDither : 
+    public AudioEffectX 
+{
+public:
+    NotJustAnotherDither(audioMasterCallback audioMaster);
+    ~NotJustAnotherDither();
+    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 bynL[13];
+	double bynR[13];
+	long double noiseShapingL;
+	long double noiseShapingR;
+
+};
+
+#endif
diff --git a/plugins/MacVST/NotJustAnotherDither/source/NotJustAnotherDitherProc.cpp b/plugins/MacVST/NotJustAnotherDither/source/NotJustAnotherDitherProc.cpp
new file mode 100755
index 0000000..46a073e
--- /dev/null
+++ b/plugins/MacVST/NotJustAnotherDither/source/NotJustAnotherDitherProc.cpp
@@ -0,0 +1,505 @@
+/* ========================================
+ *  NotJustAnotherDither - NotJustAnotherDither.h
+ *  Copyright (c) 2016 airwindows, All rights reserved
+ * ======================================== */
+
+#ifndef __NotJustAnotherDither_H
+#include "NotJustAnotherDither.h"
+#endif
+
+void NotJustAnotherDither::processReplacing(float **inputs, float **outputs, VstInt32 sampleFrames) 
+{
+    float* in1  =  inputs[0];
+    float* in2  =  inputs[1];
+    float* out1 = outputs[0];
+    float* out2 = outputs[1];
+	
+	long double inputSampleL;
+	long double inputSampleR;
+	
+	double benfordize;
+	int hotbinA;
+	int hotbinB;
+	double totalA;
+	double totalB;
+	float drySampleL;
+	float 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;
+
+		inputSampleL -= noiseShapingL;
+		inputSampleR -= noiseShapingR;
+
+		inputSampleL *= 8388608.0;
+		inputSampleR *= 8388608.0;
+		//0-1 is now one bit, now we dither
+				
+		//begin L
+		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;}
+		hotbinA = floor(benfordize);
+		//hotbin becomes the Benford bin value for this number floored
+		totalA = 0;
+		if ((hotbinA > 0) && (hotbinA < 10))
+		{
+			bynL[hotbinA] += 1;
+			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;}
+		hotbinB = floor(benfordize);
+		//hotbin becomes the Benford bin value for this number ceiled
+		totalB = 0;
+		if ((hotbinB > 0) && (hotbinB < 10))
+		{
+			bynL[hotbinB] += 1;
+			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;
+			inputSampleL = floor(inputSampleL);
+		}
+		else
+		{
+			bynL[hotbinB] += 1;
+			inputSampleL = ceil(inputSampleL);
+		}
+		//assign the relevant one to the delay line
+		//and floor/ceil signal accordingly
+		
+		totalA = bynL[1] + bynL[2] + bynL[3] + bynL[4] + bynL[5] + bynL[6] + bynL[7] + bynL[8] + bynL[9];
+		totalA /= 1000;
+		if (totalA = 0) totalA = 1;
+		bynL[1] /= totalA;
+		bynL[2] /= totalA;
+		bynL[3] /= totalA;
+		bynL[4] /= totalA;
+		bynL[5] /= totalA;
+		bynL[6] /= totalA;
+		bynL[7] /= totalA;
+		bynL[8] /= totalA;
+		bynL[9] /= totalA;
+		bynL[10] /= 2; //catchall for garbage data
+		//end L
+		
+		//begin R
+		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;
+			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;
+			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;
+			inputSampleR = floor(inputSampleR);
+		}
+		else
+		{
+			bynR[hotbinB] += 1;
+			inputSampleR = ceil(inputSampleR);
+		}
+		//assign the relevant one to the delay line
+		//and floor/ceil signal accordingly
+		
+		totalA = bynR[1] + bynR[2] + bynR[3] + bynR[4] + bynR[5] + bynR[6] + bynR[7] + bynR[8] + bynR[9];
+		totalA /= 1000;
+		if (totalA = 0) totalA = 1;
+		bynR[1] /= totalA;
+		bynR[2] /= totalA;
+		bynR[3] /= totalA;
+		bynR[4] /= totalA;
+		bynR[5] /= totalA;
+		bynR[6] /= totalA;
+		bynR[7] /= totalA;
+		bynR[8] /= totalA;
+		bynR[9] /= totalA;
+		bynR[10] /= 2; //catchall for garbage data
+		//end R
+		
+		inputSampleL /= 8388608.0;
+		inputSampleR /= 8388608.0;
+
+		noiseShapingL += inputSampleL - drySampleL;
+		noiseShapingR += inputSampleR - drySampleR;
+
+		*out1 = inputSampleL;
+		*out2 = inputSampleR;
+		
+		*in1++;
+		*in2++;
+		*out1++;
+		*out2++;
+    }
+}
+
+void NotJustAnotherDither::processDoubleReplacing(double **inputs, double **outputs, VstInt32 sampleFrames) 
+{
+    double* in1  =  inputs[0];
+    double* in2  =  inputs[1];
+    double* out1 = outputs[0];
+    double* out2 = outputs[1];
+	
+	
+	long double inputSampleL;
+	long double inputSampleR;
+	
+	double benfordize;
+	int hotbinA;
+	int hotbinB;
+	double totalA;
+	double totalB;
+	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;
+
+		inputSampleL -= noiseShapingL;
+		inputSampleR -= noiseShapingR;
+		
+		inputSampleL *= 8388608.0;
+		inputSampleR *= 8388608.0;
+		//0-1 is now one bit, now we dither
+		
+		//begin L
+		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;}
+		hotbinA = floor(benfordize);
+		//hotbin becomes the Benford bin value for this number floored
+		totalA = 0;
+		if ((hotbinA > 0) && (hotbinA < 10))
+		{
+			bynL[hotbinA] += 1;
+			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;}
+		hotbinB = floor(benfordize);
+		//hotbin becomes the Benford bin value for this number ceiled
+		totalB = 0;
+		if ((hotbinB > 0) && (hotbinB < 10))
+		{
+			bynL[hotbinB] += 1;
+			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;
+			inputSampleL = floor(inputSampleL);
+		}
+		else
+		{
+			bynL[hotbinB] += 1;
+			inputSampleL = ceil(inputSampleL);
+		}
+		//assign the relevant one to the delay line
+		//and floor/ceil signal accordingly
+		
+		totalA = bynL[1] + bynL[2] + bynL[3] + bynL[4] + bynL[5] + bynL[6] + bynL[7] + bynL[8] + bynL[9];
+		totalA /= 1000;
+		if (totalA = 0) totalA = 1;
+		bynL[1] /= totalA;
+		bynL[2] /= totalA;
+		bynL[3] /= totalA;
+		bynL[4] /= totalA;
+		bynL[5] /= totalA;
+		bynL[6] /= totalA;
+		bynL[7] /= totalA;
+		bynL[8] /= totalA;
+		bynL[9] /= totalA;
+		bynL[10] /= 2; //catchall for garbage data
+		//end L
+		
+		//begin R
+		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;
+			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;
+			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;
+			inputSampleR = floor(inputSampleR);
+		}
+		else
+		{
+			bynR[hotbinB] += 1;
+			inputSampleR = ceil(inputSampleR);
+		}
+		//assign the relevant one to the delay line
+		//and floor/ceil signal accordingly
+		
+		totalA = bynR[1] + bynR[2] + bynR[3] + bynR[4] + bynR[5] + bynR[6] + bynR[7] + bynR[8] + bynR[9];
+		totalA /= 1000;
+		if (totalA = 0) totalA = 1;
+		bynR[1] /= totalA;
+		bynR[2] /= totalA;
+		bynR[3] /= totalA;
+		bynR[4] /= totalA;
+		bynR[5] /= totalA;
+		bynR[6] /= totalA;
+		bynR[7] /= totalA;
+		bynR[8] /= totalA;
+		bynR[9] /= totalA;
+		bynR[10] /= 2; //catchall for garbage data
+		//end R
+		
+		inputSampleL /= 8388608.0;
+		inputSampleR /= 8388608.0;
+		
+		noiseShapingL += inputSampleL - drySampleL;
+		noiseShapingR += inputSampleR - drySampleR;
+		
+		*out1 = inputSampleL;
+		*out2 = inputSampleR;
+		
+		*in1++;
+		*in2++;
+		*out1++;
+		*out2++;
+    }
+}
\ No newline at end of file