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

1 files changed, 267 insertions, 0 deletions

diff --git a/src/mdaDegrade.cpp b/src/mdaDegrade.cpp
new file mode 100644
index 0000000..7ad1a82
--- /dev/null
+++ b/src/mdaDegrade.cpp
@@ -0,0 +1,267 @@
+#include "mdaDegrade.h"
+
+#include <math.h>
+#include <float.h>
+
+AudioEffect *createEffectInstance(audioMasterCallback audioMaster)
+{
+  return new mdaDegrade(audioMaster);
+}
+
+mdaDegrade::mdaDegrade(audioMasterCallback audioMaster)	: AudioEffectX(audioMaster, 1, 6)	// programs, parameters
+{
+  //inits here!
+  fParam1 = (float)0.8; //clip
+  fParam2 = (float)0.50; //bits
+  fParam3 = (float)0.65; //rate
+  fParam4 = (float)0.9; //postfilt
+  fParam5 = (float)0.58; //non-lin
+  fParam6 = (float)0.5; //level
+
+  setNumInputs(2);		  
+	setNumOutputs(2);		  
+	setUniqueID("mdaC");    // identify here
+	DECLARE_LVZ_DEPRECATED(canMono) ();				      
+	canProcessReplacing();	
+	strcpy(programName, "Degrade");
+

+  buf0 = buf1 = buf2 = buf3 = buf4 = buf5 = buf6 = buf7 = buf8 = buf9 = 0.0f;

+
+  setParameter(5, 0.5f);
+}
+
+bool  mdaDegrade::getProductString(char* text) { strcpy(text, "mda Degrade"); return true; }
+bool  mdaDegrade::getVendorString(char* text)  { strcpy(text, "mda"); return true; }
+bool  mdaDegrade::getEffectName(char* name)    { strcpy(name, "Degrade"); return true; }
+
+void mdaDegrade::setParameter(LvzInt32 index, float value)
+{
+	float f;
+  
+  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;
+  }
+  //calcs here
+  if(fParam3>0.5) { f = fParam3 - 0.5f;  mode = 1.0f; }
+             else { f = 0.5f - fParam3;  mode = 0.0f; }
+    
+  tn = (int)exp(18.0f * f);
+    
+    //tn = (int)(18.0 * fParam3 - 8.0); mode=1.f; }
+    //         else { tn = (int)(10.0 - 18.0 * fParam3); mode=0.f; }
+  tcount = 1;
+  clp = (float)(pow(10.0,(-1.5 + 1.5 * fParam1)) );
+  fo2 = filterFreq((float)pow(10.0f, 2.30104f + 2.f*fParam4));
+  fi2 = (1.f-fo2); fi2=fi2*fi2; fi2=fi2*fi2;
+  float _g1 = (float)(pow(2.0,2.0 + int(fParam2*12.0))); 
+  g2 = (float)(1.0/(2.0 * _g1));
+  if(fParam3>0.5) g1 = -_g1/(float)tn; else g1= -_g1; 
+  g3 = (float)(pow(10.0,2.0*fParam6 - 1.0));
+  if(fParam5>0.5) { lin = (float)(pow(10.0,0.3 * (0.5 - fParam5))); lin2=lin; }
+  else { lin = (float)pow(10.0,0.3 * (fParam5 - 0.5)); lin2=1.0; }
+}
+
+mdaDegrade::~mdaDegrade()
+{
+}
+
+void mdaDegrade::suspend()
+{
+}
+
+float mdaDegrade::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 mdaDegrade::setProgramName(char *name)
+{
+	strcpy(programName, name);
+}
+
+void mdaDegrade::getProgramName(char *name)
+{
+	strcpy(name, programName);
+}
+
+float mdaDegrade::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;
+  }
+  return v;
+}
+
+void mdaDegrade::getParameterName(LvzInt32 index, char *label)
+{
+	switch(index)
+  {
+    case 0: strcpy(label, "Headroom"); break;
+    case 1: strcpy(label, "Quant"); break;
+    case 2: strcpy(label, "Rate"); break;
+    case 3: strcpy(label, "PostFilt"); break;
+    case 4: strcpy(label, "Non-Lin"); break;
+    case 5: strcpy(label, "Output"); break;
+  }
+}
+
+#include <stdio.h>
+void long2string(long value, char *string) { sprintf(string, "%ld", value); }
+
+void mdaDegrade::getParameterDisplay(LvzInt32 index, char *text)
+{
+	switch(index)
+  {
+    case 0: long2string((long)(-30.0 * (1.0 - fParam1)), text); break;
+    case 1: long2string((long)(4.0 + 12.0 * fParam2), text); break;
+    case 2: long2string((long)(getSampleRate()/tn), text); break;
+    case 3: long2string((long)pow(10.0f, 2.30104f + 2.f*fParam4), text); break;
+    case 4: long2string((long)(200.0 * fabs(fParam5 - 0.5)), text); break;
+    case 5: long2string((long)(40.0 * fParam6 - 20.0), text); break;
+  }
+}
+
+void mdaDegrade::getParameterLabel(LvzInt32 index, char *label)
+{
+	switch(index)
+  {
+    case 0: strcpy(label, "dB"); break;
+    case 1: strcpy(label, "bits"); break;
+    case 2: strcpy(label, "S<>S&&H"); break;
+    case 3: strcpy(label, "Hz"); break;
+    case 4: strcpy(label, "Odd<>Eve"); break;
+  }
+}
+
+//--------------------------------------------------------------------------------
+// process
+
+void mdaDegrade::process(float **inputs, float **outputs, LvzInt32 sampleFrames)
+{
+	float *in1 = inputs[0];
+	float *in2 = inputs[1];
+	float *out1 = outputs[0];
+	float *out2 = outputs[1];
+	float b0=buf0, c, d, l=lin, l2=lin2;
+  float cl=clp, i2=fi2, o2=fo2;
+  float b1=buf1, b2=buf2, b3=buf3, b4=buf4, b5=buf5;
+  float b6=buf6, b7=buf7, b8=buf8, b9=buf9;
+  float gi=g1, go=g2, ga=g3, m=mode;  
+  int n=tn, t=tcount;   
+  
+  --in1;	
+	--in2;	
+	--out1;
+	--out2;
+	while(--sampleFrames >= 0)
+	{
+    c = out1[1];
+		d = out2[1]; //process from here...
+          
+    b0 = (*++in1 + *++in2) + m * b0;
+
+    if(t==n)
+    {
+      t=0;
+      b5 = (float)(go * int(b0 * gi));
+      if(b5>0)
+      { b5=(float)pow(b5,l2); if(b5>cl) b5=cl;}
+      else
+      { b5=-(float)pow(-b5,l); if(b5<-cl) b5=-cl; }
+      b0=0;
+    } 
+    t=t+1;
+
+    b1 = i2 * (b5 * ga) + o2 * b1;
+    b2 =      b1 + o2 * b2;
+    b3 =      b2 + o2 * b3;
+    b4 =      b3 + o2 * b4;
+    b6 = i2 * b4 + o2 * b6;
+    b7 =      b6 + o2 * b7;
+    b8 =      b7 + o2 * b8;
+    b9 =      b8 + o2 * b9;
+
+    c += b9; // output
+		d += b9;
+
+    *++out1 = c;	
+		*++out2 = d;
+	}
+  if(fabs(b1)<1.0e-10) 
+  { buf1=0.f; buf2=0.f; buf3=0.f; buf4=0.f; buf6=0.f; buf7=0.f;
+    buf8=0.f; buf9=0.f; buf0=0.f; buf5=0.f; }
+  else 
+  { buf1=b1; buf2=b2; buf3=b3; buf4=b4; buf6=b6; buf7=b7; 
+    buf8=b8, buf9=b9; buf0=b0; buf5=b5; tcount=t; }
+}
+
+void mdaDegrade::processReplacing(float **inputs, float **outputs, LvzInt32 sampleFrames)
+{
+	float *in1 = inputs[0];
+	float *in2 = inputs[1];
+	float *out1 = outputs[0];
+	float *out2 = outputs[1];
+	float b0=buf0, l=lin, l2=lin2;
+  float cl=clp, i2=fi2, o2=fo2;
+  float b1=buf1, b2=buf2, b3=buf3, b4=buf4, b5=buf5;
+  float b6=buf6, b7=buf7, b8=buf8, b9=buf9;
+  float gi=g1, go=g2, ga=g3, m=mode;  
+  int n=tn, t=tcount;   
+  
+	--in1;	
+	--in2;	
+	--out1;
+	--out2;
+	while(--sampleFrames >= 0)
+	{
+    b0 = (*++in1 + *++in2) + m * b0;
+
+    if(t==n)
+    {
+      t=0;
+      b5 = (float)(go * int(b0 * gi));
+      if(b5>0)
+      { b5=(float)pow(b5,l2); if(b5>cl) b5=cl;}
+      else
+      { b5=-(float)pow(-b5,l); if(b5<-cl) b5=-cl; }
+      b0=0;
+    } 
+    t=t+1;
+
+    b1 = i2 * (b5 * ga) + o2 * b1;
+    b2 =      b1 + o2 * b2;
+    b3 =      b2 + o2 * b3;
+    b4 =      b3 + o2 * b4;
+    b6 = i2 * b4 + o2 * b6;
+    b7 =      b6 + o2 * b7;
+    b8 =      b7 + o2 * b8;
+    b9 =      b8 + o2 * b9;
+
+		*++out1 = b9;
+		*++out2 = b9;
+	}
+  if(fabs(b1)<1.0e-10) 
+  { buf1=0.f; buf2=0.f; buf3=0.f; buf4=0.f; buf6=0.f; buf7=0.f;
+    buf8=0.f; buf9=0.f; buf0=0.f; buf5=0.f; }
+  else 
+  { buf1=b1; buf2=b2; buf3=b3; buf4=b4; buf6=b6; buf7=b7; 
+    buf8=b8, buf9=b9; buf0=b0; buf5=b5; tcount=t; }
+}