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, 323 insertions, 0 deletions

diff --git a/src/mdaVocInput.cpp b/src/mdaVocInput.cpp
new file mode 100644
index 0000000..a595158
--- /dev/null
+++ b/src/mdaVocInput.cpp
@@ -0,0 +1,323 @@
+//
+// Plug-in: "mda VocSource" v1.0  (generates carrier signal for vocoder)
+//
+// Copyright(c)1999-2000 Paul Kellett (maxim digital audio)
+//
+
+#include "mdaVocInput.h"
+
+#include <stdlib.h>
+#include <stdio.h>
+#include <float.h>
+#include <math.h>
+
+AudioEffect *createEffectInstance(audioMasterCallback audioMaster)
+{
+  return new mdaVocInput(audioMaster);
+}
+
+mdaVocInputProgram::mdaVocInputProgram() ///default program settings
+{
+  param[0] = 0.25f;  //Tracking Off / On / Quant
+  param[1] = 0.50f;  //Pitch
+  param[2] = 0.20f;  //Breath Noise
+  param[3] = 0.50f;  //Voiced/Unvoiced Thresh
+  param[4] = 0.35f;  //Max Freq
+  strcpy(name, "Vocoder Carrier Signal");
+}
+
+
+mdaVocInput::mdaVocInput(audioMasterCallback audioMaster): AudioEffectX(audioMaster, NPROGS, NPARAMS)
+{
+  setNumInputs(2);
+  setNumOutputs(2);
+  setUniqueID("mda3");  ///identify plug-in here 
+	DECLARE_LVZ_DEPRECATED(canMono) ();				      
+  canProcessReplacing();
+
+  programs = new mdaVocInputProgram[numPrograms];
+  setProgram(0);
+
+  suspend();
+}
+
+bool  mdaVocInput::getProductString(char* text) { strcpy(text, "mda VocInput"); return true; }
+bool  mdaVocInput::getVendorString(char* text)  { strcpy(text, "mda"); return true; }
+bool  mdaVocInput::getEffectName(char* name)    { strcpy(name, "VocInput"); return true; }
+
+void mdaVocInput::resume() ///update internal parameters...
+{
+  float fs, ifs;
+  fs = getSampleRate();
+  ifs = 1.0f / fs;
+
+  track = (long)(2.99f * param[0]);
+  pmult = (float)pow(1.0594631f, floor(48.0f * param[1] - 24.0f));
+  if(track==0) pstep = 110.0f * pmult * ifs;
+
+  noise = 6.0f * param[2];
+  lfreq = 660.0f * ifs;
+  minp = (float)pow(16.0f, 0.5f - param[4]) * fs / 440.0f;
+  maxp = 0.03f * fs;
+  root = log10(8.1757989f * ifs);
+  vuv = param[3] * param[3];
+}
+
+
+void mdaVocInput::suspend() ///clear any buffers...
+{
+  lbuf0 = lbuf1 = lbuf2 = lbuf3 = 0.0f;
+  pstep = sawbuf = lenv = 0.0f;
+}
+
+
+mdaVocInput::~mdaVocInput() ///destroy any buffers...
+{
+  if(programs) delete [] programs;
+}
+
+
+void mdaVocInput::setProgram(LvzInt32 program)
+{
+  int i=0;
+
+  mdaVocInputProgram *p = &programs[program];
+  curProgram = program;
+  setProgramName(p->name);
+  for(i=0; i<NPARAMS; i++) param[i] = p->param[i];
+  resume();
+}
+
+
+void  mdaVocInput::setParameter(LvzInt32 index, float value) { param[index] = value; resume(); }
+float mdaVocInput::getParameter(LvzInt32 index) { return param[index]; }
+void  mdaVocInput::setProgramName(char *name) { strcpy(programName, name); }
+void  mdaVocInput::getProgramName(char *name) { strcpy(name, programName); }
+
+
+void mdaVocInput::getParameterName(LvzInt32 index, char *label)
+{
+  switch(index)
+  {
+    case  0: strcpy(label, "Tracking"); break;
+    case  1: strcpy(label, "Pitch"); break;
+    case  2: strcpy(label, "Breath"); break;
+    case  3: strcpy(label, "S Thresh"); break;
+    default: strcpy(label, "Max Freq");
+  }
+}
+
+
+void mdaVocInput::getParameterDisplay(LvzInt32 index, char *text)
+{
+ 	char string[16];
+
+  switch(index)
+  {
+    case  0: switch(track) { case 0: strcpy(string, "OFF");   break;
+                             case 1: strcpy(string, "FREE");   break;
+                             case 2: strcpy(string, "QUANT"); } break;
+
+    case  1: if(track) sprintf(string, "%ld", (long)(48.0f * param[1] - 24.0f)); 
+                  else midi2string((long)(48.0f * param[1] + 21.0f), string); break;
+
+    case  4: midi2string((long)(48.0f * param[4] + 45.0f), string); break;
+
+    default: sprintf(string, "%.0f", 100.0f * param[index]);
+  }
+	string[8] = 0;
+	strcpy(text, (char *)string);
+}
+
+
+void mdaVocInput::getParameterLabel(LvzInt32 index, char *label)
+{
+  switch(index)
+  {
+    case  0: strcpy(label, ""); break;
+    case  1: strcpy(label, ""); break;
+    case  4: strcpy(label, "Hz"); break;
+    default: strcpy(label, "%");
+  }
+}
+
+
+void mdaVocInput::midi2string(long n, char *text) //show name of MIDI note number (60=C3)
+{
+  char t[8];
+  long o, s, p=0;
+
+  t[p++] = ' ';
+  t[p++] = ' ';
+  t[p++] = ' ';
+
+  o = int(n/12.f); s = n-(12*o); o -= 2;
+
+  switch(s)
+  {
+    case  0: t[p++]='C';             break;
+    case  1: t[p++]='C'; t[p++]='#'; break;
+    case  2: t[p++]='D';             break;
+    case  3: t[p++]='D'; t[p++]='#'; break;
+    case  4: t[p++]='E';             break;
+    case  5: t[p++]='F';             break;
+    case  6: t[p++]='F'; t[p++]='#'; break;
+    case  7: t[p++]='G';             break;
+    case  8: t[p++]='G'; t[p++]='#'; break;
+    case  9: t[p++]='A';             break;
+    case 10: t[p++]='A'; t[p++]='#'; break;
+    default: t[p++]='B';             
+  }    
+  t[p++] = ' ';
+  
+  if(o < 0) { t[p++] = '-';  o = -o; }
+  t[p++] = (char)(48 + (o % 10)); 
+  
+  t[p] = 0; 
+  strcpy(text, t);
+}
+
+
+void mdaVocInput::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;
+  float ds=pstep, s=sawbuf, n=noise;
+  float l0=lbuf0, l1=lbuf1, l2=lbuf2, l3=lbuf3; 
+  float le=lenv, he=henv, et=lfreq*0.1f, lf=lfreq, v=vuv, mn=minp, mx=maxp;
+  float rootm=39.863137f;
+  long  tr=track;
+
+  --in1;
+  --in2;
+  --out1;
+  --out2;
+  while(--sampleFrames >= 0)
+  {
+    a = *++in1; 
+    b = *++in2;
+    c = out1[1];
+    d = out2[1];
+
+    l0 -= lf * (l1 + a);       //fundamental filter (peaking 2nd-order 100Hz lpf)
+    l1 -= lf * (l1 - l0);
+
+    b = l0; if(b<0.0f) b = -b;
+    le -= et * (le - b);       //fundamental level
+
+    b = (a + 0.03f) * v;
+    if(b<0.0f) b = -b;
+    he -= et * (he - b);       //overall level (+ constant so >f0 when quiet)
+
+    l3 += 1.0f;
+    if(tr>0)                   //pitch tracking
+    {
+      if(l1>0.0f && l2<=0.0f)  //found +ve zero crossing
+      {
+        if(l3>mn && l3<mx)     //...in allowed range
+        {
+          mn = 0.6f * l3;       //new max pitch to discourage octave jumps!
+          l2 = l1 / (l1 - l2);   //fractional period...
+          ds = pmult / (l3 - l2); //new period
+
+          if(tr==2)            //quantize pitch
+          {
+            ds = rootm * (float)(log10(ds) - root);
+            ds = (float)pow(1.0594631, floor(ds + 0.5) + rootm * root);   
+          }
+        }
+        l3 = l2;               //restart period measurement
+      }
+      l2=l1;                   //remember previous sample
+    }
+  
+    b = 0.00001f * (float)((rand() & 32767) - 16384);  //sibilance
+    if(le>he) b *= s * n;                    //...or modulated breath noise
+    b += s; s += ds; if(s>0.5f) s-=1.0f;     //badly aliased sawtooth!
+
+    c += a;
+    d += b;
+
+    *++out1 = c;
+    *++out2 = d;
+  }
+  sawbuf=s;
+
+  if(fabs(he)>1.0e-10) henv = he; else henv=0.0f; //catch denormals
+  if(fabs(l1)>1.0e-10) { lbuf0=l0; lbuf1=l1; lenv=le; } else { lbuf0 = lbuf1= lenv = 0.0f; }
+  lbuf2=l2, lbuf3=l3;
+  if(tr) pstep=ds;     
+}
+
+
+void mdaVocInput::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;
+  float ds=pstep, s=sawbuf, n=noise;
+  float l0=lbuf0, l1=lbuf1, l2=lbuf2, l3=lbuf3; 
+  float le=lenv, he=henv, et=lfreq*0.1f, lf=lfreq, v=vuv, mn=minp, mx=maxp;
+  float rootm=39.863137f;
+  long  tr=track;
+
+  --in1;
+  --in2;
+  --out1;
+  --out2;
+  while(--sampleFrames >= 0)
+  {
+    a = *++in1; 
+    b = *++in2;
+
+    l0 -= lf * (l1 + a);       //fundamental filter (peaking 2nd-order 100Hz lpf)
+    l1 -= lf * (l1 - l0);
+
+    b = l0; if(b<0.0f) b = -b;
+    le -= et * (le - b);       //fundamental level
+
+    b = (a + 0.03f) * v;
+    if(b<0.0f) b = -b;
+    he -= et * (he - b);       //overall level (+ constant so >f0 when quiet)
+
+    l3 += 1.0f;
+    if(tr>0)                   //pitch tracking
+    {
+      if(l1>0.0f && l2<=0.0f)  //found +ve zero crossing
+      {
+        if(l3>mn && l3<mx)     //...in allowed range
+        {
+          mn = 0.6f * l3;       //new max pitch to discourage octave jumps!
+          l2 = l1 / (l1 - l2);   //fractional period...
+          ds = pmult / (l3 - l2); //new period
+
+          if(tr==2)            //quantize pitch
+          {
+            ds = rootm * (float)(log10(ds) - root);
+            ds = (float)pow(1.0594631, floor(ds + 0.5) + rootm * root);   
+          }
+        }
+        l3 = l2;               //restart period measurement
+      }
+      l2=l1;                   //remember previous sample
+    }
+  
+    b = 0.00001f * (float)((rand() & 32767) - 16384);  //sibilance
+    if(le>he) b *= s * n;                    //...or modulated breath noise
+    b += s; s += ds; if(s>0.5f) s-=1.0f;     //badly aliased sawtooth!
+
+    *++out1 = a;
+    *++out2 = b;
+  }
+  sawbuf=s;
+
+  if(fabs(he)>1.0e-10) henv = he; else henv=0.0f; //catch denormals
+  if(fabs(l1)>1.0e-10) { lbuf0=l0; lbuf1=l1; lenv=le; } else { lbuf0 = lbuf1= lenv = 0.0f; }
+  lbuf2=l2, lbuf3=l3;
+  if(tr) pstep=ds; 
+}