brepjs-sheetmetal 0.3.1 → 0.4.1

package/README.md

@@ -300,7 +300,10 @@ authored path. It reads only the B-rep:
   other surface type (cone, sphere, spline, …) is reported `UNSUPPORTED_FACE` and skipped (the unfold never
   silently invents a flat for an unrecognised face).
 - **Pair flats**: the two large parallel planar faces of a panel, a sheet thickness apart, are the
-  mid-surface flat. Thickness is detected as the smallest opposed-planar-pair gap.
+  mid-surface flat. Two antiparallel faces count as opposed only when they face each other (an interior
+  point of one projects inside the other), so a flange's end cap and the adjacent flange's inner face, a
+  bend radius apart across a corner, never pair. Thickness is detected as the smallest opposed-planar-pair
+  gap.
 - **Fit bends**: `fitCylinder(face)` recovers a bend's axis + radius + swept angle from sampled points and
   normals. The axis direction is the common perpendicular of the sampled normals (every cylinder normal is
   ⟂ the axis), recovered as the sign-aligned average of cross products of non-parallel normal pairs; the
@@ -320,8 +323,9 @@ non-uniform thickness, conical/spline bend faces, non-manifold / non-sheet topol
 fit a cylinder within tolerance all surface as warnings inside the `Ok` payload (or a clear `Err` when no
 sheet-metal structure is detected at all, e.g. a bare sphere). The detection is validated by an oracle:
 author a part, take only its `.solid`, run `unfoldSolid`, and confirm the detected unfold reproduces the
-authored unfold's developed area, bend count, flat count, and total flat bbox (single bend, L-bracket,
-U-channel), reading only the solid, never the feature tree.
+authored unfold's developed area, bend count, flat count, and total flat bbox with no warnings (single bend,
+plus the L-bracket, U-channel, and four-sided tray at thicknesses 1, 1.5, 2, and 3 against a 2 mm inner
+radius), reading only the solid, never the feature tree.
 
 ## Bend tables
 

package/dist/brepjs-sheetmetal.cjs

@@ -1,6 +1,6 @@
 Object.defineProperty(exports, Symbol.toStringTag, { value: "Module" });
 let brepjs = require("brepjs");
-//#region src/materials.ts
+// src/materials.ts
 var DEFAULT_K_FACTOR = .44;
 function gauge(name, thickness, kFactor = DEFAULT_K_FACTOR) {
 	return {
@@ -38,8 +38,7 @@ var MATERIALS = {
 function getMaterial(name) {
 	return MATERIALS[name];
 }
-//#endregion
-//#region src/bendTableFns.ts
+// src/bendTableFns.ts
 var DEG_TO_RAD = Math.PI / 180;
 var registry = /* @__PURE__ */ new Map();
 /**
@@ -378,8 +377,7 @@ function ensureStarterTables() {
 		});
 	}
 }
-//#endregion
-//#region src/allowanceFns.ts
+// src/allowanceFns.ts
 function validateRule(rule) {
 	if (!Number.isFinite(rule.innerRadius) || rule.innerRadius < 0) return (0, brepjs.err)((0, brepjs.validationError)("INVALID_RADIUS", `innerRadius must be a finite, non-negative number, got ${rule.innerRadius}`));
 	if (!Number.isFinite(rule.kFactor) || rule.kFactor < 0 || rule.kFactor > 1) return (0, brepjs.err)((0, brepjs.validationError)("INVALID_K_FACTOR", `kFactor must be in [0, 1], got ${rule.kFactor}`));
@@ -413,8 +411,7 @@ function neutralRadius(thickness, rule) {
 	if (!ruleCheck.ok) return ruleCheck;
 	return (0, brepjs.ok)(rule.innerRadius + rule.kFactor * thickness);
 }
-//#endregion
-//#region src/featureTreeFns.ts
+// src/featureTreeFns.ts
 var ROOT_FLAT_ID = "root";
 /**
 * Build the feature graph for an authored part. Each flange contributes a flat
@@ -559,8 +556,7 @@ function featureTree(part, rootId = ROOT_FLAT_ID) {
 	if (!graph.ok) return graph;
 	return buildFeatureTree(graph.value, rootId);
 }
-//#endregion
-//#region src/internal.ts
+// src/internal.ts
 /**
 * OCCT boolean fuse returns a compound wrapping the unioned solid. Extract the
 * single solid so a part's cached shape is a true `Solid` (clean STEP/GLB
@@ -583,8 +579,7 @@ function sideFromAxisDir(axisDir) {
 	if (Math.abs(axisDir[1]) > .5) return "xmax";
 	if (Math.abs(axisDir[0]) > .5) return "ymax";
 }
-//#endregion
-//#region src/unfoldFns.ts
+// src/unfoldFns.ts
 /**
 * Tree-driven analytic unfold of an arbitrary straight-bend part into its flat
 * pattern. Walks the feature tree in BFS order placing every flat in the flat
@@ -1291,8 +1286,62 @@ function reportFromLayout(layout) {
 		totalFlatSize: [maxX - minX, maxY - minY]
 	};
 }
-//#endregion
-//#region src/foreignUnfoldFns.ts
+// \0@oxc-project+runtime@0.155.0/helpers/esm/usingCtx.js
+function _usingCtx() {
+	var r = "function" == typeof SuppressedError ? SuppressedError : function(r, e) {
+		var n = Error();
+		return n.name = "SuppressedError", n.error = r, n.suppressed = e, n;
+	}, e = {}, n = [];
+	function using(r, e) {
+		if (null != e) {
+			if (Object(e) !== e) throw new TypeError("using declarations can only be used with objects, functions, null, or undefined.");
+			if (r) var o = e[Symbol.asyncDispose || Symbol["for"]("Symbol.asyncDispose")];
+			if (void 0 === o && (o = e[Symbol.dispose || Symbol["for"]("Symbol.dispose")], r)) var t = o;
+			if ("function" != typeof o) throw new TypeError("Object is not disposable.");
+			t && (o = function o() {
+				try {
+					t.call(e);
+				} catch (r) {
+					return Promise.reject(r);
+				}
+			}), n.push({
+				v: e,
+				d: o,
+				a: r
+			});
+		} else r && n.push({
+			d: e,
+			a: r
+		});
+		return e;
+	}
+	return {
+		e,
+		u: using.bind(null, !1),
+		a: using.bind(null, !0),
+		d: function d() {
+			var o, t = this.e, s = 0;
+			function next() {
+				for (; o = n.pop();) try {
+					if (!o.a && 1 === s) return s = 0, n.push(o), Promise.resolve().then(next);
+					if (o.d) {
+						var r = o.d.call(o.v);
+						if (o.a) return s |= 2, Promise.resolve(r).then(next, err);
+					} else s |= 1;
+				} catch (r) {
+					return err(r);
+				}
+				if (1 === s) return t !== e ? Promise.reject(t) : Promise.resolve();
+				if (t !== e) throw t;
+			}
+			function err(n) {
+				return t = t !== e ? new r(n, t) : n, next();
+			}
+			return next();
+		}
+	};
+}
+// src/foreignUnfoldFns.ts
 /** Default neutral-axis K-factor for a purely geometric unfold (mid-surface). */
 var DEFAULT_FOREIGN_K = .5;
 /** Relative tolerance for the cylinder-fit residual (fraction of the radius). */
@@ -1602,10 +1651,11 @@ function unfoldForeignSolid(solid, opts) {
 			message: `UNSUPPORTED_FACE: surface type '${t.value}' is not a planar panel or cylindrical bend; the unfold ignores it and may be incomplete`
 		});
 	}
-	const thicknessResult = detectThickness(planarFaces);
+	const planar = planarFaces.map(describePlanar);
+	const thicknessResult = detectThickness(planar);
 	if (!thicknessResult.ok) return thicknessResult;
 	const thickness = thicknessResult.value;
-	const flats = pairFlats(planarFaces, thickness, warnings);
+	const flats = pairFlats(planar, thickness, warnings);
 	if (flats.length === 0) return (0, brepjs.err)((0, brepjs.validationError)("NO_FLATS", "no planar panel pairs detected; not a recognizable sheet-metal solid"));
 	const layout = layoutForeign(flats, buildBendGraph(flats, pairBends(cylFaces, warnings), warnings), thickness, kFactor, warnings);
 	if (!layout.ok) return layout;
@@ -1615,59 +1665,106 @@ function unfoldForeignSolid(solid, opts) {
 		warnings
 	});
 }
+function describePlanar(face) {
+	const center = (0, brepjs.faceCenter)(face);
+	const n = 4;
+	const candidates = [center];
+	for (let i = 0; i < n; i += 1) for (let j = 0; j < n; j += 1) candidates.push((0, brepjs.pointOnSurface)(face, (i + .5) / n, (j + .5) / n));
+	const interior = candidates.filter((p) => liesOnFace(face, p));
+	return {
+		face,
+		normal: faceNormal(face),
+		center,
+		interior
+	};
+}
+/** True if `point` (already on the face's surface) lies within the face's trimmed boundary. */
+function liesOnFace(face, point) {
+	try {
+		var _usingCtx$1 = _usingCtx();
+		const probe = _usingCtx$1.u((0, brepjs.vertex)(point));
+		const d = (0, brepjs.measureDistance)(probe, face);
+		return d.ok && d.value < GEOM_TOL;
+	} catch (_) {
+		_usingCtx$1.e = _;
+	} finally {
+		_usingCtx$1.d();
+	}
+}
+/** Normal-direction gap between two planar faces, or `null` unless they are antiparallel. */
+function antiparallelGap(a, b) {
+	if ((0, brepjs.vecDot)(a.normal, b.normal) > -.99) return null;
+	return Math.abs((0, brepjs.vecDot)((0, brepjs.vecSub)(b.center, a.center), a.normal));
+}
 /**
-* Detect the sheet thickness as the smallest gap between an antiparallel pair of
-* planar faces (the two faces of one flat panel). For a roughly-uniform-thickness
-* part this gap is the same across panels; a spread beyond tolerance warns.
+* Antiparallel faces are the two sides of one panel only if they face each other:
+* an interior point of one projects inside the other. The normal-direction gap
+* alone also matches faces that merely sit on parallel planes, e.g. a flange's end
+* cap and the adjacent flange's inner face across a corner. Their gap is the bend's
+* inner radius, which passes for the sheet thickness as T approaches or exceeds R.
+*/
+function facesOppose(a, b) {
+	return b.interior.some((p) => liesOnFace(a.face, projectToPlane(p, a.normal, a.center))) || a.interior.some((p) => liesOnFace(b.face, projectToPlane(p, b.normal, b.center)));
+}
+/**
+* Detect the sheet thickness as the smallest gap between two antiparallel planar
+* faces that oppose each other (the two faces of one flat panel).
 */
 function detectThickness(planar) {
-	const gaps = [];
+	const pairs = [];
 	for (let i = 0; i < planar.length; i += 1) for (let j = i + 1; j < planar.length; j += 1) {
-		const fa = planar[i];
-		const fb = planar[j];
-		if (fa === void 0 || fb === void 0) continue;
-		const na = faceNormal(fa);
-		const nb = faceNormal(fb);
-		if ((0, brepjs.vecDot)(na, nb) > -.99) continue;
-		const ca = (0, brepjs.faceCenter)(fa);
-		const cb = (0, brepjs.faceCenter)(fb);
-		const gap = Math.abs((0, brepjs.vecDot)((0, brepjs.vecSub)(cb, ca), na));
-		if (gap > GEOM_TOL) gaps.push(gap);
+		const a = planar[i];
+		const b = planar[j];
+		if (a === void 0 || b === void 0) continue;
+		const gap = antiparallelGap(a, b);
+		if (gap !== null && gap > GEOM_TOL) pairs.push({
+			gap,
+			a,
+			b
+		});
+	}
+	pairs.sort((p, q) => p.gap - q.gap);
+	const thinnest = pairs.find((p) => facesOppose(p.a, p.b));
+	if (thinnest === void 0) return (0, brepjs.err)((0, brepjs.validationError)("NO_THICKNESS", "could not detect a sheet thickness (no opposed planar pair)"));
+	return (0, brepjs.ok)(thinnest.gap);
+}
+/** Unused face opposing `a` whose gap is nearest the sheet thickness, within tolerance. */
+function bestOpposedPartner(planar, i, thickness, used) {
+	const a = planar[i];
+	if (a === void 0) return null;
+	const candidates = [];
+	for (let j = 0; j < planar.length; j += 1) {
+		const b = planar[j];
+		if (j === i || used.has(j) || b === void 0) continue;
+		const gap = antiparallelGap(a, b);
+		if (gap === null) continue;
+		const e = Math.abs(gap - thickness);
+		if (e <= thickness * .5 + GEOM_TOL) candidates.push({
+			j,
+			e
+		});
 	}
-	if (gaps.length === 0) return (0, brepjs.err)((0, brepjs.validationError)("NO_THICKNESS", "could not detect a sheet thickness (no opposed planar pair)"));
-	gaps.sort((a, b) => a - b);
-	return (0, brepjs.ok)(gaps[0]);
+	candidates.sort((p, q) => p.e - q.e);
+	const match = candidates.find((c) => {
+		const b = planar[c.j];
+		return b !== void 0 && facesOppose(a, b);
+	});
+	return match === void 0 ? null : match.j;
 }
 /**
-* Pair planar faces into flat panels: each panel is two parallel faces a
-* `thickness` apart with overlapping footprints. Greedy nearest-match by gap.
+* Pair planar faces into flat panels: each panel is two antiparallel faces a
+* `thickness` apart that oppose each other. Greedy nearest-match by gap.
 */
 function pairFlats(planar, thickness, warnings) {
 	const used = /* @__PURE__ */ new Set();
 	const flats = [];
 	for (let i = 0; i < planar.length; i += 1) {
 		if (used.has(i)) continue;
-		const fa = planar[i];
-		if (fa === void 0) continue;
-		const na = faceNormal(fa);
-		const ca = (0, brepjs.faceCenter)(fa);
-		let best = -1;
-		let bestErr = Infinity;
-		for (let j = 0; j < planar.length; j += 1) {
-			if (j === i || used.has(j)) continue;
-			const fb = planar[j];
-			if (fb === void 0) continue;
-			const nb = faceNormal(fb);
-			if ((0, brepjs.vecDot)(na, nb) > -.99) continue;
-			const cb = (0, brepjs.faceCenter)(fb);
-			const gap = Math.abs((0, brepjs.vecDot)((0, brepjs.vecSub)(cb, ca), na));
-			const e = Math.abs(gap - thickness);
-			if (e < bestErr) {
-				bestErr = e;
-				best = j;
-			}
-		}
-		if (best < 0 || bestErr > thickness * .5 + GEOM_TOL) {
+		const pa = planar[i];
+		if (pa === void 0) continue;
+		const fa = pa.face;
+		const best = bestOpposedPartner(planar, i, thickness, used);
+		if (best === null) {
 			const db = (0, brepjs.getBounds)(fa);
 			const ext = [
 				db.xMax - db.xMin,
@@ -1681,7 +1778,9 @@ function pairFlats(planar, thickness, warnings) {
 			});
 			continue;
 		}
-		const fb = planar[best];
+		const pb = planar[best];
+		if (pb === void 0) continue;
+		const fb = pb.face;
 		used.add(i);
 		used.add(best);
 		const areaA = (0, brepjs.measureArea)(fa);
@@ -1704,8 +1803,8 @@ function pairFlats(planar, thickness, warnings) {
 		flats.push({
 			id: `flat-${flats.length}`,
 			faces: [fa, fb],
-			normal: na,
-			center: (0, brepjs.vecScale)((0, brepjs.vecAdd)(ca, (0, brepjs.faceCenter)(fb)), .5),
+			normal: pa.normal,
+			center: (0, brepjs.vecScale)((0, brepjs.vecAdd)(pa.center, pb.center), .5),
 			area
 		});
 	}
@@ -2181,8 +2280,7 @@ function bendDirection(bend, parent) {
 function faceNormal(face) {
 	return (0, brepjs.vecNormalize)((0, brepjs.normalAt)(face, (0, brepjs.faceCenter)(face)));
 }
-//#endregion
-//#region src/authorFns.ts
+// src/authorFns.ts
 /**
 * Stable edge reference for flange attachment. The base flat is `face-0`; every
 * flange face is `face-<n+1>` in authoring order. The leading edge of a flat is
@@ -2799,8 +2897,7 @@ function buildFlangeFlat(axisOrigin, thetaDeg, width, thickness, length, directi
 		]
 	});
 }
-//#endregion
-//#region src/reliefFns.ts
+// src/reliefFns.ts
 var EPS$6 = 1e-6;
 var Z_AXIS$1 = [
 	0,
@@ -3164,8 +3261,7 @@ function toPt2$1(v) {
 function dist2(a, b) {
 	return Math.hypot(a[0] - b[0], a[1] - b[1]);
 }
-//#endregion
-//#region src/cutoutFns.ts
+// src/cutoutFns.ts
 var EPS$5 = 1e-6;
 /** A hair of over-cut so the tool pokes through both faces cleanly. */
 var HAIR$1 = .05;
@@ -3393,8 +3489,7 @@ function polygonArea(loop) {
 	}
 	return Math.abs(sum) / 2;
 }
-//#endregion
-//#region src/tabFns.ts
+// src/tabFns.ts
 var EPS$4 = 1e-6;
 /**
 * Fuse a rectangular tab onto a named flat region's edge — additive material (the
@@ -3536,8 +3631,7 @@ function devRect(local, dev) {
 		y1
 	];
 }
-//#endregion
-//#region src/formFns.ts
+// src/formFns.ts
 var EPS$3 = 1e-6;
 var HAIR = .05;
 var CIRCLE_SEGMENTS = 48;
@@ -3838,8 +3932,7 @@ function circleLocal(c, r) {
 	}
 	return pts;
 }
-//#endregion
-//#region src/validateFns.ts
+// src/validateFns.ts
 var OVERLAP_TOL = 1e-6;
 /**
 * Manufacturability checks for an authored sheet-metal part. Returns a list of
@@ -4021,8 +4114,7 @@ function pairKey(a, b) {
 function boundsOverlap(a, b) {
 	return a.xMin < b.xMax - OVERLAP_TOL && b.xMin < a.xMax - OVERLAP_TOL && a.yMin < b.yMax - OVERLAP_TOL && b.yMin < a.yMax - OVERLAP_TOL && a.zMin < b.zMax - OVERLAP_TOL && b.zMin < a.zMax - OVERLAP_TOL;
 }
-//#endregion
-//#region src/foldFns.ts
+// src/foldFns.ts
 /**
 * Fold a flat pattern up into a 3D part. Each {@link FoldRegion} folds about its
 * fold line off its parent region, exactly inverting {@link unfold}: folding by an
@@ -4680,8 +4772,7 @@ function uniqueSorted(values) {
 	}
 	return out;
 }
-//#endregion
-//#region src/miterFns.ts
+// src/miterFns.ts
 /**
 * General miter primitive: cut a sheet-metal part by a plane, removing the
 * material on the `+normal` side. The tool is a half-space block sized to the
@@ -4803,8 +4894,7 @@ function alignZTo(shape, target) {
 		axis
 	});
 }
-//#endregion
-//#region src/contourFlangeFns.ts
+// src/contourFlangeFns.ts
 /**
 * Author a contour flange: an OPEN 2D profile (alternating line/arc segments) swept
 * along one straight edge of the base flat. Each arc is a cylindrical bend, each
@@ -5273,8 +5363,7 @@ function matrixToAxisAngle(m) {
 		angleDeg
 	};
 }
-//#endregion
-//#region src/loftedFlangeFns.ts
+// src/loftedFlangeFns.ts
 /**
 * Relative out-of-plane tolerance for the per-quad developability test: a quad is
 * twisted (non-developable) when its far corner sits off the diagonal-triangle plane
@@ -5586,8 +5675,7 @@ function closedLoopWire(loop) {
 	}
 	return (0, brepjs.wireLoop)(edges);
 }
-//#endregion
-//#region src/hemFns.ts
+// src/hemFns.ts
 /**
 * A small clearance baked into a CLOSED hem's curl radius. A true closed hem folds
 * the return flat against the parent with zero gap; modelling that exactly makes the
@@ -5780,8 +5868,7 @@ function splitAngle(totalDeg) {
 	const per = totalDeg / n;
 	return Array.from({ length: n }, () => per);
 }
-//#endregion
-//#region src/jogFns.ts
+// src/jogFns.ts
 /**
 * Author a jog (joggle): two opposite bends (`+θ` then `−θ`) that step a flat by
 * `offsetHeight` perpendicular to its plane, then continue parallel past the step.
@@ -5940,8 +6027,7 @@ function jogId(spec) {
 function resolveRegion(region) {
 	return region === "base" || region === "face-0" ? ROOT_FLAT_ID : region;
 }
-//#endregion
-//#region src/dxfFns.ts
+// src/dxfFns.ts
 var DEFAULT_TEXT_HEIGHT = 2.5;
 var LAYER_OUTLINE = "OUTLINE";
 var LAYER_BEND_UP = "BEND_UP";
@@ -6150,8 +6236,7 @@ function pathPoints(wire) {
 function toPt2(v) {
 	return [v[0], v[1]];
 }
-//#endregion
-//#region src/polygonFns.ts
+// src/polygonFns.ts
 /**
 * 2D simple-polygon geometry for true-shape (no-fit-polygon) nesting: an ordered
 * vertex loop plus the overlap predicate that is the correctness backbone of the
@@ -6394,8 +6479,7 @@ function pointSegmentDistance(p, a, b) {
 	const dy = p[1] - (a[1] + t * vy);
 	return Math.hypot(dx, dy);
 }
-//#endregion
-//#region src/nestFns.ts
+// src/nestFns.ts
 var EPS = 1e-9;
 /**
 * Bounding-box nest: arrange developed flat patterns onto stock sheets to reduce
@@ -6861,8 +6945,7 @@ function nestToDXF$1(result, patterns, sheetIndex, options) {
 	}
 	return multiPatternToDXF(placed, options ?? {});
 }
-//#endregion
-//#region src/reportFns.ts
+// src/reportFns.ts
 /**
 * Build the bend-table JSON report directly from an authored part. Walks the
 * same tree layout the unfold consumes so the per-bend values and the total flat
@@ -6905,8 +6988,7 @@ function reportFromUnfold(result) {
 function reportToJSON(report) {
 	return JSON.stringify(report, null, 2);
 }
-//#endregion
-//#region src/api.ts
+// src/api.ts
 /** Author a straight-bend part: a base flat plus folded-up flanges. */
 function author(spec) {
 	return authorPart(spec);
@@ -7068,8 +7150,7 @@ function bendTable(id) {
 function resolveAllowance(rule, angleDeg, thickness, onWarning) {
 	return resolveBendAllowance(rule, angleDeg, thickness, onWarning);
 }
-//#endregion
-//#region src/facade.ts
+// src/facade.ts
 /** Thrown by the fluent facade when an underlying `Result<T>` is an `Err`. */
 var SheetMetalError = class extends Error {
 	code;
@@ -7339,7 +7420,6 @@ function foldFlat(input) {
 function fromSolid(solid, opts) {
 	return new ForeignSolidHandle(solid, opts);
 }
-//#endregion
 exports.ALUMINUM_GAUGES = ALUMINUM_GAUGES;
 exports.DEFAULT_K_FACTOR = DEFAULT_K_FACTOR;
 exports.ForeignSolidHandle = ForeignSolidHandle;

package/dist/brepjs-sheetmetal.js

@@ -1,5 +1,5 @@
-import { box, curveEndPoint, curveStartPoint, cut, cylinder, err, extrude, face, faceCenter, fuse, getBounds, getEdges, getFaces, getSolids, getSurfaceType, intersect, isErr, isOk, isPlanarWire, isSolid, isValid, line, measureArea, measureVolume, normalAt, ok, pointOnSurface, rotate, sharedEdges, translate, validationError, vecAdd, vecCross, vecDot, vecLength, vecNormalize, vecScale, vecSub, wire, wireLoop } from "brepjs";
-//#region src/materials.ts
+import { box, curveEndPoint, curveStartPoint, cut, cylinder, err, extrude, face, faceCenter, fuse, getBounds, getEdges, getFaces, getSolids, getSurfaceType, intersect, isErr, isOk, isPlanarWire, isSolid, isValid, line, measureArea, measureDistance, measureVolume, normalAt, ok, pointOnSurface, rotate, sharedEdges, translate, validationError, vecAdd, vecCross, vecDot, vecLength, vecNormalize, vecScale, vecSub, vertex, wire, wireLoop } from "brepjs";
+// src/materials.ts
 var DEFAULT_K_FACTOR = .44;
 function gauge(name, thickness, kFactor = DEFAULT_K_FACTOR) {
 	return {
@@ -37,8 +37,7 @@ var MATERIALS = {
 function getMaterial(name) {
 	return MATERIALS[name];
 }
-//#endregion
-//#region src/bendTableFns.ts
+// src/bendTableFns.ts
 var DEG_TO_RAD = Math.PI / 180;
 var registry = /* @__PURE__ */ new Map();
 /**
@@ -377,8 +376,7 @@ function ensureStarterTables() {
 		});
 	}
 }
-//#endregion
-//#region src/allowanceFns.ts
+// src/allowanceFns.ts
 function validateRule(rule) {
 	if (!Number.isFinite(rule.innerRadius) || rule.innerRadius < 0) return err(validationError("INVALID_RADIUS", `innerRadius must be a finite, non-negative number, got ${rule.innerRadius}`));
 	if (!Number.isFinite(rule.kFactor) || rule.kFactor < 0 || rule.kFactor > 1) return err(validationError("INVALID_K_FACTOR", `kFactor must be in [0, 1], got ${rule.kFactor}`));
@@ -412,8 +410,7 @@ function neutralRadius(thickness, rule) {
 	if (!ruleCheck.ok) return ruleCheck;
 	return ok(rule.innerRadius + rule.kFactor * thickness);
 }
-//#endregion
-//#region src/featureTreeFns.ts
+// src/featureTreeFns.ts
 var ROOT_FLAT_ID = "root";
 /**
 * Build the feature graph for an authored part. Each flange contributes a flat
@@ -558,8 +555,7 @@ function featureTree(part, rootId = ROOT_FLAT_ID) {
 	if (!graph.ok) return graph;
 	return buildFeatureTree(graph.value, rootId);
 }
-//#endregion
-//#region src/internal.ts
+// src/internal.ts
 /**
 * OCCT boolean fuse returns a compound wrapping the unioned solid. Extract the
 * single solid so a part's cached shape is a true `Solid` (clean STEP/GLB
@@ -582,8 +578,7 @@ function sideFromAxisDir(axisDir) {
 	if (Math.abs(axisDir[1]) > .5) return "xmax";
 	if (Math.abs(axisDir[0]) > .5) return "ymax";
 }
-//#endregion
-//#region src/unfoldFns.ts
+// src/unfoldFns.ts
 /**
 * Tree-driven analytic unfold of an arbitrary straight-bend part into its flat
 * pattern. Walks the feature tree in BFS order placing every flat in the flat
@@ -1290,8 +1285,62 @@ function reportFromLayout(layout) {
 		totalFlatSize: [maxX - minX, maxY - minY]
 	};
 }
-//#endregion
-//#region src/foreignUnfoldFns.ts
+// \0@oxc-project+runtime@0.155.0/helpers/esm/usingCtx.js
+function _usingCtx() {
+	var r = "function" == typeof SuppressedError ? SuppressedError : function(r, e) {
+		var n = Error();
+		return n.name = "SuppressedError", n.error = r, n.suppressed = e, n;
+	}, e = {}, n = [];
+	function using(r, e) {
+		if (null != e) {
+			if (Object(e) !== e) throw new TypeError("using declarations can only be used with objects, functions, null, or undefined.");
+			if (r) var o = e[Symbol.asyncDispose || Symbol["for"]("Symbol.asyncDispose")];
+			if (void 0 === o && (o = e[Symbol.dispose || Symbol["for"]("Symbol.dispose")], r)) var t = o;
+			if ("function" != typeof o) throw new TypeError("Object is not disposable.");
+			t && (o = function o() {
+				try {
+					t.call(e);
+				} catch (r) {
+					return Promise.reject(r);
+				}
+			}), n.push({
+				v: e,
+				d: o,
+				a: r
+			});
+		} else r && n.push({
+			d: e,
+			a: r
+		});
+		return e;
+	}
+	return {
+		e,
+		u: using.bind(null, !1),
+		a: using.bind(null, !0),
+		d: function d() {
+			var o, t = this.e, s = 0;
+			function next() {
+				for (; o = n.pop();) try {
+					if (!o.a && 1 === s) return s = 0, n.push(o), Promise.resolve().then(next);
+					if (o.d) {
+						var r = o.d.call(o.v);
+						if (o.a) return s |= 2, Promise.resolve(r).then(next, err);
+					} else s |= 1;
+				} catch (r) {
+					return err(r);
+				}
+				if (1 === s) return t !== e ? Promise.reject(t) : Promise.resolve();
+				if (t !== e) throw t;
+			}
+			function err(n) {
+				return t = t !== e ? new r(n, t) : n, next();
+			}
+			return next();
+		}
+	};
+}
+// src/foreignUnfoldFns.ts
 /** Default neutral-axis K-factor for a purely geometric unfold (mid-surface). */
 var DEFAULT_FOREIGN_K = .5;
 /** Relative tolerance for the cylinder-fit residual (fraction of the radius). */
@@ -1601,10 +1650,11 @@ function unfoldForeignSolid(solid, opts) {
 			message: `UNSUPPORTED_FACE: surface type '${t.value}' is not a planar panel or cylindrical bend; the unfold ignores it and may be incomplete`
 		});
 	}
-	const thicknessResult = detectThickness(planarFaces);
+	const planar = planarFaces.map(describePlanar);
+	const thicknessResult = detectThickness(planar);
 	if (!thicknessResult.ok) return thicknessResult;
 	const thickness = thicknessResult.value;
-	const flats = pairFlats(planarFaces, thickness, warnings);
+	const flats = pairFlats(planar, thickness, warnings);
 	if (flats.length === 0) return err(validationError("NO_FLATS", "no planar panel pairs detected; not a recognizable sheet-metal solid"));
 	const layout = layoutForeign(flats, buildBendGraph(flats, pairBends(cylFaces, warnings), warnings), thickness, kFactor, warnings);
 	if (!layout.ok) return layout;
@@ -1614,59 +1664,106 @@ function unfoldForeignSolid(solid, opts) {
 		warnings
 	});
 }
+function describePlanar(face) {
+	const center = faceCenter(face);
+	const n = 4;
+	const candidates = [center];
+	for (let i = 0; i < n; i += 1) for (let j = 0; j < n; j += 1) candidates.push(pointOnSurface(face, (i + .5) / n, (j + .5) / n));
+	const interior = candidates.filter((p) => liesOnFace(face, p));
+	return {
+		face,
+		normal: faceNormal(face),
+		center,
+		interior
+	};
+}
+/** True if `point` (already on the face's surface) lies within the face's trimmed boundary. */
+function liesOnFace(face, point) {
+	try {
+		var _usingCtx$1 = _usingCtx();
+		const probe = _usingCtx$1.u(vertex(point));
+		const d = measureDistance(probe, face);
+		return d.ok && d.value < GEOM_TOL;
+	} catch (_) {
+		_usingCtx$1.e = _;
+	} finally {
+		_usingCtx$1.d();
+	}
+}
+/** Normal-direction gap between two planar faces, or `null` unless they are antiparallel. */
+function antiparallelGap(a, b) {
+	if (vecDot(a.normal, b.normal) > -.99) return null;
+	return Math.abs(vecDot(vecSub(b.center, a.center), a.normal));
+}
 /**
-* Detect the sheet thickness as the smallest gap between an antiparallel pair of
-* planar faces (the two faces of one flat panel). For a roughly-uniform-thickness
-* part this gap is the same across panels; a spread beyond tolerance warns.
+* Antiparallel faces are the two sides of one panel only if they face each other:
+* an interior point of one projects inside the other. The normal-direction gap
+* alone also matches faces that merely sit on parallel planes, e.g. a flange's end
+* cap and the adjacent flange's inner face across a corner. Their gap is the bend's
+* inner radius, which passes for the sheet thickness as T approaches or exceeds R.
+*/
+function facesOppose(a, b) {
+	return b.interior.some((p) => liesOnFace(a.face, projectToPlane(p, a.normal, a.center))) || a.interior.some((p) => liesOnFace(b.face, projectToPlane(p, b.normal, b.center)));
+}
+/**
+* Detect the sheet thickness as the smallest gap between two antiparallel planar
+* faces that oppose each other (the two faces of one flat panel).
 */
 function detectThickness(planar) {
-	const gaps = [];
+	const pairs = [];
 	for (let i = 0; i < planar.length; i += 1) for (let j = i + 1; j < planar.length; j += 1) {
-		const fa = planar[i];
-		const fb = planar[j];
-		if (fa === void 0 || fb === void 0) continue;
-		const na = faceNormal(fa);
-		const nb = faceNormal(fb);
-		if (vecDot(na, nb) > -.99) continue;
-		const ca = faceCenter(fa);
-		const cb = faceCenter(fb);
-		const gap = Math.abs(vecDot(vecSub(cb, ca), na));
-		if (gap > GEOM_TOL) gaps.push(gap);
+		const a = planar[i];
+		const b = planar[j];
+		if (a === void 0 || b === void 0) continue;
+		const gap = antiparallelGap(a, b);
+		if (gap !== null && gap > GEOM_TOL) pairs.push({
+			gap,
+			a,
+			b
+		});
+	}
+	pairs.sort((p, q) => p.gap - q.gap);
+	const thinnest = pairs.find((p) => facesOppose(p.a, p.b));
+	if (thinnest === void 0) return err(validationError("NO_THICKNESS", "could not detect a sheet thickness (no opposed planar pair)"));
+	return ok(thinnest.gap);
+}
+/** Unused face opposing `a` whose gap is nearest the sheet thickness, within tolerance. */
+function bestOpposedPartner(planar, i, thickness, used) {
+	const a = planar[i];
+	if (a === void 0) return null;
+	const candidates = [];
+	for (let j = 0; j < planar.length; j += 1) {
+		const b = planar[j];
+		if (j === i || used.has(j) || b === void 0) continue;
+		const gap = antiparallelGap(a, b);
+		if (gap === null) continue;
+		const e = Math.abs(gap - thickness);
+		if (e <= thickness * .5 + GEOM_TOL) candidates.push({
+			j,
+			e
+		});
 	}
-	if (gaps.length === 0) return err(validationError("NO_THICKNESS", "could not detect a sheet thickness (no opposed planar pair)"));
-	gaps.sort((a, b) => a - b);
-	return ok(gaps[0]);
+	candidates.sort((p, q) => p.e - q.e);
+	const match = candidates.find((c) => {
+		const b = planar[c.j];
+		return b !== void 0 && facesOppose(a, b);
+	});
+	return match === void 0 ? null : match.j;
 }
 /**
-* Pair planar faces into flat panels: each panel is two parallel faces a
-* `thickness` apart with overlapping footprints. Greedy nearest-match by gap.
+* Pair planar faces into flat panels: each panel is two antiparallel faces a
+* `thickness` apart that oppose each other. Greedy nearest-match by gap.
 */
 function pairFlats(planar, thickness, warnings) {
 	const used = /* @__PURE__ */ new Set();
 	const flats = [];
 	for (let i = 0; i < planar.length; i += 1) {
 		if (used.has(i)) continue;
-		const fa = planar[i];
-		if (fa === void 0) continue;
-		const na = faceNormal(fa);
-		const ca = faceCenter(fa);
-		let best = -1;
-		let bestErr = Infinity;
-		for (let j = 0; j < planar.length; j += 1) {
-			if (j === i || used.has(j)) continue;
-			const fb = planar[j];
-			if (fb === void 0) continue;
-			const nb = faceNormal(fb);
-			if (vecDot(na, nb) > -.99) continue;
-			const cb = faceCenter(fb);
-			const gap = Math.abs(vecDot(vecSub(cb, ca), na));
-			const e = Math.abs(gap - thickness);
-			if (e < bestErr) {
-				bestErr = e;
-				best = j;
-			}
-		}
-		if (best < 0 || bestErr > thickness * .5 + GEOM_TOL) {
+		const pa = planar[i];
+		if (pa === void 0) continue;
+		const fa = pa.face;
+		const best = bestOpposedPartner(planar, i, thickness, used);
+		if (best === null) {
 			const db = getBounds(fa);
 			const ext = [
 				db.xMax - db.xMin,
@@ -1680,7 +1777,9 @@ function pairFlats(planar, thickness, warnings) {
 			});
 			continue;
 		}
-		const fb = planar[best];
+		const pb = planar[best];
+		if (pb === void 0) continue;
+		const fb = pb.face;
 		used.add(i);
 		used.add(best);
 		const areaA = measureArea(fa);
@@ -1703,8 +1802,8 @@ function pairFlats(planar, thickness, warnings) {
 		flats.push({
 			id: `flat-${flats.length}`,
 			faces: [fa, fb],
-			normal: na,
-			center: vecScale(vecAdd(ca, faceCenter(fb)), .5),
+			normal: pa.normal,
+			center: vecScale(vecAdd(pa.center, pb.center), .5),
 			area
 		});
 	}
@@ -2180,8 +2279,7 @@ function bendDirection(bend, parent) {
 function faceNormal(face) {
 	return vecNormalize(normalAt(face, faceCenter(face)));
 }
-//#endregion
-//#region src/authorFns.ts
+// src/authorFns.ts
 /**
 * Stable edge reference for flange attachment. The base flat is `face-0`; every
 * flange face is `face-<n+1>` in authoring order. The leading edge of a flat is
@@ -2798,8 +2896,7 @@ function buildFlangeFlat(axisOrigin, thetaDeg, width, thickness, length, directi
 		]
 	});
 }
-//#endregion
-//#region src/reliefFns.ts
+// src/reliefFns.ts
 var EPS$6 = 1e-6;
 var Z_AXIS$1 = [
 	0,
@@ -3163,8 +3260,7 @@ function toPt2$1(v) {
 function dist2(a, b) {
 	return Math.hypot(a[0] - b[0], a[1] - b[1]);
 }
-//#endregion
-//#region src/cutoutFns.ts
+// src/cutoutFns.ts
 var EPS$5 = 1e-6;
 /** A hair of over-cut so the tool pokes through both faces cleanly. */
 var HAIR$1 = .05;
@@ -3392,8 +3488,7 @@ function polygonArea(loop) {
 	}
 	return Math.abs(sum) / 2;
 }
-//#endregion
-//#region src/tabFns.ts
+// src/tabFns.ts
 var EPS$4 = 1e-6;
 /**
 * Fuse a rectangular tab onto a named flat region's edge — additive material (the
@@ -3535,8 +3630,7 @@ function devRect(local, dev) {
 		y1
 	];
 }
-//#endregion
-//#region src/formFns.ts
+// src/formFns.ts
 var EPS$3 = 1e-6;
 var HAIR = .05;
 var CIRCLE_SEGMENTS = 48;
@@ -3837,8 +3931,7 @@ function circleLocal(c, r) {
 	}
 	return pts;
 }
-//#endregion
-//#region src/validateFns.ts
+// src/validateFns.ts
 var OVERLAP_TOL = 1e-6;
 /**
 * Manufacturability checks for an authored sheet-metal part. Returns a list of
@@ -4020,8 +4113,7 @@ function pairKey(a, b) {
 function boundsOverlap(a, b) {
 	return a.xMin < b.xMax - OVERLAP_TOL && b.xMin < a.xMax - OVERLAP_TOL && a.yMin < b.yMax - OVERLAP_TOL && b.yMin < a.yMax - OVERLAP_TOL && a.zMin < b.zMax - OVERLAP_TOL && b.zMin < a.zMax - OVERLAP_TOL;
 }
-//#endregion
-//#region src/foldFns.ts
+// src/foldFns.ts
 /**
 * Fold a flat pattern up into a 3D part. Each {@link FoldRegion} folds about its
 * fold line off its parent region, exactly inverting {@link unfold}: folding by an
@@ -4679,8 +4771,7 @@ function uniqueSorted(values) {
 	}
 	return out;
 }
-//#endregion
-//#region src/miterFns.ts
+// src/miterFns.ts
 /**
 * General miter primitive: cut a sheet-metal part by a plane, removing the
 * material on the `+normal` side. The tool is a half-space block sized to the
@@ -4802,8 +4893,7 @@ function alignZTo(shape, target) {
 		axis
 	});
 }
-//#endregion
-//#region src/contourFlangeFns.ts
+// src/contourFlangeFns.ts
 /**
 * Author a contour flange: an OPEN 2D profile (alternating line/arc segments) swept
 * along one straight edge of the base flat. Each arc is a cylindrical bend, each
@@ -5272,8 +5362,7 @@ function matrixToAxisAngle(m) {
 		angleDeg
 	};
 }
-//#endregion
-//#region src/loftedFlangeFns.ts
+// src/loftedFlangeFns.ts
 /**
 * Relative out-of-plane tolerance for the per-quad developability test: a quad is
 * twisted (non-developable) when its far corner sits off the diagonal-triangle plane
@@ -5585,8 +5674,7 @@ function closedLoopWire(loop) {
 	}
 	return wireLoop(edges);
 }
-//#endregion
-//#region src/hemFns.ts
+// src/hemFns.ts
 /**
 * A small clearance baked into a CLOSED hem's curl radius. A true closed hem folds
 * the return flat against the parent with zero gap; modelling that exactly makes the
@@ -5779,8 +5867,7 @@ function splitAngle(totalDeg) {
 	const per = totalDeg / n;
 	return Array.from({ length: n }, () => per);
 }
-//#endregion
-//#region src/jogFns.ts
+// src/jogFns.ts
 /**
 * Author a jog (joggle): two opposite bends (`+θ` then `−θ`) that step a flat by
 * `offsetHeight` perpendicular to its plane, then continue parallel past the step.
@@ -5939,8 +6026,7 @@ function jogId(spec) {
 function resolveRegion(region) {
 	return region === "base" || region === "face-0" ? ROOT_FLAT_ID : region;
 }
-//#endregion
-//#region src/dxfFns.ts
+// src/dxfFns.ts
 var DEFAULT_TEXT_HEIGHT = 2.5;
 var LAYER_OUTLINE = "OUTLINE";
 var LAYER_BEND_UP = "BEND_UP";
@@ -6149,8 +6235,7 @@ function pathPoints(wire) {
 function toPt2(v) {
 	return [v[0], v[1]];
 }
-//#endregion
-//#region src/polygonFns.ts
+// src/polygonFns.ts
 /**
 * 2D simple-polygon geometry for true-shape (no-fit-polygon) nesting: an ordered
 * vertex loop plus the overlap predicate that is the correctness backbone of the
@@ -6393,8 +6478,7 @@ function pointSegmentDistance(p, a, b) {
 	const dy = p[1] - (a[1] + t * vy);
 	return Math.hypot(dx, dy);
 }
-//#endregion
-//#region src/nestFns.ts
+// src/nestFns.ts
 var EPS = 1e-9;
 /**
 * Bounding-box nest: arrange developed flat patterns onto stock sheets to reduce
@@ -6860,8 +6944,7 @@ function nestToDXF$1(result, patterns, sheetIndex, options) {
 	}
 	return multiPatternToDXF(placed, options ?? {});
 }
-//#endregion
-//#region src/reportFns.ts
+// src/reportFns.ts
 /**
 * Build the bend-table JSON report directly from an authored part. Walks the
 * same tree layout the unfold consumes so the per-bend values and the total flat
@@ -6904,8 +6987,7 @@ function reportFromUnfold(result) {
 function reportToJSON(report) {
 	return JSON.stringify(report, null, 2);
 }
-//#endregion
-//#region src/api.ts
+// src/api.ts
 /** Author a straight-bend part: a base flat plus folded-up flanges. */
 function author(spec) {
 	return authorPart(spec);
@@ -7067,8 +7149,7 @@ function bendTable(id) {
 function resolveAllowance(rule, angleDeg, thickness, onWarning) {
 	return resolveBendAllowance(rule, angleDeg, thickness, onWarning);
 }
-//#endregion
-//#region src/facade.ts
+// src/facade.ts
 /** Thrown by the fluent facade when an underlying `Result<T>` is an `Err`. */
 var SheetMetalError = class extends Error {
 	code;
@@ -7338,5 +7419,4 @@ function foldFlat(input) {
 function fromSolid(solid, opts) {
 	return new ForeignSolidHandle(solid, opts);
 }
-//#endregion
 export { ALUMINUM_GAUGES, DEFAULT_K_FACTOR, ForeignSolidHandle, MATERIALS, ROOT_FLAT_ID, STEEL_GAUGES, SheetMetalBuilder, SheetMetalError, SheetMetalPartHandle, addBendRelief, addBendTable, addCutout, addForm, addHole, addPolygonCutout, addSlot, addTab, allowance, author, authorContourFlange, authorLoftedFlange, authorPart, autoBendReliefs, autoMiterCorner, autoReliefs, baseEdgeRef, bendAllowance, bendRelief, bendTable, buildFeatureGraph, buildFeatureTree, buildReport, contourFlange, cornerRelief, developed, developedLength, emboss, featureTree, fitCylinder, flatPatternToDXF, flatPatternToPolylines, fold, foldFlat, fold$1 as foldPart, foldWithWarnings, fromPart, fromSolid, getBendTable, getMaterial, hem, jog, loftedFlange, louver, miter, miterCorner, miterCut, multiPatternToDXF, nest, nestToDXF, neutralRadius, partToFlatInput, patternBbox, patternToFlatInput, pointInPolygon, polygonBounds, polygonsOverlap, polygonsOverlapWithClearance, registerBendTable, relieveCorner, report, reportFrom, reportFromUnfold, reportJSON, reportToJSON, resolveAllowance, resolveBendAllowance, segmentsIntersect, sheetMetal, signedArea, tabAndSlot, toDXF, transformPolygon, unfold, unfoldForeignSolid, unfold$1 as unfoldPart, unfoldSolid, validate, validatePart, wireToPolygon };

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "brepjs-sheetmetal",
-  "version": "0.3.1",
+  "version": "0.4.1",
   "description": "Sheet-metal CAD domain for brepjs — flange authoring, auto-miter, fold/unfold flat patterns",
   "keywords": [
     "sheet-metal",
@@ -53,13 +53,13 @@
   },
   "devDependencies": {
     "@eslint/js": "*",
-    "@types/node": "26.2.0",
+    "@types/node": "26.6.5",
     "eslint": "*",
-    "tsx": "4.23.11",
+    "tsx": "4.23.15",
     "typescript": "*",
     "typescript-eslint": "*",
     "vite": "^8.0.0",
     "vite-plugin-dts": "^5.0.0",
-    "vitest": "^4.0.0"
+    "vitest": "^5.0.0"
   }
 }