diff options
Diffstat (limited to 'plugins/WinVST/DeRez2/DeRez2Proc.cpp')
-rwxr-xr-x | plugins/WinVST/DeRez2/DeRez2Proc.cpp | 364 |
1 files changed, 364 insertions, 0 deletions
diff --git a/plugins/WinVST/DeRez2/DeRez2Proc.cpp b/plugins/WinVST/DeRez2/DeRez2Proc.cpp new file mode 100755 index 0000000..889fd79 --- /dev/null +++ b/plugins/WinVST/DeRez2/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++; + } +} |