github-starred/airwindows
2
0
Fork 0
mirror of https://github.com/airwindows/airwindows.git synced 2026-05-15 06:05:55 -06:00
Code Issues 33 Projects Packages Wiki Activity Actions
airwindows/plugins/MacSignedVST/PunchyDeluxe/source/PunchyDeluxeProc.cpp
Christopher Johnson 2335148573 BezEQ3
2026-05-09 20:09:00 -04:00

344 lines
15 KiB
C++
Executable file
Raw Blame History

/* ========================================
* PunchyDeluxe - PunchyDeluxe.h
* Copyright (c) airwindows, Airwindows uses the MIT license
* ======================================== */
#ifndef __PunchyDeluxe_H
#include "PunchyDeluxe.h"
#endif
void PunchyDeluxe::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 drive = 1.0;
double pad = 1.0;
angG[1] = sqrt(A*2.0);
if (pad > angG[1]) pad = angG[1];
if (drive < angG[1]) drive = angG[1];
angG[0] = (1.0+angG[1])*0.5; //if boost, boost half
if (angG[0] > angG[1]) angG[0] = angG[1];
angG[2] = sqrt(B*2.0);
if (pad > angG[2]) pad = angG[2];
if (drive < angG[2]) drive = angG[2];
angG[3] = sqrt(C*2.0);
if (pad > angG[3]) pad = angG[3];
if (drive < angG[3]) drive = angG[3];
angG[4] = sqrt(D*2.0);
if (pad > angG[4]) pad = angG[4];
if (drive < angG[4]) drive = angG[4];
angG[5] = sqrt(E*2.0);
if (pad > angG[5]) pad = angG[5];
if (drive < angG[5]) drive = angG[5];
angG[6] = sqrt(F*2.0);
if (pad > angG[6]) pad = angG[6];
if (drive < angG[6]) drive = angG[6];
angG[7] = sqrt(G*2.0);
if (pad > angG[7]) pad = angG[7];
if (drive < angG[7]) drive = angG[7];
angG[8] = sqrt(H*2.0);
if (pad > angG[8]) pad = angG[8];
if (drive < angG[8]) drive = angG[8];
angG[9] = sqrt(I*2.0);
if (pad > angG[9]) pad = angG[9];
if (drive < angG[9]) drive = angG[9];
angG[10] = sqrt(J*2.0);
if (pad > angG[10]) pad = angG[10];
if (drive < angG[10]) drive = angG[10];
angG[11] = 1.0;
angG[12] = 1.0;
pad = sqrt(pad);
int poles = (int)(drive*11.0);
double tune = 0.618+(overallscale*0.0055);
while (--sampleFrames >= 0)
{
double inputSampleL = *in1;
double inputSampleR = *in2;
if (fabs(inputSampleL)<1.18e-23) inputSampleL = fpdL * 1.18e-17;
if (fabs(inputSampleR)<1.18e-23) inputSampleR = fpdR * 1.18e-17;
for (int x = 0; x < poles; x++) {
double fr = (0.92/overallscale)+(overallscale*0.01);
double band = inputSampleL; inputSampleL = 0.0;
for (int y = 0; y < 12; y++) {
angAL[x][y] = (angAL[x][y]*(1.0-fr)) + ((band-angSL[x][y])*fr);
double temp = band; band = ((angSL[x][y]+(angAL[x][y]*fr)) * (1.0-fr))+(band*fr);
angSL[x][y] = ((angSL[x][y]+(angAL[x][y]*fr)) * (1.0-fr))+(band*fr);
inputSampleL += ((temp-band)*angG[y]);
fr *= tune;
}
inputSampleL += band;
inputSampleL *= drive;
inputSampleL = fmin(fmax(inputSampleL,-M_PI_2),M_PI_2);
long double X = inputSampleL; X *= X; //long double for even
long double temp = inputSampleL * X; //the initial multiplies
inputSampleL -= temp*0.16666666666666666666666666666666666; temp *= X;
inputSampleL += temp*0.00833333333333333333333333333333333; temp *= X;
inputSampleL -= temp*0.00019841269841269841269841269841269; temp *= X;
inputSampleL += temp*0.00000275573192239858906525573192239; temp *= X;
inputSampleL -= temp*0.00000002505210838544171877521083854; temp *= X;
inputSampleL += temp*0.00000000016059043836821614599392377; temp *= X;
inputSampleL -= temp*0.00000000000076471637318198164759011; temp *= X;
inputSampleL += temp*0.00000000000000281145725434552076319; temp *= X;
inputSampleL -= temp*0.00000000000000000822063524662432971; temp *= X;
inputSampleL += temp*0.00000000000000000001957294106339126;
fr = (0.92/overallscale)+(overallscale*0.01);
band = inputSampleR; inputSampleR = 0.0;
for (int y = 0; y < 12; y++) {
angAR[x][y] = (angAR[x][y]*(1.0-fr)) + ((band-angSR[x][y])*fr);
double temp = band; band = ((angSR[x][y]+(angAR[x][y]*fr)) * (1.0-fr))+(band*fr);
angSR[x][y] = ((angSR[x][y]+(angAR[x][y]*fr)) * (1.0-fr))+(band*fr);
inputSampleR += ((temp-band)*angG[y]);
fr *= tune;
}
inputSampleR += band;
inputSampleR *= drive;
inputSampleR = fmin(fmax(inputSampleR,-M_PI_2),M_PI_2);
X = inputSampleR; X *= X; //long double for even
temp = inputSampleR * X; //the initial multiplies
inputSampleR -= temp*0.16666666666666666666666666666666666; temp *= X;
inputSampleR += temp*0.00833333333333333333333333333333333; temp *= X;
inputSampleR -= temp*0.00019841269841269841269841269841269; temp *= X;
inputSampleR += temp*0.00000275573192239858906525573192239; temp *= X;
inputSampleR -= temp*0.00000002505210838544171877521083854; temp *= X;
inputSampleR += temp*0.00000000016059043836821614599392377; temp *= X;
inputSampleR -= temp*0.00000000000076471637318198164759011; temp *= X;
inputSampleR += temp*0.00000000000000281145725434552076319; temp *= X;
inputSampleR -= temp*0.00000000000000000822063524662432971; temp *= X;
inputSampleR += temp*0.00000000000000000001957294106339126;
//retain mantissa of a long double increasing power function
//long double probably doesn't handle more than 36 digits or so
}
if (pad < 1.0) {
inputSampleL *= pad;
inputSampleR *= pad;
}
inputSampleL = fmin(fmax(inputSampleL,-M_PI_2),M_PI_2);
inputSampleR = fmin(fmax(inputSampleR,-M_PI_2),M_PI_2);
long double X = inputSampleL; X *= X; //long double for even
long double temp = inputSampleL * X; //the initial multiplies
inputSampleL -= temp*0.16666666666666666666666666666666666; temp *= X;
inputSampleL += temp*0.00833333333333333333333333333333333; temp *= X;
inputSampleL -= temp*0.00019841269841269841269841269841269; temp *= X;
inputSampleL += temp*0.00000275573192239858906525573192239; temp *= X;
inputSampleL -= temp*0.00000002505210838544171877521083854; temp *= X;
inputSampleL += temp*0.00000000016059043836821614599392377; temp *= X;
inputSampleL -= temp*0.00000000000076471637318198164759011; temp *= X;
inputSampleL += temp*0.00000000000000281145725434552076319; temp *= X;
inputSampleL -= temp*0.00000000000000000822063524662432971; temp *= X;
inputSampleL += temp*0.00000000000000000001957294106339126;
X = inputSampleR; X *= X; //long double for even
temp = inputSampleR * X; //the initial multiplies
inputSampleR -= temp*0.16666666666666666666666666666666666; temp *= X;
inputSampleR += temp*0.00833333333333333333333333333333333; temp *= X;
inputSampleR -= temp*0.00019841269841269841269841269841269; temp *= X;
inputSampleR += temp*0.00000275573192239858906525573192239; temp *= X;
inputSampleR -= temp*0.00000002505210838544171877521083854; temp *= X;
inputSampleR += temp*0.00000000016059043836821614599392377; temp *= X;
inputSampleR -= temp*0.00000000000076471637318198164759011; temp *= X;
inputSampleR += temp*0.00000000000000281145725434552076319; temp *= X;
inputSampleR -= temp*0.00000000000000000822063524662432971; temp *= X;
inputSampleR += temp*0.00000000000000000001957294106339126;
//retain mantissa of a long double increasing power function
//long double probably doesn't handle more than 36 digits or so
//begin 32 bit stereo floating point dither
int expon; frexpf((float)inputSampleL, &expon);
fpdL ^= fpdL << 13; fpdL ^= fpdL >> 17; fpdL ^= fpdL << 5;
inputSampleL += ((double(fpdL)-uint32_t(0x7fffffff)) * 3.553e-44l * pow(2,expon+62));
frexpf((float)inputSampleR, &expon);
fpdR ^= fpdR << 13; fpdR ^= fpdR >> 17; fpdR ^= fpdR << 5;
if (fpdL-fpdR < 1073741824 || fpdR-fpdL < 1073741824) {
fpdR ^= fpdR << 13; fpdR ^= fpdR >> 17; fpdR ^= fpdR << 5;}
inputSampleR += ((double(fpdR)-uint32_t(0x7fffffff)) * 3.553e-44l * pow(2,expon+62));
//end 32 bit stereo floating point dither
*out1 = inputSampleL;
*out2 = inputSampleR;
in1++;
in2++;
out1++;
out2++;
}
}
void PunchyDeluxe::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 drive = 1.0;
double pad = 1.0;
angG[1] = sqrt(A*2.0);
if (pad > angG[1]) pad = angG[1];
if (drive < angG[1]) drive = angG[1];
angG[0] = (1.0+angG[1])*0.5; //if boost, boost half
if (angG[0] > angG[1]) angG[0] = angG[1];
angG[2] = sqrt(B*2.0);
if (pad > angG[2]) pad = angG[2];
if (drive < angG[2]) drive = angG[2];
angG[3] = sqrt(C*2.0);
if (pad > angG[3]) pad = angG[3];
if (drive < angG[3]) drive = angG[3];
angG[4] = sqrt(D*2.0);
if (pad > angG[4]) pad = angG[4];
if (drive < angG[4]) drive = angG[4];
angG[5] = sqrt(E*2.0);
if (pad > angG[5]) pad = angG[5];
if (drive < angG[5]) drive = angG[5];
angG[6] = sqrt(F*2.0);
if (pad > angG[6]) pad = angG[6];
if (drive < angG[6]) drive = angG[6];
angG[7] = sqrt(G*2.0);
if (pad > angG[7]) pad = angG[7];
if (drive < angG[7]) drive = angG[7];
angG[8] = sqrt(H*2.0);
if (pad > angG[8]) pad = angG[8];
if (drive < angG[8]) drive = angG[8];
angG[9] = sqrt(I*2.0);
if (pad > angG[9]) pad = angG[9];
if (drive < angG[9]) drive = angG[9];
angG[10] = sqrt(J*2.0);
if (pad > angG[10]) pad = angG[10];
if (drive < angG[10]) drive = angG[10];
angG[11] = 1.0;
angG[12] = 1.0;
pad = sqrt(pad);
int poles = (int)(drive*11.0);
double tune = 0.618+(overallscale*0.0055);
while (--sampleFrames >= 0)
{
double inputSampleL = *in1;
double inputSampleR = *in2;
if (fabs(inputSampleL)<1.18e-23) inputSampleL = fpdL * 1.18e-17;
if (fabs(inputSampleR)<1.18e-23) inputSampleR = fpdR * 1.18e-17;
for (int x = 0; x < poles; x++) {
double fr = (0.92/overallscale)+(overallscale*0.01);
double band = inputSampleL; inputSampleL = 0.0;
for (int y = 0; y < 12; y++) {
angAL[x][y] = (angAL[x][y]*(1.0-fr)) + ((band-angSL[x][y])*fr);
double temp = band; band = ((angSL[x][y]+(angAL[x][y]*fr)) * (1.0-fr))+(band*fr);
angSL[x][y] = ((angSL[x][y]+(angAL[x][y]*fr)) * (1.0-fr))+(band*fr);
inputSampleL += ((temp-band)*angG[y]);
fr *= tune;
}
inputSampleL += band;
inputSampleL *= drive;
inputSampleL = fmin(fmax(inputSampleL,-M_PI_2),M_PI_2);
long double X = inputSampleL; X *= X; //long double for even
long double temp = inputSampleL * X; //the initial multiplies
inputSampleL -= temp*0.16666666666666666666666666666666666; temp *= X;
inputSampleL += temp*0.00833333333333333333333333333333333; temp *= X;
inputSampleL -= temp*0.00019841269841269841269841269841269; temp *= X;
inputSampleL += temp*0.00000275573192239858906525573192239; temp *= X;
inputSampleL -= temp*0.00000002505210838544171877521083854; temp *= X;
inputSampleL += temp*0.00000000016059043836821614599392377; temp *= X;
inputSampleL -= temp*0.00000000000076471637318198164759011; temp *= X;
inputSampleL += temp*0.00000000000000281145725434552076319; temp *= X;
inputSampleL -= temp*0.00000000000000000822063524662432971; temp *= X;
inputSampleL += temp*0.00000000000000000001957294106339126;
fr = (0.92/overallscale)+(overallscale*0.01);
band = inputSampleR; inputSampleR = 0.0;
for (int y = 0; y < 12; y++) {
angAR[x][y] = (angAR[x][y]*(1.0-fr)) + ((band-angSR[x][y])*fr);
double temp = band; band = ((angSR[x][y]+(angAR[x][y]*fr)) * (1.0-fr))+(band*fr);
angSR[x][y] = ((angSR[x][y]+(angAR[x][y]*fr)) * (1.0-fr))+(band*fr);
inputSampleR += ((temp-band)*angG[y]);
fr *= tune;
}
inputSampleR += band;
inputSampleR *= drive;
inputSampleR = fmin(fmax(inputSampleR,-M_PI_2),M_PI_2);
X = inputSampleR; X *= X; //long double for even
temp = inputSampleR * X; //the initial multiplies
inputSampleR -= temp*0.16666666666666666666666666666666666; temp *= X;
inputSampleR += temp*0.00833333333333333333333333333333333; temp *= X;
inputSampleR -= temp*0.00019841269841269841269841269841269; temp *= X;
inputSampleR += temp*0.00000275573192239858906525573192239; temp *= X;
inputSampleR -= temp*0.00000002505210838544171877521083854; temp *= X;
inputSampleR += temp*0.00000000016059043836821614599392377; temp *= X;
inputSampleR -= temp*0.00000000000076471637318198164759011; temp *= X;
inputSampleR += temp*0.00000000000000281145725434552076319; temp *= X;
inputSampleR -= temp*0.00000000000000000822063524662432971; temp *= X;
inputSampleR += temp*0.00000000000000000001957294106339126;
//retain mantissa of a long double increasing power function
//long double probably doesn't handle more than 36 digits or so
}
if (pad < 1.0) {
inputSampleL *= pad;
inputSampleR *= pad;
}
inputSampleL = fmin(fmax(inputSampleL,-M_PI_2),M_PI_2);
inputSampleR = fmin(fmax(inputSampleR,-M_PI_2),M_PI_2);
long double X = inputSampleL; X *= X; //long double for even
long double temp = inputSampleL * X; //the initial multiplies
inputSampleL -= temp*0.16666666666666666666666666666666666; temp *= X;
inputSampleL += temp*0.00833333333333333333333333333333333; temp *= X;
inputSampleL -= temp*0.00019841269841269841269841269841269; temp *= X;
inputSampleL += temp*0.00000275573192239858906525573192239; temp *= X;
inputSampleL -= temp*0.00000002505210838544171877521083854; temp *= X;
inputSampleL += temp*0.00000000016059043836821614599392377; temp *= X;
inputSampleL -= temp*0.00000000000076471637318198164759011; temp *= X;
inputSampleL += temp*0.00000000000000281145725434552076319; temp *= X;
inputSampleL -= temp*0.00000000000000000822063524662432971; temp *= X;
inputSampleL += temp*0.00000000000000000001957294106339126;
X = inputSampleR; X *= X; //long double for even
temp = inputSampleR * X; //the initial multiplies
inputSampleR -= temp*0.16666666666666666666666666666666666; temp *= X;
inputSampleR += temp*0.00833333333333333333333333333333333; temp *= X;
inputSampleR -= temp*0.00019841269841269841269841269841269; temp *= X;
inputSampleR += temp*0.00000275573192239858906525573192239; temp *= X;
inputSampleR -= temp*0.00000002505210838544171877521083854; temp *= X;
inputSampleR += temp*0.00000000016059043836821614599392377; temp *= X;
inputSampleR -= temp*0.00000000000076471637318198164759011; temp *= X;
inputSampleR += temp*0.00000000000000281145725434552076319; temp *= X;
inputSampleR -= temp*0.00000000000000000822063524662432971; temp *= X;
inputSampleR += temp*0.00000000000000000001957294106339126;
//retain mantissa of a long double increasing power function
//long double probably doesn't handle more than 36 digits or so
//begin 64 bit stereo floating point dither
//int expon; frexp((double)inputSampleL, &expon);
fpdL ^= fpdL << 13; fpdL ^= fpdL >> 17; fpdL ^= fpdL << 5;
//inputSampleL += ((double(fpdL)-uint32_t(0x7fffffff)) * 3.553e-44l * pow(2,expon+62));
//frexp((double)inputSampleR, &expon);
fpdR ^= fpdR << 13; fpdR ^= fpdR >> 17; fpdR ^= fpdR << 5;
if (fpdL-fpdR < 1073741824 || fpdR-fpdL < 1073741824) {
fpdR ^= fpdR << 13; fpdR ^= fpdR >> 17; fpdR ^= fpdR << 5;}
//inputSampleR += ((double(fpdR)-uint32_t(0x7fffffff)) * 3.553e-44l * pow(2,expon+62));
//end 64 bit stereo floating point dither
*out1 = inputSampleL;
*out2 = inputSampleR;
in1++;
in2++;
out1++;
out2++;
}
}
Reference in a new issue View git blame Copy permalink