@vcad/mcp 0.9.4-main.15 → 0.9.4-main.17

package/index.mjs

@@ -150,6 +150,15 @@ __export(vcad_kernel_wasm_exports, {
   ecadVerifyReceipt: () => ecadVerifyReceipt,
   ecadVerifySubstitution: () => ecadVerifySubstitution,
   emSimulate: () => emSimulate,
+  embroideryDesignToMesh: () => embroideryDesignToMesh,
+  enclosure_component_extents: () => enclosure_component_extents,
+  enclosure_connectors: () => enclosure_connectors,
+  enclosure_derive_board: () => enclosure_derive_board,
+  enclosure_features: () => enclosure_features,
+  enclosure_fit: () => enclosure_fit,
+  enclosure_mounting_holes: () => enclosure_mounting_holes,
+  enclosure_outline_aabb: () => enclosure_outline_aabb,
+  enclosure_to_world: () => enclosure_to_world,
   estimatePrintCost: () => estimatePrintCost,
   estimate_cost_for_process: () => estimate_cost_for_process,
   eulerXyzDegToQuat: () => eulerXyzDegToQuat,
@@ -257,6 +266,7 @@ __export(vcad_kernel_wasm_exports, {
   toleranceAnalyze: () => toleranceAnalyze,
   topologyOptimizeBox: () => topologyOptimizeBox,
   topologyOptimizeMesh: () => topologyOptimizeMesh,
+  transformMeshBuffers: () => transformMeshBuffers,
   writeEmbroideryDst: () => writeEmbroideryDst,
   writeEmbroideryPes: () => writeEmbroideryPes
 });
@@ -1519,6 +1529,193 @@ function emSimulate(spec_json, params_json, options_json) {
   }
   return takeFromExternrefTable0(ret[0]);
 }
+function embroideryDesignToMesh(design_json) {
+  const ptr0 = passStringToWasm0(design_json, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+  const len0 = WASM_VECTOR_LEN;
+  const ret = wasm.embroideryDesignToMesh(ptr0, len0);
+  if (ret[2]) {
+    throw takeFromExternrefTable0(ret[1]);
+  }
+  return takeFromExternrefTable0(ret[0]);
+}
+function enclosure_component_extents(meshes_json, pcb_json) {
+  let deferred4_0;
+  let deferred4_1;
+  try {
+    const ptr0 = passStringToWasm0(meshes_json, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+    const len0 = WASM_VECTOR_LEN;
+    const ptr1 = passStringToWasm0(pcb_json, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+    const len1 = WASM_VECTOR_LEN;
+    const ret = wasm.enclosure_component_extents(ptr0, len0, ptr1, len1);
+    var ptr3 = ret[0];
+    var len3 = ret[1];
+    if (ret[3]) {
+      ptr3 = 0;
+      len3 = 0;
+      throw takeFromExternrefTable0(ret[2]);
+    }
+    deferred4_0 = ptr3;
+    deferred4_1 = len3;
+    return getStringFromWasm0(ptr3, len3);
+  } finally {
+    wasm.__wbindgen_free(deferred4_0, deferred4_1, 1);
+  }
+}
+function enclosure_connectors(pcb_json, outline_json) {
+  let deferred4_0;
+  let deferred4_1;
+  try {
+    const ptr0 = passStringToWasm0(pcb_json, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+    const len0 = WASM_VECTOR_LEN;
+    const ptr1 = passStringToWasm0(outline_json, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+    const len1 = WASM_VECTOR_LEN;
+    const ret = wasm.enclosure_connectors(ptr0, len0, ptr1, len1);
+    var ptr3 = ret[0];
+    var len3 = ret[1];
+    if (ret[3]) {
+      ptr3 = 0;
+      len3 = 0;
+      throw takeFromExternrefTable0(ret[2]);
+    }
+    deferred4_0 = ptr3;
+    deferred4_1 = len3;
+    return getStringFromWasm0(ptr3, len3);
+  } finally {
+    wasm.__wbindgen_free(deferred4_0, deferred4_1, 1);
+  }
+}
+function enclosure_derive_board(cavity_json, standoffs_json, options_json) {
+  let deferred5_0;
+  let deferred5_1;
+  try {
+    const ptr0 = passStringToWasm0(cavity_json, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+    const len0 = WASM_VECTOR_LEN;
+    const ptr1 = passStringToWasm0(standoffs_json, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+    const len1 = WASM_VECTOR_LEN;
+    const ptr2 = passStringToWasm0(options_json, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+    const len2 = WASM_VECTOR_LEN;
+    const ret = wasm.enclosure_derive_board(ptr0, len0, ptr1, len1, ptr2, len2);
+    var ptr4 = ret[0];
+    var len4 = ret[1];
+    if (ret[3]) {
+      ptr4 = 0;
+      len4 = 0;
+      throw takeFromExternrefTable0(ret[2]);
+    }
+    deferred5_0 = ptr4;
+    deferred5_1 = len4;
+    return getStringFromWasm0(ptr4, len4);
+  } finally {
+    wasm.__wbindgen_free(deferred5_0, deferred5_1, 1);
+  }
+}
+function enclosure_features(positions, indices) {
+  let deferred4_0;
+  let deferred4_1;
+  try {
+    const ptr0 = passArrayF64ToWasm0(positions, wasm.__wbindgen_malloc);
+    const len0 = WASM_VECTOR_LEN;
+    const ptr1 = passArray32ToWasm0(indices, wasm.__wbindgen_malloc);
+    const len1 = WASM_VECTOR_LEN;
+    const ret = wasm.enclosure_features(ptr0, len0, ptr1, len1);
+    var ptr3 = ret[0];
+    var len3 = ret[1];
+    if (ret[3]) {
+      ptr3 = 0;
+      len3 = 0;
+      throw takeFromExternrefTable0(ret[2]);
+    }
+    deferred4_0 = ptr3;
+    deferred4_1 = len3;
+    return getStringFromWasm0(ptr3, len3);
+  } finally {
+    wasm.__wbindgen_free(deferred4_0, deferred4_1, 1);
+  }
+}
+function enclosure_fit(input_json) {
+  let deferred3_0;
+  let deferred3_1;
+  try {
+    const ptr0 = passStringToWasm0(input_json, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+    const len0 = WASM_VECTOR_LEN;
+    const ret = wasm.enclosure_fit(ptr0, len0);
+    var ptr2 = ret[0];
+    var len2 = ret[1];
+    if (ret[3]) {
+      ptr2 = 0;
+      len2 = 0;
+      throw takeFromExternrefTable0(ret[2]);
+    }
+    deferred3_0 = ptr2;
+    deferred3_1 = len2;
+    return getStringFromWasm0(ptr2, len2);
+  } finally {
+    wasm.__wbindgen_free(deferred3_0, deferred3_1, 1);
+  }
+}
+function enclosure_mounting_holes(pcb_json) {
+  let deferred3_0;
+  let deferred3_1;
+  try {
+    const ptr0 = passStringToWasm0(pcb_json, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+    const len0 = WASM_VECTOR_LEN;
+    const ret = wasm.enclosure_mounting_holes(ptr0, len0);
+    var ptr2 = ret[0];
+    var len2 = ret[1];
+    if (ret[3]) {
+      ptr2 = 0;
+      len2 = 0;
+      throw takeFromExternrefTable0(ret[2]);
+    }
+    deferred3_0 = ptr2;
+    deferred3_1 = len2;
+    return getStringFromWasm0(ptr2, len2);
+  } finally {
+    wasm.__wbindgen_free(deferred3_0, deferred3_1, 1);
+  }
+}
+function enclosure_outline_aabb(outline_json) {
+  let deferred3_0;
+  let deferred3_1;
+  try {
+    const ptr0 = passStringToWasm0(outline_json, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+    const len0 = WASM_VECTOR_LEN;
+    const ret = wasm.enclosure_outline_aabb(ptr0, len0);
+    var ptr2 = ret[0];
+    var len2 = ret[1];
+    if (ret[3]) {
+      ptr2 = 0;
+      len2 = 0;
+      throw takeFromExternrefTable0(ret[2]);
+    }
+    deferred3_0 = ptr2;
+    deferred3_1 = len2;
+    return getStringFromWasm0(ptr2, len2);
+  } finally {
+    wasm.__wbindgen_free(deferred3_0, deferred3_1, 1);
+  }
+}
+function enclosure_to_world(x, y, z, placement_json) {
+  let deferred3_0;
+  let deferred3_1;
+  try {
+    const ptr0 = passStringToWasm0(placement_json, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+    const len0 = WASM_VECTOR_LEN;
+    const ret = wasm.enclosure_to_world(x, y, z, ptr0, len0);
+    var ptr2 = ret[0];
+    var len2 = ret[1];
+    if (ret[3]) {
+      ptr2 = 0;
+      len2 = 0;
+      throw takeFromExternrefTable0(ret[2]);
+    }
+    deferred3_0 = ptr2;
+    deferred3_1 = len2;
+    return getStringFromWasm0(ptr2, len2);
+  } finally {
+    wasm.__wbindgen_free(deferred3_0, deferred3_1, 1);
+  }
+}
 function estimatePrintCost(volume_mm3, infill_density, wall_count, line_width, material_name) {
   const ptr0 = passStringToWasm0(material_name, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
   const len0 = WASM_VECTOR_LEN;
@@ -2946,6 +3143,19 @@ function topologyOptimizeMesh(spec_json, positions, indices) {
   }
   return takeFromExternrefTable0(ret[0]);
 }
+function transformMeshBuffers(positions, normals, transform_json) {
+  const ptr0 = passArrayF32ToWasm0(positions, wasm.__wbindgen_malloc);
+  const len0 = WASM_VECTOR_LEN;
+  var ptr1 = isLikeNone(normals) ? 0 : passArrayF32ToWasm0(normals, wasm.__wbindgen_malloc);
+  var len1 = WASM_VECTOR_LEN;
+  const ptr2 = passStringToWasm0(transform_json, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+  const len2 = WASM_VECTOR_LEN;
+  const ret = wasm.transformMeshBuffers(ptr0, len0, ptr1, len1, ptr2, len2);
+  if (ret[2]) {
+    throw takeFromExternrefTable0(ret[1]);
+  }
+  return takeFromExternrefTable0(ret[0]);
+}
 function writeEmbroideryDst(json) {
   const ptr0 = passStringToWasm0(json, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
   const len0 = WASM_VECTOR_LEN;
@@ -10895,753 +11105,57 @@ var init_atoms = __esm({
 });
 
 // ../engine/dist/enclosure-fit.js
-function outlineAabb(outline) {
-  let minX = Infinity;
-  let maxX = -Infinity;
-  let minY = Infinity;
-  let maxY = -Infinity;
-  for (const v of outline.vertices) {
-    if (v.x < minX)
-      minX = v.x;
-    if (v.x > maxX)
-      maxX = v.x;
-    if (v.y < minY)
-      minY = v.y;
-    if (v.y > maxY)
-      maxY = v.y;
-  }
-  return { minX, maxX, minY, maxY };
+function enclosureWasm() {
+  const wasm3 = getKernelWasmSync();
+  if (!wasm3?.enclosure_fit) {
+    throw new Error("enclosure checks require the kernel WASM to be initialized \u2014 await getKernelWasm() (or Engine.init()) first");
+  }
+  return wasm3;
 }
-function toWorld(p, placement) {
-  const t2 = placement.rotationDeg * Math.PI / 180;
-  const cos = Math.cos(t2);
-  const sin = Math.sin(t2);
-  return {
-    x: placement.offset.x + p.x * cos - p.y * sin,
-    y: placement.offset.y + p.x * sin + p.y * cos,
-    z: placement.offset.z + (p.z ?? 0)
-  };
+function compact(obj) {
+  return Object.fromEntries(Object.entries(obj).filter(([, v]) => v !== void 0));
 }
-function nearestEdge(x, y, b) {
-  const d = [
-    ["minX", Math.abs(x - b.minX)],
-    ["maxX", Math.abs(x - b.maxX)],
-    ["minY", Math.abs(y - b.minY)],
-    ["maxY", Math.abs(y - b.maxY)]
-  ];
-  d.sort((a, c) => a[1] - c[1]);
-  return d[0][0];
+function outlineAabb(outline) {
+  return JSON.parse(enclosureWasm().enclosure_outline_aabb(JSON.stringify(outline)));
 }
-function padWorld(fp, pad) {
-  const t2 = (fp.rotation ?? 0) * Math.PI / 180;
-  const cos = Math.cos(t2);
-  const sin = Math.sin(t2);
-  return {
-    x: fp.position.x + pad.position.x * cos - pad.position.y * sin,
-    y: fp.position.y + pad.position.x * sin + pad.position.y * cos
-  };
+function toWorld(p, placement) {
+  return JSON.parse(enclosureWasm().enclosure_to_world(p.x, p.y, p.z ?? 0, JSON.stringify(placement)));
 }
 function mountingHolesFromPcb(pcb) {
-  const holes = [];
-  for (const fp of pcb.footprints) {
-    const isMount = MOUNT_RE.test(fp.footprintName) || MOUNT_RE.test(fp.ref);
-    if (isMount) {
-      const pad = fp.pads[0];
-      let dia = 3.2;
-      if (pad) {
-        if (pad.drill && typeof pad.drill === "object") {
-          const dd = pad.drill;
-          if (typeof dd.diameter === "number")
-            dia = dd.diameter;
-        } else if (pad.shape.type === "Circle") {
-          dia = pad.shape.diameter;
-        }
-      }
-      holes.push({ x: round2(fp.position.x), y: round2(fp.position.y), diameter: round2(dia), ref: fp.ref });
-      continue;
-    }
-    for (const pad of fp.pads) {
-      if (pad.padType !== "NPTH")
-        continue;
-      const w = padWorld(fp, pad);
-      let dia = 3.2;
-      const dd = pad.drill;
-      if (dd && typeof dd.diameter === "number")
-        dia = dd.diameter;
-      else if (pad.shape.type === "Circle")
-        dia = pad.shape.diameter;
-      holes.push({ x: round2(w.x), y: round2(w.y), diameter: round2(dia), ref: fp.ref });
-    }
-  }
-  return holes;
+  return JSON.parse(enclosureWasm().enclosure_mounting_holes(JSON.stringify(pcb)));
 }
 function componentExtentsFromMeshes(meshes, pcb) {
-  const frontByRef = new Map(pcb.footprints.map((fp) => [fp.ref, fp.front ?? true]));
-  const out = [];
-  for (const m of meshes) {
-    let minZ = Infinity;
-    let maxZ = -Infinity;
-    for (let i = 2; i < m.positions.length; i += 3) {
-      const z = m.positions[i];
-      if (z < minZ)
-        minZ = z;
-      if (z > maxZ)
-        maxZ = z;
-    }
-    if (!Number.isFinite(minZ))
-      continue;
-    out.push({
-      ref: m.footprint_ref,
-      front: frontByRef.get(m.footprint_ref) ?? true,
-      topZ: round2(maxZ),
-      bottomZ: round2(minZ)
-    });
-  }
-  return out;
+  const payload = meshes.map((m) => ({
+    footprint_ref: m.footprint_ref,
+    positions: Array.from(m.positions)
+  }));
+  return JSON.parse(enclosureWasm().enclosure_component_extents(JSON.stringify(payload), JSON.stringify(pcb)));
 }
 function connectorsFromPcb(pcb, outline) {
-  const aabb = outlineAabb(outline);
-  const out = [];
-  for (const fp of pcb.footprints) {
-    const isConn = CONNECTOR_REF_RE.test(fp.ref) || CONNECTOR_NAME_RE.test(fp.footprintName) || CONNECTOR_NAME_RE.test(fp.value ?? "");
-    if (!isConn)
-      continue;
-    out.push({
-      ref: fp.ref,
-      x: round2(fp.position.x),
-      y: round2(fp.position.y),
-      edge: nearestEdge(fp.position.x, fp.position.y, aabb),
-      height: 0
-    });
-  }
-  return out;
-}
-function autoPlacement(input, clearance) {
-  const { cavity } = input;
-  const a = outlineAabb(input.outline);
-  const boardW = a.maxX - a.minX;
-  const boardH = a.maxY - a.minY;
-  const cavW = cavity.maxX - cavity.minX;
-  const cavH = cavity.maxY - cavity.minY;
-  const offX = cavity.minX + (cavW - boardW) / 2 - a.minX;
-  const offY = cavity.minY + (cavH - boardH) / 2 - a.minY;
-  const standoffTop = input.standoffs && input.standoffs.length > 0 ? Math.max(...input.standoffs.map((s) => s.topZ)) : cavity.floorZ + (input.standoffHeight ?? 0);
-  void clearance;
-  return { offset: { x: round2(offX), y: round2(offY), z: round2(standoffTop) }, rotationDeg: 0 };
-}
-function checkBoardFit(input, placement, clearance) {
-  const { cavity } = input;
-  let minX = Infinity;
-  let maxX = -Infinity;
-  let minY = Infinity;
-  let maxY = -Infinity;
-  for (const v of input.outline.vertices) {
-    const w = toWorld(v, placement);
-    if (w.x < minX)
-      minX = w.x;
-    if (w.x > maxX)
-      maxX = w.x;
-    if (w.y < minY)
-      minY = w.y;
-    if (w.y > maxY)
-      maxY = w.y;
-  }
-  const marginMinX = minX - cavity.minX;
-  const marginMaxX = cavity.maxX - maxX;
-  const marginMinY = minY - cavity.minY;
-  const marginMaxY = cavity.maxY - maxY;
-  const worst2 = Math.min(marginMinX, marginMaxX, marginMinY, marginMaxY);
-  const sides = [
-    ["-X", marginMinX],
-    ["+X", marginMaxX],
-    ["-Y", marginMinY],
-    ["+Y", marginMaxY]
-  ];
-  const tight = sides.filter(([, m]) => m < clearance).map(([s]) => s);
-  const ok4 = worst2 >= clearance - 1e-6;
-  return {
-    id: "board_fit",
-    label: "Board fits cavity with clearance",
-    status: ok4 ? "pass" : "fail",
-    detail: ok4 ? `Board fits with ${round2(worst2)}mm worst-case clearance (need ${clearance}mm)` : worst2 < 0 ? `Board overhangs the cavity by ${round2(-worst2)}mm on ${tight.join(", ")}` : `Clearance on ${tight.join(", ")} is ${round2(worst2)}mm < required ${clearance}mm`,
-    measurements: {
-      worst_clearance_mm: round2(worst2),
-      margin_minus_x: round2(marginMinX),
-      margin_plus_x: round2(marginMaxX),
-      margin_minus_y: round2(marginMinY),
-      margin_plus_y: round2(marginMaxY),
-      board_w: round2(maxX - minX),
-      board_h: round2(maxY - minY),
-      cavity_w: round2(cavity.maxX - cavity.minX),
-      cavity_h: round2(cavity.maxY - cavity.minY)
-    }
-  };
-}
-function checkLidClearance(input, placement, clearance) {
-  const { cavity } = input;
-  const extents = input.componentExtents ?? [];
-  if (extents.length === 0) {
-    return {
-      id: "lid_clearance",
-      label: "Components clear the lid",
-      status: "skip",
-      detail: "No component heights available (kernel component meshes unavailable)"
-    };
-  }
-  const cavityDepth = cavity.ceilZ - cavity.floorZ;
-  const front = extents.filter((e) => e.front);
-  const back = extents.filter((e) => !e.front);
-  let tallest = { ref: "", top: -Infinity };
-  for (const e of front) {
-    const top = placement.offset.z + e.topZ;
-    if (top > tallest.top)
-      tallest = { ref: e.ref, top };
-  }
-  const lidGap = cavity.ceilZ - tallest.top;
-  const topOk = front.length === 0 || lidGap >= clearance - 1e-6;
-  let lowest = { ref: "", bottom: Infinity };
-  for (const e of back) {
-    const bot = placement.offset.z + e.bottomZ;
-    if (bot < lowest.bottom)
-      lowest = { ref: e.ref, bottom: bot };
-  }
-  const floorGap = back.length > 0 ? lowest.bottom - cavity.floorZ : Infinity;
-  const botOk = back.length === 0 || floorGap >= -1e-6;
-  const ok4 = topOk && botOk;
-  let detail;
-  if (ok4) {
-    detail = `Tallest part ${tallest.ref || "\u2014"} leaves ${round2(lidGap)}mm under the lid (cavity depth ${round2(cavityDepth)}mm)`;
-    if (back.length > 0 && Number.isFinite(floorGap)) {
-      detail += `; back-side ${lowest.ref} clears floor by ${round2(floorGap)}mm`;
-    }
-  } else if (!topOk) {
-    detail = `${tallest.ref} is ${round2(-lidGap + clearance)}mm too tall \u2014 it ${lidGap < 0 ? "punches through" : "is within clearance of"} the lid (cavity depth ${round2(cavityDepth)}mm)`;
-  } else {
-    detail = `Back-side ${lowest.ref} collides with the floor by ${round2(-floorGap)}mm \u2014 raise the standoffs`;
-  }
-  return {
-    id: "lid_clearance",
-    label: "Components clear the lid",
-    status: ok4 ? "pass" : "fail",
-    detail,
-    measurements: {
-      cavity_depth_mm: round2(cavityDepth),
-      tallest_ref: tallest.ref || "none",
-      lid_gap_mm: Number.isFinite(lidGap) ? round2(lidGap) : "n/a",
-      stack_top_z: Number.isFinite(tallest.top) ? round2(tallest.top) : "n/a",
-      floor_gap_mm: Number.isFinite(floorGap) ? round2(floorGap) : "n/a"
-    }
-  };
-}
-function checkMountingHoles(input, placement) {
-  const holes = input.mountingHoles ?? [];
-  const standoffs = input.standoffs ?? [];
-  const tol = input.holeTolerance ?? DEFAULT_HOLE_TOL;
-  if (holes.length === 0) {
-    return {
-      id: "mounting_holes",
-      label: "Mounting holes land on standoffs",
-      status: "skip",
-      detail: "Board declares no mounting holes"
-    };
-  }
-  if (standoffs.length === 0) {
-    return {
-      id: "mounting_holes",
-      label: "Mounting holes land on standoffs",
-      status: "skip",
-      detail: `Board has ${holes.length} mounting hole(s) but no standoffs were detected in the enclosure`
-    };
-  }
-  let matched = 0;
-  let worst2 = 0;
-  const misses = [];
-  for (const h of holes) {
-    const w = toWorld(h, placement);
-    let best = Infinity;
-    for (const s of standoffs) {
-      const d = Math.hypot(w.x - s.x, w.y - s.y);
-      if (d < best)
-        best = d;
-    }
-    if (best <= tol) {
-      matched++;
-      if (best > worst2)
-        worst2 = best;
-    } else {
-      misses.push(`${h.ref ?? "hole"}@(${round2(w.x)},${round2(w.y)}) is ${round2(best)}mm off`);
-    }
-  }
-  const ok4 = matched === holes.length;
-  return {
-    id: "mounting_holes",
-    label: "Mounting holes land on standoffs",
-    status: ok4 ? "pass" : "fail",
-    detail: ok4 ? `All ${holes.length} mounting holes align to standoffs (worst offset ${round2(worst2)}mm, tol ${tol}mm)` : `${matched}/${holes.length} holes align \u2014 ${misses.join("; ")}`,
-    measurements: {
-      holes_total: holes.length,
-      holes_matched: matched,
-      standoffs: standoffs.length,
-      tolerance_mm: tol,
-      worst_offset_mm: round2(worst2)
-    }
-  };
-}
-function checkConnectors(input, placement, clearance) {
-  const conns = input.connectors ?? [];
-  const openings = input.openings ?? [];
-  if (conns.length === 0) {
-    return {
-      id: "connector_cutouts",
-      label: "Connectors align to wall cutouts",
-      status: "skip",
-      detail: "Board declares no edge connectors"
-    };
-  }
-  const cav = input.cavity;
-  let aligned = 0;
-  const problems = [];
-  for (const c of conns) {
-    const w = toWorld(c, placement);
-    const wallEdge = nearestEdge(w.x, w.y, cav);
-    const along = wallEdge === "minX" || wallEdge === "maxX" ? w.y : w.x;
-    const onWall = openings.filter((o) => o.edge === wallEdge);
-    if (onWall.length === 0) {
-      problems.push(`${c.ref} faces the ${wallEdge} wall but it has no cutout`);
-      continue;
-    }
-    const hit = onWall.find((o) => {
-      const oc = o.edge === "minX" || o.edge === "maxX" ? o.center.y : o.center.x;
-      return Math.abs(along - oc) <= o.width / 2 + clearance;
-    });
-    if (hit) {
-      aligned++;
-    } else {
-      const nearest = onWall.reduce((best, o) => {
-        const oc = o.edge === "minX" || o.edge === "maxX" ? o.center.y : o.center.x;
-        const off = Math.abs(along - oc) - o.width / 2;
-        return off < best ? off : best;
-      }, Infinity);
-      problems.push(`${c.ref} on ${wallEdge} misses its cutout by ${round2(nearest)}mm`);
-    }
-  }
-  const ok4 = aligned === conns.length;
-  const status = ok4 ? "pass" : openings.length === 0 ? "warn" : "fail";
-  return {
-    id: "connector_cutouts",
-    label: "Connectors align to wall cutouts",
-    status,
-    detail: ok4 ? `All ${conns.length} connector(s) line up with wall cutouts` : openings.length === 0 ? `No wall cutouts detected; ${conns.length} connector(s) would be enclosed: ${problems.join("; ")}` : `${aligned}/${conns.length} connectors aligned \u2014 ${problems.join("; ")}`,
-    measurements: {
-      connectors_total: conns.length,
-      connectors_aligned: aligned,
-      wall_openings: openings.length
-    }
-  };
+  return JSON.parse(enclosureWasm().enclosure_connectors(JSON.stringify(pcb), JSON.stringify(outline)));
 }
 function checkEnclosureFit(input) {
-  const clearance = input.clearance ?? DEFAULT_CLEARANCE;
-  const placement = input.placement ?? autoPlacement(input, clearance);
-  const checks = [
-    checkBoardFit(input, placement, clearance),
-    checkLidClearance(input, placement, clearance),
-    checkMountingHoles(input, placement),
-    checkConnectors(input, placement, clearance)
-  ];
-  const failed = checks.filter((c) => c.status === "fail");
-  const warned = checks.filter((c) => c.status === "warn");
-  const passed = checks.filter((c) => c.status === "pass");
-  const ok4 = failed.length === 0;
-  const verified = ok4 && warned.length === 0;
-  let summary;
-  if (failed.length > 0) {
-    summary = `Enclosure fit: FAIL \u2014 ${failed.map((c) => c.label.toLowerCase()).join("; ")}`;
-  } else if (warned.length > 0) {
-    summary = `Enclosure fit: UNVERIFIED \u2014 ${passed.length} passed, ${warned.length} warning(s): ${warned.map((c) => c.detail).join("; ")}`;
-  } else {
-    summary = `Enclosure fit: PASS \u2014 ${passed.length}/${checks.length} checks (${checks.filter((c) => c.status === "pass").map((c) => c.label.toLowerCase()).join(", ")})`;
-  }
-  return { ok: ok4, verified, summary, clearance, placement, checks };
+  return JSON.parse(enclosureWasm().enclosure_fit(JSON.stringify(compact(input))));
 }
 function deriveBoardFromCavity(cavity, standoffs, opts = {}) {
-  const clearance = opts.clearance ?? DEFAULT_CLEARANCE;
-  const thickness = opts.thickness ?? 1.6;
-  const holeDia = opts.holeDiameter ?? 3.2;
-  const w = round2(cavity.maxX - cavity.minX - 2 * clearance);
-  const h = round2(cavity.maxY - cavity.minY - 2 * clearance);
-  const outline = {
-    vertices: [
-      { x: 0, y: 0 },
-      { x: w, y: 0 },
-      { x: w, y: h },
-      { x: 0, y: h }
-    ],
-    thickness
-  };
-  const standoffTop = standoffs.length > 0 ? Math.max(...standoffs.map((s) => s.topZ)) : cavity.floorZ + (opts.standoffHeight ?? 0);
-  const offX = cavity.minX + clearance;
-  const offY = cavity.minY + clearance;
-  const placement = {
-    offset: { x: round2(offX), y: round2(offY), z: round2(standoffTop) },
-    rotationDeg: 0
-  };
-  const mountingHoles = [];
-  for (const s of standoffs) {
-    const lx = round2(s.x - offX);
-    const ly = round2(s.y - offY);
-    if (lx >= 0 && lx <= w && ly >= 0 && ly <= h) {
-      mountingHoles.push({ x: lx, y: ly, diameter: holeDia });
-    }
-  }
-  return { outline, mountingHoles, placement };
+  return JSON.parse(enclosureWasm().enclosure_derive_board(JSON.stringify(cavity), JSON.stringify(standoffs), JSON.stringify(compact(opts))));
 }
-var DEFAULT_CLEARANCE, DEFAULT_HOLE_TOL, round2, MOUNT_RE, CONNECTOR_REF_RE, CONNECTOR_NAME_RE;
 var init_enclosure_fit = __esm({
   "../engine/dist/enclosure-fit.js"() {
     "use strict";
-    DEFAULT_CLEARANCE = 0.5;
-    DEFAULT_HOLE_TOL = 0.6;
-    round2 = (n) => Math.round(n * 100) / 100;
-    MOUNT_RE = /mount(ing)?[_-]?hole|mountingpad|mounthole/i;
-    CONNECTOR_REF_RE = /^(J|CN|CON|USB|P)\d/i;
-    CONNECTOR_NAME_RE = /usb|type[-_]?c|micro|mini|conn|header|jst|molex|rj45|hdr|terminal|socket|receptacle|barrel|dcjack/i;
+    init_wasm_singleton();
   }
 });
 
 // ../engine/dist/enclosure-mesh.js
-function gwn(positions, indices, qx, qy, qz) {
-  let w = 0;
-  const tris = indices.length;
-  for (let t2 = 0; t2 < tris; t2 += 3) {
-    const i0 = indices[t2] * 3;
-    const i1 = indices[t2 + 1] * 3;
-    const i2 = indices[t2 + 2] * 3;
-    const ax = positions[i0] - qx;
-    const ay = positions[i0 + 1] - qy;
-    const az = positions[i0 + 2] - qz;
-    const bx = positions[i1] - qx;
-    const by = positions[i1 + 1] - qy;
-    const bz = positions[i1 + 2] - qz;
-    const cx = positions[i2] - qx;
-    const cy = positions[i2 + 1] - qy;
-    const cz = positions[i2 + 2] - qz;
-    const la = Math.sqrt(ax * ax + ay * ay + az * az);
-    const lb = Math.sqrt(bx * bx + by * by + bz * bz);
-    const lc = Math.sqrt(cx * cx + cy * cy + cz * cz);
-    const cbx = by * cz - bz * cy;
-    const cby = bz * cx - bx * cz;
-    const cbz = bx * cy - by * cx;
-    const num2 = ax * cbx + ay * cby + az * cbz;
-    const den = la * lb * lc + (ax * bx + ay * by + az * bz) * lc + (bx * cx + by * cy + bz * cz) * la + (cx * ax + cy * ay + cz * az) * lb;
-    w += Math.atan2(num2, den);
-  }
-  return w / (2 * Math.PI);
-}
-function buildOccupancy(positions, indices) {
-  let minX = Infinity;
-  let maxX = -Infinity;
-  let minY = Infinity;
-  let maxY = -Infinity;
-  let minZ = Infinity;
-  let maxZ = -Infinity;
-  for (let i = 0; i < positions.length; i += 3) {
-    const x = positions[i];
-    const y = positions[i + 1];
-    const z = positions[i + 2];
-    if (x < minX)
-      minX = x;
-    if (x > maxX)
-      maxX = x;
-    if (y < minY)
-      minY = y;
-    if (y > maxY)
-      maxY = y;
-    if (z < minZ)
-      minZ = z;
-    if (z > maxZ)
-      maxZ = z;
-  }
-  if (!(maxX > minX) || !(maxY > minY) || !(maxZ > minZ))
-    return null;
-  const gw = GRID_XY;
-  const gh = GRID_XY;
-  const gz = GRID_Z;
-  const dx = (maxX - minX) / gw;
-  const dy = (maxY - minY) / gh;
-  const dz = (maxZ - minZ) / gz;
-  const occ = new Uint8Array(gw * gh * gz);
-  for (let k = 0; k < gz; k++) {
-    const qz = minZ + (k + JZ) * dz;
-    for (let j = 0; j < gh; j++) {
-      const qy = minY + (j + JY) * dy;
-      for (let i = 0; i < gw; i++) {
-        const qx = minX + (i + JX) * dx;
-        if (Math.abs(gwn(positions, indices, qx, qy, qz)) > INSIDE) {
-          occ[i + gw * (j + gh * k)] = 1;
-        }
-      }
-    }
-  }
-  return { gw, gh, gz, ox: minX, oy: minY, oz: minZ, dx, dy, dz, minZ, maxZ, occ };
-}
-function bottomRunTopZ(o, i, j) {
-  const base = i + o.gw * j;
-  if (o.occ[base] !== 1)
-    return null;
-  let k = 1;
-  while (k < o.gz && o.occ[base + o.gw * o.gh * k] === 1)
-    k++;
-  return o.oz + k * o.dz;
-}
-function openAtTop(o, i, j) {
-  return o.occ[i + o.gw * (j + o.gh * (o.gz - 1))] !== 1;
-}
-function hasGap(o, i, j, kLo, kHi) {
-  for (let k = kLo; k <= kHi; k++) {
-    if (o.occ[i + o.gw * (j + o.gh * k)] !== 1)
-      return true;
-  }
-  return false;
-}
 function extractEnclosureFeatures(positions, indices) {
-  const o = buildOccupancy(positions, indices);
-  if (!o) {
-    return {
-      outer: { minX: 0, maxX: 0, minY: 0, maxY: 0, minZ: 0, maxZ: 0 },
-      cavity: null,
-      standoffs: [],
-      openings: []
-    };
-  }
-  const outer = {
-    minX: o.ox,
-    maxX: o.ox + o.gw * o.dx,
-    minY: o.oy,
-    maxY: o.oy + o.gh * o.dy,
-    minZ: o.minZ,
-    maxZ: o.maxZ
-  };
-  const isPocket = new Uint8Array(o.gw * o.gh);
-  const firstTops = [];
-  const countX = new Int32Array(o.gw);
-  const countY = new Int32Array(o.gh);
-  let pocketCount = 0;
-  for (let j = 0; j < o.gh; j++) {
-    for (let i = 0; i < o.gw; i++) {
-      const top = bottomRunTopZ(o, i, j);
-      if (top != null && openAtTop(o, i, j)) {
-        isPocket[i + o.gw * j] = 1;
-        pocketCount++;
-        firstTops.push(top);
-        countX[i]++;
-        countY[j]++;
-      }
-    }
-  }
-  if (pocketCount === 0) {
-    return { outer, cavity: null, standoffs: [], openings: [] };
-  }
-  const core = (counts) => {
-    let peak = 0;
-    for (const c of counts)
-      if (c > peak)
-        peak = c;
-    const thr = peak * 0.5;
-    let lo = 0;
-    let hi = counts.length - 1;
-    while (lo < counts.length && counts[lo] < thr)
-      lo++;
-    while (hi >= 0 && counts[hi] < thr)
-      hi--;
-    return [lo, hi];
-  };
-  const [pminI, pmaxI] = core(countX);
-  const [pminJ, pmaxJ] = core(countY);
-  const sortedTops = [...firstTops].sort((a, b) => a - b);
-  const floorZ = sortedTops[Math.floor(sortedTops.length / 2)] ?? o.minZ;
-  const ceilZ = o.maxZ;
-  const cavity = {
-    minX: o.ox + pminI * o.dx,
-    maxX: o.ox + (pmaxI + 1) * o.dx,
-    minY: o.oy + pminJ * o.dy,
-    maxY: o.oy + (pmaxJ + 1) * o.dy,
-    floorZ: round22(floorZ),
-    ceilZ: round22(ceilZ),
-    hasLid: false
-  };
-  const standoffs = extractStandoffs(o, isPocket, floorZ);
-  const openings = extractOpenings(o, cavity, floorZ, ceilZ);
-  return { outer, cavity, standoffs, openings };
-}
-function extractStandoffs(o, isPocket, floorZ) {
-  const isPost = new Uint8Array(o.gw * o.gh);
-  for (let j = 0; j < o.gh; j++) {
-    for (let i = 0; i < o.gw; i++) {
-      if (!isPocket[i + o.gw * j])
-        continue;
-      const top = bottomRunTopZ(o, i, j);
-      if (top != null && top > floorZ + MIN_POST_HEIGHT)
-        isPost[i + o.gw * j] = 1;
-    }
-  }
-  const seen = new Uint8Array(o.gw * o.gh);
-  const standoffs = [];
-  for (let j = 0; j < o.gh; j++) {
-    for (let i = 0; i < o.gw; i++) {
-      const k0 = i + o.gw * j;
-      if (!isPost[k0] || seen[k0])
-        continue;
-      const stack = [[i, j]];
-      seen[k0] = 1;
-      let sumX = 0;
-      let sumY = 0;
-      let n = 0;
-      let topMax = -Infinity;
-      let minI = i;
-      let maxI = i;
-      let minJ = j;
-      let maxJ = j;
-      while (stack.length) {
-        const [ci, cj] = stack.pop();
-        sumX += o.ox + (ci + 0.5) * o.dx;
-        sumY += o.oy + (cj + 0.5) * o.dy;
-        n++;
-        const t2 = bottomRunTopZ(o, ci, cj);
-        if (t2 != null && t2 > topMax)
-          topMax = t2;
-        if (ci < minI)
-          minI = ci;
-        if (ci > maxI)
-          maxI = ci;
-        if (cj < minJ)
-          minJ = cj;
-        if (cj > maxJ)
-          maxJ = cj;
-        for (const [ni, nj] of [
-          [ci - 1, cj],
-          [ci + 1, cj],
-          [ci, cj - 1],
-          [ci, cj + 1]
-        ]) {
-          if (ni < 0 || nj < 0 || ni >= o.gw || nj >= o.gh)
-            continue;
-          const nk = ni + o.gw * nj;
-          if (isPost[nk] && !seen[nk]) {
-            seen[nk] = 1;
-            stack.push([ni, nj]);
-          }
-        }
-      }
-      if (n < 2)
-        continue;
-      const radius = Math.max(((maxI - minI + 1) * o.dx + (maxJ - minJ + 1) * o.dy) / 4, o.dx);
-      standoffs.push({
-        x: round22(sumX / n),
-        y: round22(sumY / n),
-        topZ: round22(topMax),
-        radius: round22(radius)
-      });
-    }
-  }
-  return standoffs;
-}
-function extractOpenings(o, cavity, floorZ, ceilZ) {
-  const pIminI = Math.round((cavity.minX - o.ox) / o.dx);
-  const pImaxI = Math.round((cavity.maxX - o.ox) / o.dx) - 1;
-  const pIminJ = Math.round((cavity.minY - o.oy) / o.dy);
-  const pImaxJ = Math.round((cavity.maxY - o.oy) / o.dy) - 1;
-  const kLo = Math.max(0, Math.floor((floorZ - o.oz) / o.dz) + 1);
-  const kHi = Math.min(o.gz - 1, Math.ceil((ceilZ - o.oz) / o.dz) - 1);
-  const openings = [];
-  const wallOpen = (i, j) => {
-    if (i < 0 || j < 0 || i >= o.gw || j >= o.gh)
-      return true;
-    return hasGap(o, i, j, kLo, kHi);
-  };
-  const scans = [];
-  for (const [edge, wi] of [
-    ["minX", pIminI - 1],
-    ["maxX", pImaxI + 1]
-  ]) {
-    const cells = [];
-    for (let j = pIminJ; j <= pImaxJ; j++) {
-      cells.push({ open: wallOpen(wi, j), x: o.ox + (wi + 0.5) * o.dx, y: o.oy + (j + 0.5) * o.dy, i: wi, j });
-    }
-    scans.push({ edge, cells });
-  }
-  for (const [edge, wj] of [
-    ["minY", pIminJ - 1],
-    ["maxY", pImaxJ + 1]
-  ]) {
-    const cells = [];
-    for (let i = pIminI; i <= pImaxI; i++) {
-      cells.push({ open: wallOpen(i, wj), x: o.ox + (i + 0.5) * o.dx, y: o.oy + (wj + 0.5) * o.dy, i, j: wj });
-    }
-    scans.push({ edge, cells });
-  }
-  for (const scan of scans) {
-    let run = [];
-    const flush = () => {
-      if (run.length >= MIN_OPENING_CELLS) {
-        const xs = run.map((p) => p.x);
-        const ys = run.map((p) => p.y);
-        const horiz = scan.edge === "minY" || scan.edge === "maxY";
-        const center = {
-          x: (Math.min(...xs) + Math.max(...xs)) / 2,
-          y: (Math.min(...ys) + Math.max(...ys)) / 2
-        };
-        const width = horiz ? Math.max(...xs) - Math.min(...xs) + o.dx : Math.max(...ys) - Math.min(...ys) + o.dy;
-        const mid = run[Math.floor(run.length / 2)];
-        const { zMin, zMax } = openingZSpan(o, mid.i, mid.j, kLo, kHi);
-        openings.push({
-          edge: scan.edge,
-          center: { x: round22(center.x), y: round22(center.y) },
-          width: round22(width),
-          zMin: round22(zMin),
-          zMax: round22(zMax)
-        });
-      }
-      run = [];
-    };
-    for (const c of scan.cells) {
-      if (c.open)
-        run.push({ x: c.x, y: c.y, i: c.i, j: c.j });
-      else
-        flush();
-    }
-    flush();
-  }
-  return openings;
-}
-function openingZSpan(o, i, j, kLo, kHi) {
-  let zMin = Infinity;
-  let zMax = -Infinity;
-  for (let k = kLo; k <= kHi; k++) {
-    if (o.occ[i + o.gw * (j + o.gh * k)] !== 1) {
-      const z = o.oz + (k + 0.5) * o.dz;
-      if (z < zMin)
-        zMin = z;
-      if (z > zMax)
-        zMax = z;
-    }
-  }
-  if (!Number.isFinite(zMin))
-    return { zMin: o.oz, zMax: o.maxZ };
-  return { zMin, zMax };
+  const json = enclosureWasm().enclosure_features(Float64Array.from(positions), Uint32Array.from(indices));
+  return JSON.parse(json);
 }
-var GRID_XY, GRID_Z, INSIDE, MIN_POST_HEIGHT, MIN_OPENING_CELLS, JX, JY, JZ, round22;
 var init_enclosure_mesh = __esm({
   "../engine/dist/enclosure-mesh.js"() {
     "use strict";
-    GRID_XY = 48;
-    GRID_Z = 28;
-    INSIDE = 0.5;
-    MIN_POST_HEIGHT = 0.8;
-    MIN_OPENING_CELLS = 2;
-    JX = 0.5 + 0.137;
-    JY = 0.5 - 0.077;
-    JZ = 0.5 + 0.041;
-    round22 = (n) => Math.round(n * 100) / 100;
+    init_enclosure_fit();
   }
 });
 
@@ -11955,6 +11469,7 @@ var init_dist = __esm({
           antennaAnalyze: wasmModule7.antennaAnalyze,
           photonicsSimulate: wasmModule7.photonicsSimulate,
           neutronicsSimulate: wasmModule7.neutronicsSimulate,
+          latticeGaugeSimulate: wasmModule7.latticeGaugeSimulate,
           analyzeStaticsBox: wasmModule7.analyzeStaticsBox,
           analyzeStaticsMesh: wasmModule7.analyzeStaticsMesh
         }, compiledWasmModule);
@@ -12315,6 +11830,19 @@ var init_dist = __esm({
         }
         return fn(specJson, paramsJson);
       }
+      /**
+       * Lattice gauge theory Monte Carlo (quenched SU(2)/SU(3) Wilson
+       * action): plaquette, Wilson loops, string tension, Polyakov order
+       * parameter, flux-tube profile, field snapshots — jackknife errors
+       * throughout; see `vcad-kernel-qcd`.
+       */
+      latticeGaugeSimulate(specJson) {
+        const fn = this.kernel.latticeGaugeSimulate;
+        if (typeof fn !== "function") {
+          throw new Error("latticeGaugeSimulate is not exported by this kernel WASM build \u2014 rebuild packages/kernel-wasm");
+        }
+        return fn(specJson);
+      }
       topologyOptimizeBox(min, max, spec) {
         const fn = this.kernel.topologyOptimizeBox;
         if (typeof fn !== "function") {
@@ -32264,13 +31792,13 @@ var init_document_store = __esm({
           const ang = (fp.rotation ?? 0) * Math.PI / 180;
           const cos = Math.cos(ang);
           const sin = Math.sin(ang);
-          const padWorld3 = fp.pads.map((p) => ({
+          const padWorld2 = fp.pads.map((p) => ({
             x: fp.position.x + (p.position.x * cos - p.position.y * sin),
             y: fp.position.y + (p.position.x * sin + p.position.y * cos)
           }));
           for (const tr of pcb.traces) {
             for (const end of [tr.start, tr.end]) {
-              if (padWorld3.some((pw) => Math.abs(end.x - pw.x) < TOL && Math.abs(end.y - pw.y) < TOL)) {
+              if (padWorld2.some((pw) => Math.abs(end.x - pw.x) < TOL && Math.abs(end.y - pw.y) < TOL)) {
                 end.x += dx;
                 end.y += dy;
               }
@@ -41281,6 +40809,21 @@ var init_CHANGELOG = __esm({
             "validate_for_fab"
           ]
         },
+        {
+          id: "2026-07-24-lattice-gauge-mcp",
+          version: "0.9.4",
+          date: "2026-07-24",
+          category: "feat",
+          title: "Lattice gauge theory simulation",
+          summary: "New simulate_lattice_gauge tool: quenched SU(2)/SU(3) Monte Carlo \u2014 confinement, string tension, deconfinement, flux tubes and instantons, with jackknife error bars on every number.",
+          features: [
+            "simulation",
+            "physics"
+          ],
+          mcpTools: [
+            "simulate_lattice_gauge"
+          ]
+        },
         {
           id: "2026-07-23-router-post-route-legalization",
           version: "0.9.4",
@@ -52309,6 +51852,10 @@ var init_tool_metadata = __esm({
         title: "Simulate Neutron Shield",
         annotations: RO
       },
+      simulate_lattice_gauge: {
+        title: "Simulate Lattice Gauge Theory",
+        annotations: RO
+      },
       analyze_tolerance_stackup: {
         title: "Analyze Tolerance Stackup",
         annotations: RO
@@ -58340,7 +57887,7 @@ function computeUtilization(footprints, vertices, cutouts, shape, innerRadius, e
   }
   const boardArea = Math.abs(loopSignedArea(vertices)) - (cutouts ?? []).reduce((s, c) => s + Math.abs(loopSignedArea(c)), 0);
   if (!Number.isFinite(occMinX) || boardArea <= 0) return void 0;
-  const round25 = (v) => Math.round(v * 100) / 100;
+  const round23 = (v) => Math.round(v * 100) / 100;
   const ceilHalf = (v) => Math.ceil(v * 2) / 2;
   const bbW = occMaxX - occMinX;
   const bbH = occMaxY - occMinY;
@@ -58353,8 +57900,8 @@ function computeUtilization(footprints, vertices, cutouts, shape, innerRadius, e
     suggested = {
       type: "circle",
       outer_diameter: od,
-      center: { x: round25(cx), y: round25(cy) },
-      ...innerRadius > 0 ? { inner_diameter: round25(2 * innerRadius) } : {},
+      center: { x: round23(cx), y: round23(cy) },
+      ...innerRadius > 0 ? { inner_diameter: round23(2 * innerRadius) } : {},
       note: `Minimum enclosing circle with ${margin}mm edge clearance`
     };
   } else {
@@ -58362,19 +57909,19 @@ function computeUtilization(footprints, vertices, cutouts, shape, innerRadius, e
       type: "rect",
       width: ceilHalf(bbW + 2 * margin),
       height: ceilHalf(bbH + 2 * margin),
-      origin: { x: round25(occMinX - margin), y: round25(occMinY - margin) },
+      origin: { x: round23(occMinX - margin), y: round23(occMinY - margin) },
       note: `Minimum enclosing rectangle with ${margin}mm edge clearance`
     };
   }
   return {
-    board_area_mm2: round25(boardArea),
-    component_area_mm2: round25(componentArea),
+    board_area_mm2: round23(boardArea),
+    component_area_mm2: round23(componentArea),
     utilization_pct: Math.round(componentArea / boardArea * 1e3) / 10,
-    bounding_box: { x: round25(occMinX), y: round25(occMinY), w: round25(bbW), h: round25(bbH) },
+    bounding_box: { x: round23(occMinX), y: round23(occMinY), w: round23(bbW), h: round23(bbH) },
     suggested_outline: suggested
   };
 }
-function padWorld2(fp, pad) {
+function padWorld(fp, pad) {
   const ang = (fp.rotation ?? 0) * Math.PI / 180;
   const cos = Math.cos(ang);
   const sin = Math.sin(ang);
@@ -58393,7 +57940,7 @@ function detectStaleNets(pcb) {
       const arr = padsByNet.get(pad.net) ?? [];
       const sh = pad.shape;
       const halfExtent = sh.diameter ? sh.diameter / 2 : Math.hypot(sh.width ?? 0, sh.height ?? 0) / 2;
-      arr.push({ pos: padWorld2(fp, pad), layers: pad.layers, halfExtent });
+      arr.push({ pos: padWorld(fp, pad), layers: pad.layers, halfExtent });
       padsByNet.set(pad.net, arr);
     }
   }
@@ -59277,7 +58824,7 @@ function layoutLint(pcb) {
   for (const fp of fps) {
     const prefix = refPrefix(fp.ref);
     for (const pad of fp.pads) {
-      const w = { ref: fp.ref, prefix, pad, pos: padWorld2(fp, pad) };
+      const w = { ref: fp.ref, prefix, pad, pos: padWorld(fp, pad) };
       allPads.push(w);
       if (pad.net) {
         const arr = byNet.get(pad.net) ?? [];
@@ -59293,7 +58840,7 @@ function layoutLint(pcb) {
     for (const pad of fp.pads) {
       const net = pad.net;
       if (!net || LINT_GROUND_NET_RE.test(net) || LINT_SUPPLY_NET_RE.test(net)) continue;
-      const from = padWorld2(fp, pad);
+      const from = padWorld(fp, pad);
       let nearest = null;
       for (const other of byNet.get(net) ?? []) {
         if (other.ref === fp.ref || other.prefix !== "U") continue;
@@ -59308,9 +58855,9 @@ function layoutLint(pcb) {
       warnings.push({
         kind: "crystal_far_from_ic",
         refs: [fp.ref, worst2.icRef],
-        distance_mm: round23(worst2.dist),
+        distance_mm: round2(worst2.dist),
         threshold_mm: LINT_CRYSTAL_MAX_MM,
-        message: `Crystal ${fp.ref} is ${round23(worst2.dist)}mm from ${worst2.icRef} pad ${worst2.icPad} (net '${worst2.net}') \u2014 oscillator loops want <${LINT_CRYSTAL_MAX_MM}mm; move ${fp.ref} next to its XIN/XOUT pins.`
+        message: `Crystal ${fp.ref} is ${round2(worst2.dist)}mm from ${worst2.icRef} pad ${worst2.icPad} (net '${worst2.net}') \u2014 oscillator loops want <${LINT_CRYSTAL_MAX_MM}mm; move ${fp.ref} next to its XIN/XOUT pins.`
       });
     }
   }
@@ -59321,7 +58868,7 @@ function layoutLint(pcb) {
     const pwrIdx = nets.findIndex((n) => n && LINT_SUPPLY_NET_RE.test(n) && !LINT_GROUND_NET_RE.test(n));
     if (gndIdx < 0 || pwrIdx < 0) continue;
     const pwrPad = fp.pads[pwrIdx];
-    const from = padWorld2(fp, pwrPad);
+    const from = padWorld(fp, pwrPad);
     let nearest = null;
     for (const other of byNet.get(pwrPad.net) ?? []) {
       if (other.ref === fp.ref || other.prefix !== "U") continue;
@@ -59334,9 +58881,9 @@ function layoutLint(pcb) {
       warnings.push({
         kind: "decoupling_cap_far_from_pin",
         refs: [fp.ref, nearest.icRef],
-        distance_mm: round23(nearest.dist),
+        distance_mm: round2(nearest.dist),
         threshold_mm: LINT_DECAP_MAX_MM,
-        message: `Decoupling cap ${fp.ref} is ${round23(nearest.dist)}mm from ${nearest.icRef} pad ${nearest.icPad} (net '${pwrPad.net}') \u2014 decouplers want <${LINT_DECAP_MAX_MM}mm to the pin they decouple.`
+        message: `Decoupling cap ${fp.ref} is ${round2(nearest.dist)}mm from ${nearest.icRef} pad ${nearest.icPad} (net '${pwrPad.net}') \u2014 decouplers want <${LINT_DECAP_MAX_MM}mm to the pin they decouple.`
       });
     }
   }
@@ -59354,15 +58901,15 @@ function layoutLint(pcb) {
     for (const fp of fps) {
       const prefix = refPrefix(fp.ref);
       if (prefix !== "J" && prefix !== "P" && prefix !== "CN" && prefix !== "USB") continue;
-      const pts = fp.pads.length > 0 ? fp.pads.map((p) => padWorld2(fp, p)) : [fp.position];
+      const pts = fp.pads.length > 0 ? fp.pads.map((p) => padWorld(fp, p)) : [fp.position];
       const d = Math.min(...pts.map(edgeDist));
       if (d > LINT_CONNECTOR_EDGE_MAX_MM) {
         warnings.push({
           kind: "connector_not_on_edge",
           refs: [fp.ref],
-          distance_mm: round23(d),
+          distance_mm: round2(d),
           threshold_mm: LINT_CONNECTOR_EDGE_MAX_MM,
-          message: `Connector ${fp.ref} sits ${round23(d)}mm from the board edge \u2014 connectors want edge access (<${LINT_CONNECTOR_EDGE_MAX_MM}mm); move it to the outline.`
+          message: `Connector ${fp.ref} sits ${round2(d)}mm from the board edge \u2014 connectors want edge access (<${LINT_CONNECTOR_EDGE_MAX_MM}mm); move it to the outline.`
         });
       }
     }
@@ -59398,9 +58945,9 @@ function layoutLint(pcb) {
       warnings.push({
         kind: "high_current_near_sensitive",
         refs: [c.hi.ref, c.lo.ref],
-        distance_mm: round23(c.dist),
+        distance_mm: round2(c.dist),
         threshold_mm: LINT_HIGH_CURRENT_SEP_MIN_MM,
-        message: `High-current net '${c.hiNet}' (${c.hi.ref} pad ${c.hi.pad.number}) is ${round23(c.dist)}mm from '${c.loNet}' (${c.lo.ref} pad ${c.lo.pad.number}) \u2014 keep \u2265${LINT_HIGH_CURRENT_SEP_MIN_MM}mm between high-current and USB/analog pads.`
+        message: `High-current net '${c.hiNet}' (${c.hi.ref} pad ${c.hi.pad.number}) is ${round2(c.dist)}mm from '${c.loNet}' (${c.lo.ref} pad ${c.lo.pad.number}) \u2014 keep \u2265${LINT_HIGH_CURRENT_SEP_MIN_MM}mm between high-current and USB/analog pads.`
       });
     }
   }
@@ -63009,7 +62556,7 @@ async function addMotorWinding(args) {
       if (!pad.net) continue;
       const arr = padsByNet.get(pad.net) ?? [];
       arr.push({
-        pos: padWorld2(fp, pad),
+        pos: padWorld(fp, pad),
         onReturnLayer: pad.layers.includes(returnLayer),
         size: padApproxRadius(pad)
       });
@@ -63352,7 +62899,7 @@ async function addMotorWinding(args) {
           for (const fp of pcb.footprints) {
             for (const pad of fp.pads) {
               if (pad.net === net) continue;
-              const pw = padWorld2(fp, pad);
+              const pw = padWorld(fp, pad);
               if (pointSegDist(pw, t2.a, t2.b) < padApproxRadius(pad) + halfW + clearance) return false;
             }
           }
@@ -65045,7 +64592,7 @@ async function fixDrc(args) {
     let corridorFail = null;
     for (const pt of movePoints) {
       const onPad = pcb.footprints.some(
-        (fp) => fp.pads.some((pad) => samePoint(padWorld2(fp, pad), pt))
+        (fp) => fp.pads.some((pad) => samePoint(padWorld(fp, pad), pt))
       );
       if (onPad) {
         corridorFail = "endpoint lands on a component pad \u2014 fixing requires moving the component";
@@ -65195,7 +64742,7 @@ async function fixDrc(args) {
   };
   return { content: [{ type: "text", text: JSON.stringify(payload) }] };
 }
-var PCB_LAYER_SET, COPPER_LAYERS, COPPER_LAYER_SET, POSITION_TOLERANCE, PIN_SEP, pinKey, createSchematicSchema, docInputProperties, placeComponentsSchema, routeNetsSchema, runDrcSchema, runErcSchema, critiqueRouteSchema, routeDiffPairSchema, lengthMatchTracesSchema, exportGerberSchema, validateForFabSchema, exportKicadSchema, addCoilSchema, boardFromSolidSchema, addCoilArraySchema, windingLayoutSchema, calcImpedanceSchema, sizeImpedanceSchema, sizePdnSchema, calcCoilSchema, sizeCoilSchema, calcRfSchema, COIL_GEOMETRY, MU0, RAIL_PIN_NAMES, PIN_TYPES, CHIP_SIZE_CODES, CONNECTIVITY_GUARD_MAX_NETS, POWER_NET_RE, DRC_CATEGORY, DRC_DELTA_SAMPLE_CAP, DRC_DELTA_FULL_BOARD_MAX, LINT_CRYSTAL_MAX_MM, LINT_DECAP_MAX_MM, LINT_CONNECTOR_EDGE_MAX_MM, LINT_HIGH_CURRENT_SEP_MIN_MM, LINT_GROUND_NET_RE, LINT_SUPPLY_NET_RE, SENSITIVE_NET_RE, HIGH_CURRENT_CLASS_RE, refPrefix, dist2d, round23, sumSq, round33, vec2Schema, setBoardOutlineSchema, FOOTPRINT_FAMILIES, listFootprintsSchema, searchFootprintsSchema, getPadPositionsSchema, getFootprintSchema, ROTATION_CONVENTION, describePcbSchema, DESCRIBE_NET_NAME_CAP, addTraceSchema, addViaSchema, OZ_TO_MM, setStackupSchema, setPlacementSchema, addZoneSchema, deleteElementProps, deleteZoneSchema, deleteTraceSchema, deleteViaSchema, COPPER_KINDS, GET_COPPER_CAP, getCopperSchema, TIE_POSITION_TOL, addNetTieSchema, deleteNetTieSchema, undoSchema, setDesignRulesSchema, sizeTraceForCurrentSchema, addViaArraySchema, MAX_VIA_ARRAY, addMotorWindingSchema, calcMotorSchema, sig6, BEARING_FRICTION_MNM, checkSelfStartSchema, solidFromBoardSchema, FR4_DENSITY_KG_M3, COPPER_DENSITY_KG_M3, DEFAULT_COPPER_THICKNESS_MM, loopArea, MOUNTING_FP_RE, KEEPOUT_ECHO_CAP, INLINE_DOC_ECHO_CAP, searchElectronicPartsSchema, resolvePartSchema, findAlternativesSchema, verifySubstitutionSchema, buildReceiptSchema, verifyReceiptSchema, FIX_EPS, FIX_DRC_NEVER, fixDrcSchema, toolDefs20;
+var PCB_LAYER_SET, COPPER_LAYERS, COPPER_LAYER_SET, POSITION_TOLERANCE, PIN_SEP, pinKey, createSchematicSchema, docInputProperties, placeComponentsSchema, routeNetsSchema, runDrcSchema, runErcSchema, critiqueRouteSchema, routeDiffPairSchema, lengthMatchTracesSchema, exportGerberSchema, validateForFabSchema, exportKicadSchema, addCoilSchema, boardFromSolidSchema, addCoilArraySchema, windingLayoutSchema, calcImpedanceSchema, sizeImpedanceSchema, sizePdnSchema, calcCoilSchema, sizeCoilSchema, calcRfSchema, COIL_GEOMETRY, MU0, RAIL_PIN_NAMES, PIN_TYPES, CHIP_SIZE_CODES, CONNECTIVITY_GUARD_MAX_NETS, POWER_NET_RE, DRC_CATEGORY, DRC_DELTA_SAMPLE_CAP, DRC_DELTA_FULL_BOARD_MAX, LINT_CRYSTAL_MAX_MM, LINT_DECAP_MAX_MM, LINT_CONNECTOR_EDGE_MAX_MM, LINT_HIGH_CURRENT_SEP_MIN_MM, LINT_GROUND_NET_RE, LINT_SUPPLY_NET_RE, SENSITIVE_NET_RE, HIGH_CURRENT_CLASS_RE, refPrefix, dist2d, round2, sumSq, round33, vec2Schema, setBoardOutlineSchema, FOOTPRINT_FAMILIES, listFootprintsSchema, searchFootprintsSchema, getPadPositionsSchema, getFootprintSchema, ROTATION_CONVENTION, describePcbSchema, DESCRIBE_NET_NAME_CAP, addTraceSchema, addViaSchema, OZ_TO_MM, setStackupSchema, setPlacementSchema, addZoneSchema, deleteElementProps, deleteZoneSchema, deleteTraceSchema, deleteViaSchema, COPPER_KINDS, GET_COPPER_CAP, getCopperSchema, TIE_POSITION_TOL, addNetTieSchema, deleteNetTieSchema, undoSchema, setDesignRulesSchema, sizeTraceForCurrentSchema, addViaArraySchema, MAX_VIA_ARRAY, addMotorWindingSchema, calcMotorSchema, sig6, BEARING_FRICTION_MNM, checkSelfStartSchema, solidFromBoardSchema, FR4_DENSITY_KG_M3, COPPER_DENSITY_KG_M3, DEFAULT_COPPER_THICKNESS_MM, loopArea, MOUNTING_FP_RE, KEEPOUT_ECHO_CAP, INLINE_DOC_ECHO_CAP, searchElectronicPartsSchema, resolvePartSchema, findAlternativesSchema, verifySubstitutionSchema, buildReceiptSchema, verifyReceiptSchema, FIX_EPS, FIX_DRC_NEVER, fixDrcSchema, toolDefs20;
 var init_ecad2 = __esm({
   "src/tools/ecad.ts"() {
     "use strict";
@@ -65997,7 +65544,7 @@ var init_ecad2 = __esm({
     HIGH_CURRENT_CLASS_RE = /pwr|power|high[-_ ]?current|motor|\bhv\b|hi[-_ ]?amp/i;
     refPrefix = (ref) => /^[A-Za-z]+/.exec(ref)?.[0]?.toUpperCase() ?? "";
     dist2d = (a, b) => Math.hypot(a.x - b.x, a.y - b.y);
-    round23 = (n) => Math.round(n * 100) / 100;
+    round2 = (n) => Math.round(n * 100) / 100;
     sumSq = (v) => v.reduce((s, x) => s + x * x, 0);
     round33 = (v) => Math.round(v * 1e3) / 1e3;
     vec2Schema = {
@@ -71074,6 +70621,128 @@ var init_neutronics = __esm({
   }
 });
 
+// src/tools/qcd.ts
+function textResult11(payload) {
+  return {
+    content: [
+      {
+        type: "text",
+        text: JSON.stringify(payload, null, 2)
+      }
+    ]
+  };
+}
+var simSpecSchema, toolDefs36;
+var init_qcd = __esm({
+  "src/tools/qcd.ts"() {
+    "use strict";
+    init_tool_def();
+    simSpecSchema = {
+      type: "object",
+      required: [
+        "dims",
+        "beta",
+        "thermalization_sweeps",
+        "measurement_sweeps",
+        "overrelax_per_heatbath",
+        "bin_size",
+        "max_wilson_extent",
+        "seed",
+        "hot_start"
+      ],
+      description: "Lattice simulation spec. Deterministic per seed. Direction 3 is time \u2014 shorten it (e.g. [6,6,6,2]) for finite temperature / deconfinement studies.",
+      properties: {
+        gauge: {
+          type: "string",
+          enum: ["Su2", "Su3"],
+          description: "Gauge group. Su2 (default): quaternion links, exact Kennedy-Pendleton heatbath, fastest. Su3: the real QCD gauge group, Cabibbo-Marinari updates, ~6x cost per link."
+        },
+        dims: {
+          type: "array",
+          items: { type: "number" },
+          minItems: 4,
+          maxItems: 4,
+          description: "Lattice extents [n0,n1,n2,n3], each >= 2. 4^4-8^4 are seconds-scale; 12^4+ is where the cost gate starts to bite."
+        },
+        beta: {
+          type: "number",
+          description: "Inverse coupling beta = 2N/g^2 (Wilson action). SU(2): strong coupling < ~1.5, crossover ~2.2-2.5, weak > 4. SU(3): strong < ~4, deconfinement at N_t=2 near 5.1, weak > 7."
+        },
+        thermalization_sweeps: {
+          type: "number",
+          description: "Discarded equilibration sweeps (>= 1). 50 is typical."
+        },
+        measurement_sweeps: {
+          type: "number",
+          description: "Measured sweeps. Must give >= 2 full jackknife bins (>= 5 bins to mint claims). 100-200 is typical."
+        },
+        overrelax_per_heatbath: {
+          type: "number",
+          description: "Overrelaxation sweeps interleaved per heatbath sweep (1-3 decorrelates much faster per unit work)."
+        },
+        bin_size: {
+          type: "number",
+          description: "Jackknife bin size in measurements (10-20 typical)."
+        },
+        max_wilson_extent: {
+          type: "number",
+          description: "Largest square Wilson loop W(r,t) to measure; 0 = plaquette only. Capped at half the smallest lattice extent."
+        },
+        seed: { type: "number" },
+        hot_start: {
+          type: "boolean",
+          description: "Start from random links instead of cold (identity)."
+        },
+        smear: {
+          type: "object",
+          required: ["alpha", "iterations"],
+          description: "APE spatial smearing applied to a measurement copy before loop observables (lifts the string-tension signal at intermediate coupling; plaquette stays unsmeared). alpha 0.5, 2-3 iterations is customary.",
+          properties: {
+            alpha: { type: "number" },
+            iterations: { type: "number" }
+          }
+        },
+        measure_temporal_loops: {
+          type: "boolean",
+          description: "Also measure spatial x temporal Wilson loops and derive the static quark potential V(r) (+ Cornell fit when >= 3 points resolve)."
+        },
+        measure_polyakov: {
+          type: "boolean",
+          description: "Measure <|L|>, the Polyakov-loop deconfinement order parameter (small = confined, O(1) = deconfined)."
+        },
+        flux_tube: {
+          type: "object",
+          required: ["separation"],
+          description: "Measure the chromoelectric flux-tube profile between a static quark-antiquark pair at this separation (lattice units, spatial axis 0): the connected (Polyakov-pair x action-density) 3D excess field plus the pair correlator whose decay with separation IS confinement. Use a short time axis (N_t 2-4).",
+          properties: { separation: { type: "number" } }
+        },
+        snapshot_cooling: {
+          type: "number",
+          description: "Export a rendering FieldSnapshot (action density per site + complex Polyakov field per spatial site) of the final configuration after this many cooling sweeps (0 = raw vacuum boil; 20-30 = classical lumps). SU(2) runs with cooling >= 1 also report the clover topological charge (near-integer = instanton content)."
+        }
+      }
+    };
+    toolDefs36 = [
+      {
+        name: "simulate_lattice_gauge",
+        pack: null,
+        description: "Lattice gauge theory Monte Carlo (quenched Wilson action, SU(2) or SU(3)): confinement from first principles on a laptop-scale 4D lattice. Returns the average plaquette, Wilson loops, Creutz ratios + static potential + Cornell fit (the string tension sigma*a^2), the Polyakov deconfinement order parameter, and optionally viewport-ready field snapshots (action density, Polyakov field, cooled topological charge) and the static-pair flux-tube profile \u2014 every observable a binned-jackknife mean +- error, bit-deterministic per seed. Fail-closed: starved statistics refuse to mint vcad.qcd-claims/1 (raise measurement_sweeps); logs require loops resolved >= 3-sigma from zero. Honesty bounds ride every claim: quenched (no dynamical fermions), lattice units at fixed coupling (no continuum limit) \u2014 model physics, not physical-QCD numbers. Cost ~ volume x sweeps (x6 for SU(3)); 6^4 with 150 sweeps runs in seconds. Recipes: confinement demo = Su2, [6,6,6,4], beta 2.2, flux_tube {separation 2}; deconfinement = [6,6,6,2], measure_polyakov, beta scan across ~1.88 (Su2) / ~5.1 (Su3); instantons = Su2, beta 2.4, snapshot_cooling 30.",
+        inputSchema: {
+          type: "object",
+          required: ["spec"],
+          properties: { spec: simSpecSchema }
+        },
+        handler: (args, ctx) => {
+          const a = args;
+          const out = ctx.engine.latticeGaugeSimulate(JSON.stringify(a.spec));
+          return textResult11(out);
+        },
+        behavior: behavior({})
+      }
+    ];
+  }
+});
+
 // src/tools/physics.ts
 function limitClaim(id, description, subject, limit, actual, unit, basis, converged) {
   if (!converged || !Number.isFinite(actual)) {
@@ -71195,7 +70864,7 @@ function predictPhysicsTool(args, engine) {
     ]
   };
 }
-var PHYSICS_DOMAIN, ORACLE, TIER_RESOLUTION, regionSchema3, predictPhysicsSchema, round52, toolDefs36;
+var PHYSICS_DOMAIN, ORACLE, TIER_RESOLUTION, regionSchema3, predictPhysicsSchema, round52, toolDefs37;
 var init_physics2 = __esm({
   "src/tools/physics.ts"() {
     "use strict";
@@ -71308,7 +70977,7 @@ var init_physics2 = __esm({
       }
     };
     round52 = (v) => Number(v.toPrecision(5));
-    toolDefs36 = [
+    toolDefs37 = [
       {
         name: "predict_physics",
         pack: null,
@@ -71550,7 +71219,7 @@ function applySetParam(doc, node, path, value) {
   }
   cursor[parts[parts.length - 1]] = value;
 }
-var lastReports, mechanicalProcessEnum, pcbProcessEnum, processEnum, dfmCheckSchema, dfmExplainSchema, dfmSuggestFixSchema, dfmApplyFixSchema, toolDefs37;
+var lastReports, mechanicalProcessEnum, pcbProcessEnum, processEnum, dfmCheckSchema, dfmExplainSchema, dfmSuggestFixSchema, dfmApplyFixSchema, toolDefs38;
 var init_dfm2 = __esm({
   "src/tools/dfm.ts"() {
     "use strict";
@@ -71647,7 +71316,7 @@ var init_dfm2 = __esm({
       },
       required: ["issue_id"]
     };
-    toolDefs37 = [
+    toolDefs38 = [
       {
         name: "dfm_check",
         pack: "dfm",
@@ -71846,7 +71515,7 @@ function renderedOf(engine, doc) {
     "document has no sheet-metal part (evaluation produced no sheetMetal bundle)"
   );
 }
-function textResult11(payload) {
+function textResult12(payload) {
   return {
     content: [{ type: "text", text: JSON.stringify(payload, null, 2) }]
   };
@@ -71888,7 +71557,7 @@ function sheetMetalCreate(input, engine) {
   const errors = rendered.violations.filter(
     (v) => v.severity === "Error"
   ).length;
-  return textResult11({
+  return textResult12({
     document_id: documentId,
     model: rendered.model,
     flat: {
@@ -71904,7 +71573,7 @@ function sheetMetalUnfold(input, engine) {
   const doc = getSession(String(a.document_id ?? ""));
   const rendered = renderedOf(engine, doc);
   const includeDxf = a.include_dxf !== false;
-  return textResult11({
+  return textResult12({
     flat_pattern: rendered.flatPattern,
     ...includeDxf ? { dxf: rendered.dxf } : {},
     note: "The DXF is a fab-ready merged silhouette: one closed exterior LWPOLYLINE plus hole loops on CUT, DASHED bend centerlines (on the allowance midline) on BEND_UP/BEND_DOWN, and any surface markings as open polylines on ENGRAVE (select the laser-marking/engraving service for that layer in the fab's UI). DXF carries no bend angles \u2014 you enter those in the fab's UI. For zero data entry, export the folded body as STEP instead (export_cad with a .step filename); bends are auto-detected by the shop, but radii/K must match their tooling at model time (use shop_profile in sheet_metal_create)."
@@ -71912,7 +71581,7 @@ function sheetMetalUnfold(input, engine) {
 }
 function sheetMetalMaterials(_input, engine) {
   const materials = engine.getSheetMetalMaterials();
-  return textResult11({
+  return textResult12({
     materials,
     note: 'Pass any `name` (e.g. `"al-soft"`, `"steel-mild"`) as the `material` field of `sheet_metal_create`. Aliases like `"aluminum"`, `"stainless"`, `"6061-T6"` also resolve.'
   });
@@ -71921,13 +71590,13 @@ function sheetMetalBendTable(input, engine) {
   const a = input ?? {};
   if (typeof a.shop_profile === "string" && a.shop_profile.length > 0) {
     const catalog = engine.getSheetMetalShopCatalog(a.shop_profile);
-    return textResult11({
+    return textResult12({
       catalog,
       note: `Published catalog for ${catalog.name} (retrieved ${catalog.retrieved}). inside_radius_mm is FIXED per material/thickness \u2014 pass shop_profile: "${catalog.id}" to sheet_metal_create and omit per-bend radius to use it.`
     });
   }
   const table = engine.getSheetMetalBendTable();
-  return textResult11({
+  return textResult12({
     table,
     note: "K-factors are interpolated by closest R/t for the chosen material. To override K for a specific bend, pass `manual_k` on that flange in sheet_metal_create."
   });
@@ -71949,7 +71618,7 @@ function sheetMetalCheck(input, engine) {
   const errors = result.violations.filter(
     (v) => v.severity === "Error"
   ).length;
-  return textResult11({
+  return textResult12({
     shop: result.shop,
     violations: result.violations,
     shop_ready: result.violations.length === 0,
@@ -72001,7 +71670,7 @@ function sheetMetalNest(input, engine) {
     footprints,
     params
   );
-  return textResult11({
+  return textResult12({
     parts: footprints,
     result,
     summary: `${result.placements.length} placements across ${result.sheets_used} sheet(s); ${result.utilization_pct.toFixed(1)}% utilization${result.unplaceable.length > 0 ? `; ${result.unplaceable.length} unplaceable (oversize)` : ""}.`
@@ -72014,7 +71683,7 @@ function sheetMetalSequence3(input, engine) {
   if (!steps) {
     throw new Error("document has no sheet-metal part");
   }
-  return textResult11({
+  return textResult12({
     count: steps.length,
     steps,
     note: "Form each bend in order. compensated_angle_rad is the angle to dial in on the brake \u2014 springs back to angle_rad once released."
@@ -72037,7 +71706,7 @@ function sheetMetalSuggestFix(input, engine) {
     violation: v,
     fix: suggestFix(v)
   }));
-  return textResult11({
+  return textResult12({
     shop: result.shop,
     count: suggestions.length,
     suggestions,
@@ -72140,13 +71809,13 @@ function sheetMetalCost(input, engine) {
   if (!result) {
     throw new Error("document has no sheet-metal part");
   }
-  return textResult11({
+  return textResult12({
     breakdown: result.breakdown,
     rates: result.rates,
     summary: `${result.breakdown.currency} ${result.breakdown.total_each.toFixed(2)} each @ qty ${result.breakdown.quantity} (mass ${result.breakdown.mass_kg_each.toFixed(3)} kg, ${result.breakdown.cut_length_m.toFixed(2)} m cut, ${result.breakdown.bends} bend(s)).`
   });
 }
-var sheetMetalCreateSchema, sheetMetalUnfoldSchema, sheetMetalMaterialsSchema, sheetMetalBendTableSchema, sheetMetalCheckSchema, sheetMetalNestSchema, sheetMetalSequenceSchema, sheetMetalSuggestFixSchema, sheetMetalCostSchema, toolDefs38;
+var sheetMetalCreateSchema, sheetMetalUnfoldSchema, sheetMetalMaterialsSchema, sheetMetalBendTableSchema, sheetMetalCheckSchema, sheetMetalNestSchema, sheetMetalSequenceSchema, sheetMetalSuggestFixSchema, sheetMetalCostSchema, toolDefs39;
 var init_sheet_metal2 = __esm({
   "src/tools/sheet-metal.ts"() {
     "use strict";
@@ -72407,7 +72076,7 @@ var init_sheet_metal2 = __esm({
       },
       required: ["document_id"]
     };
-    toolDefs38 = [
+    toolDefs39 = [
       {
         name: "sheet_metal_create",
         pack: "sheet_metal",
@@ -72618,12 +72287,12 @@ function simulateStrike(input, engine) {
     },
     modes: result.modes.map((m) => ({
       n: m.n,
-      hz: round24(m.hz),
+      hz: round22(m.hz),
       gain: round4(m.gain),
       q: Math.round(m.q),
-      t60_s: round24(m.t60_s)
+      t60_s: round22(m.t60_s)
     })),
-    spectrum_peaks: result.spectrum_peaks.map((p) => ({ hz: round24(p.hz), db: round24(p.db) })),
+    spectrum_peaks: result.spectrum_peaks.map((p) => ({ hz: round22(p.hz), db: round22(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"
   };
@@ -72652,7 +72321,7 @@ function simulateStrike(input, engine) {
   }
   return { content: [{ type: "text", text: JSON.stringify(payload, null, 2) }] };
 }
-var cents, MATERIAL_ALIASES, simulateStrikeSchema, round24, round4, toolDefs39;
+var cents, MATERIAL_ALIASES, simulateStrikeSchema, round22, round4, toolDefs40;
 var init_acoustics = __esm({
   "src/tools/acoustics.ts"() {
     "use strict";
@@ -72724,9 +72393,9 @@ var init_acoustics = __esm({
         }
       }
     };
-    round24 = (x) => Math.round(x * 100) / 100;
+    round22 = (x) => Math.round(x * 100) / 100;
     round4 = (x) => Math.round(x * 1e4) / 1e4;
-    toolDefs39 = [
+    toolDefs40 = [
       {
         name: "simulate_strike",
         pack: "sheet_metal",
@@ -72877,7 +72546,7 @@ function importStep(input, engine) {
     ]
   };
 }
-var MAX_STEP_BYTES, importStepSchema, toolDefs40;
+var MAX_STEP_BYTES, importStepSchema, toolDefs41;
 var init_import = __esm({
   "src/tools/import.ts"() {
     "use strict";
@@ -72909,7 +72578,7 @@ var init_import = __esm({
         }
       }
     };
-    toolDefs40 = [
+    toolDefs41 = [
       {
         name: "import_step",
         pack: null,
@@ -73037,7 +72706,7 @@ function importEagle(_input) {
     "Eagle (.brd) import is not yet supported. Export your board from Eagle as KiCad (File > Export > KiCad .kicad_pcb) and use import_kicad instead."
   );
 }
-var MAX_PCB_BYTES, importKicadSchema, importEagleSchema, toolDefs41;
+var MAX_PCB_BYTES, importKicadSchema, importEagleSchema, toolDefs42;
 var init_import_pcb = __esm({
   "src/tools/import-pcb.ts"() {
     "use strict";
@@ -73078,7 +72747,7 @@ var init_import_pcb = __esm({
         }
       }
     };
-    toolDefs41 = [
+    toolDefs42 = [
       {
         name: "import_kicad",
         pack: "ecad",
@@ -73128,8 +72797,8 @@ function parseDocument(input) {
   }
   throw new Error("Document must be JSON (starting with '{') or VCode (starting with '#')");
 }
-function compressForUrl(compact) {
-  const compressed = gzipSync(Buffer.from(compact, "utf-8"), { level: 9 });
+function compressForUrl(compact2) {
+  const compressed = gzipSync(Buffer.from(compact2, "utf-8"), { level: 9 });
   return base64urlEncode(compressed);
 }
 function openInBrowser(input) {
@@ -73174,7 +72843,7 @@ Compressed URL param: ${encoded.length} bytes`
     ]
   };
 }
-var openInBrowserSchema, toolDefs42;
+var openInBrowserSchema, toolDefs43;
 var init_share = __esm({
   "src/tools/share.ts"() {
     "use strict";
@@ -73198,7 +72867,7 @@ var init_share = __esm({
         }
       }
     };
-    toolDefs42 = [
+    toolDefs43 = [
       {
         name: "open_in_browser",
         pack: null,
@@ -73523,7 +73192,7 @@ function batchReset(input) {
     };
   }
 }
-var simulations, MAX_TRAJECTORY, envRecords, batchGroups, createRobotEnvSchema, gymStepSchema, gymResetSchema, gymObserveSchema, gymCloseSchema, batchCreateEnvsSchema, batchStepSchema, batchResetSchema, toolDefs43;
+var simulations, MAX_TRAJECTORY, envRecords, batchGroups, createRobotEnvSchema, gymStepSchema, gymResetSchema, gymObserveSchema, gymCloseSchema, batchCreateEnvsSchema, batchStepSchema, batchResetSchema, toolDefs44;
 var init_gym = __esm({
   "src/tools/gym.ts"() {
     "use strict";
@@ -73675,7 +73344,7 @@ var init_gym = __esm({
       },
       required: ["batch_id"]
     };
-    toolDefs43 = [
+    toolDefs44 = [
       {
         name: "create_robot_env",
         pack: "physics",
@@ -74033,7 +73702,7 @@ async function renderMolecule(input) {
     return fail(err7);
   }
 }
-var loadStructureSchema, inspectMoleculeSchema, configSchema, minimizeEnergySchema, mdRunSchema, designMaterialSchema, homogenizeMaterialSchema, renderMoleculeSchema, CPK, DEFAULT_RADIUS, toolDefs44;
+var loadStructureSchema, inspectMoleculeSchema, configSchema, minimizeEnergySchema, mdRunSchema, designMaterialSchema, homogenizeMaterialSchema, renderMoleculeSchema, CPK, DEFAULT_RADIUS, toolDefs45;
 var init_atoms2 = __esm({
   "src/tools/atoms.ts"() {
     "use strict";
@@ -74168,7 +73837,7 @@ var init_atoms2 = __esm({
       S: 1.05,
       P: 1.07
     };
-    toolDefs44 = [
+    toolDefs45 = [
       {
         name: "load_structure",
         pack: "atoms",
@@ -75376,7 +75045,7 @@ function errorResult(text) {
     isError: true
   };
 }
-var recordSimulationSchema, SVG_SCALE2, MAX_STEPS, MIN_STEPS, DEFAULT_WIDTH_PX2, MIN_WIDTH_PX, MAX_WIDTH_PX, DEFAULT_FPS, MIN_FPS, MAX_FPS, toolDefs45;
+var recordSimulationSchema, SVG_SCALE2, MAX_STEPS, MIN_STEPS, DEFAULT_WIDTH_PX2, MIN_WIDTH_PX, MAX_WIDTH_PX, DEFAULT_FPS, MIN_FPS, MAX_FPS, toolDefs46;
 var init_record2 = __esm({
   "src/tools/record.ts"() {
     "use strict";
@@ -75440,7 +75109,7 @@ var init_record2 = __esm({
     DEFAULT_FPS = 30;
     MIN_FPS = 1;
     MAX_FPS = 60;
-    toolDefs45 = [
+    toolDefs46 = [
       {
         name: "record_simulation",
         pack: "physics",
@@ -75533,7 +75202,7 @@ ${text}`
     ]
   };
 }
-var getChangelogSchema, toolDefs46;
+var getChangelogSchema, toolDefs47;
 var init_changelog2 = __esm({
   "src/tools/changelog.ts"() {
     "use strict";
@@ -75573,7 +75242,7 @@ var init_changelog2 = __esm({
         }
       }
     };
-    toolDefs46 = [
+    toolDefs47 = [
       {
         name: "get_changelog",
         pack: null,
@@ -75722,7 +75391,7 @@ function getSimVersion(input) {
     ]
   };
 }
-var getSimReplaySchema, getSimVersionSchema, toolDefs47;
+var getSimReplaySchema, getSimVersionSchema, toolDefs48;
 var init_sim_replay = __esm({
   "src/tools/sim-replay.ts"() {
     "use strict";
@@ -75750,7 +75419,7 @@ var init_sim_replay = __esm({
       },
       required: ["env_id"]
     };
-    toolDefs47 = [
+    toolDefs48 = [
       {
         name: "get_sim_replay",
         pack: null,
@@ -75869,7 +75538,7 @@ async function getOrderFeed(input, store, user) {
     version: version2
   });
 }
-var toUsd4, getOrderFeedSchema, toolDefs48;
+var toUsd4, getOrderFeedSchema, toolDefs49;
 var init_order_feed = __esm({
   "src/tools/order-feed.ts"() {
     "use strict";
@@ -75888,7 +75557,7 @@ var init_order_feed = __esm({
       },
       required: ["document_id"]
     };
-    toolDefs48 = [
+    toolDefs49 = [
       {
         name: "get_order_feed",
         pack: "fabricate",
@@ -76790,7 +76459,7 @@ ${issues.map((i) => i.problem).join("\n")}`
     note: "Rollout compiled to the document timeline. Preview with render_sequence, ship with export_video."
   });
 }
-var MAX_FRAMES, MAX_GEO_SAMPLES, MAX_VERIFY_FRAMES, INLINE_GLB_MAX_BASE64, DEFAULT_WIDTH_PX3, MIN_WIDTH_PX2, MAX_WIDTH_PX2, SVG_SCALE3, KEYFRAMES_DESCRIPTION, animateSchema, renderSequenceSchema, CAMERA_NODE, exportVideoSchema, ffmpegAvailable, timelineFromSimulationSchema, toolDefs49;
+var MAX_FRAMES, MAX_GEO_SAMPLES, MAX_VERIFY_FRAMES, INLINE_GLB_MAX_BASE64, DEFAULT_WIDTH_PX3, MIN_WIDTH_PX2, MAX_WIDTH_PX2, SVG_SCALE3, KEYFRAMES_DESCRIPTION, animateSchema, renderSequenceSchema, CAMERA_NODE, exportVideoSchema, ffmpegAvailable, timelineFromSimulationSchema, toolDefs50;
 var init_animate = __esm({
   "src/tools/animate.ts"() {
     "use strict";
@@ -76942,7 +76611,7 @@ var init_animate = __esm({
       },
       required: ["env_id", "document_id"]
     };
-    toolDefs49 = [
+    toolDefs50 = [
       {
         name: "animate",
         pack: null,
@@ -77787,6 +77456,7 @@ var init_server3 = __esm({
     init_antenna();
     init_photonics();
     init_neutronics();
+    init_qcd();
     init_physics2();
     init_loon_macros();
     init_dfm2();
@@ -77805,8 +77475,8 @@ var init_server3 = __esm({
     init_order_feed();
     init_animate();
     PKG_VERSION = (() => {
-      if ("0.9.4-main.15") {
-        return "0.9.4-main.15";
+      if ("0.9.4-main.17") {
+        return "0.9.4-main.17";
       }
       try {
         const req = createRequire2(import.meta.url);
@@ -77815,8 +77485,8 @@ var init_server3 = __esm({
         return "0.0.0";
       }
     })();
-    BUILD_SHA = "ba0bbea017e3e720ec8796167b6824f3cdf235cf";
-    BUILD_TIME = "2026-07-24T17:18:07.764Z";
+    BUILD_SHA = "cfb504c93d13b7e0339f79b4b8863ffeb1ca83b1";
+    BUILD_TIME = "2026-07-24T20:18:23.601Z";
     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);
@@ -77876,8 +77546,8 @@ var init_server3 = __esm({
       ...toolDefs34,
       ...toolDefs35,
       ...toolDefs36,
-      ...toolDefs2,
       ...toolDefs37,
+      ...toolDefs2,
       ...toolDefs38,
       ...toolDefs39,
       ...toolDefs40,
@@ -77887,11 +77557,12 @@ var init_server3 = __esm({
       ...toolDefs44,
       ...toolDefs45,
       ...toolDefs46,
+      ...toolDefs47,
       ...toolDefs20,
       ...toolDefs19,
-      ...toolDefs47,
       ...toolDefs48,
-      ...toolDefs49
+      ...toolDefs49,
+      ...toolDefs50
     ];
     LIST_TOOL_ORDER = [
       // ── Session lifecycle ──────────────────────────────────────
@@ -77962,6 +77633,7 @@ var init_server3 = __esm({
       "solve_thermal",
       "analyze_structure",
       "simulate_neutron_shield",
+      "simulate_lattice_gauge",
       "analyze_tolerance_stackup",
       "simulate_circuit",
       "tune_circuit",