brepjs-sheetmetal 0.3.0 → 0.4.0

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
+// 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
@@ -504,7 +501,7 @@ function buildFeatureTree(graph, rootId = ROOT_FLAT_ID) {
 		adjacency.get(edge.parent)?.push(edge);
 		adjacency.get(edge.child)?.push(edge);
 	}
-	const visited = new Set([rootId]);
+	const visited = /* @__PURE__ */ new Set([rootId]);
 	const treeEdges = /* @__PURE__ */ new Set();
 	const bends = [];
 	const queue = [rootId];
@@ -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
@@ -673,18 +668,20 @@ function unfold(part) {
 		direction: p.direction,
 		inward: [-p.frame.v[0], -p.frame.v[1]]
 	})).concat(contourResult.value.bendLines).concat(hemJogResult.value.bendLines);
+	const pattern = {
+		outline: outlineResult.value,
+		bendLines,
+		holes: holesResult.value,
+		formCuts: formsResult.value.cuts,
+		formMarkers: formsResult.value.markers,
+		formHinges: formsResult.value.hinges,
+		loftedDevelopments: loftedResult.value,
+		developedArea: layout.developedArea
+	};
+	const report = reportFromLayout(layout);
 	return ok({
-		pattern: {
-			outline: outlineResult.value,
-			bendLines,
-			holes: holesResult.value,
-			formCuts: formsResult.value.cuts,
-			formMarkers: formsResult.value.markers,
-			formHinges: formsResult.value.hinges,
-			loftedDevelopments: loftedResult.value,
-			developedArea: layout.developedArea
-		},
-		report: reportFromLayout(layout),
+		pattern,
+		report,
 		warnings
 	});
 }
@@ -1288,8 +1285,7 @@ function reportFromLayout(layout) {
 		totalFlatSize: [maxX - minX, maxY - minY]
 	};
 }
-//#endregion
-//#region src/foreignUnfoldFns.ts
+// 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). */
@@ -1307,7 +1303,8 @@ function sampleFace(face, n) {
 	for (let i = 0; i < n; i += 1) {
 		const u = n === 1 ? .5 : i / (n - 1);
 		for (let j = 0; j < n; j += 1) {
-			const point = pointOnSurface(face, u, n === 1 ? .5 : j / (n - 1));
+			const v = n === 1 ? .5 : j / (n - 1);
+			const point = pointOnSurface(face, u, v);
 			const normal = normalAt(face, point);
 			const len = vecLength(normal);
 			if (len < GEOM_TOL) continue;
@@ -1398,7 +1395,7 @@ function axisFromNormals(normals) {
 }
 /** An orthonormal pair spanning the plane perpendicular to `axis`. */
 function perpBasis(axis) {
-	const e1 = vecNormalize(vecCross(axis, Math.abs(axis[0]) < .9 ? [
+	const seed = Math.abs(axis[0]) < .9 ? [
 		1,
 		0,
 		0
@@ -1406,7 +1403,8 @@ function perpBasis(axis) {
 		0,
 		1,
 		0
-	]));
+	];
+	const e1 = vecNormalize(vecCross(axis, seed));
 	return {
 		e1,
 		e2: vecNormalize(vecCross(axis, e1))
@@ -1415,7 +1413,9 @@ function perpBasis(axis) {
 /** Perpendicular distance from `p` to the line through `o` with unit direction `d`. */
 function distanceToAxis(p, o, d) {
 	const w = vecSub(p, o);
-	return vecLength(vecSub(w, vecScale(d, vecDot(w, d))));
+	const along = vecDot(w, d);
+	const perp = vecSub(w, vecScale(d, along));
+	return vecLength(perp);
 }
 /** Algebraic (Kåsa) least-squares circle fit to 2D points. */
 function fitCircle2(pts) {
@@ -1620,7 +1620,8 @@ function detectThickness(planar) {
 		const fb = planar[j];
 		if (fa === void 0 || fb === void 0) continue;
 		const na = faceNormal(fa);
-		if (vecDot(na, faceNormal(fb)) > -.99) continue;
+		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));
@@ -1649,7 +1650,8 @@ function pairFlats(planar, thickness, warnings) {
 			if (j === i || used.has(j)) continue;
 			const fb = planar[j];
 			if (fb === void 0) continue;
-			if (vecDot(na, faceNormal(fb)) > -.99) 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);
@@ -1818,7 +1820,7 @@ function layoutForeign(flats, graph, thickness, kFactor, warnings) {
 		adjacency.get(e.flats[0])?.push(e);
 		adjacency.get(e.flats[1])?.push(e);
 	}
-	const visited = new Set([0]);
+	const visited = /* @__PURE__ */ new Set([0]);
 	const treeEdges = /* @__PURE__ */ new Set();
 	const treeBends = [];
 	const queue = [0];
@@ -1912,23 +1914,26 @@ function layoutForeign(flats, graph, thickness, kFactor, warnings) {
 		});
 	}
 	const bbox = aabb2(rects);
-	const outlineWire = wireLoop(buildRectOutline2(bbox));
+	const outline = buildRectOutline2(bbox);
+	const outlineWire = wireLoop(outline);
 	if (!outlineWire.ok) return outlineWire;
+	const pattern = {
+		outline: outlineWire.value,
+		bendLines,
+		holes: [],
+		formCuts: [],
+		formMarkers: [],
+		formHinges: [],
+		loftedDevelopments: [],
+		developedArea
+	};
+	const report = {
+		bends: bendReportRows,
+		totalFlatSize: [bbox.x1 - bbox.x0, bbox.y1 - bbox.y0]
+	};
 	return ok({
-		pattern: {
-			outline: outlineWire.value,
-			bendLines,
-			holes: [],
-			formCuts: [],
-			formMarkers: [],
-			formHinges: [],
-			loftedDevelopments: [],
-			developedArea
-		},
-		report: {
-			bends: bendReportRows,
-			totalFlatSize: [bbox.x1 - bbox.x0, bbox.y1 - bbox.y0]
-		}
+		pattern,
+		report
 	});
 }
 /** Map a 3D in-plane point through a placed flat's affine map. */
@@ -2018,7 +2023,9 @@ function parentFootprintCentreAlong(parent, axisInPlane, axisRef) {
 		min = Math.min(min, v);
 		max = Math.max(max, v);
 	}
-	return vecAdd(axisRef, vecScale(axisInPlane, (min + max) / 2 - vecDot(axisRef, axisInPlane)));
+	const mid = (min + max) / 2;
+	const refAlong = vecDot(axisRef, axisInPlane);
+	return vecAdd(axisRef, vecScale(axisInPlane, mid - refAlong));
 }
 /** In-plane run extent of a flat along `out` (perpendicular to the bend axis). */
 function inPlaneRun(flat, out) {
@@ -2034,7 +2041,8 @@ function inPlaneRun(flat, out) {
 }
 /** Run direction of a child flat pointing away from the bend axis. */
 function runDirectionAway(flat, out, axisOrigin) {
-	return vecDot(vecSub(flat.center, axisOrigin), out) >= 0 ? out : vecScale(out, -1);
+	const toCenter = vecSub(flat.center, axisOrigin);
+	return vecDot(toCenter, out) >= 0 ? out : vecScale(out, -1);
 }
 /**
 * The child corner that develops onto the strip-far near point: minimum along the
@@ -2160,13 +2168,13 @@ function bendAxialExtent(bend) {
 * parent flat's mid-surface, 'down' otherwise.
 */
 function bendDirection(bend, parent) {
-	return vecDot(vecSub(bend.fit.axisOrigin, parent.center), parent.normal) >= 0 ? "up" : "down";
+	const toAxis = vecSub(bend.fit.axisOrigin, parent.center);
+	return vecDot(toAxis, parent.normal) >= 0 ? "up" : "down";
 }
 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
@@ -2494,10 +2502,10 @@ function parentEdge(f, side) {
 			length = f.uLen;
 			a = o;
 			b = uTop;
-			break;
 	}
 	const dir = vecNormalize(vecCross(f.n, out));
-	const start = vecDot(vecSub(b, a), dir) >= 0 ? a : b;
+	const toB = vecSub(b, a);
+	const start = vecDot(toB, dir) >= 0 ? a : b;
 	return {
 		dir,
 		out,
@@ -2610,14 +2618,15 @@ function frameTransform$1(axisT, runT, nT, edgeStart, offset) {
 	const applyR = (v) => vecAdd(vecAdd(vecScale(runT, v[0]), vecScale(axisT, v[1])), vecScale(nT, v[2]));
 	return {
 		solid: (s) => {
-			return translate(aa.angleDeg === 0 ? s : rotate(s, aa.angleDeg, {
+			const rotated = aa.angleDeg === 0 ? s : rotate(s, aa.angleDeg, {
 				at: [
 					0,
 					0,
 					0
 				],
 				axis: aa.axis
-			}), target);
+			});
+			return translate(rotated, target);
 		},
 		point: (p) => vecAdd(applyR(p), target),
 		vector: applyR
@@ -2679,12 +2688,13 @@ function matrixToAxisAngle$2(m) {
 		};
 	}
 	const denom = 2 * Math.sin(angle);
+	const axis = [
+		(m21 - m12) / denom,
+		(m02 - m20) / denom,
+		(m10 - m01) / denom
+	];
 	return {
-		axis: vecNormalize([
-			(m21 - m12) / denom,
-			(m02 - m20) / denom,
-			(m10 - m01) / denom
-		]),
+		axis: vecNormalize(axis),
 		angleDeg
 	};
 }
@@ -2706,14 +2716,15 @@ function buildBendPatch$1(axisOrigin, thetaDeg, r, width, thickness, direction)
 	});
 	let tube = outer;
 	if (r > 1e-9) {
-		const carved = cut(outer, cylinder(r, width, {
+		const inner = cylinder(r, width, {
 			at: axisOrigin,
 			axis: [
 				0,
 				1,
 				0
 			]
-		}));
+		});
+		const carved = cut(outer, inner);
 		if (!carved.ok) return carved;
 		tube = carved.value;
 	}
@@ -2731,39 +2742,47 @@ function buildBendPatch$1(axisOrigin, thetaDeg, r, width, thickness, direction)
 function buildWedge$1(axisOrigin, thetaDeg, outerR, width, direction) {
 	const span = outerR * 2 + 2;
 	const margin = 1;
-	const halfA = translate(box(span, width + 2 * margin, 2 * span), [
+	const blockA = box(span, width + 2, 2 * span);
+	const halfA = translate(blockA, [
 		axisOrigin[0],
 		axisOrigin[1] - margin,
 		axisOrigin[2] - span
 	]);
-	return intersect(halfA, rotate(halfA, (direction === "up" ? 1 : -1) * (180 - thetaDeg), {
+	const halfB = rotate(halfA, (direction === "up" ? 1 : -1) * (180 - thetaDeg), {
 		at: axisOrigin,
 		axis: [
 			0,
 			1,
 			0
 		]
-	}));
+	});
+	return intersect(halfA, halfB);
 }
 /** The flange flat, built lying along +X past the bend, then folded to angle θ. */
 function buildFlangeFlat(axisOrigin, thetaDeg, width, thickness, length, direction) {
-	if (direction === "up") return rotate(translate(box(length, width, thickness), [
-		axisOrigin[0],
-		axisOrigin[1],
-		0
-	]), -thetaDeg, {
-		at: axisOrigin,
-		axis: [
-			0,
-			1,
+	if (direction === "up") {
+		const flat = box(length, width, thickness);
+		const positioned = translate(flat, [
+			axisOrigin[0],
+			axisOrigin[1],
 			0
-		]
-	});
-	return rotate(translate(box(length, width, thickness), [
+		]);
+		return rotate(positioned, -thetaDeg, {
+			at: axisOrigin,
+			axis: [
+				0,
+				1,
+				0
+			]
+		});
+	}
+	const flat = box(length, width, thickness);
+	const positioned = translate(flat, [
 		axisOrigin[0],
 		axisOrigin[1],
 		-thickness
-	]), thetaDeg, {
+	]);
+	return rotate(positioned, thetaDeg, {
 		at: axisOrigin,
 		axis: [
 			0,
@@ -2772,8 +2791,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,
@@ -3047,7 +3065,8 @@ function cornerPoint$1(originA, dirA, originB, dirB) {
 	const e = vecDot(dirB, w0);
 	const denom = a * c - b * b;
 	if (Math.abs(denom) < 1e-9) return originA;
-	return vecAdd(originA, vecScale(dirA, (b * e - c * d) / denom));
+	const sc = (b * e - c * d) / denom;
+	return vecAdd(originA, vecScale(dirA, sc));
 }
 /** In-plane perpendicular to `axis` pointing toward the part center (into the parent). */
 function inwardRun(axisOrigin, axis, center) {
@@ -3080,11 +3099,12 @@ function parentEdgeLength(part, flange) {
 */
 function slotTool3d(at, along, inward, width, depth, thickness) {
 	const margin = thickness + 2;
-	const y0 = vecDot([
+	const rotPerp = [
 		-along[1],
 		along[0],
 		0
-	], inward) >= 0 ? 0 : -depth;
+	];
+	const y0 = vecDot(rotPerp, inward) >= 0 ? 0 : -depth;
 	let tool = box(width, depth, thickness + 2 * margin);
 	tool = translate(tool, [
 		-width / 2,
@@ -3109,7 +3129,8 @@ function squareTool(solid, cornerXY, side) {
 	const b = getBounds(solid);
 	const margin = 2;
 	const h = side / 2;
-	return translate(box(side, side, b.zMax - b.zMin + 2 * margin), [
+	const tool = box(side, side, b.zMax - b.zMin + 4);
+	return translate(tool, [
 		cornerXY[0] - h,
 		cornerXY[1] - h,
 		b.zMin - margin
@@ -3134,8 +3155,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;
@@ -3363,8 +3383,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
@@ -3506,8 +3525,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;
@@ -3602,18 +3620,20 @@ function embossForm(part, partSolid, spec, regionId, world, dev, ext) {
 	let solid;
 	if (spec.form === "emboss") {
 		const base = vecAdd(faceCentre, vecScale(world.n, -.05));
-		const fused = fuse(partSolid, cylinder(r, spec.height + HAIR, {
+		const bump = cylinder(r, spec.height + HAIR, {
 			at: base,
 			axis: world.n
-		}));
+		});
+		const fused = fuse(partSolid, bump);
 		if (!fused.ok) return fused;
 		solid = normalizeSolid(fused.value);
 	} else {
 		const top = vecAdd(faceCentre, vecScale(world.n, HAIR));
-		const cutResult = cut(partSolid, cylinder(r, spec.height + HAIR, {
+		const tool = cylinder(r, spec.height + HAIR, {
 			at: top,
 			axis: vecScale(world.n, -1)
-		}));
+		});
+		const cutResult = cut(partSolid, tool);
 		if (!cutResult.ok) return cutResult;
 		solid = normalizeSolid(cutResult.value);
 	}
@@ -3705,7 +3725,8 @@ function buildLouverFlap(f, x0, x1, y0, width, height, thickness, sign) {
 	const n = vecNormalize(f.n);
 	const v = vecNormalize(vecCross(n, u));
 	const hingeCorner = vecAdd(vecAdd(vecAdd(f.origin, vecScale(f.u, x0)), vecScale(f.v, y0)), vecScale(f.n, thickness));
-	return ok(mapLocalToWorld(flap, u, v, n, hingeCorner));
+	const placed = mapLocalToWorld(flap, u, v, n, hingeCorner);
+	return ok(placed);
 }
 /** Map a solid built in a canonical XYZ frame onto world axes (uT, vT, nT) at origin. */
 function mapLocalToWorld(s, uT, vT, nT, origin) {
@@ -3720,14 +3741,15 @@ function mapLocalToWorld(s, uT, vT, nT, origin) {
 		vT[2],
 		nT[2]
 	]);
-	return translate(aa.angleDeg === 0 ? s : rotate(s, aa.angleDeg, {
+	const rotated = aa.angleDeg === 0 ? s : rotate(s, aa.angleDeg, {
 		at: [
 			0,
 			0,
 			0
 		],
 		axis: aa.axis
-	}), origin);
+	});
+	return translate(rotated, origin);
 }
 /** Axis-angle of a 3×3 row-major rotation matrix (columns = target basis). */
 function matrixToAxisAngle$1(m) {
@@ -3785,12 +3807,13 @@ function matrixToAxisAngle$1(m) {
 		};
 	}
 	const denom = 2 * Math.sin(angle);
+	const axis = [
+		(m21 - m12) / denom,
+		(m02 - m20) / denom,
+		(m10 - m01) / denom
+	];
 	return {
-		axis: vecNormalize([
-			(m21 - m12) / denom,
-			(m02 - m20) / denom,
-			(m10 - m01) / denom
-		]),
+		axis: vecNormalize(axis),
 		angleDeg
 	};
 }
@@ -3803,8 +3826,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
@@ -3986,8 +4008,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
@@ -4645,8 +4666,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
@@ -4705,7 +4725,8 @@ function cornerPoint(originA, dirA, originB, dirB) {
 	const e = vecDot(dirB, w0);
 	const denom = a * c - b * b;
 	if (Math.abs(denom) < 1e-9) return originA;
-	return vecAdd(originA, vecScale(dirA, (b * e - c * d) / denom));
+	const sc = (b * e - c * d) / denom;
+	return vecAdd(originA, vecScale(dirA, sc));
 }
 /**
 * Normal of the miter plane between two flanges whose bend axes are `dirA`/`dirB`.
@@ -4737,13 +4758,14 @@ function halfSpaceTool(solid, plane) {
 		-size / 2,
 		0
 	]);
-	return ok(translate(alignZTo(tool, normal), plane.origin));
+	const aligned = alignZTo(tool, normal);
+	return ok(translate(aligned, plane.origin));
 }
 /** Rotate a shape so its local +Z axis points along `target` (about the origin). */
 function alignZTo(shape, target) {
 	const dot = vecDot(Z_AXIS, target);
-	if (dot > .999999999) return shape;
-	if (dot < -.999999999) return rotate(shape, 180, {
+	if (dot > 1 - 1e-9) return shape;
+	if (dot < -1 + 1e-9) return rotate(shape, 180, {
 		at: [
 			0,
 			0,
@@ -4756,7 +4778,8 @@ function alignZTo(shape, target) {
 		]
 	});
 	const axis = vecNormalize(vecCross(Z_AXIS, target));
-	return rotate(shape, Math.acos(Math.max(-1, Math.min(1, dot))) * 180 / Math.PI, {
+	const angleDeg = Math.acos(Math.max(-1, Math.min(1, dot))) * 180 / Math.PI;
+	return rotate(shape, angleDeg, {
 		at: [
 			0,
 			0,
@@ -4765,8 +4788,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
@@ -4887,7 +4909,8 @@ function authorContourFlange(part, spec) {
 function buildLineLeg(base, frame, span, thickness, legLength) {
 	const place = frameTransform(frame.run, frame.axis, frame.n, frame.origin);
 	const canon = box(legLength, span, thickness);
-	const fused = fuse(base, place.solid(canon));
+	const placed = place.solid(canon);
+	const fused = fuse(base, placed);
 	if (!fused.ok) return fused;
 	const nextOrigin = vecAdd(frame.origin, vecScale(frame.run, legLength));
 	return ok({
@@ -4919,7 +4942,8 @@ function buildArcBend(base, frame, span, thickness, seg) {
 	const place = frameTransform(frame.run, frame.axis, frame.n, frame.origin);
 	const patch = buildBendPatch(canonOrigin, seg.angleDeg, r, span, thickness, seg.direction);
 	if (!patch.ok) return patch;
-	const fused = fuse(base, place.solid(patch.value));
+	const placed = place.solid(patch.value);
+	const fused = fuse(base, placed);
 	if (!fused.ok) return fused;
 	const theta = seg.angleDeg * Math.PI / 180;
 	const a = -sign * theta;
@@ -4976,14 +5000,15 @@ function frameTransform(runT, axisT, nT, origin) {
 	const applyR = (v) => vecAdd(vecAdd(vecScale(runT, v[0]), vecScale(axisT, v[1])), vecScale(nT, v[2]));
 	return {
 		solid: (s) => {
-			return translate(aa.angleDeg === 0 ? s : rotate(s, aa.angleDeg, {
+			const rotated = aa.angleDeg === 0 ? s : rotate(s, aa.angleDeg, {
 				at: [
 					0,
 					0,
 					0
 				],
 				axis: aa.axis
-			}), origin);
+			});
+			return translate(rotated, origin);
 		},
 		point: (p) => vecAdd(applyR(p), origin),
 		vector: applyR
@@ -5044,10 +5069,10 @@ function baseEdge(part, side) {
 			length = uLen;
 			a = o;
 			b = uTop;
-			break;
 	}
 	const dir = vecNormalize(vecCross(n, out));
-	const start = vecDot(vecSub(b, a), dir) >= 0 ? a : b;
+	const toB = vecSub(b, a);
+	const start = vecDot(toB, dir) >= 0 ? a : b;
 	return {
 		dir,
 		out,
@@ -5094,10 +5119,10 @@ function regionEdge(frame, side) {
 			length = frame.uLen;
 			a = o;
 			b = uTop;
-			break;
 	}
 	const dir = vecNormalize(vecCross(frame.n, out));
-	const start = vecDot(vecSub(b, a), dir) >= 0 ? a : b;
+	const toB = vecSub(b, a);
+	const start = vecDot(toB, dir) >= 0 ? a : b;
 	return {
 		dir,
 		out,
@@ -5131,14 +5156,15 @@ function buildBendPatch(axisOrigin, thetaDeg, r, width, thickness, direction) {
 	});
 	let tube = outer;
 	if (r > 1e-9) {
-		const carved = cut(outer, cylinder(r, width, {
+		const inner = cylinder(r, width, {
 			at: axisOrigin,
 			axis: [
 				0,
 				1,
 				0
 			]
-		}));
+		});
+		const carved = cut(outer, inner);
 		if (!carved.ok) return carved;
 		tube = carved.value;
 	}
@@ -5149,19 +5175,21 @@ function buildBendPatch(axisOrigin, thetaDeg, r, width, thickness, direction) {
 function buildWedge(axisOrigin, thetaDeg, outerR, width, direction) {
 	const span = outerR * 2 + 2;
 	const margin = 1;
-	const halfA = translate(box(span, width + 2 * margin, 2 * span), [
+	const blockA = box(span, width + 2, 2 * span);
+	const halfA = translate(blockA, [
 		axisOrigin[0],
 		axisOrigin[1] - margin,
 		axisOrigin[2] - span
 	]);
-	return intersect(halfA, rotate(halfA, (direction === "up" ? 1 : -1) * (180 - thetaDeg), {
+	const halfB = rotate(halfA, (direction === "up" ? 1 : -1) * (180 - thetaDeg), {
 		at: axisOrigin,
 		axis: [
 			0,
 			1,
 			0
 		]
-	}));
+	});
+	return intersect(halfA, halfB);
 }
 /** Axis-angle of a 3×3 rotation matrix (row-major). */
 function matrixToAxisAngle(m) {
@@ -5219,17 +5247,17 @@ function matrixToAxisAngle(m) {
 		};
 	}
 	const denom = 2 * Math.sin(angle);
+	const axis = [
+		(m21 - m12) / denom,
+		(m02 - m20) / denom,
+		(m10 - m01) / denom
+	];
 	return {
-		axis: vecNormalize([
-			(m21 - m12) / denom,
-			(m02 - m20) / denom,
-			(m10 - m01) / denom
-		]),
+		axis: vecNormalize(axis),
 		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
@@ -5349,10 +5377,13 @@ function buildQuadPlate(id, a0, a1, b1, b0, thickness) {
 }
 /** A flat triangle (p0,p1,p2) extruded by `thickness` along its own normal. */
 function buildTrianglePlate(id, p0, p1, p2, thickness) {
+	const e0 = line(p0, p1);
+	const e1 = line(p1, p2);
+	const e2 = line(p2, p0);
 	const loopResult = wireLoop([
-		line(p0, p1),
-		line(p1, p2),
-		line(p2, p0)
+		e0,
+		e1,
+		e2
 	]);
 	if (!loopResult.ok) return loopResult;
 	const loop = loopResult.value;
@@ -5538,8 +5569,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
@@ -5732,8 +5762,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.
@@ -5892,8 +5921,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";
@@ -6102,8 +6130,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
@@ -6346,8 +6373,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
@@ -6813,8 +6839,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
@@ -6857,8 +6882,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);
@@ -7020,8 +7044,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;
@@ -7291,5 +7314,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 };