@vessel-dsp/runtime 0.4.1 → 0.6.0

package/dist/worklet/v2-audio-worklet.js

@@ -585,37 +585,38 @@ var require_v2_dsp = __commonJS({
         var _v2_testonly_resample_reset = Module["_v2_testonly_resample_reset"] = (a0) => (_v2_testonly_resample_reset = Module["_v2_testonly_resample_reset"] = wasmExports["q"])(a0);
         var _v2_testonly_resample_process = Module["_v2_testonly_resample_process"] = (a0, a1) => (_v2_testonly_resample_process = Module["_v2_testonly_resample_process"] = wasmExports["r"])(a0, a1);
         var _v2_testonly_resample_destroy = Module["_v2_testonly_resample_destroy"] = (a0) => (_v2_testonly_resample_destroy = Module["_v2_testonly_resample_destroy"] = wasmExports["s"])(a0);
-        var _v2_engine_reset = Module["_v2_engine_reset"] = (a0) => (_v2_engine_reset = Module["_v2_engine_reset"] = wasmExports["t"])(a0);
-        var _v2_engine_set_control = Module["_v2_engine_set_control"] = (a0, a1, a2) => (_v2_engine_set_control = Module["_v2_engine_set_control"] = wasmExports["u"])(a0, a1, a2);
-        var _v2_engine_get_control = Module["_v2_engine_get_control"] = (a0, a1) => (_v2_engine_get_control = Module["_v2_engine_get_control"] = wasmExports["v"])(a0, a1);
-        var _v2_engine_set_supply = Module["_v2_engine_set_supply"] = (a0, a1, a2, a3, a4) => (_v2_engine_set_supply = Module["_v2_engine_set_supply"] = wasmExports["w"])(a0, a1, a2, a3, a4);
-        var _v2_engine_process_sample = Module["_v2_engine_process_sample"] = (a0, a1) => (_v2_engine_process_sample = Module["_v2_engine_process_sample"] = wasmExports["x"])(a0, a1);
-        var _v2_engine_process_block = Module["_v2_engine_process_block"] = (a0, a1, a2, a3) => (_v2_engine_process_block = Module["_v2_engine_process_block"] = wasmExports["y"])(a0, a1, a2, a3);
-        var _v2_engine_get_input_buffer = Module["_v2_engine_get_input_buffer"] = (a0, a1) => (_v2_engine_get_input_buffer = Module["_v2_engine_get_input_buffer"] = wasmExports["z"])(a0, a1);
-        var _v2_engine_get_output_buffer = Module["_v2_engine_get_output_buffer"] = (a0, a1) => (_v2_engine_get_output_buffer = Module["_v2_engine_get_output_buffer"] = wasmExports["A"])(a0, a1);
-        var _v2_engine_process_internal = Module["_v2_engine_process_internal"] = (a0, a1) => (_v2_engine_process_internal = Module["_v2_engine_process_internal"] = wasmExports["B"])(a0, a1);
-        var _v2_engine_get_last_iteration_count = Module["_v2_engine_get_last_iteration_count"] = (a0) => (_v2_engine_get_last_iteration_count = Module["_v2_engine_get_last_iteration_count"] = wasmExports["C"])(a0);
-        var _v2_engine_get_last_converged = Module["_v2_engine_get_last_converged"] = (a0) => (_v2_engine_get_last_converged = Module["_v2_engine_get_last_converged"] = wasmExports["D"])(a0);
-        var _v2_engine_get_max_iterations = Module["_v2_engine_get_max_iterations"] = (a0) => (_v2_engine_get_max_iterations = Module["_v2_engine_get_max_iterations"] = wasmExports["E"])(a0);
-        var _v2_engine_get_operating_point = Module["_v2_engine_get_operating_point"] = (a0, a1, a2) => (_v2_engine_get_operating_point = Module["_v2_engine_get_operating_point"] = wasmExports["F"])(a0, a1, a2);
-        var _v2_engine_get_state = Module["_v2_engine_get_state"] = (a0, a1, a2) => (_v2_engine_get_state = Module["_v2_engine_get_state"] = wasmExports["G"])(a0, a1, a2);
-        var _v2_engine_get_last_error = Module["_v2_engine_get_last_error"] = (a0) => (_v2_engine_get_last_error = Module["_v2_engine_get_last_error"] = wasmExports["H"])(a0);
-        var _v2_engine_get_schedule_solves = Module["_v2_engine_get_schedule_solves"] = (a0) => (_v2_engine_get_schedule_solves = Module["_v2_engine_get_schedule_solves"] = wasmExports["I"])(a0);
-        var _v2_engine_get_schedule_fallbacks = Module["_v2_engine_get_schedule_fallbacks"] = (a0) => (_v2_engine_get_schedule_fallbacks = Module["_v2_engine_get_schedule_fallbacks"] = wasmExports["J"])(a0);
-        var _v2_engine_get_predictor_seeds = Module["_v2_engine_get_predictor_seeds"] = (a0) => (_v2_engine_get_predictor_seeds = Module["_v2_engine_get_predictor_seeds"] = wasmExports["K"])(a0);
-        var _v2_engine_get_one_iteration_solves = Module["_v2_engine_get_one_iteration_solves"] = (a0) => (_v2_engine_get_one_iteration_solves = Module["_v2_engine_get_one_iteration_solves"] = wasmExports["L"])(a0);
-        var _v2_engine_get_total_iterations = Module["_v2_engine_get_total_iterations"] = (a0) => (_v2_engine_get_total_iterations = Module["_v2_engine_get_total_iterations"] = wasmExports["M"])(a0);
-        var _v2_engine_get_predictor_order = Module["_v2_engine_get_predictor_order"] = (a0, a1) => (_v2_engine_get_predictor_order = Module["_v2_engine_get_predictor_order"] = wasmExports["N"])(a0, a1);
-        var _v2_engine_get_kernel_solves = Module["_v2_engine_get_kernel_solves"] = (a0) => (_v2_engine_get_kernel_solves = Module["_v2_engine_get_kernel_solves"] = wasmExports["O"])(a0);
-        var _v2_engine_get_repivoted_schedule_blocks = Module["_v2_engine_get_repivoted_schedule_blocks"] = (a0) => (_v2_engine_get_repivoted_schedule_blocks = Module["_v2_engine_get_repivoted_schedule_blocks"] = wasmExports["P"])(a0);
-        var _v2_engine_get_abandoned_schedule_blocks = Module["_v2_engine_get_abandoned_schedule_blocks"] = (a0) => (_v2_engine_get_abandoned_schedule_blocks = Module["_v2_engine_get_abandoned_schedule_blocks"] = wasmExports["Q"])(a0);
-        var _v2_engine_get_dropped_schedule_blocks = Module["_v2_engine_get_dropped_schedule_blocks"] = (a0) => (_v2_engine_get_dropped_schedule_blocks = Module["_v2_engine_get_dropped_schedule_blocks"] = wasmExports["R"])(a0);
-        var _malloc = Module["_malloc"] = (a0) => (_malloc = Module["_malloc"] = wasmExports["S"])(a0);
-        var _free = Module["_free"] = (a0) => (_free = Module["_free"] = wasmExports["T"])(a0);
-        var ___trap = () => (___trap = wasmExports["U"])();
-        var __emscripten_stack_restore = (a0) => (__emscripten_stack_restore = wasmExports["V"])(a0);
-        var __emscripten_stack_alloc = (a0) => (__emscripten_stack_alloc = wasmExports["W"])(a0);
-        var _emscripten_stack_get_current = () => (_emscripten_stack_get_current = wasmExports["X"])();
+        var _v2_testonly_taper_fraction = Module["_v2_testonly_taper_fraction"] = (a0, a1) => (_v2_testonly_taper_fraction = Module["_v2_testonly_taper_fraction"] = wasmExports["t"])(a0, a1);
+        var _v2_engine_reset = Module["_v2_engine_reset"] = (a0) => (_v2_engine_reset = Module["_v2_engine_reset"] = wasmExports["u"])(a0);
+        var _v2_engine_set_control = Module["_v2_engine_set_control"] = (a0, a1, a2) => (_v2_engine_set_control = Module["_v2_engine_set_control"] = wasmExports["v"])(a0, a1, a2);
+        var _v2_engine_get_control = Module["_v2_engine_get_control"] = (a0, a1) => (_v2_engine_get_control = Module["_v2_engine_get_control"] = wasmExports["w"])(a0, a1);
+        var _v2_engine_set_supply = Module["_v2_engine_set_supply"] = (a0, a1, a2, a3, a4) => (_v2_engine_set_supply = Module["_v2_engine_set_supply"] = wasmExports["x"])(a0, a1, a2, a3, a4);
+        var _v2_engine_process_sample = Module["_v2_engine_process_sample"] = (a0, a1) => (_v2_engine_process_sample = Module["_v2_engine_process_sample"] = wasmExports["y"])(a0, a1);
+        var _v2_engine_process_block = Module["_v2_engine_process_block"] = (a0, a1, a2, a3) => (_v2_engine_process_block = Module["_v2_engine_process_block"] = wasmExports["z"])(a0, a1, a2, a3);
+        var _v2_engine_get_input_buffer = Module["_v2_engine_get_input_buffer"] = (a0, a1) => (_v2_engine_get_input_buffer = Module["_v2_engine_get_input_buffer"] = wasmExports["A"])(a0, a1);
+        var _v2_engine_get_output_buffer = Module["_v2_engine_get_output_buffer"] = (a0, a1) => (_v2_engine_get_output_buffer = Module["_v2_engine_get_output_buffer"] = wasmExports["B"])(a0, a1);
+        var _v2_engine_process_internal = Module["_v2_engine_process_internal"] = (a0, a1) => (_v2_engine_process_internal = Module["_v2_engine_process_internal"] = wasmExports["C"])(a0, a1);
+        var _v2_engine_get_last_iteration_count = Module["_v2_engine_get_last_iteration_count"] = (a0) => (_v2_engine_get_last_iteration_count = Module["_v2_engine_get_last_iteration_count"] = wasmExports["D"])(a0);
+        var _v2_engine_get_last_converged = Module["_v2_engine_get_last_converged"] = (a0) => (_v2_engine_get_last_converged = Module["_v2_engine_get_last_converged"] = wasmExports["E"])(a0);
+        var _v2_engine_get_max_iterations = Module["_v2_engine_get_max_iterations"] = (a0) => (_v2_engine_get_max_iterations = Module["_v2_engine_get_max_iterations"] = wasmExports["F"])(a0);
+        var _v2_engine_get_operating_point = Module["_v2_engine_get_operating_point"] = (a0, a1, a2) => (_v2_engine_get_operating_point = Module["_v2_engine_get_operating_point"] = wasmExports["G"])(a0, a1, a2);
+        var _v2_engine_get_state = Module["_v2_engine_get_state"] = (a0, a1, a2) => (_v2_engine_get_state = Module["_v2_engine_get_state"] = wasmExports["H"])(a0, a1, a2);
+        var _v2_engine_get_last_error = Module["_v2_engine_get_last_error"] = (a0) => (_v2_engine_get_last_error = Module["_v2_engine_get_last_error"] = wasmExports["I"])(a0);
+        var _v2_engine_get_schedule_solves = Module["_v2_engine_get_schedule_solves"] = (a0) => (_v2_engine_get_schedule_solves = Module["_v2_engine_get_schedule_solves"] = wasmExports["J"])(a0);
+        var _v2_engine_get_schedule_fallbacks = Module["_v2_engine_get_schedule_fallbacks"] = (a0) => (_v2_engine_get_schedule_fallbacks = Module["_v2_engine_get_schedule_fallbacks"] = wasmExports["K"])(a0);
+        var _v2_engine_get_predictor_seeds = Module["_v2_engine_get_predictor_seeds"] = (a0) => (_v2_engine_get_predictor_seeds = Module["_v2_engine_get_predictor_seeds"] = wasmExports["L"])(a0);
+        var _v2_engine_get_one_iteration_solves = Module["_v2_engine_get_one_iteration_solves"] = (a0) => (_v2_engine_get_one_iteration_solves = Module["_v2_engine_get_one_iteration_solves"] = wasmExports["M"])(a0);
+        var _v2_engine_get_total_iterations = Module["_v2_engine_get_total_iterations"] = (a0) => (_v2_engine_get_total_iterations = Module["_v2_engine_get_total_iterations"] = wasmExports["N"])(a0);
+        var _v2_engine_get_predictor_order = Module["_v2_engine_get_predictor_order"] = (a0, a1) => (_v2_engine_get_predictor_order = Module["_v2_engine_get_predictor_order"] = wasmExports["O"])(a0, a1);
+        var _v2_engine_get_kernel_solves = Module["_v2_engine_get_kernel_solves"] = (a0) => (_v2_engine_get_kernel_solves = Module["_v2_engine_get_kernel_solves"] = wasmExports["P"])(a0);
+        var _v2_engine_get_repivoted_schedule_blocks = Module["_v2_engine_get_repivoted_schedule_blocks"] = (a0) => (_v2_engine_get_repivoted_schedule_blocks = Module["_v2_engine_get_repivoted_schedule_blocks"] = wasmExports["Q"])(a0);
+        var _v2_engine_get_abandoned_schedule_blocks = Module["_v2_engine_get_abandoned_schedule_blocks"] = (a0) => (_v2_engine_get_abandoned_schedule_blocks = Module["_v2_engine_get_abandoned_schedule_blocks"] = wasmExports["R"])(a0);
+        var _v2_engine_get_dropped_schedule_blocks = Module["_v2_engine_get_dropped_schedule_blocks"] = (a0) => (_v2_engine_get_dropped_schedule_blocks = Module["_v2_engine_get_dropped_schedule_blocks"] = wasmExports["S"])(a0);
+        var _malloc = Module["_malloc"] = (a0) => (_malloc = Module["_malloc"] = wasmExports["T"])(a0);
+        var _free = Module["_free"] = (a0) => (_free = Module["_free"] = wasmExports["U"])(a0);
+        var ___trap = () => (___trap = wasmExports["V"])();
+        var __emscripten_stack_restore = (a0) => (__emscripten_stack_restore = wasmExports["W"])(a0);
+        var __emscripten_stack_alloc = (a0) => (__emscripten_stack_alloc = wasmExports["X"])(a0);
+        var _emscripten_stack_get_current = () => (_emscripten_stack_get_current = wasmExports["Y"])();
         Module["ccall"] = ccall;
         Module["cwrap"] = cwrap;
         Module["setValue"] = setValue;
@@ -1167,6 +1168,58 @@ function cascadeLatencyHostSamples(stages) {
 }
 
 // src/taper.ts
+var IMPLEMENTED_TAPERS = {
+  linear: true,
+  logarithmic: true,
+  "reverse-logarithmic": true,
+  "reverse-linear": true,
+  "s-curve": true
+};
+var IMPLEMENTED_TAPER_NAMES = Object.keys(IMPLEMENTED_TAPERS);
+function unrecognisedTapers(program) {
+  const found = [];
+  const check = (where, taper) => {
+    if (taper === void 0) return;
+    if (typeof taper !== "string" || !Object.hasOwn(IMPLEMENTED_TAPERS, taper)) {
+      found.push({ where, taper: String(taper) });
+    }
+  };
+  for (const control of program.controls) {
+    check(`control "${control.id}"`, control.taper);
+  }
+  for (const block of program.blocks) {
+    if (block.kind === "mna") {
+      for (const stamp of block.stamps) {
+        if ("taper" in stamp) {
+          check(`${stamp.kind} stamp in block "${block.id}"`, stamp.taper);
+        }
+      }
+    } else if (block.clockControl != null) {
+      check(`clockControl of block "${block.id}"`, block.clockControl.taper);
+    }
+  }
+  return found;
+}
+var S_CURVE_SEGMENTS = [
+  { x0: 0, h: 0.3, y0: 0, y1: 0.1, m0: 0, m1: 0.6 },
+  { x0: 0.3, h: 0.2, y0: 0.1, y1: 0.5, m0: 0.6, m1: 2 },
+  { x0: 0.5, h: 0.2, y0: 0.5, y1: 0.9, m0: 2, m1: 0.6 },
+  { x0: 0.7, h: 0.3, y0: 0.9, y1: 1, m0: 0.6, m1: 0 }
+];
+function sCurveFraction(x) {
+  if (x <= 0) return 0;
+  if (x >= 1) return 1;
+  const segment = x < 0.3 ? 0 : x < 0.5 ? 1 : x < 0.7 ? 2 : 3;
+  const { x0, h, y0, y1, m0, m1 } = S_CURVE_SEGMENTS[segment];
+  const t = (x - x0) / h;
+  const t2 = t * t;
+  const t3 = t2 * t;
+  const h00 = 2 * t3 - 3 * t2 + 1;
+  const h10 = t3 - 2 * t2 + t;
+  const h01 = -2 * t3 + 3 * t2;
+  const h11 = t3 - t2;
+  return h00 * y0 + h10 * h * m0 + h01 * y1 + h11 * h * m1;
+}
 function taperFraction(taper, position) {
   const x = Math.min(1, Math.max(0, position));
   const audioBase = 1089 / 49;
@@ -1179,6 +1232,9 @@ function taperFraction(taper, position) {
   if (taper === "reverse-linear") {
     return 1 - x;
   }
+  if (taper === "s-curve") {
+    return sCurveFraction(x);
+  }
   return x;
 }
 
@@ -1339,7 +1395,10 @@ var OPAMP_FOLD_STREAK = 6;
 function opAmpHalfSwing(railHigh, railLow) {
   return Math.max((railHigh - railLow) / 2, 1e-9);
 }
-var OPAMP_GAIN_BANDWIDTH_HZ = 1e6;
+var OPAMP_DEFAULT_GAIN_BANDWIDTH_HZ = 1e6;
+function opAmpGainBandwidthHz(stamp) {
+  return stamp.gainBandwidthHz ?? OPAMP_DEFAULT_GAIN_BANDWIDTH_HZ;
+}
 var OPAMP_POLE_BAND = 8;
 function boundOpAmpRaw(raw, halfSwing) {
   const limit = OPAMP_POLE_BAND * halfSwing;
@@ -1850,6 +1909,8 @@ var ReferenceRuntime = class {
   limiterTrace = typeof process !== "undefined" && process.env?.VESSEL_LIMITER_TRACE === "1";
   limitedIterationCount = 0;
   limitedIterationPerSample = [];
+  /** Fold reseeds taken so far (diagnostic, read past `private` by the stall study's tracer). */
+  foldReseedCount = 0;
   /** Series ohms of whatever drives the input jack -- see `prepare()`'s option. */
   inputSourceOhms = 0;
   /**
@@ -2162,6 +2223,12 @@ var ReferenceRuntime = class {
         `program requires ${missingModels.length === 1 ? "a DSP model" : "DSP models"} this runtime does not implement: ${missingModels.join(", ")}`
       );
     }
+    const unrecognised = unrecognisedTapers(this.program);
+    if (unrecognised.length > 0) {
+      throw new RuntimeError(
+        `program names ${unrecognised.length === 1 ? "a taper" : "tapers"} this runtime does not implement: ${unrecognised.map((entry) => `"${entry.taper}" (${entry.where})`).join(", ")}; it implements ${IMPLEMENTED_TAPER_NAMES.join(", ")}`
+      );
+    }
     const dispatched = this.program.blocks.filter(
       (block) => block.kind === "macro"
     );
@@ -4025,7 +4092,9 @@ var ReferenceRuntime = class {
     let bestDelta = Number.POSITIVE_INFINITY;
     let noImprovement = 0;
     let foldStreak = 0;
-    let foldReseeded = false;
+    let foldReseeds = 0;
+    let lastWorstNode = -1;
+    let lastWorstSign = 0;
     for (let iteration = 0; iteration < iterations; iteration += 1) {
       used = iteration + 1;
       this.limitedIterate = false;
@@ -4153,6 +4222,12 @@ var ReferenceRuntime = class {
           }
         }
       }
+      const damping = dc ? 1 : this.opAmpStepDamping(block, blockIndex, current, next, dt, dc);
+      if (damping < 1) {
+        for (let i = 0; i < size; i++) {
+          next[i] = current[i] + damping * (next[i] - current[i]);
+        }
+      }
       if (relaxing) {
         for (let i = 0; i < size; i++) {
           next[i] = current[i] + alpha * (next[i] - current[i]);
@@ -4161,6 +4236,10 @@ var ReferenceRuntime = class {
       const delta = maxAbsDifference(next, current);
       worstNode = worstDifferenceIndex(next, current);
       worstDelta = delta;
+      const worstSign = Math.sign((next[worstNode] ?? 0) - (current[worstNode] ?? 0));
+      const reversed = dc || worstNode === lastWorstNode && worstSign !== 0 && lastWorstSign !== 0 && worstSign !== lastWorstSign;
+      lastWorstNode = worstNode;
+      lastWorstSign = worstSign;
       if (delta < bestDelta * 0.999) {
         bestDelta = delta;
         noImprovement = 0;
@@ -4173,7 +4252,7 @@ var ReferenceRuntime = class {
           }
         }
       } else {
-        noImprovement += 1;
+        noImprovement += reversed ? 1 : 0;
         if (!relaxing) {
           if (noImprovement >= NEWTON_NON_CONTRACTING_LIMIT && iteration >= NEWTON_RELAXATION_EARLIEST_ITERATION) {
             relaxing = true;
@@ -4204,13 +4283,15 @@ var ReferenceRuntime = class {
       }
       const limitedOpamp = this.limitedOpamp;
       foldStreak = limitedOpamp?.folded === true ? foldStreak + 1 : 0;
-      if (!dc && !foldReseeded && limitedOpamp !== null && foldStreak >= OPAMP_FOLD_STREAK) {
-        foldReseeded = true;
+      if (!dc && limitedOpamp !== null && foldStreak >= OPAMP_FOLD_STREAK) {
+        foldReseeds += 1;
+        this.foldReseedCount += 1;
         foldStreak = 0;
         const wasHigh = (start[limitedOpamp.output] ?? limitedOpamp.centre) >= limitedOpamp.centre;
-        current[limitedOpamp.output] = wasHigh ? limitedOpamp.railLow : limitedOpamp.railHigh;
+        const toHigh = foldReseeds % 2 === 1 ? wasHigh : !wasHigh;
+        current[limitedOpamp.output] = toHigh ? limitedOpamp.railHigh : limitedOpamp.railLow;
         this.opampHistory.set(limitedOpamp.key, {
-          differential: wasHigh ? -limitedOpamp.band : limitedOpamp.band,
+          differential: toHigh ? limitedOpamp.band : -limitedOpamp.band,
           step: 0,
           cap: limitedOpamp.maxStep
         });
@@ -4250,6 +4331,53 @@ var ReferenceRuntime = class {
     }
     return { solution: current, converged, used, worstNode, worstDelta };
   }
+  /**
+   * The factor the Newton update is scaled by so that no railed op-amp's differential moves
+   * further than its step cap and no op-amp output node moves further than one full swing:
+   * the smallest `cap / |requested differential step|` and `2 * halfSwing / |requested output
+   * step|` over the block's railed op-amps, or 1. `current` is the iterate the rows were
+   * stamped at and `next` the undamped solve. The history read here is the one the stamps
+   * wrote this iteration, so `differential` is the point each row's tangent was taken at and
+   * `cap` the step the next evaluation will allow; a request inside both bounds scales nothing.
+   * Mirrored in `Engine.cpp` (`opAmpStepDamping`); the two must agree to the operation.
+   */
+  opAmpStepDamping(block, blockIndex, current, next, dt, dc) {
+    let damping = 1;
+    for (const stamp of block.stamps) {
+      if (stamp.kind !== "ideal-opamp" || stamp.railHigh === null || stamp.railLow === null) {
+        continue;
+      }
+      const history = this.opampHistory.get(
+        packHistoryKey3(blockIndex, stamp.plus, stamp.minus, stamp.output)
+      );
+      if (history === void 0) continue;
+      const halfSwing = opAmpHalfSwing(stamp.railHigh, stamp.railLow);
+      const alpha = dc ? 1 : opAmpPoleAlpha(
+        stamp.openLoopGain,
+        opAmpGainBandwidthHz(stamp),
+        dt,
+        this.sampleRate
+      );
+      const effectiveGain = alpha * stamp.openLoopGain;
+      const linearWidth = halfSwing / effectiveGain;
+      const band = OPAMP_DIFFERENTIAL_BAND * linearWidth;
+      const nextDifferential = clamp(
+        (next[stamp.plus] ?? 0) - (next[stamp.minus] ?? 0),
+        -band,
+        band
+      );
+      const magnitude = Math.abs(nextDifferential - history.differential);
+      if (magnitude > history.cap) {
+        damping = Math.min(damping, history.cap / magnitude);
+      }
+      const swing = 2 * halfSwing;
+      const outputStep = Math.abs((next[stamp.output] ?? 0) - (current[stamp.output] ?? 0));
+      if (outputStep > swing) {
+        damping = Math.min(damping, swing / outputStep);
+      }
+    }
+    return damping;
+  }
   /**
    * Build-order step 6, vertical slice: the part of a block's matrix that is the same for
    * every Newton iteration of this sample -- every stamp `newtonPortRowsOf` calls linear,
@@ -5073,7 +5201,7 @@ var ReferenceRuntime = class {
       const key = packHistoryKey3(blockIndex, stamp.plus, stamp.minus, stamp.output);
       const alpha = opAmpPoleAlpha(
         stamp.openLoopGain,
-        OPAMP_GAIN_BANDWIDTH_HZ,
+        opAmpGainBandwidthHz(stamp),
         dt,
         this.sampleRate
       );
@@ -5718,26 +5846,29 @@ var ReferenceRuntime = class {
         );
         const alpha = dc ? 1 : opAmpPoleAlpha(
           stamp.openLoopGain,
-          OPAMP_GAIN_BANDWIDTH_HZ,
+          opAmpGainBandwidthHz(stamp),
           dt,
           this.sampleRate
         );
         const effectiveGain = alpha * stamp.openLoopGain;
         const linearWidth = halfSwing / effectiveGain;
         const previous = this.opampHistory.get(opampKey);
-        const previousDifferential = previous?.differential ?? 0;
         const rawDifferential = (solution[stamp.plus] ?? 0) - (solution[stamp.minus] ?? 0);
         const band = OPAMP_DIFFERENTIAL_BAND * linearWidth;
         const bounded = clamp(rawDifferential, -band, band);
+        const previousDifferential = previous?.differential ?? bounded;
         const step = bounded - previousDifferential;
         const maxStep = OPAMP_MAX_DIFFERENTIAL_STEP * linearWidth;
         const reversing = previous !== void 0 && previous.step !== 0 && step !== 0 && Math.sign(step) !== Math.sign(previous.step);
         const previousCap = previous?.cap ?? maxStep;
         const cap = reversing ? Math.max(previousCap / 2, linearWidth / 64) : Math.min(maxStep, previousCap * 1.2);
         const limited = Math.abs(step) > cap;
+        const folded = cap <= linearWidth / 8 && Math.abs(step) >= linearWidth / 256;
         if (limited) {
           this.limitedIterate = true;
           this.limitedBy = `opamp +=${stamp.plus} -=${stamp.minus} out=${stamp.output}`;
+        }
+        if (limited || folded) {
           this.limitedOpamp = {
             key: opampKey,
             output: stamp.output,
@@ -5746,7 +5877,7 @@ var ReferenceRuntime = class {
             railLow: stamp.railLow,
             band,
             maxStep,
-            folded: cap <= linearWidth / 8
+            folded
           };
         }
         const differential = limited ? previousDifferential + Math.sign(step) * cap : bounded;
@@ -5763,10 +5894,11 @@ var ReferenceRuntime = class {
           clamp(scaled, -OPAMP_SATURATION_BAND, OPAMP_SATURATION_BAND)
         );
         const slope = effectiveGain * (1 - derivativeShape * derivativeShape);
+        const linearisedAt = limited ? differential : rawDifferential;
         matrix[row][stamp.output] += 1;
         matrix[row][stamp.plus] -= slope;
         matrix[row][stamp.minus] += slope;
-        rhs[row] = output - slope * differential;
+        rhs[row] = output - slope * linearisedAt;
         break;
       }
       case "vccs": {

package/package.json

@@ -1,6 +1,6 @@
 {
 	"name": "@vessel-dsp/runtime",
-	"version": "0.4.1",
+	"version": "0.6.0",
 	"description": "Real-time MNA simulation runtime and solver console for VesselDSP Program ROMs.",
 	"keywords": [
 		"guitar-pedal",
@@ -58,7 +58,7 @@
 		"generate-kernels": "bun scripts/generate-kernels.ts"
 	},
 	"dependencies": {
-		"@vessel-dsp/compiler": "0.4.0"
+		"@vessel-dsp/compiler": "0.6.0"
 	},
 	"publishConfig": {
 		"access": "public"