band-limited wavetable Sqr and Tri implementation finished for now in…

… W.purs

src/AudioWorklet.purs

@@ -60,14 +60,33 @@ function clamp(min,max,x) { return Math.max(Math.min(max,x),min); }
 function ain(input,n) { return (n >= input.length ? 0.0 : input[n]); }
 function genSaw() {
   var r = new Float32Array(4096).fill(0.5);
-  for(var k=1;k<=84;k++) { -- with 84 harmonics, highest harmonic of middle C is just below 22050 Hz
+  for(var k=1;k<=84;k++) { // with 84 harmonics, highest harmonic of middle C is just below 22050 Hz
     var x = Math.pow(-1,k);
     for(var t=0;t<4096;t++) {
       r[t] -= Math.sin(Math.PI * 2.0 * k * t / 4096) * x / (k * Math.PI);
     }
   }
   return r;
 }
+function genSqr() {
+  var r = new Float32Array(4096).fill(0.5);
+  for(var k=1;k<=83;k+=2) { // with 84 harmonics, highest harmonic of middle C is just below 22050 Hz
+    for(var t=0;t<4096;t++) {
+      r[t] += Math.sin(Math.PI * 2.0 * k * t / 4096) * 4 / (k * Math.PI);
+    }
+  }
+  return r;
+}
+function genTri() {
+  var r = new Float32Array(4096).fill(0.5);
+  for(var k=1;k<=83;k+=2) { // with 84 harmonics, highest harmonic of middle C is just below 22050 Hz
+    var x = Math.pow(-1,(k-1)*0.5);
+    for(var t=0;t<4096;t++) {
+      r[t] += Math.sin(Math.PI * 2.0 * k * t / 4096) * 8 * x / (k * k * Math.PI * Math.PI);
+    }
+  }
+  return r;
+}
 
 """
     classHeader = "class " <> name <> " extends AudioWorkletProcessor {\n\n"
@@ -81,7 +100,7 @@ return [
 ];}
 
 """
-    constructor = "constructor() { super(); this.saw = genSaw(); this.framesOut = 0; this.runTime = currentTime; this.f = new Float32Array(" <> show wState.allocatedFloats <> ").fill(0); this.i = new Int32Array(" <> show wState.allocatedInts <> ").fill(0);}\n\n"
+    constructor = "constructor() { super(); this.saw=genSaw(); this.sqr=genSqr(); this.tri=genTri(); this.framesOut=0; this.runTime=currentTime; this.f=new Float32Array(" <> show wState.allocatedFloats <> ").fill(0); this.i=new Int32Array(" <> show wState.allocatedInts <> ").fill(0);}\n\n"
     innerLoopPrefix = """process(inputs,outputs,parameters) {
 const input = inputs[0];
 const output = outputs[0];
@@ -101,6 +120,8 @@ const eTime = t - evalTimeAudio;
 const eBeat = eTime * cps;
 const fOut = fOutEnd == -1.0 ? 1.0 : clamp(0,1,(fOutEnd-t)/fOutDur);
 const saw = this.saw;
+const sqr = this.sqr;
+const tri = this.tri;
 """
     fadeCalculations = "const fIn = clamp(0,1,(t-" <> show fInStart <> ")/" <> show fInDur <> ");\nconst fade = Math.min(fIn,fOut);\n"
     outputIndices = range 0 (length frame - 1)

src/W.purs

@@ -209,6 +209,12 @@ wavetable tableName tableLength phase = do
 saw :: Sample -> W Sample
 saw f = phasor f >>= wavetable "saw" 4096
 
+sqr :: Sample -> W Sample
+sqr f = phasor f >>= wavetable "sqr" 4096
+
+tri :: Sample -> W Sample
+tri f = phasor f >>= wavetable "tri" 4096
+
 type Frame = Matrix Sample
 
 signalToFrame :: Signal -> W Frame
@@ -251,9 +257,9 @@ signalToFrame (Bipolar x) = signalToFrame x >>= traverse bipolar
 signalToFrame (Unipolar x) = signalToFrame x >>= traverse unipolar
 
 signalToFrame (Osc f) = signalToFrame f >>= traverse osc
--- Tri Signal
+signalToFrame (Tri f) = signalToFrame f >>= traverse tri
 signalToFrame (Saw f) = signalToFrame f >>= traverse saw
--- Sqr Signal
+signalToFrame (Sqr f) = signalToFrame f >>= traverse sqr
 
 signalToFrame (LFTri f) = do
   xs <- signalToFrame f >>= traverse phasor