diff options
Diffstat (limited to 'plugins/MacVST/GrooveWear/source/GrooveWearProc.cpp')
-rwxr-xr-x | plugins/MacVST/GrooveWear/source/GrooveWearProc.cpp | 920 |
1 files changed, 920 insertions, 0 deletions
diff --git a/plugins/MacVST/GrooveWear/source/GrooveWearProc.cpp b/plugins/MacVST/GrooveWear/source/GrooveWearProc.cpp new file mode 100755 index 0000000..a22f6df --- /dev/null +++ b/plugins/MacVST/GrooveWear/source/GrooveWearProc.cpp @@ -0,0 +1,920 @@ +/* ======================================== + * GrooveWear - GrooveWear.h + * Copyright (c) 2016 airwindows, All rights reserved + * ======================================== */ + +#ifndef __GrooveWear_H +#include "GrooveWear.h" +#endif + +void GrooveWear::processReplacing(float **inputs, float **outputs, VstInt32 sampleFrames) +{ + float* in1 = inputs[0]; + float* in2 = inputs[1]; + float* out1 = outputs[0]; + float* out2 = outputs[1]; + + float fpTemp; + long double fpOld = 0.618033988749894848204586; //golden ratio! + long double fpNew = 1.0 - fpOld; + + double overallscale = (pow(A,2)*19.0)+1.0; + double gain = overallscale; + //mid groove wear + if (gain > 1.0) {fMid[0] = 1.0; gain -= 1.0;} else {fMid[0] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[1] = 1.0; gain -= 1.0;} else {fMid[1] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[2] = 1.0; gain -= 1.0;} else {fMid[2] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[3] = 1.0; gain -= 1.0;} else {fMid[3] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[4] = 1.0; gain -= 1.0;} else {fMid[4] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[5] = 1.0; gain -= 1.0;} else {fMid[5] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[6] = 1.0; gain -= 1.0;} else {fMid[6] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[7] = 1.0; gain -= 1.0;} else {fMid[7] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[8] = 1.0; gain -= 1.0;} else {fMid[8] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[9] = 1.0; gain -= 1.0;} else {fMid[9] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[10] = 1.0; gain -= 1.0;} else {fMid[10] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[11] = 1.0; gain -= 1.0;} else {fMid[11] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[12] = 1.0; gain -= 1.0;} else {fMid[12] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[13] = 1.0; gain -= 1.0;} else {fMid[13] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[14] = 1.0; gain -= 1.0;} else {fMid[14] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[15] = 1.0; gain -= 1.0;} else {fMid[15] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[16] = 1.0; gain -= 1.0;} else {fMid[16] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[17] = 1.0; gain -= 1.0;} else {fMid[17] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[18] = 1.0; gain -= 1.0;} else {fMid[18] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[19] = 1.0; gain -= 1.0;} else {fMid[19] = gain; gain = 0.0;} + //there, now we have a neat little moving average with remainders, in stereo + + if (overallscale < 1.0) overallscale = 1.0; + fMid[0] /= overallscale; + fMid[1] /= overallscale; + fMid[2] /= overallscale; + fMid[3] /= overallscale; + fMid[4] /= overallscale; + fMid[5] /= overallscale; + fMid[6] /= overallscale; + fMid[7] /= overallscale; + fMid[8] /= overallscale; + fMid[9] /= overallscale; + fMid[10] /= overallscale; + fMid[11] /= overallscale; + fMid[12] /= overallscale; + fMid[13] /= overallscale; + fMid[14] /= overallscale; + fMid[15] /= overallscale; + fMid[16] /= overallscale; + fMid[17] /= overallscale; + fMid[18] /= overallscale; + fMid[19] /= overallscale; + //and now it's neatly scaled, too + + double accumulatorSampleL; + double correctionL; + double accumulatorSampleR; + double correctionR; + + + double aWet = 1.0; + double bWet = 1.0; + double cWet = 1.0; + double dWet = B*4.0; + //four-stage wet/dry control using progressive stages that bypass when not engaged + if (dWet < 1.0) {aWet = dWet; bWet = 0.0; cWet = 0.0; dWet = 0.0;} + else if (dWet < 2.0) {bWet = dWet - 1.0; cWet = 0.0; dWet = 0.0;} + else if (dWet < 3.0) {cWet = dWet - 2.0; dWet = 0.0;} + else {dWet -= 3.0;} + //this is one way to make a little set of dry/wet stages that are successively added to the + //output as the control is turned up. Each one independently goes from 0-1 and stays at 1 + //beyond that point: this is a way to progressively add a 'black box' sound processing + //which lets you fall through to simpler processing at lower settings. + + //now we set them up so each full intensity one is blended evenly with dry for each stage. + //That's because the GrooveWear algorithm works best combined with dry. + aWet *= 0.5; + bWet *= 0.5; + cWet *= 0.5; + dWet *= 0.5; + //if you are using a more typical algorithm (like a sin() or something) you won't use this part + + double aDry = 1.0 - aWet; + double bDry = 1.0 - bWet; + double cDry = 1.0 - cWet; + double dDry = 1.0 - dWet; + + double drySampleL; + double drySampleR; + long double inputSampleL; + long double inputSampleR; + + 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 (aWet > 0.0) { + aMidL[19] = aMidL[18]; aMidL[18] = aMidL[17]; aMidL[17] = aMidL[16]; aMidL[16] = aMidL[15]; + aMidL[15] = aMidL[14]; aMidL[14] = aMidL[13]; aMidL[13] = aMidL[12]; aMidL[12] = aMidL[11]; + aMidL[11] = aMidL[10]; aMidL[10] = aMidL[9]; + aMidL[9] = aMidL[8]; aMidL[8] = aMidL[7]; aMidL[7] = aMidL[6]; aMidL[6] = aMidL[5]; + aMidL[5] = aMidL[4]; aMidL[4] = aMidL[3]; aMidL[3] = aMidL[2]; aMidL[2] = aMidL[1]; + aMidL[1] = aMidL[0]; aMidL[0] = accumulatorSampleL = (inputSampleL-aMidPrevL); + //this is different from Aura because that is accumulating rates of change OF the rate of change + //this is just averaging slews directly, and we have two stages of it. + + accumulatorSampleL *= fMid[0]; + accumulatorSampleL += (aMidL[1] * fMid[1]); + accumulatorSampleL += (aMidL[2] * fMid[2]); + accumulatorSampleL += (aMidL[3] * fMid[3]); + accumulatorSampleL += (aMidL[4] * fMid[4]); + accumulatorSampleL += (aMidL[5] * fMid[5]); + accumulatorSampleL += (aMidL[6] * fMid[6]); + accumulatorSampleL += (aMidL[7] * fMid[7]); + accumulatorSampleL += (aMidL[8] * fMid[8]); + accumulatorSampleL += (aMidL[9] * fMid[9]); + accumulatorSampleL += (aMidL[10] * fMid[10]); + accumulatorSampleL += (aMidL[11] * fMid[11]); + accumulatorSampleL += (aMidL[12] * fMid[12]); + accumulatorSampleL += (aMidL[13] * fMid[13]); + accumulatorSampleL += (aMidL[14] * fMid[14]); + accumulatorSampleL += (aMidL[15] * fMid[15]); + accumulatorSampleL += (aMidL[16] * fMid[16]); + accumulatorSampleL += (aMidL[17] * fMid[17]); + accumulatorSampleL += (aMidL[18] * fMid[18]); + accumulatorSampleL += (aMidL[19] * fMid[19]); + //we are doing our repetitive calculations on a separate value + correctionL = (inputSampleL-aMidPrevL) - accumulatorSampleL; + aMidPrevL = inputSampleL; + inputSampleL -= correctionL; + inputSampleL = (inputSampleL * aWet) + (drySampleL * aDry); + drySampleL = inputSampleL; + + aMidR[19] = aMidR[18]; aMidR[18] = aMidR[17]; aMidR[17] = aMidR[16]; aMidR[16] = aMidR[15]; + aMidR[15] = aMidR[14]; aMidR[14] = aMidR[13]; aMidR[13] = aMidR[12]; aMidR[12] = aMidR[11]; + aMidR[11] = aMidR[10]; aMidR[10] = aMidR[9]; + aMidR[9] = aMidR[8]; aMidR[8] = aMidR[7]; aMidR[7] = aMidR[6]; aMidR[6] = aMidR[5]; + aMidR[5] = aMidR[4]; aMidR[4] = aMidR[3]; aMidR[3] = aMidR[2]; aMidR[2] = aMidR[1]; + aMidR[1] = aMidR[0]; aMidR[0] = accumulatorSampleR = (inputSampleR-aMidPrevR); + //this is different from Aura because that is accumulating rates of change OF the rate of change + //this is just averaging slews directly, and we have two stages of it. + + accumulatorSampleR *= fMid[0]; + accumulatorSampleR += (aMidR[1] * fMid[1]); + accumulatorSampleR += (aMidR[2] * fMid[2]); + accumulatorSampleR += (aMidR[3] * fMid[3]); + accumulatorSampleR += (aMidR[4] * fMid[4]); + accumulatorSampleR += (aMidR[5] * fMid[5]); + accumulatorSampleR += (aMidR[6] * fMid[6]); + accumulatorSampleR += (aMidR[7] * fMid[7]); + accumulatorSampleR += (aMidR[8] * fMid[8]); + accumulatorSampleR += (aMidR[9] * fMid[9]); + accumulatorSampleR += (aMidR[10] * fMid[10]); + accumulatorSampleR += (aMidR[11] * fMid[11]); + accumulatorSampleR += (aMidR[12] * fMid[12]); + accumulatorSampleR += (aMidR[13] * fMid[13]); + accumulatorSampleR += (aMidR[14] * fMid[14]); + accumulatorSampleR += (aMidR[15] * fMid[15]); + accumulatorSampleR += (aMidR[16] * fMid[16]); + accumulatorSampleR += (aMidR[17] * fMid[17]); + accumulatorSampleR += (aMidR[18] * fMid[18]); + accumulatorSampleR += (aMidR[19] * fMid[19]); + //we are doing our repetitive calculations on a separate value + correctionR = (inputSampleR-aMidPrevR) - accumulatorSampleR; + aMidPrevR = inputSampleR; + inputSampleR -= correctionR; + inputSampleR = (inputSampleR * aWet) + (drySampleR * aDry); + drySampleR = inputSampleR; + } + + if (bWet > 0.0) { + bMidL[19] = bMidL[18]; bMidL[18] = bMidL[17]; bMidL[17] = bMidL[16]; bMidL[16] = bMidL[15]; + bMidL[15] = bMidL[14]; bMidL[14] = bMidL[13]; bMidL[13] = bMidL[12]; bMidL[12] = bMidL[11]; + bMidL[11] = bMidL[10]; bMidL[10] = bMidL[9]; + bMidL[9] = bMidL[8]; bMidL[8] = bMidL[7]; bMidL[7] = bMidL[6]; bMidL[6] = bMidL[5]; + bMidL[5] = bMidL[4]; bMidL[4] = bMidL[3]; bMidL[3] = bMidL[2]; bMidL[2] = bMidL[1]; + bMidL[1] = bMidL[0]; bMidL[0] = accumulatorSampleL = (inputSampleL-bMidPrevL); + + accumulatorSampleL *= fMid[0]; + accumulatorSampleL += (bMidL[1] * fMid[1]); + accumulatorSampleL += (bMidL[2] * fMid[2]); + accumulatorSampleL += (bMidL[3] * fMid[3]); + accumulatorSampleL += (bMidL[4] * fMid[4]); + accumulatorSampleL += (bMidL[5] * fMid[5]); + accumulatorSampleL += (bMidL[6] * fMid[6]); + accumulatorSampleL += (bMidL[7] * fMid[7]); + accumulatorSampleL += (bMidL[8] * fMid[8]); + accumulatorSampleL += (bMidL[9] * fMid[9]); + accumulatorSampleL += (bMidL[10] * fMid[10]); + accumulatorSampleL += (bMidL[11] * fMid[11]); + accumulatorSampleL += (bMidL[12] * fMid[12]); + accumulatorSampleL += (bMidL[13] * fMid[13]); + accumulatorSampleL += (bMidL[14] * fMid[14]); + accumulatorSampleL += (bMidL[15] * fMid[15]); + accumulatorSampleL += (bMidL[16] * fMid[16]); + accumulatorSampleL += (bMidL[17] * fMid[17]); + accumulatorSampleL += (bMidL[18] * fMid[18]); + accumulatorSampleL += (bMidL[19] * fMid[19]); + //we are doing our repetitive calculations on a separate value + correctionL = (inputSampleL-bMidPrevL) - accumulatorSampleL; + bMidPrevL = inputSampleL; + inputSampleL -= correctionL; + inputSampleL = (inputSampleL * bWet) + (drySampleL * bDry); + drySampleL = inputSampleL; + + bMidR[19] = bMidR[18]; bMidR[18] = bMidR[17]; bMidR[17] = bMidR[16]; bMidR[16] = bMidR[15]; + bMidR[15] = bMidR[14]; bMidR[14] = bMidR[13]; bMidR[13] = bMidR[12]; bMidR[12] = bMidR[11]; + bMidR[11] = bMidR[10]; bMidR[10] = bMidR[9]; + bMidR[9] = bMidR[8]; bMidR[8] = bMidR[7]; bMidR[7] = bMidR[6]; bMidR[6] = bMidR[5]; + bMidR[5] = bMidR[4]; bMidR[4] = bMidR[3]; bMidR[3] = bMidR[2]; bMidR[2] = bMidR[1]; + bMidR[1] = bMidR[0]; bMidR[0] = accumulatorSampleR = (inputSampleR-bMidPrevR); + + accumulatorSampleR *= fMid[0]; + accumulatorSampleR += (bMidR[1] * fMid[1]); + accumulatorSampleR += (bMidR[2] * fMid[2]); + accumulatorSampleR += (bMidR[3] * fMid[3]); + accumulatorSampleR += (bMidR[4] * fMid[4]); + accumulatorSampleR += (bMidR[5] * fMid[5]); + accumulatorSampleR += (bMidR[6] * fMid[6]); + accumulatorSampleR += (bMidR[7] * fMid[7]); + accumulatorSampleR += (bMidR[8] * fMid[8]); + accumulatorSampleR += (bMidR[9] * fMid[9]); + accumulatorSampleR += (bMidR[10] * fMid[10]); + accumulatorSampleR += (bMidR[11] * fMid[11]); + accumulatorSampleR += (bMidR[12] * fMid[12]); + accumulatorSampleR += (bMidR[13] * fMid[13]); + accumulatorSampleR += (bMidR[14] * fMid[14]); + accumulatorSampleR += (bMidR[15] * fMid[15]); + accumulatorSampleR += (bMidR[16] * fMid[16]); + accumulatorSampleR += (bMidR[17] * fMid[17]); + accumulatorSampleR += (bMidR[18] * fMid[18]); + accumulatorSampleR += (bMidR[19] * fMid[19]); + //we are doing our repetitive calculations on a separate value + correctionR = (inputSampleR-bMidPrevR) - accumulatorSampleR; + bMidPrevR = inputSampleR; + inputSampleR -= correctionR; + inputSampleR = (inputSampleR * bWet) + (drySampleR * bDry); + drySampleR = inputSampleR; + } + + if (cWet > 0.0) { + cMidL[19] = cMidL[18]; cMidL[18] = cMidL[17]; cMidL[17] = cMidL[16]; cMidL[16] = cMidL[15]; + cMidL[15] = cMidL[14]; cMidL[14] = cMidL[13]; cMidL[13] = cMidL[12]; cMidL[12] = cMidL[11]; + cMidL[11] = cMidL[10]; cMidL[10] = cMidL[9]; + cMidL[9] = cMidL[8]; cMidL[8] = cMidL[7]; cMidL[7] = cMidL[6]; cMidL[6] = cMidL[5]; + cMidL[5] = cMidL[4]; cMidL[4] = cMidL[3]; cMidL[3] = cMidL[2]; cMidL[2] = cMidL[1]; + cMidL[1] = cMidL[0]; cMidL[0] = accumulatorSampleL = (inputSampleL-cMidPrevL); + + accumulatorSampleL *= fMid[0]; + accumulatorSampleL += (cMidL[1] * fMid[1]); + accumulatorSampleL += (cMidL[2] * fMid[2]); + accumulatorSampleL += (cMidL[3] * fMid[3]); + accumulatorSampleL += (cMidL[4] * fMid[4]); + accumulatorSampleL += (cMidL[5] * fMid[5]); + accumulatorSampleL += (cMidL[6] * fMid[6]); + accumulatorSampleL += (cMidL[7] * fMid[7]); + accumulatorSampleL += (cMidL[8] * fMid[8]); + accumulatorSampleL += (cMidL[9] * fMid[9]); + accumulatorSampleL += (cMidL[10] * fMid[10]); + accumulatorSampleL += (cMidL[11] * fMid[11]); + accumulatorSampleL += (cMidL[12] * fMid[12]); + accumulatorSampleL += (cMidL[13] * fMid[13]); + accumulatorSampleL += (cMidL[14] * fMid[14]); + accumulatorSampleL += (cMidL[15] * fMid[15]); + accumulatorSampleL += (cMidL[16] * fMid[16]); + accumulatorSampleL += (cMidL[17] * fMid[17]); + accumulatorSampleL += (cMidL[18] * fMid[18]); + accumulatorSampleL += (cMidL[19] * fMid[19]); + //we are doing our repetitive calculations on a separate value + correctionL = (inputSampleL-cMidPrevL) - accumulatorSampleL; + cMidPrevL = inputSampleL; + inputSampleL -= correctionL; + inputSampleL = (inputSampleL * cWet) + (drySampleL * cDry); + drySampleL = inputSampleL; + + cMidR[19] = cMidR[18]; cMidR[18] = cMidR[17]; cMidR[17] = cMidR[16]; cMidR[16] = cMidR[15]; + cMidR[15] = cMidR[14]; cMidR[14] = cMidR[13]; cMidR[13] = cMidR[12]; cMidR[12] = cMidR[11]; + cMidR[11] = cMidR[10]; cMidR[10] = cMidR[9]; + cMidR[9] = cMidR[8]; cMidR[8] = cMidR[7]; cMidR[7] = cMidR[6]; cMidR[6] = cMidR[5]; + cMidR[5] = cMidR[4]; cMidR[4] = cMidR[3]; cMidR[3] = cMidR[2]; cMidR[2] = cMidR[1]; + cMidR[1] = cMidR[0]; cMidR[0] = accumulatorSampleR = (inputSampleR-cMidPrevR); + + accumulatorSampleR *= fMid[0]; + accumulatorSampleR += (cMidR[1] * fMid[1]); + accumulatorSampleR += (cMidR[2] * fMid[2]); + accumulatorSampleR += (cMidR[3] * fMid[3]); + accumulatorSampleR += (cMidR[4] * fMid[4]); + accumulatorSampleR += (cMidR[5] * fMid[5]); + accumulatorSampleR += (cMidR[6] * fMid[6]); + accumulatorSampleR += (cMidR[7] * fMid[7]); + accumulatorSampleR += (cMidR[8] * fMid[8]); + accumulatorSampleR += (cMidR[9] * fMid[9]); + accumulatorSampleR += (cMidR[10] * fMid[10]); + accumulatorSampleR += (cMidR[11] * fMid[11]); + accumulatorSampleR += (cMidR[12] * fMid[12]); + accumulatorSampleR += (cMidR[13] * fMid[13]); + accumulatorSampleR += (cMidR[14] * fMid[14]); + accumulatorSampleR += (cMidR[15] * fMid[15]); + accumulatorSampleR += (cMidR[16] * fMid[16]); + accumulatorSampleR += (cMidR[17] * fMid[17]); + accumulatorSampleR += (cMidR[18] * fMid[18]); + accumulatorSampleR += (cMidR[19] * fMid[19]); + //we are doing our repetitive calculations on a separate value + correctionR = (inputSampleR-cMidPrevR) - accumulatorSampleR; + cMidPrevR = inputSampleR; + inputSampleR -= correctionR; + inputSampleR = (inputSampleR * cWet) + (drySampleR * cDry); + drySampleR = inputSampleR; + } + + if (dWet > 0.0) { + dMidL[19] = dMidL[18]; dMidL[18] = dMidL[17]; dMidL[17] = dMidL[16]; dMidL[16] = dMidL[15]; + dMidL[15] = dMidL[14]; dMidL[14] = dMidL[13]; dMidL[13] = dMidL[12]; dMidL[12] = dMidL[11]; + dMidL[11] = dMidL[10]; dMidL[10] = dMidL[9]; + dMidL[9] = dMidL[8]; dMidL[8] = dMidL[7]; dMidL[7] = dMidL[6]; dMidL[6] = dMidL[5]; + dMidL[5] = dMidL[4]; dMidL[4] = dMidL[3]; dMidL[3] = dMidL[2]; dMidL[2] = dMidL[1]; + dMidL[1] = dMidL[0]; dMidL[0] = accumulatorSampleL = (inputSampleL-dMidPrevL); + + accumulatorSampleL *= fMid[0]; + accumulatorSampleL += (dMidL[1] * fMid[1]); + accumulatorSampleL += (dMidL[2] * fMid[2]); + accumulatorSampleL += (dMidL[3] * fMid[3]); + accumulatorSampleL += (dMidL[4] * fMid[4]); + accumulatorSampleL += (dMidL[5] * fMid[5]); + accumulatorSampleL += (dMidL[6] * fMid[6]); + accumulatorSampleL += (dMidL[7] * fMid[7]); + accumulatorSampleL += (dMidL[8] * fMid[8]); + accumulatorSampleL += (dMidL[9] * fMid[9]); + accumulatorSampleL += (dMidL[10] * fMid[10]); + accumulatorSampleL += (dMidL[11] * fMid[11]); + accumulatorSampleL += (dMidL[12] * fMid[12]); + accumulatorSampleL += (dMidL[13] * fMid[13]); + accumulatorSampleL += (dMidL[14] * fMid[14]); + accumulatorSampleL += (dMidL[15] * fMid[15]); + accumulatorSampleL += (dMidL[16] * fMid[16]); + accumulatorSampleL += (dMidL[17] * fMid[17]); + accumulatorSampleL += (dMidL[18] * fMid[18]); + accumulatorSampleL += (dMidL[19] * fMid[19]); + //we are doing our repetitive calculations on a separate value + correctionL = (inputSampleL-dMidPrevL) - accumulatorSampleL; + dMidPrevL = inputSampleL; + inputSampleL -= correctionL; + inputSampleL = (inputSampleL * dWet) + (drySampleL * dDry); + drySampleL = inputSampleL; + + dMidR[19] = dMidR[18]; dMidR[18] = dMidR[17]; dMidR[17] = dMidR[16]; dMidR[16] = dMidR[15]; + dMidR[15] = dMidR[14]; dMidR[14] = dMidR[13]; dMidR[13] = dMidR[12]; dMidR[12] = dMidR[11]; + dMidR[11] = dMidR[10]; dMidR[10] = dMidR[9]; + dMidR[9] = dMidR[8]; dMidR[8] = dMidR[7]; dMidR[7] = dMidR[6]; dMidR[6] = dMidR[5]; + dMidR[5] = dMidR[4]; dMidR[4] = dMidR[3]; dMidR[3] = dMidR[2]; dMidR[2] = dMidR[1]; + dMidR[1] = dMidR[0]; dMidR[0] = accumulatorSampleR = (inputSampleR-dMidPrevR); + + accumulatorSampleR *= fMid[0]; + accumulatorSampleR += (dMidR[1] * fMid[1]); + accumulatorSampleR += (dMidR[2] * fMid[2]); + accumulatorSampleR += (dMidR[3] * fMid[3]); + accumulatorSampleR += (dMidR[4] * fMid[4]); + accumulatorSampleR += (dMidR[5] * fMid[5]); + accumulatorSampleR += (dMidR[6] * fMid[6]); + accumulatorSampleR += (dMidR[7] * fMid[7]); + accumulatorSampleR += (dMidR[8] * fMid[8]); + accumulatorSampleR += (dMidR[9] * fMid[9]); + accumulatorSampleR += (dMidR[10] * fMid[10]); + accumulatorSampleR += (dMidR[11] * fMid[11]); + accumulatorSampleR += (dMidR[12] * fMid[12]); + accumulatorSampleR += (dMidR[13] * fMid[13]); + accumulatorSampleR += (dMidR[14] * fMid[14]); + accumulatorSampleR += (dMidR[15] * fMid[15]); + accumulatorSampleR += (dMidR[16] * fMid[16]); + accumulatorSampleR += (dMidR[17] * fMid[17]); + accumulatorSampleR += (dMidR[18] * fMid[18]); + accumulatorSampleR += (dMidR[19] * fMid[19]); + //we are doing our repetitive calculations on a separate value + correctionR = (inputSampleR-dMidPrevR) - accumulatorSampleR; + dMidPrevR = inputSampleR; + inputSampleR -= correctionR; + inputSampleR = (inputSampleR * dWet) + (drySampleR * dDry); + drySampleR = inputSampleR; + } + + //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 GrooveWear::processDoubleReplacing(double **inputs, double **outputs, VstInt32 sampleFrames) +{ + double* in1 = inputs[0]; + double* in2 = inputs[1]; + double* out1 = outputs[0]; + double* out2 = outputs[1]; + + double fpTemp; + long double fpOld = 0.618033988749894848204586; //golden ratio! + long double fpNew = 1.0 - fpOld; + + double overallscale = (pow(A,2)*19.0)+1.0; + double gain = overallscale; + //mid groove wear + if (gain > 1.0) {fMid[0] = 1.0; gain -= 1.0;} else {fMid[0] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[1] = 1.0; gain -= 1.0;} else {fMid[1] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[2] = 1.0; gain -= 1.0;} else {fMid[2] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[3] = 1.0; gain -= 1.0;} else {fMid[3] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[4] = 1.0; gain -= 1.0;} else {fMid[4] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[5] = 1.0; gain -= 1.0;} else {fMid[5] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[6] = 1.0; gain -= 1.0;} else {fMid[6] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[7] = 1.0; gain -= 1.0;} else {fMid[7] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[8] = 1.0; gain -= 1.0;} else {fMid[8] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[9] = 1.0; gain -= 1.0;} else {fMid[9] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[10] = 1.0; gain -= 1.0;} else {fMid[10] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[11] = 1.0; gain -= 1.0;} else {fMid[11] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[12] = 1.0; gain -= 1.0;} else {fMid[12] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[13] = 1.0; gain -= 1.0;} else {fMid[13] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[14] = 1.0; gain -= 1.0;} else {fMid[14] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[15] = 1.0; gain -= 1.0;} else {fMid[15] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[16] = 1.0; gain -= 1.0;} else {fMid[16] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[17] = 1.0; gain -= 1.0;} else {fMid[17] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[18] = 1.0; gain -= 1.0;} else {fMid[18] = gain; gain = 0.0;} + if (gain > 1.0) {fMid[19] = 1.0; gain -= 1.0;} else {fMid[19] = gain; gain = 0.0;} + //there, now we have a neat little moving average with remainders, in stereo + + if (overallscale < 1.0) overallscale = 1.0; + fMid[0] /= overallscale; + fMid[1] /= overallscale; + fMid[2] /= overallscale; + fMid[3] /= overallscale; + fMid[4] /= overallscale; + fMid[5] /= overallscale; + fMid[6] /= overallscale; + fMid[7] /= overallscale; + fMid[8] /= overallscale; + fMid[9] /= overallscale; + fMid[10] /= overallscale; + fMid[11] /= overallscale; + fMid[12] /= overallscale; + fMid[13] /= overallscale; + fMid[14] /= overallscale; + fMid[15] /= overallscale; + fMid[16] /= overallscale; + fMid[17] /= overallscale; + fMid[18] /= overallscale; + fMid[19] /= overallscale; + //and now it's neatly scaled, too + + double accumulatorSampleL; + double correctionL; + double accumulatorSampleR; + double correctionR; + + + double aWet = 1.0; + double bWet = 1.0; + double cWet = 1.0; + double dWet = B*4.0; + //four-stage wet/dry control using progressive stages that bypass when not engaged + if (dWet < 1.0) {aWet = dWet; bWet = 0.0; cWet = 0.0; dWet = 0.0;} + else if (dWet < 2.0) {bWet = dWet - 1.0; cWet = 0.0; dWet = 0.0;} + else if (dWet < 3.0) {cWet = dWet - 2.0; dWet = 0.0;} + else {dWet -= 3.0;} + //this is one way to make a little set of dry/wet stages that are successively added to the + //output as the control is turned up. Each one independently goes from 0-1 and stays at 1 + //beyond that point: this is a way to progressively add a 'black box' sound processing + //which lets you fall through to simpler processing at lower settings. + + //now we set them up so each full intensity one is blended evenly with dry for each stage. + //That's because the GrooveWear algorithm works best combined with dry. + aWet *= 0.5; + bWet *= 0.5; + cWet *= 0.5; + dWet *= 0.5; + //if you are using a more typical algorithm (like a sin() or something) you won't use this part + + double aDry = 1.0 - aWet; + double bDry = 1.0 - bWet; + double cDry = 1.0 - cWet; + double dDry = 1.0 - dWet; + + double drySampleL; + double drySampleR; + long double inputSampleL; + long double inputSampleR; + + 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 (aWet > 0.0) { + aMidL[19] = aMidL[18]; aMidL[18] = aMidL[17]; aMidL[17] = aMidL[16]; aMidL[16] = aMidL[15]; + aMidL[15] = aMidL[14]; aMidL[14] = aMidL[13]; aMidL[13] = aMidL[12]; aMidL[12] = aMidL[11]; + aMidL[11] = aMidL[10]; aMidL[10] = aMidL[9]; + aMidL[9] = aMidL[8]; aMidL[8] = aMidL[7]; aMidL[7] = aMidL[6]; aMidL[6] = aMidL[5]; + aMidL[5] = aMidL[4]; aMidL[4] = aMidL[3]; aMidL[3] = aMidL[2]; aMidL[2] = aMidL[1]; + aMidL[1] = aMidL[0]; aMidL[0] = accumulatorSampleL = (inputSampleL-aMidPrevL); + //this is different from Aura because that is accumulating rates of change OF the rate of change + //this is just averaging slews directly, and we have two stages of it. + + accumulatorSampleL *= fMid[0]; + accumulatorSampleL += (aMidL[1] * fMid[1]); + accumulatorSampleL += (aMidL[2] * fMid[2]); + accumulatorSampleL += (aMidL[3] * fMid[3]); + accumulatorSampleL += (aMidL[4] * fMid[4]); + accumulatorSampleL += (aMidL[5] * fMid[5]); + accumulatorSampleL += (aMidL[6] * fMid[6]); + accumulatorSampleL += (aMidL[7] * fMid[7]); + accumulatorSampleL += (aMidL[8] * fMid[8]); + accumulatorSampleL += (aMidL[9] * fMid[9]); + accumulatorSampleL += (aMidL[10] * fMid[10]); + accumulatorSampleL += (aMidL[11] * fMid[11]); + accumulatorSampleL += (aMidL[12] * fMid[12]); + accumulatorSampleL += (aMidL[13] * fMid[13]); + accumulatorSampleL += (aMidL[14] * fMid[14]); + accumulatorSampleL += (aMidL[15] * fMid[15]); + accumulatorSampleL += (aMidL[16] * fMid[16]); + accumulatorSampleL += (aMidL[17] * fMid[17]); + accumulatorSampleL += (aMidL[18] * fMid[18]); + accumulatorSampleL += (aMidL[19] * fMid[19]); + //we are doing our repetitive calculations on a separate value + correctionL = (inputSampleL-aMidPrevL) - accumulatorSampleL; + aMidPrevL = inputSampleL; + inputSampleL -= correctionL; + inputSampleL = (inputSampleL * aWet) + (drySampleL * aDry); + drySampleL = inputSampleL; + + aMidR[19] = aMidR[18]; aMidR[18] = aMidR[17]; aMidR[17] = aMidR[16]; aMidR[16] = aMidR[15]; + aMidR[15] = aMidR[14]; aMidR[14] = aMidR[13]; aMidR[13] = aMidR[12]; aMidR[12] = aMidR[11]; + aMidR[11] = aMidR[10]; aMidR[10] = aMidR[9]; + aMidR[9] = aMidR[8]; aMidR[8] = aMidR[7]; aMidR[7] = aMidR[6]; aMidR[6] = aMidR[5]; + aMidR[5] = aMidR[4]; aMidR[4] = aMidR[3]; aMidR[3] = aMidR[2]; aMidR[2] = aMidR[1]; + aMidR[1] = aMidR[0]; aMidR[0] = accumulatorSampleR = (inputSampleR-aMidPrevR); + //this is different from Aura because that is accumulating rates of change OF the rate of change + //this is just averaging slews directly, and we have two stages of it. + + accumulatorSampleR *= fMid[0]; + accumulatorSampleR += (aMidR[1] * fMid[1]); + accumulatorSampleR += (aMidR[2] * fMid[2]); + accumulatorSampleR += (aMidR[3] * fMid[3]); + accumulatorSampleR += (aMidR[4] * fMid[4]); + accumulatorSampleR += (aMidR[5] * fMid[5]); + accumulatorSampleR += (aMidR[6] * fMid[6]); + accumulatorSampleR += (aMidR[7] * fMid[7]); + accumulatorSampleR += (aMidR[8] * fMid[8]); + accumulatorSampleR += (aMidR[9] * fMid[9]); + accumulatorSampleR += (aMidR[10] * fMid[10]); + accumulatorSampleR += (aMidR[11] * fMid[11]); + accumulatorSampleR += (aMidR[12] * fMid[12]); + accumulatorSampleR += (aMidR[13] * fMid[13]); + accumulatorSampleR += (aMidR[14] * fMid[14]); + accumulatorSampleR += (aMidR[15] * fMid[15]); + accumulatorSampleR += (aMidR[16] * fMid[16]); + accumulatorSampleR += (aMidR[17] * fMid[17]); + accumulatorSampleR += (aMidR[18] * fMid[18]); + accumulatorSampleR += (aMidR[19] * fMid[19]); + //we are doing our repetitive calculations on a separate value + correctionR = (inputSampleR-aMidPrevR) - accumulatorSampleR; + aMidPrevR = inputSampleR; + inputSampleR -= correctionR; + inputSampleR = (inputSampleR * aWet) + (drySampleR * aDry); + drySampleR = inputSampleR; + } + + if (bWet > 0.0) { + bMidL[19] = bMidL[18]; bMidL[18] = bMidL[17]; bMidL[17] = bMidL[16]; bMidL[16] = bMidL[15]; + bMidL[15] = bMidL[14]; bMidL[14] = bMidL[13]; bMidL[13] = bMidL[12]; bMidL[12] = bMidL[11]; + bMidL[11] = bMidL[10]; bMidL[10] = bMidL[9]; + bMidL[9] = bMidL[8]; bMidL[8] = bMidL[7]; bMidL[7] = bMidL[6]; bMidL[6] = bMidL[5]; + bMidL[5] = bMidL[4]; bMidL[4] = bMidL[3]; bMidL[3] = bMidL[2]; bMidL[2] = bMidL[1]; + bMidL[1] = bMidL[0]; bMidL[0] = accumulatorSampleL = (inputSampleL-bMidPrevL); + + accumulatorSampleL *= fMid[0]; + accumulatorSampleL += (bMidL[1] * fMid[1]); + accumulatorSampleL += (bMidL[2] * fMid[2]); + accumulatorSampleL += (bMidL[3] * fMid[3]); + accumulatorSampleL += (bMidL[4] * fMid[4]); + accumulatorSampleL += (bMidL[5] * fMid[5]); + accumulatorSampleL += (bMidL[6] * fMid[6]); + accumulatorSampleL += (bMidL[7] * fMid[7]); + accumulatorSampleL += (bMidL[8] * fMid[8]); + accumulatorSampleL += (bMidL[9] * fMid[9]); + accumulatorSampleL += (bMidL[10] * fMid[10]); + accumulatorSampleL += (bMidL[11] * fMid[11]); + accumulatorSampleL += (bMidL[12] * fMid[12]); + accumulatorSampleL += (bMidL[13] * fMid[13]); + accumulatorSampleL += (bMidL[14] * fMid[14]); + accumulatorSampleL += (bMidL[15] * fMid[15]); + accumulatorSampleL += (bMidL[16] * fMid[16]); + accumulatorSampleL += (bMidL[17] * fMid[17]); + accumulatorSampleL += (bMidL[18] * fMid[18]); + accumulatorSampleL += (bMidL[19] * fMid[19]); + //we are doing our repetitive calculations on a separate value + correctionL = (inputSampleL-bMidPrevL) - accumulatorSampleL; + bMidPrevL = inputSampleL; + inputSampleL -= correctionL; + inputSampleL = (inputSampleL * bWet) + (drySampleL * bDry); + drySampleL = inputSampleL; + + bMidR[19] = bMidR[18]; bMidR[18] = bMidR[17]; bMidR[17] = bMidR[16]; bMidR[16] = bMidR[15]; + bMidR[15] = bMidR[14]; bMidR[14] = bMidR[13]; bMidR[13] = bMidR[12]; bMidR[12] = bMidR[11]; + bMidR[11] = bMidR[10]; bMidR[10] = bMidR[9]; + bMidR[9] = bMidR[8]; bMidR[8] = bMidR[7]; bMidR[7] = bMidR[6]; bMidR[6] = bMidR[5]; + bMidR[5] = bMidR[4]; bMidR[4] = bMidR[3]; bMidR[3] = bMidR[2]; bMidR[2] = bMidR[1]; + bMidR[1] = bMidR[0]; bMidR[0] = accumulatorSampleR = (inputSampleR-bMidPrevR); + + accumulatorSampleR *= fMid[0]; + accumulatorSampleR += (bMidR[1] * fMid[1]); + accumulatorSampleR += (bMidR[2] * fMid[2]); + accumulatorSampleR += (bMidR[3] * fMid[3]); + accumulatorSampleR += (bMidR[4] * fMid[4]); + accumulatorSampleR += (bMidR[5] * fMid[5]); + accumulatorSampleR += (bMidR[6] * fMid[6]); + accumulatorSampleR += (bMidR[7] * fMid[7]); + accumulatorSampleR += (bMidR[8] * fMid[8]); + accumulatorSampleR += (bMidR[9] * fMid[9]); + accumulatorSampleR += (bMidR[10] * fMid[10]); + accumulatorSampleR += (bMidR[11] * fMid[11]); + accumulatorSampleR += (bMidR[12] * fMid[12]); + accumulatorSampleR += (bMidR[13] * fMid[13]); + accumulatorSampleR += (bMidR[14] * fMid[14]); + accumulatorSampleR += (bMidR[15] * fMid[15]); + accumulatorSampleR += (bMidR[16] * fMid[16]); + accumulatorSampleR += (bMidR[17] * fMid[17]); + accumulatorSampleR += (bMidR[18] * fMid[18]); + accumulatorSampleR += (bMidR[19] * fMid[19]); + //we are doing our repetitive calculations on a separate value + correctionR = (inputSampleR-bMidPrevR) - accumulatorSampleR; + bMidPrevR = inputSampleR; + inputSampleR -= correctionR; + inputSampleR = (inputSampleR * bWet) + (drySampleR * bDry); + drySampleR = inputSampleR; + } + + if (cWet > 0.0) { + cMidL[19] = cMidL[18]; cMidL[18] = cMidL[17]; cMidL[17] = cMidL[16]; cMidL[16] = cMidL[15]; + cMidL[15] = cMidL[14]; cMidL[14] = cMidL[13]; cMidL[13] = cMidL[12]; cMidL[12] = cMidL[11]; + cMidL[11] = cMidL[10]; cMidL[10] = cMidL[9]; + cMidL[9] = cMidL[8]; cMidL[8] = cMidL[7]; cMidL[7] = cMidL[6]; cMidL[6] = cMidL[5]; + cMidL[5] = cMidL[4]; cMidL[4] = cMidL[3]; cMidL[3] = cMidL[2]; cMidL[2] = cMidL[1]; + cMidL[1] = cMidL[0]; cMidL[0] = accumulatorSampleL = (inputSampleL-cMidPrevL); + + accumulatorSampleL *= fMid[0]; + accumulatorSampleL += (cMidL[1] * fMid[1]); + accumulatorSampleL += (cMidL[2] * fMid[2]); + accumulatorSampleL += (cMidL[3] * fMid[3]); + accumulatorSampleL += (cMidL[4] * fMid[4]); + accumulatorSampleL += (cMidL[5] * fMid[5]); + accumulatorSampleL += (cMidL[6] * fMid[6]); + accumulatorSampleL += (cMidL[7] * fMid[7]); + accumulatorSampleL += (cMidL[8] * fMid[8]); + accumulatorSampleL += (cMidL[9] * fMid[9]); + accumulatorSampleL += (cMidL[10] * fMid[10]); + accumulatorSampleL += (cMidL[11] * fMid[11]); + accumulatorSampleL += (cMidL[12] * fMid[12]); + accumulatorSampleL += (cMidL[13] * fMid[13]); + accumulatorSampleL += (cMidL[14] * fMid[14]); + accumulatorSampleL += (cMidL[15] * fMid[15]); + accumulatorSampleL += (cMidL[16] * fMid[16]); + accumulatorSampleL += (cMidL[17] * fMid[17]); + accumulatorSampleL += (cMidL[18] * fMid[18]); + accumulatorSampleL += (cMidL[19] * fMid[19]); + //we are doing our repetitive calculations on a separate value + correctionL = (inputSampleL-cMidPrevL) - accumulatorSampleL; + cMidPrevL = inputSampleL; + inputSampleL -= correctionL; + inputSampleL = (inputSampleL * cWet) + (drySampleL * cDry); + drySampleL = inputSampleL; + + cMidR[19] = cMidR[18]; cMidR[18] = cMidR[17]; cMidR[17] = cMidR[16]; cMidR[16] = cMidR[15]; + cMidR[15] = cMidR[14]; cMidR[14] = cMidR[13]; cMidR[13] = cMidR[12]; cMidR[12] = cMidR[11]; + cMidR[11] = cMidR[10]; cMidR[10] = cMidR[9]; + cMidR[9] = cMidR[8]; cMidR[8] = cMidR[7]; cMidR[7] = cMidR[6]; cMidR[6] = cMidR[5]; + cMidR[5] = cMidR[4]; cMidR[4] = cMidR[3]; cMidR[3] = cMidR[2]; cMidR[2] = cMidR[1]; + cMidR[1] = cMidR[0]; cMidR[0] = accumulatorSampleR = (inputSampleR-cMidPrevR); + + accumulatorSampleR *= fMid[0]; + accumulatorSampleR += (cMidR[1] * fMid[1]); + accumulatorSampleR += (cMidR[2] * fMid[2]); + accumulatorSampleR += (cMidR[3] * fMid[3]); + accumulatorSampleR += (cMidR[4] * fMid[4]); + accumulatorSampleR += (cMidR[5] * fMid[5]); + accumulatorSampleR += (cMidR[6] * fMid[6]); + accumulatorSampleR += (cMidR[7] * fMid[7]); + accumulatorSampleR += (cMidR[8] * fMid[8]); + accumulatorSampleR += (cMidR[9] * fMid[9]); + accumulatorSampleR += (cMidR[10] * fMid[10]); + accumulatorSampleR += (cMidR[11] * fMid[11]); + accumulatorSampleR += (cMidR[12] * fMid[12]); + accumulatorSampleR += (cMidR[13] * fMid[13]); + accumulatorSampleR += (cMidR[14] * fMid[14]); + accumulatorSampleR += (cMidR[15] * fMid[15]); + accumulatorSampleR += (cMidR[16] * fMid[16]); + accumulatorSampleR += (cMidR[17] * fMid[17]); + accumulatorSampleR += (cMidR[18] * fMid[18]); + accumulatorSampleR += (cMidR[19] * fMid[19]); + //we are doing our repetitive calculations on a separate value + correctionR = (inputSampleR-cMidPrevR) - accumulatorSampleR; + cMidPrevR = inputSampleR; + inputSampleR -= correctionR; + inputSampleR = (inputSampleR * cWet) + (drySampleR * cDry); + drySampleR = inputSampleR; + } + + if (dWet > 0.0) { + dMidL[19] = dMidL[18]; dMidL[18] = dMidL[17]; dMidL[17] = dMidL[16]; dMidL[16] = dMidL[15]; + dMidL[15] = dMidL[14]; dMidL[14] = dMidL[13]; dMidL[13] = dMidL[12]; dMidL[12] = dMidL[11]; + dMidL[11] = dMidL[10]; dMidL[10] = dMidL[9]; + dMidL[9] = dMidL[8]; dMidL[8] = dMidL[7]; dMidL[7] = dMidL[6]; dMidL[6] = dMidL[5]; + dMidL[5] = dMidL[4]; dMidL[4] = dMidL[3]; dMidL[3] = dMidL[2]; dMidL[2] = dMidL[1]; + dMidL[1] = dMidL[0]; dMidL[0] = accumulatorSampleL = (inputSampleL-dMidPrevL); + + accumulatorSampleL *= fMid[0]; + accumulatorSampleL += (dMidL[1] * fMid[1]); + accumulatorSampleL += (dMidL[2] * fMid[2]); + accumulatorSampleL += (dMidL[3] * fMid[3]); + accumulatorSampleL += (dMidL[4] * fMid[4]); + accumulatorSampleL += (dMidL[5] * fMid[5]); + accumulatorSampleL += (dMidL[6] * fMid[6]); + accumulatorSampleL += (dMidL[7] * fMid[7]); + accumulatorSampleL += (dMidL[8] * fMid[8]); + accumulatorSampleL += (dMidL[9] * fMid[9]); + accumulatorSampleL += (dMidL[10] * fMid[10]); + accumulatorSampleL += (dMidL[11] * fMid[11]); + accumulatorSampleL += (dMidL[12] * fMid[12]); + accumulatorSampleL += (dMidL[13] * fMid[13]); + accumulatorSampleL += (dMidL[14] * fMid[14]); + accumulatorSampleL += (dMidL[15] * fMid[15]); + accumulatorSampleL += (dMidL[16] * fMid[16]); + accumulatorSampleL += (dMidL[17] * fMid[17]); + accumulatorSampleL += (dMidL[18] * fMid[18]); + accumulatorSampleL += (dMidL[19] * fMid[19]); + //we are doing our repetitive calculations on a separate value + correctionL = (inputSampleL-dMidPrevL) - accumulatorSampleL; + dMidPrevL = inputSampleL; + inputSampleL -= correctionL; + inputSampleL = (inputSampleL * dWet) + (drySampleL * dDry); + drySampleL = inputSampleL; + + dMidR[19] = dMidR[18]; dMidR[18] = dMidR[17]; dMidR[17] = dMidR[16]; dMidR[16] = dMidR[15]; + dMidR[15] = dMidR[14]; dMidR[14] = dMidR[13]; dMidR[13] = dMidR[12]; dMidR[12] = dMidR[11]; + dMidR[11] = dMidR[10]; dMidR[10] = dMidR[9]; + dMidR[9] = dMidR[8]; dMidR[8] = dMidR[7]; dMidR[7] = dMidR[6]; dMidR[6] = dMidR[5]; + dMidR[5] = dMidR[4]; dMidR[4] = dMidR[3]; dMidR[3] = dMidR[2]; dMidR[2] = dMidR[1]; + dMidR[1] = dMidR[0]; dMidR[0] = accumulatorSampleR = (inputSampleR-dMidPrevR); + + accumulatorSampleR *= fMid[0]; + accumulatorSampleR += (dMidR[1] * fMid[1]); + accumulatorSampleR += (dMidR[2] * fMid[2]); + accumulatorSampleR += (dMidR[3] * fMid[3]); + accumulatorSampleR += (dMidR[4] * fMid[4]); + accumulatorSampleR += (dMidR[5] * fMid[5]); + accumulatorSampleR += (dMidR[6] * fMid[6]); + accumulatorSampleR += (dMidR[7] * fMid[7]); + accumulatorSampleR += (dMidR[8] * fMid[8]); + accumulatorSampleR += (dMidR[9] * fMid[9]); + accumulatorSampleR += (dMidR[10] * fMid[10]); + accumulatorSampleR += (dMidR[11] * fMid[11]); + accumulatorSampleR += (dMidR[12] * fMid[12]); + accumulatorSampleR += (dMidR[13] * fMid[13]); + accumulatorSampleR += (dMidR[14] * fMid[14]); + accumulatorSampleR += (dMidR[15] * fMid[15]); + accumulatorSampleR += (dMidR[16] * fMid[16]); + accumulatorSampleR += (dMidR[17] * fMid[17]); + accumulatorSampleR += (dMidR[18] * fMid[18]); + accumulatorSampleR += (dMidR[19] * fMid[19]); + //we are doing our repetitive calculations on a separate value + correctionR = (inputSampleR-dMidPrevR) - accumulatorSampleR; + dMidPrevR = inputSampleR; + inputSampleR -= correctionR; + inputSampleR = (inputSampleR * dWet) + (drySampleR * dDry); + drySampleR = inputSampleR; + } + + //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 |