From 04af2f23dd59cfb33c1cc3d8760ecb8b92bddce3 Mon Sep 17 00:00:00 2001 From: Chris Johnson Date: Sun, 3 Mar 2019 21:38:54 -0500 Subject: PurestSquish --- plugins/LinuxVST/CMakeLists.txt | 1 + plugins/LinuxVST/src/PurestSquish/PurestSquish.cpp | 190 ++++ plugins/LinuxVST/src/PurestSquish/PurestSquish.h | 118 +++ .../LinuxVST/src/PurestSquish/PurestSquishProc.cpp | 1016 ++++++++++++++++++++ 4 files changed, 1325 insertions(+) create mode 100755 plugins/LinuxVST/src/PurestSquish/PurestSquish.cpp create mode 100755 plugins/LinuxVST/src/PurestSquish/PurestSquish.h create mode 100755 plugins/LinuxVST/src/PurestSquish/PurestSquishProc.cpp (limited to 'plugins/LinuxVST') diff --git a/plugins/LinuxVST/CMakeLists.txt b/plugins/LinuxVST/CMakeLists.txt index d8764e4..02cf1f0 100755 --- a/plugins/LinuxVST/CMakeLists.txt +++ b/plugins/LinuxVST/CMakeLists.txt @@ -112,6 +112,7 @@ add_airwindows_plugin(PurestConsoleChannel) add_airwindows_plugin(PurestDrive) add_airwindows_plugin(PurestEcho) add_airwindows_plugin(PurestGain) +add_airwindows_plugin(PurestSquish) add_airwindows_plugin(PurestWarm) add_airwindows_plugin(Pyewacket) add_airwindows_plugin(RawGlitters) diff --git a/plugins/LinuxVST/src/PurestSquish/PurestSquish.cpp b/plugins/LinuxVST/src/PurestSquish/PurestSquish.cpp new file mode 100755 index 0000000..16543c4 --- /dev/null +++ b/plugins/LinuxVST/src/PurestSquish/PurestSquish.cpp @@ -0,0 +1,190 @@ +/* ======================================== + * PurestSquish - PurestSquish.h + * Copyright (c) 2016 airwindows, All rights reserved + * ======================================== */ + +#ifndef __PurestSquish_H +#include "PurestSquish.h" +#endif + +AudioEffect* createEffectInstance(audioMasterCallback audioMaster) {return new PurestSquish(audioMaster);} + +PurestSquish::PurestSquish(audioMasterCallback audioMaster) : + AudioEffectX(audioMaster, kNumPrograms, kNumParameters) +{ + A = 0.0; + B = 0.0; + C = 1.0; + D = 1.0; + + muSpeedAL = 10000; + muSpeedBL = 10000; + muSpeedCL = 10000; + muSpeedDL = 10000; + muSpeedEL = 10000; + muCoefficientAL = 1; + muCoefficientBL = 1; + muCoefficientCL = 1; + muCoefficientDL = 1; + muCoefficientEL = 1; + iirSampleAL = 0.0; + iirSampleBL = 0.0; + iirSampleCL = 0.0; + iirSampleDL = 0.0; + iirSampleEL = 0.0; + lastCoefficientAL = 1; + lastCoefficientBL = 1; + lastCoefficientCL = 1; + lastCoefficientDL = 1; + mergedCoefficientsL = 1; + muVaryL = 1; + + muSpeedAR = 10000; + muSpeedBR = 10000; + muSpeedCR = 10000; + muSpeedDR = 10000; + muSpeedER = 10000; + muCoefficientAR = 1; + muCoefficientBR = 1; + muCoefficientCR = 1; + muCoefficientDR = 1; + muCoefficientER = 1; + iirSampleAR = 0.0; + iirSampleBR = 0.0; + iirSampleCR = 0.0; + iirSampleDR = 0.0; + iirSampleER = 0.0; + lastCoefficientAR = 1; + lastCoefficientBR = 1; + lastCoefficientCR = 1; + lastCoefficientDR = 1; + mergedCoefficientsR = 1; + muVaryR = 1; + + count = 1; + fpFlip = true; + + fpd = 17; + //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 +} + +PurestSquish::~PurestSquish() {} +VstInt32 PurestSquish::getVendorVersion () {return 1000;} +void PurestSquish::setProgramName(char *name) {vst_strncpy (_programName, name, kVstMaxProgNameLen);} +void PurestSquish::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 PurestSquish::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 PurestSquish::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 PurestSquish::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 PurestSquish::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 PurestSquish::getParameterName(VstInt32 index, char *text) { + switch (index) { + case kParamA: vst_strncpy (text, "Squish", kVstMaxParamStrLen); break; + case kParamB: vst_strncpy (text, "BassBlm", 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 PurestSquish::getParameterDisplay(VstInt32 index, char *text) { + switch (index) { + case kParamA: float2string (A, text, kVstMaxParamStrLen); break; + case kParamB: float2string (B, text, kVstMaxParamStrLen); break; + case kParamC: float2string (C, text, kVstMaxParamStrLen); break; + case kParamD: float2string (D, text, kVstMaxParamStrLen); break; + default: break; // unknown parameter, shouldn't happen! + } //this displays the values and handles 'popups' where it's discrete choices +} + +void PurestSquish::getParameterLabel(VstInt32 index, char *text) { + switch (index) { + case kParamA: vst_strncpy (text, "", kVstMaxParamStrLen); break; + case kParamB: vst_strncpy (text, "", kVstMaxParamStrLen); break; + case kParamC: vst_strncpy (text, "", kVstMaxParamStrLen); break; + case kParamD: vst_strncpy (text, "", kVstMaxParamStrLen); break; + default: break; // unknown parameter, shouldn't happen! + } +} + +VstInt32 PurestSquish::canDo(char *text) +{ return (_canDo.find(text) == _canDo.end()) ? -1: 1; } // 1 = yes, -1 = no, 0 = don't know + +bool PurestSquish::getEffectName(char* name) { + vst_strncpy(name, "PurestSquish", kVstMaxProductStrLen); return true; +} + +VstPlugCategory PurestSquish::getPlugCategory() {return kPlugCategEffect;} + +bool PurestSquish::getProductString(char* text) { + vst_strncpy (text, "airwindows PurestSquish", kVstMaxProductStrLen); return true; +} + +bool PurestSquish::getVendorString(char* text) { + vst_strncpy (text, "airwindows", kVstMaxVendorStrLen); return true; +} diff --git a/plugins/LinuxVST/src/PurestSquish/PurestSquish.h b/plugins/LinuxVST/src/PurestSquish/PurestSquish.h new file mode 100755 index 0000000..d5b93a9 --- /dev/null +++ b/plugins/LinuxVST/src/PurestSquish/PurestSquish.h @@ -0,0 +1,118 @@ +/* ======================================== + * PurestSquish - PurestSquish.h + * Created 8/12/11 by SPIAdmin + * Copyright (c) 2011 __MyCompanyName__, All rights reserved + * ======================================== */ + +#ifndef __PurestSquish_H +#define __PurestSquish_H + +#ifndef __audioeffect__ +#include "audioeffectx.h" +#endif + +#include +#include +#include + +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 = 'pusq'; //Change this to what the AU identity is! + +class PurestSquish : + public AudioEffectX +{ +public: + PurestSquish(audioMasterCallback audioMaster); + ~PurestSquish(); + 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; + + uint32_t fpd; + //default stuff + + float A; + float B; + float C; + float D; + + double muVaryL; + double muAttackL; + double muNewSpeedL; + double muSpeedAL; + double muSpeedBL; + double muSpeedCL; + double muSpeedDL; + double muSpeedEL; + double muCoefficientAL; + double muCoefficientBL; + double muCoefficientCL; + double muCoefficientDL; + double muCoefficientEL; + double iirSampleAL; + double iirSampleBL; + double iirSampleCL; + double iirSampleDL; + double iirSampleEL; + double lastCoefficientAL; + double lastCoefficientBL; + double lastCoefficientCL; + double lastCoefficientDL; + double mergedCoefficientsL; + + double muVaryR; + double muAttackR; + double muNewSpeedR; + double muSpeedAR; + double muSpeedBR; + double muSpeedCR; + double muSpeedDR; + double muSpeedER; + double muCoefficientAR; + double muCoefficientBR; + double muCoefficientCR; + double muCoefficientDR; + double muCoefficientER; + double iirSampleAR; + double iirSampleBR; + double iirSampleCR; + double iirSampleDR; + double iirSampleER; + double lastCoefficientAR; + double lastCoefficientBR; + double lastCoefficientCR; + double lastCoefficientDR; + double mergedCoefficientsR; + + int count; + bool fpFlip; +}; + +#endif diff --git a/plugins/LinuxVST/src/PurestSquish/PurestSquishProc.cpp b/plugins/LinuxVST/src/PurestSquish/PurestSquishProc.cpp new file mode 100755 index 0000000..1be8537 --- /dev/null +++ b/plugins/LinuxVST/src/PurestSquish/PurestSquishProc.cpp @@ -0,0 +1,1016 @@ +/* ======================================== + * PurestSquish - PurestSquish.h + * Copyright (c) 2016 airwindows, All rights reserved + * ======================================== */ + +#ifndef __PurestSquish_H +#include "PurestSquish.h" +#endif + +void PurestSquish::processReplacing(float **inputs, float **outputs, VstInt32 sampleFrames) +{ + float* in1 = inputs[0]; + float* in2 = inputs[1]; + float* out1 = outputs[0]; + float* out2 = outputs[1]; + + double overallscale = 1.0; + overallscale /= 44100.0; + overallscale *= getSampleRate(); + + long double highpassSampleL; + long double highpassSampleBL; + long double highpassSampleR; + long double highpassSampleBR; + double threshold = 1.01 - (1.0-pow(1.0-(A*0.5),4)); + double iirAmount = pow(B,4)/overallscale; + double output = C; + double wet = D; + + while (--sampleFrames >= 0) + { + long double inputSampleL = *in1; + long double inputSampleR = *in2; + + static int noisesourceL = 0; + static int noisesourceR = 850010; + int residue; + double applyresidue; + + noisesourceL = noisesourceL % 1700021; noisesourceL++; + residue = noisesourceL * noisesourceL; + residue = residue % 170003; residue *= residue; + residue = residue % 17011; residue *= residue; + residue = residue % 1709; residue *= residue; + residue = residue % 173; residue *= residue; + residue = residue % 17; + applyresidue = residue; + applyresidue *= 0.00000001; + applyresidue *= 0.00000001; + inputSampleL += applyresidue; + if (inputSampleL<1.2e-38 && -inputSampleL<1.2e-38) { + inputSampleL -= applyresidue; + } + + noisesourceR = noisesourceR % 1700021; noisesourceR++; + residue = noisesourceR * noisesourceR; + residue = residue % 170003; residue *= residue; + residue = residue % 17011; residue *= residue; + residue = residue % 1709; residue *= residue; + residue = residue % 173; residue *= residue; + residue = residue % 17; + applyresidue = residue; + applyresidue *= 0.00000001; + applyresidue *= 0.00000001; + inputSampleR += applyresidue; + if (inputSampleR<1.2e-38 && -inputSampleR<1.2e-38) { + inputSampleR -= applyresidue; + } + //for live air, we always apply the dither noise. Then, if our result is + //effectively digital black, we'll subtract it again. We want a 'air' hiss + long double drySampleL = inputSampleL; + long double drySampleR = inputSampleR; + + double muMakeupGainL = 1.0 / threshold; + double outMakeupGainL = sqrt(muMakeupGainL); + muMakeupGainL += outMakeupGainL; + muMakeupGainL *= 0.5; + outMakeupGainL *= 0.5; + //gain settings around threshold + double releaseL = mergedCoefficientsL * 32768.0; + releaseL /= overallscale; + double fastestL = sqrt(releaseL); + //speed settings around release + double lastCorrectionL = mergedCoefficientsL; + // µ µ µ µ µ µ µ µ µ µ µ µ is the kitten song o/~ + + double muMakeupGainR = 1.0 / threshold; + double outMakeupGainR = sqrt(muMakeupGainR); + muMakeupGainR += outMakeupGainR; + muMakeupGainR *= 0.5; + outMakeupGainR *= 0.5; + //gain settings around threshold + double releaseR = mergedCoefficientsR * 32768.0; + releaseR /= overallscale; + double fastestR = sqrt(releaseR); + //speed settings around release + double lastCorrectionR = mergedCoefficientsR; + // µ µ µ µ µ µ µ µ µ µ µ µ is the kitten song o/~ + + if (muMakeupGainL != 1.0) inputSampleL = inputSampleL * muMakeupGainL; + if (muMakeupGainR != 1.0) inputSampleR = inputSampleR * muMakeupGainR; + + highpassSampleBL = highpassSampleL = inputSampleL; + highpassSampleBR = highpassSampleR = inputSampleR; + + if (count < 1 || count > 3) count = 1; + switch (count) + { + case 1: + //begin L + iirSampleAL = (iirSampleAL * (1 - iirAmount)) + (highpassSampleL * iirAmount); + highpassSampleL -= iirSampleAL; + if (fabs(highpassSampleL) > threshold) + { + muVaryL = threshold / fabs(highpassSampleL); + muAttackL = sqrt(fabs(muSpeedAL)); + muCoefficientAL = muCoefficientAL * (muAttackL-1.0); + if (muVaryL < threshold) + { + muCoefficientAL = muCoefficientAL + threshold; + } + else + { + muCoefficientAL = muCoefficientAL + muVaryL; + } + muCoefficientAL = muCoefficientAL / muAttackL; + } + else + { + muCoefficientAL = muCoefficientAL * ((muSpeedAL * muSpeedAL)-1.0); + muCoefficientAL = muCoefficientAL + 1.0; + muCoefficientAL = muCoefficientAL / (muSpeedAL * muSpeedAL); + } + muNewSpeedL = muSpeedAL * (muSpeedAL-1); + muNewSpeedL = muNewSpeedL + fabs(highpassSampleL*releaseL)+fastestL; + muSpeedAL = muNewSpeedL / muSpeedAL; + lastCoefficientAL = pow(muCoefficientAL,2); + mergedCoefficientsL = lastCoefficientBL; + mergedCoefficientsL += lastCoefficientAL; + lastCoefficientAL *= (1.0-lastCorrectionL); + lastCoefficientAL += (muCoefficientAL * lastCorrectionL); + lastCoefficientBL = lastCoefficientAL; + //end L + + //begin R + iirSampleAR = (iirSampleAR * (1 - iirAmount)) + (highpassSampleR * iirAmount); + highpassSampleR -= iirSampleAR; + if (fabs(highpassSampleR) > threshold) + { + muVaryR = threshold / fabs(highpassSampleR); + muAttackR = sqrt(fabs(muSpeedAR)); + muCoefficientAR = muCoefficientAR * (muAttackR-1.0); + if (muVaryR < threshold) + { + muCoefficientAR = muCoefficientAR + threshold; + } + else + { + muCoefficientAR = muCoefficientAR + muVaryR; + } + muCoefficientAR = muCoefficientAR / muAttackR; + } + else + { + muCoefficientAR = muCoefficientAR * ((muSpeedAR * muSpeedAR)-1.0); + muCoefficientAR = muCoefficientAR + 1.0; + muCoefficientAR = muCoefficientAR / (muSpeedAR * muSpeedAR); + } + muNewSpeedR = muSpeedAR * (muSpeedAR-1); + muNewSpeedR = muNewSpeedR + fabs(highpassSampleR*releaseR)+fastestR; + muSpeedAR = muNewSpeedR / muSpeedAR; + lastCoefficientAR = pow(muCoefficientAR,2); + mergedCoefficientsR = lastCoefficientBR; + mergedCoefficientsR += lastCoefficientAR; + lastCoefficientAR *= (1.0-lastCorrectionR); + lastCoefficientAR += (muCoefficientAR * lastCorrectionR); + lastCoefficientBR = lastCoefficientAR; + //end R + + break; + case 2: + //begin L + iirSampleBL = (iirSampleBL * (1 - iirAmount)) + (highpassSampleL * iirAmount); + highpassSampleL -= iirSampleBL; + if (fabs(highpassSampleL) > threshold) + { + muVaryL = threshold / fabs(highpassSampleL); + muAttackL = sqrt(fabs(muSpeedBL)); + muCoefficientBL = muCoefficientBL * (muAttackL-1); + if (muVaryL < threshold) + { + muCoefficientBL = muCoefficientBL + threshold; + } + else + { + muCoefficientBL = muCoefficientBL + muVaryL; + } + muCoefficientBL = muCoefficientBL / muAttackL; + } + else + { + muCoefficientBL = muCoefficientBL * ((muSpeedBL * muSpeedBL)-1.0); + muCoefficientBL = muCoefficientBL + 1.0; + muCoefficientBL = muCoefficientBL / (muSpeedBL * muSpeedBL); + } + muNewSpeedL = muSpeedBL * (muSpeedBL-1); + muNewSpeedL = muNewSpeedL + fabs(highpassSampleL*releaseL)+fastestL; + muSpeedBL = muNewSpeedL / muSpeedBL; + lastCoefficientAL = pow(muCoefficientBL,2); + mergedCoefficientsL = lastCoefficientBL; + mergedCoefficientsL += lastCoefficientAL; + lastCoefficientAL *= (1.0-lastCorrectionL); + lastCoefficientAL += (muCoefficientBL * lastCorrectionL); + lastCoefficientBL = lastCoefficientAL; + //end L + + //begin R + iirSampleBR = (iirSampleBR * (1 - iirAmount)) + (highpassSampleR * iirAmount); + highpassSampleR -= iirSampleBR; + if (fabs(highpassSampleR) > threshold) + { + muVaryR = threshold / fabs(highpassSampleR); + muAttackR = sqrt(fabs(muSpeedBR)); + muCoefficientBR = muCoefficientBR * (muAttackR-1); + if (muVaryR < threshold) + { + muCoefficientBR = muCoefficientBR + threshold; + } + else + { + muCoefficientBR = muCoefficientBR + muVaryR; + } + muCoefficientBR = muCoefficientBR / muAttackR; + } + else + { + muCoefficientBR = muCoefficientBR * ((muSpeedBR * muSpeedBR)-1.0); + muCoefficientBR = muCoefficientBR + 1.0; + muCoefficientBR = muCoefficientBR / (muSpeedBR * muSpeedBR); + } + muNewSpeedR = muSpeedBR * (muSpeedBR-1); + muNewSpeedR = muNewSpeedR + fabs(highpassSampleR*releaseR)+fastestR; + muSpeedBR = muNewSpeedR / muSpeedBR; + lastCoefficientAR = pow(muCoefficientBR,2); + mergedCoefficientsR = lastCoefficientBR; + mergedCoefficientsR += lastCoefficientAR; + lastCoefficientAR *= (1.0-lastCorrectionR); + lastCoefficientAR += (muCoefficientBR * lastCorrectionR); + lastCoefficientBR = lastCoefficientAR; + //end R + + break; + case 3: + //begin L + iirSampleCL = (iirSampleCL * (1 - iirAmount)) + (highpassSampleL * iirAmount); + highpassSampleL -= iirSampleCL; + if (fabs(highpassSampleL) > threshold) + { + muVaryL = threshold / fabs(highpassSampleL); + muAttackL = sqrt(fabs(muSpeedCL)); + muCoefficientCL = muCoefficientCL * (muAttackL-1); + if (muVaryL < threshold) + { + muCoefficientCL = muCoefficientCL + threshold; + } + else + { + muCoefficientCL = muCoefficientCL + muVaryL; + } + muCoefficientCL = muCoefficientCL / muAttackL; + } + else + { + muCoefficientCL = muCoefficientCL * ((muSpeedCL * muSpeedCL)-1.0); + muCoefficientCL = muCoefficientCL + 1.0; + muCoefficientCL = muCoefficientCL / (muSpeedCL * muSpeedCL); + } + muNewSpeedL = muSpeedCL * (muSpeedCL-1); + muNewSpeedL = muNewSpeedL + fabs(highpassSampleL*releaseL)+fastestL; + muSpeedCL = muNewSpeedL / muSpeedCL; + lastCoefficientAL = pow(muCoefficientCL,2); + mergedCoefficientsL = lastCoefficientBL; + mergedCoefficientsL += lastCoefficientAL; + lastCoefficientAL *= (1.0-lastCorrectionL); + lastCoefficientAL += (muCoefficientCL * lastCorrectionL); + lastCoefficientBL = lastCoefficientAL; + //end L + + //begin R + iirSampleCR = (iirSampleCR * (1 - iirAmount)) + (highpassSampleR * iirAmount); + highpassSampleR -= iirSampleCR; + if (fabs(highpassSampleR) > threshold) + { + muVaryR = threshold / fabs(highpassSampleR); + muAttackR = sqrt(fabs(muSpeedCR)); + muCoefficientCR = muCoefficientCR * (muAttackR-1); + if (muVaryR < threshold) + { + muCoefficientCR = muCoefficientCR + threshold; + } + else + { + muCoefficientCR = muCoefficientCR + muVaryR; + } + muCoefficientCR = muCoefficientCR / muAttackR; + } + else + { + muCoefficientCR = muCoefficientCR * ((muSpeedCR * muSpeedCR)-1.0); + muCoefficientCR = muCoefficientCR + 1.0; + muCoefficientCR = muCoefficientCR / (muSpeedCR * muSpeedCR); + } + muNewSpeedR = muSpeedCR * (muSpeedCR-1); + muNewSpeedR = muNewSpeedR + fabs(highpassSampleR*releaseR)+fastestR; + muSpeedCR = muNewSpeedR / muSpeedCR; + lastCoefficientAR = pow(muCoefficientCR,2); + mergedCoefficientsR = lastCoefficientBR; + mergedCoefficientsR += lastCoefficientAR; + lastCoefficientAR *= (1.0-lastCorrectionR); + lastCoefficientAR += (muCoefficientCR * lastCorrectionR); + lastCoefficientBR = lastCoefficientAR; + //end R + + break; + } + count++; + + //applied compression with vari-vari-µ-µ-µ-µ-µ-µ-is-the-kitten-song o/~ + //applied gain correction to control output level- tends to constrain sound rather than inflate it + + if (fpFlip) { + //begin L + iirSampleDL = (iirSampleDL * (1 - iirAmount)) + (highpassSampleBL * iirAmount); + highpassSampleBL -= iirSampleDL; + if (fabs(highpassSampleBL) > threshold) + { + muVaryL = threshold / fabs(highpassSampleBL); + muAttackL = sqrt(fabs(muSpeedDL)); + muCoefficientDL = muCoefficientDL * (muAttackL-1.0); + if (muVaryL < threshold) + { + muCoefficientDL = muCoefficientDL + threshold; + } + else + { + muCoefficientDL = muCoefficientDL + muVaryL; + } + muCoefficientDL = muCoefficientDL / muAttackL; + } + else + { + muCoefficientDL = muCoefficientDL * ((muSpeedDL * muSpeedDL)-1.0); + muCoefficientDL = muCoefficientDL + 1.0; + muCoefficientDL = muCoefficientDL / (muSpeedDL * muSpeedDL); + } + muNewSpeedL = muSpeedDL * (muSpeedDL-1); + muNewSpeedL = muNewSpeedL + fabs(highpassSampleBL*releaseL)+fastestL; + muSpeedDL = muNewSpeedL / muSpeedDL; + lastCoefficientCL = pow(muCoefficientEL,2); + mergedCoefficientsL += lastCoefficientDL; + mergedCoefficientsL += lastCoefficientCL; + lastCoefficientCL *= (1.0-lastCorrectionL); + lastCoefficientCL += (muCoefficientDL * lastCorrectionL); + lastCoefficientDL = lastCoefficientCL; + //end L + + //begin R + iirSampleDR = (iirSampleDR * (1 - iirAmount)) + (highpassSampleBR * iirAmount); + highpassSampleBR -= iirSampleDR; + if (fabs(highpassSampleBR) > threshold) + { + muVaryR = threshold / fabs(highpassSampleBR); + muAttackR = sqrt(fabs(muSpeedDR)); + muCoefficientDR = muCoefficientDR * (muAttackR-1.0); + if (muVaryR < threshold) + { + muCoefficientDR = muCoefficientDR + threshold; + } + else + { + muCoefficientDR = muCoefficientDR + muVaryR; + } + muCoefficientDR = muCoefficientDR / muAttackR; + } + else + { + muCoefficientDR = muCoefficientDR * ((muSpeedDR * muSpeedDR)-1.0); + muCoefficientDR = muCoefficientDR + 1.0; + muCoefficientDR = muCoefficientDR / (muSpeedDR * muSpeedDR); + } + muNewSpeedR = muSpeedDR * (muSpeedDR-1); + muNewSpeedR = muNewSpeedR + fabs(highpassSampleBR*releaseR)+fastestR; + muSpeedDR = muNewSpeedR / muSpeedDR; + lastCoefficientCR = pow(muCoefficientER,2); + mergedCoefficientsR += lastCoefficientDR; + mergedCoefficientsR += lastCoefficientCR; + lastCoefficientCR *= (1.0-lastCorrectionR); + lastCoefficientCR += (muCoefficientDR * lastCorrectionR); + lastCoefficientDR = lastCoefficientCR; + //end R + + } else { + //begin L + iirSampleEL = (iirSampleEL * (1 - iirAmount)) + (highpassSampleBL * iirAmount); + highpassSampleBL -= iirSampleEL; + if (fabs(highpassSampleBL) > threshold) + { + muVaryL = threshold / fabs(highpassSampleBL); + muAttackL = sqrt(fabs(muSpeedEL)); + muCoefficientEL = muCoefficientEL * (muAttackL-1.0); + if (muVaryL < threshold) + { + muCoefficientEL = muCoefficientEL + threshold; + } + else + { + muCoefficientEL = muCoefficientEL + muVaryL; + } + muCoefficientEL = muCoefficientEL / muAttackL; + } + else + { + muCoefficientEL = muCoefficientEL * ((muSpeedEL * muSpeedEL)-1.0); + muCoefficientEL = muCoefficientEL + 1.0; + muCoefficientEL = muCoefficientEL / (muSpeedEL * muSpeedEL); + } + muNewSpeedL = muSpeedEL * (muSpeedEL-1); + muNewSpeedL = muNewSpeedL + fabs(highpassSampleBL*releaseL)+fastestL; + muSpeedEL = muNewSpeedL / muSpeedEL; + lastCoefficientCL = pow(muCoefficientEL,2); + mergedCoefficientsL += lastCoefficientDL; + mergedCoefficientsL += lastCoefficientCL; + lastCoefficientCL *= (1.0-lastCorrectionL); + lastCoefficientCL += (muCoefficientEL * lastCorrectionL); + lastCoefficientDL = lastCoefficientCL; + //end L + + //begin R + iirSampleER = (iirSampleER * (1 - iirAmount)) + (highpassSampleBR * iirAmount); + highpassSampleBR -= iirSampleER; + if (fabs(highpassSampleBR) > threshold) + { + muVaryR = threshold / fabs(highpassSampleBR); + muAttackR = sqrt(fabs(muSpeedER)); + muCoefficientER = muCoefficientER * (muAttackR-1.0); + if (muVaryR < threshold) + { + muCoefficientER = muCoefficientER + threshold; + } + else + { + muCoefficientER = muCoefficientER + muVaryR; + } + muCoefficientER = muCoefficientER / muAttackR; + } + else + { + muCoefficientER = muCoefficientER * ((muSpeedER * muSpeedER)-1.0); + muCoefficientER = muCoefficientER + 1.0; + muCoefficientER = muCoefficientER / (muSpeedER * muSpeedER); + } + muNewSpeedR = muSpeedER * (muSpeedER-1); + muNewSpeedR = muNewSpeedR + fabs(highpassSampleBR*releaseR)+fastestR; + muSpeedER = muNewSpeedR / muSpeedER; + lastCoefficientCR = pow(muCoefficientER,2); + mergedCoefficientsR += lastCoefficientDR; + mergedCoefficientsR += lastCoefficientCR; + lastCoefficientCR *= (1.0-lastCorrectionR); + lastCoefficientCR += (muCoefficientER * lastCorrectionR); + lastCoefficientDR = lastCoefficientCR; + //end R + + } + mergedCoefficientsL *= 0.25; + inputSampleL *= mergedCoefficientsL; + + mergedCoefficientsR *= 0.25; + inputSampleR *= mergedCoefficientsR; + + if (outMakeupGainL != 1.0) inputSampleL = inputSampleL * outMakeupGainL; + if (outMakeupGainR != 1.0) inputSampleR = inputSampleR * outMakeupGainR; + + fpFlip = !fpFlip; + + if (output < 1.0) { + inputSampleL *= output; + inputSampleR *= output; + } + + if (wet < 1.0) { + inputSampleL = (inputSampleL * wet) + (drySampleL * (1.0-wet)); + inputSampleR = (inputSampleR * wet) + (drySampleR * (1.0-wet)); + } + + //begin 32 bit stereo floating point dither + int expon; frexpf((float)inputSampleL, &expon); + fpd ^= fpd << 13; fpd ^= fpd >> 17; fpd ^= fpd << 5; + inputSampleL += static_cast(fpd) * 5.960464655174751e-36L * pow(2,expon+62); + frexpf((float)inputSampleR, &expon); + fpd ^= fpd << 13; fpd ^= fpd >> 17; fpd ^= fpd << 5; + inputSampleR += static_cast(fpd) * 5.960464655174751e-36L * pow(2,expon+62); + //end 32 bit stereo floating point dither + + *out1 = inputSampleL; + *out2 = inputSampleR; + + *in1++; + *in2++; + *out1++; + *out2++; + } +} + +void PurestSquish::processDoubleReplacing(double **inputs, double **outputs, VstInt32 sampleFrames) +{ + double* in1 = inputs[0]; + double* in2 = inputs[1]; + double* out1 = outputs[0]; + double* out2 = outputs[1]; + + double overallscale = 1.0; + overallscale /= 44100.0; + overallscale *= getSampleRate(); + + long double highpassSampleL; + long double highpassSampleBL; + long double highpassSampleR; + long double highpassSampleBR; + double threshold = 1.01 - (1.0-pow(1.0-(A*0.5),4)); + double iirAmount = pow(B,4)/overallscale; + double output = C; + double wet = D; + + while (--sampleFrames >= 0) + { + long double inputSampleL = *in1; + long double inputSampleR = *in2; + + static int noisesourceL = 0; + static int noisesourceR = 850010; + int residue; + double applyresidue; + + noisesourceL = noisesourceL % 1700021; noisesourceL++; + residue = noisesourceL * noisesourceL; + residue = residue % 170003; residue *= residue; + residue = residue % 17011; residue *= residue; + residue = residue % 1709; residue *= residue; + residue = residue % 173; residue *= residue; + residue = residue % 17; + applyresidue = residue; + applyresidue *= 0.00000001; + applyresidue *= 0.00000001; + inputSampleL += applyresidue; + if (inputSampleL<1.2e-38 && -inputSampleL<1.2e-38) { + inputSampleL -= applyresidue; + } + + noisesourceR = noisesourceR % 1700021; noisesourceR++; + residue = noisesourceR * noisesourceR; + residue = residue % 170003; residue *= residue; + residue = residue % 17011; residue *= residue; + residue = residue % 1709; residue *= residue; + residue = residue % 173; residue *= residue; + residue = residue % 17; + applyresidue = residue; + applyresidue *= 0.00000001; + applyresidue *= 0.00000001; + inputSampleR += applyresidue; + if (inputSampleR<1.2e-38 && -inputSampleR<1.2e-38) { + inputSampleR -= applyresidue; + } + //for live air, we always apply the dither noise. Then, if our result is + //effectively digital black, we'll subtract it again. We want a 'air' hiss + long double drySampleL = inputSampleL; + long double drySampleR = inputSampleR; + + double muMakeupGainL = 1.0 / threshold; + double outMakeupGainL = sqrt(muMakeupGainL); + muMakeupGainL += outMakeupGainL; + muMakeupGainL *= 0.5; + outMakeupGainL *= 0.5; + //gain settings around threshold + double releaseL = mergedCoefficientsL * 32768.0; + releaseL /= overallscale; + double fastestL = sqrt(releaseL); + //speed settings around release + double lastCorrectionL = mergedCoefficientsL; + // µ µ µ µ µ µ µ µ µ µ µ µ is the kitten song o/~ + + double muMakeupGainR = 1.0 / threshold; + double outMakeupGainR = sqrt(muMakeupGainR); + muMakeupGainR += outMakeupGainR; + muMakeupGainR *= 0.5; + outMakeupGainR *= 0.5; + //gain settings around threshold + double releaseR = mergedCoefficientsR * 32768.0; + releaseR /= overallscale; + double fastestR = sqrt(releaseR); + //speed settings around release + double lastCorrectionR = mergedCoefficientsR; + // µ µ µ µ µ µ µ µ µ µ µ µ is the kitten song o/~ + + if (muMakeupGainL != 1.0) inputSampleL = inputSampleL * muMakeupGainL; + if (muMakeupGainR != 1.0) inputSampleR = inputSampleR * muMakeupGainR; + + highpassSampleBL = highpassSampleL = inputSampleL; + highpassSampleBR = highpassSampleR = inputSampleR; + + if (count < 1 || count > 3) count = 1; + switch (count) + { + case 1: + //begin L + iirSampleAL = (iirSampleAL * (1 - iirAmount)) + (highpassSampleL * iirAmount); + highpassSampleL -= iirSampleAL; + if (fabs(highpassSampleL) > threshold) + { + muVaryL = threshold / fabs(highpassSampleL); + muAttackL = sqrt(fabs(muSpeedAL)); + muCoefficientAL = muCoefficientAL * (muAttackL-1.0); + if (muVaryL < threshold) + { + muCoefficientAL = muCoefficientAL + threshold; + } + else + { + muCoefficientAL = muCoefficientAL + muVaryL; + } + muCoefficientAL = muCoefficientAL / muAttackL; + } + else + { + muCoefficientAL = muCoefficientAL * ((muSpeedAL * muSpeedAL)-1.0); + muCoefficientAL = muCoefficientAL + 1.0; + muCoefficientAL = muCoefficientAL / (muSpeedAL * muSpeedAL); + } + muNewSpeedL = muSpeedAL * (muSpeedAL-1); + muNewSpeedL = muNewSpeedL + fabs(highpassSampleL*releaseL)+fastestL; + muSpeedAL = muNewSpeedL / muSpeedAL; + lastCoefficientAL = pow(muCoefficientAL,2); + mergedCoefficientsL = lastCoefficientBL; + mergedCoefficientsL += lastCoefficientAL; + lastCoefficientAL *= (1.0-lastCorrectionL); + lastCoefficientAL += (muCoefficientAL * lastCorrectionL); + lastCoefficientBL = lastCoefficientAL; + //end L + + //begin R + iirSampleAR = (iirSampleAR * (1 - iirAmount)) + (highpassSampleR * iirAmount); + highpassSampleR -= iirSampleAR; + if (fabs(highpassSampleR) > threshold) + { + muVaryR = threshold / fabs(highpassSampleR); + muAttackR = sqrt(fabs(muSpeedAR)); + muCoefficientAR = muCoefficientAR * (muAttackR-1.0); + if (muVaryR < threshold) + { + muCoefficientAR = muCoefficientAR + threshold; + } + else + { + muCoefficientAR = muCoefficientAR + muVaryR; + } + muCoefficientAR = muCoefficientAR / muAttackR; + } + else + { + muCoefficientAR = muCoefficientAR * ((muSpeedAR * muSpeedAR)-1.0); + muCoefficientAR = muCoefficientAR + 1.0; + muCoefficientAR = muCoefficientAR / (muSpeedAR * muSpeedAR); + } + muNewSpeedR = muSpeedAR * (muSpeedAR-1); + muNewSpeedR = muNewSpeedR + fabs(highpassSampleR*releaseR)+fastestR; + muSpeedAR = muNewSpeedR / muSpeedAR; + lastCoefficientAR = pow(muCoefficientAR,2); + mergedCoefficientsR = lastCoefficientBR; + mergedCoefficientsR += lastCoefficientAR; + lastCoefficientAR *= (1.0-lastCorrectionR); + lastCoefficientAR += (muCoefficientAR * lastCorrectionR); + lastCoefficientBR = lastCoefficientAR; + //end R + + break; + case 2: + //begin L + iirSampleBL = (iirSampleBL * (1 - iirAmount)) + (highpassSampleL * iirAmount); + highpassSampleL -= iirSampleBL; + if (fabs(highpassSampleL) > threshold) + { + muVaryL = threshold / fabs(highpassSampleL); + muAttackL = sqrt(fabs(muSpeedBL)); + muCoefficientBL = muCoefficientBL * (muAttackL-1); + if (muVaryL < threshold) + { + muCoefficientBL = muCoefficientBL + threshold; + } + else + { + muCoefficientBL = muCoefficientBL + muVaryL; + } + muCoefficientBL = muCoefficientBL / muAttackL; + } + else + { + muCoefficientBL = muCoefficientBL * ((muSpeedBL * muSpeedBL)-1.0); + muCoefficientBL = muCoefficientBL + 1.0; + muCoefficientBL = muCoefficientBL / (muSpeedBL * muSpeedBL); + } + muNewSpeedL = muSpeedBL * (muSpeedBL-1); + muNewSpeedL = muNewSpeedL + fabs(highpassSampleL*releaseL)+fastestL; + muSpeedBL = muNewSpeedL / muSpeedBL; + lastCoefficientAL = pow(muCoefficientBL,2); + mergedCoefficientsL = lastCoefficientBL; + mergedCoefficientsL += lastCoefficientAL; + lastCoefficientAL *= (1.0-lastCorrectionL); + lastCoefficientAL += (muCoefficientBL * lastCorrectionL); + lastCoefficientBL = lastCoefficientAL; + //end L + + //begin R + iirSampleBR = (iirSampleBR * (1 - iirAmount)) + (highpassSampleR * iirAmount); + highpassSampleR -= iirSampleBR; + if (fabs(highpassSampleR) > threshold) + { + muVaryR = threshold / fabs(highpassSampleR); + muAttackR = sqrt(fabs(muSpeedBR)); + muCoefficientBR = muCoefficientBR * (muAttackR-1); + if (muVaryR < threshold) + { + muCoefficientBR = muCoefficientBR + threshold; + } + else + { + muCoefficientBR = muCoefficientBR + muVaryR; + } + muCoefficientBR = muCoefficientBR / muAttackR; + } + else + { + muCoefficientBR = muCoefficientBR * ((muSpeedBR * muSpeedBR)-1.0); + muCoefficientBR = muCoefficientBR + 1.0; + muCoefficientBR = muCoefficientBR / (muSpeedBR * muSpeedBR); + } + muNewSpeedR = muSpeedBR * (muSpeedBR-1); + muNewSpeedR = muNewSpeedR + fabs(highpassSampleR*releaseR)+fastestR; + muSpeedBR = muNewSpeedR / muSpeedBR; + lastCoefficientAR = pow(muCoefficientBR,2); + mergedCoefficientsR = lastCoefficientBR; + mergedCoefficientsR += lastCoefficientAR; + lastCoefficientAR *= (1.0-lastCorrectionR); + lastCoefficientAR += (muCoefficientBR * lastCorrectionR); + lastCoefficientBR = lastCoefficientAR; + //end R + + break; + case 3: + //begin L + iirSampleCL = (iirSampleCL * (1 - iirAmount)) + (highpassSampleL * iirAmount); + highpassSampleL -= iirSampleCL; + if (fabs(highpassSampleL) > threshold) + { + muVaryL = threshold / fabs(highpassSampleL); + muAttackL = sqrt(fabs(muSpeedCL)); + muCoefficientCL = muCoefficientCL * (muAttackL-1); + if (muVaryL < threshold) + { + muCoefficientCL = muCoefficientCL + threshold; + } + else + { + muCoefficientCL = muCoefficientCL + muVaryL; + } + muCoefficientCL = muCoefficientCL / muAttackL; + } + else + { + muCoefficientCL = muCoefficientCL * ((muSpeedCL * muSpeedCL)-1.0); + muCoefficientCL = muCoefficientCL + 1.0; + muCoefficientCL = muCoefficientCL / (muSpeedCL * muSpeedCL); + } + muNewSpeedL = muSpeedCL * (muSpeedCL-1); + muNewSpeedL = muNewSpeedL + fabs(highpassSampleL*releaseL)+fastestL; + muSpeedCL = muNewSpeedL / muSpeedCL; + lastCoefficientAL = pow(muCoefficientCL,2); + mergedCoefficientsL = lastCoefficientBL; + mergedCoefficientsL += lastCoefficientAL; + lastCoefficientAL *= (1.0-lastCorrectionL); + lastCoefficientAL += (muCoefficientCL * lastCorrectionL); + lastCoefficientBL = lastCoefficientAL; + //end L + + //begin R + iirSampleCR = (iirSampleCR * (1 - iirAmount)) + (highpassSampleR * iirAmount); + highpassSampleR -= iirSampleCR; + if (fabs(highpassSampleR) > threshold) + { + muVaryR = threshold / fabs(highpassSampleR); + muAttackR = sqrt(fabs(muSpeedCR)); + muCoefficientCR = muCoefficientCR * (muAttackR-1); + if (muVaryR < threshold) + { + muCoefficientCR = muCoefficientCR + threshold; + } + else + { + muCoefficientCR = muCoefficientCR + muVaryR; + } + muCoefficientCR = muCoefficientCR / muAttackR; + } + else + { + muCoefficientCR = muCoefficientCR * ((muSpeedCR * muSpeedCR)-1.0); + muCoefficientCR = muCoefficientCR + 1.0; + muCoefficientCR = muCoefficientCR / (muSpeedCR * muSpeedCR); + } + muNewSpeedR = muSpeedCR * (muSpeedCR-1); + muNewSpeedR = muNewSpeedR + fabs(highpassSampleR*releaseR)+fastestR; + muSpeedCR = muNewSpeedR / muSpeedCR; + lastCoefficientAR = pow(muCoefficientCR,2); + mergedCoefficientsR = lastCoefficientBR; + mergedCoefficientsR += lastCoefficientAR; + lastCoefficientAR *= (1.0-lastCorrectionR); + lastCoefficientAR += (muCoefficientCR * lastCorrectionR); + lastCoefficientBR = lastCoefficientAR; + //end R + + break; + } + count++; + + //applied compression with vari-vari-µ-µ-µ-µ-µ-µ-is-the-kitten-song o/~ + //applied gain correction to control output level- tends to constrain sound rather than inflate it + + if (fpFlip) { + //begin L + iirSampleDL = (iirSampleDL * (1 - iirAmount)) + (highpassSampleBL * iirAmount); + highpassSampleBL -= iirSampleDL; + if (fabs(highpassSampleBL) > threshold) + { + muVaryL = threshold / fabs(highpassSampleBL); + muAttackL = sqrt(fabs(muSpeedDL)); + muCoefficientDL = muCoefficientDL * (muAttackL-1.0); + if (muVaryL < threshold) + { + muCoefficientDL = muCoefficientDL + threshold; + } + else + { + muCoefficientDL = muCoefficientDL + muVaryL; + } + muCoefficientDL = muCoefficientDL / muAttackL; + } + else + { + muCoefficientDL = muCoefficientDL * ((muSpeedDL * muSpeedDL)-1.0); + muCoefficientDL = muCoefficientDL + 1.0; + muCoefficientDL = muCoefficientDL / (muSpeedDL * muSpeedDL); + } + muNewSpeedL = muSpeedDL * (muSpeedDL-1); + muNewSpeedL = muNewSpeedL + fabs(highpassSampleBL*releaseL)+fastestL; + muSpeedDL = muNewSpeedL / muSpeedDL; + lastCoefficientCL = pow(muCoefficientEL,2); + mergedCoefficientsL += lastCoefficientDL; + mergedCoefficientsL += lastCoefficientCL; + lastCoefficientCL *= (1.0-lastCorrectionL); + lastCoefficientCL += (muCoefficientDL * lastCorrectionL); + lastCoefficientDL = lastCoefficientCL; + //end L + + //begin R + iirSampleDR = (iirSampleDR * (1 - iirAmount)) + (highpassSampleBR * iirAmount); + highpassSampleBR -= iirSampleDR; + if (fabs(highpassSampleBR) > threshold) + { + muVaryR = threshold / fabs(highpassSampleBR); + muAttackR = sqrt(fabs(muSpeedDR)); + muCoefficientDR = muCoefficientDR * (muAttackR-1.0); + if (muVaryR < threshold) + { + muCoefficientDR = muCoefficientDR + threshold; + } + else + { + muCoefficientDR = muCoefficientDR + muVaryR; + } + muCoefficientDR = muCoefficientDR / muAttackR; + } + else + { + muCoefficientDR = muCoefficientDR * ((muSpeedDR * muSpeedDR)-1.0); + muCoefficientDR = muCoefficientDR + 1.0; + muCoefficientDR = muCoefficientDR / (muSpeedDR * muSpeedDR); + } + muNewSpeedR = muSpeedDR * (muSpeedDR-1); + muNewSpeedR = muNewSpeedR + fabs(highpassSampleBR*releaseR)+fastestR; + muSpeedDR = muNewSpeedR / muSpeedDR; + lastCoefficientCR = pow(muCoefficientER,2); + mergedCoefficientsR += lastCoefficientDR; + mergedCoefficientsR += lastCoefficientCR; + lastCoefficientCR *= (1.0-lastCorrectionR); + lastCoefficientCR += (muCoefficientDR * lastCorrectionR); + lastCoefficientDR = lastCoefficientCR; + //end R + + } else { + //begin L + iirSampleEL = (iirSampleEL * (1 - iirAmount)) + (highpassSampleBL * iirAmount); + highpassSampleBL -= iirSampleEL; + if (fabs(highpassSampleBL) > threshold) + { + muVaryL = threshold / fabs(highpassSampleBL); + muAttackL = sqrt(fabs(muSpeedEL)); + muCoefficientEL = muCoefficientEL * (muAttackL-1.0); + if (muVaryL < threshold) + { + muCoefficientEL = muCoefficientEL + threshold; + } + else + { + muCoefficientEL = muCoefficientEL + muVaryL; + } + muCoefficientEL = muCoefficientEL / muAttackL; + } + else + { + muCoefficientEL = muCoefficientEL * ((muSpeedEL * muSpeedEL)-1.0); + muCoefficientEL = muCoefficientEL + 1.0; + muCoefficientEL = muCoefficientEL / (muSpeedEL * muSpeedEL); + } + muNewSpeedL = muSpeedEL * (muSpeedEL-1); + muNewSpeedL = muNewSpeedL + fabs(highpassSampleBL*releaseL)+fastestL; + muSpeedEL = muNewSpeedL / muSpeedEL; + lastCoefficientCL = pow(muCoefficientEL,2); + mergedCoefficientsL += lastCoefficientDL; + mergedCoefficientsL += lastCoefficientCL; + lastCoefficientCL *= (1.0-lastCorrectionL); + lastCoefficientCL += (muCoefficientEL * lastCorrectionL); + lastCoefficientDL = lastCoefficientCL; + //end L + + //begin R + iirSampleER = (iirSampleER * (1 - iirAmount)) + (highpassSampleBR * iirAmount); + highpassSampleBR -= iirSampleER; + if (fabs(highpassSampleBR) > threshold) + { + muVaryR = threshold / fabs(highpassSampleBR); + muAttackR = sqrt(fabs(muSpeedER)); + muCoefficientER = muCoefficientER * (muAttackR-1.0); + if (muVaryR < threshold) + { + muCoefficientER = muCoefficientER + threshold; + } + else + { + muCoefficientER = muCoefficientER + muVaryR; + } + muCoefficientER = muCoefficientER / muAttackR; + } + else + { + muCoefficientER = muCoefficientER * ((muSpeedER * muSpeedER)-1.0); + muCoefficientER = muCoefficientER + 1.0; + muCoefficientER = muCoefficientER / (muSpeedER * muSpeedER); + } + muNewSpeedR = muSpeedER * (muSpeedER-1); + muNewSpeedR = muNewSpeedR + fabs(highpassSampleBR*releaseR)+fastestR; + muSpeedER = muNewSpeedR / muSpeedER; + lastCoefficientCR = pow(muCoefficientER,2); + mergedCoefficientsR += lastCoefficientDR; + mergedCoefficientsR += lastCoefficientCR; + lastCoefficientCR *= (1.0-lastCorrectionR); + lastCoefficientCR += (muCoefficientER * lastCorrectionR); + lastCoefficientDR = lastCoefficientCR; + //end R + + } + mergedCoefficientsL *= 0.25; + inputSampleL *= mergedCoefficientsL; + + mergedCoefficientsR *= 0.25; + inputSampleR *= mergedCoefficientsR; + + if (outMakeupGainL != 1.0) inputSampleL = inputSampleL * outMakeupGainL; + if (outMakeupGainR != 1.0) inputSampleR = inputSampleR * outMakeupGainR; + + fpFlip = !fpFlip; + + if (output < 1.0) { + inputSampleL *= output; + inputSampleR *= output; + } + + if (wet < 1.0) { + inputSampleL = (inputSampleL * wet) + (drySampleL * (1.0-wet)); + inputSampleR = (inputSampleR * wet) + (drySampleR * (1.0-wet)); + } + + //begin 64 bit stereo floating point dither + int expon; frexp((double)inputSampleL, &expon); + fpd ^= fpd << 13; fpd ^= fpd >> 17; fpd ^= fpd << 5; + inputSampleL += static_cast(fpd) * 1.110223024625156e-44L * pow(2,expon+62); + frexp((double)inputSampleR, &expon); + fpd ^= fpd << 13; fpd ^= fpd >> 17; fpd ^= fpd << 5; + inputSampleR += static_cast(fpd) * 1.110223024625156e-44L * pow(2,expon+62); + //end 64 bit stereo floating point dither + + *out1 = inputSampleL; + *out2 = inputSampleR; + + *in1++; + *in2++; + *out1++; + *out2++; + } +} -- cgit v1.2.3