@code3d/core 0.0.1-alpha.0

package/bld/library/sketch-curve-intersections.js

@@ -0,0 +1,117 @@
+import { sketchCurveClosestParameter, sketchCurvePosition, sketchCurveTolerance, } from './sketch-curves.js';
+const subtract = (a, b) => [
+    a[0] - b[0],
+    a[1] - b[1],
+];
+const cross = (a, b) => a[0] * b[1] - a[1] * b[0];
+const distance = (a, b) => Math.hypot(...subtract(a, b));
+const endpoints = (curve) => curve.kind === 'circle'
+    ? []
+    : [sketchCurvePosition(curve, 0), sketchCurvePosition(curve, 1)];
+/**
+ * Isolated intersections and finite overlap boundaries, ordered along first.
+ * A shared continuous interior does not enumerate points: coincident full
+ * circles have no boundaries. Parameters are local to this evaluation, not IDs.
+ */
+export function sketchCurveIntersections(first, second) {
+    const tolerances = [
+        sketchCurveTolerance(first),
+        sketchCurveTolerance(second),
+    ];
+    const tolerance = Math.min(...tolerances);
+    const contacts = [];
+    for (const position of candidates(first, second, tolerance)) {
+        const parameters = [
+            sketchCurveClosestParameter(first, position),
+            sketchCurveClosestParameter(second, position),
+        ];
+        if ([first, second].some((curve, i) => distance(position, sketchCurvePosition(curve, parameters[i])) >
+            tolerances[i]) ||
+            contacts.some(contact => distance(position, contact.position) <= tolerance))
+            continue;
+        contacts.push({ position, parameters });
+    }
+    return contacts.sort((a, b) => a.parameters[0] - b.parameters[0]);
+}
+function candidates(first, second, tolerance) {
+    if (first.kind === 'line') {
+        const [a, b] = first.points;
+        const length = distance(a, b);
+        if (!length)
+            return [a];
+        const unit = [
+            (b[0] - a[0]) / length,
+            (b[1] - a[1]) / length,
+        ];
+        if (second.kind === 'line') {
+            const [c, d] = second.points;
+            const otherLength = distance(c, d);
+            if (!otherLength)
+                return [c];
+            // Collinear endpoints delimit overlap. Check them before dividing by
+            // the angle, which can be roundoff-sized for rotated/translated lines.
+            if (second.points.every(p => Math.abs(cross(unit, subtract(p, a))) <= tolerance))
+                return [a, b, c, d];
+            const otherUnit = [
+                (d[0] - c[0]) / otherLength,
+                (d[1] - c[1]) / otherLength,
+            ];
+            const denominator = cross(unit, otherUnit);
+            if (Math.abs(denominator) <= Number.EPSILON * 16)
+                return [];
+            const along = cross(subtract(c, a), otherUnit) / denominator;
+            return [[a[0] + along * unit[0], a[1] + along * unit[1]]];
+        }
+        const perpendicular = cross(unit, subtract(second.center, a));
+        const gap = Math.abs(perpendicular);
+        if (gap > second.radius + tolerance)
+            return [];
+        const normalized = gap / second.radius;
+        // Classify tangency before sqrt magnifies coordinate roundoff into two
+        // artificial cuts. Use the same model-space tolerance as finite contacts.
+        const offset = Math.abs(gap - second.radius) <= tolerance
+            ? 0
+            : second.radius *
+                Math.sqrt(Math.max(0, (1 - normalized) * (1 + normalized)));
+        // Start at the circular center's projection, not a far-away line endpoint.
+        const foot = [
+            second.center[0] + perpendicular * unit[1],
+            second.center[1] - perpendicular * unit[0],
+        ];
+        return [-offset, offset].map(along => [
+            foot[0] + along * unit[0],
+            foot[1] + along * unit[1],
+        ]);
+    }
+    if (second.kind === 'line')
+        return candidates(second, first, tolerance);
+    const delta = subtract(second.center, first.center);
+    const separation = Math.hypot(...delta);
+    if (separation <= tolerance &&
+        Math.abs(first.radius - second.radius) <= tolerance)
+        return [...endpoints(first), ...endpoints(second)];
+    if (!separation ||
+        separation > first.radius + second.radius + tolerance ||
+        separation < Math.abs(first.radius - second.radius) - tolerance)
+        return [];
+    // Normalize before squaring; the factored radius difference also avoids
+    // subtracting nearly equal squared radii.
+    const scale = Math.max(separation, first.radius, second.radius);
+    const d = separation / scale, r = first.radius / scale, s = second.radius / scale;
+    const along = (d * d + (r - s) * (r + s)) / (2 * d);
+    const tangent = Math.abs(separation - first.radius - second.radius) <= tolerance ||
+        Math.abs(separation - Math.abs(first.radius - second.radius)) <= tolerance;
+    const height = tangent
+        ? 0
+        : Math.sqrt(Math.max(0, (r - along) * (r + along))) * scale;
+    const unit = [delta[0] / separation, delta[1] / separation];
+    const foot = [
+        first.center[0] + along * scale * unit[0],
+        first.center[1] + along * scale * unit[1],
+    ];
+    return [-height, height].map(offset => [
+        foot[0] - offset * unit[1],
+        foot[1] + offset * unit[0],
+    ]);
+}
+//# sourceMappingURL=sketch-curve-intersections.js.map

package/bld/library/sketch-curve-intersections.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"sketch-curve-intersections.js","sourceRoot":"","sources":["../../src/library/sketch-curve-intersections.ts"],"names":[],"mappings":"AACA,OAAO,EACL,2BAA2B,EAC3B,mBAAmB,EACnB,oBAAoB,GAErB,MAAM,oBAAoB,CAAC;AAO5B,MAAM,QAAQ,GAAG,CAAC,CAAiB,EAAE,CAAiB,EAAkB,EAAE,CAAC;IACzE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACX,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;CACZ,CAAC;AACF,MAAM,KAAK,GAAG,CAAC,CAAiB,EAAE,CAAiB,EAAE,EAAE,CACrD,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;AAC5B,MAAM,QAAQ,GAAG,CAAC,CAAiB,EAAE,CAAiB,EAAE,EAAE,CACxD,IAAI,CAAC,KAAK,CAAC,GAAG,QAAQ,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;AAChC,MAAM,SAAS,GAAG,CAAC,KAAkB,EAA6B,EAAE,CAClE,KAAK,CAAC,IAAI,KAAK,QAAQ;IACrB,CAAC,CAAC,EAAE;IACJ,CAAC,CAAC,CAAC,mBAAmB,CAAC,KAAK,EAAE,CAAC,CAAC,EAAE,mBAAmB,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC;AAErE;;;;GAIG;AACH,MAAM,UAAU,wBAAwB,CACtC,KAAkB,EAClB,MAAmB;IAEnB,MAAM,UAAU,GAAG;QACjB,oBAAoB,CAAC,KAAK,CAAC;QAC3B,oBAAoB,CAAC,MAAM,CAAC;KAC7B,CAAC;IACF,MAAM,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,UAAU,CAAC,CAAC;IAC1C,MAAM,QAAQ,GAA8B,EAAE,CAAC;IAC/C,KAAK,MAAM,QAAQ,IAAI,UAAU,CAAC,KAAK,EAAE,MAAM,EAAE,SAAS,CAAC,EAAE,CAAC;QAC5D,MAAM,UAAU,GAAG;YACjB,2BAA2B,CAAC,KAAK,EAAE,QAAQ,CAAC;YAC5C,2BAA2B,CAAC,MAAM,EAAE,QAAQ,CAAC;SACrC,CAAC;QACX,IACE,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC,IAAI,CAClB,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE,CACX,QAAQ,CAAC,QAAQ,EAAE,mBAAmB,CAAC,KAAK,EAAE,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC;YAC7D,UAAU,CAAC,CAAC,CAAC,CAChB;YACD,QAAQ,CAAC,IAAI,CACX,OAAO,CAAC,EAAE,CAAC,QAAQ,CAAC,QAAQ,EAAE,OAAO,CAAC,QAAQ,CAAC,IAAI,SAAS,CAC7D;YAED,SAAS;QACX,QAAQ,CAAC,IAAI,CAAC,EAAC,QAAQ,EAAE,UAAU,EAAC,CAAC,CAAC;IACxC,CAAC;IACD,OAAO,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,UAAU,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC;AACpE,CAAC;AAED,SAAS,UAAU,CACjB,KAAkB,EAClB,MAAmB,EACnB,SAAiB;IAEjB,IAAI,KAAK,CAAC,IAAI,KAAK,MAAM,EAAE,CAAC;QAC1B,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC;QAC5B,MAAM,MAAM,GAAG,QAAQ,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAC9B,IAAI,CAAC,MAAM;YAAE,OAAO,CAAC,CAAC,CAAC,CAAC;QACxB,MAAM,IAAI,GAAmB;YAC3B,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM;YACtB,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM;SACvB,CAAC;QACF,IAAI,MAAM,CAAC,IAAI,KAAK,MAAM,EAAE,CAAC;YAC3B,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;YAC7B,MAAM,WAAW,GAAG,QAAQ,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YACnC,IAAI,CAAC,WAAW;gBAAE,OAAO,CAAC,CAAC,CAAC,CAAC;YAC7B,qEAAqE;YACrE,uEAAuE;YACvE,IACE,MAAM,CAAC,MAAM,CAAC,KAAK,CACjB,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI,SAAS,CACxD;gBAED,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YACtB,MAAM,SAAS,GAAmB;gBAChC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,WAAW;gBAC3B,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,WAAW;aAC5B,CAAC;YACF,MAAM,WAAW,GAAG,KAAK,CAAC,IAAI,EAAE,SAAS,CAAC,CAAC;YAC3C,IAAI,IAAI,CAAC,GAAG,CAAC,WAAW,CAAC,IAAI,MAAM,CAAC,OAAO,GAAG,EAAE;gBAAE,OAAO,EAAE,CAAC;YAC5D,MAAM,KAAK,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,SAAS,CAAC,GAAG,WAAW,CAAC;YAC7D,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC5D,CAAC;QACD,MAAM,aAAa,GAAG,KAAK,CAAC,IAAI,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC,CAAC;QAC9D,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,aAAa,CAAC,CAAC;QACpC,IAAI,GAAG,GAAG,MAAM,CAAC,MAAM,GAAG,SAAS;YAAE,OAAO,EAAE,CAAC;QAC/C,MAAM,UAAU,GAAG,GAAG,GAAG,MAAM,CAAC,MAAM,CAAC;QACvC,uEAAuE;QACvE,0EAA0E;QAC1E,MAAM,MAAM,GACV,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,MAAM,CAAC,MAAM,CAAC,IAAI,SAAS;YACxC,CAAC,CAAC,CAAC;YACH,CAAC,CAAC,MAAM,CAAC,MAAM;gBACb,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,UAAU,CAAC,GAAG,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;QAClE,2EAA2E;QAC3E,MAAM,IAAI,GAAmB;YAC3B,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,aAAa,GAAG,IAAI,CAAC,CAAC,CAAC;YAC1C,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,aAAa,GAAG,IAAI,CAAC,CAAC,CAAC;SAC3C,CAAC;QACF,OAAO,CAAC,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,CAAC;YACpC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,IAAI,CAAC,CAAC,CAAC;YACzB,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,IAAI,CAAC,CAAC,CAAC;SAC1B,CAAC,CAAC;IACL,CAAC;IACD,IAAI,MAAM,CAAC,IAAI,KAAK,MAAM;QAAE,OAAO,UAAU,CAAC,MAAM,EAAE,KAAK,EAAE,SAAS,CAAC,CAAC;IACxE,MAAM,KAAK,GAAG,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,KAAK,CAAC,MAAM,CAAC,CAAC;IACpD,MAAM,UAAU,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,KAAK,CAAC,CAAC;IACxC,IACE,UAAU,IAAI,SAAS;QACvB,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,MAAM,GAAG,MAAM,CAAC,MAAM,CAAC,IAAI,SAAS;QAEnD,OAAO,CAAC,GAAG,SAAS,CAAC,KAAK,CAAC,EAAE,GAAG,SAAS,CAAC,MAAM,CAAC,CAAC,CAAC;IACrD,IACE,CAAC,UAAU;QACX,UAAU,GAAG,KAAK,CAAC,MAAM,GAAG,MAAM,CAAC,MAAM,GAAG,SAAS;QACrD,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,MAAM,GAAG,MAAM,CAAC,MAAM,CAAC,GAAG,SAAS;QAE/D,OAAO,EAAE,CAAC;IACZ,wEAAwE;IACxE,0CAA0C;IAC1C,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,UAAU,EAAE,KAAK,CAAC,MAAM,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC;IAChE,MAAM,CAAC,GAAG,UAAU,GAAG,KAAK,EAC1B,CAAC,GAAG,KAAK,CAAC,MAAM,GAAG,KAAK,EACxB,CAAC,GAAG,MAAM,CAAC,MAAM,GAAG,KAAK,CAAC;IAC5B,MAAM,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACpD,MAAM,OAAO,GACX,IAAI,CAAC,GAAG,CAAC,UAAU,GAAG,KAAK,CAAC,MAAM,GAAG,MAAM,CAAC,MAAM,CAAC,IAAI,SAAS;QAChE,IAAI,CAAC,GAAG,CAAC,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,MAAM,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,IAAI,SAAS,CAAC;IAC7E,MAAM,MAAM,GAAG,OAAO;QACpB,CAAC,CAAC,CAAC;QACH,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;IAC9D,MAAM,IAAI,GAAmB,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,UAAU,EAAE,KAAK,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC;IAC5E,MAAM,IAAI,GAAmB;QAC3B,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,KAAK,GAAG,IAAI,CAAC,CAAC,CAAC;QACzC,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,KAAK,GAAG,IAAI,CAAC,CAAC,CAAC;KAC1C,CAAC;IACF,OAAO,CAAC,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,EAAE,CAAC;QACrC,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC;QAC1B,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC;KAC3B,CAAC,CAAC;AACL,CAAC"}

package/bld/library/sketch-curves.d.ts

@@ -0,0 +1,31 @@
+import type { SketchArcDirection, SketchEntitySnapshot, SketchPointAddress, SketchPosition } from './sketch.js';
+/** Analytic geometry only. Parameters are evaluation-local, never entity identity. */
+export type SketchCurve = Readonly<{
+    kind: 'line';
+    points: readonly [SketchPosition, SketchPosition];
+}> | Readonly<{
+    kind: 'circle';
+    center: SketchPosition;
+    radius: number;
+}> | Readonly<{
+    kind: 'arc';
+    center: SketchPosition;
+    radius: number;
+    start: number;
+    sweep: number;
+}>;
+export declare const sketchPositiveAngle: (angle: number) => number;
+export declare function sketchArcGeometry(center: SketchPosition, start: SketchPosition, end: SketchPosition, direction: SketchArcDirection): Extract<SketchCurve, {
+    kind: 'arc';
+}>;
+export declare function sketchCurveGeometry<E extends SketchEntitySnapshot>(entity: E, point: (address: SketchPointAddress) => SketchPosition): E extends {
+    kind: 'point';
+} ? undefined : Extract<SketchCurve, {
+    kind: E['kind'];
+}>;
+export declare function sketchCurvePosition(curve: SketchCurve, t: number): SketchPosition;
+/** Nearest point on the finite curve, including arc endpoints (not its full circle). */
+export declare function sketchCurveClosestParameter(curve: SketchCurve, position: SketchPosition): number;
+export declare function sketchCurveBounds(curve: SketchCurve): readonly SketchPosition[];
+export declare function sketchCurveTolerance(curve: SketchCurve): number;
+//# sourceMappingURL=sketch-curves.d.ts.map

package/bld/library/sketch-curves.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"sketch-curves.d.ts","sourceRoot":"","sources":["../../src/library/sketch-curves.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EACV,kBAAkB,EAClB,oBAAoB,EACpB,kBAAkB,EAClB,cAAc,EACf,MAAM,aAAa,CAAC;AAErB,sFAAsF;AACtF,MAAM,MAAM,WAAW,GACnB,QAAQ,CAAC;IAAC,IAAI,EAAE,MAAM,CAAC;IAAC,MAAM,EAAE,SAAS,CAAC,cAAc,EAAE,cAAc,CAAC,CAAA;CAAC,CAAC,GAC3E,QAAQ,CAAC;IAAC,IAAI,EAAE,QAAQ,CAAC;IAAC,MAAM,EAAE,cAAc,CAAC;IAAC,MAAM,EAAE,MAAM,CAAA;CAAC,CAAC,GAClE,QAAQ,CAAC;IACP,IAAI,EAAE,KAAK,CAAC;IACZ,MAAM,EAAE,cAAc,CAAC;IACvB,MAAM,EAAE,MAAM,CAAC;IACf,KAAK,EAAE,MAAM,CAAC;IACd,KAAK,EAAE,MAAM,CAAC;CACf,CAAC,CAAC;AAGP,eAAO,MAAM,mBAAmB,UAAW,MAAM,WACpB,CAAC;AAE9B,wBAAgB,iBAAiB,CAC/B,MAAM,EAAE,cAAc,EACtB,KAAK,EAAE,cAAc,EACrB,GAAG,EAAE,cAAc,EACnB,SAAS,EAAE,kBAAkB,GAC5B,OAAO,CAAC,WAAW,EAAE;IAAC,IAAI,EAAE,KAAK,CAAA;CAAC,CAAC,CAWrC;AAED,wBAAgB,mBAAmB,CAAC,CAAC,SAAS,oBAAoB,EAChE,MAAM,EAAE,CAAC,EACT,KAAK,EAAE,CAAC,OAAO,EAAE,kBAAkB,KAAK,cAAc,GACrD,CAAC,SAAS;IAAC,IAAI,EAAE,OAAO,CAAA;CAAC,GACxB,SAAS,GACT,OAAO,CAAC,WAAW,EAAE;IAAC,IAAI,EAAE,CAAC,CAAC,MAAM,CAAC,CAAA;CAAC,CAAC,CAAC;AAgC5C,wBAAgB,mBAAmB,CACjC,KAAK,EAAE,WAAW,EAClB,CAAC,EAAE,MAAM,GACR,cAAc,CAWhB;AAED,wFAAwF;AACxF,wBAAgB,2BAA2B,CACzC,KAAK,EAAE,WAAW,EAClB,QAAQ,EAAE,cAAc,GACvB,MAAM,CAwBR;AAED,wBAAgB,iBAAiB,CAC/B,KAAK,EAAE,WAAW,GACjB,SAAS,cAAc,EAAE,CAiB3B;AAED,wBAAgB,oBAAoB,CAAC,KAAK,EAAE,WAAW,GAAG,MAAM,CAc/D"}

package/bld/library/sketch-curves.js

@@ -0,0 +1,89 @@
+import { sketchRelativePrecision } from './sketch-precision.js';
+const tau = Math.PI * 2;
+export const sketchPositiveAngle = (angle) => ((angle % tau) + tau) % tau;
+export function sketchArcGeometry(center, start, end, direction) {
+    const angle = Math.atan2(start[1] - center[1], start[0] - center[0]);
+    const endAngle = Math.atan2(end[1] - center[1], end[0] - center[0]);
+    const sign = direction === 'ccw' ? 1 : -1;
+    return {
+        kind: 'arc',
+        center,
+        radius: Math.hypot(start[0] - center[0], start[1] - center[1]),
+        start: angle,
+        sweep: sign * sketchPositiveAngle(sign * (endAngle - angle)),
+    };
+}
+export function sketchCurveGeometry(entity, point) {
+    switch (entity.kind) {
+        case 'point':
+            return undefined;
+        case 'line':
+            return {
+                kind: 'line',
+                points: [point(entity.points[0]), point(entity.points[1])],
+            };
+        case 'circle':
+            return {
+                kind: 'circle',
+                center: point(entity.center),
+                radius: entity.radius,
+            };
+        case 'arc':
+            return {
+                ...sketchArcGeometry(point(entity.center), point(entity.points[0]), point(entity.points[1]), entity.direction),
+                radius: entity.radius,
+            };
+    }
+}
+export function sketchCurvePosition(curve, t) {
+    if (curve.kind === 'line') {
+        const [a, b] = curve.points;
+        return [a[0] + t * (b[0] - a[0]), a[1] + t * (b[1] - a[1])];
+    }
+    const angle = curve.kind === 'circle' ? t * tau : curve.start + t * curve.sweep;
+    return [
+        curve.center[0] + curve.radius * Math.cos(angle),
+        curve.center[1] + curve.radius * Math.sin(angle),
+    ];
+}
+/** Nearest point on the finite curve, including arc endpoints (not its full circle). */
+export function sketchCurveClosestParameter(curve, position) {
+    if (curve.kind === 'line') {
+        const [a, b] = curve.points, dx = b[0] - a[0], dy = b[1] - a[1];
+        return Math.max(0, Math.min(1, ((position[0] - a[0]) * dx + (position[1] - a[1]) * dy) /
+            (dx * dx + dy * dy) || 0));
+    }
+    const angle = Math.atan2(position[1] - curve.center[1], position[0] - curve.center[0]);
+    if (curve.kind === 'circle')
+        return sketchPositiveAngle(angle) / tau;
+    const sign = Math.sign(curve.sweep), sweep = Math.abs(curve.sweep);
+    const distance = sketchPositiveAngle(sign * (angle - curve.start));
+    if (distance <= sweep)
+        return sweep ? distance / sweep : 0;
+    return distance - sweep < tau - distance ? 1 : 0;
+}
+export function sketchCurveBounds(curve) {
+    if (curve.kind === 'line')
+        return curve.points;
+    const candidates = [0, Math.PI / 2, Math.PI, (3 * Math.PI) / 2].filter(angle => curve.kind === 'circle' ||
+        sketchPositiveAngle(Math.sign(curve.sweep) * (angle - curve.start)) <=
+            Math.abs(curve.sweep));
+    return [
+        ...(curve.kind === 'arc'
+            ? [sketchCurvePosition(curve, 0), sketchCurvePosition(curve, 1)]
+            : []),
+        ...candidates.map((angle) => [
+            curve.center[0] + curve.radius * Math.cos(angle),
+            curve.center[1] + curve.radius * Math.sin(angle),
+        ]),
+    ];
+}
+export function sketchCurveTolerance(curve) {
+    const size = curve.kind === 'line'
+        ? Math.hypot(curve.points[1][0] - curve.points[0][0], curve.points[1][1] - curve.points[0][1])
+        : curve.radius;
+    return Math.max(size * sketchRelativePrecision, Number.EPSILON *
+        16 *
+        Math.max(...sketchCurveBounds(curve).flatMap(p => p.map(Math.abs))));
+}
+//# sourceMappingURL=sketch-curves.js.map

package/bld/library/sketch-curves.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"sketch-curves.js","sourceRoot":"","sources":["../../src/library/sketch-curves.ts"],"names":[],"mappings":"AAAA,OAAO,EAAC,uBAAuB,EAAC,MAAM,uBAAuB,CAAC;AAoB9D,MAAM,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;AACxB,MAAM,CAAC,MAAM,mBAAmB,GAAG,CAAC,KAAa,EAAE,EAAE,CACnD,CAAC,CAAC,KAAK,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;AAE9B,MAAM,UAAU,iBAAiB,CAC/B,MAAsB,EACtB,KAAqB,EACrB,GAAmB,EACnB,SAA6B;IAE7B,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;IACrE,MAAM,QAAQ,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;IACpE,MAAM,IAAI,GAAG,SAAS,KAAK,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1C,OAAO;QACL,IAAI,EAAE,KAAK;QACX,MAAM;QACN,MAAM,EAAE,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QAC9D,KAAK,EAAE,KAAK;QACZ,KAAK,EAAE,IAAI,GAAG,mBAAmB,CAAC,IAAI,GAAG,CAAC,QAAQ,GAAG,KAAK,CAAC,CAAC;KAC7D,CAAC;AACJ,CAAC;AAQD,MAAM,UAAU,mBAAmB,CACjC,MAA4B,EAC5B,KAAsD;IAEtD,QAAQ,MAAM,CAAC,IAAI,EAAE,CAAC;QACpB,KAAK,OAAO;YACV,OAAO,SAAS,CAAC;QACnB,KAAK,MAAM;YACT,OAAO;gBACL,IAAI,EAAE,MAAM;gBACZ,MAAM,EAAE,CAAC,KAAK,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;aAC3D,CAAC;QACJ,KAAK,QAAQ;YACX,OAAO;gBACL,IAAI,EAAE,QAAQ;gBACd,MAAM,EAAE,KAAK,CAAC,MAAM,CAAC,MAAM,CAAC;gBAC5B,MAAM,EAAE,MAAM,CAAC,MAAM;aACtB,CAAC;QACJ,KAAK,KAAK;YACR,OAAO;gBACL,GAAG,iBAAiB,CAClB,KAAK,CAAC,MAAM,CAAC,MAAM,CAAC,EACpB,KAAK,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,EACvB,KAAK,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,EACvB,MAAM,CAAC,SAAS,CACjB;gBACD,MAAM,EAAE,MAAM,CAAC,MAAM;aACtB,CAAC;IACN,CAAC;AACH,CAAC;AAED,MAAM,UAAU,mBAAmB,CACjC,KAAkB,EAClB,CAAS;IAET,IAAI,KAAK,CAAC,IAAI,KAAK,MAAM,EAAE,CAAC;QAC1B,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC;QAC5B,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC9D,CAAC;IACD,MAAM,KAAK,GACT,KAAK,CAAC,IAAI,KAAK,QAAQ,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,KAAK,GAAG,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC;IACpE,OAAO;QACL,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC;QAChD,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC;KACjD,CAAC;AACJ,CAAC;AAED,wFAAwF;AACxF,MAAM,UAAU,2BAA2B,CACzC,KAAkB,EAClB,QAAwB;IAExB,IAAI,KAAK,CAAC,IAAI,KAAK,MAAM,EAAE,CAAC;QAC1B,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,EACzB,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAChB,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACnB,OAAO,IAAI,CAAC,GAAG,CACb,CAAC,EACD,IAAI,CAAC,GAAG,CACN,CAAC,EACD,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;YACrD,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,IAAI,CAAC,CAC3B,CACF,CAAC;IACJ,CAAC;IACD,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CACtB,QAAQ,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,EAC7B,QAAQ,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,CAC9B,CAAC;IACF,IAAI,KAAK,CAAC,IAAI,KAAK,QAAQ;QAAE,OAAO,mBAAmB,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC;IACrE,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,EACjC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC;IAChC,MAAM,QAAQ,GAAG,mBAAmB,CAAC,IAAI,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC;IACnE,IAAI,QAAQ,IAAI,KAAK;QAAE,OAAO,KAAK,CAAC,CAAC,CAAC,QAAQ,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IAC3D,OAAO,QAAQ,GAAG,KAAK,GAAG,GAAG,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AACnD,CAAC;AAED,MAAM,UAAU,iBAAiB,CAC/B,KAAkB;IAElB,IAAI,KAAK,CAAC,IAAI,KAAK,MAAM;QAAE,OAAO,KAAK,CAAC,MAAM,CAAC;IAC/C,MAAM,UAAU,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,GAAG,CAAC,EAAE,IAAI,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,CACpE,KAAK,CAAC,EAAE,CACN,KAAK,CAAC,IAAI,KAAK,QAAQ;QACvB,mBAAmB,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC;YACjE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,KAAK,CAAC,CAC1B,CAAC;IACF,OAAO;QACL,GAAG,CAAC,KAAK,CAAC,IAAI,KAAK,KAAK;YACtB,CAAC,CAAC,CAAC,mBAAmB,CAAC,KAAK,EAAE,CAAC,CAAC,EAAE,mBAAmB,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;YAChE,CAAC,CAAC,EAAE,CAAC;QACP,GAAG,UAAU,CAAC,GAAG,CAAC,CAAC,KAAK,EAAkB,EAAE,CAAC;YAC3C,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC;YAChD,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC;SACjD,CAAC;KACH,CAAC;AACJ,CAAC;AAED,MAAM,UAAU,oBAAoB,CAAC,KAAkB;IACrD,MAAM,IAAI,GACR,KAAK,CAAC,IAAI,KAAK,MAAM;QACnB,CAAC,CAAC,IAAI,CAAC,KAAK,CACR,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EACvC,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CACxC;QACH,CAAC,CAAC,KAAK,CAAC,MAAM,CAAC;IACnB,OAAO,IAAI,CAAC,GAAG,CACb,IAAI,GAAG,uBAAuB,EAC9B,MAAM,CAAC,OAAO;QACZ,EAAE;QACF,IAAI,CAAC,GAAG,CAAC,GAAG,iBAAiB,CAAC,KAAK,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CACtE,CAAC;AACJ,CAAC"}

package/bld/library/sketch-drag-rules.d.ts

@@ -0,0 +1,35 @@
+import { type SketchSolveObjective, type SketchSolveProblem, type SketchSolveTarget, type SketchSolveResult } from './sketch-solver.js';
+import type { SketchPosition } from './sketch.js';
+import { type SketchIncidence, type SketchIncidenceGeometry } from './sketch-incidence.js';
+export type SketchDragContext = Readonly<{
+    reference: SketchSolveProblem;
+    target: SketchSolveTarget;
+    /** Displayed gesture-start geometry, before applying AST coordinate locks. */
+    geometry?: SketchIncidenceGeometry;
+}>;
+export type SketchDragPlan = Readonly<{
+    problem: SketchSolveProblem;
+    stages: readonly [SketchDragStage, ...SketchDragStage[]];
+    incidences?: readonly SketchIncidence[];
+}>;
+/** A rule may prepare a stage's seed from the preceding feasible geometry.
+ * The seed is only an initial guess; all accumulated locks still apply. */
+export type SketchDragStage = Readonly<{
+    objectives: (current: SketchSolveProblem) => readonly SketchSolveObjective[];
+    seed?: (current: SketchSolveProblem, objectives: readonly SketchSolveObjective[]) => SketchSolveProblem;
+}>;
+export type SketchDragSession = (current: SketchSolveProblem, target: SketchSolveTarget) => SketchDragPlan;
+export type SketchDragRule = (context: SketchDragContext) => SketchDragSession | undefined;
+/** Dispatch knows neither point roles nor connected components. */
+export declare function createSketchDragSession(context: SketchDragContext, rules?: readonly SketchDragRule[]): SketchDragSession;
+export declare function solveSketchDrag(current: SketchSolveProblem, target: SketchSolveTarget, reference: SketchSolveProblem, geometry?: SketchDragContext['geometry']): Readonly<{
+    positions: readonly SketchPosition[];
+    radii: readonly number[];
+    arcRadii: readonly number[];
+    degreesOfFreedom: number;
+    redundant: readonly number[];
+}>;
+/** Each stage keeps the preceding stage's chosen parameter values, not all
+ * equivalent optima. These locks belong to this solve, never the next frame. */
+export declare function solveSketchDragPlan(plan: SketchDragPlan): SketchSolveResult;
+//# sourceMappingURL=sketch-drag-rules.d.ts.map

package/bld/library/sketch-drag-rules.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"sketch-drag-rules.d.ts","sourceRoot":"","sources":["../../src/library/sketch-drag-rules.ts"],"names":[],"mappings":"AAAA,OAAO,EAEL,KAAK,oBAAoB,EACzB,KAAK,kBAAkB,EACvB,KAAK,iBAAiB,EAEtB,KAAK,iBAAiB,EACvB,MAAM,oBAAoB,CAAC;AAC5B,OAAO,KAAK,EAAC,cAAc,EAAC,MAAM,aAAa,CAAC;AAChD,OAAO,EAML,KAAK,eAAe,EACpB,KAAK,uBAAuB,EAC7B,MAAM,uBAAuB,CAAC;AAM/B,MAAM,MAAM,iBAAiB,GAAG,QAAQ,CAAC;IACvC,SAAS,EAAE,kBAAkB,CAAC;IAC9B,MAAM,EAAE,iBAAiB,CAAC;IAC1B,8EAA8E;IAC9E,QAAQ,CAAC,EAAE,uBAAuB,CAAC;CACpC,CAAC,CAAC;AACH,MAAM,MAAM,cAAc,GAAG,QAAQ,CAAC;IACpC,OAAO,EAAE,kBAAkB,CAAC;IAC5B,MAAM,EAAE,SAAS,CAAC,eAAe,EAAE,GAAG,eAAe,EAAE,CAAC,CAAC;IACzD,UAAU,CAAC,EAAE,SAAS,eAAe,EAAE,CAAC;CACzC,CAAC,CAAC;AACH;2EAC2E;AAC3E,MAAM,MAAM,eAAe,GAAG,QAAQ,CAAC;IACrC,UAAU,EAAE,CAAC,OAAO,EAAE,kBAAkB,KAAK,SAAS,oBAAoB,EAAE,CAAC;IAC7E,IAAI,CAAC,EAAE,CACL,OAAO,EAAE,kBAAkB,EAC3B,UAAU,EAAE,SAAS,oBAAoB,EAAE,KACxC,kBAAkB,CAAC;CACzB,CAAC,CAAC;AACH,MAAM,MAAM,iBAAiB,GAAG,CAC9B,OAAO,EAAE,kBAAkB,EAC3B,MAAM,EAAE,iBAAiB,KACtB,cAAc,CAAC;AACpB,MAAM,MAAM,cAAc,GAAG,CAC3B,OAAO,EAAE,iBAAiB,KACvB,iBAAiB,GAAG,SAAS,CAAC;AAEnC,mEAAmE;AACnE,wBAAgB,uBAAuB,CACrC,OAAO,EAAE,iBAAiB,EAC1B,KAAK,GAAE,SAAS,cAAc,EAAc,GAC3C,iBAAiB,CAMnB;AAED,wBAAgB,eAAe,CAC7B,OAAO,EAAE,kBAAkB,EAC3B,MAAM,EAAE,iBAAiB,EACzB,SAAS,EAAE,kBAAkB,EAC7B,QAAQ,CAAC,EAAE,iBAAiB,CAAC,UAAU,CAAC;;;;;;GA4CzC;AAED;gFACgF;AAChF,wBAAgB,mBAAmB,CAAC,IAAI,EAAE,cAAc,GAAG,iBAAiB,CAwC3E"}