aboutsummaryrefslogtreecommitdiffstats
path: root/plugins/MacVST/DeRez2/source/DeRez2Proc.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'plugins/MacVST/DeRez2/source/DeRez2Proc.cpp')
-rwxr-xr-xplugins/MacVST/DeRez2/source/DeRez2Proc.cpp364
1 files changed, 364 insertions, 0 deletions
diff --git a/plugins/MacVST/DeRez2/source/DeRez2Proc.cpp b/plugins/MacVST/DeRez2/source/DeRez2Proc.cpp
new file mode 100755
index 0000000..889fd79
--- /dev/null
+++ b/plugins/MacVST/DeRez2/source/DeRez2Proc.cpp
@@ -0,0 +1,364 @@
+/* ========================================
+ * DeRez2 - DeRez2.h
+ * Copyright (c) 2016 airwindows, All rights reserved
+ * ======================================== */
+
+#ifndef __DeRez2_H
+#include "DeRez2.h"
+#endif
+
+void DeRez2::processReplacing(float **inputs, float **outputs, VstInt32 sampleFrames)
+{
+ float* in1 = inputs[0];
+ float* in2 = inputs[1];
+ float* out1 = outputs[0];
+ float* out2 = outputs[1];
+
+ double targetA = pow(A,3)+0.0005;
+ if (targetA > 1.0) targetA = 1.0;
+ double soften = (1.0 + targetA)/2;
+ double targetB = pow(1.0-B,3) / 3;
+ double hard = C;
+ double wet = D;
+
+ double overallscale = 1.0;
+ overallscale /= 44100.0;
+ overallscale *= getSampleRate();
+ targetA /= overallscale;
+
+ 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;
+
+
+ incrementA = ((incrementA*999.0)+targetA)/1000.0;
+ incrementB = ((incrementB*999.0)+targetB)/1000.0;
+ //incrementA is the frequency derez
+ //incrementB is the bit depth derez
+ position += incrementA;
+
+ long double outputSampleL = heldSampleL;
+ long double outputSampleR = heldSampleR;
+ if (position > 1.0)
+ {
+ position -= 1.0;
+ heldSampleL = (lastSampleL * position) + (inputSampleL * (1.0-position));
+ outputSampleL = (outputSampleL * (1.0-soften)) + (heldSampleL * soften);
+ //softens the edge of the derez
+ heldSampleR = (lastSampleR * position) + (inputSampleR * (1.0-position));
+ outputSampleR = (outputSampleR * (1.0-soften)) + (heldSampleR * soften);
+ //softens the edge of the derez
+ }
+ inputSampleL = outputSampleL;
+ inputSampleR = outputSampleR;
+
+ long double tempL = inputSampleL;
+ long double tempR = inputSampleR;
+
+ if (inputSampleL != lastOutputSampleL) {
+ tempL = inputSampleL;
+ inputSampleL = (inputSampleL * hard) + (lastDrySampleL * (1.0-hard));
+ //transitions get an intermediate dry sample
+ lastOutputSampleL = tempL; //only one intermediate sample
+ } else {
+ lastOutputSampleL = inputSampleL;
+ }
+
+ if (inputSampleR != lastOutputSampleR) {
+ tempR = inputSampleR;
+ inputSampleR = (inputSampleR * hard) + (lastDrySampleR * (1.0-hard));
+ //transitions get an intermediate dry sample
+ lastOutputSampleR = tempR; //only one intermediate sample
+ } else {
+ lastOutputSampleR = inputSampleR;
+ }
+
+ lastDrySampleL = drySampleL;
+ lastDrySampleR = drySampleR;
+ //freq section of soft/hard interpolates dry samples
+
+ tempL = inputSampleL;
+ tempR = inputSampleR;
+
+ if (inputSampleL > 1.0) inputSampleL = 1.0;
+ if (inputSampleL < -1.0) inputSampleL = -1.0;
+
+ if (inputSampleR > 1.0) inputSampleR = 1.0;
+ if (inputSampleR < -1.0) inputSampleR = -1.0;
+
+ if (inputSampleL > 0) inputSampleL = log(1.0+(255*fabs(inputSampleL))) / log(256);
+ if (inputSampleL < 0) inputSampleL = -log(1.0+(255*fabs(inputSampleL))) / log(256);
+
+ if (inputSampleR > 0) inputSampleR = log(1.0+(255*fabs(inputSampleR))) / log(256);
+ if (inputSampleR < 0) inputSampleR = -log(1.0+(255*fabs(inputSampleR))) / log(256);
+
+ inputSampleL = (tempL * hard) + (inputSampleL * (1.0-hard));
+ inputSampleR = (tempR * hard) + (inputSampleR * (1.0-hard)); //uLaw encode as part of soft/hard
+
+ tempL = inputSampleL;
+ tempR = inputSampleR;
+
+ if (incrementB > 0.0005)
+ {
+ if (inputSampleL > 0)
+ {
+ tempL = inputSampleL;
+ while (tempL > 0) {tempL -= incrementB;}
+ inputSampleL -= tempL;
+ //it's below 0 so subtracting adds the remainder
+ }
+ if (inputSampleR > 0)
+ {
+ tempR = inputSampleR;
+ while (tempR > 0) {tempR -= incrementB;}
+ inputSampleR -= tempR;
+ //it's below 0 so subtracting adds the remainder
+ }
+
+ if (inputSampleL < 0)
+ {
+ tempL = inputSampleL;
+ while (tempL < 0) {tempL += incrementB;}
+ inputSampleL -= tempL;
+ //it's above 0 so subtracting subtracts the remainder
+ }
+ if (inputSampleR < 0)
+ {
+ tempR = inputSampleR;
+ while (tempR < 0) {tempR += incrementB;}
+ inputSampleR -= tempR;
+ //it's above 0 so subtracting subtracts the remainder
+ }
+
+ inputSampleL *= (1.0 - incrementB);
+ inputSampleR *= (1.0 - incrementB);
+ }
+
+ tempL = inputSampleL;
+ tempR = inputSampleR;
+
+ if (inputSampleL > 1.0) inputSampleL = 1.0;
+ if (inputSampleL < -1.0) inputSampleL = -1.0;
+
+ if (inputSampleR > 1.0) inputSampleR = 1.0;
+ if (inputSampleR < -1.0) inputSampleR = -1.0;
+
+ if (inputSampleL > 0) inputSampleL = (pow(256,fabs(inputSampleL))-1.0) / 255;
+ if (inputSampleL < 0) inputSampleL = -(pow(256,fabs(inputSampleL))-1.0) / 255;
+
+ if (inputSampleR > 0) inputSampleR = (pow(256,fabs(inputSampleR))-1.0) / 255;
+ if (inputSampleR < 0) inputSampleR = -(pow(256,fabs(inputSampleR))-1.0) / 255;
+
+ inputSampleL = (tempL * hard) + (inputSampleL * (1.0-hard));
+ inputSampleR = (tempR * hard) + (inputSampleR * (1.0-hard)); //uLaw decode as part of soft/hard
+
+ if (wet !=1.0) {
+ inputSampleL = (inputSampleL * wet) + (drySampleL * (1.0-wet));
+ inputSampleR = (inputSampleR * wet) + (drySampleR * (1.0-wet));
+ }
+ //Dry/Wet control, defaults to the last slider
+
+ lastSampleL = drySampleL;
+ lastSampleR = drySampleR;
+
+ //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 DeRez2::processDoubleReplacing(double **inputs, double **outputs, VstInt32 sampleFrames)
+{
+ double* in1 = inputs[0];
+ double* in2 = inputs[1];
+ double* out1 = outputs[0];
+ double* out2 = outputs[1];
+
+ double targetA = pow(A,3)+0.0005;
+ if (targetA > 1.0) targetA = 1.0;
+ double soften = (1.0 + targetA)/2;
+ double targetB = pow(1.0-B,3) / 3;
+ double hard = C;
+ double wet = D;
+
+ double overallscale = 1.0;
+ overallscale /= 44100.0;
+ overallscale *= getSampleRate();
+ targetA /= overallscale;
+
+ 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;
+
+
+ incrementA = ((incrementA*999.0)+targetA)/1000.0;
+ incrementB = ((incrementB*999.0)+targetB)/1000.0;
+ //incrementA is the frequency derez
+ //incrementB is the bit depth derez
+ position += incrementA;
+
+ long double outputSampleL = heldSampleL;
+ long double outputSampleR = heldSampleR;
+ if (position > 1.0)
+ {
+ position -= 1.0;
+ heldSampleL = (lastSampleL * position) + (inputSampleL * (1.0-position));
+ outputSampleL = (outputSampleL * (1.0-soften)) + (heldSampleL * soften);
+ //softens the edge of the derez
+ heldSampleR = (lastSampleR * position) + (inputSampleR * (1.0-position));
+ outputSampleR = (outputSampleR * (1.0-soften)) + (heldSampleR * soften);
+ //softens the edge of the derez
+ }
+ inputSampleL = outputSampleL;
+ inputSampleR = outputSampleR;
+
+ long double tempL = inputSampleL;
+ long double tempR = inputSampleR;
+
+ if (inputSampleL != lastOutputSampleL) {
+ tempL = inputSampleL;
+ inputSampleL = (inputSampleL * hard) + (lastDrySampleL * (1.0-hard));
+ //transitions get an intermediate dry sample
+ lastOutputSampleL = tempL; //only one intermediate sample
+ } else {
+ lastOutputSampleL = inputSampleL;
+ }
+
+ if (inputSampleR != lastOutputSampleR) {
+ tempR = inputSampleR;
+ inputSampleR = (inputSampleR * hard) + (lastDrySampleR * (1.0-hard));
+ //transitions get an intermediate dry sample
+ lastOutputSampleR = tempR; //only one intermediate sample
+ } else {
+ lastOutputSampleR = inputSampleR;
+ }
+
+ lastDrySampleL = drySampleL;
+ lastDrySampleR = drySampleR;
+ //freq section of soft/hard interpolates dry samples
+
+ tempL = inputSampleL;
+ tempR = inputSampleR;
+
+ if (inputSampleL > 1.0) inputSampleL = 1.0;
+ if (inputSampleL < -1.0) inputSampleL = -1.0;
+
+ if (inputSampleR > 1.0) inputSampleR = 1.0;
+ if (inputSampleR < -1.0) inputSampleR = -1.0;
+
+ if (inputSampleL > 0) inputSampleL = log(1.0+(255*fabs(inputSampleL))) / log(256);
+ if (inputSampleL < 0) inputSampleL = -log(1.0+(255*fabs(inputSampleL))) / log(256);
+
+ if (inputSampleR > 0) inputSampleR = log(1.0+(255*fabs(inputSampleR))) / log(256);
+ if (inputSampleR < 0) inputSampleR = -log(1.0+(255*fabs(inputSampleR))) / log(256);
+
+ inputSampleL = (tempL * hard) + (inputSampleL * (1.0-hard));
+ inputSampleR = (tempR * hard) + (inputSampleR * (1.0-hard)); //uLaw encode as part of soft/hard
+
+ tempL = inputSampleL;
+ tempR = inputSampleR;
+
+ if (incrementB > 0.0005)
+ {
+ if (inputSampleL > 0)
+ {
+ tempL = inputSampleL;
+ while (tempL > 0) {tempL -= incrementB;}
+ inputSampleL -= tempL;
+ //it's below 0 so subtracting adds the remainder
+ }
+ if (inputSampleR > 0)
+ {
+ tempR = inputSampleR;
+ while (tempR > 0) {tempR -= incrementB;}
+ inputSampleR -= tempR;
+ //it's below 0 so subtracting adds the remainder
+ }
+
+ if (inputSampleL < 0)
+ {
+ tempL = inputSampleL;
+ while (tempL < 0) {tempL += incrementB;}
+ inputSampleL -= tempL;
+ //it's above 0 so subtracting subtracts the remainder
+ }
+ if (inputSampleR < 0)
+ {
+ tempR = inputSampleR;
+ while (tempR < 0) {tempR += incrementB;}
+ inputSampleR -= tempR;
+ //it's above 0 so subtracting subtracts the remainder
+ }
+
+ inputSampleL *= (1.0 - incrementB);
+ inputSampleR *= (1.0 - incrementB);
+ }
+
+ tempL = inputSampleL;
+ tempR = inputSampleR;
+
+ if (inputSampleL > 1.0) inputSampleL = 1.0;
+ if (inputSampleL < -1.0) inputSampleL = -1.0;
+
+ if (inputSampleR > 1.0) inputSampleR = 1.0;
+ if (inputSampleR < -1.0) inputSampleR = -1.0;
+
+ if (inputSampleL > 0) inputSampleL = (pow(256,fabs(inputSampleL))-1.0) / 255;
+ if (inputSampleL < 0) inputSampleL = -(pow(256,fabs(inputSampleL))-1.0) / 255;
+
+ if (inputSampleR > 0) inputSampleR = (pow(256,fabs(inputSampleR))-1.0) / 255;
+ if (inputSampleR < 0) inputSampleR = -(pow(256,fabs(inputSampleR))-1.0) / 255;
+
+ inputSampleL = (tempL * hard) + (inputSampleL * (1.0-hard));
+ inputSampleR = (tempR * hard) + (inputSampleR * (1.0-hard)); //uLaw decode as part of soft/hard
+
+ if (wet !=1.0) {
+ inputSampleL = (inputSampleL * wet) + (drySampleL * (1.0-wet));
+ inputSampleR = (inputSampleR * wet) + (drySampleR * (1.0-wet));
+ }
+ //Dry/Wet control, defaults to the last slider
+
+ lastSampleL = drySampleL;
+ lastSampleR = drySampleR;
+
+ //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++;
+ }
+}
generated by cgit v1.2.3 (git 2.25.1) at 2024-06-26 13:33:07 +0000