From 95d0811f50332fe99cdb52eed0bd929faec2d43c Mon Sep 17 00:00:00 2001 From: Chris Johnson Date: Sun, 26 Jan 2020 17:47:49 -0500 Subject: Tape. Just 'Tape'. --- plugins/LinuxVST/src/Tape/Tape.cpp | 133 ++++++++ plugins/LinuxVST/src/Tape/Tape.h | 85 +++++ plugins/LinuxVST/src/Tape/TapeProc.cpp | 564 +++++++++++++++++++++++++++++++++ 3 files changed, 782 insertions(+) create mode 100755 plugins/LinuxVST/src/Tape/Tape.cpp create mode 100755 plugins/LinuxVST/src/Tape/Tape.h create mode 100755 plugins/LinuxVST/src/Tape/TapeProc.cpp (limited to 'plugins/LinuxVST/src/Tape') diff --git a/plugins/LinuxVST/src/Tape/Tape.cpp b/plugins/LinuxVST/src/Tape/Tape.cpp new file mode 100755 index 0000000..eb3712c --- /dev/null +++ b/plugins/LinuxVST/src/Tape/Tape.cpp @@ -0,0 +1,133 @@ +/* ======================================== + * Tape - Tape.h + * Copyright (c) 2016 airwindows, All rights reserved + * ======================================== */ + +#ifndef __Tape_H +#include "Tape.h" +#endif + +AudioEffect* createEffectInstance(audioMasterCallback audioMaster) {return new Tape(audioMaster);} + +Tape::Tape(audioMasterCallback audioMaster) : + AudioEffectX(audioMaster, kNumPrograms, kNumParameters) +{ + A = 0.5; + iirMidRollerAL = 0.0; + iirMidRollerBL = 0.0; + iirHeadBumpAL = 0.0; + iirHeadBumpBL = 0.0; + iirMidRollerAR = 0.0; + iirMidRollerBR = 0.0; + iirHeadBumpAR = 0.0; + iirHeadBumpBR = 0.0; + for (int x = 0; x < 9; x++) { + biquadAL[x] = 0.0;biquadBL[x] = 0.0;biquadCL[x] = 0.0;biquadDL[x] = 0.0; + biquadAR[x] = 0.0;biquadBR[x] = 0.0;biquadCR[x] = 0.0;biquadDR[x] = 0.0; + } + flip = false; + lastSampleL = 0.0; + lastSampleR = 0.0; + 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 +} + +Tape::~Tape() {} +VstInt32 Tape::getVendorVersion () {return 1000;} +void Tape::setProgramName(char *name) {vst_strncpy (_programName, name, kVstMaxProgNameLen);} +void Tape::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 Tape::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 Tape::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 Tape::setParameter(VstInt32 index, float value) { + switch (index) { + case kParamA: A = value; break; + default: throw; // unknown parameter, shouldn't happen! + } +} + +float Tape::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 Tape::getParameterName(VstInt32 index, char *text) { + switch (index) { + case kParamA: vst_strncpy (text, "Slam", kVstMaxParamStrLen); break; + default: break; // unknown parameter, shouldn't happen! + } //this is our labels for displaying in the VST host +} + +void Tape::getParameterDisplay(VstInt32 index, char *text) { + switch (index) { + case kParamA: float2string ((A-0.5)*24.0, text, kVstMaxParamStrLen); break; + default: break; // unknown parameter, shouldn't happen! + } //this displays the values and handles 'popups' where it's discrete choices +} + +void Tape::getParameterLabel(VstInt32 index, char *text) { + switch (index) { + case kParamA: vst_strncpy (text, "dB", kVstMaxParamStrLen); break; + default: break; // unknown parameter, shouldn't happen! + } +} + +VstInt32 Tape::canDo(char *text) +{ return (_canDo.find(text) == _canDo.end()) ? -1: 1; } // 1 = yes, -1 = no, 0 = don't know + +bool Tape::getEffectName(char* name) { + vst_strncpy(name, "Tape", kVstMaxProductStrLen); return true; +} + +VstPlugCategory Tape::getPlugCategory() {return kPlugCategEffect;} + +bool Tape::getProductString(char* text) { + vst_strncpy (text, "airwindows Tape", kVstMaxProductStrLen); return true; +} + +bool Tape::getVendorString(char* text) { + vst_strncpy (text, "airwindows", kVstMaxVendorStrLen); return true; +} diff --git a/plugins/LinuxVST/src/Tape/Tape.h b/plugins/LinuxVST/src/Tape/Tape.h new file mode 100755 index 0000000..6c1dc00 --- /dev/null +++ b/plugins/LinuxVST/src/Tape/Tape.h @@ -0,0 +1,85 @@ +/* ======================================== + * Tape - Tape.h + * Created 8/12/11 by SPIAdmin + * Copyright (c) 2011 __MyCompanyName__, All rights reserved + * ======================================== */ + +#ifndef __Tape_H +#define __Tape_H + +#ifndef __audioeffect__ +#include "audioeffectx.h" +#endif + +#include +#include +#include + +enum { + kParamA = 0, + kNumParameters = 1 +}; // + +const int kNumPrograms = 0; +const int kNumInputs = 2; +const int kNumOutputs = 2; +const unsigned long kUniqueId = 'tape'; //Change this to what the AU identity is! + +class Tape : + public AudioEffectX +{ +public: + Tape(audioMasterCallback audioMaster); + ~Tape(); + virtual bool getEffectName(char* name); // The plug-in name + virtual VstPlugCategory getPlugCategory(); // The general category for the plug-in + virtual bool getProductString(char* text); // This is a unique plug-in string provided by Steinberg + virtual bool getVendorString(char* text); // Vendor info + virtual VstInt32 getVendorVersion(); // Version number + virtual void processReplacing (float** inputs, float** outputs, VstInt32 sampleFrames); + virtual void processDoubleReplacing (double** inputs, double** outputs, VstInt32 sampleFrames); + virtual void getProgramName(char *name); // read the name from the host + virtual void setProgramName(char *name); // changes the name of the preset displayed in the host + virtual VstInt32 getChunk (void** data, bool isPreset); + virtual VstInt32 setChunk (void* data, VstInt32 byteSize, bool isPreset); + virtual float getParameter(VstInt32 index); // get the parameter value at the specified index + virtual void setParameter(VstInt32 index, float value); // set the parameter at index to value + virtual void getParameterLabel(VstInt32 index, char *text); // label for the parameter (eg dB) + virtual void getParameterName(VstInt32 index, char *text); // name of the parameter + virtual void getParameterDisplay(VstInt32 index, char *text); // text description of the current value + virtual VstInt32 canDo(char *text); +private: + char _programName[kVstMaxProgNameLen + 1]; + std::set< std::string > _canDo; + + double iirMidRollerAL; + double iirMidRollerBL; + double iirHeadBumpAL; + double iirHeadBumpBL; + + double iirMidRollerAR; + double iirMidRollerBR; + double iirHeadBumpAR; + double iirHeadBumpBR; + + long double biquadAL[9]; + long double biquadBL[9]; + long double biquadCL[9]; + long double biquadDL[9]; + + long double biquadAR[9]; + long double biquadBR[9]; + long double biquadCR[9]; + long double biquadDR[9]; + bool flip; + + long double lastSampleL; + long double lastSampleR; + + uint32_t fpd; + //default stuff + + float A; +}; + +#endif diff --git a/plugins/LinuxVST/src/Tape/TapeProc.cpp b/plugins/LinuxVST/src/Tape/TapeProc.cpp new file mode 100755 index 0000000..1d1050c --- /dev/null +++ b/plugins/LinuxVST/src/Tape/TapeProc.cpp @@ -0,0 +1,564 @@ +/* ======================================== + * Tape - Tape.h + * Copyright (c) 2016 airwindows, All rights reserved + * ======================================== */ + +#ifndef __Tape_H +#include "Tape.h" +#endif + +void Tape::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(); + + double inputgain = pow(10.0,((A-0.5)*24.0)/20.0); + double HeadBumpFreq = 0.12/overallscale; + double softness = 0.618033988749894848204586; + double RollAmount = (1.0 - softness) / overallscale; + //[0] is frequency: 0.000001 to 0.499999 is near-zero to near-Nyquist + //[1] is resonance, 0.7071 is Butterworth. Also can't be zero + biquadAL[0] = biquadBL[0] = biquadAR[0] = biquadBR[0] = 0.0072/overallscale; + biquadAL[1] = biquadBL[1] = biquadAR[1] = biquadBR[1] = 0.0009; + double K = tan(M_PI * biquadBR[0]); + double norm = 1.0 / (1.0 + K / biquadBR[1] + K * K); + biquadAL[2] = biquadBL[2] = biquadAR[2] = biquadBR[2] = K / biquadBR[1] * norm; + biquadAL[4] = biquadBL[4] = biquadAR[4] = biquadBR[4] = -biquadBR[2]; + biquadAL[5] = biquadBL[5] = biquadAR[5] = biquadBR[5] = 2.0 * (K * K - 1.0) * norm; + biquadAL[6] = biquadBL[6] = biquadAR[6] = biquadBR[6] = (1.0 - K / biquadBR[1] + K * K) * norm; + + biquadCL[0] = biquadDL[0] = biquadCR[0] = biquadDR[0] = 0.032/overallscale; + biquadCL[1] = biquadDL[1] = biquadCR[1] = biquadDR[1] = 0.0007; + K = tan(M_PI * biquadDR[0]); + norm = 1.0 / (1.0 + K / biquadDR[1] + K * K); + biquadCL[2] = biquadDL[2] = biquadCR[2] = biquadDR[2] = K / biquadDR[1] * norm; + biquadCL[4] = biquadDL[4] = biquadCR[4] = biquadDR[4] = -biquadDR[2]; + biquadCL[5] = biquadDL[5] = biquadCR[5] = biquadDR[5] = 2.0 * (K * K - 1.0) * norm; + biquadCL[6] = biquadDL[6] = biquadCR[6] = biquadDR[6] = (1.0 - K / biquadDR[1] + K * K) * norm; + + while (--sampleFrames >= 0) + { + long double inputSampleL = *in1; + long double inputSampleR = *in2; + if (fabs(inputSampleL)<1.18e-37) inputSampleL = fpd * 1.18e-37; + if (fabs(inputSampleR)<1.18e-37) inputSampleR = fpd * 1.18e-37; + long double drySampleL = inputSampleL; + long double drySampleR = inputSampleR; + + long double HighsSampleL = 0.0; + long double HighsSampleR = 0.0; + long double NonHighsSampleL = 0.0; + long double NonHighsSampleR = 0.0; + long double tempSample; + + if (flip) + { + iirMidRollerAL = (iirMidRollerAL * (1.0 - RollAmount)) + (inputSampleL * RollAmount); + iirMidRollerAR = (iirMidRollerAR * (1.0 - RollAmount)) + (inputSampleR * RollAmount); + HighsSampleL = inputSampleL - iirMidRollerAL; + HighsSampleR = inputSampleR - iirMidRollerAR; + NonHighsSampleL = iirMidRollerAL; + NonHighsSampleR = iirMidRollerAR; + + iirHeadBumpAL += (inputSampleL * 0.05); + iirHeadBumpAR += (inputSampleR * 0.05); + iirHeadBumpAL -= (iirHeadBumpAL * iirHeadBumpAL * iirHeadBumpAL * HeadBumpFreq); + iirHeadBumpAR -= (iirHeadBumpAR * iirHeadBumpAR * iirHeadBumpAR * HeadBumpFreq); + iirHeadBumpAL = sin(iirHeadBumpAL); + iirHeadBumpAR = sin(iirHeadBumpAR); + + tempSample = (iirHeadBumpAL * biquadAL[2]) + biquadAL[7]; + biquadAL[7] = (iirHeadBumpAL * biquadAL[3]) - (tempSample * biquadAL[5]) + biquadAL[8]; + biquadAL[8] = (iirHeadBumpAL * biquadAL[4]) - (tempSample * biquadAL[6]); + iirHeadBumpAL = tempSample; //interleaved biquad + if (iirHeadBumpAL > 1.0) iirHeadBumpAL = 1.0; + if (iirHeadBumpAL < -1.0) iirHeadBumpAL = -1.0; + iirHeadBumpAL = asin(iirHeadBumpAL); + + tempSample = (iirHeadBumpAR * biquadAR[2]) + biquadAR[7]; + biquadAR[7] = (iirHeadBumpAR * biquadAR[3]) - (tempSample * biquadAR[5]) + biquadAR[8]; + biquadAR[8] = (iirHeadBumpAR * biquadAR[4]) - (tempSample * biquadAR[6]); + iirHeadBumpAR = tempSample; //interleaved biquad + if (iirHeadBumpAR > 1.0) iirHeadBumpAR = 1.0; + if (iirHeadBumpAR < -1.0) iirHeadBumpAR = -1.0; + iirHeadBumpAR = asin(iirHeadBumpAR); + + inputSampleL = sin(inputSampleL); + tempSample = (inputSampleL * biquadCL[2]) + biquadCL[7]; + biquadCL[7] = (inputSampleL * biquadCL[3]) - (tempSample * biquadCL[5]) + biquadCL[8]; + biquadCL[8] = (inputSampleL * biquadCL[4]) - (tempSample * biquadCL[6]); + inputSampleL = tempSample; //interleaved biquad + if (inputSampleL > 1.0) inputSampleL = 1.0; + if (inputSampleL < -1.0) inputSampleL = -1.0; + inputSampleL = asin(inputSampleL); + + inputSampleR = sin(inputSampleR); + tempSample = (inputSampleR * biquadCR[2]) + biquadCR[7]; + biquadCR[7] = (inputSampleR * biquadCR[3]) - (tempSample * biquadCR[5]) + biquadCR[8]; + biquadCR[8] = (inputSampleR * biquadCR[4]) - (tempSample * biquadCR[6]); + inputSampleR = tempSample; //interleaved biquad + if (inputSampleR > 1.0) inputSampleR = 1.0; + if (inputSampleR < -1.0) inputSampleR = -1.0; + inputSampleR = asin(inputSampleR); + } else { + iirMidRollerBL = (iirMidRollerBL * (1.0 - RollAmount)) + (inputSampleL * RollAmount); + iirMidRollerBR = (iirMidRollerBR * (1.0 - RollAmount)) + (inputSampleR * RollAmount); + HighsSampleL = inputSampleL - iirMidRollerBL; + HighsSampleR = inputSampleR - iirMidRollerBR; + NonHighsSampleL = iirMidRollerBL; + NonHighsSampleR = iirMidRollerBR; + + iirHeadBumpBL += (inputSampleL * 0.05); + iirHeadBumpBR += (inputSampleR * 0.05); + iirHeadBumpBL -= (iirHeadBumpBL * iirHeadBumpBL * iirHeadBumpBL * HeadBumpFreq); + iirHeadBumpBR -= (iirHeadBumpBR * iirHeadBumpBR * iirHeadBumpBR * HeadBumpFreq); + iirHeadBumpBL = sin(iirHeadBumpBL); + iirHeadBumpBR = sin(iirHeadBumpBR); + + tempSample = (iirHeadBumpBL * biquadBL[2]) + biquadBL[7]; + biquadBL[7] = (iirHeadBumpBL * biquadBL[3]) - (tempSample * biquadBL[5]) + biquadBL[8]; + biquadBL[8] = (iirHeadBumpBL * biquadBL[4]) - (tempSample * biquadBL[6]); + iirHeadBumpBL = tempSample; //interleaved biquad + if (iirHeadBumpBL > 1.0) iirHeadBumpBL = 1.0; + if (iirHeadBumpBL < -1.0) iirHeadBumpBL = -1.0; + iirHeadBumpBL = asin(iirHeadBumpBL); + + tempSample = (iirHeadBumpBR * biquadBR[2]) + biquadBR[7]; + biquadBR[7] = (iirHeadBumpBR * biquadBR[3]) - (tempSample * biquadBR[5]) + biquadBR[8]; + biquadBR[8] = (iirHeadBumpBR * biquadBR[4]) - (tempSample * biquadBR[6]); + iirHeadBumpBR = tempSample; //interleaved biquad + if (iirHeadBumpBR > 1.0) iirHeadBumpBR = 1.0; + if (iirHeadBumpBR < -1.0) iirHeadBumpBR = -1.0; + iirHeadBumpBR = asin(iirHeadBumpBR); + + inputSampleL = sin(inputSampleL); + tempSample = (inputSampleL * biquadDL[2]) + biquadDL[7]; + biquadDL[7] = (inputSampleL * biquadDL[3]) - (tempSample * biquadDL[5]) + biquadDL[8]; + biquadDL[8] = (inputSampleL * biquadDL[4]) - (tempSample * biquadDL[6]); + inputSampleL = tempSample; //interleaved biquad + if (inputSampleL > 1.0) inputSampleL = 1.0; + if (inputSampleL < -1.0) inputSampleL = -1.0; + inputSampleL = asin(inputSampleL); + + inputSampleR = sin(inputSampleR); + tempSample = (inputSampleR * biquadDR[2]) + biquadDR[7]; + biquadDR[7] = (inputSampleR * biquadDR[3]) - (tempSample * biquadDR[5]) + biquadDR[8]; + biquadDR[8] = (inputSampleR * biquadDR[4]) - (tempSample * biquadDR[6]); + inputSampleR = tempSample; //interleaved biquad + if (inputSampleR > 1.0) inputSampleR = 1.0; + if (inputSampleR < -1.0) inputSampleR = -1.0; + inputSampleR = asin(inputSampleR); + } + flip = !flip; + + long double groundSampleL = drySampleL - inputSampleL; //set up UnBox + long double groundSampleR = drySampleR - inputSampleR; //set up UnBox + + if (inputgain != 1.0) { + inputSampleL *= inputgain; + inputSampleR *= inputgain; + } //gain boost inside UnBox: do not boost fringe audio + + long double applySoften = fabs(HighsSampleL)*1.57079633; + if (applySoften > 1.57079633) applySoften = 1.57079633; + applySoften = 1-cos(applySoften); + if (HighsSampleL > 0) inputSampleL -= applySoften; + if (HighsSampleL < 0) inputSampleL += applySoften; + //apply Soften depending on polarity + applySoften = fabs(HighsSampleR)*1.57079633; + if (applySoften > 1.57079633) applySoften = 1.57079633; + applySoften = 1-cos(applySoften); + if (HighsSampleR > 0) inputSampleR -= applySoften; + if (HighsSampleR < 0) inputSampleR += applySoften; + //apply Soften depending on polarity + + if (inputSampleL > 1.2533141373155) inputSampleL = 1.2533141373155; + if (inputSampleL < -1.2533141373155) inputSampleL = -1.2533141373155; + //clip to 1.2533141373155 to reach maximum output + inputSampleL = sin(inputSampleL * fabs(inputSampleL)) / ((inputSampleL == 0.0) ?1:fabs(inputSampleL)); + //Spiral, for cleanest most optimal tape effect + if (inputSampleR > 1.2533141373155) inputSampleR = 1.2533141373155; + if (inputSampleR < -1.2533141373155) inputSampleR = -1.2533141373155; + //clip to 1.2533141373155 to reach maximum output + inputSampleR = sin(inputSampleR * fabs(inputSampleR)) / ((inputSampleR == 0.0) ?1:fabs(inputSampleR)); + //Spiral, for cleanest most optimal tape effect + + double suppress = (1.0-fabs(inputSampleL)) * 0.00013; + if (iirHeadBumpAL > suppress) iirHeadBumpAL -= suppress; + if (iirHeadBumpAL < -suppress) iirHeadBumpAL += suppress; + if (iirHeadBumpBL > suppress) iirHeadBumpBL -= suppress; + if (iirHeadBumpBL < -suppress) iirHeadBumpBL += suppress; + //restrain resonant quality of head bump algorithm + suppress = (1.0-fabs(inputSampleR)) * 0.00013; + if (iirHeadBumpAR > suppress) iirHeadBumpAR -= suppress; + if (iirHeadBumpAR < -suppress) iirHeadBumpAR += suppress; + if (iirHeadBumpBR > suppress) iirHeadBumpBR -= suppress; + if (iirHeadBumpBR < -suppress) iirHeadBumpBR += suppress; + //restrain resonant quality of head bump algorithm + + inputSampleL += groundSampleL; //apply UnBox processing + inputSampleR += groundSampleR; //apply UnBox processing + + inputSampleL += ((iirHeadBumpAL + iirHeadBumpBL) * 0.1);//and head bump + inputSampleR += ((iirHeadBumpAR + iirHeadBumpBR) * 0.1);//and head bump + + if (lastSampleL >= 0.99) + { + if (inputSampleL < 0.99) lastSampleL = ((0.99*softness) + (inputSampleL * (1.0-softness))); + else lastSampleL = 0.99; + } + + if (lastSampleL <= -0.99) + { + if (inputSampleL > -0.99) lastSampleL = ((-0.99*softness) + (inputSampleL * (1.0-softness))); + else lastSampleL = -0.99; + } + + if (inputSampleL > 0.99) + { + if (lastSampleL < 0.99) inputSampleL = ((0.99*softness) + (lastSampleL * (1.0-softness))); + else inputSampleL = 0.99; + } + + if (inputSampleL < -0.99) + { + if (lastSampleL > -0.99) inputSampleL = ((-0.99*softness) + (lastSampleL * (1.0-softness))); + else inputSampleL = -0.99; + } + lastSampleL = inputSampleL; //end ADClip L + + + if (lastSampleR >= 0.99) + { + if (inputSampleR < 0.99) lastSampleR = ((0.99*softness) + (inputSampleR * (1.0-softness))); + else lastSampleR = 0.99; + } + + if (lastSampleR <= -0.99) + { + if (inputSampleR > -0.99) lastSampleR = ((-0.99*softness) + (inputSampleR * (1.0-softness))); + else lastSampleR = -0.99; + } + + if (inputSampleR > 0.99) + { + if (lastSampleR < 0.99) inputSampleR = ((0.99*softness) + (lastSampleR * (1.0-softness))); + else inputSampleR = 0.99; + } + + if (inputSampleR < -0.99) + { + if (lastSampleR > -0.99) inputSampleR = ((-0.99*softness) + (lastSampleR * (1.0-softness))); + else inputSampleR = -0.99; + } + lastSampleR = inputSampleR; //end ADClip R + + if (inputSampleL > 0.99) inputSampleL = 0.99; + if (inputSampleL < -0.99) inputSampleL = -0.99; + //final iron bar + if (inputSampleR > 0.99) inputSampleR = 0.99; + if (inputSampleR < -0.99) inputSampleR = -0.99; + //final iron bar + + //begin 32 bit stereo floating point dither + int expon; frexpf((float)inputSampleL, &expon); + fpd ^= fpd << 13; fpd ^= fpd >> 17; fpd ^= fpd << 5; + inputSampleL += ((double(fpd)-uint32_t(0x7fffffff)) * 5.5e-36l * pow(2,expon+62)); + frexpf((float)inputSampleR, &expon); + fpd ^= fpd << 13; fpd ^= fpd >> 17; fpd ^= fpd << 5; + inputSampleR += ((double(fpd)-uint32_t(0x7fffffff)) * 5.5e-36l * pow(2,expon+62)); + //end 32 bit stereo floating point dither + + *out1 = inputSampleL; + *out2 = inputSampleR; + + *in1++; + *in2++; + *out1++; + *out2++; + } +} + +void Tape::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(); + + double inputgain = pow(10.0,((A-0.5)*24.0)/20.0); + double HeadBumpFreq = 0.12/overallscale; + double softness = 0.618033988749894848204586; + double RollAmount = (1.0 - softness) / overallscale; + //[0] is frequency: 0.000001 to 0.499999 is near-zero to near-Nyquist + //[1] is resonance, 0.7071 is Butterworth. Also can't be zero + biquadAL[0] = biquadBL[0] = biquadAR[0] = biquadBR[0] = 0.0072/overallscale; + biquadAL[1] = biquadBL[1] = biquadAR[1] = biquadBR[1] = 0.0009; + double K = tan(M_PI * biquadBR[0]); + double norm = 1.0 / (1.0 + K / biquadBR[1] + K * K); + biquadAL[2] = biquadBL[2] = biquadAR[2] = biquadBR[2] = K / biquadBR[1] * norm; + biquadAL[4] = biquadBL[4] = biquadAR[4] = biquadBR[4] = -biquadBR[2]; + biquadAL[5] = biquadBL[5] = biquadAR[5] = biquadBR[5] = 2.0 * (K * K - 1.0) * norm; + biquadAL[6] = biquadBL[6] = biquadAR[6] = biquadBR[6] = (1.0 - K / biquadBR[1] + K * K) * norm; + + biquadCL[0] = biquadDL[0] = biquadCR[0] = biquadDR[0] = 0.032/overallscale; + biquadCL[1] = biquadDL[1] = biquadCR[1] = biquadDR[1] = 0.0007; + K = tan(M_PI * biquadDR[0]); + norm = 1.0 / (1.0 + K / biquadDR[1] + K * K); + biquadCL[2] = biquadDL[2] = biquadCR[2] = biquadDR[2] = K / biquadDR[1] * norm; + biquadCL[4] = biquadDL[4] = biquadCR[4] = biquadDR[4] = -biquadDR[2]; + biquadCL[5] = biquadDL[5] = biquadCR[5] = biquadDR[5] = 2.0 * (K * K - 1.0) * norm; + biquadCL[6] = biquadDL[6] = biquadCR[6] = biquadDR[6] = (1.0 - K / biquadDR[1] + K * K) * norm; + + while (--sampleFrames >= 0) + { + long double inputSampleL = *in1; + long double inputSampleR = *in2; + if (fabs(inputSampleL)<1.18e-43) inputSampleL = fpd * 1.18e-43; + if (fabs(inputSampleR)<1.18e-43) inputSampleR = fpd * 1.18e-43; + long double drySampleL = inputSampleL; + long double drySampleR = inputSampleR; + + long double HighsSampleL = 0.0; + long double HighsSampleR = 0.0; + long double NonHighsSampleL = 0.0; + long double NonHighsSampleR = 0.0; + long double tempSample; + + if (flip) + { + iirMidRollerAL = (iirMidRollerAL * (1.0 - RollAmount)) + (inputSampleL * RollAmount); + iirMidRollerAR = (iirMidRollerAR * (1.0 - RollAmount)) + (inputSampleR * RollAmount); + HighsSampleL = inputSampleL - iirMidRollerAL; + HighsSampleR = inputSampleR - iirMidRollerAR; + NonHighsSampleL = iirMidRollerAL; + NonHighsSampleR = iirMidRollerAR; + + iirHeadBumpAL += (inputSampleL * 0.05); + iirHeadBumpAR += (inputSampleR * 0.05); + iirHeadBumpAL -= (iirHeadBumpAL * iirHeadBumpAL * iirHeadBumpAL * HeadBumpFreq); + iirHeadBumpAR -= (iirHeadBumpAR * iirHeadBumpAR * iirHeadBumpAR * HeadBumpFreq); + iirHeadBumpAL = sin(iirHeadBumpAL); + iirHeadBumpAR = sin(iirHeadBumpAR); + + tempSample = (iirHeadBumpAL * biquadAL[2]) + biquadAL[7]; + biquadAL[7] = (iirHeadBumpAL * biquadAL[3]) - (tempSample * biquadAL[5]) + biquadAL[8]; + biquadAL[8] = (iirHeadBumpAL * biquadAL[4]) - (tempSample * biquadAL[6]); + iirHeadBumpAL = tempSample; //interleaved biquad + if (iirHeadBumpAL > 1.0) iirHeadBumpAL = 1.0; + if (iirHeadBumpAL < -1.0) iirHeadBumpAL = -1.0; + iirHeadBumpAL = asin(iirHeadBumpAL); + + tempSample = (iirHeadBumpAR * biquadAR[2]) + biquadAR[7]; + biquadAR[7] = (iirHeadBumpAR * biquadAR[3]) - (tempSample * biquadAR[5]) + biquadAR[8]; + biquadAR[8] = (iirHeadBumpAR * biquadAR[4]) - (tempSample * biquadAR[6]); + iirHeadBumpAR = tempSample; //interleaved biquad + if (iirHeadBumpAR > 1.0) iirHeadBumpAR = 1.0; + if (iirHeadBumpAR < -1.0) iirHeadBumpAR = -1.0; + iirHeadBumpAR = asin(iirHeadBumpAR); + + inputSampleL = sin(inputSampleL); + tempSample = (inputSampleL * biquadCL[2]) + biquadCL[7]; + biquadCL[7] = (inputSampleL * biquadCL[3]) - (tempSample * biquadCL[5]) + biquadCL[8]; + biquadCL[8] = (inputSampleL * biquadCL[4]) - (tempSample * biquadCL[6]); + inputSampleL = tempSample; //interleaved biquad + if (inputSampleL > 1.0) inputSampleL = 1.0; + if (inputSampleL < -1.0) inputSampleL = -1.0; + inputSampleL = asin(inputSampleL); + + inputSampleR = sin(inputSampleR); + tempSample = (inputSampleR * biquadCR[2]) + biquadCR[7]; + biquadCR[7] = (inputSampleR * biquadCR[3]) - (tempSample * biquadCR[5]) + biquadCR[8]; + biquadCR[8] = (inputSampleR * biquadCR[4]) - (tempSample * biquadCR[6]); + inputSampleR = tempSample; //interleaved biquad + if (inputSampleR > 1.0) inputSampleR = 1.0; + if (inputSampleR < -1.0) inputSampleR = -1.0; + inputSampleR = asin(inputSampleR); + } else { + iirMidRollerBL = (iirMidRollerBL * (1.0 - RollAmount)) + (inputSampleL * RollAmount); + iirMidRollerBR = (iirMidRollerBR * (1.0 - RollAmount)) + (inputSampleR * RollAmount); + HighsSampleL = inputSampleL - iirMidRollerBL; + HighsSampleR = inputSampleR - iirMidRollerBR; + NonHighsSampleL = iirMidRollerBL; + NonHighsSampleR = iirMidRollerBR; + + iirHeadBumpBL += (inputSampleL * 0.05); + iirHeadBumpBR += (inputSampleR * 0.05); + iirHeadBumpBL -= (iirHeadBumpBL * iirHeadBumpBL * iirHeadBumpBL * HeadBumpFreq); + iirHeadBumpBR -= (iirHeadBumpBR * iirHeadBumpBR * iirHeadBumpBR * HeadBumpFreq); + iirHeadBumpBL = sin(iirHeadBumpBL); + iirHeadBumpBR = sin(iirHeadBumpBR); + + tempSample = (iirHeadBumpBL * biquadBL[2]) + biquadBL[7]; + biquadBL[7] = (iirHeadBumpBL * biquadBL[3]) - (tempSample * biquadBL[5]) + biquadBL[8]; + biquadBL[8] = (iirHeadBumpBL * biquadBL[4]) - (tempSample * biquadBL[6]); + iirHeadBumpBL = tempSample; //interleaved biquad + if (iirHeadBumpBL > 1.0) iirHeadBumpBL = 1.0; + if (iirHeadBumpBL < -1.0) iirHeadBumpBL = -1.0; + iirHeadBumpBL = asin(iirHeadBumpBL); + + tempSample = (iirHeadBumpBR * biquadBR[2]) + biquadBR[7]; + biquadBR[7] = (iirHeadBumpBR * biquadBR[3]) - (tempSample * biquadBR[5]) + biquadBR[8]; + biquadBR[8] = (iirHeadBumpBR * biquadBR[4]) - (tempSample * biquadBR[6]); + iirHeadBumpBR = tempSample; //interleaved biquad + if (iirHeadBumpBR > 1.0) iirHeadBumpBR = 1.0; + if (iirHeadBumpBR < -1.0) iirHeadBumpBR = -1.0; + iirHeadBumpBR = asin(iirHeadBumpBR); + + inputSampleL = sin(inputSampleL); + tempSample = (inputSampleL * biquadDL[2]) + biquadDL[7]; + biquadDL[7] = (inputSampleL * biquadDL[3]) - (tempSample * biquadDL[5]) + biquadDL[8]; + biquadDL[8] = (inputSampleL * biquadDL[4]) - (tempSample * biquadDL[6]); + inputSampleL = tempSample; //interleaved biquad + if (inputSampleL > 1.0) inputSampleL = 1.0; + if (inputSampleL < -1.0) inputSampleL = -1.0; + inputSampleL = asin(inputSampleL); + + inputSampleR = sin(inputSampleR); + tempSample = (inputSampleR * biquadDR[2]) + biquadDR[7]; + biquadDR[7] = (inputSampleR * biquadDR[3]) - (tempSample * biquadDR[5]) + biquadDR[8]; + biquadDR[8] = (inputSampleR * biquadDR[4]) - (tempSample * biquadDR[6]); + inputSampleR = tempSample; //interleaved biquad + if (inputSampleR > 1.0) inputSampleR = 1.0; + if (inputSampleR < -1.0) inputSampleR = -1.0; + inputSampleR = asin(inputSampleR); + } + flip = !flip; + + long double groundSampleL = drySampleL - inputSampleL; //set up UnBox + long double groundSampleR = drySampleR - inputSampleR; //set up UnBox + + if (inputgain != 1.0) { + inputSampleL *= inputgain; + inputSampleR *= inputgain; + } //gain boost inside UnBox: do not boost fringe audio + + long double applySoften = fabs(HighsSampleL)*1.57079633; + if (applySoften > 1.57079633) applySoften = 1.57079633; + applySoften = 1-cos(applySoften); + if (HighsSampleL > 0) inputSampleL -= applySoften; + if (HighsSampleL < 0) inputSampleL += applySoften; + //apply Soften depending on polarity + applySoften = fabs(HighsSampleR)*1.57079633; + if (applySoften > 1.57079633) applySoften = 1.57079633; + applySoften = 1-cos(applySoften); + if (HighsSampleR > 0) inputSampleR -= applySoften; + if (HighsSampleR < 0) inputSampleR += applySoften; + //apply Soften depending on polarity + + if (inputSampleL > 1.2533141373155) inputSampleL = 1.2533141373155; + if (inputSampleL < -1.2533141373155) inputSampleL = -1.2533141373155; + //clip to 1.2533141373155 to reach maximum output + inputSampleL = sin(inputSampleL * fabs(inputSampleL)) / ((inputSampleL == 0.0) ?1:fabs(inputSampleL)); + //Spiral, for cleanest most optimal tape effect + if (inputSampleR > 1.2533141373155) inputSampleR = 1.2533141373155; + if (inputSampleR < -1.2533141373155) inputSampleR = -1.2533141373155; + //clip to 1.2533141373155 to reach maximum output + inputSampleR = sin(inputSampleR * fabs(inputSampleR)) / ((inputSampleR == 0.0) ?1:fabs(inputSampleR)); + //Spiral, for cleanest most optimal tape effect + + double suppress = (1.0-fabs(inputSampleL)) * 0.00013; + if (iirHeadBumpAL > suppress) iirHeadBumpAL -= suppress; + if (iirHeadBumpAL < -suppress) iirHeadBumpAL += suppress; + if (iirHeadBumpBL > suppress) iirHeadBumpBL -= suppress; + if (iirHeadBumpBL < -suppress) iirHeadBumpBL += suppress; + //restrain resonant quality of head bump algorithm + suppress = (1.0-fabs(inputSampleR)) * 0.00013; + if (iirHeadBumpAR > suppress) iirHeadBumpAR -= suppress; + if (iirHeadBumpAR < -suppress) iirHeadBumpAR += suppress; + if (iirHeadBumpBR > suppress) iirHeadBumpBR -= suppress; + if (iirHeadBumpBR < -suppress) iirHeadBumpBR += suppress; + //restrain resonant quality of head bump algorithm + + inputSampleL += groundSampleL; //apply UnBox processing + inputSampleR += groundSampleR; //apply UnBox processing + + inputSampleL += ((iirHeadBumpAL + iirHeadBumpBL) * 0.1);//and head bump + inputSampleR += ((iirHeadBumpAR + iirHeadBumpBR) * 0.1);//and head bump + + if (lastSampleL >= 0.99) + { + if (inputSampleL < 0.99) lastSampleL = ((0.99*softness) + (inputSampleL * (1.0-softness))); + else lastSampleL = 0.99; + } + + if (lastSampleL <= -0.99) + { + if (inputSampleL > -0.99) lastSampleL = ((-0.99*softness) + (inputSampleL * (1.0-softness))); + else lastSampleL = -0.99; + } + + if (inputSampleL > 0.99) + { + if (lastSampleL < 0.99) inputSampleL = ((0.99*softness) + (lastSampleL * (1.0-softness))); + else inputSampleL = 0.99; + } + + if (inputSampleL < -0.99) + { + if (lastSampleL > -0.99) inputSampleL = ((-0.99*softness) + (lastSampleL * (1.0-softness))); + else inputSampleL = -0.99; + } + lastSampleL = inputSampleL; //end ADClip L + + + if (lastSampleR >= 0.99) + { + if (inputSampleR < 0.99) lastSampleR = ((0.99*softness) + (inputSampleR * (1.0-softness))); + else lastSampleR = 0.99; + } + + if (lastSampleR <= -0.99) + { + if (inputSampleR > -0.99) lastSampleR = ((-0.99*softness) + (inputSampleR * (1.0-softness))); + else lastSampleR = -0.99; + } + + if (inputSampleR > 0.99) + { + if (lastSampleR < 0.99) inputSampleR = ((0.99*softness) + (lastSampleR * (1.0-softness))); + else inputSampleR = 0.99; + } + + if (inputSampleR < -0.99) + { + if (lastSampleR > -0.99) inputSampleR = ((-0.99*softness) + (lastSampleR * (1.0-softness))); + else inputSampleR = -0.99; + } + lastSampleR = inputSampleR; //end ADClip R + + if (inputSampleL > 0.99) inputSampleL = 0.99; + if (inputSampleL < -0.99) inputSampleL = -0.99; + //final iron bar + if (inputSampleR > 0.99) inputSampleR = 0.99; + if (inputSampleR < -0.99) inputSampleR = -0.99; + //final iron bar + + //begin 64 bit stereo floating point dither + int expon; frexp((double)inputSampleL, &expon); + fpd ^= fpd << 13; fpd ^= fpd >> 17; fpd ^= fpd << 5; + inputSampleL += ((double(fpd)-uint32_t(0x7fffffff)) * 1.1e-44l * pow(2,expon+62)); + frexp((double)inputSampleR, &expon); + fpd ^= fpd << 13; fpd ^= fpd >> 17; fpd ^= fpd << 5; + inputSampleR += ((double(fpd)-uint32_t(0x7fffffff)) * 1.1e-44l * pow(2,expon+62)); + //end 64 bit stereo floating point dither + + *out1 = inputSampleL; + *out2 = inputSampleR; + + *in1++; + *in2++; + *out1++; + *out2++; + } +} -- cgit v1.2.3