diff options
author | Chris Johnson <jinx6568@sover.net> | 2018-10-22 18:04:06 -0400 |
---|---|---|
committer | Chris Johnson <jinx6568@sover.net> | 2018-10-22 18:04:06 -0400 |
commit | 633be2e22c6648c901f08f3b4cd4e8e14ea86443 (patch) | |
tree | 1e272c3d2b5bd29636b9f9f521af62734e4df012 /plugins/MacVST/Melt/source | |
parent | 057757aa8eb0a463caf0cdfdb5894ac5f723ff3f (diff) | |
download | airwindows-lv2-port-633be2e22c6648c901f08f3b4cd4e8e14ea86443.tar.gz airwindows-lv2-port-633be2e22c6648c901f08f3b4cd4e8e14ea86443.tar.bz2 airwindows-lv2-port-633be2e22c6648c901f08f3b4cd4e8e14ea86443.zip |
Updates (in case my plane crashes)
Diffstat (limited to 'plugins/MacVST/Melt/source')
-rwxr-xr-x | plugins/MacVST/Melt/source/Melt.cpp | 156 | ||||
-rwxr-xr-x | plugins/MacVST/Melt/source/Melt.h | 87 | ||||
-rwxr-xr-x | plugins/MacVST/Melt/source/MeltProc.cpp | 596 |
3 files changed, 839 insertions, 0 deletions
diff --git a/plugins/MacVST/Melt/source/Melt.cpp b/plugins/MacVST/Melt/source/Melt.cpp new file mode 100755 index 0000000..9722e1f --- /dev/null +++ b/plugins/MacVST/Melt/source/Melt.cpp @@ -0,0 +1,156 @@ +/* ======================================== + * Melt - Melt.h + * Copyright (c) 2016 airwindows, All rights reserved + * ======================================== */ + +#ifndef __Melt_H +#include "Melt.h" +#endif + +AudioEffect* createEffectInstance(audioMasterCallback audioMaster) {return new Melt(audioMaster);} + +Melt::Melt(audioMasterCallback audioMaster) : + AudioEffectX(audioMaster, kNumPrograms, kNumParameters) +{ + A = 0.25; + B = 0.75; + C = 1.0; + D = 1.0; + + for(int count = 0; count < 32001; count++) {dL[count] = 0;dR[count] = 0;} + combineL = 0; combineR = 0; + scalefactorL = 0.999; scalefactorR = 0.999; + + for(int count = 0; count < 31; count++) {minTap[count] = 0; maxTap[count] = 0; position[count] = 1; stepTap[count] = 1;} + stepCount = 0; + slowCount = 0; + gcount = 0; + + fpNShapeLA = 0.0; + fpNShapeLB = 0.0; + fpNShapeRA = 0.0; + fpNShapeRB = 0.0; + fpFlip = true; + //this is reset: values being initialized only once. Startup values, whatever they are. + + _canDo.insert("plugAsChannelInsert"); // plug-in can be used as a channel insert effect. + _canDo.insert("plugAsSend"); // plug-in can be used as a send effect. + _canDo.insert("x2in2out"); + setNumInputs(kNumInputs); + setNumOutputs(kNumOutputs); + setUniqueID(kUniqueId); + canProcessReplacing(); // supports output replacing + canDoubleReplacing(); // supports double precision processing + programsAreChunks(true); + vst_strncpy (_programName, "Default", kVstMaxProgNameLen); // default program name +} + +Melt::~Melt() {} +VstInt32 Melt::getVendorVersion () {return 1000;} +void Melt::setProgramName(char *name) {vst_strncpy (_programName, name, kVstMaxProgNameLen);} +void Melt::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 Melt::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 Melt::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 Melt::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 Melt::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 Melt::getParameterName(VstInt32 index, char *text) { + switch (index) { + case kParamA: vst_strncpy (text, "Depth", kVstMaxParamStrLen); break; + case kParamB: vst_strncpy (text, "Range", 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 Melt::getParameterDisplay(VstInt32 index, char *text) { + switch (index) { + case kParamA: float2string (A, text, kVstMaxParamStrLen); break; + case kParamB: float2string (B, text, kVstMaxParamStrLen); break; //also display 0-1 as percent + 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 Melt::getParameterLabel(VstInt32 index, char *text) { + switch (index) { + case kParamA: vst_strncpy (text, " ", kVstMaxParamStrLen); break; + case kParamB: vst_strncpy (text, " ", kVstMaxParamStrLen); break; //the percent + case kParamC: vst_strncpy (text, " ", kVstMaxParamStrLen); break; + case kParamD: vst_strncpy (text, " ", kVstMaxParamStrLen); break; //the popup + default: break; // unknown parameter, shouldn't happen! + } +} + +VstInt32 Melt::canDo(char *text) +{ return (_canDo.find(text) == _canDo.end()) ? -1: 1; } // 1 = yes, -1 = no, 0 = don't know + +bool Melt::getEffectName(char* name) { + vst_strncpy(name, "Melt", kVstMaxProductStrLen); return true; +} + +VstPlugCategory Melt::getPlugCategory() {return kPlugCategEffect;} + +bool Melt::getProductString(char* text) { + vst_strncpy (text, "airwindows Melt", kVstMaxProductStrLen); return true; +} + +bool Melt::getVendorString(char* text) { + vst_strncpy (text, "airwindows", kVstMaxVendorStrLen); return true; +} diff --git a/plugins/MacVST/Melt/source/Melt.h b/plugins/MacVST/Melt/source/Melt.h new file mode 100755 index 0000000..a8e0f38 --- /dev/null +++ b/plugins/MacVST/Melt/source/Melt.h @@ -0,0 +1,87 @@ +/* ======================================== + * Melt - Melt.h + * Created 8/12/11 by SPIAdmin + * Copyright (c) 2011 __MyCompanyName__, All rights reserved + * ======================================== */ + +#ifndef __Melt_H +#define __Melt_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 = 'melt'; //Change this to what the AU identity is! + +class Melt : + public AudioEffectX +{ +public: + Melt(audioMasterCallback audioMaster); + ~Melt(); + 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; + + float dL[32002]; + float dR[32002]; + float combineL; + float combineR; + float scalefactorL; + float scalefactorR; + + int minTap[32]; + int maxTap[32]; + int position[32]; + int stepTap[32]; + int stepCount; + int slowCount; + int gcount; + + long double fpNShapeLA; + long double fpNShapeLB; + long double fpNShapeRA; + long double fpNShapeRB; + bool fpFlip; + //default stuff + + float A; + float B; + float C; + float D; //parameters. Always 0-1, and we scale/alter them elsewhere. + +}; + +#endif diff --git a/plugins/MacVST/Melt/source/MeltProc.cpp b/plugins/MacVST/Melt/source/MeltProc.cpp new file mode 100755 index 0000000..48a1c11 --- /dev/null +++ b/plugins/MacVST/Melt/source/MeltProc.cpp @@ -0,0 +1,596 @@ +/* ======================================== + * Melt - Melt.h + * Copyright (c) 2016 airwindows, All rights reserved + * ======================================== */ + +#ifndef __Melt_H +#include "Melt.h" +#endif + +void Melt::processReplacing(float **inputs, float **outputs, VstInt32 sampleFrames) +{ + float* in1 = inputs[0]; + float* in2 = inputs[1]; + float* out1 = outputs[0]; + float* out2 = outputs[1]; + + double rate = 1 / (pow(A,2) + 0.001); + double depthB = (B * 139.5)+2; + double depthA = depthB * (1.0 - A); + double output = C * 0.05; + double wet = D; + double dry = 1.0-wet; + + float fpTemp; + long double fpOld = 0.618033988749894848204586; //golden ratio! + long double fpNew = 1.0 - fpOld; + + long double inputSampleL; + long double inputSampleR; + long double drySampleL; + long double drySampleR; + + minTap[0] = floor(2 * depthA); maxTap[0] = floor(2 * depthB); + minTap[1] = floor(3 * depthA); maxTap[1] = floor(3 * depthB); + minTap[2] = floor(5 * depthA); maxTap[2] = floor(5 * depthB); + minTap[3] = floor(7 * depthA); maxTap[3] = floor(7 * depthB); + minTap[4] = floor(11 * depthA); maxTap[4] = floor(11 * depthB); + minTap[5] = floor(13 * depthA); maxTap[5] = floor(13 * depthB); + minTap[6] = floor(17 * depthA); maxTap[6] = floor(17 * depthB); + minTap[7] = floor(19 * depthA); maxTap[7] = floor(19 * depthB); + minTap[8] = floor(23 * depthA); maxTap[8] = floor(23 * depthB); + minTap[9] = floor(29 * depthA); maxTap[9] = floor(29 * depthB); + minTap[10] = floor(31 * depthA); maxTap[10] = floor(31 * depthB); + minTap[11] = floor(37 * depthA); maxTap[11] = floor(37 * depthB); + minTap[12] = floor(41 * depthA); maxTap[12] = floor(41 * depthB); + minTap[13] = floor(43 * depthA); maxTap[13] = floor(43 * depthB); + minTap[14] = floor(47 * depthA); maxTap[14] = floor(47 * depthB); + minTap[15] = floor(53 * depthA); maxTap[15] = floor(53 * depthB); + minTap[16] = floor(59 * depthA); maxTap[16] = floor(59 * depthB); + minTap[17] = floor(61 * depthA); maxTap[17] = floor(61 * depthB); + minTap[18] = floor(67 * depthA); maxTap[18] = floor(67 * depthB); + minTap[19] = floor(71 * depthA); maxTap[19] = floor(71 * depthB); + minTap[20] = floor(73 * depthA); maxTap[20] = floor(73 * depthB); + minTap[21] = floor(79 * depthA); maxTap[21] = floor(79 * depthB); + minTap[22] = floor(83 * depthA); maxTap[22] = floor(83 * depthB); + minTap[23] = floor(89 * depthA); maxTap[23] = floor(89 * depthB); + minTap[24] = floor(97 * depthA); maxTap[24] = floor(97 * depthB); + minTap[25] = floor(101 * depthA); maxTap[25] = floor(101 * depthB); + minTap[26] = floor(103 * depthA); maxTap[26] = floor(103 * depthB); + minTap[27] = floor(107 * depthA); maxTap[27] = floor(107 * depthB); + minTap[28] = floor(109 * depthA); maxTap[28] = floor(109 * depthB); + minTap[29] = floor(113 * depthA); maxTap[29] = floor(113 * depthB); + minTap[30] = floor(117 * depthA); maxTap[30] = floor(117 * depthB); + + 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; + + if (gcount < 0 || gcount > 16000) {gcount = 16000;} + dL[gcount+16000] = dL[gcount] = inputSampleL; + dR[gcount+16000] = dR[gcount] = inputSampleR; + + if (slowCount > rate || slowCount < 0) { + slowCount = 0; + stepCount++; + if (stepCount > 29 || stepCount < 0) {stepCount = 0;} + position[stepCount] += stepTap[stepCount]; + if (position[stepCount] < minTap[stepCount]) { + position[stepCount] = minTap[stepCount]; + stepTap[stepCount] = 1; + } + if (position[stepCount] > maxTap[stepCount]) { + position[stepCount] = maxTap[stepCount]; + stepTap[stepCount] = -1; + } + } + + //begin L + scalefactorL *= 0.9999; + scalefactorL += (100.0 - fabs(combineL)) * 0.000001; + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[29]]); + combineL += (dL[gcount+position[28]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[27]]); + combineL += (dL[gcount+position[26]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[25]]); + combineL += (dL[gcount+position[24]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[23]]); + combineL += (dL[gcount+position[22]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[21]]); + combineL += (dL[gcount+position[20]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[19]]); + combineL += (dL[gcount+position[18]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[17]]); + combineL += (dL[gcount+position[16]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[15]]); + combineL += (dL[gcount+position[14]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[13]]); + combineL += (dL[gcount+position[12]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[11]]); + combineL += (dL[gcount+position[10]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[9]]); + combineL += (dL[gcount+position[8]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[7]]); + combineL += (dL[gcount+position[6]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[5]]); + combineL += (dL[gcount+position[4]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[3]]); + combineL += (dL[gcount+position[2]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[1]]); + combineL += (dL[gcount+position[0]]); + + inputSampleL = combineL; + //done with L + + //begin R + scalefactorR *= 0.9999; + scalefactorR += (100.0 - fabs(combineR)) * 0.000001; + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[29]]); + combineR += (dR[gcount+position[28]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[27]]); + combineR += (dR[gcount+position[26]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[25]]); + combineR += (dR[gcount+position[24]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[23]]); + combineR += (dR[gcount+position[22]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[21]]); + combineR += (dR[gcount+position[20]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[19]]); + combineR += (dR[gcount+position[18]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[17]]); + combineR += (dR[gcount+position[16]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[15]]); + combineR += (dR[gcount+position[14]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[13]]); + combineR += (dR[gcount+position[12]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[11]]); + combineR += (dR[gcount+position[10]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[9]]); + combineR += (dR[gcount+position[8]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[7]]); + combineR += (dR[gcount+position[6]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[5]]); + combineR += (dR[gcount+position[4]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[3]]); + combineR += (dR[gcount+position[2]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[1]]); + combineR += (dR[gcount+position[0]]); + + inputSampleR = combineR; + //done with R + + gcount--; + slowCount++; + + if (output < 1.0) {inputSampleL *= output; inputSampleR *= output;} + if (wet < 1.0) { + inputSampleL = (drySampleL * dry)+(inputSampleL*wet); + inputSampleR = (drySampleR * dry)+(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 + if (fpFlip) { + fpTemp = inputSampleL; + fpNShapeLA = (fpNShapeLA*fpOld)+((inputSampleL-fpTemp)*fpNew); + inputSampleL += fpNShapeLA; + fpTemp = inputSampleR; + fpNShapeRA = (fpNShapeRA*fpOld)+((inputSampleR-fpTemp)*fpNew); + inputSampleR += fpNShapeRA; + } + else { + fpTemp = inputSampleL; + fpNShapeLB = (fpNShapeLB*fpOld)+((inputSampleL-fpTemp)*fpNew); + inputSampleL += fpNShapeLB; + fpTemp = inputSampleR; + fpNShapeRB = (fpNShapeRB*fpOld)+((inputSampleR-fpTemp)*fpNew); + inputSampleR += fpNShapeRB; + } + fpFlip = !fpFlip; + //end noise shaping on 32 bit output + + *out1 = inputSampleL; + *out2 = inputSampleR; + + *in1++; + *in2++; + *out1++; + *out2++; + } +} + +void Melt::processDoubleReplacing(double **inputs, double **outputs, VstInt32 sampleFrames) +{ + double* in1 = inputs[0]; + double* in2 = inputs[1]; + double* out1 = outputs[0]; + double* out2 = outputs[1]; + + double rate = 1 / (pow(A,2) + 0.001); + double depthB = (B * 139.5)+2; + double depthA = depthB * (1.0 - A); + double output = C * 0.05; + double wet = D; + double dry = 1.0-wet; + + double fpTemp; //this is different from singlereplacing + long double fpOld = 0.618033988749894848204586; //golden ratio! + long double fpNew = 1.0 - fpOld; + + long double inputSampleL; + long double inputSampleR; + long double drySampleL; + long double drySampleR; + + minTap[0] = floor(2 * depthA); maxTap[0] = floor(2 * depthB); + minTap[1] = floor(3 * depthA); maxTap[1] = floor(3 * depthB); + minTap[2] = floor(5 * depthA); maxTap[2] = floor(5 * depthB); + minTap[3] = floor(7 * depthA); maxTap[3] = floor(7 * depthB); + minTap[4] = floor(11 * depthA); maxTap[4] = floor(11 * depthB); + minTap[5] = floor(13 * depthA); maxTap[5] = floor(13 * depthB); + minTap[6] = floor(17 * depthA); maxTap[6] = floor(17 * depthB); + minTap[7] = floor(19 * depthA); maxTap[7] = floor(19 * depthB); + minTap[8] = floor(23 * depthA); maxTap[8] = floor(23 * depthB); + minTap[9] = floor(29 * depthA); maxTap[9] = floor(29 * depthB); + minTap[10] = floor(31 * depthA); maxTap[10] = floor(31 * depthB); + minTap[11] = floor(37 * depthA); maxTap[11] = floor(37 * depthB); + minTap[12] = floor(41 * depthA); maxTap[12] = floor(41 * depthB); + minTap[13] = floor(43 * depthA); maxTap[13] = floor(43 * depthB); + minTap[14] = floor(47 * depthA); maxTap[14] = floor(47 * depthB); + minTap[15] = floor(53 * depthA); maxTap[15] = floor(53 * depthB); + minTap[16] = floor(59 * depthA); maxTap[16] = floor(59 * depthB); + minTap[17] = floor(61 * depthA); maxTap[17] = floor(61 * depthB); + minTap[18] = floor(67 * depthA); maxTap[18] = floor(67 * depthB); + minTap[19] = floor(71 * depthA); maxTap[19] = floor(71 * depthB); + minTap[20] = floor(73 * depthA); maxTap[20] = floor(73 * depthB); + minTap[21] = floor(79 * depthA); maxTap[21] = floor(79 * depthB); + minTap[22] = floor(83 * depthA); maxTap[22] = floor(83 * depthB); + minTap[23] = floor(89 * depthA); maxTap[23] = floor(89 * depthB); + minTap[24] = floor(97 * depthA); maxTap[24] = floor(97 * depthB); + minTap[25] = floor(101 * depthA); maxTap[25] = floor(101 * depthB); + minTap[26] = floor(103 * depthA); maxTap[26] = floor(103 * depthB); + minTap[27] = floor(107 * depthA); maxTap[27] = floor(107 * depthB); + minTap[28] = floor(109 * depthA); maxTap[28] = floor(109 * depthB); + minTap[29] = floor(113 * depthA); maxTap[29] = floor(113 * depthB); + minTap[30] = floor(117 * depthA); maxTap[30] = floor(117 * depthB); + + 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; + + if (gcount < 0 || gcount > 16000) {gcount = 16000;} + dL[gcount+16000] = dL[gcount] = inputSampleL; + dR[gcount+16000] = dR[gcount] = inputSampleR; + + if (slowCount > rate || slowCount < 0) { + slowCount = 0; + stepCount++; + if (stepCount > 29 || stepCount < 0) {stepCount = 0;} + position[stepCount] += stepTap[stepCount]; + if (position[stepCount] < minTap[stepCount]) { + position[stepCount] = minTap[stepCount]; + stepTap[stepCount] = 1; + } + if (position[stepCount] > maxTap[stepCount]) { + position[stepCount] = maxTap[stepCount]; + stepTap[stepCount] = -1; + } + } + + //begin L + scalefactorL *= 0.9999; + scalefactorL += (100.0 - fabs(combineL)) * 0.000001; + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[29]]); + combineL += (dL[gcount+position[28]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[27]]); + combineL += (dL[gcount+position[26]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[25]]); + combineL += (dL[gcount+position[24]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[23]]); + combineL += (dL[gcount+position[22]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[21]]); + combineL += (dL[gcount+position[20]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[19]]); + combineL += (dL[gcount+position[18]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[17]]); + combineL += (dL[gcount+position[16]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[15]]); + combineL += (dL[gcount+position[14]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[13]]); + combineL += (dL[gcount+position[12]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[11]]); + combineL += (dL[gcount+position[10]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[9]]); + combineL += (dL[gcount+position[8]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[7]]); + combineL += (dL[gcount+position[6]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[5]]); + combineL += (dL[gcount+position[4]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[3]]); + combineL += (dL[gcount+position[2]]); + + combineL *= scalefactorL; + combineL -= (dL[gcount+position[1]]); + combineL += (dL[gcount+position[0]]); + + inputSampleL = combineL; + //done with L + + //begin R + scalefactorR *= 0.9999; + scalefactorR += (100.0 - fabs(combineR)) * 0.000001; + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[29]]); + combineR += (dR[gcount+position[28]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[27]]); + combineR += (dR[gcount+position[26]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[25]]); + combineR += (dR[gcount+position[24]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[23]]); + combineR += (dR[gcount+position[22]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[21]]); + combineR += (dR[gcount+position[20]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[19]]); + combineR += (dR[gcount+position[18]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[17]]); + combineR += (dR[gcount+position[16]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[15]]); + combineR += (dR[gcount+position[14]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[13]]); + combineR += (dR[gcount+position[12]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[11]]); + combineR += (dR[gcount+position[10]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[9]]); + combineR += (dR[gcount+position[8]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[7]]); + combineR += (dR[gcount+position[6]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[5]]); + combineR += (dR[gcount+position[4]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[3]]); + combineR += (dR[gcount+position[2]]); + + combineR *= scalefactorR; + combineR -= (dR[gcount+position[1]]); + combineR += (dR[gcount+position[0]]); + + inputSampleR = combineR; + //done with R + + gcount--; + slowCount++; + + if (output < 1.0) {inputSampleL *= output; inputSampleR *= output;} + if (wet < 1.0) { + inputSampleL = (drySampleL * dry)+(inputSampleL*wet); + inputSampleR = (drySampleR * dry)+(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 + if (fpFlip) { + fpTemp = inputSampleL; + fpNShapeLA = (fpNShapeLA*fpOld)+((inputSampleL-fpTemp)*fpNew); + inputSampleL += fpNShapeLA; + fpTemp = inputSampleR; + fpNShapeRA = (fpNShapeRA*fpOld)+((inputSampleR-fpTemp)*fpNew); + inputSampleR += fpNShapeRA; + } + else { + fpTemp = inputSampleL; + fpNShapeLB = (fpNShapeLB*fpOld)+((inputSampleL-fpTemp)*fpNew); + inputSampleL += fpNShapeLB; + fpTemp = inputSampleR; + fpNShapeRB = (fpNShapeRB*fpOld)+((inputSampleR-fpTemp)*fpNew); + inputSampleR += fpNShapeRB; + } + fpFlip = !fpFlip; + //end noise shaping on 64 bit output + + *out1 = inputSampleL; + *out2 = inputSampleR; + + *in1++; + *in2++; + *out1++; + *out2++; + } +}
\ No newline at end of file |