/* ========================================
* BitGlitter - BitGlitter.h
* Copyright (c) 2016 airwindows, All rights reserved
* ======================================== */
#ifndef __BitGlitter_H
#include "BitGlitter.h"
#endif
void BitGlitter::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 factor = B+1.0;
factor = pow(factor,7)+2.0;
int divvy = (int)(factor*overallscale);
double rateA = 1.0 / divvy;
double rezA = 0.0016666666666667; //looks to be a fixed bitcrush
double rateB = 1.61803398875 / divvy;
double rezB = 0.0026666666666667; //looks to be a fixed bitcrush
double offset;
double ingain = pow(10.0,((A * 36.0)-18.0)/14.0); //add adjustment factor
double outgain = pow(10.0,((C * 36.0)-18.0)/14.0); //add adjustment factor
double wet = D;
while (--sampleFrames >= 0)
{
long double inputSampleL = *in1;
long double inputSampleR = *in2;
long double drySampleL = inputSampleL;
long double drySampleR = inputSampleR;
//first, the distortion section
inputSampleL *= ingain;
inputSampleR *= ingain;
if (inputSampleL > 1.0) inputSampleL = 1.0;
if (inputSampleL < -1.0) inputSampleL = -1.0;
inputSampleL *= 1.2533141373155;
//clip to 1.2533141373155 to reach maximum output
inputSampleL = sin(inputSampleL * fabs(inputSampleL)) / ((inputSampleL == 0.0) ?1:fabs(inputSampleL));
if (inputSampleR > 1.0) inputSampleR = 1.0;
if (inputSampleR < -1.0) inputSampleR = -1.0;
inputSampleR *= 1.2533141373155;
//clip to 1.2533141373155 to reach maximum output
inputSampleR = sin(inputSampleR * fabs(inputSampleR)) / ((inputSampleR == 0.0) ?1:fabs(inputSampleR));
ataDrySampleL = inputSampleL;
ataHalfwaySampleL = (inputSampleL + ataLastSampleL ) / 2.0;
ataLastSampleL = inputSampleL;
//setting up crude oversampling
ataDrySampleR = inputSampleR;
ataHalfwaySampleR = (inputSampleR + ataLastSampleR ) / 2.0;
ataLastSampleR = inputSampleR;
//setting up crude oversampling
//begin raw sample L
positionAL += rateA;
long double outputSampleL = heldSampleAL;
if (positionAL > 1.0)
{
positionAL -= 1.0;
heldSampleAL = (lastSampleL * positionAL) + (inputSampleL * (1-positionAL));
outputSampleL = (outputSampleL * 0.5) + (heldSampleAL * 0.5);
//softens the edge of the derez
}
if (outputSampleL > 0)
{
offset = outputSampleL;
while (offset > 0) {offset -= rezA;}
outputSampleL -= offset;
//it's below 0 so subtracting adds the remainder
}
if (outputSampleL < 0)
{
offset = outputSampleL;
while (offset < 0) {offset += rezA;}
outputSampleL -= offset;
//it's above 0 so subtracting subtracts the remainder
}
outputSampleL *= (1.0 - rezA);
if (fabs(outputSampleL) < rezA) outputSampleL = 0.0;
inputSampleL = outputSampleL;
//end raw sample L
//begin raw sample R
positionAR += rateA;
long double outputSampleR = heldSampleAR;
if (positionAR > 1.0)
{
positionAR -= 1.0;
heldSampleAR = (lastSampleR * positionAR) + (inputSampleR * (1-positionAR));
outputSampleR = (outputSampleR * 0.5) + (heldSampleAR * 0.5);
//softens the edge of the derez
}
if (outputSampleR > 0)
{
offset = outputSampleR;
while (offset > 0) {offset -= rezA;}
outputSampleR -= offset;
//it's below 0 so subtracting adds the remainder
}
if (outputSampleR < 0)
{
offset = outputSampleR;
while (offset < 0) {offset += rezA;}
outputSampleR -= offset;
//it's above 0 so subtracting subtracts the remainder
}
outputSampleR *= (1.0 - rezA);
if (fabs(outputSampleR) < rezA) outputSampleR = 0.0;
inputSampleR = outputSampleR;
//end raw sample R
//begin interpolated sample L
positionBL += rateB;
outputSampleL = heldSampleBL;
if (positionBL > 1.0)
{
positionBL -= 1.0;
heldSampleBL = (lastSampleL * positionBL) + (ataHalfwaySampleL * (1-positionBL));
outputSampleL = (outputSampleL * 0.5) + (heldSampleBL * 0.5);
//softens the edge of the derez
}
if (outputSampleL > 0)
{
offset = outputSampleL;
while (offset > 0) {offset -= rezB;}
outputSampleL -= offset;
//it's below 0 so subtracting adds the remainder
}
if (outputSampleL < 0)
{
offset = outputSampleL;
while (offset < 0) {offset += rezB;}
outputSampleL -= offset;
//it's above 0 so subtracting subtracts the remainder
}
outputSampleL *= (1.0 - rezB);
if (fabs(outputSampleL) < rezB) outputSampleL = 0.0;
ataHalfwaySampleL = outputSampleL;
//end interpolated sample L
//begin interpolated sample R
positionBR += rateB;
outputSampleR = heldSampleBR;
if (positionBR > 1.0)
{
positionBR -= 1.0;
heldSampleBR = (lastSampleR * positionBR) + (ataHalfwaySampleR * (1-positionBR));
outputSampleR = (outputSampleR * 0.5) + (heldSampleBR * 0.5);
//softens the edge of the derez
}
if (outputSampleR > 0)
{
offset = outputSampleR;
while (offset > 0) {offset -= rezB;}
outputSampleR -= offset;
//it's below 0 so subtracting adds the remainder
}
if (outputSampleR < 0)
{
offset = outputSampleR;
while (offset < 0) {offset += rezB;}
outputSampleR -= offset;
//it's above 0 so subtracting subtracts the remainder
}
outputSampleR *= (1.0 - rezB);
if (fabs(outputSampleR) < rezB) outputSampleR = 0.0;
ataHalfwaySampleR = outputSampleR;
//end interpolated sample R
inputSampleL += ataHalfwaySampleL;
inputSampleL /= 2.0;
//plain old blend the two
inputSampleR += ataHalfwaySampleR;
inputSampleR /= 2.0;
//plain old blend the two
outputSampleL = (inputSampleL * (1.0-(wet/2))) + (lastOutputSampleL*(wet/2));
//darken to extent of wet in wet/dry, maximum 50%
lastOutputSampleL = inputSampleL;
outputSampleL *= outgain;
outputSampleR = (inputSampleR * (1.0-(wet/2))) + (lastOutputSampleR*(wet/2));
//darken to extent of wet in wet/dry, maximum 50%
lastOutputSampleR = inputSampleR;
outputSampleR *= outgain;
if (wet < 1.0) {
outputSampleL = (drySampleL * (1.0-wet)) + (outputSampleL * wet);
outputSampleR = (drySampleR * (1.0-wet)) + (outputSampleR * wet);
}
//stereo 32 bit dither, made small and tidy.
int expon; frexpf((float)outputSampleL, &expon);
long double dither = (rand()/(RAND_MAX*7.737125245533627e+25))*pow(2,expon+62);
outputSampleL += (dither-fpNShapeL); fpNShapeL = dither;
frexpf((float)outputSampleR, &expon);
dither = (rand()/(RAND_MAX*7.737125245533627e+25))*pow(2,expon+62);
outputSampleR += (dither-fpNShapeR); fpNShapeR = dither;
//end 32 bit dither
*out1 = outputSampleL;
*out2 = outputSampleR;
*in1++;
*in2++;
*out1++;
*out2++;
}
}
void BitGlitter::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 factor = B+1.0;
factor = pow(factor,7)+2.0;
int divvy = (int)(factor*overallscale);
double rateA = 1.0 / divvy;
double rezA = 0.0016666666666667; //looks to be a fixed bitcrush
double rateB = 1.61803398875 / divvy;
double rezB = 0.0026666666666667; //looks to be a fixed bitcrush
double offset;
double ingain = pow(10.0,((A * 36.0)-18.0)/14.0); //add adjustment factor
double outgain = pow(10.0,((C * 36.0)-18.0)/14.0); //add adjustment factor
double wet = D;
while (--sampleFrames >= 0)
{
long double inputSampleL = *in1;
long double inputSampleR = *in2;
long double drySampleL = inputSampleL;
long double drySampleR = inputSampleR;
//first, the distortion section
inputSampleL *= ingain;
inputSampleR *= ingain;
if (inputSampleL > 1.0) inputSampleL = 1.0;
if (inputSampleL < -1.0) inputSampleL = -1.0;
inputSampleL *= 1.2533141373155;
//clip to 1.2533141373155 to reach maximum output
inputSampleL = sin(inputSampleL * fabs(inputSampleL)) / ((inputSampleL == 0.0) ?1:fabs(inputSampleL));
if (inputSampleR > 1.0) inputSampleR = 1.0;
if (inputSampleR < -1.0) inputSampleR = -1.0;
inputSampleR *= 1.2533141373155;
//clip to 1.2533141373155 to reach maximum output
inputSampleR = sin(inputSampleR * fabs(inputSampleR)) / ((inputSampleR == 0.0) ?1:fabs(inputSampleR));
ataDrySampleL = inputSampleL;
ataHalfwaySampleL = (inputSampleL + ataLastSampleL ) / 2.0;
ataLastSampleL = inputSampleL;
//setting up crude oversampling
ataDrySampleR = inputSampleR;
ataHalfwaySampleR = (inputSampleR + ataLastSampleR ) / 2.0;
ataLastSampleR = inputSampleR;
//setting up crude oversampling
//begin raw sample L
positionAL += rateA;
long double outputSampleL = heldSampleAL;
if (positionAL > 1.0)
{
positionAL -= 1.0;
heldSampleAL = (lastSampleL * positionAL) + (inputSampleL * (1-positionAL));
outputSampleL = (outputSampleL * 0.5) + (heldSampleAL * 0.5);
//softens the edge of the derez
}
if (outputSampleL > 0)
{
offset = outputSampleL;
while (offset > 0) {offset -= rezA;}
outputSampleL -= offset;
//it's below 0 so subtracting adds the remainder
}
if (outputSampleL < 0)
{
offset = outputSampleL;
while (offset < 0) {offset += rezA;}
outputSampleL -= offset;
//it's above 0 so subtracting subtracts the remainder
}
outputSampleL *= (1.0 - rezA);
if (fabs(outputSampleL) < rezA) outputSampleL = 0.0;
inputSampleL = outputSampleL;
//end raw sample L
//begin raw sample R
positionAR += rateA;
long double outputSampleR = heldSampleAR;
if (positionAR > 1.0)
{
positionAR -= 1.0;
heldSampleAR = (lastSampleR * positionAR) + (inputSampleR * (1-positionAR));
outputSampleR = (outputSampleR * 0.5) + (heldSampleAR * 0.5);
//softens the edge of the derez
}
if (outputSampleR > 0)
{
offset = outputSampleR;
while (offset > 0) {offset -= rezA;}
outputSampleR -= offset;
//it's below 0 so subtracting adds the remainder
}
if (outputSampleR < 0)
{
offset = outputSampleR;
while (offset < 0) {offset += rezA;}
outputSampleR -= offset;
//it's above 0 so subtracting subtracts the remainder
}
outputSampleR *= (1.0 - rezA);
if (fabs(outputSampleR) < rezA) outputSampleR = 0.0;
inputSampleR = outputSampleR;
//end raw sample R
//begin interpolated sample L
positionBL += rateB;
outputSampleL = heldSampleBL;
if (positionBL > 1.0)
{
positionBL -= 1.0;
heldSampleBL = (lastSampleL * positionBL) + (ataHalfwaySampleL * (1-positionBL));
outputSampleL = (outputSampleL * 0.5) + (heldSampleBL * 0.5);
//softens the edge of the derez
}
if (outputSampleL > 0)
{
offset = outputSampleL;
while (offset > 0) {offset -= rezB;}
outputSampleL -= offset;
//it's below 0 so subtracting adds the remainder
}
if (outputSampleL < 0)
{
offset = outputSampleL;
while (offset < 0) {offset += rezB;}
outputSampleL -= offset;
//it's above 0 so subtracting subtracts the remainder
}
outputSampleL *= (1.0 - rezB);
if (fabs(outputSampleL) < rezB) outputSampleL = 0.0;
ataHalfwaySampleL = outputSampleL;
//end interpolated sample L
//begin interpolated sample R
positionBR += rateB;
outputSampleR = heldSampleBR;
if (positionBR > 1.0)
{
positionBR -= 1.0;
heldSampleBR = (lastSampleR * positionBR) + (ataHalfwaySampleR * (1-positionBR));
outputSampleR = (outputSampleR * 0.5) + (heldSampleBR * 0.5);
//softens the edge of the derez
}
if (outputSampleR > 0)
{
offset = outputSampleR;
while (offset > 0) {offset -= rezB;}
outputSampleR -= offset;
//it's below 0 so subtracting adds the remainder
}
if (outputSampleR < 0)
{
offset = outputSampleR;
while (offset < 0) {offset += rezB;}
outputSampleR -= offset;
//it's above 0 so subtracting subtracts the remainder
}
outputSampleR *= (1.0 - rezB);
if (fabs(outputSampleR) < rezB) outputSampleR = 0.0;
ataHalfwaySampleR = outputSampleR;
//end interpolated sample R
inputSampleL += ataHalfwaySampleL;
inputSampleL /= 2.0;
//plain old blend the two
inputSampleR += ataHalfwaySampleR;
inputSampleR /= 2.0;
//plain old blend the two
outputSampleL = (inputSampleL * (1.0-(wet/2))) + (lastOutputSampleL*(wet/2));
//darken to extent of wet in wet/dry, maximum 50%
lastOutputSampleL = inputSampleL;
outputSampleL *= outgain;
outputSampleR = (inputSampleR * (1.0-(wet/2))) + (lastOutputSampleR*(wet/2));
//darken to extent of wet in wet/dry, maximum 50%
lastOutputSampleR = inputSampleR;
outputSampleR *= outgain;
if (wet < 1.0) {
outputSampleL = (drySampleL * (1.0-wet)) + (outputSampleL * wet);
outputSampleR = (drySampleR * (1.0-wet)) + (outputSampleR * wet);
}
//stereo 64 bit dither, made small and tidy.
int expon; frexp((double)outputSampleL, &expon);
long double dither = (rand()/(RAND_MAX*7.737125245533627e+25))*pow(2,expon+62);
dither /= 536870912.0; //needs this to scale to 64 bit zone
outputSampleL += (dither-fpNShapeL); fpNShapeL = dither;
frexp((double)outputSampleR, &expon);
dither = (rand()/(RAND_MAX*7.737125245533627e+25))*pow(2,expon+62);
dither /= 536870912.0; //needs this to scale to 64 bit zone
outputSampleR += (dither-fpNShapeR); fpNShapeR = dither;
//end 64 bit dither
*out1 = outputSampleL;
*out2 = outputSampleR;
*in1++;
*in2++;
*out1++;
*out2++;
}
}
