brepjs 18.67.0 → 18.68.0

package/dist/brepjs.cjs

@@ -182,7 +182,7 @@ function validateInputs(mesh, queries) {
 	return null;
 }
 /** Resolve a registered voxel engine, mapping an unregistered id to an error. */
-function resolveEngine(id) {
+function resolveEngine$1(id) {
 	try {
 		return require_errors.ok(getVoxel(id));
 	} catch (cause) {
@@ -199,7 +199,7 @@ function resolveEngine(id) {
 function windingNumbers(mesh, queries, id) {
 	const invalid = validateInputs(mesh, queries);
 	if (invalid) return require_errors.err(invalid);
-	const engine = resolveEngine(id);
+	const engine = resolveEngine$1(id);
 	if (require_errors.isErr(engine)) return engine;
 	return require_errors.ok(engine.value.winding_numbers(mesh.vertices, mesh.triangles, queries));
 }
@@ -211,15 +211,15 @@ function windingNumbers(mesh, queries, id) {
 function pointsInside(mesh, queries, id) {
 	const invalid = validateInputs(mesh, queries);
 	if (invalid) return require_errors.err(invalid);
-	const engine = resolveEngine(id);
+	const engine = resolveEngine$1(id);
 	if (require_errors.isErr(engine)) return engine;
 	const flags = engine.value.points_inside(mesh.vertices, mesh.triangles, queries);
 	return require_errors.ok(Array.from(flags, (flag) => flag === 1));
 }
 //#endregion
 //#region src/voxel/repairFns.ts
-var DEFAULT_RESOLUTION$3 = 48;
-var DEFAULT_PADDING$3 = 2;
+var DEFAULT_RESOLUTION$4 = 48;
+var DEFAULT_PADDING$4 = 2;
 /**
 * Repair a (possibly non-watertight) triangle-soup mesh into a closed surface:
 * voxelize an FWN-signed SDF, then Surface-Nets contour it back to triangles.
@@ -232,11 +232,11 @@ function repairMesh(mesh, opts, id) {
 	const invalid = validateMesh(mesh);
 	if (invalid) return require_errors.err(invalid);
 	if (mesh.vertices.length === 0 || mesh.triangles.length === 0) return require_errors.err(require_errors.validationError("VOXEL_EMPTY_MESH", "repairMesh requires a non-empty triangle mesh."));
-	const resolution = opts?.resolution ?? DEFAULT_RESOLUTION$3;
-	const padding = opts?.padding ?? DEFAULT_PADDING$3;
+	const resolution = opts?.resolution ?? DEFAULT_RESOLUTION$4;
+	const padding = opts?.padding ?? DEFAULT_PADDING$4;
 	if (!Number.isInteger(resolution) || resolution < 1) return require_errors.err(require_errors.validationError("VOXEL_INVALID_RESOLUTION", "resolution must be an integer >= 1."));
 	if (!Number.isInteger(padding) || padding < 1) return require_errors.err(require_errors.validationError("VOXEL_INVALID_PADDING", "padding must be an integer >= 1."));
-	const engine = resolveEngine(id);
+	const engine = resolveEngine$1(id);
 	if (require_errors.isErr(engine)) return engine;
 	try {
 		try {
@@ -286,16 +286,16 @@ function shapeToMeshInput(shape, deflection = DEFAULT_DEFLECTION) {
 }
 //#endregion
 //#region src/voxel/meshOpsFns.ts
-var DEFAULT_RESOLUTION$2 = 48;
-var DEFAULT_PADDING$2 = 2;
+var DEFAULT_RESOLUTION$3 = 48;
+var DEFAULT_PADDING$3 = 2;
 var BOOLEAN_OP_CODES$1 = {
 	union: 0,
 	intersection: 1,
 	difference: 2
 };
-function resolveGridParams$1(opts) {
-	const resolution = opts?.resolution ?? DEFAULT_RESOLUTION$2;
-	const padding = opts?.padding ?? DEFAULT_PADDING$2;
+function resolveGridParams$2(opts) {
+	const resolution = opts?.resolution ?? DEFAULT_RESOLUTION$3;
+	const padding = opts?.padding ?? DEFAULT_PADDING$3;
 	if (!Number.isInteger(resolution) || resolution < 1) return require_errors.err(require_errors.validationError("VOXEL_INVALID_RESOLUTION", "resolution must be an integer >= 1."));
 	if (!Number.isInteger(padding) || padding < 1) return require_errors.err(require_errors.validationError("VOXEL_INVALID_PADDING", "padding must be an integer >= 1."));
 	return require_errors.ok({
@@ -330,9 +330,9 @@ function offsetMesh(mesh, distance, opts, id) {
 	if (invalid) return require_errors.err(invalid);
 	if (mesh.vertices.length === 0 || mesh.triangles.length === 0) return require_errors.err(require_errors.validationError("VOXEL_EMPTY_MESH", "offsetMesh requires a non-empty triangle mesh."));
 	if (!Number.isFinite(distance)) return require_errors.err(require_errors.validationError("VOXEL_INVALID_DISTANCE", "distance must be a finite number."));
-	const params = resolveGridParams$1(opts);
+	const params = resolveGridParams$2(opts);
 	if (require_errors.isErr(params)) return params;
-	const engine = resolveEngine(id);
+	const engine = resolveEngine$1(id);
 	if (require_errors.isErr(engine)) return engine;
 	try {
 		try {
@@ -359,9 +359,9 @@ function shellMesh(mesh, thickness, opts, id) {
 	if (invalid) return require_errors.err(invalid);
 	if (mesh.vertices.length === 0 || mesh.triangles.length === 0) return require_errors.err(require_errors.validationError("VOXEL_EMPTY_MESH", "shellMesh requires a non-empty triangle mesh."));
 	if (!Number.isFinite(thickness) || thickness <= 0) return require_errors.err(require_errors.validationError("VOXEL_INVALID_THICKNESS", "thickness must be a finite number > 0."));
-	const params = resolveGridParams$1(opts);
+	const params = resolveGridParams$2(opts);
 	if (require_errors.isErr(params)) return params;
-	const engine = resolveEngine(id);
+	const engine = resolveEngine$1(id);
 	if (require_errors.isErr(engine)) return engine;
 	try {
 		try {
@@ -391,9 +391,9 @@ function voxelBoolean(a, b, op, opts, id) {
 	if (a.vertices.length === 0 || a.triangles.length === 0) return require_errors.err(require_errors.validationError("VOXEL_EMPTY_MESH", "voxelBoolean requires a non-empty mesh for operand A."));
 	if (b.vertices.length === 0 || b.triangles.length === 0) return require_errors.err(require_errors.validationError("VOXEL_EMPTY_MESH", "voxelBoolean requires a non-empty mesh for operand B."));
 	const opCode = BOOLEAN_OP_CODES$1[op];
-	const params = resolveGridParams$1(opts);
+	const params = resolveGridParams$2(opts);
 	if (require_errors.isErr(params)) return params;
-	const engine = resolveEngine(id);
+	const engine = resolveEngine$1(id);
 	if (require_errors.isErr(engine)) return engine;
 	try {
 		try {
@@ -442,16 +442,16 @@ function voxelBooleanShapes(a, b, op, opts, id) {
 }
 //#endregion
 //#region src/voxel/fieldFns.ts
-var DEFAULT_RESOLUTION$1 = 48;
-var DEFAULT_PADDING$1 = 2;
+var DEFAULT_RESOLUTION$2 = 48;
+var DEFAULT_PADDING$2 = 2;
 var BOOLEAN_OP_CODES = {
 	union: 0,
 	intersection: 1,
 	difference: 2
 };
-function resolveGridParams(opts) {
-	const resolution = opts?.resolution ?? DEFAULT_RESOLUTION$1;
-	const padding = opts?.padding ?? DEFAULT_PADDING$1;
+function resolveGridParams$1(opts) {
+	const resolution = opts?.resolution ?? DEFAULT_RESOLUTION$2;
+	const padding = opts?.padding ?? DEFAULT_PADDING$2;
 	if (!Number.isInteger(resolution) || resolution < 1) return require_errors.err(require_errors.validationError("VOXEL_INVALID_RESOLUTION", "resolution must be an integer >= 1."));
 	if (!Number.isInteger(padding) || padding < 1) return require_errors.err(require_errors.validationError("VOXEL_INVALID_PADDING", "padding must be an integer >= 1."));
 	return require_errors.ok({
@@ -501,6 +501,9 @@ function fieldDeletable(raw) {
 		}
 	};
 }
+function makeFieldHandle(raw) {
+	return makeHandle(raw);
+}
 function makeHandle(raw) {
 	const inner = require_shapeTypes.createKernelHandle(fieldDeletable(raw));
 	const handle = {
@@ -557,9 +560,9 @@ function voxelField(mesh, opts, id) {
 	const invalid = validateMesh(mesh);
 	if (invalid) return require_errors.err(invalid);
 	if (mesh.vertices.length === 0 || mesh.triangles.length === 0) return require_errors.err(require_errors.validationError("VOXEL_EMPTY_MESH", "voxelField requires a non-empty triangle mesh."));
-	const params = resolveGridParams(opts);
+	const params = resolveGridParams$1(opts);
 	if (require_errors.isErr(params)) return params;
-	const engine = resolveEngine(id);
+	const engine = resolveEngine$1(id);
 	if (require_errors.isErr(engine)) return engine;
 	try {
 		return require_errors.ok(makeHandle(new engine.value.VoxelField(mesh.vertices, mesh.triangles, params.value.resolution, params.value.padding)));
@@ -586,9 +589,9 @@ function voxelBooleanField(a, b, op, opts, id) {
 	if (invalidB) return require_errors.err(invalidB);
 	if (a.vertices.length === 0 || a.triangles.length === 0) return require_errors.err(require_errors.validationError("VOXEL_EMPTY_MESH", "voxelBooleanField requires a non-empty mesh for operand A."));
 	if (b.vertices.length === 0 || b.triangles.length === 0) return require_errors.err(require_errors.validationError("VOXEL_EMPTY_MESH", "voxelBooleanField requires a non-empty mesh for operand B."));
-	const params = resolveGridParams(opts);
+	const params = resolveGridParams$1(opts);
 	if (require_errors.isErr(params)) return params;
-	const engine = resolveEngine(id);
+	const engine = resolveEngine$1(id);
 	if (require_errors.isErr(engine)) return engine;
 	try {
 		return require_errors.ok(makeHandle(engine.value.VoxelField.boolean_of(a.vertices, a.triangles, b.vertices, b.triangles, BOOLEAN_OP_CODES[op], params.value.resolution, params.value.padding)));
@@ -698,6 +701,163 @@ function voxelBooleanFieldShapes(a, b, op, opts, id) {
 	return voxelBooleanField(meshA.value, meshB.value, op, opts, id);
 }
 //#endregion
+//#region src/implicit/sdfFns.ts
+var DEFAULT_RESOLUTION$1 = 48;
+var DEFAULT_PADDING$1 = 2;
+function resolveEngine(id) {
+	try {
+		return require_errors.ok(getVoxel(id));
+	} catch (cause) {
+		return require_errors.err(require_errors.moduleInitError("VOXEL_NOT_INITIALIZED", cause instanceof Error ? cause.message : "voxel engine not initialized", cause));
+	}
+}
+function resolveGridParams(opts) {
+	const resolution = opts?.resolution ?? DEFAULT_RESOLUTION$1;
+	const padding = opts?.padding ?? DEFAULT_PADDING$1;
+	if (!Number.isInteger(resolution) || resolution < 1) return require_errors.err(require_errors.validationError("SDF_INVALID_RESOLUTION", "resolution must be an integer >= 1."));
+	if (!Number.isInteger(padding) || padding < 1) return require_errors.err(require_errors.validationError("SDF_INVALID_PADDING", "padding must be an integer >= 1."));
+	return require_errors.ok({
+		resolution,
+		padding
+	});
+}
+function sdfDeletable(raw) {
+	return {
+		raw,
+		delete() {
+			raw.free();
+		}
+	};
+}
+function makeSdfHandle(raw) {
+	const inner = require_shapeTypes.createKernelHandle(sdfDeletable(raw));
+	return {
+		get value() {
+			return inner.value.raw;
+		},
+		get disposed() {
+			return inner.disposed;
+		},
+		[Symbol.dispose]() {
+			inner[Symbol.dispose]();
+		},
+		union(other) {
+			return makeSdfHandle(this.value.union(other.value));
+		},
+		intersection(other) {
+			return makeSdfHandle(this.value.intersection(other.value));
+		},
+		difference(other) {
+			return makeSdfHandle(this.value.difference(other.value));
+		},
+		smoothUnion(other, k) {
+			return makeSdfHandle(this.value.smooth_union(other.value, k));
+		},
+		smoothIntersection(other, k) {
+			return makeSdfHandle(this.value.smooth_intersection(other.value, k));
+		},
+		smoothDifference(other, k) {
+			return makeSdfHandle(this.value.smooth_difference(other.value, k));
+		},
+		offset(distance) {
+			return makeSdfHandle(this.value.offset(distance));
+		},
+		round(radius) {
+			return makeSdfHandle(this.value.round(radius));
+		},
+		shell(thickness) {
+			return makeSdfHandle(this.value.shell(thickness));
+		},
+		onion(thickness) {
+			return makeSdfHandle(this.value.onion(thickness));
+		},
+		translate(x, y, z) {
+			return makeSdfHandle(this.value.translate(x, y, z));
+		},
+		rotate(ax, ay, az, angle) {
+			return makeSdfHandle(this.value.rotate(ax, ay, az, angle));
+		},
+		scale(s) {
+			return makeSdfHandle(this.value.scale(s));
+		},
+		rasterize(opts) {
+			return rasterizeField(this.value, opts);
+		},
+		rasterizeIn(bounds, opts) {
+			return rasterizeFieldIn(this.value, bounds, opts);
+		}
+	};
+}
+function rasterizeField(sdf, opts) {
+	const params = resolveGridParams(opts);
+	if (require_errors.isErr(params)) return params;
+	try {
+		return require_errors.ok(makeFieldHandle(sdf.rasterize(params.value.resolution, params.value.padding)));
+	} catch (cause) {
+		return rasterizeError(cause);
+	}
+}
+function rasterizeFieldIn(sdf, bounds, opts) {
+	const params = resolveGridParams(opts);
+	if (require_errors.isErr(params)) return params;
+	const invalid = validateBounds(bounds);
+	if (invalid) return require_errors.err(invalid);
+	try {
+		return require_errors.ok(makeFieldHandle(sdf.rasterize_in(bounds.min[0], bounds.min[1], bounds.min[2], bounds.max[0], bounds.max[1], bounds.max[2], params.value.resolution, params.value.padding)));
+	} catch (cause) {
+		return rasterizeError(cause);
+	}
+}
+function rasterizeError(cause) {
+	return require_errors.err(require_errors.computationError("SDF_RASTERIZE_FAILED", cause instanceof Error ? cause.message : "SDF rasterization failed (grid too large?).", cause));
+}
+function validateBounds(bounds) {
+	if (![...bounds.min, ...bounds.max].every((v) => Number.isFinite(v))) return require_errors.validationError("SDF_INVALID_BOUNDS", "bounds must be finite numbers.");
+	for (let axis = 0; axis < 3; axis++) if (bounds.max[axis] <= bounds.min[axis]) return require_errors.validationError("SDF_INVALID_BOUNDS", "bounds max must exceed min on every axis.");
+	return null;
+}
+function build(make, id) {
+	const engine = resolveEngine(id);
+	if (require_errors.isErr(engine)) return engine;
+	try {
+		return require_errors.ok(makeSdfHandle(make(engine.value)));
+	} catch (cause) {
+		return require_errors.err(require_errors.computationError("SDF_BUILD_FAILED", cause instanceof Error ? cause.message : "SDF primitive construction failed.", cause));
+	}
+}
+/** A sphere of radius `r`, centered at the origin. */
+function sphere(r, id) {
+	return build((e) => e.Sdf.sphere(r), id);
+}
+/** An axis-aligned box of half-extents `(hx, hy, hz)`, centered at the origin. */
+function box$1(hx, hy, hz, id) {
+	return build((e) => e.Sdf.box_(hx, hy, hz), id);
+}
+/** A box with rounded edges of radius `r`. */
+function roundedBox(hx, hy, hz, r, id) {
+	return build((e) => e.Sdf.rounded_box(hx, hy, hz, r), id);
+}
+/** A capped cylinder, axis +Z, radius `r`, total height `h`, centered at origin. */
+function cylinder$1(r, h, id) {
+	return build((e) => e.Sdf.cylinder(r, h), id);
+}
+/** A capped cone centered at the origin: base radius `r` at z = −h/2 tapering to an apex at z = +h/2. */
+function cone$1(r, h, id) {
+	return build((e) => e.Sdf.cone(r, h), id);
+}
+/** A capsule: a line segment `a`→`b` of radius `r`. */
+function capsule(a, b, r, id) {
+	return build((e) => e.Sdf.capsule(a[0], a[1], a[2], b[0], b[1], b[2], r), id);
+}
+/** A torus in the XY plane (axis +Z): a `minor`-radius tube on a `major` circle. */
+function torus(major, minor, id) {
+	return build((e) => e.Sdf.torus(major, minor), id);
+}
+/** A half-space: the plane through `h·n` with outward normal `n` (normalized). */
+function plane(n, h, id) {
+	return build((e) => e.Sdf.plane(n[0], n[1], n[2], h), id);
+}
+//#endregion
 //#region src/lattice/latticeFns.ts
 var LATTICE_TAGS = {
 	gyroid: 0,
@@ -751,7 +911,7 @@ function latticeInfill(mesh, opts, id) {
 	if (mesh.vertices.length === 0 || mesh.triangles.length === 0) return require_errors.err(require_errors.validationError("LATTICE_EMPTY_MESH", "latticeInfill requires a non-empty triangle mesh."));
 	const resolved = validateOptions(opts);
 	if (require_errors.isErr(resolved)) return resolved;
-	const engine = resolveEngine(id);
+	const engine = resolveEngine$1(id);
 	if (require_errors.isErr(engine)) return engine;
 	const { tag, period, thickness, resolution, padding } = resolved.value;
 	try {
@@ -792,7 +952,7 @@ function tpmsLattice(bounds, opts, id) {
 	}
 	const resolved = validateOptions(opts);
 	if (require_errors.isErr(resolved)) return resolved;
-	const engine = resolveEngine(id);
+	const engine = resolveEngine$1(id);
 	if (require_errors.isErr(engine)) return engine;
 	const { tag, period, thickness, resolution, padding } = resolved.value;
 	const [minX, minY, minZ] = bounds.min;
@@ -4461,7 +4621,7 @@ function depsOf(...sources) {
 	}
 	return acc.size === 0 ? EMPTY_DEPS : acc;
 }
-function box$1(x, y, z) {
+function box$2(x, y, z) {
 	const xe = asScalarExpr(x);
 	const ye = asScalarExpr(y);
 	const ze = asScalarExpr(z);
@@ -4476,7 +4636,7 @@ function box$1(x, y, z) {
 		freeParams: depsOf(xe, ye, ze)
 	};
 }
-function sphere$1(radius) {
+function sphere$2(radius) {
 	const re = asScalarExpr(radius);
 	return {
 		kind: "Sphere",
@@ -4485,7 +4645,7 @@ function sphere$1(radius) {
 		freeParams: re.freeParams
 	};
 }
-function cylinder$1(radius, height) {
+function cylinder$2(radius, height) {
 	const re = asScalarExpr(radius);
 	const he = asScalarExpr(height);
 	return {
@@ -4496,7 +4656,7 @@ function cylinder$1(radius, height) {
 		freeParams: depsOf(re, he)
 	};
 }
-function cone$1(radius1, radius2, height) {
+function cone$2(radius1, radius2, height) {
 	const r1 = asScalarExpr(radius1);
 	const r2 = asScalarExpr(radius2);
 	const he = asScalarExpr(height);
@@ -4511,7 +4671,7 @@ function cone$1(radius1, radius2, height) {
 		freeParams: depsOf(r1, r2, he)
 	};
 }
-function torus$1(majorRadius, minorRadius) {
+function torus$2(majorRadius, minorRadius) {
 	const ma = asScalarExpr(majorRadius);
 	const mi = asScalarExpr(minorRadius);
 	return {
@@ -5416,7 +5576,7 @@ function readSingleExpr(j, key, build) {
 function readPrimitive(kind, j) {
 	switch (kind) {
 		case "Box": return readBox(j);
-		case "Sphere": return readSingleExpr(j, "radius", sphere$1);
+		case "Sphere": return readSingleExpr(j, "radius", sphere$2);
 		case "Cylinder": return readCylinder(j);
 		case "Cone": return readCone(j);
 		case "Torus": return readTorus(j);
@@ -5435,14 +5595,14 @@ function readBox(j) {
 	if (!y.ok) return y;
 	const z = readExpr(j["z"]);
 	if (!z.ok) return z;
-	return require_errors.ok(box$1(x.value, y.value, z.value));
+	return require_errors.ok(box$2(x.value, y.value, z.value));
 }
 function readCylinder(j) {
 	const r = readExpr(j["radius"]);
 	if (!r.ok) return r;
 	const h = readExpr(j["height"]);
 	if (!h.ok) return h;
-	return require_errors.ok(cylinder$1(r.value, h.value));
+	return require_errors.ok(cylinder$2(r.value, h.value));
 }
 function readCone(j) {
 	const r1 = readExpr(j["radius1"]);
@@ -5451,14 +5611,14 @@ function readCone(j) {
 	if (!r2.ok) return r2;
 	const h = readExpr(j["height"]);
 	if (!h.ok) return h;
-	return require_errors.ok(cone$1(r1.value, r2.value, h.value));
+	return require_errors.ok(cone$2(r1.value, r2.value, h.value));
 }
 function readTorus(j) {
 	const ma = readExpr(j["majorRadius"]);
 	if (!ma.ok) return ma;
 	const mi = readExpr(j["minorRadius"]);
 	if (!mi.ok) return mi;
-	return require_errors.ok(torus$1(ma.value, mi.value));
+	return require_errors.ok(torus$2(ma.value, mi.value));
 }
 function readPolygon(j) {
 	const pts = j["points"];
@@ -5634,11 +5794,11 @@ function foldBuildVec(dim, comps) {
 }
 function optimizeNode(n) {
 	switch (n.kind) {
-		case "Box": return box$1(foldExpr(n.x), foldExpr(n.y), foldExpr(n.z));
-		case "Sphere": return sphere$1(foldExpr(n.radius));
-		case "Cylinder": return cylinder$1(foldExpr(n.radius), foldExpr(n.height));
-		case "Cone": return cone$1(foldExpr(n.radius1), foldExpr(n.radius2), foldExpr(n.height));
-		case "Torus": return torus$1(foldExpr(n.majorRadius), foldExpr(n.minorRadius));
+		case "Box": return box$2(foldExpr(n.x), foldExpr(n.y), foldExpr(n.z));
+		case "Sphere": return sphere$2(foldExpr(n.radius));
+		case "Cylinder": return cylinder$2(foldExpr(n.radius), foldExpr(n.height));
+		case "Cone": return cone$2(foldExpr(n.radius1), foldExpr(n.radius2), foldExpr(n.height));
+		case "Torus": return torus$2(foldExpr(n.majorRadius), foldExpr(n.minorRadius));
 		case "Polygon": return polygon$1(n.points.map(foldExpr));
 		case "Circle": return circle$1(foldExpr(n.radius));
 		case "Line": return line$1(foldExpr(n.from), foldExpr(n.to));
@@ -5796,15 +5956,15 @@ var csg_exports = /* @__PURE__ */ require_textBlueprints.__exportAll({
 	asVec2Expr: () => asVec2Expr,
 	asVec3Expr: () => asVec3Expr,
 	binOp: () => binOp,
-	box: () => box$1,
+	box: () => box$2,
 	buildVec: () => buildVec,
 	circle: () => circle$1,
 	component: () => component,
 	compound: () => compound$1,
-	cone: () => cone$1,
+	cone: () => cone$2,
 	cut: () => cut$1,
 	cutAll: () => cutAll$1,
-	cylinder: () => cylinder$1,
+	cylinder: () => cylinder$2,
 	emptyFace: () => emptyFace,
 	emptySolid: () => emptySolid,
 	emptyWire: () => emptyWire,
@@ -5826,9 +5986,9 @@ var csg_exports = /* @__PURE__ */ require_textBlueprints.__exportAll({
 	replaceNode: () => replaceNode,
 	rotate: () => rotate$1,
 	scale: () => scale$1,
-	sphere: () => sphere$1,
+	sphere: () => sphere$2,
 	toJSON: () => toJSON,
-	torus: () => torus$1,
+	torus: () => torus$2,
 	translate: () => translate$1,
 	unaryOp: () => unaryOp,
 	vec2Lit: () => vec2Lit,
@@ -6308,6 +6468,14 @@ exports.roundedRectangleBlueprint = require_boolean2D.roundedRectangleBlueprint;
 exports.scale = scale;
 exports.scale2D = require_blueprintFns.scale2D;
 exports.scaleDrawing = require_drawFns.scaleDrawing;
+exports.sdfBox = box$1;
+exports.sdfCapsule = capsule;
+exports.sdfCone = cone$1;
+exports.sdfCylinder = cylinder$1;
+exports.sdfPlane = plane;
+exports.sdfRoundedBox = roundedBox;
+exports.sdfSphere = sphere;
+exports.sdfTorus = torus;
 exports.section = section;
 exports.sectionToFace = sectionToFace;
 exports.serializeHistory = require_historyFns.serializeHistory;

package/dist/brepjs.js

@@ -14,7 +14,7 @@ import { a as meshEdges$1, c as createMeshCache, i as mesh$1, n as exportSTEP, o
 import { n as getAtOrThrow, r as lastOrThrow, t as firstOrThrow } from "./arrayAccess-DrUGPADn.js";
 import { _ as makeThreePointArc, d as makeCircle, h as makeLine, l as makeBSplineInterpolation, n as fill, r as makeFace, s as assembleWire } from "./surfaceBuilders-jx81G_YJ.js";
 import { a as fuseAll, c as sectionToFace$1, i as fuse$2, l as slice$1, n as cut$2, o as intersect$2, r as cutAll, s as section$1, t as booleanPipeline, u as split$1 } from "./booleanFns-CaKngiuu.js";
-import { $ as fuseAllBisect, A as fixShape, B as offset$1, C as threePointArc, D as wireLoop, E as wire, F as isValid$1, G as chamferWithEvolution, H as thicken$1, I as solidFromShell, J as fuseWithEvolution, K as cutWithEvolution, L as chamfer$1, M as healFace, N as healSolid, O as autoHeal, P as healWire, Q as cutAllBisect, R as draft$1, S as tangentArc, T as vertex, U as variableFillet, V as shell$1, W as positionOnCurve, X as shellWithEvolution, Y as intersectWithEvolution, Z as checkBoolean, _ as polygon, a as circle, at as sharedEdges, b as sphere, c as cylinder, ct as chamferDistAngle, d as ellipsoid, dt as toLODGeometryData, et as getNurbsCurveData, f as face, ft as toLineGeometryData, g as offsetFace, h as line, i as bsplineApprox, it as facesOfEdge, j as heal$1, k as fixSelfIntersection, l as ellipse, lt as toBufferGeometryData, m as helix, n as bezier, nt as adjacentFaces, o as compound, ot as verticesOfEdge, p as filledFace, q as filletWithEvolution, r as box, rt as edgesOfFace, s as cone, st as wiresOfFace, t as addHoles, tt as getNurbsSurfaceData, u as ellipseArc, ut as toGroupedBufferGeometryData, v as sewShells, w as torus, x as subFace, y as solid, z as fillet$1 } from "./primitiveFns-ecKWNC5k.js";
+import { $ as fuseAllBisect, A as fixShape, B as offset$1, C as threePointArc, D as wireLoop, E as wire, F as isValid$1, G as chamferWithEvolution, H as thicken$1, I as solidFromShell, J as fuseWithEvolution, K as cutWithEvolution, L as chamfer$1, M as healFace, N as healSolid, O as autoHeal, P as healWire, Q as cutAllBisect, R as draft$1, S as tangentArc, T as vertex, U as variableFillet, V as shell$1, W as positionOnCurve, X as shellWithEvolution, Y as intersectWithEvolution, Z as checkBoolean, _ as polygon, a as circle, at as sharedEdges, b as sphere$1, c as cylinder, ct as chamferDistAngle, d as ellipsoid, dt as toLODGeometryData, et as getNurbsCurveData, f as face, ft as toLineGeometryData, g as offsetFace, h as line, i as bsplineApprox, it as facesOfEdge, j as heal$1, k as fixSelfIntersection, l as ellipse, lt as toBufferGeometryData, m as helix, n as bezier, nt as adjacentFaces, o as compound, ot as verticesOfEdge, p as filledFace, q as filletWithEvolution, r as box, rt as edgesOfFace, s as cone, st as wiresOfFace, t as addHoles, tt as getNurbsSurfaceData, u as ellipseArc, ut as toGroupedBufferGeometryData, v as sewShells, w as torus$1, x as subFace, y as solid, z as fillet$1 } from "./primitiveFns-ecKWNC5k.js";
 import { C as walkAssembly, D as exportAssemblySTEP, E as linearPattern, O as createAssembly, S as updateNode, _ as collectShapes, a as findStep, b as findNode, c as registerOperation, d as replayHistory, f as serializeHistory, g as addChild, h as undoLast, i as deserializeHistory, l as registerShape, m as stepsFrom, n as createHistory, o as getShape, p as stepCount, r as createRegistry, s as modifyStep, t as addStep, u as replayFrom, v as countNodes, w as circularPattern, x as removeChild, y as createAssemblyNode } from "./historyFns-CZ9oNL7j.js";
 import { n as BaseSketcher2d, r as organiseBlueprints, t as BlueprintSketcher } from "./blueprintSketcher-0DeTMXwj.js";
 import { a as createTypedFinder, i as wireFinder, n as edgeFinder, r as faceFinder, t as getSingleFace } from "./helpers-BX-0e71G.js";
@@ -193,7 +193,7 @@ function validateInputs(mesh, queries) {
 	return null;
 }
 /** Resolve a registered voxel engine, mapping an unregistered id to an error. */
-function resolveEngine(id) {
+function resolveEngine$1(id) {
 	try {
 		return ok(getVoxel(id));
 	} catch (cause) {
@@ -210,7 +210,7 @@ function resolveEngine(id) {
 function windingNumbers(mesh, queries, id) {
 	const invalid = validateInputs(mesh, queries);
 	if (invalid) return err(invalid);
-	const engine = resolveEngine(id);
+	const engine = resolveEngine$1(id);
 	if (isErr(engine)) return engine;
 	return ok(engine.value.winding_numbers(mesh.vertices, mesh.triangles, queries));
 }
@@ -222,15 +222,15 @@ function windingNumbers(mesh, queries, id) {
 function pointsInside(mesh, queries, id) {
 	const invalid = validateInputs(mesh, queries);
 	if (invalid) return err(invalid);
-	const engine = resolveEngine(id);
+	const engine = resolveEngine$1(id);
 	if (isErr(engine)) return engine;
 	const flags = engine.value.points_inside(mesh.vertices, mesh.triangles, queries);
 	return ok(Array.from(flags, (flag) => flag === 1));
 }
 //#endregion
 //#region src/voxel/repairFns.ts
-var DEFAULT_RESOLUTION$3 = 48;
-var DEFAULT_PADDING$3 = 2;
+var DEFAULT_RESOLUTION$4 = 48;
+var DEFAULT_PADDING$4 = 2;
 /**
 * Repair a (possibly non-watertight) triangle-soup mesh into a closed surface:
 * voxelize an FWN-signed SDF, then Surface-Nets contour it back to triangles.
@@ -243,11 +243,11 @@ function repairMesh(mesh, opts, id) {
 	const invalid = validateMesh(mesh);
 	if (invalid) return err(invalid);
 	if (mesh.vertices.length === 0 || mesh.triangles.length === 0) return err(validationError("VOXEL_EMPTY_MESH", "repairMesh requires a non-empty triangle mesh."));
-	const resolution = opts?.resolution ?? DEFAULT_RESOLUTION$3;
-	const padding = opts?.padding ?? DEFAULT_PADDING$3;
+	const resolution = opts?.resolution ?? DEFAULT_RESOLUTION$4;
+	const padding = opts?.padding ?? DEFAULT_PADDING$4;
 	if (!Number.isInteger(resolution) || resolution < 1) return err(validationError("VOXEL_INVALID_RESOLUTION", "resolution must be an integer >= 1."));
 	if (!Number.isInteger(padding) || padding < 1) return err(validationError("VOXEL_INVALID_PADDING", "padding must be an integer >= 1."));
-	const engine = resolveEngine(id);
+	const engine = resolveEngine$1(id);
 	if (isErr(engine)) return engine;
 	try {
 		try {
@@ -297,16 +297,16 @@ function shapeToMeshInput(shape, deflection = DEFAULT_DEFLECTION) {
 }
 //#endregion
 //#region src/voxel/meshOpsFns.ts
-var DEFAULT_RESOLUTION$2 = 48;
-var DEFAULT_PADDING$2 = 2;
+var DEFAULT_RESOLUTION$3 = 48;
+var DEFAULT_PADDING$3 = 2;
 var BOOLEAN_OP_CODES$1 = {
 	union: 0,
 	intersection: 1,
 	difference: 2
 };
-function resolveGridParams$1(opts) {
-	const resolution = opts?.resolution ?? DEFAULT_RESOLUTION$2;
-	const padding = opts?.padding ?? DEFAULT_PADDING$2;
+function resolveGridParams$2(opts) {
+	const resolution = opts?.resolution ?? DEFAULT_RESOLUTION$3;
+	const padding = opts?.padding ?? DEFAULT_PADDING$3;
 	if (!Number.isInteger(resolution) || resolution < 1) return err(validationError("VOXEL_INVALID_RESOLUTION", "resolution must be an integer >= 1."));
 	if (!Number.isInteger(padding) || padding < 1) return err(validationError("VOXEL_INVALID_PADDING", "padding must be an integer >= 1."));
 	return ok({
@@ -341,9 +341,9 @@ function offsetMesh(mesh, distance, opts, id) {
 	if (invalid) return err(invalid);
 	if (mesh.vertices.length === 0 || mesh.triangles.length === 0) return err(validationError("VOXEL_EMPTY_MESH", "offsetMesh requires a non-empty triangle mesh."));
 	if (!Number.isFinite(distance)) return err(validationError("VOXEL_INVALID_DISTANCE", "distance must be a finite number."));
-	const params = resolveGridParams$1(opts);
+	const params = resolveGridParams$2(opts);
 	if (isErr(params)) return params;
-	const engine = resolveEngine(id);
+	const engine = resolveEngine$1(id);
 	if (isErr(engine)) return engine;
 	try {
 		try {
@@ -370,9 +370,9 @@ function shellMesh(mesh, thickness, opts, id) {
 	if (invalid) return err(invalid);
 	if (mesh.vertices.length === 0 || mesh.triangles.length === 0) return err(validationError("VOXEL_EMPTY_MESH", "shellMesh requires a non-empty triangle mesh."));
 	if (!Number.isFinite(thickness) || thickness <= 0) return err(validationError("VOXEL_INVALID_THICKNESS", "thickness must be a finite number > 0."));
-	const params = resolveGridParams$1(opts);
+	const params = resolveGridParams$2(opts);
 	if (isErr(params)) return params;
-	const engine = resolveEngine(id);
+	const engine = resolveEngine$1(id);
 	if (isErr(engine)) return engine;
 	try {
 		try {
@@ -402,9 +402,9 @@ function voxelBoolean(a, b, op, opts, id) {
 	if (a.vertices.length === 0 || a.triangles.length === 0) return err(validationError("VOXEL_EMPTY_MESH", "voxelBoolean requires a non-empty mesh for operand A."));
 	if (b.vertices.length === 0 || b.triangles.length === 0) return err(validationError("VOXEL_EMPTY_MESH", "voxelBoolean requires a non-empty mesh for operand B."));
 	const opCode = BOOLEAN_OP_CODES$1[op];
-	const params = resolveGridParams$1(opts);
+	const params = resolveGridParams$2(opts);
 	if (isErr(params)) return params;
-	const engine = resolveEngine(id);
+	const engine = resolveEngine$1(id);
 	if (isErr(engine)) return engine;
 	try {
 		try {
@@ -453,16 +453,16 @@ function voxelBooleanShapes(a, b, op, opts, id) {
 }
 //#endregion
 //#region src/voxel/fieldFns.ts
-var DEFAULT_RESOLUTION$1 = 48;
-var DEFAULT_PADDING$1 = 2;
+var DEFAULT_RESOLUTION$2 = 48;
+var DEFAULT_PADDING$2 = 2;
 var BOOLEAN_OP_CODES = {
 	union: 0,
 	intersection: 1,
 	difference: 2
 };
-function resolveGridParams(opts) {
-	const resolution = opts?.resolution ?? DEFAULT_RESOLUTION$1;
-	const padding = opts?.padding ?? DEFAULT_PADDING$1;
+function resolveGridParams$1(opts) {
+	const resolution = opts?.resolution ?? DEFAULT_RESOLUTION$2;
+	const padding = opts?.padding ?? DEFAULT_PADDING$2;
 	if (!Number.isInteger(resolution) || resolution < 1) return err(validationError("VOXEL_INVALID_RESOLUTION", "resolution must be an integer >= 1."));
 	if (!Number.isInteger(padding) || padding < 1) return err(validationError("VOXEL_INVALID_PADDING", "padding must be an integer >= 1."));
 	return ok({
@@ -512,6 +512,9 @@ function fieldDeletable(raw) {
 		}
 	};
 }
+function makeFieldHandle(raw) {
+	return makeHandle(raw);
+}
 function makeHandle(raw) {
 	const inner = createKernelHandle(fieldDeletable(raw));
 	const handle = {
@@ -568,9 +571,9 @@ function voxelField(mesh, opts, id) {
 	const invalid = validateMesh(mesh);
 	if (invalid) return err(invalid);
 	if (mesh.vertices.length === 0 || mesh.triangles.length === 0) return err(validationError("VOXEL_EMPTY_MESH", "voxelField requires a non-empty triangle mesh."));
-	const params = resolveGridParams(opts);
+	const params = resolveGridParams$1(opts);
 	if (isErr(params)) return params;
-	const engine = resolveEngine(id);
+	const engine = resolveEngine$1(id);
 	if (isErr(engine)) return engine;
 	try {
 		return ok(makeHandle(new engine.value.VoxelField(mesh.vertices, mesh.triangles, params.value.resolution, params.value.padding)));
@@ -597,9 +600,9 @@ function voxelBooleanField(a, b, op, opts, id) {
 	if (invalidB) return err(invalidB);
 	if (a.vertices.length === 0 || a.triangles.length === 0) return err(validationError("VOXEL_EMPTY_MESH", "voxelBooleanField requires a non-empty mesh for operand A."));
 	if (b.vertices.length === 0 || b.triangles.length === 0) return err(validationError("VOXEL_EMPTY_MESH", "voxelBooleanField requires a non-empty mesh for operand B."));
-	const params = resolveGridParams(opts);
+	const params = resolveGridParams$1(opts);
 	if (isErr(params)) return params;
-	const engine = resolveEngine(id);
+	const engine = resolveEngine$1(id);
 	if (isErr(engine)) return engine;
 	try {
 		return ok(makeHandle(engine.value.VoxelField.boolean_of(a.vertices, a.triangles, b.vertices, b.triangles, BOOLEAN_OP_CODES[op], params.value.resolution, params.value.padding)));
@@ -709,6 +712,163 @@ function voxelBooleanFieldShapes(a, b, op, opts, id) {
 	return voxelBooleanField(meshA.value, meshB.value, op, opts, id);
 }
 //#endregion
+//#region src/implicit/sdfFns.ts
+var DEFAULT_RESOLUTION$1 = 48;
+var DEFAULT_PADDING$1 = 2;
+function resolveEngine(id) {
+	try {
+		return ok(getVoxel(id));
+	} catch (cause) {
+		return err(moduleInitError("VOXEL_NOT_INITIALIZED", cause instanceof Error ? cause.message : "voxel engine not initialized", cause));
+	}
+}
+function resolveGridParams(opts) {
+	const resolution = opts?.resolution ?? DEFAULT_RESOLUTION$1;
+	const padding = opts?.padding ?? DEFAULT_PADDING$1;
+	if (!Number.isInteger(resolution) || resolution < 1) return err(validationError("SDF_INVALID_RESOLUTION", "resolution must be an integer >= 1."));
+	if (!Number.isInteger(padding) || padding < 1) return err(validationError("SDF_INVALID_PADDING", "padding must be an integer >= 1."));
+	return ok({
+		resolution,
+		padding
+	});
+}
+function sdfDeletable(raw) {
+	return {
+		raw,
+		delete() {
+			raw.free();
+		}
+	};
+}
+function makeSdfHandle(raw) {
+	const inner = createKernelHandle(sdfDeletable(raw));
+	return {
+		get value() {
+			return inner.value.raw;
+		},
+		get disposed() {
+			return inner.disposed;
+		},
+		[Symbol.dispose]() {
+			inner[Symbol.dispose]();
+		},
+		union(other) {
+			return makeSdfHandle(this.value.union(other.value));
+		},
+		intersection(other) {
+			return makeSdfHandle(this.value.intersection(other.value));
+		},
+		difference(other) {
+			return makeSdfHandle(this.value.difference(other.value));
+		},
+		smoothUnion(other, k) {
+			return makeSdfHandle(this.value.smooth_union(other.value, k));
+		},
+		smoothIntersection(other, k) {
+			return makeSdfHandle(this.value.smooth_intersection(other.value, k));
+		},
+		smoothDifference(other, k) {
+			return makeSdfHandle(this.value.smooth_difference(other.value, k));
+		},
+		offset(distance) {
+			return makeSdfHandle(this.value.offset(distance));
+		},
+		round(radius) {
+			return makeSdfHandle(this.value.round(radius));
+		},
+		shell(thickness) {
+			return makeSdfHandle(this.value.shell(thickness));
+		},
+		onion(thickness) {
+			return makeSdfHandle(this.value.onion(thickness));
+		},
+		translate(x, y, z) {
+			return makeSdfHandle(this.value.translate(x, y, z));
+		},
+		rotate(ax, ay, az, angle) {
+			return makeSdfHandle(this.value.rotate(ax, ay, az, angle));
+		},
+		scale(s) {
+			return makeSdfHandle(this.value.scale(s));
+		},
+		rasterize(opts) {
+			return rasterizeField(this.value, opts);
+		},
+		rasterizeIn(bounds, opts) {
+			return rasterizeFieldIn(this.value, bounds, opts);
+		}
+	};
+}
+function rasterizeField(sdf, opts) {
+	const params = resolveGridParams(opts);
+	if (isErr(params)) return params;
+	try {
+		return ok(makeFieldHandle(sdf.rasterize(params.value.resolution, params.value.padding)));
+	} catch (cause) {
+		return rasterizeError(cause);
+	}
+}
+function rasterizeFieldIn(sdf, bounds, opts) {
+	const params = resolveGridParams(opts);
+	if (isErr(params)) return params;
+	const invalid = validateBounds(bounds);
+	if (invalid) return err(invalid);
+	try {
+		return ok(makeFieldHandle(sdf.rasterize_in(bounds.min[0], bounds.min[1], bounds.min[2], bounds.max[0], bounds.max[1], bounds.max[2], params.value.resolution, params.value.padding)));
+	} catch (cause) {
+		return rasterizeError(cause);
+	}
+}
+function rasterizeError(cause) {
+	return err(computationError("SDF_RASTERIZE_FAILED", cause instanceof Error ? cause.message : "SDF rasterization failed (grid too large?).", cause));
+}
+function validateBounds(bounds) {
+	if (![...bounds.min, ...bounds.max].every((v) => Number.isFinite(v))) return validationError("SDF_INVALID_BOUNDS", "bounds must be finite numbers.");
+	for (let axis = 0; axis < 3; axis++) if (bounds.max[axis] <= bounds.min[axis]) return validationError("SDF_INVALID_BOUNDS", "bounds max must exceed min on every axis.");
+	return null;
+}
+function build(make, id) {
+	const engine = resolveEngine(id);
+	if (isErr(engine)) return engine;
+	try {
+		return ok(makeSdfHandle(make(engine.value)));
+	} catch (cause) {
+		return err(computationError("SDF_BUILD_FAILED", cause instanceof Error ? cause.message : "SDF primitive construction failed.", cause));
+	}
+}
+/** A sphere of radius `r`, centered at the origin. */
+function sphere(r, id) {
+	return build((e) => e.Sdf.sphere(r), id);
+}
+/** An axis-aligned box of half-extents `(hx, hy, hz)`, centered at the origin. */
+function box$1(hx, hy, hz, id) {
+	return build((e) => e.Sdf.box_(hx, hy, hz), id);
+}
+/** A box with rounded edges of radius `r`. */
+function roundedBox(hx, hy, hz, r, id) {
+	return build((e) => e.Sdf.rounded_box(hx, hy, hz, r), id);
+}
+/** A capped cylinder, axis +Z, radius `r`, total height `h`, centered at origin. */
+function cylinder$1(r, h, id) {
+	return build((e) => e.Sdf.cylinder(r, h), id);
+}
+/** A capped cone centered at the origin: base radius `r` at z = −h/2 tapering to an apex at z = +h/2. */
+function cone$1(r, h, id) {
+	return build((e) => e.Sdf.cone(r, h), id);
+}
+/** A capsule: a line segment `a`→`b` of radius `r`. */
+function capsule(a, b, r, id) {
+	return build((e) => e.Sdf.capsule(a[0], a[1], a[2], b[0], b[1], b[2], r), id);
+}
+/** A torus in the XY plane (axis +Z): a `minor`-radius tube on a `major` circle. */
+function torus(major, minor, id) {
+	return build((e) => e.Sdf.torus(major, minor), id);
+}
+/** A half-space: the plane through `h·n` with outward normal `n` (normalized). */
+function plane(n, h, id) {
+	return build((e) => e.Sdf.plane(n[0], n[1], n[2], h), id);
+}
+//#endregion
 //#region src/lattice/latticeFns.ts
 var LATTICE_TAGS = {
 	gyroid: 0,
@@ -762,7 +922,7 @@ function latticeInfill(mesh, opts, id) {
 	if (mesh.vertices.length === 0 || mesh.triangles.length === 0) return err(validationError("LATTICE_EMPTY_MESH", "latticeInfill requires a non-empty triangle mesh."));
 	const resolved = validateOptions(opts);
 	if (isErr(resolved)) return resolved;
-	const engine = resolveEngine(id);
+	const engine = resolveEngine$1(id);
 	if (isErr(engine)) return engine;
 	const { tag, period, thickness, resolution, padding } = resolved.value;
 	try {
@@ -803,7 +963,7 @@ function tpmsLattice(bounds, opts, id) {
 	}
 	const resolved = validateOptions(opts);
 	if (isErr(resolved)) return resolved;
-	const engine = resolveEngine(id);
+	const engine = resolveEngine$1(id);
 	if (isErr(engine)) return engine;
 	const { tag, period, thickness, resolution, padding } = resolved.value;
 	const [minX, minY, minZ] = bounds.min;
@@ -4032,11 +4192,11 @@ var primitives_exports = /* @__PURE__ */ __exportAll({
 	polyhedron: () => polyhedron,
 	sewShells: () => sewShells,
 	solid: () => solid,
-	sphere: () => sphere,
+	sphere: () => sphere$1,
 	subFace: () => subFace,
 	tangentArc: () => tangentArc,
 	threePointArc: () => threePointArc,
-	torus: () => torus,
+	torus: () => torus$1,
 	vertex: () => vertex,
 	wire: () => wire,
 	wireLoop: () => wireLoop
@@ -4464,7 +4624,7 @@ function depsOf(...sources) {
 	}
 	return acc.size === 0 ? EMPTY_DEPS : acc;
 }
-function box$1(x, y, z) {
+function box$2(x, y, z) {
 	const xe = asScalarExpr(x);
 	const ye = asScalarExpr(y);
 	const ze = asScalarExpr(z);
@@ -4479,7 +4639,7 @@ function box$1(x, y, z) {
 		freeParams: depsOf(xe, ye, ze)
 	};
 }
-function sphere$1(radius) {
+function sphere$2(radius) {
 	const re = asScalarExpr(radius);
 	return {
 		kind: "Sphere",
@@ -4488,7 +4648,7 @@ function sphere$1(radius) {
 		freeParams: re.freeParams
 	};
 }
-function cylinder$1(radius, height) {
+function cylinder$2(radius, height) {
 	const re = asScalarExpr(radius);
 	const he = asScalarExpr(height);
 	return {
@@ -4499,7 +4659,7 @@ function cylinder$1(radius, height) {
 		freeParams: depsOf(re, he)
 	};
 }
-function cone$1(radius1, radius2, height) {
+function cone$2(radius1, radius2, height) {
 	const r1 = asScalarExpr(radius1);
 	const r2 = asScalarExpr(radius2);
 	const he = asScalarExpr(height);
@@ -4514,7 +4674,7 @@ function cone$1(radius1, radius2, height) {
 		freeParams: depsOf(r1, r2, he)
 	};
 }
-function torus$1(majorRadius, minorRadius) {
+function torus$2(majorRadius, minorRadius) {
 	const ma = asScalarExpr(majorRadius);
 	const mi = asScalarExpr(minorRadius);
 	return {
@@ -4750,7 +4910,7 @@ function evalBox(node, ctx) {
 function evalSphere(node, ctx) {
 	const r = evalScalar(node.radius, ctx.env, "Sphere.radius");
 	if (!r.ok) return r;
-	return ok(sphere(r.value));
+	return ok(sphere$1(r.value));
 }
 function evalCylinder(node, ctx) {
 	const r = evalScalar(node.radius, ctx.env, "Cylinder.radius");
@@ -4773,7 +4933,7 @@ function evalTorus(node, ctx) {
 	if (!ma.ok) return ma;
 	const mi = evalScalar(node.minorRadius, ctx.env, "Torus.minorRadius");
 	if (!mi.ok) return mi;
-	return ok(torus(ma.value, mi.value));
+	return ok(torus$1(ma.value, mi.value));
 }
 function evalPolygon(node, ctx) {
 	const pts = [];
@@ -5419,7 +5579,7 @@ function readSingleExpr(j, key, build) {
 function readPrimitive(kind, j) {
 	switch (kind) {
 		case "Box": return readBox(j);
-		case "Sphere": return readSingleExpr(j, "radius", sphere$1);
+		case "Sphere": return readSingleExpr(j, "radius", sphere$2);
 		case "Cylinder": return readCylinder(j);
 		case "Cone": return readCone(j);
 		case "Torus": return readTorus(j);
@@ -5438,14 +5598,14 @@ function readBox(j) {
 	if (!y.ok) return y;
 	const z = readExpr(j["z"]);
 	if (!z.ok) return z;
-	return ok(box$1(x.value, y.value, z.value));
+	return ok(box$2(x.value, y.value, z.value));
 }
 function readCylinder(j) {
 	const r = readExpr(j["radius"]);
 	if (!r.ok) return r;
 	const h = readExpr(j["height"]);
 	if (!h.ok) return h;
-	return ok(cylinder$1(r.value, h.value));
+	return ok(cylinder$2(r.value, h.value));
 }
 function readCone(j) {
 	const r1 = readExpr(j["radius1"]);
@@ -5454,14 +5614,14 @@ function readCone(j) {
 	if (!r2.ok) return r2;
 	const h = readExpr(j["height"]);
 	if (!h.ok) return h;
-	return ok(cone$1(r1.value, r2.value, h.value));
+	return ok(cone$2(r1.value, r2.value, h.value));
 }
 function readTorus(j) {
 	const ma = readExpr(j["majorRadius"]);
 	if (!ma.ok) return ma;
 	const mi = readExpr(j["minorRadius"]);
 	if (!mi.ok) return mi;
-	return ok(torus$1(ma.value, mi.value));
+	return ok(torus$2(ma.value, mi.value));
 }
 function readPolygon(j) {
 	const pts = j["points"];
@@ -5637,11 +5797,11 @@ function foldBuildVec(dim, comps) {
 }
 function optimizeNode(n) {
 	switch (n.kind) {
-		case "Box": return box$1(foldExpr(n.x), foldExpr(n.y), foldExpr(n.z));
-		case "Sphere": return sphere$1(foldExpr(n.radius));
-		case "Cylinder": return cylinder$1(foldExpr(n.radius), foldExpr(n.height));
-		case "Cone": return cone$1(foldExpr(n.radius1), foldExpr(n.radius2), foldExpr(n.height));
-		case "Torus": return torus$1(foldExpr(n.majorRadius), foldExpr(n.minorRadius));
+		case "Box": return box$2(foldExpr(n.x), foldExpr(n.y), foldExpr(n.z));
+		case "Sphere": return sphere$2(foldExpr(n.radius));
+		case "Cylinder": return cylinder$2(foldExpr(n.radius), foldExpr(n.height));
+		case "Cone": return cone$2(foldExpr(n.radius1), foldExpr(n.radius2), foldExpr(n.height));
+		case "Torus": return torus$2(foldExpr(n.majorRadius), foldExpr(n.minorRadius));
 		case "Polygon": return polygon$1(n.points.map(foldExpr));
 		case "Circle": return circle$1(foldExpr(n.radius));
 		case "Line": return line$1(foldExpr(n.from), foldExpr(n.to));
@@ -5799,15 +5959,15 @@ var csg_exports = /* @__PURE__ */ __exportAll({
 	asVec2Expr: () => asVec2Expr,
 	asVec3Expr: () => asVec3Expr,
 	binOp: () => binOp,
-	box: () => box$1,
+	box: () => box$2,
 	buildVec: () => buildVec,
 	circle: () => circle$1,
 	component: () => component,
 	compound: () => compound$1,
-	cone: () => cone$1,
+	cone: () => cone$2,
 	cut: () => cut$1,
 	cutAll: () => cutAll$1,
-	cylinder: () => cylinder$1,
+	cylinder: () => cylinder$2,
 	emptyFace: () => emptyFace,
 	emptySolid: () => emptySolid,
 	emptyWire: () => emptyWire,
@@ -5829,9 +5989,9 @@ var csg_exports = /* @__PURE__ */ __exportAll({
 	replaceNode: () => replaceNode,
 	rotate: () => rotate$1,
 	scale: () => scale$1,
-	sphere: () => sphere$1,
+	sphere: () => sphere$2,
 	toJSON: () => toJSON,
-	torus: () => torus$1,
+	torus: () => torus$2,
 	translate: () => translate$1,
 	unaryOp: () => unaryOp,
 	vec2Lit: () => vec2Lit,
@@ -5840,4 +6000,4 @@ var csg_exports = /* @__PURE__ */ __exportAll({
 	withEvaluator: () => withEvaluator
 });
 //#endregion
-export { BaseSketcher2d, BlueprintSketcher, BrepBugError, BrepErrorCode, BrepWrapperError, BrepkitAdapter, CompoundSketch, DEFAULT_CAPABILITIES, DEG2RAD, DisposalScope, EXACT_BREP_CAPABILITIES, FaceSketcher, HASH_CODE_MAX, OK, OcctWasmAdapter, RAD2DEG, Sketch, Sketcher, Sketches, addChild, addHoles, addMate, addStep, adjacentFaces, all, andThen, applyGlue, applyMatrix, approximateCurve, as2D, as3D, asTopo, assignRoles, autoHeal, bezier, blueprintToDXF, booleanPipeline, booleans_exports as booleans, boss, box, bsplineApprox, bug, cameraFromPlane, cameraLookAt, captureHint, cast, castShape, castShape3D, chamfer, chamferDistAngle as chamferDistAngleShape, chamferWithEvolution, checkAllInterferences, checkBoolean, checkInterference, circle, circularPattern, classifyPointOnFace, clearMeshCache, clone, closedWire, collect, collectShapes, colorFaces, colorShape, complexExtrude, composeTransforms, compound, compoundSketchExtrude, compoundSketchFace, compoundSketchLoft, compoundSketchRevolve, computationError, computeStraightSkeleton, cone, construction_exports as construction, convexHull, cornerFinder, countNodes, createAssembly, createAssemblyNode, createBlueprint, createCamera, createCompound, createCompoundBlueprint, createDistanceQuery, createEdge, createFace, createHandle, createHistory, createKernelHandle, createMeshCache, createNamedPlane, createOperationRegistry, createPlane, createRef, createRegistry, createShell, createSolid, createTaskQueue, createVertex, createWire, createWorkerClient, createWorkerHandler, csg_exports as csg, currentQuality, curve2dBoundingBox, curve2dDistanceFrom, curve2dFirstPoint, curve2dIsOnCurve, curve2dLastPoint, curve2dParameter, curve2dSplitAt, curve2dTangentAt, curveEndPoint, curveIsClosed, curveIsPeriodic, curveLength, curvePeriod, curvePointAt, curveStartPoint, curveTangentAt, cut, cut2D, cutAll, cutAllBisect, cutBlueprints, cutWithEvolution, cylinder, defaultScorer, dequeueTask, describe, deserializeDrawing, deserializeHistory, fromBREP as deserializeShape, downcast, draft, draw, drawCircle, drawEllipse, drawFaceOutline, drawParametricFunction, drawPointsInterpolation, drawPolysides, drawProjection, drawRectangle, drawRoundedRectangle, drawSingleCircle, drawSingleEllipse, drawText, drawingChamfer, drawingCut, drawingFillet, drawingFuse, drawingIntersect, drawingToSketchOnPlane, drill, edgeFinder, edgesOfFace, ellipse, ellipseArc, ellipsoid, enqueueTask, err, exportAssemblySTEP, exportDXF, exportGlb, exportGltf, exportIGES, exportOBJ, exportSTEP, exportSTEPConfigured, exportSTL, exportThreeMF, extrude, extrudeAll, face, faceCenter, faceFinder, faceGeomType, faceOrientation, facesOfEdge, fieldBoolean, fieldContour, fieldOffset, fieldReinit, fieldShell, fill, filledFace, fillet, filletWithEvolution, findFacesByTag, findNode, findStep, fixSelfIntersection, fixShape, flatMap, flatten, flipFaceOrientation, flipOrientation, fontMetrics, fromBREP$1 as fromBREP, fromKernelDir, fromKernelPnt, fromKernelVec, fromNullable, fuse, fuse2D, fuseAll, fuseAllBisect, fuseBlueprints, fuseWithEvolution, gearGeometry, getActiveVoxelId, getBounds, getBounds2D, getCompSolids, getCurveType, getDisposalStats, getEdges, getFaceColor, getFaceOrigins, getFaceTags, getFaces, getFont, getHashCode, getShape as getHistoryShape, getKernel, getKernelCapabilities, getKernelTier, getNurbsCurveData, getNurbsSurfaceData, getOrientation, getOrientation2D, getPerformanceStats, getShapeColor, getShapeKind, getShells, getSingleFace, getSolids, getSurfaceType, getTagMetadata, getVertices, getVoxel, getWires, guidedSweep, heal, healFace, healSolid, healWire, helix, hull, importDXF, importGLB, importIGES, importOBJ, importSTEP, importSTL, importSVG, importSVGPathD, importThreeMF, init, initFromManifold, initFromOC, initVoxel, innerWires, interpolateCurve, intersect, intersect2D, intersectBlueprints, intersectWithEvolution, invalidateShapeCache, ioNs_exports as io, ioError, is2D, is3D, isChamferRadius, isClosedWire, isCompSolid, isCompound, isDisposeRequest, isEdge, isEmpty, isEqualShape, isErr, isErrorResponse, isFace, isFilletRadius, isInitRequest, isInside2D, isLive, isManifoldShell, isNumber, isOk, isOperationRequest, isOrientedFace, isPlanarFace, isPlanarWire, isProjectionPlane, isEmpty$1 as isQueueEmpty, isSameShape, isShape1D, isShape3D, isShell, isSolid, isSuccessResponse, isValid, isValidSolid, isVertex, isWire, iterCompSolids, iterEdges, iterFaces, iterShells, iterSolids, iterTopo, iterVertices, iterWires, kernelCall, kernelCallRaw, kernelCallScoped, kernelError, latticeInfill, latticeInfillShape, line, linearPattern, loadFont, loft, loftAll, makeBaseBox, makeExternalGear, makeInternalGear, makePlane, makePlanetaryGear, makeProjectedEdges, manifoldShell, map, mapBoth, mapErr, match, measureArea, measureCurvatureAt, measureCurvatureAtMid, measureDistance, measureDistanceProps, measureLength, measureLinearProps, measureSurfaceProps, measureVolume, measureVolumeProps, measurement_exports as measurement, mesh, meshEdges, meshMultiLOD, minkowski, mirror, mirror2D, mirrorDrawing, mirrorJoin, modifiers_exports as modifiers, modifyStep, moduleInitError, multiSectionSweep, normalAt, offset, offsetFace, offsetMesh, offsetShape, offsetWire2D, ok, or, orElse, organiseBlueprints, orientedFace, outerWire, patterns_exports as patterns, pendingCount, pipeline, pivotPlane, planarFace, planarWire, planetPlacements, pocket, pointOnSurface, pointsInside, polygon, polyhedron, polysideInnerRadius, polysidesBlueprint, positionOnCurve, prewarm, primitives_exports as primitives, projectEdges, projectPointOnFace, query_exports as query, queryError, rectangularPattern, registerHandler, registerKernel, registerKernelTier, registerOperation, registerShape, registerVoxel, rejectAll, removeChild, removeHolesFromFace, repairMesh, replayFrom, replayHistory, resetDisposalStats, resetPerformanceStats, resize, resolve, resolve3D, resolveDirection, resolvePlane, resolveRef, reverseCurve, revolve, roof, rotate, rotate2D, rotateDrawing, roundedRectangleBlueprint, scale, scale2D, scaleDrawing, section, sectionToFace, serializeHistory, setShapeOrigin, setTagMetadata, sewShells, shape, shapeToMeshInput, shapeType, sharedEdges, shell, shellMesh, shellShape, shellWithEvolution, simplify, sketchCircle, sketchEllipse, sketchExtrude, sketchFace, sketchFaceOffset, sketchHelix, sketchLoft, sketchOnFace2D, sketchOnPlane2D, sketchParametricFunction, sketchPolysides, sketchRectangle, sketchRevolve, sketchRoundedRectangle, sketchSweep, sketchText, sketchWires, sketcherStateError, slice, solid, solidFromShell, solveAssembly, sphere, split, stepCount, stepsFrom, stretch2D, subFace, supportExtrude, supportsConstraintSketch, supportsProjection, surfaceFromGrid, surfaceFromImage, sweep, tagFaces, tangentArc, tap, tapErr, textBlueprints, textMetrics, thicken, threePointArc, toBREP, toBufferGeometryData, toGroupedBufferGeometryData, toKernelVec, toLODGeometryData, toLineGeometryData, toSVGPathD, toVec2, toVec3, torus, tpmsLattice, transformCopy, transforms_exports as transforms, translate, translate2D, translateDrawing, translatePlane, tryCatch, tryCatchAsync, twistExtrude, typeCastError, undoLast, unsupportedError, unwrap, unwrapErr, unwrapOr, unwrapOrElse, updateNode, updateRoles, uvBounds, uvCoordinates, validSolid, validatePlanetary, validationError, variableFillet, vecAdd, vecAngle, vecCross, vecDistance, vecDot, vecEquals, vecIsZero, vecLength, vecLengthSq, vecNegate, vecNormalize, vecProjectToPlane, vecRepr, vecRotate, vecScale, vecSub, vertex, vertexFinder, vertexPosition, verticesOfEdge, voxelBoolean, voxelBooleanField, voxelBooleanFieldShapes, voxelBooleanShapes, voxelField, voxelFieldFromShape, walkAssembly, windingNumbers, wire, wireFinder, wireLoop, wiresOfFace, withKernel, withKernelDir, withKernelPnt, withKernelVec, withQuality, withScope, withScopeResult, withScopeResultAsync, withTier, zip as zipResults };
+export { BaseSketcher2d, BlueprintSketcher, BrepBugError, BrepErrorCode, BrepWrapperError, BrepkitAdapter, CompoundSketch, DEFAULT_CAPABILITIES, DEG2RAD, DisposalScope, EXACT_BREP_CAPABILITIES, FaceSketcher, HASH_CODE_MAX, OK, OcctWasmAdapter, RAD2DEG, Sketch, Sketcher, Sketches, addChild, addHoles, addMate, addStep, adjacentFaces, all, andThen, applyGlue, applyMatrix, approximateCurve, as2D, as3D, asTopo, assignRoles, autoHeal, bezier, blueprintToDXF, booleanPipeline, booleans_exports as booleans, boss, box, bsplineApprox, bug, cameraFromPlane, cameraLookAt, captureHint, cast, castShape, castShape3D, chamfer, chamferDistAngle as chamferDistAngleShape, chamferWithEvolution, checkAllInterferences, checkBoolean, checkInterference, circle, circularPattern, classifyPointOnFace, clearMeshCache, clone, closedWire, collect, collectShapes, colorFaces, colorShape, complexExtrude, composeTransforms, compound, compoundSketchExtrude, compoundSketchFace, compoundSketchLoft, compoundSketchRevolve, computationError, computeStraightSkeleton, cone, construction_exports as construction, convexHull, cornerFinder, countNodes, createAssembly, createAssemblyNode, createBlueprint, createCamera, createCompound, createCompoundBlueprint, createDistanceQuery, createEdge, createFace, createHandle, createHistory, createKernelHandle, createMeshCache, createNamedPlane, createOperationRegistry, createPlane, createRef, createRegistry, createShell, createSolid, createTaskQueue, createVertex, createWire, createWorkerClient, createWorkerHandler, csg_exports as csg, currentQuality, curve2dBoundingBox, curve2dDistanceFrom, curve2dFirstPoint, curve2dIsOnCurve, curve2dLastPoint, curve2dParameter, curve2dSplitAt, curve2dTangentAt, curveEndPoint, curveIsClosed, curveIsPeriodic, curveLength, curvePeriod, curvePointAt, curveStartPoint, curveTangentAt, cut, cut2D, cutAll, cutAllBisect, cutBlueprints, cutWithEvolution, cylinder, defaultScorer, dequeueTask, describe, deserializeDrawing, deserializeHistory, fromBREP as deserializeShape, downcast, draft, draw, drawCircle, drawEllipse, drawFaceOutline, drawParametricFunction, drawPointsInterpolation, drawPolysides, drawProjection, drawRectangle, drawRoundedRectangle, drawSingleCircle, drawSingleEllipse, drawText, drawingChamfer, drawingCut, drawingFillet, drawingFuse, drawingIntersect, drawingToSketchOnPlane, drill, edgeFinder, edgesOfFace, ellipse, ellipseArc, ellipsoid, enqueueTask, err, exportAssemblySTEP, exportDXF, exportGlb, exportGltf, exportIGES, exportOBJ, exportSTEP, exportSTEPConfigured, exportSTL, exportThreeMF, extrude, extrudeAll, face, faceCenter, faceFinder, faceGeomType, faceOrientation, facesOfEdge, fieldBoolean, fieldContour, fieldOffset, fieldReinit, fieldShell, fill, filledFace, fillet, filletWithEvolution, findFacesByTag, findNode, findStep, fixSelfIntersection, fixShape, flatMap, flatten, flipFaceOrientation, flipOrientation, fontMetrics, fromBREP$1 as fromBREP, fromKernelDir, fromKernelPnt, fromKernelVec, fromNullable, fuse, fuse2D, fuseAll, fuseAllBisect, fuseBlueprints, fuseWithEvolution, gearGeometry, getActiveVoxelId, getBounds, getBounds2D, getCompSolids, getCurveType, getDisposalStats, getEdges, getFaceColor, getFaceOrigins, getFaceTags, getFaces, getFont, getHashCode, getShape as getHistoryShape, getKernel, getKernelCapabilities, getKernelTier, getNurbsCurveData, getNurbsSurfaceData, getOrientation, getOrientation2D, getPerformanceStats, getShapeColor, getShapeKind, getShells, getSingleFace, getSolids, getSurfaceType, getTagMetadata, getVertices, getVoxel, getWires, guidedSweep, heal, healFace, healSolid, healWire, helix, hull, importDXF, importGLB, importIGES, importOBJ, importSTEP, importSTL, importSVG, importSVGPathD, importThreeMF, init, initFromManifold, initFromOC, initVoxel, innerWires, interpolateCurve, intersect, intersect2D, intersectBlueprints, intersectWithEvolution, invalidateShapeCache, ioNs_exports as io, ioError, is2D, is3D, isChamferRadius, isClosedWire, isCompSolid, isCompound, isDisposeRequest, isEdge, isEmpty, isEqualShape, isErr, isErrorResponse, isFace, isFilletRadius, isInitRequest, isInside2D, isLive, isManifoldShell, isNumber, isOk, isOperationRequest, isOrientedFace, isPlanarFace, isPlanarWire, isProjectionPlane, isEmpty$1 as isQueueEmpty, isSameShape, isShape1D, isShape3D, isShell, isSolid, isSuccessResponse, isValid, isValidSolid, isVertex, isWire, iterCompSolids, iterEdges, iterFaces, iterShells, iterSolids, iterTopo, iterVertices, iterWires, kernelCall, kernelCallRaw, kernelCallScoped, kernelError, latticeInfill, latticeInfillShape, line, linearPattern, loadFont, loft, loftAll, makeBaseBox, makeExternalGear, makeInternalGear, makePlane, makePlanetaryGear, makeProjectedEdges, manifoldShell, map, mapBoth, mapErr, match, measureArea, measureCurvatureAt, measureCurvatureAtMid, measureDistance, measureDistanceProps, measureLength, measureLinearProps, measureSurfaceProps, measureVolume, measureVolumeProps, measurement_exports as measurement, mesh, meshEdges, meshMultiLOD, minkowski, mirror, mirror2D, mirrorDrawing, mirrorJoin, modifiers_exports as modifiers, modifyStep, moduleInitError, multiSectionSweep, normalAt, offset, offsetFace, offsetMesh, offsetShape, offsetWire2D, ok, or, orElse, organiseBlueprints, orientedFace, outerWire, patterns_exports as patterns, pendingCount, pipeline, pivotPlane, planarFace, planarWire, planetPlacements, pocket, pointOnSurface, pointsInside, polygon, polyhedron, polysideInnerRadius, polysidesBlueprint, positionOnCurve, prewarm, primitives_exports as primitives, projectEdges, projectPointOnFace, query_exports as query, queryError, rectangularPattern, registerHandler, registerKernel, registerKernelTier, registerOperation, registerShape, registerVoxel, rejectAll, removeChild, removeHolesFromFace, repairMesh, replayFrom, replayHistory, resetDisposalStats, resetPerformanceStats, resize, resolve, resolve3D, resolveDirection, resolvePlane, resolveRef, reverseCurve, revolve, roof, rotate, rotate2D, rotateDrawing, roundedRectangleBlueprint, scale, scale2D, scaleDrawing, box$1 as sdfBox, capsule as sdfCapsule, cone$1 as sdfCone, cylinder$1 as sdfCylinder, plane as sdfPlane, roundedBox as sdfRoundedBox, sphere as sdfSphere, torus as sdfTorus, section, sectionToFace, serializeHistory, setShapeOrigin, setTagMetadata, sewShells, shape, shapeToMeshInput, shapeType, sharedEdges, shell, shellMesh, shellShape, shellWithEvolution, simplify, sketchCircle, sketchEllipse, sketchExtrude, sketchFace, sketchFaceOffset, sketchHelix, sketchLoft, sketchOnFace2D, sketchOnPlane2D, sketchParametricFunction, sketchPolysides, sketchRectangle, sketchRevolve, sketchRoundedRectangle, sketchSweep, sketchText, sketchWires, sketcherStateError, slice, solid, solidFromShell, solveAssembly, sphere$1 as sphere, split, stepCount, stepsFrom, stretch2D, subFace, supportExtrude, supportsConstraintSketch, supportsProjection, surfaceFromGrid, surfaceFromImage, sweep, tagFaces, tangentArc, tap, tapErr, textBlueprints, textMetrics, thicken, threePointArc, toBREP, toBufferGeometryData, toGroupedBufferGeometryData, toKernelVec, toLODGeometryData, toLineGeometryData, toSVGPathD, toVec2, toVec3, torus$1 as torus, tpmsLattice, transformCopy, transforms_exports as transforms, translate, translate2D, translateDrawing, translatePlane, tryCatch, tryCatchAsync, twistExtrude, typeCastError, undoLast, unsupportedError, unwrap, unwrapErr, unwrapOr, unwrapOrElse, updateNode, updateRoles, uvBounds, uvCoordinates, validSolid, validatePlanetary, validationError, variableFillet, vecAdd, vecAngle, vecCross, vecDistance, vecDot, vecEquals, vecIsZero, vecLength, vecLengthSq, vecNegate, vecNormalize, vecProjectToPlane, vecRepr, vecRotate, vecScale, vecSub, vertex, vertexFinder, vertexPosition, verticesOfEdge, voxelBoolean, voxelBooleanField, voxelBooleanFieldShapes, voxelBooleanShapes, voxelField, voxelFieldFromShape, walkAssembly, windingNumbers, wire, wireFinder, wireLoop, wiresOfFace, withKernel, withKernelDir, withKernelPnt, withKernelVec, withQuality, withScope, withScopeResult, withScopeResultAsync, withTier, zip as zipResults };

package/dist/implicit/index.d.ts

@@ -0,0 +1,11 @@
+/**
+ * Field-first implicit CAD domain (ADR-0013, brepjs-implicit Phase 1).
+ *
+ * An analytic SDF expression tree that rasterizes DIRECTLY into the voxel
+ * substrate's dense grid with no input mesh — the field-first twin of the
+ * mesh-first voxel path. Primitives compose through CSG / smooth / domain
+ * combinators into an {@link SdfHandle}, which rasterizes to a banded-SDF
+ * {@link VoxelFieldHandle} for contour / offset / shell.
+ */
+export type { SdfHandle, SdfBounds } from './sdfFns.js';
+export { sphere, box, roundedBox, cylinder, cone, capsule, torus, plane } from './sdfFns.js';

package/dist/implicit/sdfFns.d.ts

@@ -0,0 +1,59 @@
+import { Result } from '../core/result.js';
+import { WasmSdf } from '../voxel/engine.js';
+import { VoxelFieldHandle, VoxelFieldOptions } from '../voxel/fieldFns.js';
+/** Explicit world bounds `[min..max]` for {@link SdfHandle.rasterizeIn}. */
+export interface SdfBounds {
+    min: [number, number, number];
+    max: [number, number, number];
+}
+/**
+ * A disposable handle around an analytic SDF expression. Every combinator clones
+ * into a NEW wasm node and returns a NEW handle (the wasm `Sdf` is a value, not a
+ * mutable builder), so an `SdfHandle` is immutable — the receiver stays valid and
+ * must be disposed independently. `rasterize` produces a banded-SDF
+ * {@link VoxelFieldHandle} you can boolean / offset / shell / contour.
+ *
+ * Dispose is mandatory: use `using`, or call `[Symbol.dispose]()`, to free the
+ * WASM expression tree. Intermediate handles in a chain (`a.union(b).shell(t)`)
+ * each own a wasm allocation; bind them to `using` or dispose them explicitly.
+ */
+export interface SdfHandle {
+    /** The wrapped WASM expression. Throws if the handle has been disposed. */
+    readonly value: WasmSdf;
+    /** Whether the backing WASM expression has been freed. */
+    readonly disposed: boolean;
+    [Symbol.dispose](): void;
+    union(other: SdfHandle): SdfHandle;
+    intersection(other: SdfHandle): SdfHandle;
+    difference(other: SdfHandle): SdfHandle;
+    smoothUnion(other: SdfHandle, k: number): SdfHandle;
+    smoothIntersection(other: SdfHandle, k: number): SdfHandle;
+    smoothDifference(other: SdfHandle, k: number): SdfHandle;
+    offset(distance: number): SdfHandle;
+    round(radius: number): SdfHandle;
+    shell(thickness: number): SdfHandle;
+    onion(thickness: number): SdfHandle;
+    translate(x: number, y: number, z: number): SdfHandle;
+    rotate(ax: number, ay: number, az: number, angle: number): SdfHandle;
+    scale(s: number): SdfHandle;
+    /** Rasterize into a banded-SDF field over the expression's analytic bounds. */
+    rasterize(opts?: VoxelFieldOptions): Result<VoxelFieldHandle>;
+    /** Rasterize over explicit world bounds (clips unbounded primitives). */
+    rasterizeIn(bounds: SdfBounds, opts?: VoxelFieldOptions): Result<VoxelFieldHandle>;
+}
+/** A sphere of radius `r`, centered at the origin. */
+export declare function sphere(r: number, id?: string): Result<SdfHandle>;
+/** An axis-aligned box of half-extents `(hx, hy, hz)`, centered at the origin. */
+export declare function box(hx: number, hy: number, hz: number, id?: string): Result<SdfHandle>;
+/** A box with rounded edges of radius `r`. */
+export declare function roundedBox(hx: number, hy: number, hz: number, r: number, id?: string): Result<SdfHandle>;
+/** A capped cylinder, axis +Z, radius `r`, total height `h`, centered at origin. */
+export declare function cylinder(r: number, h: number, id?: string): Result<SdfHandle>;
+/** A capped cone centered at the origin: base radius `r` at z = −h/2 tapering to an apex at z = +h/2. */
+export declare function cone(r: number, h: number, id?: string): Result<SdfHandle>;
+/** A capsule: a line segment `a`→`b` of radius `r`. */
+export declare function capsule(a: [number, number, number], b: [number, number, number], r: number, id?: string): Result<SdfHandle>;
+/** A torus in the XY plane (axis +Z): a `minor`-radius tube on a `major` circle. */
+export declare function torus(major: number, minor: number, id?: string): Result<SdfHandle>;
+/** A half-space: the plane through `h·n` with outward normal `n` (normalized). */
+export declare function plane(n: [number, number, number], h: number, id?: string): Result<SdfHandle>;

package/dist/index.d.ts

@@ -40,6 +40,8 @@ export { importSVGPathD, importSVG, type SVGImportOptions } from './io/svgImport
 export { exportSTEPConfigured, type StepExportOptions, type StepExportPart, } from './io/stepConfigFns.js';
 export { initVoxel, registerVoxel, getVoxel, getActiveVoxelId, windingNumbers, pointsInside, repairMesh, offsetMesh, shellMesh, voxelBoolean, offsetShape, shellShape, voxelBooleanShapes, voxelField, voxelBooleanField, fieldBoolean, fieldOffset, fieldShell, fieldReinit, fieldContour, voxelFieldFromShape, voxelBooleanFieldShapes, shapeToMeshInput, } from './voxel/index.js';
 export type { VoxelEngine, VoxelMeshInput, VoxelRepairResult, RepairOptions, VoxelOpOptions, VoxelFieldHandle, VoxelFieldOptions, VoxelBooleanOp, } from './voxel/index.js';
+export { sphere as sdfSphere, box as sdfBox, roundedBox as sdfRoundedBox, cylinder as sdfCylinder, cone as sdfCone, capsule as sdfCapsule, torus as sdfTorus, plane as sdfPlane, } from './implicit/index.js';
+export type { SdfHandle, SdfBounds } from './implicit/index.js';
 export { latticeInfill, latticeInfillShape, tpmsLattice } from './lattice/index.js';
 export type { LatticeType, LatticeOptions, LatticeBounds } from './lattice/index.js';
 export { default as Sketcher } from './sketching/sketcher.js';

package/dist/voxel/engine.d.ts

@@ -29,9 +29,60 @@ export interface VoxelEngine {
      * satisfied by the generated `VoxelField` wasm-bindgen class.
      */
     VoxelField: WasmVoxelFieldConstructor;
+    /**
+     * Field-first analytic SDF builder (ADR-0013). Static primitive constructors
+     * and combinator methods compose an opaque expression tree that rasterizes
+     * directly into a {@link WasmVoxelField} with no input mesh. Structurally
+     * satisfied by the generated `Sdf` wasm-bindgen class.
+     */
+    Sdf: WasmSdfConstructor;
     /** Engine artifact version, for loader/artifact compatibility checks. */
     version(): string;
 }
+/**
+ * Static surface of the wasm `Sdf` class: the primitive constructors that seed an
+ * expression tree (centered at the origin unless noted). Each returns a fresh
+ * {@link WasmSdf}. Structurally satisfied by the generated `Sdf` class.
+ */
+export interface WasmSdfConstructor {
+    sphere(r: number): WasmSdf;
+    box_(hx: number, hy: number, hz: number): WasmSdf;
+    rounded_box(hx: number, hy: number, hz: number, r: number): WasmSdf;
+    cylinder(r: number, h: number): WasmSdf;
+    cone(r: number, h: number): WasmSdf;
+    capsule(ax: number, ay: number, az: number, bx: number, by: number, bz: number, r: number): WasmSdf;
+    torus(major: number, minor: number): WasmSdf;
+    plane(nx: number, ny: number, nz: number, h: number): WasmSdf;
+}
+/**
+ * An opaque analytic SDF expression. Every combinator CLONES into a new node and
+ * returns a fresh `WasmSdf` (wasm-bindgen has no shared borrow across calls), so an
+ * `Sdf` is a value, not a mutable builder. `rasterize` builds a banded-SDF
+ * {@link WasmVoxelField}; `free()` releases the backing WASM expression tree.
+ * Structurally satisfied by the generated `Sdf` wasm-bindgen class.
+ */
+export interface WasmSdf {
+    union(other: WasmSdf): WasmSdf;
+    intersection(other: WasmSdf): WasmSdf;
+    difference(other: WasmSdf): WasmSdf;
+    smooth_union(other: WasmSdf, k: number): WasmSdf;
+    smooth_intersection(other: WasmSdf, k: number): WasmSdf;
+    smooth_difference(other: WasmSdf, k: number): WasmSdf;
+    offset(d: number): WasmSdf;
+    round(r: number): WasmSdf;
+    shell(t: number): WasmSdf;
+    onion(t: number): WasmSdf;
+    translate(x: number, y: number, z: number): WasmSdf;
+    rotate(ax: number, ay: number, az: number, angle: number): WasmSdf;
+    scale(s: number): WasmSdf;
+    /** Rasterize into a banded-SDF dense field over the expression's analytic bounds. */
+    rasterize(resolution: number, padding: number): WasmVoxelField;
+    /** Rasterize over explicit `[min..max]` bounds (clips unbounded primitives). */
+    rasterize_in(min_x: number, min_y: number, min_z: number, max_x: number, max_y: number, max_z: number, resolution: number, padding: number): WasmVoxelField;
+    /** Release the backing WASM expression-tree allocation (wasm-bindgen lifecycle). */
+    free(): void;
+    [Symbol.dispose](): void;
+}
 /**
  * Constructor of the wasm `VoxelField` class. `new VoxelField(verts, tris, res,
  * padding)` voxelizes a mesh into a persistent dense field. Throws (as a JS

package/dist/voxel/fieldFns.d.ts

@@ -45,6 +45,7 @@ export interface VoxelFieldHandle {
     /** Surface-Nets contour the current field to a mesh (the field stays alive). */
     contour(): KernelMeshResult;
 }
+export declare function makeFieldHandle(raw: WasmVoxelField): VoxelFieldHandle;
 /**
  * Voxelize a mesh into a persistent dense {@link VoxelFieldHandle}: one grid you
  * can boolean / offset / shell / reinit in place, then contour once. The handle

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "brepjs",
-  "version": "18.67.0",
+  "version": "18.68.0",
   "description": "Web CAD library with pluggable geometry kernel",
   "keywords": [
     "cad",