@vcad/mcp 0.9.4-main.11 → 0.9.4-main.13

package/README.md

@@ -6,4 +6,4 @@ vcad's MCP server as a self-contained bundle (server + BRep kernel WASM).
 { "mcpServers": { "vcad": { "command": "npx", "args": ["-y", "@vcad/mcp"] } } }
 ```
 
-Built from 681857ca2f62f173237de79ff2621cd1609480a7 at 2026-07-24T16:25:54.297Z. Source: https://github.com/ecto/vcad
+Built from 70552964e3760fb5b103235917b00364f6d45763 at 2026-07-24T16:29:07.170Z. Source: https://github.com/ecto/vcad

package/index.mjs

@@ -8022,14 +8022,22 @@ async function loadAndInit() {
     }
     await mod.default(hint ? { module_or_path: hint } : void 0);
   }
+  assertKernelBaseline(mod);
   wasmModule2 = mod;
   generation++;
   return mod;
 }
+function assertKernelBaseline(mod) {
+  const missing = REQUIRED_BINDINGS.filter((name) => typeof mod[name] !== "function");
+  if (missing.length > 0) {
+    const version2 = typeof mod.get_kernel_version === "function" ? mod.get_kernel_version() : "unknown";
+    throw new Error(`stale @vcad/kernel-wasm bundle (kernel version ${version2}): missing required bindings ${missing.join(", ")} \u2014 rebuild the workspace (VCAD_WASM_SKIP unset) or update to a current bundle`);
+  }
+}
 function getKernelWasmSync() {
   return wasmModule2;
 }
-var wasmPromise, wasmModule2, wasmInputHint, bindgenStarted, lastModuleBuffer, lastTrap, generation;
+var wasmPromise, wasmModule2, wasmInputHint, bindgenStarted, lastModuleBuffer, lastTrap, generation, REQUIRED_BINDINGS;
 var init_wasm_singleton = __esm({
   "../engine/dist/wasm-singleton.js"() {
     "use strict";
@@ -8038,6 +8046,13 @@ var init_wasm_singleton = __esm({
     bindgenStarted = false;
     lastTrap = null;
     generation = 0;
+    REQUIRED_BINDINGS = [
+      "evaluateDocument",
+      "documentDiff",
+      "documentMerge",
+      "documentToLoon",
+      "documentToLoonChecked"
+    ];
   }
 });
 
@@ -8509,7 +8524,10 @@ function evaluateDocument2(doc, kernel2, options = {}) {
   const resolved = resolveDocument(doc);
   doc = resolved.doc;
   doc = expandPartInstances(doc, kernel2);
-  if (kernel2.evaluateDocument) {
+  if (!kernel2.evaluateDocument) {
+    throw new Error("kernel WASM bundle is missing the evaluateDocument binding \u2014 rebuild @vcad/kernel-wasm (stale bundle); old bundles are unsupported");
+  }
+  {
     try {
       const docJson = JSON.stringify(doc);
       const result = kernel2.evaluateDocument(docJson, options.skipClashDetection ?? false);
@@ -9637,114 +9655,11 @@ var init_gpu = __esm({
 });
 
 // ../engine/dist/diff.js
-function decompose(doc) {
-  const obj = doc;
-  const entities = /* @__PURE__ */ new Map();
-  for (const [field, kind] of MAP_TABLES) {
-    const table = obj[field];
-    if (table && typeof table === "object" && !Array.isArray(table)) {
-      for (const [id, v] of Object.entries(table)) {
-        entities.set(`${kind}\0${id}`, v);
-      }
-    }
-  }
-  for (const [field, kind, idField] of VEC_TABLES) {
-    const items = obj[field];
-    if (Array.isArray(items)) {
-      for (const item of items) {
-        const raw = item?.[idField];
-        const id = typeof raw === "string" ? raw : JSON.stringify(raw);
-        entities.set(`${kind}\0${id}`, item);
-      }
-    }
-  }
-  for (const field of DOC_FIELDS) {
-    const v = obj[field];
-    if (v !== void 0 && v !== null)
-      entities.set(`doc\0${field}`, v);
-  }
-  return entities;
-}
-function deepEqual(a, b) {
-  if (a === b)
-    return true;
-  if (typeof a !== typeof b || a === null || b === null)
-    return false;
-  if (typeof a !== "object")
-    return false;
-  if (Array.isArray(a) !== Array.isArray(b))
-    return false;
-  if (Array.isArray(a) && Array.isArray(b)) {
-    if (a.length !== b.length)
-      return false;
-    return a.every((v, i) => deepEqual(v, b[i]));
-  }
-  const ka = Object.keys(a).filter((k) => a[k] !== void 0);
-  const kb = Object.keys(b).filter((k) => b[k] !== void 0);
-  if (ka.length !== kb.length)
-    return false;
-  return ka.every((k) => deepEqual(a[k], b[k]));
-}
-function valueDiff(prefix, oldV, newV, out) {
-  if (deepEqual(oldV, newV))
-    return;
-  const join4 = (k) => prefix ? `${prefix}.${k}` : k;
-  const isObj = (v) => typeof v === "object" && v !== null && !Array.isArray(v);
-  if (isObj(oldV) && isObj(newV)) {
-    const keys = /* @__PURE__ */ new Set([...Object.keys(oldV), ...Object.keys(newV)]);
-    for (const k of keys) {
-      valueDiff(join4(k), oldV[k] ?? null, newV[k] ?? null, out);
-    }
-    return;
-  }
-  if (Array.isArray(oldV) && Array.isArray(newV) && oldV.length === newV.length) {
-    oldV.forEach((v, i) => valueDiff(join4(String(i)), v, newV[i], out));
-    return;
-  }
-  out.push({ path: prefix, old: oldV, new: newV });
-}
-function entityName(v) {
-  const n = v?.name;
-  return typeof n === "string" ? n : void 0;
-}
-function semanticDiffFallback(oldDoc, newDoc) {
-  const a = decompose(oldDoc);
-  const b = decompose(newDoc);
-  const keys = [.../* @__PURE__ */ new Set([...a.keys(), ...b.keys()])].sort();
-  const changes = [];
-  for (const key of keys) {
-    const [kind, id] = key.split("\0");
-    const va = a.get(key);
-    const vb = b.get(key);
-    if (va === void 0 && vb !== void 0) {
-      changes.push({ kind, id, name: entityName(vb), type: "added", value: vb });
-    } else if (va !== void 0 && vb === void 0) {
-      changes.push({ kind, id, name: entityName(va), type: "removed", value: va });
-    } else if (!deepEqual(va, vb)) {
-      const fields = [];
-      valueDiff("", va, vb, fields);
-      if (fields.length > 0) {
-        changes.push({
-          kind,
-          id,
-          name: entityName(vb) ?? entityName(va),
-          type: "modified",
-          fields
-        });
-      }
-    }
-  }
-  return { changes };
-}
 function semanticDiff(oldDoc, newDoc, kernel2) {
-  if (kernel2?.documentDiff) {
-    try {
-      return kernel2.documentDiff(JSON.stringify(oldDoc), JSON.stringify(newDoc));
-    } catch (e) {
-      console.warn("[diff] kernel documentDiff failed, using TS fallback:", e);
-    }
+  if (!kernel2?.documentDiff) {
+    throw new Error("semantic diff needs a kernel build with the documentDiff binding \u2014 ensure the kernel WASM is initialized (getKernelWasm) and current");
   }
-  return semanticDiffFallback(oldDoc, newDoc);
+  return kernel2.documentDiff(JSON.stringify(oldDoc), JSON.stringify(newDoc));
 }
 function mergeAvailable(kernel2) {
   return typeof kernel2?.documentMerge === "function";
@@ -9755,33 +9670,9 @@ function threeWayMerge(base, ours, theirs, resolutions, kernel2) {
   }
   return kernel2.documentMerge(JSON.stringify(base), JSON.stringify(ours), JSON.stringify(theirs), resolutions.length > 0 ? JSON.stringify(resolutions) : null);
 }
-var MAP_TABLES, VEC_TABLES, DOC_FIELDS;
 var init_diff = __esm({
   "../engine/dist/diff.js"() {
     "use strict";
-    MAP_TABLES = [
-      ["nodes", "node"],
-      ["materials", "material"],
-      ["part_materials", "part-material"],
-      ["partDefs", "part-def"],
-      ["parameters", "parameter"],
-      ["bindings", "binding"]
-    ];
-    VEC_TABLES = [
-      ["roots", "root", "root"],
-      ["instances", "instance", "id"],
-      ["joints", "joint", "id"],
-      ["clearance_specs", "clearance", "label"]
-    ];
-    DOC_FIELDS = [
-      "version",
-      "scene",
-      "groundInstanceId",
-      "schematic",
-      "pcb",
-      "molecule",
-      "timeline"
-    ];
   }
 });
 
@@ -11631,7 +11522,6 @@ __export(dist_exports, {
   searchEcadParts: () => searchEcadParts,
   searchParts: () => searchParts,
   semanticDiff: () => semanticDiff,
-  semanticDiffFallback: () => semanticDiffFallback,
   solveDesignConstraints: () => solveDesignConstraints2,
   solveForwardKinematics: () => solveForwardKinematics2,
   threeWayMerge: () => threeWayMerge,
@@ -11801,6 +11691,8 @@ var init_dist = __esm({
           sheetMetalSequence: wasmModule7.sheetMetalSequence,
           nestSheetMetalParts: wasmModule7.nestSheetMetalParts,
           getSheetMetalMaterials: wasmModule7.getSheetMetalMaterials,
+          simulateStrikeKernel: wasmModule7.simulateStrikeKernel,
+          noteToHz: wasmModule7.noteToHz,
           getSheetMetalBendTable: wasmModule7.getSheetMetalBendTable,
           getSheetMetalShopCatalog: wasmModule7.getSheetMetalShopCatalog,
           sheetMetalFoldedStep: wasmModule7.sheetMetalFoldedStep,
@@ -12342,6 +12234,27 @@ var init_dist = __esm({
           return null;
         return sheetMetalSequence2(chain, this.kernel);
       }
+      /**
+       * Run the kernel's mallet-strike pipeline (free-free bar modal analysis +
+       * synthesis + FFT verdict). Input/output shapes are JSON mirrors of
+       * `vcad_kernel_acoustics::strike::{StrikeInput, StrikeResult}` with the
+       * WAV base64-encoded as `wav_base64`.
+       */
+      simulateStrike(input) {
+        const kernel2 = this.kernel;
+        if (!kernel2.simulateStrikeKernel) {
+          throw new Error("kernel WASM build lacks simulateStrikeKernel \u2014 rebuild @vcad/kernel-wasm");
+        }
+        return JSON.parse(kernel2.simulateStrikeKernel(JSON.stringify(input)));
+      }
+      /** Parse a note name ("C6", "F#4") to Hz via the kernel. */
+      noteToHz(note) {
+        const kernel2 = this.kernel;
+        if (!kernel2.noteToHz) {
+          throw new Error("kernel WASM build lacks noteToHz \u2014 rebuild @vcad/kernel-wasm");
+        }
+        return kernel2.noteToHz(note);
+      }
       /** Return the kernel's curated sheet-metal materials registry. */
       getSheetMetalMaterials() {
         return getSheetMetalMaterials2(this.kernel);
@@ -26615,33 +26528,33 @@ var init_server = __esm({
 });
 
 // node_modules/@modelcontextprotocol/sdk/dist/esm/experimental/tasks/helpers.js
-function assertToolsCallTaskCapability(requests, method, entityName2) {
+function assertToolsCallTaskCapability(requests, method, entityName) {
   if (!requests) {
-    throw new Error(`${entityName2} does not support task creation (required for ${method})`);
+    throw new Error(`${entityName} does not support task creation (required for ${method})`);
   }
   switch (method) {
     case "tools/call":
       if (!requests.tools?.call) {
-        throw new Error(`${entityName2} does not support task creation for tools/call (required for ${method})`);
+        throw new Error(`${entityName} does not support task creation for tools/call (required for ${method})`);
       }
       break;
     default:
       break;
   }
 }
-function assertClientRequestTaskCapability(requests, method, entityName2) {
+function assertClientRequestTaskCapability(requests, method, entityName) {
   if (!requests) {
-    throw new Error(`${entityName2} does not support task creation (required for ${method})`);
+    throw new Error(`${entityName} does not support task creation (required for ${method})`);
   }
   switch (method) {
     case "sampling/createMessage":
       if (!requests.sampling?.createMessage) {
-        throw new Error(`${entityName2} does not support task creation for sampling/createMessage (required for ${method})`);
+        throw new Error(`${entityName} does not support task creation for sampling/createMessage (required for ${method})`);
       }
       break;
     case "elicitation/create":
       if (!requests.elicitation?.create) {
-        throw new Error(`${entityName2} does not support task creation for elicitation/create (required for ${method})`);
+        throw new Error(`${entityName} does not support task creation for elicitation/create (required for ${method})`);
       }
       break;
     default:
@@ -40930,6 +40843,7 @@ var init_save_load = __esm({
 var init_document_to_loon = __esm({
   "../core/dist/utils/document-to-loon.js"() {
     "use strict";
+    init_wasm_singleton2();
   }
 });
 
@@ -72670,335 +72584,6 @@ var init_sheet_metal2 = __esm({
 
 // src/tools/acoustics.ts
 import { writeFileSync as writeFileSync5 } from "node:fs";
-function freeFreeBetaL(count) {
-  const roots = [];
-  for (let n = 1; n <= count; n++) {
-    let x = (2 * n + 1) * Math.PI / 2;
-    for (let i = 0; i < 50; i++) {
-      const f = Math.cosh(x) * Math.cos(x) - 1;
-      const df = Math.sinh(x) * Math.cos(x) - Math.cosh(x) * Math.sin(x);
-      const step = f / df;
-      x -= step;
-      if (Math.abs(step) < 1e-13) break;
-    }
-    roots.push(x);
-  }
-  return roots;
-}
-function modeShape(betaL, xi) {
-  const sigma = (Math.cosh(betaL) - Math.cos(betaL)) / (Math.sinh(betaL) - Math.sin(betaL));
-  const raw = (t2) => Math.cosh(betaL * t2) + Math.cos(betaL * t2) - sigma * (Math.sinh(betaL * t2) + Math.sin(betaL * t2));
-  return raw(xi) / Math.abs(raw(0));
-}
-function closedFormHz(bar, count) {
-  const L = bar.length_mm / 1e3;
-  const t2 = bar.thickness_mm / 1e3;
-  const E = bar.modulus_gpa * 1e9;
-  const c = t2 * Math.sqrt(E / (12 * bar.density_kg_m3));
-  return freeFreeBetaL(count).map((bl) => bl * bl * c / (2 * Math.PI * L * L));
-}
-function effectiveWidthMm(bar, xMm) {
-  let w = bar.width_mm;
-  const r = bar.hole_dia_mm / 2;
-  for (const h of bar.holes_mm) {
-    const d = Math.abs(xMm - h);
-    if (d < r) w -= 2 * Math.sqrt(r * r - d * d);
-  }
-  return Math.max(w, 1e-6);
-}
-function femHz(bar, count, nel = 96) {
-  const L = bar.length_mm / 1e3;
-  const t2 = bar.thickness_mm / 1e3;
-  const E = bar.modulus_gpa * 1e9;
-  const rho = bar.density_kg_m3;
-  const nDof = 2 * (nel + 1);
-  const K = zeros(nDof, nDof);
-  const M = zeros(nDof, nDof);
-  const le = L / nel;
-  const gp = [-0.8611363115940526, -0.3399810435848563, 0.3399810435848563, 0.8611363115940526];
-  const gw = [0.3478548451374538, 0.6521451548625461, 0.6521451548625461, 0.3478548451374538];
-  for (let e = 0; e < nel; e++) {
-    const x0 = e * le;
-    const ke = zeros(4, 4);
-    const me = zeros(4, 4);
-    for (let g = 0; g < 4; g++) {
-      const xi = gp[g];
-      const x = x0 + (xi + 1) / 2 * le;
-      const wEff = effectiveWidthMm(bar, x * 1e3) / 1e3;
-      const A = wEff * t2;
-      const I = wEff * t2 * t2 * t2 / 12;
-      const s = le / 2;
-      const N = [
-        0.25 * (1 - xi) * (1 - xi) * (2 + xi),
-        s * 0.25 * (1 - xi) * (1 - xi) * (1 + xi),
-        0.25 * (1 + xi) * (1 + xi) * (2 - xi),
-        s * 0.25 * (1 + xi) * (1 + xi) * (xi - 1)
-      ];
-      const c2 = 2 / le * (2 / le);
-      const B = [
-        1.5 * xi * c2,
-        s * (1.5 * xi - 0.5) * c2,
-        -1.5 * xi * c2,
-        s * (1.5 * xi + 0.5) * c2
-      ];
-      const wJ = gw[g] * s;
-      for (let i = 0; i < 4; i++) {
-        for (let j = 0; j < 4; j++) {
-          ke[i][j] += E * I * B[i] * B[j] * wJ;
-          me[i][j] += rho * A * N[i] * N[j] * wJ;
-        }
-      }
-    }
-    const dof = [2 * e, 2 * e + 1, 2 * e + 2, 2 * e + 3];
-    for (let i = 0; i < 4; i++) {
-      for (let j = 0; j < 4; j++) {
-        K[dof[i]][dof[j]] += ke[i][j];
-        M[dof[i]][dof[j]] += me[i][j];
-      }
-    }
-  }
-  const Lc = cholesky(M);
-  const C = congruence(K, Lc);
-  const lambda = jacobiEigenvalues(C);
-  lambda.sort((a, b) => a - b);
-  const elastic = lambda.filter((l) => l > (2 * Math.PI) ** 2);
-  return elastic.slice(0, count).map((l) => Math.sqrt(l) / (2 * Math.PI));
-}
-function zeros(r, c) {
-  return Array.from({ length: r }, () => new Array(c).fill(0));
-}
-function cholesky(A) {
-  const n = A.length;
-  const L = zeros(n, n);
-  for (let i = 0; i < n; i++) {
-    for (let j = 0; j <= i; j++) {
-      let s = A[i][j];
-      for (let k = 0; k < j; k++) s -= L[i][k] * L[j][k];
-      if (i === j) {
-        if (s <= 0) throw new Error("mass matrix not positive definite");
-        L[i][i] = Math.sqrt(s);
-      } else {
-        L[i][j] = s / L[j][j];
-      }
-    }
-  }
-  return L;
-}
-function congruence(K, L) {
-  const n = K.length;
-  const Y = zeros(n, n);
-  for (let c = 0; c < n; c++) {
-    for (let i = 0; i < n; i++) {
-      let s = K[i][c];
-      for (let k = 0; k < i; k++) s -= L[i][k] * Y[k][c];
-      Y[i][c] = s / L[i][i];
-    }
-  }
-  const C = zeros(n, n);
-  for (let r = 0; r < n; r++) {
-    for (let i = 0; i < n; i++) {
-      let s = Y[r][i];
-      for (let k = 0; k < i; k++) s -= L[i][k] * C[r][k];
-      C[r][i] = s / L[i][i];
-    }
-  }
-  for (let i = 0; i < n; i++)
-    for (let j = 0; j < i; j++) {
-      const v = 0.5 * (C[i][j] + C[j][i]);
-      C[i][j] = v;
-      C[j][i] = v;
-    }
-  return C;
-}
-function jacobiEigenvalues(A) {
-  const n = A.length;
-  const a = A.map((row) => row.slice());
-  for (let sweep = 0; sweep < 30; sweep++) {
-    let off = 0;
-    for (let i = 0; i < n; i++)
-      for (let j = i + 1; j < n; j++) off += a[i][j] * a[i][j];
-    if (off < 1e-18 * n * n) break;
-    for (let p = 0; p < n - 1; p++) {
-      for (let q = p + 1; q < n; q++) {
-        if (Math.abs(a[p][q]) < 1e-30) continue;
-        const theta = (a[q][q] - a[p][p]) / (2 * a[p][q]);
-        const t2 = Math.sign(theta) / (Math.abs(theta) + Math.sqrt(theta * theta + 1));
-        const c = 1 / Math.sqrt(t2 * t2 + 1);
-        const s = t2 * c;
-        for (let k = 0; k < n; k++) {
-          const akp = a[k][p];
-          const akq = a[k][q];
-          a[k][p] = c * akp - s * akq;
-          a[k][q] = s * akp + c * akq;
-        }
-        for (let k = 0; k < n; k++) {
-          const apk = a[p][k];
-          const aqk = a[q][k];
-          a[p][k] = c * apk - s * aqk;
-          a[q][k] = s * apk + c * aqk;
-        }
-      }
-    }
-  }
-  return a.map((row, i) => row[i]);
-}
-function malletSpectrum(hz, contactMs) {
-  const tc = contactMs / 1e3;
-  const u = 2 * hz * tc;
-  const denom = Math.abs(1 - u * u);
-  if (denom < 1e-9) return Math.PI / 4;
-  return Math.abs(Math.cos(Math.PI * hz * tc)) / denom;
-}
-function strikeModes(bar, femFreqs, strikeFrac, contactMs, materialQ = 2500, suspensionQ0 = 150) {
-  const betas = freeFreeBetaL(femFreqs.length);
-  const modes = femFreqs.map((hz, i) => {
-    const phiStrike = modeShape(betas[i], strikeFrac);
-    const gain = Math.abs(phiStrike) * malletSpectrum(hz, contactMs);
-    let phiSq = 0;
-    for (const h of bar.holes_mm) {
-      const phi = modeShape(betas[i], h / bar.length_mm);
-      phiSq += phi * phi;
-    }
-    const qSusp = suspensionQ0 / Math.max(phiSq, 1e-9);
-    const q = 1 / (1 / materialQ + 1 / qSusp);
-    return {
-      n: i + 1,
-      hz,
-      gain,
-      q,
-      t60_s: Math.log(1e3) * q / (Math.PI * hz)
-    };
-  });
-  const peak = Math.max(...modes.map((m) => m.gain), 1e-12);
-  for (const m of modes) m.gain /= peak;
-  return modes;
-}
-function synthesize(modes, durationS, sampleRate) {
-  const n = Math.floor(durationS * sampleRate);
-  const out = new Float64Array(n);
-  for (const m of modes) {
-    if (m.hz >= 0.45 * sampleRate) continue;
-    const w = 2 * Math.PI * m.hz / sampleRate;
-    const tau = m.q / (Math.PI * m.hz);
-    const decay = Math.exp(-1 / (tau * sampleRate));
-    let amp = m.gain;
-    for (let k = 0; k < n; k++) {
-      out[k] += amp * Math.sin(w * k);
-      amp *= decay;
-    }
-  }
-  let peak = 1e-12;
-  for (let k = 0; k < n; k++) peak = Math.max(peak, Math.abs(out[k]));
-  const norm = 0.891 / peak;
-  for (let k = 0; k < n; k++) out[k] *= norm;
-  return out;
-}
-function encodeWav(samples, sampleRate) {
-  const n = samples.length;
-  const buf = new ArrayBuffer(44 + n * 2);
-  const v = new DataView(buf);
-  const str2 = (off, s) => {
-    for (let i = 0; i < s.length; i++) v.setUint8(off + i, s.charCodeAt(i));
-  };
-  str2(0, "RIFF");
-  v.setUint32(4, 36 + n * 2, true);
-  str2(8, "WAVE");
-  str2(12, "fmt ");
-  v.setUint32(16, 16, true);
-  v.setUint16(20, 1, true);
-  v.setUint16(22, 1, true);
-  v.setUint32(24, sampleRate, true);
-  v.setUint32(28, sampleRate * 2, true);
-  v.setUint16(32, 2, true);
-  v.setUint16(34, 16, true);
-  str2(36, "data");
-  v.setUint32(40, n * 2, true);
-  for (let i = 0; i < n; i++) {
-    const s = Math.max(-1, Math.min(1, samples[i]));
-    v.setInt16(44 + i * 2, Math.round(s * 32767), true);
-  }
-  return new Uint8Array(buf);
-}
-function fft(re, im) {
-  const n = re.length;
-  for (let i = 1, j = 0; i < n; i++) {
-    let bit = n >> 1;
-    for (; j & bit; bit >>= 1) j ^= bit;
-    j ^= bit;
-    if (i < j) {
-      [re[i], re[j]] = [re[j], re[i]];
-      [im[i], im[j]] = [im[j], im[i]];
-    }
-  }
-  for (let len = 2; len <= n; len <<= 1) {
-    const ang = -2 * Math.PI / len;
-    const wr = Math.cos(ang);
-    const wi = Math.sin(ang);
-    for (let i = 0; i < n; i += len) {
-      let cr = 1;
-      let ci = 0;
-      for (let k = 0; k < len / 2; k++) {
-        const ur = re[i + k];
-        const ui = im[i + k];
-        const vr = re[i + k + len / 2] * cr - im[i + k + len / 2] * ci;
-        const vi = re[i + k + len / 2] * ci + im[i + k + len / 2] * cr;
-        re[i + k] = ur + vr;
-        im[i + k] = ui + vi;
-        re[i + k + len / 2] = ur - vr;
-        im[i + k + len / 2] = ui - vi;
-        const ncr = cr * wr - ci * wi;
-        ci = cr * wi + ci * wr;
-        cr = ncr;
-      }
-    }
-  }
-}
-function spectrumPeaks(samples, sampleRate, maxPeaks = 8) {
-  const n = Math.min(65536, 1 << Math.floor(Math.log2(samples.length)));
-  const re = new Float64Array(n);
-  const im = new Float64Array(n);
-  for (let i = 0; i < n; i++) {
-    const w = 0.5 * (1 - Math.cos(2 * Math.PI * i / (n - 1)));
-    re[i] = samples[i] * w;
-  }
-  fft(re, im);
-  const half = n / 2;
-  const mag = new Float64Array(half);
-  let maxMag = 1e-30;
-  for (let i = 0; i < half; i++) {
-    mag[i] = Math.hypot(re[i], im[i]);
-    maxMag = Math.max(maxMag, mag[i]);
-  }
-  const peaks = [];
-  const floorDb = -60;
-  for (let i = 2; i < half - 2; i++) {
-    if (mag[i] <= mag[i - 1] || mag[i] < mag[i + 1]) continue;
-    const db = 20 * Math.log10(mag[i] / maxMag);
-    if (db < floorDb) continue;
-    const a = Math.log(mag[i - 1] + 1e-30);
-    const b = Math.log(mag[i] + 1e-30);
-    const c = Math.log(mag[i + 1] + 1e-30);
-    const delta = 0.5 * (a - c) / (a - 2 * b + c || 1e-30);
-    peaks.push({ hz: (i + delta) * sampleRate / n, db });
-  }
-  peaks.sort((a, b) => b.db - a.db);
-  const kept = [];
-  for (const p of peaks) {
-    if (kept.every((k) => Math.abs(p.hz - k.hz) / k.hz > 0.03)) kept.push(p);
-    if (kept.length >= maxPeaks) break;
-  }
-  return kept.sort((a, b) => a.hz - b.hz);
-}
-function noteToHz(note) {
-  const m = /^([A-Ga-g])([#b]?)(-?\d)$/.exec(note.trim());
-  if (!m) throw new Error(`unparseable note ${JSON.stringify(note)} \u2014 use e.g. "C6", "F#4"`);
-  const base = { C: 0, D: 2, E: 4, F: 5, G: 7, A: 9, B: 11 };
-  let semis = base[m[1].toUpperCase()];
-  if (m[2] === "#") semis += 1;
-  if (m[2] === "b") semis -= 1;
-  const midi = (Number(m[3]) + 1) * 12 + semis;
-  return 440 * 2 ** ((midi - 69) / 12);
-}
 function lookupMaterial(engine, name) {
   const key = name.trim().toLowerCase().replace(/[_ ]/g, "-");
   const registry4 = engine.getSheetMetalMaterials();
@@ -73045,10 +72630,8 @@ function barFromDocument(engine, documentId) {
   let holeDia = 0;
   for (const loop of holeLoops) {
     const cs = loop.map((p) => axis === "x" ? p.x : p.y);
-    const os = loop.map((p) => axis === "x" ? p.y : p.x);
     holesMm.push((Math.min(...cs) + Math.max(...cs)) / 2);
     holeDia = Math.max(holeDia, Math.max(...cs) - Math.min(...cs));
-    void os;
   }
   return {
     bar: {
@@ -73095,28 +72678,26 @@ function simulateStrike(input, engine) {
   const durationS = Math.min(Math.max(Number(a.duration_s ?? 2.5), 0.5), 10);
   const strikeFrac = Math.min(Math.max(Number(a.strike_position ?? 0.5), 0.02), 0.98);
   const contactMs = Math.min(Math.max(Number(a.mallet_contact_ms ?? 0.5), 0.05), 10);
-  const nModes = 6;
-  const closed = closedFormHz(bar, nModes);
-  const fem = femHz(bar, nModes);
-  const modes = strikeModes(bar, fem, strikeFrac, contactMs);
-  const audible = modes.filter((m) => m.hz < 0.45 * sampleRate && m.gain > 1e-4);
-  const samples = synthesize(audible, durationS, sampleRate);
-  const peaks = spectrumPeaks(samples, sampleRate);
   let expected;
   if (typeof a.expect_hz === "number") expected = a.expect_hz;
-  else if (typeof a.note === "string") expected = noteToHz(a.note);
-  const toleranceCents = Number(a.tolerance_cents ?? 10);
-  const dominant = peaks.reduce(
-    (best, p) => p.db > best.db ? p : best,
-    { hz: 0, db: -Infinity }
-  );
-  const verdict = expected !== void 0 ? {
-    expected_hz: expected,
-    measured_hz: dominant.hz,
-    cents_error: cents(dominant.hz, expected),
-    tolerance_cents: toleranceCents,
-    pass: Math.abs(cents(dominant.hz, expected)) <= toleranceCents
-  } : void 0;
+  else if (typeof a.note === "string") expected = engine.noteToHz(a.note);
+  const wantWav = typeof a.wav_filename === "string" && a.wav_filename.length > 0;
+  if (wantWav && !/\.wav$/i.test(a.wav_filename)) {
+    throw new Error("wav_filename must end in .wav");
+  }
+  const result = engine.simulateStrike({
+    bar,
+    strike_position: strikeFrac,
+    mallet_contact_ms: contactMs,
+    duration_s: durationS,
+    sample_rate: sampleRate,
+    n_modes: 6,
+    expected_hz: expected ?? null,
+    tolerance_cents: Number(a.tolerance_cents ?? 10),
+    include_wav: wantWav
+  });
+  const closed = result.closed_form_hz;
+  const fem = result.fem_hz;
   const payload = {
     bar: { ...bar, material: mat.name, material_display: mat.display_name },
     physics: {
@@ -73133,22 +72714,19 @@ function simulateStrike(input, engine) {
       duration_s: durationS,
       sample_rate: sampleRate
     },
-    modes: audible.map((m) => ({
+    modes: result.modes.map((m) => ({
       n: m.n,
       hz: round24(m.hz),
       gain: round4(m.gain),
       q: Math.round(m.q),
       t60_s: round24(m.t60_s)
     })),
-    spectrum_peaks: peaks.map((p) => ({ hz: round24(p.hz), db: round24(p.db) })),
-    ...verdict ? { verdict } : {},
+    spectrum_peaks: result.spectrum_peaks.map((p) => ({ hz: round24(p.hz), db: round24(p.db) })),
+    ...result.verdict ? { verdict: result.verdict } : {},
     limits: "1-D transverse bending only (no torsional/lateral modes); decay Q is a heuristic (material + cord-at-holes), frequencies are not"
   };
-  if (typeof a.wav_filename === "string" && a.wav_filename.length > 0) {
-    if (!/\.wav$/i.test(a.wav_filename)) {
-      throw new Error("wav_filename must end in .wav");
-    }
-    const bytes = encodeWav(samples, sampleRate);
+  if (wantWav && result.wav_base64) {
+    const bytes = Buffer.from(result.wav_base64, "base64");
     if (isRemoteDeployment()) {
       const cap = maxInlineExportBytes();
       if (bytes.length > cap) {
@@ -73159,7 +72737,7 @@ function simulateStrike(input, engine) {
       payload.wav = {
         filename: a.wav_filename,
         bytes: bytes.length,
-        data_base64: Buffer.from(bytes).toString("base64")
+        data_base64: result.wav_base64
       };
     } else {
       const path = resolveWithinRoot(
@@ -78325,8 +77903,8 @@ var init_server3 = __esm({
     init_order_feed();
     init_animate();
     PKG_VERSION = (() => {
-      if ("0.9.4-main.11") {
-        return "0.9.4-main.11";
+      if ("0.9.4-main.13") {
+        return "0.9.4-main.13";
       }
       try {
         const req = createRequire2(import.meta.url);
@@ -78335,8 +77913,8 @@ var init_server3 = __esm({
         return "0.0.0";
       }
     })();
-    BUILD_SHA = "681857ca2f62f173237de79ff2621cd1609480a7";
-    BUILD_TIME = "2026-07-24T16:25:54.297Z";
+    BUILD_SHA = "70552964e3760fb5b103235917b00364f6d45763";
+    BUILD_TIME = "2026-07-24T16:29:07.170Z";
     SHORT_SHA = BUILD_SHA === "unknown" ? "unknown" : BUILD_SHA.slice(0, 7);
     VERSION_WITH_BUILD = SHORT_SHA === "unknown" ? PKG_VERSION : `${PKG_VERSION}+${SHORT_SHA}`;
     INSTANCE_ID = randomUUID6().slice(0, 8);

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@vcad/mcp",
-  "version": "0.9.4-main.11",
+  "version": "0.9.4-main.13",
   "description": "vcad MCP server — parametric CAD + PCB design tools for AI agents (self-contained: bundled server + kernel WASM)",
   "license": "MIT",
   "type": "module",
@@ -19,7 +19,7 @@
   },
   "homepage": "https://vcad.io",
   "vcadBuild": {
-    "sha": "681857ca2f62f173237de79ff2621cd1609480a7",
-    "builtAt": "2026-07-24T16:25:54.297Z"
+    "sha": "70552964e3760fb5b103235917b00364f6d45763",
+    "builtAt": "2026-07-24T16:29:07.170Z"
   }
 }