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/mdaLeslie.cpp | 400 ++++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 400 insertions(+) create mode 100644 src/mdaLeslie.cpp (limited to 'src/mdaLeslie.cpp') diff --git a/src/mdaLeslie.cpp b/src/mdaLeslie.cpp new file mode 100644 index 0000000..97d2718 --- /dev/null +++ b/src/mdaLeslie.cpp @@ -0,0 +1,400 @@ +//Remove optimization /Og else wavelab crashes!!! + +#include "mdaLeslie.h" + +#include +#include +#include + +AudioEffect *createEffectInstance(audioMasterCallback audioMaster) +{ + return new mdaLeslie(audioMaster); +} + +mdaLeslieProgram::mdaLeslieProgram() +{ + fParam1 = 0.66f; + fParam7 = 0.50f; + fParam9 = 0.60f; + fParam4 = 0.70f; + fParam5 = 0.60f; + fParam6 = 0.70f; + fParam3 = 0.48f; + fParam2 = 0.50f; + fParam8 = 0.50f; + strcpy(name, "Leslie Simulator"); +} + +mdaLeslie::mdaLeslie(audioMasterCallback audioMaster) : AudioEffectX(audioMaster, 3, 9) // programs, parameters +{ + + size = 256; hpos = 0; + hbuf = new float[size]; + fbuf1 = fbuf2 = 0.0f; + twopi = 6.2831853f; + + setNumInputs(2); + setNumOutputs(2); + setUniqueID("mdaH"); // identify here + DECLARE_LVZ_DEPRECATED(canMono) (); + canProcessReplacing(); + suspend(); + + programs = new mdaLeslieProgram[numPrograms]; + if(programs) + { + programs[1].fParam1 = 0.33f; + programs[1].fParam5 = 0.75f; + programs[1].fParam6 = 0.57f; + strcpy(programs[1].name,"Slow"); + programs[2].fParam1 = 0.66f; + programs[2].fParam5 = 0.60f; + programs[2].fParam6 = 0.70f; + strcpy(programs[2].name,"Fast"); + setProgram(0); + } + + chp = dchp = clp = dclp = shp = dshp = slp = dslp = 0.0f; + + lspd = 0.0f; hspd = 0.0f; + lphi = 0.0f; hphi = 1.6f; + setParameter(0, 0.66f); +} + +bool mdaLeslie::getProductString(char* text) { strcpy(text, "mda Leslie"); return true; } +bool mdaLeslie::getVendorString(char* text) { strcpy(text, "mda"); return true; } +bool mdaLeslie::getEffectName(char* name) { strcpy(name, "Leslie"); return true; } + +void mdaLeslie::setParameter(LvzInt32 index, float value) +{ + + float ifs = 1.0f / getSampleRate(); + float spd = twopi * ifs * 2.0f * fParam8; + + switch(index) + { + case 0: programs[curProgram].fParam1 = fParam1 = value; break; + case 1: programs[curProgram].fParam7 = fParam7 = value; break; + case 2: programs[curProgram].fParam9 = fParam9 = value; break; + case 3: programs[curProgram].fParam4 = fParam4 = value; break; + case 4: programs[curProgram].fParam5 = fParam5 = value; break; + case 5: programs[curProgram].fParam6 = fParam6 = value; break; + case 6: programs[curProgram].fParam3 = fParam3 = value; break; + case 7: programs[curProgram].fParam2 = fParam2 = value; break; + case 8: programs[curProgram].fParam8 = fParam8 = value; break; + } + //calcs here! + filo = 1.f - (float)pow(10.0f, fParam3 * (2.27f - 0.54f * fParam3) - 1.92f); + + if(fParam1<0.50f) + { + if(fParam1<0.1f) //stop + { + lset = 0.00f; hset = 0.00f; + lmom = 0.12f; hmom = 0.10f; + } + else //low speed + { + lset = 0.49f; hset = 0.66f; + lmom = 0.27f; hmom = 0.18f; + } + } + else //high speed + { + lset = 5.31f; hset = 6.40f; + lmom = 0.14f; hmom = 0.09f; + } + hmom = (float)pow(10.0f, -ifs / hmom); + lmom = (float)pow(10.0f, -ifs / lmom); + hset *= spd; + lset *= spd; + + gain = 0.4f * (float)pow(10.0f, 2.0f * fParam2 - 1.0f); + lwid = fParam7 * fParam7; + llev = gain * 0.9f * fParam9 * fParam9; + hwid = fParam4 * fParam4; + hdep = fParam5 * fParam5 * getSampleRate() / 760.0f; + hlev = gain * 0.9f * fParam6 * fParam6; +} + +mdaLeslie::~mdaLeslie() +{ + if(hbuf) delete [] hbuf; + if(programs) delete [] programs; +} + +void mdaLeslie::setProgram(LvzInt32 program) +{ + mdaLeslieProgram *p = &programs[program]; + + curProgram = program; + setParameter(0, p->fParam1); + setParameter(1, p->fParam7); + setParameter(2, p->fParam9); + setParameter(3, p->fParam4); + setParameter(4, p->fParam5); + setParameter(5, p->fParam6); + setParameter(6, p->fParam3); + setParameter(7, p->fParam2); + setParameter(8, p->fParam8); + setProgramName(p->name); +} + +void mdaLeslie::suspend() +{ + memset(hbuf, 0, size * sizeof(float)); +} + +void mdaLeslie::setProgramName(char *name) +{ + strcpy(programName, name); +} + +void mdaLeslie::getProgramName(char *name) +{ + strcpy(name, programName); +} + +float mdaLeslie::getParameter(LvzInt32 index) +{ + float v=0; + + switch(index) + { + case 0: v = fParam1; break; + case 1: v = fParam7; break; + case 2: v = fParam9; break; + case 3: v = fParam4; break; + case 4: v = fParam5; break; + case 5: v = fParam6; break; + case 6: v = fParam3; break; + case 7: v = fParam2; break; + case 8: v = fParam8; break; + } + return v; +} + +void mdaLeslie::getParameterName(LvzInt32 index, char *label) +{ + switch(index) + { + case 0: strcpy(label, "Speed:"); break; + case 1: strcpy(label, "Lo Width"); break; + case 2: strcpy(label, "Lo Throb"); break; + case 3: strcpy(label, "Hi Width"); break; + case 4: strcpy(label, "Hi Depth"); break; + case 5: strcpy(label, "Hi Throb"); break; + case 6: strcpy(label, "X-Over"); break; + case 7: strcpy(label, "Output"); break; + case 8: strcpy(label, "Speed "); break; + } +} + +#include +void long2string(long value, char *string) { sprintf(string, "%ld", value); } + +void mdaLeslie::getParameterDisplay(LvzInt32 index, char *text) +{ + switch(index) + { + case 0: + if(fParam1<0.5f) + { + if(fParam1 < 0.1f) strcpy(text, "STOP"); + else strcpy(text, "SLOW"); + } else strcpy(text, "FAST"); break; + case 1: long2string((long)(100 * fParam7), text); break; + case 2: long2string((long)(100 * fParam9), text); break; + case 3: long2string((long)(100 * fParam4), text); break; + case 4: long2string((long)(100 * fParam5), text); break; + case 5: long2string((long)(100 * fParam6), text); break; + case 6: long2string((long)(10*int((float)pow(10.0f,1.179f + fParam3))), text); break; + case 7: long2string((long)(40 * fParam2 - 20), text); break; + case 8: long2string((long)(200 * fParam8), text); break; + } +} + +void mdaLeslie::getParameterLabel(LvzInt32 index, char *label) +{ + switch(index) + { + case 0: strcpy(label, ""); break; + case 6: strcpy(label, "Hz"); break; + case 7: strcpy(label, "dB"); break; + default: strcpy(label, "%"); break; + } +} + +//-------------------------------------------------------------------------------- + +void mdaLeslie::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, c, d, g=gain, h, l; + float fo=filo, fb1=fbuf1, fb2=fbuf2; + float hl=hlev, hs=hspd, ht, hm=hmom, hp=hphi, hw=hwid, hd=hdep; + float ll=llev, ls=lspd, lt, lm=lmom, lp=lphi, lw=lwid; + float hint, k0=0.03125f, k1=32.f; + long hdd, hdd2, k=0, hps=hpos; + + ht=hset*(1.f-hm); + lt=lset*(1.f-lm); + + chp = (float)cos(hp); chp *= chp * chp; + clp = (float)cos(lp); + shp = (float)sin(hp); + slp = (float)sin(lp); + + --in1; + --in2; + --out1; + --out2; + while(--sampleFrames >= 0) + { + a = *++in1 + *++in2; + c = out1[1]; + d = out2[1]; //see processReplacing() for comments + + if(k) k--; else + { + ls = (lm * ls) + lt; + hs = (hm * hs) + ht; + lp += k1 * ls; + hp += k1 * hs; + + dchp = (float)cos(hp + k1*hs); + dchp = k0 * (dchp * dchp * dchp - chp); + dclp = k0 * ((float)cos(lp + k1*ls) - clp); + dshp = k0 * ((float)sin(hp + k1*hs) - shp); + dslp = k0 * ((float)sin(lp + k1*ls) - slp); + + k=(long)k1; + } + + fb1 = fo * (fb1 - a) + a; + fb2 = fo * (fb2 - fb1) + fb1; + h = (g - hl * chp) * (a - fb2); + l = (g - ll * clp) * fb2; + + if(hps>0) hps--; else hps=200; + hint = hps + hd * (1.0f + chp); + hdd = (int)hint; + hint = hint - hdd; + hdd2 = hdd + 1; + if(hdd>199) { if(hdd>200) hdd -= 201; hdd2 -= 201; } + + *(hbuf + hps) = h; + a = *(hbuf + hdd); + h += a + hint * ( *(hbuf + hdd2) - a); + + c += l + h; + d += l + h; + h *= hw * shp; + l *= lw * slp; + d += l - h; + c += h - l; + + *++out1 = c; + *++out2 = d; + + chp += dchp; + clp += dclp; + shp += dshp; + slp += dslp; + } + lspd = ls; + hspd = hs; + hpos = hps; + lphi = (float)fmod(lp+(k1-k)*ls,twopi); + hphi = (float)fmod(hp+(k1-k)*hs,twopi); + if(fabs(fb1)>1.0e-10) fbuf1=fb1; else fbuf1=0.f; + if(fabs(fb2)>1.0e-10) fbuf2=fb2; else fbuf2=0.f; +} + +void mdaLeslie::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, c, d, g=gain, h, l; + float fo=filo, fb1=fbuf1, fb2=fbuf2; + float hl=hlev, hs=hspd, ht, hm=hmom, hp=hphi, hw=hwid, hd=hdep; + float ll=llev, ls=lspd, lt, lm=lmom, lp=lphi, lw=lwid; + float hint, k0=0.03125f, k1=32.f; //k0 = 1/k1 + long hdd, hdd2, k=0, hps=hpos; + + ht=hset*(1.f-hm); //target speeds + lt=lset*(1.f-lm); + + chp = (float)cos(hp); chp *= chp * chp; //set LFO values + clp = (float)cos(lp); + shp = (float)sin(hp); + slp = (float)sin(lp); + + --in1; + --in2; + --out1; + --out2; + while(--sampleFrames >= 0) + { + a = *++in1 + *++in2; //mono input + + if(k) k--; else //linear piecewise approx to LFO waveforms + { + ls = (lm * ls) + lt; //tend to required speed + hs = (hm * hs) + ht; + lp += k1 * ls; + hp += k1 * hs; + + dchp = (float)cos(hp + k1*hs); + dchp = k0 * (dchp * dchp * dchp - chp); //sin^3 level mod + dclp = k0 * ((float)cos(lp + k1*ls) - clp); + dshp = k0 * ((float)sin(hp + k1*hs) - shp); + dslp = k0 * ((float)sin(lp + k1*ls) - slp); + + k=(long)k1; + } + + fb1 = fo * (fb1 - a) + a; //crossover + fb2 = fo * (fb2 - fb1) + fb1; + h = (g - hl * chp) * (a - fb2); //volume + l = (g - ll * clp) * fb2; + + if(hps>0) hps--; else hps=200; //delay input pos + hint = hps + hd * (1.0f + chp); //delay output pos + hdd = (int)hint; + hint = hint - hdd; //linear intrpolation + hdd2 = hdd + 1; + if(hdd>199) { if(hdd>200) hdd -= 201; hdd2 -= 201; } + + *(hbuf + hps) = h; //delay input + a = *(hbuf + hdd); + h += a + hint * ( *(hbuf + hdd2) - a); //delay output + + c = l + h; + d = l + h; + h *= hw * shp; + l *= lw * slp; + d += l - h; + c += h - l; + + *++out1 = c; //output + *++out2 = d; + + chp += dchp; + clp += dclp; + shp += dshp; + slp += dslp; + } + lspd = ls; + hspd = hs; + hpos = hps; + lphi = (float)fmod(lp+(k1-k)*ls,twopi); + hphi = (float)fmod(hp+(k1-k)*hs,twopi); + if(fabs(fb1)>1.0e-10) fbuf1=fb1; else fbuf1=0.0f; //catch denormals + if(fabs(fb2)>1.0e-10) fbuf2=fb2; else fbuf2=0.0f; +} -- cgit v1.2.1