From e360047054117d63fb579ec9231e9dc77c99f12a Mon Sep 17 00:00:00 2001 From: David Robillard Date: Fri, 8 Aug 2008 22:45:58 +0000 Subject: Add preliminary (library side only) LV2 port of MDA (open-sourced VST plugins). git-svn-id: http://svn.drobilla.net/lad/mda-lv2@1321 a436a847-0d15-0410-975c-d299462d15a1 --- src/mdaCombo.cpp | 492 +++++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 492 insertions(+) create mode 100644 src/mdaCombo.cpp (limited to 'src/mdaCombo.cpp') diff --git a/src/mdaCombo.cpp b/src/mdaCombo.cpp new file mode 100644 index 0000000..944ee62 --- /dev/null +++ b/src/mdaCombo.cpp @@ -0,0 +1,492 @@ +#include "mdaCombo.h" + +#include +#include + +AudioEffect *createEffectInstance(audioMasterCallback audioMaster) +{ + return new mdaCombo(audioMaster); +} + +mdaCombo::mdaCombo(audioMasterCallback audioMaster) : AudioEffectX(audioMaster, 1, 7) // programs, parameters +{ + //inits here! + fParam1 = 1.00f; //select + fParam2 = 0.50f; //drive + fParam3 = 0.50f; //bias + fParam4 = 0.50f; //output + fParam5 = 0.40f; //stereo + fParam6 = 0.00f; //hpf freq + fParam7 = 0.50f; //hpf reso + + size = 1024; + bufpos = 0; + buffer = new float[size]; + buffe2 = new float[size]; + + setNumInputs(2); + setNumOutputs(2); + setUniqueID("mdaX"); // identify here + DECLARE_LVZ_DEPRECATED(canMono) (); + canProcessReplacing(); + strcpy(programName, "Amp & Speaker Simulator"); + suspend(); // flush buffer + + setParameter(0, 0.f); //go and set initial values! +} + +bool mdaCombo::getProductString(char* text) { strcpy(text, "mda Combo"); return true; } +bool mdaCombo::getVendorString(char* text) { strcpy(text, "mda"); return true; } +bool mdaCombo::getEffectName(char* name) { strcpy(name, "Combo"); return true; } + +void mdaCombo::setParameter(LvzInt32 index, float value) +{ + switch(index) + { + case 0: fParam1 = value; break; + case 1: fParam2 = value; break; + case 2: fParam3 = value; break; + case 3: fParam4 = value; break; + case 4: fParam5 = value; break; + case 5: fParam6 = value; break; + case 6: fParam7 = value; break; + } + //calcs here + ster=0; if(fParam5>0.5) ster=1; + hpf = filterFreq(25.f); + switch(int(fParam1*6.9)) + { + case 0: trim = 0.5f; lpf = 0.f; //DI + mix1 = (float)0.0; mix2 = (float)0.0; + del1 = 0; del2 = 0; + break; + + case 1: trim = 0.53f; lpf = filterFreq(2700.f); //speaker sim + mix1 = (float)0.0; mix2 = (float)0.0; + del1 = 0; del2 = 0; + hpf = filterFreq(382.f); + break; + + case 2: trim = 1.10f; lpf = filterFreq(1685.f); //radio + mix1 = -1.70f; mix2 = 0.82f; + del1 = int(getSampleRate() / 6546.f); + del2 = int(getSampleRate() / 4315.f); + break; + + case 3: trim = 0.98f; lpf = filterFreq(1385.f); //mesa boogie 1" + mix1 = -0.53f; mix2 = 0.21f; + del1 = int(getSampleRate() / 7345.f); + del2 = int(getSampleRate() / 1193.f); + break; + + case 4: trim = 0.96f; lpf = filterFreq(1685.f); //mesa boogie 8" + mix1 = -0.85f; mix2 = 0.41f; + del1 = int(getSampleRate() / 6546.f); + del2 = int(getSampleRate() / 3315.f); + break; + + case 5: trim = 0.59f; lpf = lpf = filterFreq(2795.f); + mix1 = -0.29f; mix2 = 0.38f; //Marshall 4x12" celestion + del1 = int(getSampleRate() / 982.f); + del2 = int(getSampleRate() / 2402.f); + hpf = filterFreq(459.f); + break; + + case 6: trim = 0.30f; lpf = filterFreq(1744.f); //scooped-out metal + mix1 = -0.96f; mix2 = 1.6f; + del1 = int(getSampleRate() / 356.f); + del2 = int(getSampleRate() / 1263.f); + hpf = filterFreq(382.f); + break; + } + + mode = (fParam2<0.5)? 1 : 0; + if(mode) //soft clipping + { + drive = (float)pow(10.f, 2.f - 6.f * fParam2); + trim *= 0.55f + 150.f * (float)pow(fParam2,4.0f); + } + else //hard clipping + { + drive = 1.f; + clip = 11.7f - 16.f*fParam2; + if(fParam2>0.7) + { + drive = (float)pow(10.0f, 7.f*fParam2 - 4.9f); + clip = 0.5f; + } + } + bias = 1.2f * fParam3 - 0.6f; + if(fParam2>0.5) bias /= (1.f + 3.f * (fParam2-0.5f)); + else bias /= (1.f + 3.f * (0.5f-fParam2)); + trim *= (float)pow(10.f, 2.f*fParam4 - 1.f); + if(ster) trim *=2.f; + + hhf = fParam6; + hhq = 1.1f - fParam7; + if(fParam6>0.05f) drive = drive * (1 + 0.1f * drive); +} + +mdaCombo::~mdaCombo() +{ + if(buffer) delete [] buffer; + if(buffe2) delete [] buffe2; +} + +void mdaCombo::suspend() +{ + memset(buffer, 0, size * sizeof(float)); + memset(buffe2, 0, size * sizeof(float)); + + ff1=0.f; ff2=0.f; ff3=0.f; ff4=0.f; ff5=0.f; + ff6=0.f; ff7=0.f; ff8=0.f; ff9=0.f; ff10=0.f; + hh0 = hh1 = 0.0f; +} + +float mdaCombo::filterFreq(float hz) +{ + float j, k, r=0.999f; + + j = r * r - 1; + k = (float)(2.f - 2.f * r * r * cos(0.647f * hz / getSampleRate() )); + return (float)((sqrt(k*k - 4.f*j*j) - k) / (2.f*j)); +} + +void mdaCombo::setProgramName(char *name) +{ + strcpy(programName, name); +} + +void mdaCombo::getProgramName(char *name) +{ + strcpy(name, programName); +} + +float mdaCombo::getParameter(LvzInt32 index) +{ + float v=0; + + switch(index) + { + case 0: v = fParam1; break; + case 1: v = fParam2; break; + case 2: v = fParam3; break; + case 3: v = fParam4; break; + case 4: v = fParam5; break; + case 5: v = fParam6; break; + case 6: v = fParam7; break; + } + return v; +} + +void mdaCombo::getParameterName(LvzInt32 index, char *label) +{ + switch(index) + { + case 0: strcpy(label, "Model"); break; + case 1: strcpy(label, "Drive"); break; + case 2: strcpy(label, "Bias"); break; + case 3: strcpy(label, "Output"); break; + case 4: strcpy(label, "Process:"); break; + case 5: strcpy(label, "HPF Freq"); break; + case 6: strcpy(label, "HPF Reso"); break; + } +} + +#include +void long2string(long value, char *string) { sprintf(string, "%ld", value); } +void float2strng(float value, char *string) { sprintf(string, "%.2f", value); } + +void mdaCombo::getParameterDisplay(LvzInt32 index, char *text) +{ + switch(index) + { + case 0: switch(int(fParam1*6.9)) + { + case 0: strcpy(text, "D.I."); break; + case 1: strcpy(text, "Spkr Sim"); break; + //case 2: strcpy(text, "Phone"); break; + case 2: strcpy(text, "Radio"); break; + case 3: strcpy(text, "MB 1\""); break; + case 4: strcpy(text, "MB 8\""); break; + case 5: strcpy(text, "4x12 ^"); break; + case 6: strcpy(text, "4x12 >"); break; + } break; + + case 1: long2string((long)(200 * fParam2 - 100), text); break; + case 2: long2string((long)(200 * fParam3 - 100), text); break; + case 3: long2string((long)(40 * fParam4 - 20), text); break; + case 4: if(fParam5>0.5) strcpy(text, "STEREO"); + else strcpy(text, "MONO"); break; + case 5: long2string((long)(100 * fParam6), text); break; + case 6: long2string((long)(100 * fParam7), text); break; + } +} + +void mdaCombo::getParameterLabel(LvzInt32 index, char *label) +{ + switch(index) + { + case 0: strcpy(label, ""); break; + case 1: strcpy(label, "S <> H"); break; + case 2: strcpy(label, ""); break; + case 3: strcpy(label, "dB"); break; + case 4: strcpy(label, ""); break; + case 5: strcpy(label, "%"); break; + case 6: strcpy(label, "%"); break; + } +} + +//-------------------------------------------------------------------------------- +// process + +void mdaCombo::process(float **inputs, float **outputs, LvzInt32 sampleFrames) +{ + float *in1 = inputs[0]; + float *in2 = inputs[1]; + float *out1 = outputs[0]; + float *out2 = outputs[1]; + float a, b, c, d, trm, m1=mix1, m2=mix2, clp=clip; + float o=lpf, i=1.f-lpf, o2=hpf, i2=1.f-hpf, bi=bias, drv=drive; + float f1=ff1, f2=ff2, f3=ff3, f4=ff4, f5=ff5; + float a2, b2, f6=ff6, f7=ff7, f8=ff8, f9=ff9, f10=ff10; + float h0=hh0, h1=hh1; + long d1=del1, d2=del2, bp = bufpos; + + trm = trim * i * i * i * i; + + --in1; + --in2; + --out1; + --out2; + + if(fParam5>0.5) //stereo + { + while(--sampleFrames >= 0) + { + a = drv * (*++in1 + bi); a2 = drv * (*++in2 + bi); + c = out1[1]; + d = out2[1]; //process from here... + + if(mode) + { + b = (a>0.f)? a : -a; b2 = (a2>0.f)? a2 : -a2; + b = a / (1.f + b); b2 = a2 / (1.f + b2); + } + else + { + b = (a> clp)? clp : a; b2 = (a2> clp)? clp : a2; //distort + b = (a<-clp)? -clp : b; b2 = (a2<-clp)? -clp : b2; + } + + *(buffer + bp) = b; *(buffe2 + bp) = b2; + b += (m1* *(buffer+((bp+d1)%1000))) + (m2* *(buffer+((bp+d2)%1000))); + b2+= (m1* *(buffe2+((bp+d1)%1000))) + (m2* *(buffe2+((bp+d2)%1000))); + + f1 = o * f1 + trm * b; f6 = o * f6 + trm * b2; + f2 = o * f2 + f1; f7 = o * f7 + f6; + f3 = o * f3 + f2; f8 = o * f8 + f7; + f4 = o * f4 + f3; f9 = o * f9 + f8; //-24dB/oct filter + + f5 = o2 * f5 + i2 * f4; f10 = o2 * f10 + i2 * f9;//high pass + b = f4 - f5; b2 = f9 - f10; + + bp = (bp==0)? 999 : bp - 1; //buffer position + + *++out1 = c + b; + *++out2 = d + b2; + bp=bufpos; + } + } + else //mono + { + if(mode) //soft clip + { + while(--sampleFrames >= 0) + { + a = drv * (*++in1 + *++in2 + bi); + c = out1[1]; + d = out2[1]; //process from here... + + b = (a>0.f)? a : -a; + b = a / (1.f + b); + + *(buffer + bp) = b; + b += (m1* *(buffer+((bp+d1)%1000))) + (m2* *(buffer+((bp+d2)%1000))); + + f1 = o * f1 + trm * b; + f2 = o * f2 + f1; + f3 = o * f3 + f2; + f4 = o * f4 + f3; //-24dB/oct filter + + f5 = o2 * f5 + i2 * f4; //high pass + b = f4 - f5; + + bp = (bp==0)? 999 : bp - 1; //buffer position + + *++out1 = c + b; + *++out2 = d + b; + } + } + else //hard clip + { + while(--sampleFrames >= 0) + { + a = drv * (*++in1 + *++in2 + bi); + c = out1[1]; + d = out2[1]; //process from here... + + b = (a>clp)? clp : a; //distort + b = (a<-clp)? -clp : b; + + *(buffer + bp) = b; + b += (m1* *(buffer+((bp+d1)%1000))) + (m2* *(buffer+((bp+d2)%1000))); + + f1 = o * f1 + trm * b; + f2 = o * f2 + f1; + f3 = o * f3 + f2; + f4 = o * f4 + f3; //-24dB/oct filter + + f5 = o2 * f5 + i2 * f4; //high pass + b = f4 - f5; + + bp = (bp==0)? 999 : bp - 1; //buffer position + + *++out1 = c + b; + *++out2 = d + b; + } + } + } + bufpos = bp; + if(fabs(f1)<1.0e-10) { ff1=0.f; ff2=0.f; ff3=0.f; ff4=0.f; ff5=0.f; } + else { ff1=f1; ff2=f2; ff3=f3; ff4=f4; ff5=f5; } + if(fabs(f6)<1.0e-10) { ff6=0.f; ff7=0.f; ff8=0.f; ff9=0.f; ff10=0.f; } + else { ff6=f6; ff7=f7; ff8=f8; ff9=f9; ff10=f10; } + if(fabs(h0)<1.0e-10) { hh0 = hh1 = 0.0f; } else { hh0=h0; hh1=h1; } +} + +void mdaCombo::processReplacing(float **inputs, float **outputs, LvzInt32 sampleFrames) +{ + float *in1 = inputs[0]; + float *in2 = inputs[1]; + float *out1 = outputs[0]; + float *out2 = outputs[1]; + float a, b, trm, m1=mix1, m2=mix2, clp=clip; + float o=lpf, i=1.f-lpf, o2=hpf, bi=bias, drv=drive; + float f1=ff1, f2=ff2, f3=ff3, f4=ff4, f5=ff5; + float a2, b2, f6=ff6, f7=ff7, f8=ff8, f9=ff9, f10=ff10; + float hf=hhf, hq=hhq, h0=hh0, h1=hh1; + long d1=del1, d2=del2, bp = bufpos; + + trm = trim * i * i * i * i; + + --in1; + --in2; + --out1; + --out2; + + if(ster) //stereo + { + while(--sampleFrames >= 0) + { + a = drv * (*++in1 + bi); a2 = drv * (*++in2 + bi); + + if(mode) + { + b = (a>0.f)? a : -a; b2 = (a2>0.f)? a2 : -a2; + b = a / (1.f + b); b2 = a2 / (1.f + b2); + } + else + { + b = (a> clp)? clp : a; b2 = (a2> clp)? clp : a2; //distort + b = (a<-clp)? -clp : b; b2 = (a2<-clp)? -clp : b2; + } + + *(buffer + bp) = b; *(buffe2 + bp) = b2; + b += (m1* *(buffer+((bp+d1)%1000))) + (m2* *(buffer+((bp+d2)%1000))); + b2+= (m1* *(buffe2+((bp+d1)%1000))) + (m2* *(buffe2+((bp+d2)%1000))); + + f1 = o * f1 + trm * b; f6 = o * f6 + trm * b2; + f2 = o * f2 + f1; f7 = o * f7 + f6; + f3 = o * f3 + f2; f8 = o * f8 + f7; + f4 = o * f4 + f3; f9 = o * f9 + f8; //-24dB/oct filter + + f5 = o2 * (f5 - f4) + f4; f10 = o2 * (f10 - f9) + f9;//high pass + b = f4 - f5; b2 = f9 - f10; + + if(bp==0) bufpos=999; else bufpos=bp-1; + + *++out1 = b; + *++out2 = b2; + } + } + else //mono + { + if(mode) //soft clip + { + while(--sampleFrames >= 0) + { + a = drv * (*++in1 + *++in2 + bi); + + h0 += hf * (h1 + a); //resonant highpass (Chamberlin SVF) + h1 -= hf * (h0 + hq * h1); + a += h1; + + b = (a>0.f)? a : -a; + b = a / (1.f + b); + + *(buffer + bp) = b; + b += (m1* *(buffer+((bp+d1)%1000))) + (m2* *(buffer+((bp+d2)%1000))); + + f1 = o * f1 + trm * b; + f2 = o * f2 + f1; + f3 = o * f3 + f2; + f4 = o * f4 + f3; //-24dB/oct filter + + f5 = o2 * (f5 - f4) + f4; //high pass + b = f4 - f5; + + bp = (bp==0)? 999 : bp - 1; //buffer position + + *++out1 = b; + *++out2 = b; + } + } + else //hard clip + { + while(--sampleFrames >= 0) + { + a = drv * (*++in1 + *++in2 + bi); + + h0 += hf * (h1 + a); //resonant highpass (Chamberlin SVF) + h1 -= hf * (h0 + hq * h1); + a += h1; + + b = (a>clp)? clp : a; //distort + b = (a<-clp)? -clp : b; + + *(buffer + bp) = b; + b += (m1* *(buffer+((bp+d1)%1000))) + (m2* *(buffer+((bp+d2)%1000))); + + f1 = o * f1 + trm * b; + f2 = o * f2 + f1; + f3 = o * f3 + f2; + f4 = o * f4 + f3; //-24dB/oct filter + + f5 = o2 * (f5 - f4) + f4; //high pass //also want smile curve here... + b = f4 - f5; + + bp = (bp==0)? 999 : bp - 1; //buffer position + + *++out1 = b; + *++out2 = b; + } + } + } + bufpos = bp; + if(fabs(f1)<1.0e-10) { ff1=0.f; ff2=0.f; ff3=0.f; ff4=0.f; ff5=0.f; } + else { ff1=f1; ff2=f2; ff3=f3; ff4=f4; ff5=f5; } + if(fabs(f6)<1.0e-10) { ff6=0.f; ff7=0.f; ff8=0.f; ff9=0.f; ff10=0.f; } + else { ff6=f6; ff7=f7; ff8=f8; ff9=f9; ff10=f10; } + if(fabs(h0)<1.0e-10) { hh0 = hh1 = 0.0f; } else { hh0=h0; hh1=h1; } +} -- cgit v1.2.1