@vcad/mcp 0.9.4-main.16 → 0.9.4-main.18

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 63c1a040cf3352ca640144ec217f3f96271180c9 at 2026-07-24T17:21:48.535Z. Source: https://github.com/ecto/vcad
+Built from 48276771946e67b08cb31a66d123bdc2d7af926d at 2026-07-24T22:30:13.865Z. Source: https://github.com/ecto/vcad

package/index.mjs

@@ -151,6 +151,14 @@ __export(vcad_kernel_wasm_exports, {
   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,
@@ -190,6 +198,7 @@ __export(vcad_kernel_wasm_exports, {
   isGpuAvailable: () => isGpuAvailable,
   isPhysicsAvailable: () => isPhysicsAvailable,
   isSlicerAvailable: () => isSlicerAvailable,
+  latticeGaugeSimulate: () => latticeGaugeSimulate,
   mesh_clearance: () => mesh_clearance,
   nestSheetMetalParts: () => nestSheetMetalParts,
   nestedSheetMetalDxf: () => nestedSheetMetalDxf,
@@ -1530,6 +1539,184 @@ function embroideryDesignToMesh(design_json) {
   }
   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;
@@ -1940,6 +2127,15 @@ function isSlicerAvailable() {
   const ret = wasm.isSlicerAvailable();
   return ret !== 0;
 }
+function latticeGaugeSimulate(spec_json) {
+  const ptr0 = passStringToWasm0(spec_json, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+  const len0 = WASM_VECTOR_LEN;
+  const ret = wasm.latticeGaugeSimulate(ptr0, len0);
+  if (ret[2]) {
+    throw takeFromExternrefTable0(ret[1]);
+  }
+  return takeFromExternrefTable0(ret[0]);
+}
 function mesh_clearance(positions_a, indices_a, positions_b, indices_b) {
   const ptr0 = passArrayF32ToWasm0(positions_a, wasm.__wbindgen_malloc);
   const len0 = WASM_VECTOR_LEN;
@@ -11269,6 +11465,7 @@ var init_dist = __esm({
           toleranceAnalyze: wasmModule7.toleranceAnalyze,
           thermalSolve: wasmModule7.thermalSolve,
           thermalSolveTransient: wasmModule7.thermalSolveTransient,
+          simulateFlow: wasmModule7.simulateFlow,
           documentDiff: wasmModule7.documentDiff,
           documentDiffApply: wasmModule7.documentDiffApply,
           documentMerge: wasmModule7.documentMerge,
@@ -11283,6 +11480,7 @@ var init_dist = __esm({
           antennaAnalyze: wasmModule7.antennaAnalyze,
           photonicsSimulate: wasmModule7.photonicsSimulate,
           neutronicsSimulate: wasmModule7.neutronicsSimulate,
+          latticeGaugeSimulate: wasmModule7.latticeGaugeSimulate,
           analyzeStaticsBox: wasmModule7.analyzeStaticsBox,
           analyzeStaticsMesh: wasmModule7.analyzeStaticsMesh
         }, compiledWasmModule);
@@ -11546,6 +11744,20 @@ var init_dist = __esm({
         }
         return fn(specJson, transientJson, paramsJson, optionsJson);
       }
+      /**
+       * Steady laminar flow solve (D3Q19 BGK lattice Boltzmann): pressure
+       * drop, flow rates, mass audit, optional thermal pickup, and predicted
+       * receipt claims. Per-voxel fields are only returned when
+       * `includeFields` is true. Inputs are JSON strings (FlowSpec, options);
+       * see `vcad-kernel-flow`.
+       */
+      simulateFlow(specJson, optionsJson, includeFields) {
+        const fn = this.kernel.simulateFlow;
+        if (typeof fn !== "function") {
+          throw new Error("simulateFlow is not exported by this kernel WASM build \u2014 rebuild packages/kernel-wasm");
+        }
+        return fn(specJson, optionsJson, includeFields);
+      }
       circuitFn(name) {
         const fn = this.kernel[name];
         if (typeof fn !== "function") {
@@ -11643,6 +11855,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") {
@@ -40609,6 +40834,37 @@ 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-24-flow-cfd-lbm",
+          version: "0.9.4",
+          date: "2026-07-24",
+          category: "feat",
+          title: "Laminar CFD: lattice-Boltzmann flow solver",
+          summary: "New simulate_flow tool: steady laminar D3Q19 LBM on a voxel grid with pressure drop, flow rates, mass audit, thermal pickup, and predicted vcad.flow-claims/1 receipts.",
+          features: [
+            "flow",
+            "cfd",
+            "simulation"
+          ],
+          mcpTools: [
+            "simulate_flow"
+          ]
+        },
         {
           id: "2026-07-23-router-post-route-legalization",
           version: "0.9.4",
@@ -51632,11 +51888,16 @@ var init_tool_metadata = __esm({
       simulate_photonics: { title: "Simulate Photonics", annotations: RO },
       analyze_antenna: { title: "Analyze Antenna", annotations: RO },
       solve_thermal: { title: "Solve Thermal", annotations: RO },
+      simulate_flow: { title: "Simulate Flow", annotations: RO },
       analyze_structure: { title: "Analyze Structure", annotations: RO },
       simulate_neutron_shield: {
         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
@@ -69352,7 +69613,7 @@ var init_thermal = __esm({
   }
 });
 
-// src/tools/circuit.ts
+// src/tools/flow.ts
 function textResult5(payload) {
   return {
     content: [
@@ -69363,6 +69624,207 @@ function textResult5(payload) {
     ]
   };
 }
+var vec32, shapeSchema2, flowSpecSchema, toolDefs30;
+var init_flow = __esm({
+  "src/tools/flow.ts"() {
+    "use strict";
+    init_tool_def();
+    vec32 = {
+      type: "array",
+      minItems: 3,
+      maxItems: 3,
+      items: { type: "number" }
+    };
+    shapeSchema2 = {
+      type: "object",
+      required: ["type"],
+      description: 'Region shape: `{"type":"Box","min_mm":[..],"size_mm":[..]}` or `{"type":"Tube","axis":"Z","center_mm":[..2],"span_mm":[..2],"outer_radius_mm":..,"inner_radius_mm":..}`.',
+      properties: {
+        type: { type: "string", enum: ["Box", "Tube"] },
+        min_mm: vec32,
+        size_mm: vec32,
+        axis: { type: "string", enum: ["X", "Y", "Z"] },
+        center_mm: {
+          type: "array",
+          minItems: 2,
+          maxItems: 2,
+          items: { type: "number" }
+        },
+        span_mm: {
+          type: "array",
+          minItems: 2,
+          maxItems: 2,
+          items: { type: "number" }
+        },
+        outer_radius_mm: { type: "number" },
+        inner_radius_mm: { type: "number" }
+      }
+    };
+    flowSpecSchema = {
+      type: "object",
+      required: ["origin_mm", "size_mm", "divisions"],
+      description: "Voxelized flow domain: an axis-aligned box discretized into cubic voxels (`size_mm[a]/divisions[a]` must agree across axes), painted with regions in order (`background` first, then `regions`). Units mm / kg / m / s / Pa / degC.",
+      properties: {
+        origin_mm: vec32,
+        size_mm: vec32,
+        divisions: {
+          type: "array",
+          minItems: 3,
+          maxItems: 3,
+          items: { type: "number" },
+          description: "Voxel counts per axis (voxels must come out cubic; the cost knob)."
+        },
+        background: {
+          type: "string",
+          enum: ["solid", "fluid"],
+          description: "What unpainted voxels are. Default solid."
+        },
+        regions: {
+          type: "array",
+          items: {
+            type: "object",
+            required: ["shape", "role"],
+            properties: {
+              shape: shapeSchema2,
+              role: {
+                type: "string",
+                enum: ["solid", "fluid", "inlet", "outlet"]
+              }
+            }
+          },
+          description: "Painted in order; later regions override earlier ones. Inlet and outlet patches must sit on the domain boundary."
+        },
+        fluid: {
+          type: "object",
+          properties: {
+            density_kg_m3: { type: "number" },
+            viscosity_pa_s: { type: "number" }
+          },
+          description: "Fluid properties. Default water (998, 1.002e-3)."
+        },
+        inlet_velocity_m_s: {
+          ...vec32,
+          description: "Plug-flow inlet velocity vector, m/s."
+        },
+        outlet_gauge_pa: {
+          type: "number",
+          description: "Outlet gauge pressure, Pa. Default 0."
+        },
+        body_force_n_m3: {
+          ...vec32,
+          description: "Volumetric body force, N/m^3 (periodic channel drive; needs `options.u_ref_m_s`)."
+        },
+        periodic: {
+          type: "array",
+          minItems: 3,
+          maxItems: 3,
+          items: { type: "boolean" },
+          description: "Per-axis periodic boundaries."
+        },
+        re_envelope: {
+          type: "number",
+          description: "Reynolds gate: the solve refuses Re above this (laminar only). Default 2300."
+        },
+        thermal: {
+          type: "object",
+          properties: {
+            inlet_temp_c: { type: "number" },
+            initial_temp_c: { type: "number" },
+            wall_temp_c: { type: "number" },
+            diffusivity_m2_s: { type: "number" },
+            heat_capacity_j_kg_k: { type: "number" },
+            buoyancy: {
+              type: "object",
+              properties: {
+                beta_per_k: { type: "number" },
+                t_ref_c: { type: "number" },
+                gravity_m_s2: { type: "number" }
+              }
+            }
+          },
+          description: "Enable the advected temperature field (film-averaged conjugate seam; optional Boussinesq buoyancy, gated at Ra <= 1e8)."
+        },
+        hot_walls: {
+          type: "array",
+          items: {
+            type: "object",
+            required: ["shape", "temp_c"],
+            properties: {
+              shape: shapeSchema2,
+              temp_c: { type: "number" }
+            }
+          },
+          description: "Isothermal wall patches (Dirichlet) for thermal runs, painted over solid voxels adjacent to fluid."
+        }
+      }
+    };
+    toolDefs30 = [
+      {
+        name: "simulate_flow",
+        pack: null,
+        description: 'Steady laminar CFD on a voxel grid: D3Q19 BGK lattice Boltzmann with half-way bounce-back walls, velocity inlets, pressure outlets, and an optional advected temperature field (isothermal walls, Boussinesq buoyancy). Returns pressure drop, inlet/outlet flow rates, the mass-balance audit, max speed, and thermal pickup \u2014 as data plus `vcad.flow-claims/1` and unified-receipt claims. Predictions carry basis "predicted" and roll up Provisional until hardware is measured. Laminar only: fail-closed gates at Re <= 2300, the stable tau window, and Ra <= 1e8; weakly compressible (pressure noise O(Ma^2)); walls are voxel staircases. Per-voxel velocity/pressure/temperature fields are only returned with `include_fields: true` \u2014 they are grid-sized (voxel cap 2M). Steadiness is detected, not assumed: a run that never goes steady is an error, not a result.',
+        inputSchema: {
+          type: "object",
+          required: ["spec"],
+          properties: {
+            spec: flowSpecSchema,
+            options: {
+              type: "object",
+              properties: {
+                u_ref_m_s: {
+                  type: "number",
+                  description: "Reference speed for unit scaling, m/s. Defaults to the inlet speed; body-force-driven cases must supply it."
+                },
+                max_steps: {
+                  type: "number",
+                  description: "Step budget before failing. Default 400000."
+                },
+                check_every: {
+                  type: "number",
+                  description: "Steadiness check interval, steps. Default 200."
+                },
+                steady_tol: {
+                  type: "number",
+                  description: "Relative L-inf velocity change per check. Default 1e-6."
+                },
+                ramp_steps: {
+                  type: "number",
+                  description: "Inlet velocity ramp length, steps. Default 1000."
+                }
+              }
+            },
+            include_fields: {
+              type: "boolean",
+              description: "Return the per-voxel velocity/pressure/temperature fields (grid-sized). Default false."
+            }
+          }
+        },
+        handler: (args, ctx) => {
+          const a = args;
+          const out = ctx.engine.simulateFlow(
+            JSON.stringify(a.spec),
+            JSON.stringify(a.options ?? {}),
+            !!a.include_fields
+          );
+          return textResult5(out);
+        },
+        behavior: behavior({})
+      }
+    ];
+  }
+});
+
+// src/tools/circuit.ts
+function textResult6(payload) {
+  return {
+    content: [
+      {
+        type: "text",
+        text: JSON.stringify(payload, null, 2)
+      }
+    ]
+  };
+}
 function frequencyGrid(ac) {
   const list = ac.frequenciesHz;
   if (Array.isArray(list) && list.length > 0) return list;
@@ -69378,7 +69840,7 @@ function frequencyGrid(ac) {
   }
   return out;
 }
-var deviceSchema, devicesSchema, toolDefs30;
+var deviceSchema, devicesSchema, toolDefs31;
 var init_circuit = __esm({
   "src/tools/circuit.ts"() {
     "use strict";
@@ -69415,7 +69877,7 @@ var init_circuit = __esm({
       items: deviceSchema,
       description: "Netlist as data. Node ids are dense integers with 0 = ground; the device's index in this array is its device id (used by sourceId / freeDevices / outNode references)."
     };
-    toolDefs30 = [
+    toolDefs31 = [
       {
         name: "simulate_circuit",
         pack: null,
@@ -69524,7 +69986,7 @@ var init_circuit = __esm({
               tr.sampleEvery ?? 1
             );
           }
-          return textResult5(out);
+          return textResult6(out);
         },
         behavior: behavior({})
       },
@@ -69587,7 +70049,7 @@ var init_circuit = __esm({
             ...typeof a.maxIters === "number" ? { maxIters: a.maxIters } : {}
           };
           const specJson = JSON.stringify({ devices: a.devices });
-          return textResult5(ctx.engine.circuitTune(specJson, JSON.stringify(tune)));
+          return textResult6(ctx.engine.circuitTune(specJson, JSON.stringify(tune)));
         },
         behavior: behavior({})
       }
@@ -69596,7 +70058,7 @@ var init_circuit = __esm({
 });
 
 // src/tools/structure.ts
-function textResult6(payload) {
+function textResult7(payload) {
   return {
     content: [
       {
@@ -69606,7 +70068,7 @@ function textResult6(payload) {
     ]
   };
 }
-var regionSchema2, toolDefs31;
+var regionSchema2, toolDefs32;
 var init_structure = __esm({
   "src/tools/structure.ts"() {
     "use strict";
@@ -69633,7 +70095,7 @@ var init_structure = __esm({
         }
       }
     };
-    toolDefs31 = [
+    toolDefs32 = [
       {
         name: "analyze_structure",
         pack: null,
@@ -69745,7 +70207,7 @@ var init_structure = __esm({
             resolved.mesh.indices
           );
           const solveMs = Math.round(performance.now() - started);
-          return textResult6({
+          return textResult7({
             document_id: a.document_id,
             part: resolved.name ?? a.part,
             solve_ms: solveMs,
@@ -69759,7 +70221,7 @@ var init_structure = __esm({
 });
 
 // src/tools/em.ts
-function textResult7(payload) {
+function textResult8(payload) {
   return {
     content: [
       {
@@ -69769,7 +70231,7 @@ function textResult7(payload) {
     ]
   };
 }
-var paramValue4, rectRegion, bcSchema, electroShape, emSpecSchema, toolDefs32;
+var paramValue4, rectRegion, bcSchema, electroShape, emSpecSchema, toolDefs33;
 var init_em = __esm({
   "src/tools/em.ts"() {
     "use strict";
@@ -69925,7 +70387,7 @@ var init_em = __esm({
         }
       }
     };
-    toolDefs32 = [
+    toolDefs33 = [
       {
         name: "simulate_em",
         pack: null,
@@ -70003,7 +70465,7 @@ var init_em = __esm({
             JSON.stringify(a.parameters ?? {}),
             JSON.stringify(a.options ?? {})
           );
-          return textResult7(out);
+          return textResult8(out);
         },
         behavior: behavior({})
       }
@@ -70012,7 +70474,7 @@ var init_em = __esm({
 });
 
 // src/tools/antenna.ts
-function textResult8(payload) {
+function textResult9(payload) {
   return {
     content: [
       {
@@ -70022,7 +70484,7 @@ function textResult8(payload) {
     ]
   };
 }
-var paramValue5, pointSpec, antennaSpecSchema, toolDefs33;
+var paramValue5, pointSpec, antennaSpecSchema, toolDefs34;
 var init_antenna = __esm({
   "src/tools/antenna.ts"() {
     "use strict";
@@ -70071,7 +70533,7 @@ var init_antenna = __esm({
         feed_mm: pointSpec
       }
     };
-    toolDefs33 = [
+    toolDefs34 = [
       {
         name: "analyze_antenna",
         pack: null,
@@ -70129,7 +70591,7 @@ var init_antenna = __esm({
             JSON.stringify(a.parameters ?? {}),
             JSON.stringify({ band: a.band, ...a.options })
           );
-          return textResult8(out);
+          return textResult9(out);
         },
         behavior: behavior({})
       }
@@ -70138,7 +70600,7 @@ var init_antenna = __esm({
 });
 
 // src/tools/photonics.ts
-function textResult9(payload) {
+function textResult10(payload) {
   return {
     content: [
       {
@@ -70148,7 +70610,7 @@ function textResult9(payload) {
     ]
   };
 }
-var photonicsSpecSchema, toolDefs34;
+var photonicsSpecSchema, toolDefs35;
 var init_photonics = __esm({
   "src/tools/photonics.ts"() {
     "use strict";
@@ -70234,7 +70696,7 @@ var init_photonics = __esm({
         }
       }
     };
-    toolDefs34 = [
+    toolDefs35 = [
       {
         name: "simulate_photonics",
         pack: null,
@@ -70281,7 +70743,7 @@ var init_photonics = __esm({
             JSON.stringify(a.spec),
             JSON.stringify(a.options ?? {})
           );
-          return textResult9(out);
+          return textResult10(out);
         },
         behavior: behavior({})
       }
@@ -70290,7 +70752,7 @@ var init_photonics = __esm({
 });
 
 // src/tools/neutronics.ts
-function textResult10(payload) {
+function textResult11(payload) {
   return {
     content: [
       {
@@ -70300,7 +70762,7 @@ function textResult10(payload) {
     ]
   };
 }
-var paramValue6, shieldSpecSchema, toolDefs35;
+var paramValue6, shieldSpecSchema, toolDefs36;
 var init_neutronics = __esm({
   "src/tools/neutronics.ts"() {
     "use strict";
@@ -70371,7 +70833,7 @@ var init_neutronics = __esm({
         }
       }
     };
-    toolDefs35 = [
+    toolDefs36 = [
       {
         name: "simulate_neutron_shield",
         pack: null,
@@ -70394,7 +70856,129 @@ var init_neutronics = __esm({
             JSON.stringify(a.spec),
             JSON.stringify(a.parameters ?? {})
           );
-          return textResult10(out);
+          return textResult11(out);
+        },
+        behavior: behavior({})
+      }
+    ];
+  }
+});
+
+// src/tools/qcd.ts
+function textResult12(payload) {
+  return {
+    content: [
+      {
+        type: "text",
+        text: JSON.stringify(payload, null, 2)
+      }
+    ]
+  };
+}
+var simSpecSchema, toolDefs37;
+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)."
+        }
+      }
+    };
+    toolDefs37 = [
+      {
+        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 textResult12(out);
         },
         behavior: behavior({})
       }
@@ -70523,7 +71107,7 @@ function predictPhysicsTool(args, engine) {
     ]
   };
 }
-var PHYSICS_DOMAIN, ORACLE, TIER_RESOLUTION, regionSchema3, predictPhysicsSchema, round52, toolDefs36;
+var PHYSICS_DOMAIN, ORACLE, TIER_RESOLUTION, regionSchema3, predictPhysicsSchema, round52, toolDefs38;
 var init_physics2 = __esm({
   "src/tools/physics.ts"() {
     "use strict";
@@ -70636,7 +71220,7 @@ var init_physics2 = __esm({
       }
     };
     round52 = (v) => Number(v.toPrecision(5));
-    toolDefs36 = [
+    toolDefs38 = [
       {
         name: "predict_physics",
         pack: null,
@@ -70878,7 +71462,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, toolDefs39;
 var init_dfm2 = __esm({
   "src/tools/dfm.ts"() {
     "use strict";
@@ -70975,7 +71559,7 @@ var init_dfm2 = __esm({
       },
       required: ["issue_id"]
     };
-    toolDefs37 = [
+    toolDefs39 = [
       {
         name: "dfm_check",
         pack: "dfm",
@@ -71174,7 +71758,7 @@ function renderedOf(engine, doc) {
     "document has no sheet-metal part (evaluation produced no sheetMetal bundle)"
   );
 }
-function textResult11(payload) {
+function textResult13(payload) {
   return {
     content: [{ type: "text", text: JSON.stringify(payload, null, 2) }]
   };
@@ -71216,7 +71800,7 @@ function sheetMetalCreate(input, engine) {
   const errors = rendered.violations.filter(
     (v) => v.severity === "Error"
   ).length;
-  return textResult11({
+  return textResult13({
     document_id: documentId,
     model: rendered.model,
     flat: {
@@ -71232,7 +71816,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 textResult13({
     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)."
@@ -71240,7 +71824,7 @@ function sheetMetalUnfold(input, engine) {
 }
 function sheetMetalMaterials(_input, engine) {
   const materials = engine.getSheetMetalMaterials();
-  return textResult11({
+  return textResult13({
     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.'
   });
@@ -71249,13 +71833,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 textResult13({
       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 textResult13({
     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."
   });
@@ -71277,7 +71861,7 @@ function sheetMetalCheck(input, engine) {
   const errors = result.violations.filter(
     (v) => v.severity === "Error"
   ).length;
-  return textResult11({
+  return textResult13({
     shop: result.shop,
     violations: result.violations,
     shop_ready: result.violations.length === 0,
@@ -71329,7 +71913,7 @@ function sheetMetalNest(input, engine) {
     footprints,
     params
   );
-  return textResult11({
+  return textResult13({
     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)` : ""}.`
@@ -71342,7 +71926,7 @@ function sheetMetalSequence3(input, engine) {
   if (!steps) {
     throw new Error("document has no sheet-metal part");
   }
-  return textResult11({
+  return textResult13({
     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."
@@ -71365,7 +71949,7 @@ function sheetMetalSuggestFix(input, engine) {
     violation: v,
     fix: suggestFix(v)
   }));
-  return textResult11({
+  return textResult13({
     shop: result.shop,
     count: suggestions.length,
     suggestions,
@@ -71468,13 +72052,13 @@ function sheetMetalCost(input, engine) {
   if (!result) {
     throw new Error("document has no sheet-metal part");
   }
-  return textResult11({
+  return textResult13({
     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, toolDefs40;
 var init_sheet_metal2 = __esm({
   "src/tools/sheet-metal.ts"() {
     "use strict";
@@ -71735,7 +72319,7 @@ var init_sheet_metal2 = __esm({
       },
       required: ["document_id"]
     };
-    toolDefs38 = [
+    toolDefs40 = [
       {
         name: "sheet_metal_create",
         pack: "sheet_metal",
@@ -71980,7 +72564,7 @@ function simulateStrike(input, engine) {
   }
   return { content: [{ type: "text", text: JSON.stringify(payload, null, 2) }] };
 }
-var cents, MATERIAL_ALIASES, simulateStrikeSchema, round22, round4, toolDefs39;
+var cents, MATERIAL_ALIASES, simulateStrikeSchema, round22, round4, toolDefs41;
 var init_acoustics = __esm({
   "src/tools/acoustics.ts"() {
     "use strict";
@@ -72054,7 +72638,7 @@ var init_acoustics = __esm({
     };
     round22 = (x) => Math.round(x * 100) / 100;
     round4 = (x) => Math.round(x * 1e4) / 1e4;
-    toolDefs39 = [
+    toolDefs41 = [
       {
         name: "simulate_strike",
         pack: "sheet_metal",
@@ -72205,7 +72789,7 @@ function importStep(input, engine) {
     ]
   };
 }
-var MAX_STEP_BYTES, importStepSchema, toolDefs40;
+var MAX_STEP_BYTES, importStepSchema, toolDefs42;
 var init_import = __esm({
   "src/tools/import.ts"() {
     "use strict";
@@ -72237,7 +72821,7 @@ var init_import = __esm({
         }
       }
     };
-    toolDefs40 = [
+    toolDefs42 = [
       {
         name: "import_step",
         pack: null,
@@ -72365,7 +72949,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, toolDefs43;
 var init_import_pcb = __esm({
   "src/tools/import-pcb.ts"() {
     "use strict";
@@ -72406,7 +72990,7 @@ var init_import_pcb = __esm({
         }
       }
     };
-    toolDefs41 = [
+    toolDefs43 = [
       {
         name: "import_kicad",
         pack: "ecad",
@@ -72502,7 +73086,7 @@ Compressed URL param: ${encoded.length} bytes`
     ]
   };
 }
-var openInBrowserSchema, toolDefs42;
+var openInBrowserSchema, toolDefs44;
 var init_share = __esm({
   "src/tools/share.ts"() {
     "use strict";
@@ -72526,7 +73110,7 @@ var init_share = __esm({
         }
       }
     };
-    toolDefs42 = [
+    toolDefs44 = [
       {
         name: "open_in_browser",
         pack: null,
@@ -72851,7 +73435,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, toolDefs45;
 var init_gym = __esm({
   "src/tools/gym.ts"() {
     "use strict";
@@ -73003,7 +73587,7 @@ var init_gym = __esm({
       },
       required: ["batch_id"]
     };
-    toolDefs43 = [
+    toolDefs45 = [
       {
         name: "create_robot_env",
         pack: "physics",
@@ -73361,7 +73945,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, toolDefs46;
 var init_atoms2 = __esm({
   "src/tools/atoms.ts"() {
     "use strict";
@@ -73496,7 +74080,7 @@ var init_atoms2 = __esm({
       S: 1.05,
       P: 1.07
     };
-    toolDefs44 = [
+    toolDefs46 = [
       {
         name: "load_structure",
         pack: "atoms",
@@ -74704,7 +75288,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, toolDefs47;
 var init_record2 = __esm({
   "src/tools/record.ts"() {
     "use strict";
@@ -74768,7 +75352,7 @@ var init_record2 = __esm({
     DEFAULT_FPS = 30;
     MIN_FPS = 1;
     MAX_FPS = 60;
-    toolDefs45 = [
+    toolDefs47 = [
       {
         name: "record_simulation",
         pack: "physics",
@@ -74861,7 +75445,7 @@ ${text}`
     ]
   };
 }
-var getChangelogSchema, toolDefs46;
+var getChangelogSchema, toolDefs48;
 var init_changelog2 = __esm({
   "src/tools/changelog.ts"() {
     "use strict";
@@ -74901,7 +75485,7 @@ var init_changelog2 = __esm({
         }
       }
     };
-    toolDefs46 = [
+    toolDefs48 = [
       {
         name: "get_changelog",
         pack: null,
@@ -75050,7 +75634,7 @@ function getSimVersion(input) {
     ]
   };
 }
-var getSimReplaySchema, getSimVersionSchema, toolDefs47;
+var getSimReplaySchema, getSimVersionSchema, toolDefs49;
 var init_sim_replay = __esm({
   "src/tools/sim-replay.ts"() {
     "use strict";
@@ -75078,7 +75662,7 @@ var init_sim_replay = __esm({
       },
       required: ["env_id"]
     };
-    toolDefs47 = [
+    toolDefs49 = [
       {
         name: "get_sim_replay",
         pack: null,
@@ -75197,7 +75781,7 @@ async function getOrderFeed(input, store, user) {
     version: version2
   });
 }
-var toUsd4, getOrderFeedSchema, toolDefs48;
+var toUsd4, getOrderFeedSchema, toolDefs50;
 var init_order_feed = __esm({
   "src/tools/order-feed.ts"() {
     "use strict";
@@ -75216,7 +75800,7 @@ var init_order_feed = __esm({
       },
       required: ["document_id"]
     };
-    toolDefs48 = [
+    toolDefs50 = [
       {
         name: "get_order_feed",
         pack: "fabricate",
@@ -75421,7 +76005,7 @@ ${issues.map((i) => i.track >= 0 ? `track[${i.track}]: ${i.problem}` : i.problem
     note: "Timeline set. Preview with render_sequence, ship with export_video."
   });
 }
-function vec32(v) {
+function vec33(v) {
   return v ? [v.x, v.y, v.z] : [0, 0, 0];
 }
 function hasParamTracks(timeline) {
@@ -75464,7 +76048,7 @@ function cameraChannels(frames) {
 }
 function explodeDirections(doc, transforms) {
   const ids = [...transforms.keys()];
-  const pts = ids.map((id) => vec32(transforms.get(id).translation));
+  const pts = ids.map((id) => vec33(transforms.get(id).translation));
   const c = pts.reduce((acc, p) => [acc[0] + p[0], acc[1] + p[1], acc[2] + p[2]], [0, 0, 0]).map((v) => v / Math.max(1, pts.length));
   let extent = 0;
   for (const p of pts) {
@@ -75522,7 +76106,7 @@ function buildSequenceGlb(doc, timeline, frames, engine) {
         for (const inst of instances) {
           const rec = perInstance.get(inst.instanceId);
           const t2 = world.get(inst.instanceId) ?? inst.transform;
-          const translation = vec32(t2?.translation);
+          const translation = vec33(t2?.translation);
           if (frame.explode > 0) {
             const d = dirs.get(inst.instanceId) ?? [0, 0, 0];
             translation[0] += d[0] * frame.explode;
@@ -76118,7 +76702,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, toolDefs51;
 var init_animate = __esm({
   "src/tools/animate.ts"() {
     "use strict";
@@ -76270,7 +76854,7 @@ var init_animate = __esm({
       },
       required: ["env_id", "document_id"]
     };
-    toolDefs49 = [
+    toolDefs51 = [
       {
         name: "animate",
         pack: null,
@@ -77109,12 +77693,14 @@ var init_server3 = __esm({
     init_particle();
     init_tolerance();
     init_thermal();
+    init_flow();
     init_circuit();
     init_structure();
     init_em();
     init_antenna();
     init_photonics();
     init_neutronics();
+    init_qcd();
     init_physics2();
     init_loon_macros();
     init_dfm2();
@@ -77133,8 +77719,8 @@ var init_server3 = __esm({
     init_order_feed();
     init_animate();
     PKG_VERSION = (() => {
-      if ("0.9.4-main.16") {
-        return "0.9.4-main.16";
+      if ("0.9.4-main.18") {
+        return "0.9.4-main.18";
       }
       try {
         const req = createRequire2(import.meta.url);
@@ -77143,8 +77729,8 @@ var init_server3 = __esm({
         return "0.0.0";
       }
     })();
-    BUILD_SHA = "63c1a040cf3352ca640144ec217f3f96271180c9";
-    BUILD_TIME = "2026-07-24T17:21:48.535Z";
+    BUILD_SHA = "48276771946e67b08cb31a66d123bdc2d7af926d";
+    BUILD_TIME = "2026-07-24T22:30:13.865Z";
     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);
@@ -77204,9 +77790,9 @@ var init_server3 = __esm({
       ...toolDefs34,
       ...toolDefs35,
       ...toolDefs36,
-      ...toolDefs2,
       ...toolDefs37,
       ...toolDefs38,
+      ...toolDefs2,
       ...toolDefs39,
       ...toolDefs40,
       ...toolDefs41,
@@ -77215,11 +77801,13 @@ var init_server3 = __esm({
       ...toolDefs44,
       ...toolDefs45,
       ...toolDefs46,
-      ...toolDefs20,
-      ...toolDefs19,
       ...toolDefs47,
       ...toolDefs48,
-      ...toolDefs49
+      ...toolDefs20,
+      ...toolDefs19,
+      ...toolDefs49,
+      ...toolDefs50,
+      ...toolDefs51
     ];
     LIST_TOOL_ORDER = [
       // ── Session lifecycle ──────────────────────────────────────
@@ -77288,8 +77876,10 @@ var init_server3 = __esm({
       "simulate_photonics",
       "analyze_antenna",
       "solve_thermal",
+      "simulate_flow",
       "analyze_structure",
       "simulate_neutron_shield",
+      "simulate_lattice_gauge",
       "analyze_tolerance_stackup",
       "simulate_circuit",
       "tune_circuit",

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@vcad/mcp",
-  "version": "0.9.4-main.16",
+  "version": "0.9.4-main.18",
   "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": "63c1a040cf3352ca640144ec217f3f96271180c9",
-    "builtAt": "2026-07-24T17:21:48.535Z"
+    "sha": "48276771946e67b08cb31a66d123bdc2d7af926d",
+    "builtAt": "2026-07-24T22:30:13.865Z"
   }
 }