@vcad/mcp 0.9.4-main.27 → 0.9.4-main.29

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 cfda229f9a317556e3c9f64a43700549cd5265ef at 2026-07-25T23:24:45.817Z. Source: https://github.com/ecto/vcad
+Built from d80a65e336d1c8b8ebf75dd1bebc553b41f7ca93 at 2026-07-26T00:45:54.095Z. Source: https://github.com/ecto/vcad

package/index.mjs

@@ -113,6 +113,7 @@ __export(vcad_kernel_wasm_exports, {
   documentToLoonChecked: () => documentToLoonChecked,
   drawingSheetToPdf: () => drawingSheetToPdf,
   ecadAirgapFluxDensity: () => ecadAirgapFluxDensity,
+  ecadAirgapSolve: () => ecadAirgapSolve,
   ecadBuildReceipt: () => ecadBuildReceipt,
   ecadBuiltinSymbols: () => ecadBuiltinSymbols,
   ecadCheckDrc: () => ecadCheckDrc,
@@ -1142,6 +1143,15 @@ function ecadAirgapFluxDensity(spec_json) {
   }
   return ret[0];
 }
+function ecadAirgapSolve(spec_json) {
+  const ptr0 = passStringToWasm0(spec_json, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+  const len0 = WASM_VECTOR_LEN;
+  const ret = wasm.ecadAirgapSolve(ptr0, len0);
+  if (ret[2]) {
+    throw takeFromExternrefTable0(ret[1]);
+  }
+  return takeFromExternrefTable0(ret[0]);
+}
 function ecadBuildReceipt(pcb_json) {
   const ptr0 = passStringToWasm0(pcb_json, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
   const len0 = WASM_VECTOR_LEN;
@@ -10367,6 +10377,20 @@ async function airgapFluxDensity(spec) {
     return null;
   }
 }
+async function airgapSolve(spec) {
+  const wasm3 = await loadEcadWasm();
+  if (!wasm3)
+    return null;
+  try {
+    const fn = wasm3.ecadAirgapSolve;
+    if (typeof fn !== "function")
+      return null;
+    return fn(JSON.stringify(spec));
+  } catch (e) {
+    console.warn("[ECAD] airgapSolve failed:", e);
+    return null;
+  }
+}
 async function generateNetlist(sheet) {
   const wasm3 = await loadEcadWasm();
   if (!wasm3)
@@ -11229,6 +11253,7 @@ __export(dist_exports, {
   PhysicsEnv: () => PhysicsEnv,
   SolidCache: () => SolidCache,
   airgapFluxDensity: () => airgapFluxDensity,
+  airgapSolve: () => airgapSolve,
   applyForwardKinematics: () => applyForwardKinematics,
   buildMoleculeReceipt: () => buildReceipt2,
   buildPartDocument: () => buildPartDocument,
@@ -40898,6 +40923,26 @@ var init_CHANGELOG = __esm({
             "sheet_metal_create"
           ]
         },
+        {
+          id: "2026-07-25-shell-and-booleans-on-filleted-solids",
+          version: "0.9.4",
+          date: "2026-07-25",
+          category: "fix",
+          title: "Shelling and cutting a filleted solid no longer leaks",
+          summary: "Hollowing a filleted part scattered its cavity and a pocket cut collapsed the neighboring blends flat; both are now exact and watertight, and the kernel is reproducible run to run.",
+          features: [
+            "kernel",
+            "fillet",
+            "shell",
+            "booleans",
+            "tessellation"
+          ],
+          mcpTools: [
+            "create_cad_loon",
+            "inspect_cad",
+            "export_cad"
+          ]
+        },
         {
           id: "2026-07-25-router-window-repair-commit-parity",
           version: "0.9.4",
@@ -41169,6 +41214,57 @@ var init_CHANGELOG = __esm({
             "quote_manufacturing"
           ]
         },
+        {
+          id: "2026-07-25-em-saturated-picard",
+          version: "0.9.4",
+          date: "2026-07-25",
+          category: "fix",
+          title: "simulate_em solves strongly saturating magnetostatics",
+          summary: "Saturated motor cross-sections no longer fail with an unfixable 'Picard not converged': the nonlinear loop's cap, tolerance and relaxation are now caller-settable, and the default cap was raised.",
+          features: [
+            "em",
+            "magnetostatics",
+            "motors"
+          ],
+          mcpTools: [
+            "simulate_em"
+          ]
+        },
+        {
+          id: "2026-07-25-curved-trim-seams-weld",
+          version: "0.9.4",
+          date: "2026-07-25",
+          category: "fix",
+          title: "Bores and cut curved faces close their seams",
+          summary: "A through-hole left the shell open along the bore wall \u2014 visible as hairline cracks and unusable for 3D print or STL export; bores, holed plates and cut fillet corners now come out watertight.",
+          features: [
+            "kernel",
+            "booleans",
+            "tessellation",
+            "export"
+          ],
+          mcpTools: [
+            "create_cad_loon",
+            "inspect_cad",
+            "export_cad"
+          ]
+        },
+        {
+          id: "2026-07-25-calc-motor-tooth-saturation",
+          version: "0.9.4",
+          date: "2026-07-25",
+          category: "fix",
+          title: "calc_motor sees saturating stator teeth",
+          summary: "calc_motor's MEC model treated iron as linear and silently over-predicted Kt on machines whose teeth saturate; it now reports tooth flux, warns past the knee, and can solve the saturating network.",
+          features: [
+            "motor",
+            "magnetics",
+            "simulation"
+          ],
+          mcpTools: [
+            "calc_motor"
+          ]
+        },
         {
           id: "2026-07-24-visionos-app",
           version: "0.9.4",
@@ -64125,6 +64221,7 @@ async function calcMotor(args) {
   let bGap = num2(args.airgap_flux_tesla);
   let bGapSource = "supplied";
   let magnetSpec;
+  let solution = null;
   let fringing;
   if (!Number.isFinite(bGap)) {
     const m = args.magnet ?? {};
@@ -64138,15 +64235,55 @@ async function calcMotor(args) {
       gapAreaMm2: mnum(m.gap_area_mm2, 1),
       ironMuRel: typeof m.iron_mu_rel === "number" ? m.iron_mu_rel : null,
       ironPathMm: mnum(m.iron_path_mm, 0),
-      ironAreaMm2: mnum(m.iron_area_mm2, 1)
-    };
-    const computed = await airgapFluxDensity(magnetSpec);
-    if (computed == null) {
-      return fail2(
-        "air-gap flux is required: pass airgap_flux_tesla, or `magnet` params (ECAD WASM must be available to compute B_gap)."
-      );
+      ironAreaMm2: mnum(m.iron_area_mm2, 1),
+      ironJsT: typeof m.iron_js_t === "number" ? m.iron_js_t : null
+    };
+    if (m.slots !== void 0) {
+      const slots = num2(m.slots);
+      if (!(slots > 0)) return fail2("magnet.slots must be > 0");
+      const meanR = Number.isFinite(num2(m.mean_radius_mm)) ? num2(m.mean_radius_mm) : (innerR + outerR) / 2;
+      const pitch = 2 * Math.PI * meanR / slots;
+      let toothW = num2(m.tooth_width_mm);
+      if (!Number.isFinite(toothW)) {
+        const frac = num2(m.tooth_fraction);
+        if (!Number.isFinite(frac)) {
+          return fail2("magnet.slots requires magnet.tooth_width_mm or magnet.tooth_fraction");
+        }
+        if (!(frac > 0 && frac <= 1)) return fail2("magnet.tooth_fraction must be in (0, 1]");
+        toothW = frac * pitch;
+      }
+      if (!(toothW > 0)) return fail2("magnet.tooth_width_mm must be > 0");
+      magnetSpec.teeth = {
+        slots,
+        toothWidthMm: toothW,
+        meanRadiusMm: meanR,
+        toothPathMm: mnum(m.tooth_path_mm, 0)
+      };
+    }
+    solution = await airgapSolve(magnetSpec);
+    if (solution == null) {
+      const computed = await airgapFluxDensity(magnetSpec);
+      if (computed == null) {
+        return fail2(
+          "air-gap flux is required: pass airgap_flux_tesla, or `magnet` params (ECAD WASM must be available to compute B_gap)."
+        );
+      }
+      const t2 = magnetSpec.teeth;
+      const k = t2 ? Math.max(1, 2 * Math.PI * t2.meanRadiusMm / t2.slots / t2.toothWidthMm) : null;
+      solution = {
+        bGapTesla: computed,
+        bToothTesla: k == null ? null : computed * k,
+        bIronTesla: null,
+        toothConcentration: k,
+        nonlinear: false,
+        iterations: 0,
+        converged: true,
+        warnings: k != null && computed * k > SILICON_STEEL_KNEE_T ? [
+          `tooth flux density ${(computed * k).toFixed(2)} T exceeds the ${SILICON_STEEL_KNEE_T.toFixed(2)} T knee \u2014 the LINEAR iron model over-predicts B_gap and Kt here; pass magnet.iron_js_t to solve the saturating network`
+        ] : []
+      };
     }
-    bGap = computed;
+    bGap = solution.bGapTesla;
     bGapSource = "computed";
     if (m.pole_width_mm !== void 0) {
       const poleW = num2(m.pole_width_mm);
@@ -64218,6 +64355,24 @@ async function calcMotor(args) {
         })
       );
     }
+    if (solution?.bToothTesla != null && magnetSpec.teeth) {
+      claims.push(
+        emClaim(
+          "tooth_flux_density",
+          r4(solution.bToothTesla),
+          "T",
+          solution.nonlinear ? "mec-saturating-iron" : "mec-tooth-concentration",
+          {
+            ...mecInputs,
+            slots: magnetSpec.teeth.slots,
+            tooth_width_mm: r4(magnetSpec.teeth.toothWidthMm),
+            mean_radius_mm: r4(magnetSpec.teeth.meanRadiusMm),
+            tooth_concentration: r4(solution.toothConcentration ?? 1),
+            ...magnetSpec.ironJsT != null ? { iron_js_t: magnetSpec.ironJsT } : {}
+          }
+        )
+      );
+    }
   }
   return {
     content: [
@@ -64233,6 +64388,16 @@ async function calcMotor(args) {
             fringing_derate: r4(fringing.derate),
             fringing_note: "Carter-like first-order derate: B under the pole drops by w/(w+2g) as flux fringes ~one gap length past each pole edge. Kt/Ke and the curve use the DERATED flux. Honest for pole width \u2273 2\xD7 gap; below that treat it as a lower bound."
           } : {},
+          ...solution?.bToothTesla != null ? {
+            tooth_flux_tesla: r4(solution.bToothTesla),
+            tooth_concentration: r4(solution.toothConcentration ?? 1),
+            // Fringing redistributes flux under the pole; it does not remove
+            // it, so the tooth carries the RAW gap flux either way.
+            tooth_note: "B_tooth = raw B_gap \xB7 (tooth pitch / tooth width) \u2014 the gap flux of a whole tooth pitch funnels into one tooth body. Compare against your steel's knee (~1.5-1.7 T for M19-class silicon steel)."
+          } : {},
+          ...solution?.bIronTesla != null ? { yoke_flux_tesla: r4(solution.bIronTesla) } : {},
+          ...solution?.nonlinear ? { iron_model: "saturating (arctangent B-H)" } : {},
+          ...solution && solution.warnings.length > 0 ? { warnings: solution.warnings } : {},
           winding_factor: windingFactor,
           kt_nm_per_a: r4(perf.ktNmPerA),
           ke_v_s_per_rad: r4(perf.keVSPerRad),
@@ -64243,7 +64408,7 @@ async function calcMotor(args) {
             speed_rad_s: r4(p.speedRadS),
             torque_nm: r4(p.torqueNm)
           })),
-          note: fringing ? "First-order steady-state estimate (no slotting/saturation/losses; pole-edge fringing derated via magnet.pole_width_mm)." : "First-order steady-state estimate (no slotting/fringing/saturation/losses; pass magnet.pole_width_mm to derate for fringing).",
+          note: "First-order steady-state estimate (no slotting, no losses)" + (fringing ? "; pole-edge fringing derated via magnet.pole_width_mm" : "; pass magnet.pole_width_mm to derate for fringing") + (solution?.nonlinear ? "; iron solved with its saturating B-H law" : "; iron is LINEAR, so Kt is OPTIMISTIC on a machine whose teeth saturate \u2014 " + (solution?.bToothTesla != null ? "see tooth_flux_tesla, and pass magnet.iron_js_t to solve the saturating network" : "pass magnet.slots + tooth_width_mm to see the tooth field")) + ".",
           claims
         })
       }
@@ -65715,7 +65880,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, fabPrepSchema, 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 PCB_LAYER_SET, COPPER_LAYERS, COPPER_LAYER_SET, POSITION_TOLERANCE, PIN_SEP, pinKey, createSchematicSchema, docInputProperties, placeComponentsSchema, routeNetsSchema, runDrcSchema, runErcSchema, critiqueRouteSchema, routeDiffPairSchema, lengthMatchTracesSchema, exportGerberSchema, validateForFabSchema, fabPrepSchema, 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, SILICON_STEEL_KNEE_T, 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";
@@ -67127,7 +67292,7 @@ var init_ecad2 = __esm({
         },
         magnet: {
           type: "object",
-          description: "PM mode: optional magnet/geometry to compute B_gap when airgap_flux_tesla is omitted (NdFeB defaults). Fields: remanence_tesla, magnet_thickness_mm, airgap_mm, recoil_mu_rel, magnet_area_mm2, gap_area_mm2, iron_mu_rel. Add pole_width_mm (magnet pole face width across the fringing direction) to also apply the first-order Carter-like fringing derate w/(w+2g) \u2014 the tool then reports raw AND derated B and uses the derated value for Kt."
+          description: "PM mode: optional magnet/geometry to compute B_gap when airgap_flux_tesla is omitted (NdFeB defaults). Fields: remanence_tesla, magnet_thickness_mm, airgap_mm, recoil_mu_rel, magnet_area_mm2, gap_area_mm2, iron_mu_rel. Add pole_width_mm (magnet pole face width across the fringing direction) to also apply the first-order Carter-like fringing derate w/(w+2g) \u2014 the tool then reports raw AND derated B and uses the derated value for Kt. SATURATION: the MEC iron is LINEAR unless you say otherwise, so it over-predicts Kt on a machine whose teeth saturate. Pass slots + (tooth_width_mm or tooth_fraction) to make the teeth visible \u2014 the tool then reports tooth_flux_density (B_gap \xB7 pitch/width, routinely 2x B_gap) and warns when it passes the ~1.5 T silicon-steel knee. Add iron_js_t (saturation polarization, silicon steel \u2248 2.0 T) with iron_mu_rel to actually SOLVE the saturating network instead of merely warning; tooth_path_mm puts the teeth in the reluctance loop."
         },
         phase_current_a: { type: "number", description: "Induction mode (required): phase current, A RMS (balanced 3-phase drive)." },
         electrical_freq_hz: { type: "number", description: "Induction mode (required): electrical drive frequency, Hz." },
@@ -67146,6 +67311,7 @@ var init_ecad2 = __esm({
       },
       required: ["pole_pairs", "turns_per_phase"]
     };
+    SILICON_STEEL_KNEE_T = 1.5;
     sig6 = (v) => v === 0 || !Number.isFinite(v) ? v : Number(v.toPrecision(6));
     BEARING_FRICTION_MNM = {
       "608-2RS": {
@@ -71400,6 +71566,22 @@ var init_em = __esm({
                   type: "number",
                   description: "SOR sweep cap. Default 200000."
                 },
+                picard_max_iters: {
+                  type: "number",
+                  description: "Saturable iron only: cap on the NONLINEAR outer loop over the B-H law. Default 300. Separate from max_sweeps, which caps the inner SOR solve and does nothing for a non-converging material state. A deeply saturated motor section contracts steadily but needs 150-250 iterations; raise this when the error says the residual is still falling."
+                },
+                picard_tol: {
+                  type: "number",
+                  description: "Saturable iron only: nonlinear convergence tolerance on the largest relative reluctivity update. Default 1e-4. Raise it to accept a looser material state deliberately."
+                },
+                picard_relax: {
+                  type: "number",
+                  description: "Saturable iron only: under-relaxation of the nonlinear loop, in (0, 1]. Default 0.7. Lower it ONLY when the error reports the residual has stopped falling (a limit cycle); when the residual is still falling, lowering this makes convergence strictly slower \u2014 raise picard_max_iters instead."
+                },
+                picard_adaptive: {
+                  type: "boolean",
+                  description: "Saturable iron only: back the relaxation off automatically when the nonlinear residual makes no progress over a 10-iteration window. Default FALSE \u2014 on a healthy saturated solve it reacts to ordinary residual jitter and costs ~1.7x the iterations. Turn it on for a device that is demonstrably oscillating."
+                },
                 drive_coil: {
                   type: "number",
                   description: "Axisym: coil index the inductance claim is priced for. Default 0."
@@ -78869,8 +79051,8 @@ var init_server3 = __esm({
     init_order_feed();
     init_animate();
     PKG_VERSION = (() => {
-      if ("0.9.4-main.27") {
-        return "0.9.4-main.27";
+      if ("0.9.4-main.29") {
+        return "0.9.4-main.29";
       }
       try {
         const req = createRequire2(import.meta.url);
@@ -78879,8 +79061,8 @@ var init_server3 = __esm({
         return "0.0.0";
       }
     })();
-    BUILD_SHA = "cfda229f9a317556e3c9f64a43700549cd5265ef";
-    BUILD_TIME = "2026-07-25T23:24:45.817Z";
+    BUILD_SHA = "d80a65e336d1c8b8ebf75dd1bebc553b41f7ca93";
+    BUILD_TIME = "2026-07-26T00:45:54.095Z";
     SHORT_SHA = BUILD_SHA === "unknown" ? "unknown" : BUILD_SHA.slice(0, 7);
     VERSION_WITH_BUILD = SHORT_SHA === "unknown" ? PKG_VERSION : `${PKG_VERSION}+${SHORT_SHA}`;
     INSTANCE_ID = randomUUID6().slice(0, 8);

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@vcad/mcp",
-  "version": "0.9.4-main.27",
+  "version": "0.9.4-main.29",
   "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": "cfda229f9a317556e3c9f64a43700549cd5265ef",
-    "builtAt": "2026-07-25T23:24:45.817Z"
+    "sha": "d80a65e336d1c8b8ebf75dd1bebc553b41f7ca93",
+    "builtAt": "2026-07-26T00:45:54.095Z"
   }
 }