diff options
author | Chris Johnson <jinx6568@sover.net> | 2019-01-14 04:29:08 -0500 |
---|---|---|
committer | Chris Johnson <jinx6568@sover.net> | 2019-01-14 04:29:08 -0500 |
commit | 0887543349dbbec0721a1fc8b1c7deba9afefa8b (patch) | |
tree | 40c650e606e37edb73221625b31d43bff4fb1013 /plugins/LinuxVST | |
parent | c27e202a05f5ff3bd1d0f060215e5024e5c7d099 (diff) | |
download | airwindows-lv2-port-0887543349dbbec0721a1fc8b1c7deba9afefa8b.tar.gz airwindows-lv2-port-0887543349dbbec0721a1fc8b1c7deba9afefa8b.tar.bz2 airwindows-lv2-port-0887543349dbbec0721a1fc8b1c7deba9afefa8b.zip |
StudioTan
Diffstat (limited to 'plugins/LinuxVST')
-rwxr-xr-x | plugins/LinuxVST/CMakeLists.txt | 1 | ||||
-rwxr-xr-x | plugins/LinuxVST/src/StudioTan/StudioTan.cpp | 156 | ||||
-rwxr-xr-x | plugins/LinuxVST/src/StudioTan/StudioTan.h | 69 | ||||
-rwxr-xr-x | plugins/LinuxVST/src/StudioTan/StudioTanProc.cpp | 564 |
4 files changed, 790 insertions, 0 deletions
diff --git a/plugins/LinuxVST/CMakeLists.txt b/plugins/LinuxVST/CMakeLists.txt index 90b707d..ea7f837 100755 --- a/plugins/LinuxVST/CMakeLists.txt +++ b/plugins/LinuxVST/CMakeLists.txt @@ -123,6 +123,7 @@ add_airwindows_plugin(Spiral) add_airwindows_plugin(Spiral2) add_airwindows_plugin(StarChild) add_airwindows_plugin(StereoFX) +add_airwindows_plugin(StudioTan) add_airwindows_plugin(SubsOnly) add_airwindows_plugin(Surge) add_airwindows_plugin(SurgeTide) diff --git a/plugins/LinuxVST/src/StudioTan/StudioTan.cpp b/plugins/LinuxVST/src/StudioTan/StudioTan.cpp new file mode 100755 index 0000000..31242f7 --- /dev/null +++ b/plugins/LinuxVST/src/StudioTan/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/LinuxVST/src/StudioTan/StudioTan.h b/plugins/LinuxVST/src/StudioTan/StudioTan.h new file mode 100755 index 0000000..66caf48 --- /dev/null +++ b/plugins/LinuxVST/src/StudioTan/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/LinuxVST/src/StudioTan/StudioTanProc.cpp b/plugins/LinuxVST/src/StudioTan/StudioTanProc.cpp new file mode 100755 index 0000000..3258f7b --- /dev/null +++ b/plugins/LinuxVST/src/StudioTan/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++; + } +} |