flo-mat 4.1.5 → 5.0.1

package/node/reconstruct-from-mat/reconstruct-from-mat.js

@@ -1,9 +1,10 @@
-import { distanceBetween, fromTo as fromToVec, dot, rotate90Degrees, cross, interpolate, translate, lineLineIntersection } from 'flo-vector2d';
-import { tangent as getTangent } from 'flo-bezier3';
+import { distanceBetween } from 'flo-vector2d';
 import { mapDeep } from '../utils/map-deep.js';
 import { removeSatInducedTrivial3Prongs } from './remove-sat-induced-trivial-3prongs.js';
 import { getCpNodessForSat } from './get-cp-nodess-for-sat.js';
-import { getRealProngCount } from '../cp-node/fs/get-real-prong-count.js';
+import { clone } from '../cp-node/fs/clone.js';
+import { getBoundaryTangentLine } from '../point-on-shape/get-boundary-tangent-line.js';
+import { getBeziersFromTangents } from '../utils/get-beziers-from-tangents.js';
 const { max, abs, round, log2 } = Math;
 // let ii = 0;
 /**
@@ -14,18 +15,25 @@ const { max, abs, round, log2 } = Math;
  * @internal
  */
 function reconstructFromMats(mats, minNumAccurateBits = undefined) {
-    const maxPO2 = max(...mats.map(mat => mat.meta.maxCoordPowerOf2));
+    let maxPO2 = -Infinity;
+    for (const mat of mats) {
+        maxPO2 = max(maxPO2, mat.meta.maxCoordPowerOf2);
+    }
     const recoverBoundary_ = recoverBoundary(maxPO2, minNumAccurateBits);
+    // `removeSatInducedTrivial3Prongs` relinks nodes in place so work on 
+    // clones - the given `Mat`s must not be modified. (Only 
+    // `meta.maxCoordPowerOf2` is used from the metas so they can be shared.)
     // Really only required if getSats was called on the `Mat`s
     const cpNodessPerSat = mats
+        .map(mat => ({ cpNode: clone(mat.cpNode), meta: mat.meta }))
         .map(removeSatInducedTrivial3Prongs)
         .map(getCpNodessForSat);
     const r = mapDeep(3, cpNodessPerSat, recoverBoundary_, true);
     // ii;//?
-    return r;
+    // Flatten the (1 or 2) beziers per boundary piece into a single array of
+    // beziers per loop.
+    return r.map(loops => loops.map(pieces => pieces.flat()));
 }
-const FROM = 0;
-const TO = 1;
 /**
  *
  * @param cpNode
@@ -48,8 +56,6 @@ function recoverBoundary(maxPO2, minNumAccurateBits) {
             return undefined;
         }
         const d = distanceBetween(circle.center, p);
-        if (getRealProngCount(cpNode)) {
-        }
         // DON'T remove these - they tell us the accuracy of the algorithm
         // if (!closeTo(2**34)(circle.radius, d) && getRealProngCount(cpNode) === 3) {
         //     ii++;
@@ -80,81 +86,15 @@ function recoverBoundary(maxPO2, minNumAccurateBits) {
         return getBoundaryPiece(cpNode);
     };
 }
-/** @internal */
-const TOLERANCE_ADD_2PRONG = 0.01;
-/** @internal */
-const TOLERANCE_USE_LINE = 0.0001; // else cubic
-function getBoundaryPiece(cpNode) {
-    const cpNode_ = cpNode.next;
-    const pos = cpNode.pointOnShape;
-    const pos_ = cpNode_.pointOnShape;
-    const { p } = pos;
-    const { p: p_ } = pos_;
-    const fromCc = p;
-    const fromL = getEdgeDirection(cpNode, FROM);
-    const toCc = p_;
-    const toL = getEdgeDirection(cpNode_, TO);
-    const mid = lineLineIntersection(fromL, toL);
-    const c = fromToVec(fromCc, toCc);
-    let twisted;
-    if (!mid) {
-        twisted = true;
-    }
-    else {
-        const a = fromToVec(fromCc, mid);
-        const b = fromToVec(toCc, mid);
-        twisted = dot(a, c) < 0 || dot(b, c) > 0;
-    }
-    if (!twisted) {
-        return [fromCc, mid, toCc];
-    }
-    const r = rotate90Degrees(c);
-    const w1 = fromToVec(fromL[0], fromL[1]); // This is a unit vector
-    const w2 = fromToVec(toL[0], toL[1]); // This is a unit vector
-    const d1 = abs(cross(c, w1)) / (3 * 3);
-    const d2 = abs(cross(c, w2)) / (3 * 3);
-    if (d1 > TOLERANCE_ADD_2PRONG || d2 > TOLERANCE_ADD_2PRONG) {
-        // FUTURE - not within tolerance - must add additional 2-prong
-        return [fromCc, toCc];
-    }
-    if (d1 > TOLERANCE_USE_LINE || d2 > TOLERANCE_USE_LINE) {
-        // approximate with cubic bezier
-        const m1 = interpolate(fromCc, toCc, 1 / 3);
-        const m2 = interpolate(fromCc, toCc, 2 / 3);
-        const v1 = translate(r, m1);
-        const v2 = translate(r, m2);
-        const l1 = [m1, v1];
-        const l2 = [m2, v2];
-        const mid1 = lineLineIntersection(fromL, l1);
-        const mid2 = lineLineIntersection(toL, l2);
-        return [fromCc, mid1, mid2, toCc];
-    }
-    // Within tolerance - approximate with a straight line.
-    return [fromCc, toCc];
-}
 /**
- * Returns a line segment of unit length starting at the given point and
- * being tangent to the shape boundary.
+ * Returns the bezier curves (1 or 2 lines or quadratics) approximating the
+ * shape boundary from the contact point of the given `CpNode` to that of the
+ * next `CpNode`.
  *
  * @internal
  */
-function getEdgeDirection(cpNode, fromOrTo) {
-    let { pointOnShape } = cpNode;
-    let { curve, t, p } = pointOnShape;
-    let { ps } = curve;
-    if (fromOrTo === FROM && t === 1) {
-        t = 0;
-        curve = curve.next;
-        ps = curve.ps;
-    }
-    if (fromOrTo === TO && t === 0) {
-        t = 1;
-        curve = curve.prev;
-        ps = curve.ps;
-    }
-    const tangent = getTangent(ps, t);
-    const v = translate(p, tangent);
-    return [p, v];
+function getBoundaryPiece(cpNode) {
+    return getBeziersFromTangents(getBoundaryTangentLine(cpNode.pointOnShape, true), getBoundaryTangentLine(cpNode.next.pointOnShape, false));
 }
 export { reconstructFromMats };
 //# sourceMappingURL=reconstruct-from-mat.js.map

package/node/reconstruct-from-mat/reconstruct-from-mat.js.map

@@ -1 +1 @@
-{"version":3,"file":"reconstruct-from-mat.js","sourceRoot":"","sources":["../../src/reconstruct-from-mat/reconstruct-from-mat.ts"],"names":[],"mappings":"AAGA,OAAO,EACW,eAAe,EAAE,MAAM,IAAI,SAAS,EAAE,GAAG,EAAE,eAAe,EACxE,KAAK,EAAE,WAAW,EAAE,SAAS,EAAE,oBAAoB,EACtD,MAAM,cAAc,CAAC;AACtB,OAAO,EAAE,OAAO,IAAI,UAAU,EAAE,MAAM,aAAa,CAAC;AACpD,OAAO,EAAE,OAAO,EAAE,MAAM,sBAAsB,CAAC;AAC/C,OAAO,EAAE,8BAA8B,EAAE,MAAM,yCAAyC,CAAC;AACzF,OAAO,EAAE,iBAAiB,EAAE,MAAM,4BAA4B,CAAC;AAC/D,OAAO,EAAE,iBAAiB,EAAE,MAAM,uCAAuC,CAAC;AAG1E,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,KAAK,EAAE,IAAI,EAAE,GAAG,IAAI,CAAC;AAGvC,cAAc;AAGd;;;;;;GAMG;AACH,SAAS,mBAAmB,CACpB,IAAW,EACX,kBAAkB,GAAuB,SAAS;IAEtD,MAAM,MAAM,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,gBAAgB,CAAC,CAAC,CAAC;IAClE,MAAM,gBAAgB,GAAG,eAAe,CAAC,MAAM,EAAE,kBAAkB,CAAC,CAAC;IAErE,2DAA2D;IAC3D,MAAM,cAAc,GAAG,IAAI;SACtB,GAAG,CAAC,8BAA8B,CAAC;SACnC,GAAG,CAAC,iBAAiB,CAAC,CAAC;IAE5B,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,EAAE,cAAc,EAAE,gBAAgB,EAAE,IAAI,CAAC,CAAC;IAC7D,SAAS;IAET,OAAO,CAAC,CAAC;AACb,CAAC;AAGD,MAAM,IAAI,GAAG,CAAC,CAAC;AACf,MAAM,EAAE,GAAG,CAAC,CAAC;AAEb;;;;;GAKG;AACH,SAAS,eAAe,CAChB,MAAc,EACd,kBAAsC;IAE1C,OAAO,UAAU,MAAc;QAC3B,MAAM,OAAO,GAAG,MAAM,CAAC,IAAI,CAAC;QAE5B,MAAM,GAAG,GAAG,MAAM,CAAC,YAAY,CAAC;QAChC,MAAM,IAAI,GAAG,OAAO,CAAC,YAAY,CAAC;QAElC,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE,GAAG,GAAG,CAAC;QACjC,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,EAAE,GAAG,IAAI,CAAC;QAEtC,IAAI,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;YACnC,OAAO,SAAS,CAAC;QACrB,CAAC;QAED,IAAI,KAAK,CAAC,IAAI,KAAK,MAAM,CAAC,IAAI,EAAE,CAAC;YAC7B,mDAAmD;YACnD,oBAAoB;YACpB,OAAO,SAAS,CAAC;QACrB,CAAC;QAED,MAAM,CAAC,GAAG,eAAe,CAAC,MAAM,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;QAE5C,IAAI,iBAAiB,CAAC,MAAM,CAAC,EAAE,CAAC;QAEhC,CAAC;QAED,kEAAkE;QAClE,8EAA8E;QAC9E,YAAY;QACZ,wBAAwB;QACxB,wBAAwB;QACxB,wBAAwB;QACxB,YAAY;QACZ,uBAAuB;QACvB,oCAAoC;QACpC,IAAI;QACJ,8EAA8E;QAC9E,YAAY;QACZ,wBAAwB;QACxB,wBAAwB;QACxB,wBAAwB;QACxB,YAAY;QACZ,uBAAuB;QACvB,oCAAoC;QACpC,IAAI;QACJ,IAAI,kBAAkB,KAAK,SAAS,EAAE,CAAC;YACnC,MAAM,GAAG,GAAG,CAAC,IAAE,CAAC,MAAM,GAAG,kBAAkB,CAAC,CAAC;YAE7C,IAAI,GAAG,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC;gBAC/B,MAAM,IAAI,KAAK,CACX,mDAAmD,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI;oBACvE,2DAA2D;oBAC3D,wBAAwB,MAAM,CAAC,MAAM,eAAe,CAAC,EAAE,CAC1D,CAAC;YACN,CAAC;QACL,CAAC;QAED,OAAO,gBAAgB,CAAC,MAAM,CAAC,CAAC;IACpC,CAAC,CAAA;AACL,CAAC;AAGD,gBAAgB;AAChB,MAAM,oBAAoB,GAAG,IAAI,CAAC;AAClC,gBAAgB;AAChB,MAAM,kBAAkB,GAAG,MAAM,CAAC,CAAC,aAAa;AAGhD,SAAS,gBAAgB,CACjB,MAAc;IAElB,MAAM,OAAO,GAAG,MAAM,CAAC,IAAI,CAAC;IAE5B,MAAM,GAAG,GAAG,MAAM,CAAC,YAAY,CAAC;IAChC,MAAM,IAAI,GAAG,OAAO,CAAC,YAAY,CAAC;IAElC,MAAM,EAAE,CAAC,EAAE,GAAG,GAAG,CAAC;IAClB,MAAM,EAAE,CAAC,EAAE,EAAE,EAAE,GAAG,IAAI,CAAC;IAEvB,MAAM,MAAM,GAAG,CAAC,CAAC;IACjB,MAAM,KAAK,GAAG,gBAAgB,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC;IAE7C,MAAM,IAAI,GAAG,EAAE,CAAC;IAChB,MAAM,GAAG,GAAG,gBAAgB,CAAC,OAAO,EAAE,EAAE,CAAC,CAAC;IAE1C,MAAM,GAAG,GAAG,oBAAoB,CAAC,KAAK,EAAE,GAAG,CAAC,CAAC;IAE7C,MAAM,CAAC,GAAG,SAAS,CAAC,MAAM,EAAC,IAAI,CAAC,CAAC;IAEjC,IAAI,OAAgB,CAAC;IACrB,IAAI,CAAC,GAAG,EAAE,CAAC;QACP,OAAO,GAAG,IAAI,CAAC;IACnB,CAAC;SAAM,CAAC;QACJ,MAAM,CAAC,GAAG,SAAS,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC;QACjC,MAAM,CAAC,GAAG,SAAS,CAAC,IAAI,EAAI,GAAG,CAAC,CAAC;QAEjC,OAAO,GAAG,GAAG,CAAC,CAAC,EAAC,CAAC,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,EAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAC3C,CAAC;IAED,IAAI,CAAC,OAAO,EAAE,CAAC;QACX,OAAO,CAAC,MAAM,EAAE,GAAI,EAAE,IAAI,CAAC,CAAC;IAChC,CAAC;IAED,MAAM,CAAC,GAAG,eAAe,CAAC,CAAC,CAAC,CAAC;IAC7B,MAAM,EAAE,GAAG,SAAS,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,wBAAwB;IAClE,MAAM,EAAE,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC,CAAC,EAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAE,wBAAwB;IAEhE,MAAM,EAAE,GAAG,GAAG,CAAC,KAAK,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAC,CAAC,CAAC,CAAC;IACrC,MAAM,EAAE,GAAG,GAAG,CAAC,KAAK,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAC,CAAC,CAAC,CAAC;IAErC,IAAI,EAAE,GAAG,oBAAoB,IAAI,EAAE,GAAG,oBAAoB,EAAE,CAAC;QACzD,8DAA8D;QAC9D,OAAO,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC;IAC1B,CAAC;IAED,IAAI,EAAE,GAAG,kBAAkB,IAAI,EAAE,GAAG,kBAAkB,EAAE,CAAC;QACrD,gCAAgC;QAChC,MAAM,EAAE,GAAG,WAAW,CAAC,MAAM,EAAC,IAAI,EAAC,CAAC,GAAC,CAAC,CAAC,CAAC;QACxC,MAAM,EAAE,GAAG,WAAW,CAAC,MAAM,EAAC,IAAI,EAAC,CAAC,GAAC,CAAC,CAAC,CAAC;QACxC,MAAM,EAAE,GAAG,SAAS,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QAC5B,MAAM,EAAE,GAAG,SAAS,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QAC5B,MAAM,EAAE,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;QACnB,MAAM,EAAE,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;QACnB,MAAM,IAAI,GAAG,oBAAoB,CAAC,KAAK,EAAE,EAAE,CAAE,CAAC;QAC9C,MAAM,IAAI,GAAG,oBAAoB,CAAC,GAAG,EAAI,EAAE,CAAE,CAAC;QAE9C,OAAO,CAAC,MAAM,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;IACtC,CAAC;IAED,uDAAuD;IACvD,OAAO,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC;AAC1B,CAAC;AAGD;;;;;GAKG;AACH,SAAS,gBAAgB,CACjB,MAAc,EACd,QAAiC;IAErC,IAAI,EAAE,YAAY,EAAE,GAAG,MAAM,CAAC;IAC9B,IAAI,EAAE,KAAK,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,YAAY,CAAC;IACnC,IAAI,EAAE,EAAE,EAAE,GAAG,KAAK,CAAC;IAEnB,IAAI,QAAQ,KAAK,IAAI,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QAC/B,CAAC,GAAG,CAAC,CAAC;QACN,KAAK,GAAG,KAAK,CAAC,IAAI,CAAC;QACnB,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC;IAClB,CAAC;IACD,IAAI,QAAQ,KAAK,EAAE,IAAI,CAAC,KAAK,CAAC,EAAG,CAAC;QAC9B,CAAC,GAAG,CAAC,CAAC;QACN,KAAK,GAAG,KAAK,CAAC,IAAI,CAAC;QACnB,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC;IAClB,CAAC;IAED,MAAM,OAAO,GAAG,UAAU,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;IAElC,MAAM,CAAC,GAAG,SAAS,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;IAEhC,OAAO,CAAC,CAAC,EAAC,CAAC,CAAC,CAAC;AACjB,CAAC;AAGD,OAAO,EAAE,mBAAmB,EAAE,CAAA"}
+{"version":3,"file":"reconstruct-from-mat.js","sourceRoot":"","sources":["../../src/reconstruct-from-mat/reconstruct-from-mat.ts"],"names":[],"mappings":"AAGA,OAAO,EAAE,eAAe,EAAE,MAAM,cAAc,CAAC;AAC/C,OAAO,EAAE,OAAO,EAAE,MAAM,sBAAsB,CAAC;AAC/C,OAAO,EAAE,8BAA8B,EAAE,MAAM,yCAAyC,CAAC;AACzF,OAAO,EAAE,iBAAiB,EAAE,MAAM,4BAA4B,CAAC;AAC/D,OAAO,EAAE,KAAK,EAAE,MAAM,wBAAwB,CAAC;AAC/C,OAAO,EAAE,sBAAsB,EAAE,MAAM,gDAAgD,CAAC;AACxF,OAAO,EAAE,sBAAsB,EAAE,MAAM,uCAAuC,CAAC;AAG/E,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,KAAK,EAAE,IAAI,EAAE,GAAG,IAAI,CAAC;AAGvC,cAAc;AAGd;;;;;;GAMG;AACH,SAAS,mBAAmB,CACpB,IAAW,EACX,kBAAkB,GAAuB,SAAS;IAEtD,IAAI,MAAM,GAAG,CAAC,QAAQ,CAAC;IACvB,KAAK,MAAM,GAAG,IAAI,IAAI,EAAE,CAAC;QAAC,MAAM,GAAG,GAAG,CAAC,MAAM,EAAE,GAAG,CAAC,IAAI,CAAC,gBAAgB,CAAC,CAAC;IAAC,CAAC;IAC5E,MAAM,gBAAgB,GAAG,eAAe,CAAC,MAAM,EAAE,kBAAkB,CAAC,CAAC;IAErE,sEAAsE;IACtE,yDAAyD;IACzD,yEAAyE;IACzE,2DAA2D;IAC3D,MAAM,cAAc,GAAG,IAAI;SACtB,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,KAAK,CAAC,GAAG,CAAC,MAAM,CAAC,EAAE,IAAI,EAAE,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC;SAC3D,GAAG,CAAC,8BAA8B,CAAC;SACnC,GAAG,CAAC,iBAAiB,CAAC,CAAC;IAE5B,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,EAAE,cAAc,EAAE,gBAAgB,EAAE,IAAI,CAAC,CAAC;IAC7D,SAAS;IAET,yEAAyE;IACzE,oBAAoB;IACpB,OAAO,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,CAAC,KAAK,CAAC,GAAG,CAAC,MAAM,CAAC,EAAE,CAAC,MAAM,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;AAC9D,CAAC;AAGD;;;;;GAKG;AACH,SAAS,eAAe,CAChB,MAAc,EACd,kBAAsC;IAE1C,OAAO,UAAU,MAAc;QAC3B,MAAM,OAAO,GAAG,MAAM,CAAC,IAAI,CAAC;QAE5B,MAAM,GAAG,GAAG,MAAM,CAAC,YAAY,CAAC;QAChC,MAAM,IAAI,GAAG,OAAO,CAAC,YAAY,CAAC;QAElC,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE,GAAG,GAAG,CAAC;QACjC,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,EAAE,GAAG,IAAI,CAAC;QAEtC,IAAI,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;YACnC,OAAO,SAAS,CAAC;QACrB,CAAC;QAED,IAAI,KAAK,CAAC,IAAI,KAAK,MAAM,CAAC,IAAI,EAAE,CAAC;YAC7B,mDAAmD;YACnD,oBAAoB;YACpB,OAAO,SAAS,CAAC;QACrB,CAAC;QAED,MAAM,CAAC,GAAG,eAAe,CAAC,MAAM,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;QAG5C,kEAAkE;QAClE,8EAA8E;QAC9E,YAAY;QACZ,wBAAwB;QACxB,wBAAwB;QACxB,wBAAwB;QACxB,YAAY;QACZ,uBAAuB;QACvB,oCAAoC;QACpC,IAAI;QACJ,8EAA8E;QAC9E,YAAY;QACZ,wBAAwB;QACxB,wBAAwB;QACxB,wBAAwB;QACxB,YAAY;QACZ,uBAAuB;QACvB,oCAAoC;QACpC,IAAI;QACJ,IAAI,kBAAkB,KAAK,SAAS,EAAE,CAAC;YACnC,MAAM,GAAG,GAAG,CAAC,IAAE,CAAC,MAAM,GAAG,kBAAkB,CAAC,CAAC;YAE7C,IAAI,GAAG,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC;gBAC/B,MAAM,IAAI,KAAK,CACX,mDAAmD,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI;oBACvE,2DAA2D;oBAC3D,wBAAwB,MAAM,CAAC,MAAM,eAAe,CAAC,EAAE,CAC1D,CAAC;YACN,CAAC;QACL,CAAC;QAED,OAAO,gBAAgB,CAAC,MAAM,CAAC,CAAC;IACpC,CAAC,CAAA;AACL,CAAC;AAGD;;;;;;GAMG;AACH,SAAS,gBAAgB,CACjB,MAAc;IAElB,OAAO,sBAAsB,CACzB,sBAAsB,CAAC,MAAM,CAAC,YAAY,EAAE,IAAI,CAAC,EACjD,sBAAsB,CAAC,MAAM,CAAC,IAAI,CAAC,YAAY,EAAE,KAAK,CAAC,CAC1D,CAAC;AACN,CAAC;AAGD,OAAO,EAAE,mBAAmB,EAAE,CAAA"}

package/node/reconstruct-from-mat/remove-sat-induced-trivial-3prongs.d.ts

@@ -3,6 +3,10 @@ import type { Mat } from '../mat/mat.js';
  * In-place removes trivial 3-prongs from the SAT graph; just to simplify
  * subsequent processing.
  *
+ * * **mutates** the given `Mat` (its `CpNode`s are relinked and `cpNode` is
+ *   reassigned; `meta.cpTrees` is left stale) - pass a clone if the original
+ *   must be kept
+ *
  * These are 3-prongs that are induced by the SAT construction, and do
  * not correspond to actual 3-prongs in the shape.
  *

package/node/reconstruct-from-mat/remove-sat-induced-trivial-3prongs.js

@@ -3,6 +3,10 @@ import { isFullyTerminating } from '../cp-node/fs/is-fully-terminating.js';
  * In-place removes trivial 3-prongs from the SAT graph; just to simplify
  * subsequent processing.
  *
+ * * **mutates** the given `Mat` (its `CpNode`s are relinked and `cpNode` is
+ *   reassigned; `meta.cpTrees` is left stale) - pass a clone if the original
+ *   must be kept
+ *
  * These are 3-prongs that are induced by the SAT construction, and do
  * not correspond to actual 3-prongs in the shape.
  *

package/node/reconstruct-from-mat/remove-sat-induced-trivial-3prongs.js.map

@@ -1 +1 @@
-{"version":3,"file":"remove-sat-induced-trivial-3prongs.js","sourceRoot":"","sources":["../../src/reconstruct-from-mat/remove-sat-induced-trivial-3prongs.ts"],"names":[],"mappings":"AAEA,OAAO,EAAE,kBAAkB,EAAE,MAAM,uCAAuC,CAAC;AAG3E;;;;;;;;GAQG;AACH,SAAS,8BAA8B,CAC/B,GAAQ;IAEZ,IAAI,EAAE,MAAM,EAAE,GAAG,GAAG,CAAC;IAErB,IAAI,OAAO,GAAG,MAAM,CAAC;IACrB,GAAG,CAAC;QACA,IAAI,kBAAkB,CAAC,MAAM,CAAC;YAC1B,iEAAiE;YACjE,MAAM,CAAC,IAAI,KAAK,MAAM,CAAC,YAAY,EAAE,CAAC;YAEtC,MAAM,GAAG,MAAM,CAAC,IAAI,CAAC;YACrB,SAAS;QACb,CAAC;QAED,+BAA+B;QAC/B,6BAA6B;QAC7B,+BAA+B;QAC/B,uDAAuD;QAEvD,MAAM,IAAI,GAAG,MAAM,CAAC,IAAI,CAAC;QACzB,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC;QAC3B,MAAM,UAAU,GAAG,MAAM,CAAC,IAAI,CAAC,YAAY,CAAC;QAE5C,+DAA+D;QAC/D,8DAA8D;QAC9D,+DAA+D;QAC9D,MAAiC,CAAC,IAAI,GAAG,QAAQ,CAAC;QAClD,MAAiC,CAAC,YAAY,GAAG,UAAU,CAAC;QAC5D,QAAqC,CAAC,IAAI,GAAG,MAAM,CAAC;QACpD,UAAyC,CAAC,YAAY,GAAG,MAAM,CAAC;QACjE,+DAA+D;QAE/D,IAAI,OAAO,KAAK,IAAI,EAAE,CAAC;YACnB,OAAO,GAAG,MAAM,CAAC,CAAE,6CAA6C;QACpE,CAAC;QAED,MAAM,GAAG,MAAM,CAAC,IAAI,CAAC;IACzB,CAAC,QAAQ,MAAM,KAAK,OAAO,EAAC;IAE5B,+FAA+F;IAC9F,GAAoB,CAAC,MAAM,GAAG,MAAM,CAAC;IAEtC,OAAO,GAAG,CAAC;AACf,CAAC;AAGD,OAAO,EAAE,8BAA8B,EAAE,CAAA"}
+{"version":3,"file":"remove-sat-induced-trivial-3prongs.js","sourceRoot":"","sources":["../../src/reconstruct-from-mat/remove-sat-induced-trivial-3prongs.ts"],"names":[],"mappings":"AAEA,OAAO,EAAE,kBAAkB,EAAE,MAAM,uCAAuC,CAAC;AAG3E;;;;;;;;;;;;GAYG;AACH,SAAS,8BAA8B,CAC/B,GAAQ;IAEZ,IAAI,EAAE,MAAM,EAAE,GAAG,GAAG,CAAC;IAErB,IAAI,OAAO,GAAG,MAAM,CAAC;IACrB,GAAG,CAAC;QACA,IAAI,kBAAkB,CAAC,MAAM,CAAC;YAC1B,iEAAiE;YACjE,MAAM,CAAC,IAAI,KAAK,MAAM,CAAC,YAAY,EAAE,CAAC;YAEtC,MAAM,GAAG,MAAM,CAAC,IAAI,CAAC;YACrB,SAAS;QACb,CAAC;QAED,+BAA+B;QAC/B,6BAA6B;QAC7B,+BAA+B;QAC/B,uDAAuD;QAEvD,MAAM,IAAI,GAAG,MAAM,CAAC,IAAI,CAAC;QACzB,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC;QAC3B,MAAM,UAAU,GAAG,MAAM,CAAC,IAAI,CAAC,YAAY,CAAC;QAE5C,+DAA+D;QAC/D,8DAA8D;QAC9D,+DAA+D;QAC9D,MAAiC,CAAC,IAAI,GAAG,QAAQ,CAAC;QAClD,MAAiC,CAAC,YAAY,GAAG,UAAU,CAAC;QAC5D,QAAqC,CAAC,IAAI,GAAG,MAAM,CAAC;QACpD,UAAyC,CAAC,YAAY,GAAG,MAAM,CAAC;QACjE,+DAA+D;QAE/D,IAAI,OAAO,KAAK,IAAI,EAAE,CAAC;YACnB,OAAO,GAAG,MAAM,CAAC,CAAE,6CAA6C;QACpE,CAAC;QAED,MAAM,GAAG,MAAM,CAAC,IAAI,CAAC;IACzB,CAAC,QAAQ,MAAM,KAAK,OAAO,EAAC;IAE5B,+FAA+F;IAC9F,GAAoB,CAAC,MAAM,GAAG,MAAM,CAAC;IAEtC,OAAO,GAAG,CAAC;AACf,CAAC;AAGD,OAAO,EAAE,8BAA8B,EAAE,CAAA"}

package/node/sat/cull.d.ts

@@ -1,14 +1,15 @@
 import type { CpNode } from '../cp-node/cp-node.js';
-import type { DualSet } from '../utils/dual-set.js';
+import type { Circle } from '../geometry/circle.js';
 /**
- * Returns the set of Vertices passing the following test: walk the MAT tree and
- * keep all Vertices not in the current cull set and any Vertices that have a
- * non-culled node further down the line toward the tree leaves.
+ * Culls (in place) the vertices of the given MAT tree: walking the tree from
+ * the maximal vertex, all vertices in the cull set are removed except those
+ * that have a non-culled vertex further down the line toward the tree leaves
+ * (which are kept to preserve the topology).
  *
- * @param culls the `CpNode`s (referred to by circles) that should be culled
+ * @param culls the vertices (identified by their `Circle`) that should be culled
  * @param maxCpNode the start `CpNode` which must represent the maximal vertex
  *
  * @internal
  */
-declare function cull(culls: DualSet<number, number>, maxCpNode: CpNode): void;
+declare function cull(culls: Set<Circle>, maxCpNode: CpNode): void;
 export { cull };

package/node/sat/cull.js

@@ -1,6 +1,5 @@
 import { getLeaves } from './get-leaves.js';
 import { isTerminating } from '../cp-node/fs/is-terminating.js';
-import { DualSetFs } from '../utils/dual-set.js';
 import { isOnSameCircle } from '../cp-node/fs/is-on-same-circle.js';
 /**
  * Returns the set of non-trivial forward edges starting from the given `CpNode`.
@@ -23,11 +22,12 @@ function getNonTrivialForwardEdges(cpStart) {
     return cps;
 }
 /**
- * Returns the set of Vertices passing the following test: walk the MAT tree and
- * keep all Vertices not in the current cull set and any Vertices that have a
- * non-culled node further down the line toward the tree leaves.
+ * Culls (in place) the vertices of the given MAT tree: walking the tree from
+ * the maximal vertex, all vertices in the cull set are removed except those
+ * that have a non-culled vertex further down the line toward the tree leaves
+ * (which are kept to preserve the topology).
  *
- * @param culls the `CpNode`s (referred to by circles) that should be culled
+ * @param culls the vertices (identified by their `Circle`) that should be culled
  * @param maxCpNode the start `CpNode` which must represent the maximal vertex
  *
  * @internal
@@ -36,8 +36,7 @@ function cull(culls, maxCpNode) {
     const leaves = getLeaves(maxCpNode);
     while (leaves.length) {
         const leaf = leaves.pop();
-        const { center: c } = leaf.pointOnShape.circle;
-        if (!DualSetFs.has(culls, c[0], c[1]) ||
+        if (!culls.has(leaf.pointOnShape.circle) ||
             // Preserve topology.
             leaf.isIntersection) {
             continue;
@@ -46,8 +45,7 @@ function cull(culls, maxCpNode) {
         let cut = false;
         while (!cut) {
             cpNode = cpNode.next;
-            const { center: c } = cpNode.pointOnShape.circle;
-            if (!DualSetFs.has(culls, c[0], c[1])) {
+            if (!culls.has(cpNode.pointOnShape.circle)) {
                 // Cut off the edge once a non-cull has been reached.
                 while (isTerminating(cpNode.prevOnCircle)) {
                     cpNode = cpNode.prevOnCircle;
@@ -67,13 +65,12 @@ function cull(culls, maxCpNode) {
                 }
             }
         }
-        if (cut) {
-            const cp1 = cpNode.prevOnCircle;
-            cp1.next = cpNode;
-            cpNode.prev = cp1;
-            cp1.nextOnCircle = cpNode;
-            cpNode.prevOnCircle = cp1;
-        }
+        // Cut the edge: link the node we stopped at directly to the one before
+        // it on its circle (`cut` is necessarily `true` here; the on-circle
+        // links are already consistent by the ring invariant)
+        const cp1 = cpNode.prevOnCircle;
+        cp1.next = cpNode;
+        cpNode.prev = cp1;
     }
 }
 export { cull };

package/node/sat/cull.js.map

@@ -1 +1 @@
-{"version":3,"file":"cull.js","sourceRoot":"","sources":["../../src/sat/cull.ts"],"names":[],"mappings":"AAGA,OAAO,EAAE,SAAS,EAAE,MAAM,iBAAiB,CAAC;AAC5C,OAAO,EAAE,aAAa,EAAE,MAAM,iCAAiC,CAAC;AAChE,OAAO,EAAE,SAAS,EAAE,MAAM,sBAAsB,CAAC;AACjD,OAAO,EAAE,cAAc,EAAE,MAAM,oCAAoC,CAAC;AAGpE;;;;GAIG;AACH,SAAS,yBAAyB,CAAC,OAAe;IAC9C,IAAI,MAAM,GAAG,OAAO,CAAC;IACrB,IAAI,KAAK,GAAG,OAAO,CAAC,YAAY,CAAC;IACjC,OAAO,aAAa,CAAC,KAAK,CAAC,EAAE,CAAC;QAC1B,KAAK,GAAG,KAAK,CAAC,YAAY,CAAC;IAC/B,CAAC;IAED,MAAM,GAAG,GAAa,EAAE,CAAC;IACzB,GAAG,CAAC;QACA,IAAI,CAAC,aAAa,CAAC,MAAM,CAAC,EAAE,CAAC;YACzB,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QACrB,CAAC;QACD,MAAM,GAAG,MAAM,CAAC,YAAY,CAAC;IACjC,CAAC,QAAQ,MAAM,KAAK,KAAK,EAAC;IAE1B,OAAO,GAAG,CAAC;AACf,CAAC;AAGD;;;;;;;;;GASG;AACH,SAAS,IAAI,CACL,KAA8B,EAC9B,SAAiB;IAErB,MAAM,MAAM,GAAG,SAAS,CAAC,SAAS,CAAC,CAAC;IAEpC,OAAO,MAAM,CAAC,MAAM,EAAE,CAAC;QACnB,MAAM,IAAI,GAAG,MAAM,CAAC,GAAG,EAAG,CAAC;QAE3B,MAAM,EAAE,MAAM,EAAE,CAAC,EAAE,GAAG,IAAI,CAAC,YAAY,CAAC,MAAM,CAAC;QAC/C,IAAI,CAAC,SAAS,CAAC,GAAG,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;YACjC,qBAAqB;YACrB,IAAI,CAAC,cAAc,EAAE,CAAC;YAEtB,SAAS;QACb,CAAC;QAED,IAAI,MAAM,GAAG,IAAI,CAAC,YAAY,CAAC,CAAC,cAAc;QAC9C,IAAI,GAAG,GAAG,KAAK,CAAC;QAChB,OAAO,CAAC,GAAG,EAAE,CAAC;YACV,MAAM,GAAG,MAAM,CAAC,IAAI,CAAC;YAErB,MAAM,EAAE,MAAM,EAAE,CAAC,EAAE,GAAG,MAAM,CAAC,YAAY,CAAC,MAAM,CAAC;YACjD,IAAI,CAAC,SAAS,CAAC,GAAG,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;gBACpC,qDAAqD;gBACrD,OAAO,aAAa,CAAC,MAAM,CAAC,YAAY,CAAC,EAAE,CAAC;oBACxC,MAAM,GAAG,MAAM,CAAC,YAAY,CAAC;gBACjC,CAAC;gBACD,GAAG,GAAG,IAAI,CAAC;YACf,CAAC;iBAAM,IAAI,cAAc,CAAC,MAAM,EAAE,SAAS,CAAC,EAAE,CAAC;gBAC3C,GAAG,GAAG,IAAI,CAAC,CAAC,0CAA0C;YAC1D,CAAC;iBAAM,CAAC;gBACJ,MAAM,GAAG,GAAG,yBAAyB,CAAC,MAAM,CAAC,CAAC;gBAE9C,IAAI,GAAG,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;oBACnB,SAAS;gBACb,CAAC;qBAAM,CAAC;oBACJ,GAAG,GAAG,IAAI,CAAC;gBACf,CAAC;YACL,CAAC;QACL,CAAC;QAED,IAAI,GAAG,EAAE,CAAC;YACN,MAAM,GAAG,GAAG,MAAM,CAAC,YAA+B,CAAC;YACnD,GAAG,CAAC,IAAI,GAAG,MAAM,CAAC;YACjB,MAA0B,CAAC,IAAI,GAAG,GAAG,CAAC;YAEvC,GAAG,CAAC,YAAY,GAAG,MAAM,CAAC;YACzB,MAA0B,CAAC,YAAY,GAAG,GAAG,CAAC;QACnD,CAAC;IACL,CAAC;AACL,CAAC;AAGD,OAAO,EAAE,IAAI,EAAE,CAAA"}
+{"version":3,"file":"cull.js","sourceRoot":"","sources":["../../src/sat/cull.ts"],"names":[],"mappings":"AAGA,OAAO,EAAE,SAAS,EAAE,MAAM,iBAAiB,CAAC;AAC5C,OAAO,EAAE,aAAa,EAAE,MAAM,iCAAiC,CAAC;AAChE,OAAO,EAAE,cAAc,EAAE,MAAM,oCAAoC,CAAC;AAGpE;;;;GAIG;AACH,SAAS,yBAAyB,CAAC,OAAe;IAC9C,IAAI,MAAM,GAAG,OAAO,CAAC;IACrB,IAAI,KAAK,GAAG,OAAO,CAAC,YAAY,CAAC;IACjC,OAAO,aAAa,CAAC,KAAK,CAAC,EAAE,CAAC;QAC1B,KAAK,GAAG,KAAK,CAAC,YAAY,CAAC;IAC/B,CAAC;IAED,MAAM,GAAG,GAAa,EAAE,CAAC;IACzB,GAAG,CAAC;QACA,IAAI,CAAC,aAAa,CAAC,MAAM,CAAC,EAAE,CAAC;YACzB,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QACrB,CAAC;QACD,MAAM,GAAG,MAAM,CAAC,YAAY,CAAC;IACjC,CAAC,QAAQ,MAAM,KAAK,KAAK,EAAC;IAE1B,OAAO,GAAG,CAAC;AACf,CAAC;AAGD;;;;;;;;;;GAUG;AACH,SAAS,IAAI,CACL,KAAkB,EAClB,SAAiB;IAErB,MAAM,MAAM,GAAG,SAAS,CAAC,SAAS,CAAC,CAAC;IAEpC,OAAO,MAAM,CAAC,MAAM,EAAE,CAAC;QACnB,MAAM,IAAI,GAAG,MAAM,CAAC,GAAG,EAAG,CAAC;QAE3B,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,IAAI,CAAC,YAAY,CAAC,MAAM,CAAC;YACpC,qBAAqB;YACrB,IAAI,CAAC,cAAc,EAAE,CAAC;YAEtB,SAAS;QACb,CAAC;QAED,IAAI,MAAM,GAAG,IAAI,CAAC,YAAY,CAAC,CAAC,cAAc;QAC9C,IAAI,GAAG,GAAG,KAAK,CAAC;QAChB,OAAO,CAAC,GAAG,EAAE,CAAC;YACV,MAAM,GAAG,MAAM,CAAC,IAAI,CAAC;YAErB,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,MAAM,CAAC,YAAY,CAAC,MAAM,CAAC,EAAE,CAAC;gBACzC,qDAAqD;gBACrD,OAAO,aAAa,CAAC,MAAM,CAAC,YAAY,CAAC,EAAE,CAAC;oBACxC,MAAM,GAAG,MAAM,CAAC,YAAY,CAAC;gBACjC,CAAC;gBACD,GAAG,GAAG,IAAI,CAAC;YACf,CAAC;iBAAM,IAAI,cAAc,CAAC,MAAM,EAAE,SAAS,CAAC,EAAE,CAAC;gBAC3C,GAAG,GAAG,IAAI,CAAC,CAAC,0CAA0C;YAC1D,CAAC;iBAAM,CAAC;gBACJ,MAAM,GAAG,GAAG,yBAAyB,CAAC,MAAM,CAAC,CAAC;gBAE9C,IAAI,GAAG,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;oBACnB,SAAS;gBACb,CAAC;qBAAM,CAAC;oBACJ,GAAG,GAAG,IAAI,CAAC;gBACf,CAAC;YACL,CAAC;QACL,CAAC;QAED,uEAAuE;QACvE,oEAAoE;QACpE,sDAAsD;QACtD,MAAM,GAAG,GAAG,MAAM,CAAC,YAA+B,CAAC;QACnD,GAAG,CAAC,IAAI,GAAG,MAAM,CAAC;QACjB,MAA0B,CAAC,IAAI,GAAG,GAAG,CAAC;IAC3C,CAAC;AACL,CAAC;AAGD,OAAO,EAAE,IAAI,EAAE,CAAA"}

package/node/sat/get-sat-culls.d.ts

@@ -1,4 +1,14 @@
-import type { DualSet } from "../utils/dual-set.js";
 import type { CpNode } from "../cp-node/cp-node.js";
-declare function getSatCulls(cpNode: CpNode, s: number): DualSet<number, number>;
+import type { Circle } from "../geometry/circle.js";
+/**
+ * Returns the set of MAT vertices (identified by their `Circle` - all `CpNode`s
+ * on a vertex share the same `Circle` object) that are to be culled by the
+ * Scale Axis Transform with the given scale factor.
+ *
+ * @param cpNode the `CpNode` representing the maximal vertex
+ * @param s the scale factor
+ *
+ * @internal
+ */
+declare function getSatCulls(cpNode: CpNode, s: number): Set<Circle>;
 export { getSatCulls };

package/node/sat/get-sat-culls.js

@@ -1,29 +1,34 @@
-import { length } from "flo-bezier3";
-import { getAllOnCircle } from "../cp-node/fs/get-all-on-circle.js";
-import { getMatCurveToNext } from "../cp-node/fs/get-mat-curve-to-next.js";
+import { getMatDistanceToNext } from "../cp-node/fs/get-distance-to-next.js";
 import { traverseEdges } from "../cp-node/fs/traverse-edges.js";
-import { DualSetFs } from "../utils/dual-set.js";
+/**
+ * Returns the set of MAT vertices (identified by their `Circle` - all `CpNode`s
+ * on a vertex share the same `Circle` object) that are to be culled by the
+ * Scale Axis Transform with the given scale factor.
+ *
+ * @param cpNode the `CpNode` representing the maximal vertex
+ * @param s the scale factor
+ *
+ * @internal
+ */
 function getSatCulls(cpNode, s) {
     /**
      * All vertices that are set to be culled initially. This may change later
      * in order to preserve topology.
      */
-    const culls = new Map();
+    const culls = new Set();
+    /** The osculating radius stored with each vertex. */
     const rMap = new Map();
     traverseEdges(cpNode, function (cpNode) {
-        /** The occulating radius stored with this vertex. */
-        const R = rMap.get(cpNode) || (s * cpNode.pointOnShape.circle.radius);
+        const circle = cpNode.pointOnShape.circle;
+        const R = rMap.get(circle) ?? s * circle.radius;
         const cpNode_ = cpNode.next;
-        const l = length([0, 1], getMatCurveToNext(cpNode));
-        const r = cpNode_.pointOnShape.circle.radius;
-        const r_ = s * r;
+        const circle_ = cpNode_.pointOnShape.circle;
+        const l = getMatDistanceToNext(cpNode);
+        const r_ = s * circle_.radius;
         if (R - l > r_) {
-            for (const cpNode of getAllOnCircle(cpNode_)) {
-                rMap.set(cpNode, R - l); // Update osculating radii
-                if (!cpNode_.isHoleClosing) {
-                    const { center: c } = cpNode_.pointOnShape.circle;
-                    DualSetFs.add(culls, c[0], c[1]);
-                }
+            rMap.set(circle_, R - l); // Update osculating radius
+            if (!cpNode_.isHoleClosing) {
+                culls.add(circle_);
             }
         }
     });

package/node/sat/get-sat-culls.js.map

@@ -1 +1 @@
-{"version":3,"file":"get-sat-culls.js","sourceRoot":"","sources":["../../src/sat/get-sat-culls.ts"],"names":[],"mappings":"AAEA,OAAO,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AACrC,OAAO,EAAE,cAAc,EAAE,MAAM,oCAAoC,CAAC;AACpE,OAAO,EAAE,iBAAiB,EAAE,MAAM,wCAAwC,CAAC;AAC3E,OAAO,EAAE,aAAa,EAAE,MAAM,iCAAiC,CAAC;AAChE,OAAO,EAAE,SAAS,EAAE,MAAM,sBAAsB,CAAC;AAGjD,SAAS,WAAW,CACZ,MAAc,EACd,CAAS;IAEb;;;OAGG;IACH,MAAM,KAAK,GAA2B,IAAI,GAAG,EAAE,CAAC;IAEhD,MAAM,IAAI,GAAuB,IAAI,GAAG,EAAE,CAAC;IAE3C,aAAa,CAAC,MAAM,EAAE,UAAS,MAAM;QACjC,qDAAqD;QACrD,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC,GAAC,MAAM,CAAC,YAAY,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC;QAEpE,MAAM,OAAO,GAAG,MAAM,CAAC,IAAI,CAAC;QAE5B,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,EAAC,CAAC,CAAC,EAAE,iBAAiB,CAAC,MAAM,CAAC,CAAC,CAAC;QAEnD,MAAM,CAAC,GAAG,OAAO,CAAC,YAAY,CAAC,MAAM,CAAC,MAAM,CAAC;QAC7C,MAAM,EAAE,GAAG,CAAC,GAAC,CAAC,CAAC;QACf,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,EAAE,CAAC;YACb,KAAK,MAAM,MAAM,IAAI,cAAc,CAAC,OAAO,CAAC,EAAE,CAAC;gBAC3C,IAAI,CAAC,GAAG,CAAC,MAAM,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,0BAA0B;gBACnD,IAAI,CAAC,OAAO,CAAC,aAAa,EAAE,CAAC;oBACzB,MAAM,EAAE,MAAM,EAAE,CAAC,EAAE,GAAG,OAAO,CAAC,YAAY,CAAC,MAAM,CAAC;oBAClD,SAAS,CAAC,GAAG,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACrC,CAAC;YACL,CAAC;QACL,CAAC;IACL,CAAC,CAAC,CAAC;IAEH,OAAO,KAAK,CAAC;AACjB,CAAC;AAGD,OAAO,EAAE,WAAW,EAAE,CAAA"}
+{"version":3,"file":"get-sat-culls.js","sourceRoot":"","sources":["../../src/sat/get-sat-culls.ts"],"names":[],"mappings":"AAEA,OAAO,EAAE,oBAAoB,EAAE,MAAM,uCAAuC,CAAC;AAC7E,OAAO,EAAE,aAAa,EAAE,MAAM,iCAAiC,CAAC;AAGhE;;;;;;;;;GASG;AACH,SAAS,WAAW,CACZ,MAAc,EACd,CAAS;IAEb;;;OAGG;IACH,MAAM,KAAK,GAAG,IAAI,GAAG,EAAU,CAAC;IAEhC,qDAAqD;IACrD,MAAM,IAAI,GAAG,IAAI,GAAG,EAAkB,CAAC;IAEvC,aAAa,CAAC,MAAM,EAAE,UAAS,MAAM;QACjC,MAAM,MAAM,GAAG,MAAM,CAAC,YAAY,CAAC,MAAM,CAAC;QAC1C,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,IAAI,CAAC,GAAC,MAAM,CAAC,MAAM,CAAC;QAE9C,MAAM,OAAO,GAAG,MAAM,CAAC,IAAI,CAAC;QAC5B,MAAM,OAAO,GAAG,OAAO,CAAC,YAAY,CAAC,MAAM,CAAC;QAE5C,MAAM,CAAC,GAAG,oBAAoB,CAAC,MAAM,CAAC,CAAC;QAEvC,MAAM,EAAE,GAAG,CAAC,GAAC,OAAO,CAAC,MAAM,CAAC;QAC5B,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,EAAE,CAAC;YACb,IAAI,CAAC,GAAG,CAAC,OAAO,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,2BAA2B;YACrD,IAAI,CAAC,OAAO,CAAC,aAAa,EAAE,CAAC;gBACzB,KAAK,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;YACvB,CAAC;QACL,CAAC;IACL,CAAC,CAAC,CAAC;IAEH,OAAO,KAAK,CAAC;AACjB,CAAC;AAGD,OAAO,EAAE,WAAW,EAAE,CAAA"}

package/node/svg/get-shape-bounds.d.ts

@@ -1,21 +1,10 @@
 import type { Loop } from 'flo-boolean';
 import type { PrePointOnShape } from '../point-on-shape/point-on-shape.js';
 /** @internal */
-declare const getShapeBounds: ((loops: Loop[]) => {
+declare function getShapeBounds(loops: Loop[]): {
     minX: PrePointOnShape;
     minY: PrePointOnShape;
     maxX: PrePointOnShape;
     maxY: PrePointOnShape;
-}) & {
-    readonly weakMapS: WeakMap<object, {
-        readonly weakMap: WeakMap<object, any>;
-        readonly map: Map<object, any>;
-    }>;
-    readonly mapS: Map<object, {
-        readonly weakMap: WeakMap<object, any>;
-        readonly map: Map<object, any>;
-    }>;
-    readonly clearCache: () => void;
-    readonly addToCache: (r: unknown, ...args: any) => void;
 };
 export { getShapeBounds };

package/node/svg/get-shape-bounds.js

@@ -1,7 +1,6 @@
-import { memoize } from 'flo-memoize';
 import { getLoopBounds$ } from './get-loop-bounds.js';
 /** @internal */
-const getShapeBounds = memoize(function (loops) {
+function getShapeBounds(loops) {
     let minX_ = +Infinity;
     let maxX_ = -Infinity;
     let minY_ = +Infinity;
@@ -35,6 +34,6 @@ const getShapeBounds = memoize(function (loops) {
         maxX: maxX,
         maxY: maxY
     };
-});
+}
 export { getShapeBounds };
 //# sourceMappingURL=get-shape-bounds.js.map

package/node/svg/get-shape-bounds.js.map

@@ -1 +1 @@
-{"version":3,"file":"get-shape-bounds.js","sourceRoot":"","sources":["../../src/svg/get-shape-bounds.ts"],"names":[],"mappings":"AAEA,OAAO,EAAE,OAAO,EAAE,MAAM,aAAa,CAAC;AACtC,OAAO,EAAE,cAAc,EAAE,MAAM,sBAAsB,CAAC;AAGtD,gBAAgB;AAChB,MAAM,cAAc,GAAG,OAAO,CAAC,UAAS,KAAa;IACjD,IAAI,KAAK,GAAG,CAAC,QAAQ,CAAC;IACtB,IAAI,KAAK,GAAG,CAAC,QAAQ,CAAC;IACtB,IAAI,KAAK,GAAG,CAAC,QAAQ,CAAC;IACtB,IAAI,KAAK,GAAG,CAAC,QAAQ,CAAC;IAEtB,IAAI,IAAqB,CAAC;IAC1B,IAAI,IAAqB,CAAC;IAC1B,IAAI,IAAqB,CAAC;IAC1B,IAAI,IAAqB,CAAC;IAE1B,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACvB,MAAM,MAAM,GAAG,cAAc,CAAC,IAAI,CAAC,CAAC;QACpC,IAAI,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,EAAE,CAAC;YAC3B,IAAI,GAAG,MAAM,CAAC,IAAI,CAAC;YACnB,KAAK,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7B,CAAC;QACD,IAAI,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,EAAE,CAAC;YAC3B,IAAI,GAAG,MAAM,CAAC,IAAI,CAAC;YACnB,KAAK,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7B,CAAC;QACD,IAAI,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,EAAE,CAAC;YAC3B,IAAI,GAAG,MAAM,CAAC,IAAI,CAAC;YACnB,KAAK,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7B,CAAC;QACD,IAAI,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,EAAE,CAAC;YAC3B,IAAI,GAAG,MAAM,CAAC,IAAI,CAAC;YACnB,KAAK,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7B,CAAC;IACL,CAAC;IAED,OAAO;QACH,IAAI,EAAE,IAAK;QACX,IAAI,EAAE,IAAK;QACX,IAAI,EAAE,IAAK;QACX,IAAI,EAAE,IAAK;KACd,CAAC;AACN,CAAC,CAAC,CAAC;AAGH,OAAO,EAAE,cAAc,EAAE,CAAA"}
+{"version":3,"file":"get-shape-bounds.js","sourceRoot":"","sources":["../../src/svg/get-shape-bounds.ts"],"names":[],"mappings":"AAEA,OAAO,EAAE,cAAc,EAAE,MAAM,sBAAsB,CAAC;AAGtD,gBAAgB;AAChB,SAAS,cAAc,CAAC,KAAa;IACjC,IAAI,KAAK,GAAG,CAAC,QAAQ,CAAC;IACtB,IAAI,KAAK,GAAG,CAAC,QAAQ,CAAC;IACtB,IAAI,KAAK,GAAG,CAAC,QAAQ,CAAC;IACtB,IAAI,KAAK,GAAG,CAAC,QAAQ,CAAC;IAEtB,IAAI,IAAqB,CAAC;IAC1B,IAAI,IAAqB,CAAC;IAC1B,IAAI,IAAqB,CAAC;IAC1B,IAAI,IAAqB,CAAC;IAE1B,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACvB,MAAM,MAAM,GAAG,cAAc,CAAC,IAAI,CAAC,CAAC;QACpC,IAAI,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,EAAE,CAAC;YAC3B,IAAI,GAAG,MAAM,CAAC,IAAI,CAAC;YACnB,KAAK,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7B,CAAC;QACD,IAAI,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,EAAE,CAAC;YAC3B,IAAI,GAAG,MAAM,CAAC,IAAI,CAAC;YACnB,KAAK,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7B,CAAC;QACD,IAAI,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,EAAE,CAAC;YAC3B,IAAI,GAAG,MAAM,CAAC,IAAI,CAAC;YACnB,KAAK,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7B,CAAC;QACD,IAAI,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,EAAE,CAAC;YAC3B,IAAI,GAAG,MAAM,CAAC,IAAI,CAAC;YACnB,KAAK,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7B,CAAC;IACL,CAAC;IAED,OAAO;QACH,IAAI,EAAE,IAAK;QACX,IAAI,EAAE,IAAK;QACX,IAAI,EAAE,IAAK;QACX,IAAI,EAAE,IAAK;KACd,CAAC;AACN,CAAC;AAGD,OAAO,EAAE,cAAc,EAAE,CAAA"}

package/node/sweep-line/sweep-line.js.map

@@ -1 +1 @@
-{"version":3,"file":"sweep-line.js","sourceRoot":"","sources":["../../src/sweep-line/sweep-line.ts"],"names":[],"mappings":"AACA,gBAAgB;AAChB,MAAM,UAAU,GAAI,CAAC,CAAC;AACtB,gBAAgB;AAChB,MAAM,WAAW,GAAG,CAAC,CAAC;AAOtB;;;;;;;;;;;;;;;GAeG;AACH,SAAS,SAAS,CACV,KAAU,EACV,WAAgC,EAChC,YAAiC,EACjC,SAA0C;IAE9C,mDAAmD;IACnD,MAAM,MAAM,GAAgB,EAAE,CAAC;IAC/B,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACvB,MAAM,CAAC,IAAI,CAAC;YACR,IAAI,EAAE,UAAU;YAChB,IAAI;YACJ,CAAC,EAAE,WAAW,CAAC,IAAI,CAAC;SACvB,CAAC,CAAC;QACH,MAAM,CAAC,IAAI,CAAC;YACR,IAAI,EAAE,WAAW;YACjB,IAAI;YACJ,CAAC,EAAE,YAAY,CAAC,IAAI,CAAC;SACxB,CAAC,CAAC;IACP,CAAC;IAED,MAAM,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;IAErB,MAAM,WAAW,GAAG,IAAI,GAAG,EAAK,CAAC;IAEjC,yDAAyD;IACzD,MAAM,WAAW,GAAY,EAAE,CAAC;IAChC,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;QACzB,MAAM,EAAE,IAAI,EAAE,GAAG,KAAK,CAAC;QAEpB,IAAI,KAAK,CAAC,IAAI,KAAK,UAAU,EAAE,CAAC;YAC5B,KAAK,MAAM,UAAU,IAAI,WAAW,CAAC,MAAM,EAAE,EAAE,CAAC;gBAC5C,IAAI,SAAS,CAAC,IAAI,EAAE,UAAU,CAAC,EAAE,CAAC;oBACjC,WAAW,CAAC,IAAI,CAAC,CAAC,IAAI,EAAE,UAAU,CAAC,CAAC,CAAC;gBACtC,CAAC;YACL,CAAC;YAED,WAAW,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QAC1B,CAAC;aAAM,IAAI,KAAK,CAAC,IAAI,KAAK,WAAW,EAAE,CAAC;YACpC,WAAW,CAAC,MAAM,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;QACnC,CAAC;IACR,CAAC;IAED,OAAO,WAAW,CAAC;AACvB,CAAC;AAeD;;;;;;;GAOG;AACH,SAAS,OAAO,CAAI,CAAY,EAAE,CAAY;IAC1C,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAEtB,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;QAAC,OAAO,GAAG,CAAC;IAAC,CAAC;IAE9B,2CAA2C;IAC3C,OAAO,CAAC,CAAC,IAAI,KAAK,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAC3C,CAAC;AAGD,OAAO,EAAE,SAAS,EAAE,CAAA"}
+{"version":3,"file":"sweep-line.js","sourceRoot":"","sources":["../../src/sweep-line/sweep-line.ts"],"names":[],"mappings":"AACA,gBAAgB;AAChB,MAAM,UAAU,GAAI,CAAC,CAAC;AACtB,gBAAgB;AAChB,MAAM,WAAW,GAAG,CAAC,CAAC;AAOtB;;;;;;;;;;;;;;;GAeG;AACH,SAAS,SAAS,CACV,KAAU,EACV,WAAgC,EAChC,YAAiC,EACjC,SAA0C;IAE9C,mDAAmD;IACnD,MAAM,MAAM,GAAgB,EAAE,CAAC;IAC/B,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACvB,MAAM,CAAC,IAAI,CAAC;YACR,IAAI,EAAE,UAAU;YAChB,IAAI;YACJ,CAAC,EAAE,WAAW,CAAC,IAAI,CAAC;SACvB,CAAC,CAAC;QACH,MAAM,CAAC,IAAI,CAAC;YACR,IAAI,EAAE,WAAW;YACjB,IAAI;YACJ,CAAC,EAAE,YAAY,CAAC,IAAI,CAAC;SACxB,CAAC,CAAC;IACP,CAAC;IAED,MAAM,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;IAErB,MAAM,WAAW,GAAG,IAAI,GAAG,EAAK,CAAC;IAEjC,yDAAyD;IACzD,MAAM,WAAW,GAAY,EAAE,CAAC;IAChC,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;QACzB,MAAM,EAAE,IAAI,EAAE,GAAG,KAAK,CAAC;QAEpB,IAAI,KAAK,CAAC,IAAI,KAAK,UAAU,EAAE,CAAC;YAC5B,KAAK,MAAM,UAAU,IAAI,WAAW,CAAC,MAAM,EAAE,EAAE,CAAC;gBAC5C,IAAI,SAAS,CAAC,IAAI,EAAE,UAAU,CAAC,EAAE,CAAC;oBAC7B,WAAW,CAAC,IAAI,CAAC,CAAC,IAAI,EAAE,UAAU,CAAC,CAAC,CAAC;gBAC1C,CAAC;YACL,CAAC;YAED,WAAW,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QAC1B,CAAC;aAAM,IAAI,KAAK,CAAC,IAAI,KAAK,WAAW,EAAE,CAAC;YACpC,WAAW,CAAC,MAAM,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;QACnC,CAAC;IACR,CAAC;IAED,OAAO,WAAW,CAAC;AACvB,CAAC;AAeD;;;;;;;GAOG;AACH,SAAS,OAAO,CAAI,CAAY,EAAE,CAAY;IAC1C,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAEtB,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;QAAC,OAAO,GAAG,CAAC;IAAC,CAAC;IAE9B,2CAA2C;IAC3C,OAAO,CAAC,CAAC,IAAI,KAAK,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAC3C,CAAC;AAGD,OAAO,EAAE,SAAS,EAAE,CAAA"}

package/node/utils/get-beziers-from-tangents.d.ts

@@ -0,0 +1,28 @@
+/**
+ * Returns an array of 1 or 2 bezier curves (each a line or a quadratic) that
+ * together form a smooth curve from the start point of `l0` to the start point
+ * of `l1` such that the curve is tangent to `l0` at its start and tangent to
+ * `l1` at its end.
+ *
+ * * `l0` and `l1` are given as `[p, p + dir]` where `p` is the curve endpoint
+ * and `dir` is the tangent direction (of any magnitude; its orientation does
+ * not matter)
+ * * if the tangent lines are not 'twisted' (they meet strictly between the
+ * endpoints along the chord) a single quadratic is returned
+ * * if the tangent lines are 'twisted' (the curve must have an inflection)
+ * 2 quadratics joined smoothly at the inflection are returned as an
+ * approximation - or, if both tangents deviate from the chord by less than
+ * `SIN_TOLERANCE_USE_LINE`, a single line
+ * * if the endpoints coincide (or a tangent is degenerate) a (degenerate) line
+ * is returned
+ *
+ * FUTURE: a strongly twisted MAT edge indicates that an additional 2-prong
+ * should be added so that the edge is approximated by two less twisted ones.
+ *
+ * @param l0 the tangent line at the start of the curve
+ * @param l1 the tangent line at the end of the curve
+ *
+ * @internal
+ */
+declare function getBeziersFromTangents(l0: number[][], l1: number[][]): number[][][];
+export { getBeziersFromTangents };

package/node/utils/get-beziers-from-tangents.js

@@ -0,0 +1,92 @@
+import { fromTo, dot, cross, lineLineIntersection } from 'flo-vector2d';
+const { abs, sqrt } = Math;
+/**
+ * The maximum sine of the angle between the chord and either tangent for
+ * which a 'twisted' (inflected) curve is approximated by a straight line -
+ * the line then deviates from the true curve by only about a millionth of the
+ * chord length.
+ *
+ * * a relative (angle) tolerance and thus independent of the scale of the
+ *   shape and of the magnitude of the given tangents
+ *
+ * @internal
+ */
+const SIN_TOLERANCE_USE_LINE = 2 ** -16;
+/**
+ * Returns an array of 1 or 2 bezier curves (each a line or a quadratic) that
+ * together form a smooth curve from the start point of `l0` to the start point
+ * of `l1` such that the curve is tangent to `l0` at its start and tangent to
+ * `l1` at its end.
+ *
+ * * `l0` and `l1` are given as `[p, p + dir]` where `p` is the curve endpoint
+ * and `dir` is the tangent direction (of any magnitude; its orientation does
+ * not matter)
+ * * if the tangent lines are not 'twisted' (they meet strictly between the
+ * endpoints along the chord) a single quadratic is returned
+ * * if the tangent lines are 'twisted' (the curve must have an inflection)
+ * 2 quadratics joined smoothly at the inflection are returned as an
+ * approximation - or, if both tangents deviate from the chord by less than
+ * `SIN_TOLERANCE_USE_LINE`, a single line
+ * * if the endpoints coincide (or a tangent is degenerate) a (degenerate) line
+ * is returned
+ *
+ * FUTURE: a strongly twisted MAT edge indicates that an additional 2-prong
+ * should be added so that the edge is approximated by two less twisted ones.
+ *
+ * @param l0 the tangent line at the start of the curve
+ * @param l1 the tangent line at the end of the curve
+ *
+ * @internal
+ */
+function getBeziersFromTangents(l0, l1) {
+    const p0 = l0[0];
+    const p1 = l1[0];
+    const c = fromTo(p0, p1);
+    if (c[0] === 0 && c[1] === 0) {
+        return [[p0, p1]]; // coinciding endpoints - a degenerate line
+    }
+    const w0 = fromTo(l0[0], l0[1]);
+    const w1 = fromTo(l1[0], l1[1]);
+    const lc = sqrt(dot(c, c));
+    const lw0 = sqrt(dot(w0, w0));
+    const lw1 = sqrt(dot(w1, w1));
+    if (lw0 === 0 || lw1 === 0) {
+        return [[p0, p1]]; // a degenerate tangent - fall back to a line
+    }
+    const mid = lineLineIntersection(l0, l1);
+    let twisted;
+    if (!mid) {
+        twisted = true;
+    }
+    else {
+        const a = fromTo(p0, mid);
+        const b = fromTo(p1, mid);
+        // (`<= 0` / `>= 0`: an intersection *at* an endpoint, which happens for
+        // collinear tangents, would give a degenerate quadratic with a zero
+        // end tangent - treat it as twisted so a line results instead)
+        twisted = dot(a, c) <= 0 || dot(b, c) >= 0;
+    }
+    if (!twisted) {
+        return [[p0, mid, p1]];
+    }
+    // the sines of the angles between the chord and the tangents
+    const sin0 = abs(cross(c, w0)) / (lc * lw0);
+    const sin1 = abs(cross(c, w1)) / (lc * lw1);
+    if (sin0 <= SIN_TOLERANCE_USE_LINE && sin1 <= SIN_TOLERANCE_USE_LINE) {
+        return [[p0, p1]]; // within tolerance - a straight line
+    }
+    // Approximate with 2 quadratic beziers joined smoothly at the inflection.
+    // The control points lie on the tangent lines a quarter of the chord
+    // length away from the endpoints, on the side facing the other endpoint
+    // (for small angles this is the midpoint degree reduction of the cubic
+    // having control points at 1/3 and 2/3 along the chord; unlike placing
+    // them at 1/4 and 3/4 *along* the chord it stays bounded for any angle).
+    const q0 = (dot(w0, c) < 0 ? -1 : 1) * lc / (4 * lw0);
+    const q1 = (dot(w1, c) < 0 ? 1 : -1) * lc / (4 * lw1);
+    const c0 = [p0[0] + q0 * w0[0], p0[1] + q0 * w0[1]];
+    const c1 = [p1[0] + q1 * w1[0], p1[1] + q1 * w1[1]];
+    const j = [(c0[0] + c1[0]) / 2, (c0[1] + c1[1]) / 2];
+    return [[p0, c0, j], [j, c1, p1]];
+}
+export { getBeziersFromTangents };
+//# sourceMappingURL=get-beziers-from-tangents.js.map

package/node/utils/get-beziers-from-tangents.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"get-beziers-from-tangents.js","sourceRoot":"","sources":["../../src/utils/get-beziers-from-tangents.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,MAAM,EAAE,GAAG,EAAE,KAAK,EAAE,oBAAoB,EAAE,MAAM,cAAc,CAAC;AAExE,MAAM,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,IAAI,CAAC;AAG3B;;;;;;;;;;GAUG;AACH,MAAM,sBAAsB,GAAG,CAAC,IAAE,CAAC,EAAE,CAAC;AAGtC;;;;;;;;;;;;;;;;;;;;;;;;;GAyBG;AACH,SAAS,sBAAsB,CACvB,EAAc,EACd,EAAc;IAElB,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IACjB,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IAEjB,MAAM,CAAC,GAAG,MAAM,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;IACxB,IAAI,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC;QAC3B,OAAO,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,CAAE,2CAA2C;IACnE,CAAC;IAED,MAAM,EAAE,GAAG,MAAM,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAChC,MAAM,EAAE,GAAG,MAAM,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAChC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAC,CAAC,CAAC,CAAC,CAAC;IAC1B,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC,CAAC;IAC7B,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC,CAAC;IAC7B,IAAI,GAAG,KAAK,CAAC,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;QACzB,OAAO,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,CAAE,6CAA6C;IACrE,CAAC;IAED,MAAM,GAAG,GAAG,oBAAoB,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IAEzC,IAAI,OAAgB,CAAC;IACrB,IAAI,CAAC,GAAG,EAAE,CAAC;QACP,OAAO,GAAG,IAAI,CAAC;IACnB,CAAC;SAAM,CAAC;QACJ,MAAM,CAAC,GAAG,MAAM,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;QAC1B,MAAM,CAAC,GAAG,MAAM,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;QAE1B,wEAAwE;QACxE,oEAAoE;QACpE,+DAA+D;QAC/D,OAAO,GAAG,GAAG,CAAC,CAAC,EAAC,CAAC,CAAC,IAAI,CAAC,IAAI,GAAG,CAAC,CAAC,EAAC,CAAC,CAAC,IAAI,CAAC,CAAC;IAC7C,CAAC;IAED,IAAI,CAAC,OAAO,EAAE,CAAC;QACX,OAAO,CAAC,CAAC,EAAE,EAAE,GAAI,EAAE,EAAE,CAAC,CAAC,CAAC;IAC5B,CAAC;IAED,6DAA6D;IAC7D,MAAM,IAAI,GAAG,GAAG,CAAC,KAAK,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,GAAC,GAAG,CAAC,CAAC;IAC1C,MAAM,IAAI,GAAG,GAAG,CAAC,KAAK,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,GAAC,GAAG,CAAC,CAAC;IAE1C,IAAI,IAAI,IAAI,sBAAsB,IAAI,IAAI,IAAI,sBAAsB,EAAE,CAAC;QACnE,OAAO,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,CAAE,qCAAqC;IAC7D,CAAC;IAED,0EAA0E;IAC1E,qEAAqE;IACrE,wEAAwE;IACxE,uEAAuE;IACvE,uEAAuE;IACvE,yEAAyE;IACzE,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,GAAC,CAAC,CAAC,GAAC,GAAG,CAAC,CAAC;IACjD,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,GAAC,CAAC,CAAC,GAAC,GAAG,CAAC,CAAC;IACjD,MAAM,EAAE,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,GAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,GAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAChD,MAAM,EAAE,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,GAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,GAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAChD,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,GAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,GAAC,CAAC,CAAC,CAAC;IAEjD,OAAO,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;AACtC,CAAC;AAGD,OAAO,EAAE,sBAAsB,EAAE,CAAA"}

package/node/vertex/remove-vertex.d.ts

@@ -1,10 +1,12 @@
 import type { CpNode } from "../cp-node/cp-node.js";
 import type { MatMeta } from '../mat/mat-meta.js';
 /**
- * Removes a cpNode from the MAT.
+ * Removes the (2-prong) vertex represented by the given `CpNode` from the MAT,
+ * i.e. both `CpNode`s on its circle are unlinked from the boundary and removed
+ * from `meta.cpTrees`.
  *
- * @param cpTree The tree graph holding the `CpNodes` of the MAT.
- * @param cpNode The `CpNode` to remove.
+ * @param cpNode a `CpNode` on the vertex to remove
+ * @param meta the `MatMeta` whose `cpTrees` hold the `CpNode`s of the MAT
  */
 declare function removeVertex(cpNode: CpNode, meta: MatMeta): void;
 export { removeVertex };

package/node/vertex/remove-vertex.js

@@ -2,10 +2,12 @@ import { getAllOnCircle } from '../cp-node/fs/get-all-on-circle.js';
 import { getProngCount } from '../cp-node/fs/get-prong-count.js';
 import { removeCpNode } from '../cp-node/fs/remove-cp-node.js';
 /**
- * Removes a cpNode from the MAT.
+ * Removes the (2-prong) vertex represented by the given `CpNode` from the MAT,
+ * i.e. both `CpNode`s on its circle are unlinked from the boundary and removed
+ * from `meta.cpTrees`.
  *
- * @param cpTree The tree graph holding the `CpNodes` of the MAT.
- * @param cpNode The `CpNode` to remove.
+ * @param cpNode a `CpNode` on the vertex to remove
+ * @param meta the `MatMeta` whose `cpTrees` hold the `CpNode`s of the MAT
  */
 function removeVertex(cpNode, meta) {
     const prongCount = getProngCount(cpNode);

package/node/vertex/remove-vertex.js.map

@@ -1 +1 @@
-{"version":3,"file":"remove-vertex.js","sourceRoot":"","sources":["../../src/vertex/remove-vertex.ts"],"names":[],"mappings":"AAEA,OAAO,EAAE,cAAc,EAAE,MAAM,oCAAoC,CAAC;AACpE,OAAO,EAAE,aAAa,EAAE,MAAM,kCAAkC,CAAC;AACjE,OAAO,EAAE,YAAY,EAAE,MAAM,iCAAiC,CAAC;AAG/D;;;;;GAKG;AACH,SAAS,YAAY,CACb,MAAc,EACd,IAAa;IAEjB,MAAM,UAAU,GAAG,aAAa,CAAC,MAAM,CAAC,CAAC;IACzC,IAAI,UAAU,KAAK,CAAC,EAAE,CAAC;QACnB,MAAM,IAAI,KAAK,CAAC,+BAA+B,UAAU,QAAQ,CAAC,CAAA;IACtE,CAAC;IAED,cAAc,CAAC,MAAM,CAAC,CAAC,OAAO,CAAC,MAAM,CAAC,EAAE,CAAC,YAAY,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC,CAAC;AACzE,CAAC;AAGD,OAAO,EAAE,YAAY,EAAE,CAAA"}
+{"version":3,"file":"remove-vertex.js","sourceRoot":"","sources":["../../src/vertex/remove-vertex.ts"],"names":[],"mappings":"AAEA,OAAO,EAAE,cAAc,EAAE,MAAM,oCAAoC,CAAC;AACpE,OAAO,EAAE,aAAa,EAAE,MAAM,kCAAkC,CAAC;AACjE,OAAO,EAAE,YAAY,EAAE,MAAM,iCAAiC,CAAC;AAG/D;;;;;;;GAOG;AACH,SAAS,YAAY,CACb,MAAc,EACd,IAAa;IAEjB,MAAM,UAAU,GAAG,aAAa,CAAC,MAAM,CAAC,CAAC;IACzC,IAAI,UAAU,KAAK,CAAC,EAAE,CAAC;QACnB,MAAM,IAAI,KAAK,CAAC,+BAA+B,UAAU,QAAQ,CAAC,CAAA;IACtE,CAAC;IAED,cAAc,CAAC,MAAM,CAAC,CAAC,OAAO,CAAC,MAAM,CAAC,EAAE,CAAC,YAAY,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC,CAAC;AACzE,CAAC;AAGD,OAAO,EAAE,YAAY,EAAE,CAAA"}