aboutsummaryrefslogtreecommitdiffstats
path: root/plugins/WinVST/Distortion/DistortionProc.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'plugins/WinVST/Distortion/DistortionProc.cpp')
-rwxr-xr-xplugins/WinVST/Distortion/DistortionProc.cpp250
1 files changed, 250 insertions, 0 deletions
diff --git a/plugins/WinVST/Distortion/DistortionProc.cpp b/plugins/WinVST/Distortion/DistortionProc.cpp
new file mode 100755
index 0000000..a31dcba
--- /dev/null
+++ b/plugins/WinVST/Distortion/DistortionProc.cpp
@@ -0,0 +1,250 @@
+/* ========================================
+ * Distortion - Distortion.h
+ * Copyright (c) 2016 airwindows, All rights reserved
+ * ======================================== */
+
+#ifndef __Distortion_H
+#include "Distortion.h"
+#endif
+
+void Distortion::processReplacing(float **inputs, float **outputs, VstInt32 sampleFrames)
+{
+ float* in1 = inputs[0];
+ float* in2 = inputs[1];
+ float* out1 = outputs[0];
+ float* out2 = outputs[1];
+
+ int stages = (int)floor(getSampleRate()/25000.0);
+ if (stages > 8) stages = 8;
+ double input = pow(10.0,((A-0.5)*24.0)/20.0);
+ int mode = (int) ( B * 4.999 );
+ double output = pow(10.0,((C-0.5)*24.0)/20.0);
+ double wet = D;
+
+ 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;
+
+ inputSampleL *= input;
+ inputSampleR *= input;
+
+ for (int x = 0; x < stages; x++) {
+ long double temp;
+ temp = (inputSampleL+previousInL[x])*0.5;
+ previousInL[x] = inputSampleL;
+ inputSampleL = temp;
+ temp = (inputSampleR+previousInR[x])*0.5;
+ previousInR[x] = inputSampleR;
+ inputSampleR = temp;
+ }
+
+ switch (mode)
+ {
+ case 0: //Density
+ if (inputSampleL > 1.570796326794897) inputSampleL = 1.570796326794897;
+ if (inputSampleL < -1.570796326794897) inputSampleL = -1.570796326794897;
+ if (inputSampleR > 1.570796326794897) inputSampleR = 1.570796326794897;
+ if (inputSampleR < -1.570796326794897) inputSampleR = -1.570796326794897;
+ //clip to 1.570796326794897 to reach maximum output
+ inputSampleL = sin(inputSampleL);
+ inputSampleR = sin(inputSampleR);
+ break;
+ case 1: //Drive
+ 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;
+ inputSampleL -= (inputSampleL * (fabs(inputSampleL) * 0.6) * (fabs(inputSampleL) * 0.6));
+ inputSampleR -= (inputSampleR * (fabs(inputSampleR) * 0.6) * (fabs(inputSampleR) * 0.6));
+ inputSampleL *= 1.5;
+ inputSampleR *= 1.5;
+ break;
+ case 2: //Spiral
+ if (inputSampleL > 1.2533141373155) inputSampleL = 1.2533141373155;
+ if (inputSampleL < -1.2533141373155) inputSampleL = -1.2533141373155;
+ if (inputSampleR > 1.2533141373155) inputSampleR = 1.2533141373155;
+ if (inputSampleR < -1.2533141373155) inputSampleR = -1.2533141373155;
+ //clip to 1.2533141373155 to reach maximum output
+ inputSampleL = sin(inputSampleL * fabs(inputSampleL)) / ((inputSampleL == 0.0) ?1:fabs(inputSampleL));
+ inputSampleR = sin(inputSampleR * fabs(inputSampleR)) / ((inputSampleR == 0.0) ?1:fabs(inputSampleR));
+ break;
+ case 3: //Mojo
+ long double mojo; mojo = pow(fabs(inputSampleL),0.25);
+ if (mojo > 0.0) inputSampleL = (sin(inputSampleL * mojo * M_PI * 0.5) / mojo) * 0.987654321;
+ mojo = pow(fabs(inputSampleR),0.25);
+ if (mojo > 0.0) inputSampleR = (sin(inputSampleR * mojo * M_PI * 0.5) / mojo) * 0.987654321;
+ //mojo is the one that flattens WAAAAY out very softly before wavefolding
+ break;
+ case 4: //Dyno
+ long double dyno; dyno = pow(fabs(inputSampleL),4);
+ if (dyno > 0.0) inputSampleL = (sin(inputSampleL * dyno) / dyno) * 1.1654321;
+ dyno = pow(fabs(inputSampleR),4);
+ if (dyno > 0.0) inputSampleR = (sin(inputSampleR * dyno) / dyno) * 1.1654321;
+ //dyno is the one that tries to raise peak energy
+ break;
+ }
+
+ for (int x = 1; x < (stages/2); x++) {
+ long double temp;
+ temp = (inputSampleL+previousOutL[x])*0.5;
+ previousOutL[x] = inputSampleL;
+ inputSampleL = temp;
+ temp = (inputSampleR+previousOutR[x])*0.5;
+ previousOutR[x] = inputSampleR;
+ inputSampleR = temp;
+ }
+
+ 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 += ((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 Distortion::processDoubleReplacing(double **inputs, double **outputs, VstInt32 sampleFrames)
+{
+ double* in1 = inputs[0];
+ double* in2 = inputs[1];
+ double* out1 = outputs[0];
+ double* out2 = outputs[1];
+
+ int stages = (int)floor(getSampleRate()/25000.0);
+ if (stages > 8) stages = 8;
+ double input = pow(10.0,((A-0.5)*24.0)/20.0);
+ int mode = (int) ( B * 4.999 );
+ double output = pow(10.0,((C-0.5)*24.0)/20.0);
+ double wet = D;
+
+ 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;
+
+ inputSampleL *= input;
+ inputSampleR *= input;
+
+ for (int x = 0; x < stages; x++) {
+ long double temp;
+ temp = (inputSampleL+previousInL[x])*0.5;
+ previousInL[x] = inputSampleL;
+ inputSampleL = temp;
+ temp = (inputSampleR+previousInR[x])*0.5;
+ previousInR[x] = inputSampleR;
+ inputSampleR = temp;
+ }
+
+ switch (mode)
+ {
+ case 0: //Density
+ if (inputSampleL > 1.570796326794897) inputSampleL = 1.570796326794897;
+ if (inputSampleL < -1.570796326794897) inputSampleL = -1.570796326794897;
+ if (inputSampleR > 1.570796326794897) inputSampleR = 1.570796326794897;
+ if (inputSampleR < -1.570796326794897) inputSampleR = -1.570796326794897;
+ //clip to 1.570796326794897 to reach maximum output
+ inputSampleL = sin(inputSampleL);
+ inputSampleR = sin(inputSampleR);
+ break;
+ case 1: //Drive
+ 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;
+ inputSampleL -= (inputSampleL * (fabs(inputSampleL) * 0.6) * (fabs(inputSampleL) * 0.6));
+ inputSampleR -= (inputSampleR * (fabs(inputSampleR) * 0.6) * (fabs(inputSampleR) * 0.6));
+ inputSampleL *= 1.5;
+ inputSampleR *= 1.5;
+ break;
+ case 2: //Spiral
+ if (inputSampleL > 1.2533141373155) inputSampleL = 1.2533141373155;
+ if (inputSampleL < -1.2533141373155) inputSampleL = -1.2533141373155;
+ if (inputSampleR > 1.2533141373155) inputSampleR = 1.2533141373155;
+ if (inputSampleR < -1.2533141373155) inputSampleR = -1.2533141373155;
+ //clip to 1.2533141373155 to reach maximum output
+ inputSampleL = sin(inputSampleL * fabs(inputSampleL)) / ((inputSampleL == 0.0) ?1:fabs(inputSampleL));
+ inputSampleR = sin(inputSampleR * fabs(inputSampleR)) / ((inputSampleR == 0.0) ?1:fabs(inputSampleR));
+ break;
+ case 3: //Mojo
+ long double mojo; mojo = pow(fabs(inputSampleL),0.25);
+ if (mojo > 0.0) inputSampleL = (sin(inputSampleL * mojo * M_PI * 0.5) / mojo) * 0.987654321;
+ mojo = pow(fabs(inputSampleR),0.25);
+ if (mojo > 0.0) inputSampleR = (sin(inputSampleR * mojo * M_PI * 0.5) / mojo) * 0.987654321;
+ //mojo is the one that flattens WAAAAY out very softly before wavefolding
+ break;
+ case 4: //Dyno
+ long double dyno; dyno = pow(fabs(inputSampleL),4);
+ if (dyno > 0.0) inputSampleL = (sin(inputSampleL * dyno) / dyno) * 1.1654321;
+ dyno = pow(fabs(inputSampleR),4);
+ if (dyno > 0.0) inputSampleR = (sin(inputSampleR * dyno) / dyno) * 1.1654321;
+ //dyno is the one that tries to raise peak energy
+ break;
+ }
+
+ for (int x = 1; x < (stages/2); x++) {
+ long double temp;
+ temp = (inputSampleL+previousOutL[x])*0.5;
+ previousOutL[x] = inputSampleL;
+ inputSampleL = temp;
+ temp = (inputSampleR+previousOutR[x])*0.5;
+ previousOutR[x] = inputSampleR;
+ inputSampleR = temp;
+ }
+
+ 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 += ((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-05-14 18:19:52 +0000