flo-mat 4.0.2 → 4.1.0

package/node/find-mat/get-meta.js

@@ -1,25 +1,8 @@
-// import { TriMapFs } from '../utils/tri-map.js';
 import { getBoundingHull, getBoundingBoxTight } from 'flo-bezier3';
 import { getBoundingBox$ } from '../geometry/get-bounding-box-.js';
 import { getCorner } from '../corner/get-corner.js';
 let timingStart;
 /** @internal */
-function getPointToCpNode(loops, cpTrees) {
-    timingStart = performance.now();
-    // const pointToCpNode: TriMap<Loop,number,number,CpNode> = new Map();
-    // for (const loop of loops) {
-    //     // Populate `posMap`
-    //     const cpTree = cpTrees.get(loop)!;
-    //     const cpNodes = cpTree.toArrayInOrder();
-    //     for (const cpNode of cpNodes) {
-    //         const { p, t } = cpNode.pointOnShape;
-    //         TriMapFs.set(pointToCpNode,loop,p[0],p[1],cpNode);
-    //     }
-    // }
-    // return pointToCpNode;
-    return undefined;
-}
-/** @internal */
 function getPartialMeta(loops) {
     const looseBoundingBoxes = [];
     const tightBoundingBoxes = [];
@@ -52,33 +35,36 @@ function getPartialMeta(loops) {
         dullCorners
     };
 }
+function addDebugInfo1() {
+    if (typeof _debug_ === 'undefined') {
+        return;
+    }
+    timingStart = performance.now();
+}
 /** @internal */
 function addDebugInfo2() {
     if (typeof _debug_ === 'undefined') {
         return;
     }
-    const timing = _debug_.generated.timing;
     const now = performance.now();
-    timing.holeClosers += now - timingStart;
+    _debug_.timing.holeClosers += now - timingStart;
     timingStart = now;
 }
+/** @internal */
 function addDebugInfo3() {
     if (typeof _debug_ === 'undefined') {
         return;
     }
-    const generated = _debug_.generated;
-    const timing = generated.timing;
     const now = performance.now();
-    timing.oneAnd2Prongs += now - timingStart;
+    _debug_.timing.oneAnd2Prongs += now - timingStart;
     timingStart = now;
 }
+/** @internal */
 function addDebugInfo4() {
     if (typeof _debug_ === 'undefined') {
         return;
     }
-    const generated = _debug_.generated;
-    const timing = generated.timing;
-    timing.threeProngs += performance.now() - timingStart;
+    _debug_.timing.threeProngs += performance.now() - timingStart;
 }
-export { getPointToCpNode, getPartialMeta, addDebugInfo2, addDebugInfo3, addDebugInfo4 };
+export { getPartialMeta, addDebugInfo1, addDebugInfo2, addDebugInfo3, addDebugInfo4 };
 //# sourceMappingURL=get-meta.js.map

package/node/find-mat/get-meta.js.map

@@ -1 +1 @@
-{"version":3,"file":"get-meta.js","sourceRoot":"","sources":["../../src/find-mat/get-meta.ts"],"names":[],"mappings":"AAMA,kDAAkD;AAClD,OAAO,EAAE,eAAe,EAAE,mBAAmB,EAAE,MAAM,aAAa,CAAC;AACnE,OAAO,EAAE,eAAe,EAAE,MAAM,kCAAkC,CAAC;AAEnE,OAAO,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AAOpD,IAAI,WAAmB,CAAC;AAGxB,gBAAgB;AAChB,SAAS,gBAAgB,CACjB,KAAa,EACb,OAAkC;IAEtC,WAAW,GAAG,WAAW,CAAC,GAAG,EAAE,CAAC;IAEhC,sEAAsE;IAEtE,8BAA8B;IAC9B,2BAA2B;IAC3B,yCAAyC;IACzC,+CAA+C;IAC/C,sCAAsC;IACtC,gDAAgD;IAChD,6DAA6D;IAC7D,QAAQ;IACR,IAAI;IAEJ,wBAAwB;IACxB,OAAO,SAAU,CAAA;AACrB,CAAC;AAGD,gBAAgB;AAChB,SAAS,cAAc,CACf,KAAa;IAUjB,MAAM,kBAAkB,GAAiB,EAAE,CAAC;IAC5C,MAAM,kBAAkB,GAAiB,EAAE,CAAC;IAC5C,MAAM,aAAa,GAAiB,EAAE,CAAC;IACvC,MAAM,YAAY,GAAY,EAAE,CAAC;IACjC,MAAM,WAAW,GAAY,EAAE,CAAC;IAEhC,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACvB,KAAK,MAAM,KAAK,IAAI,IAAI,CAAC,MAAM,EAAE,CAAC;YAC9B,MAAM,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,MAAM,IAAI,GAAG,eAAe,CAAC,EAAE,EAAE,IAAI,CAAE,CAAC;YAExC,aAAa,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;YAEzB,MAAM,gBAAgB,GAAG,eAAe,CAAC,EAAE,CAAC,CAAC;YAC7C,kBAAkB,CAAC,IAAI,CAAC,gBAAgB,CAAC,CAAC;YAE1C,MAAM,gBAAgB,GAAG,mBAAmB,CAAC,EAAE,CAAC,CAAC;YACjD,kBAAkB,CAAC,IAAI,CAAC,gBAAgB,CAAC,CAAC;YAE1C,MAAM,MAAM,GAAG,SAAS,CAAC,KAAK,CAAC,EAAE,EAAE,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;YAClD,IAAI,MAAM,CAAC,YAAY,EAAE,CAAC;gBACtB,YAAY,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;YAC7B,CAAC;iBAAM,IAAI,MAAM,CAAC,WAAW,EAAE,CAAC;gBAC5B,WAAW,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;YAC5B,CAAC;QACL,CAAC;IACL,CAAC;IAED,OAAO;QACH,kBAAkB;QAClB,kBAAkB;QAClB,aAAa;QACb,YAAY;QACZ,WAAW;KACd,CAAA;AACL,CAAC;AAGD,gBAAgB;AAChB,SAAS,aAAa;IAClB,IAAI,OAAO,OAAO,KAAK,WAAW,EAAE,CAAC;QAAC,OAAO;IAAC,CAAC;IAE/C,MAAM,MAAM,GAAG,OAAO,CAAC,SAAS,CAAC,MAAM,CAAC;IACxC,MAAM,GAAG,GAAG,WAAW,CAAC,GAAG,EAAE,CAAC;IAC9B,MAAM,CAAC,WAAW,IAAI,GAAG,GAAG,WAAW,CAAC;IACxC,WAAW,GAAG,GAAG,CAAC;AACtB,CAAC;AAGD,SAAS,aAAa;IAClB,IAAI,OAAO,OAAO,KAAK,WAAW,EAAE,CAAC;QAAC,OAAO;IAAC,CAAC;IAE/C,MAAM,SAAS,GAAG,OAAO,CAAC,SAAS,CAAC;IACpC,MAAM,MAAM,GAAG,SAAS,CAAC,MAAM,CAAC;IAChC,MAAM,GAAG,GAAG,WAAW,CAAC,GAAG,EAAE,CAAC;IAE9B,MAAM,CAAC,aAAa,IAAI,GAAG,GAAG,WAAW,CAAC;IAC1C,WAAW,GAAG,GAAG,CAAC;AACtB,CAAC;AAGD,SAAS,aAAa;IAClB,IAAI,OAAO,OAAO,KAAK,WAAW,EAAE,CAAC;QAAC,OAAO;IAAC,CAAC;IAE/C,MAAM,SAAS,GAAG,OAAO,CAAC,SAAS,CAAC;IACpC,MAAM,MAAM,GAAG,SAAS,CAAC,MAAM,CAAC;IAEhC,MAAM,CAAC,WAAW,IAAI,WAAW,CAAC,GAAG,EAAE,GAAG,WAAW,CAAC;AAC1D,CAAC;AAGD,OAAO,EACH,gBAAgB,EAChB,cAAc,EACd,aAAa,EACb,aAAa,EACb,aAAa,EAChB,CAAA"}
+{"version":3,"file":"get-meta.js","sourceRoot":"","sources":["../../src/find-mat/get-meta.ts"],"names":[],"mappings":"AAIA,OAAO,EAAE,eAAe,EAAE,mBAAmB,EAAE,MAAM,aAAa,CAAC;AACnE,OAAO,EAAE,eAAe,EAAE,MAAM,kCAAkC,CAAC;AACnE,OAAO,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AAOpD,IAAI,WAAmB,CAAC;AAGxB,gBAAgB;AAChB,SAAS,cAAc,CACf,KAAa;IAUjB,MAAM,kBAAkB,GAAiB,EAAE,CAAC;IAC5C,MAAM,kBAAkB,GAAiB,EAAE,CAAC;IAC5C,MAAM,aAAa,GAAiB,EAAE,CAAC;IACvC,MAAM,YAAY,GAAY,EAAE,CAAC;IACjC,MAAM,WAAW,GAAY,EAAE,CAAC;IAEhC,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACvB,KAAK,MAAM,KAAK,IAAI,IAAI,CAAC,MAAM,EAAE,CAAC;YAC9B,MAAM,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC;YACpB,MAAM,IAAI,GAAG,eAAe,CAAC,EAAE,EAAE,IAAI,CAAE,CAAC;YAExC,aAAa,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;YAEzB,MAAM,gBAAgB,GAAG,eAAe,CAAC,EAAE,CAAC,CAAC;YAC7C,kBAAkB,CAAC,IAAI,CAAC,gBAAgB,CAAC,CAAC;YAE1C,MAAM,gBAAgB,GAAG,mBAAmB,CAAC,EAAE,CAAC,CAAC;YACjD,kBAAkB,CAAC,IAAI,CAAC,gBAAgB,CAAC,CAAC;YAE1C,MAAM,MAAM,GAAG,SAAS,CAAC,KAAK,CAAC,EAAE,EAAE,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;YAClD,IAAI,MAAM,CAAC,YAAY,EAAE,CAAC;gBACtB,YAAY,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;YAC7B,CAAC;iBAAM,IAAI,MAAM,CAAC,WAAW,EAAE,CAAC;gBAC5B,WAAW,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;YAC5B,CAAC;QACL,CAAC;IACL,CAAC;IAED,OAAO;QACH,kBAAkB;QAClB,kBAAkB;QAClB,aAAa;QACb,YAAY;QACZ,WAAW;KACd,CAAA;AACL,CAAC;AAGD,SAAS,aAAa;IAClB,IAAI,OAAO,OAAO,KAAK,WAAW,EAAE,CAAC;QAAC,OAAO;IAAC,CAAC;IAE/C,WAAW,GAAG,WAAW,CAAC,GAAG,EAAE,CAAC;AACpC,CAAC;AAID,gBAAgB;AAChB,SAAS,aAAa;IAClB,IAAI,OAAO,OAAO,KAAK,WAAW,EAAE,CAAC;QAAC,OAAO;IAAC,CAAC;IAE/C,MAAM,GAAG,GAAG,WAAW,CAAC,GAAG,EAAE,CAAC;IAC9B,OAAO,CAAC,MAAM,CAAC,WAAW,IAAI,GAAG,GAAG,WAAW,CAAC;IAChD,WAAW,GAAG,GAAG,CAAC;AACtB,CAAC;AAGD,gBAAgB;AAChB,SAAS,aAAa;IAClB,IAAI,OAAO,OAAO,KAAK,WAAW,EAAE,CAAC;QAAC,OAAO;IAAC,CAAC;IAE/C,MAAM,GAAG,GAAG,WAAW,CAAC,GAAG,EAAE,CAAC;IAE9B,OAAO,CAAC,MAAM,CAAC,aAAa,IAAI,GAAG,GAAG,WAAW,CAAC;IAClD,WAAW,GAAG,GAAG,CAAC;AACtB,CAAC;AAGD,gBAAgB;AAChB,SAAS,aAAa;IAClB,IAAI,OAAO,OAAO,KAAK,WAAW,EAAE,CAAC;QAAC,OAAO;IAAC,CAAC;IAE/C,OAAO,CAAC,MAAM,CAAC,WAAW,IAAI,WAAW,CAAC,GAAG,EAAE,GAAG,WAAW,CAAC;AAClE,CAAC;AAGD,OAAO,EACH,cAAc,EACd,aAAa,EACb,aAAa,EACb,aAAa,EACb,aAAa,EAChB,CAAA"}

package/node/find-mat/get-sharp-corners.js

@@ -1,5 +1,4 @@
 import { controlPointLinesLength } from 'flo-bezier3';
-import { toP } from '../point-on-shape/to-p.js';
 import { getCorner } from '../corner/get-corner.js';
 /** @internal */
 function getSharpCornersOnLoop(minBezLength) {
@@ -11,7 +10,10 @@ function getSharpCornersOnLoop(minBezLength) {
                 continue;
             }
             if (getCorner(curve.ps, curve.next.ps).isQuiteSharp) {
-                sharpCorners.push({ curve, t: 1, isSource: true, p: toP(curve.ps, 1) });
+                const { ps } = curve;
+                const p = ps[ps.length - 1];
+                const pDd = [[0, p[0]], [0, p[1]]];
+                sharpCorners.push({ curve, t: 1, isSource: true, p, pDd });
             }
         }
         return sharpCorners;

package/node/find-mat/get-sharp-corners.js.map

@@ -1 +1 @@
-{"version":3,"file":"get-sharp-corners.js","sourceRoot":"","sources":["../../src/find-mat/get-sharp-corners.ts"],"names":[],"mappings":"AAEA,OAAO,EAAE,uBAAuB,EAAE,MAAM,aAAa,CAAC;AACtD,OAAO,EAAE,GAAG,EAAE,MAAM,2BAA2B,CAAC;AAChD,OAAO,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AAGpD,gBAAgB;AAChB,SAAS,qBAAqB,CACtB,YAAoB;IAExB,OAAO,CAAC,IAAU,EAAqB,EAAE;QACrC,MAAM,YAAY,GAAsB,EAAE,CAAC;QAE3C,KAAK,IAAI,CAAC,GAAC,CAAC,EAAE,CAAC,GAAC,IAAI,CAAC,MAAM,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YACtC,MAAM,KAAK,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;YAE7B,IAAI,uBAAuB,CAAC,KAAK,CAAC,EAAE,CAAC,GAAG,YAAY,EAAE,CAAC;gBACnD,SAAS;YACb,CAAC;YAED,IAAI,SAAS,CAAC,KAAK,CAAC,EAAE,EAAE,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,YAAY,EAAE,CAAC;gBAClD,YAAY,CAAC,IAAI,CAAC,EAAE,KAAK,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,EAAE,IAAI,EAAE,CAAC,EAAE,GAAG,CAAC,KAAK,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;YAC5E,CAAC;QACL,CAAC;QAED,OAAO,YAAY,CAAC;IACxB,CAAC,CAAA;AACL,CAAC;AAID,OAAO,EAAE,qBAAqB,EAAE,CAAA"}
+{"version":3,"file":"get-sharp-corners.js","sourceRoot":"","sources":["../../src/find-mat/get-sharp-corners.ts"],"names":[],"mappings":"AAEA,OAAO,EAAE,uBAAuB,EAAE,MAAM,aAAa,CAAC;AAEtD,OAAO,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AAGpD,gBAAgB;AAChB,SAAS,qBAAqB,CACtB,YAAoB;IAExB,OAAO,CAAC,IAAU,EAAqB,EAAE;QACrC,MAAM,YAAY,GAAsB,EAAE,CAAC;QAE3C,KAAK,IAAI,CAAC,GAAC,CAAC,EAAE,CAAC,GAAC,IAAI,CAAC,MAAM,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YACtC,MAAM,KAAK,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;YAE7B,IAAI,uBAAuB,CAAC,KAAK,CAAC,EAAE,CAAC,GAAG,YAAY,EAAE,CAAC;gBACnD,SAAS;YACb,CAAC;YAED,IAAI,SAAS,CAAC,KAAK,CAAC,EAAE,EAAE,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,YAAY,EAAE,CAAC;gBAClD,MAAM,EAAE,EAAE,EAAE,GAAG,KAAK,CAAC;gBAErB,MAAM,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;gBAC5B,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,EAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAC,CAAC,CAAC,EAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAChC,YAAY,CAAC,IAAI,CAAC,EAAE,KAAK,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,EAAE,IAAI,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;YAC/D,CAAC;QACL,CAAC;QAED,OAAO,YAAY,CAAC;IACxB,CAAC,CAAA;AACL,CAAC;AAID,OAAO,EAAE,qBAAqB,EAAE,CAAA"}

package/node/get-medial-points/double-double/dd-get-medial-point-coeffs-bez0.d.ts

@@ -0,0 +1,21 @@
+/**
+ * Returns the coefficients `a0` and `b0` for the linear equation
+ * `a0⋅t + b0 = 0`, so the ray parameter can be recovered as `t = -b0/a0`.
+ * That parameter value makes `q(t) = p + t⋅v` equidistant from `p` and `P`.
+ *
+ * Let `p` be a fixed point in the plane.\
+ * Let `v` be a direction vector defining the ray `q(t) = p + t⋅v`.\
+ * Let `P` be another point.\
+ *
+ * In other words, this function returns the coefficients needed to solve for
+ * the ray parameter of the medial point on the ray.
+ *
+ * @param p base point
+ * @param v ray direction from `p`
+ * @param P another "target" point, e.g. `[1,2]`
+ */
+declare function ddGetMedialPointCoeffsBez0(p: number[][], v: number[][], P: number[]): {
+    a0: number[];
+    b0: number[];
+};
+export { ddGetMedialPointCoeffsBez0 };

package/node/get-medial-points/double-double/dd-get-medial-point-coeffs-bez0.js

@@ -0,0 +1,46 @@
+import { ddAddDd, ddMultDd, ddMultBy2, ddDiffDouble } from 'double-double';
+const qaq = ddAddDd;
+const qmq = ddMultDd;
+const qm2 = ddMultBy2;
+const qdd = ddDiffDouble;
+// TODO - move all these upstream to flo-bezier3
+/**
+ * Returns the coefficients `a0` and `b0` for the linear equation
+ * `a0⋅t + b0 = 0`, so the ray parameter can be recovered as `t = -b0/a0`.
+ * That parameter value makes `q(t) = p + t⋅v` equidistant from `p` and `P`.
+ *
+ * Let `p` be a fixed point in the plane.\
+ * Let `v` be a direction vector defining the ray `q(t) = p + t⋅v`.\
+ * Let `P` be another point.\
+ *
+ * In other words, this function returns the coefficients needed to solve for
+ * the ray parameter of the medial point on the ray.
+ *
+ * @param p base point
+ * @param v ray direction from `p`
+ * @param P another "target" point, e.g. `[1,2]`
+ */
+function ddGetMedialPointCoeffsBez0(p, v, P) {
+    // -----------------------------------------------------
+    // See get-medial-points.md for implementation details.
+    // -----------------------------------------------------
+    const [px, py] = p;
+    const [vx, vy] = v;
+    const [x0, y0] = P;
+    // Constant bezier in power basis: b(s) = c
+    // u(s) = p - b(s) = u0
+    // const u0x = px - x0;
+    // const u0y = py - y0;
+    const u0x = qdd(px, x0);
+    const u0y = qdd(py, y0);
+    // E2(s,t): |t⋅v|² - |u(s) + t⋅v|² = 0
+    //         => 2⋅(v⋅u(s))⋅t + |u(s)|² = 0
+    //         => A(s)⋅t + B(s) = 0
+    // const a0 = 2*(vx*u0x + vy*u0y);
+    // const b0 = u0x*u0x + u0y*u0y;
+    const a0 = qaq(qmq(qm2(vx), u0x), qmq(qm2(vy), u0y));
+    const b0 = qaq(qmq(u0x, u0x), qmq(u0y, u0y));
+    return { a0, b0 };
+}
+export { ddGetMedialPointCoeffsBez0 };
+//# sourceMappingURL=dd-get-medial-point-coeffs-bez0.js.map

package/node/get-medial-points/double-double/dd-get-medial-point-coeffs-bez0.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"dd-get-medial-point-coeffs-bez0.js","sourceRoot":"","sources":["../../../src/get-medial-points/double-double/dd-get-medial-point-coeffs-bez0.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,SAAS,EAAE,YAAY,EAAE,MAAM,eAAe,CAAC;AAE3E,MAAM,GAAG,GAAG,OAAO,CAAC;AACpB,MAAM,GAAG,GAAG,QAAQ,CAAC;AACrB,MAAM,GAAG,GAAG,SAAS,CAAC;AACtB,MAAM,GAAG,GAAG,YAAY,CAAC;AAGzB,gDAAgD;AAEhD;;;;;;;;;;;;;;;GAeG;AACH,SAAS,0BAA0B,CAC3B,CAAa,EACb,CAAa,EACb,CAAW;IAEf,wDAAwD;IACxD,uDAAuD;IACvD,wDAAwD;IAExD,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC;IACnB,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC;IACnB,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC;IAEnB,2CAA2C;IAC3C,uBAAuB;IAEvB,uBAAuB;IACvB,uBAAuB;IACvB,MAAM,GAAG,GAAG,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;IACvB,MAAM,GAAG,GAAG,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;IAEvB,sCAAsC;IACtC,wCAAwC;IACxC,+BAA+B;IAE/B,kCAAkC;IAClC,gCAAgC;IAChC,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,EAAC,GAAG,CAAC,EAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,EAAC,GAAG,CAAC,CAAC,CAAC;IAClD,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,EAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,EAAC,GAAG,CAAC,CAAC,CAAC;IAE3C,OAAO,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;AACtB,CAAC;AAGD,OAAO,EAAE,0BAA0B,EAAE,CAAA"}

package/node/get-medial-points/double-double/dd-get-medial-point-coeffs-bez1.d.ts

@@ -0,0 +1,28 @@
+/**
+ * Returns the polynomial coefficients used to recover ray parameter values `t`
+ * for a ray `q(t) = p + t⋅v` and an order-1 bezier curve `ps`.
+ *
+ * The returned polynomials encode the medial-condition equations in `s` and
+ * `t`. For any valid solution, `t` is recovered by eliminating `s` and solving
+ * the resulting ray equation, or equivalently by using the linear form
+ * `A(s)⋅t + B(s) = 0`, which gives `t = -B(s)/A(s)`.
+ *
+ * More specifically, the returned values represent:
+ * * `A` and `B`: the coefficients of `E2(s,t) = A(s)⋅t + B(s)`
+ * * `C` and `D`: the coefficients of `E1(s,t) = C(s)⋅t + D(s)`
+ * * `H`: the eliminated polynomial `A(s)⋅D(s) - B(s)⋅C(s)` whose roots are
+ *   candidate `s` values for medial points
+ *
+ * @param p base point
+ * @param v ray direction vector starting from `p`
+ * @param ps order 1 bezier control points, i.e. a line segment
+ * given as an array of control points, e.g. `[[0,0],[2,1]]`
+ */
+declare function ddGetMedialPointCoeffsBez1(p: number[][], v: number[][], ps: number[][]): {
+    A: number[][];
+    B: number[][];
+    C: number[][];
+    D: number[][];
+    H: number[][];
+};
+export { ddGetMedialPointCoeffsBez1 };

package/node/get-medial-points/double-double/dd-get-medial-point-coeffs-bez1.js

@@ -0,0 +1,97 @@
+import { twoDiff, ddAddDd, ddMultDd, ddNegativeOf, ddMultByNeg2, ddDiffDouble, ddMultBy2 } from 'double-double';
+const td = twoDiff;
+const qaq = ddAddDd;
+const qmq = ddMultDd;
+const qno = ddNegativeOf;
+const qmn2 = ddMultByNeg2;
+const qm2 = ddMultBy2;
+const qdd = ddDiffDouble;
+/**
+ * Returns the polynomial coefficients used to recover ray parameter values `t`
+ * for a ray `q(t) = p + t⋅v` and an order-1 bezier curve `ps`.
+ *
+ * The returned polynomials encode the medial-condition equations in `s` and
+ * `t`. For any valid solution, `t` is recovered by eliminating `s` and solving
+ * the resulting ray equation, or equivalently by using the linear form
+ * `A(s)⋅t + B(s) = 0`, which gives `t = -B(s)/A(s)`.
+ *
+ * More specifically, the returned values represent:
+ * * `A` and `B`: the coefficients of `E2(s,t) = A(s)⋅t + B(s)`
+ * * `C` and `D`: the coefficients of `E1(s,t) = C(s)⋅t + D(s)`
+ * * `H`: the eliminated polynomial `A(s)⋅D(s) - B(s)⋅C(s)` whose roots are
+ *   candidate `s` values for medial points
+ *
+ * @param p base point
+ * @param v ray direction vector starting from `p`
+ * @param ps order 1 bezier control points, i.e. a line segment
+ * given as an array of control points, e.g. `[[0,0],[2,1]]`
+ */
+function ddGetMedialPointCoeffsBez1(p, v, ps) {
+    // -----------------------------------------------------
+    // See get-medial-points.md for implementation details.
+    // -----------------------------------------------------
+    const [px, py] = p;
+    const [vx, vy] = v;
+    const [[x0, y0], [x1, y1]] = ps;
+    // Linear bezier in power basis: b(s) = b*s + c
+    // We can also use: `const [[bx,by], [cx,cy]] = toPowerBasis1(ps)`
+    // const bx = x1 - x0;
+    // const by = y1 - y0;
+    const bx = td(x1, x0);
+    const by = td(y1, y0);
+    // u(s) = p - b(s) = u1*s + u0
+    // const u0x = px - x0;
+    // const u0y = py - y0;
+    const u0x = qdd(px, x0);
+    const u0y = qdd(py, y0);
+    // b'(s) = w(s) = w0
+    // -----------------------------------------------------
+    // E1(s,t): (u(s) + t⋅v) ⋅ b'(s) = 0
+    // => C(s)⋅t + D(s) = 0
+    // const c0 = vx*bx + vy*by;
+    const c0 = qaq(qmq(vx, bx), qmq(vy, by));
+    // -----------------------------------------------------
+    // -----------------------------------------------------
+    // const d1 = -(bx*bx + by*by);
+    // const d0 = bx*u0x + by*u0y;
+    // const d1 = qmd(-1,qaq(qmq(bx,bx),qmq(by,by)));
+    const d1 = qno(qaq(qmq(bx, bx), qmq(by, by)));
+    const d0 = qaq(qmq(bx, u0x), qmq(by, u0y));
+    // -----------------------------------------------------
+    // -----------------------------------------------------
+    // E2(s,t): |t⋅v|² - |u(s) + t⋅v|² = 0
+    //         => 2⋅(v⋅u(s))⋅t + |u(s)|² = 0
+    //         => A(s)⋅t + B(s) = 0
+    // const a1 = -2*c0;
+    // const a0 = 2*(vx*u0x + vy*u0y);
+    const a1 = qmn2(c0);
+    const a0 = qaq(qmq(qm2(vx), u0x), qmq(qm2(vy), u0y));
+    // -----------------------------------------------------
+    // -----------------------------------------------------
+    // const b2 = -d1;
+    // const b1 = -2*d0;
+    // const b0 = u0x*u0x + u0y*u0y;
+    const b2 = qno(d1);
+    const b1 = qmn2(d0);
+    const b0 = qaq(qmq(u0x, u0x), qmq(u0y, u0y));
+    // -----------------------------------------------------
+    // Eliminate t from:
+    //   A(s)⋅t + B(s) = 0
+    //   C(s)⋅t + D(s) = 0
+    // by taking A(s)⋅D(s) - B(s)⋅C(s) = 0 (degree ≤ 2 in s)
+    // const H2 = b2*c0;
+    // const H1 = a0*d1;
+    // const H0 = a0*d0 - b0*c0;
+    const H2 = qmq(b2, c0);
+    const H1 = qmq(a0, d1);
+    const H0 = qaq(qmq(a0, d0), qno(qmq(b0, c0)));
+    return {
+        A: [a1, a0],
+        B: [b2, b1, b0],
+        C: [c0],
+        D: [d1, d0],
+        H: [H2, H1, H0]
+    };
+}
+export { ddGetMedialPointCoeffsBez1 };
+//# sourceMappingURL=dd-get-medial-point-coeffs-bez1.js.map

package/node/get-medial-points/double-double/dd-get-medial-point-coeffs-bez1.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"dd-get-medial-point-coeffs-bez1.js","sourceRoot":"","sources":["../../../src/get-medial-points/double-double/dd-get-medial-point-coeffs-bez1.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,YAAY,EAAE,YAAY,EAAE,YAAY,EAAE,SAAS,EAAE,MAAM,eAAe,CAAC;AAEhH,MAAM,EAAE,GAAG,OAAO,CAAC;AACnB,MAAM,GAAG,GAAG,OAAO,CAAC;AACpB,MAAM,GAAG,GAAG,QAAQ,CAAC;AACrB,MAAM,GAAG,GAAG,YAAY,CAAC;AACzB,MAAM,IAAI,GAAG,YAAY,CAAC;AAC1B,MAAM,GAAG,GAAG,SAAS,CAAC;AACtB,MAAM,GAAG,GAAG,YAAY,CAAC;AAGzB;;;;;;;;;;;;;;;;;;;GAmBG;AACH,SAAS,0BAA0B,CAC3B,CAAa,EACb,CAAa,EACb,EAAc;IAElB,wDAAwD;IACxD,uDAAuD;IACvD,wDAAwD;IACxD,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC;IACnB,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC;IACnB,MAAM,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC;IAEhC,+CAA+C;IAC/C,kEAAkE;IAClE,sBAAsB;IACtB,sBAAsB;IACtB,MAAM,EAAE,GAAG,EAAE,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;IACrB,MAAM,EAAE,GAAG,EAAE,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;IAErB,8BAA8B;IAC9B,uBAAuB;IACvB,uBAAuB;IACvB,MAAM,GAAG,GAAG,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;IACvB,MAAM,GAAG,GAAG,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;IAEvB,oBAAoB;IAEpB,wDAAwD;IACxD,oCAAoC;IACpC,uBAAuB;IACvB,4BAA4B;IAC5B,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,EAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC,CAAC;IACtC,wDAAwD;IAExD,wDAAwD;IACxD,+BAA+B;IAC/B,8BAA8B;IAC9B,iDAAiD;IACjD,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,EAAE,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAC5C,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,GAAG,CAAC,EAAE,GAAG,CAAC,EAAE,EAAC,GAAG,CAAC,CAAC,CAAC;IACzC,wDAAwD;IAExD,wDAAwD;IACxD,sCAAsC;IACtC,wCAAwC;IACxC,+BAA+B;IAC/B,oBAAoB;IACpB,kCAAkC;IAClC,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;IACpB,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,EAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,EAAC,GAAG,CAAC,CAAC,CAAC;IACnD,wDAAwD;IAExD,wDAAwD;IACxD,kBAAkB;IAClB,oBAAoB;IACpB,gCAAgC;IAChC,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC;IACnB,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;IACpB,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,EAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,EAAC,GAAG,CAAC,CAAC,CAAC;IAC3C,wDAAwD;IAGxD,oBAAoB;IACpB,sBAAsB;IACtB,sBAAsB;IACtB,wDAAwD;IAExD,oBAAoB;IACpB,oBAAoB;IACpB,4BAA4B;IAC5B,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;IACtB,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;IACtB,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAE5C,OAAO;QACH,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC;QACX,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;QACf,CAAC,EAAE,CAAC,EAAE,CAAC;QACP,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC;QACX,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;KAClB,CAAC;AACN,CAAC;AAGD,OAAO,EAAE,0BAA0B,EAAE,CAAA"}

package/node/get-medial-points/double-double/dd-get-medial-point-coeffs-bez2-same-curve.d.ts

@@ -0,0 +1,28 @@
+/**
+ * Returns polynomial coefficients for the same-curve quadratic medial-point
+ * case, where `p = b(t)` on the same quadratic bezier `ps` and `v` is normal
+ * to the curve at that parameter `t`.
+ *
+ * The returned coefficients encode:
+ * * `A`, `B`: `E2(s, τ) = (s - t)^2⋅(A(s)⋅τ + B(s))`
+//  * * `C`, `D`: `E1(s, τ) = C(s)⋅τ + D(s)`
+ * * `H`: the reduced eliminant in `s`
+ *
+ * Here `τ` is the ray parameter in `q(τ) = p + τ⋅v`.
+ * For a candidate root `s` of `H`, recover `τ` from
+ * `A(s)⋅τ + B(s) = 0`, i.e. `τ = -B(s)/A(s)` when `A(s) ≠ 0`.
+ *
+ * In this same-curve setup, the full eliminant has the form
+ * `H_full(s) = (s - t)^4*(l1*s + l0)`. This function returns only the reduced
+ * factor `H(s) = l1*s + l0` (the repeated `(s - t)^4` factor is omitted).
+ *
+ * @param t parameter on `ps` where `p = b(t)`
+ * @param v ray direction from `p`; assumed normal to `ps` at `t`
+ * @param ps order 2 bezier control points, e.g. `[[0,0],[1,1],[2,1]]`
+ */
+declare function ddGetMedialPointCoeffsBez2_SameCurve(t: number, v: number[][], ps: number[][]): {
+    A: number[][];
+    B: number[][];
+    H: number[][];
+};
+export { ddGetMedialPointCoeffsBez2_SameCurve };

package/node/get-medial-points/double-double/dd-get-medial-point-coeffs-bez2-same-curve.js

@@ -0,0 +1,109 @@
+import { toPowerBasisDd } from 'flo-bezier3';
+import { ddMultDouble2, ddAddDd, ddMultDd, ddNegativeOf, ddMultBy2 } from 'double-double';
+const qmd = ddMultDouble2;
+const qaq = ddAddDd;
+const qmq = ddMultDd;
+const qno = ddNegativeOf;
+const qm2 = ddMultBy2;
+/**
+ * Returns polynomial coefficients for the same-curve quadratic medial-point
+ * case, where `p = b(t)` on the same quadratic bezier `ps` and `v` is normal
+ * to the curve at that parameter `t`.
+ *
+ * The returned coefficients encode:
+ * * `A`, `B`: `E2(s, τ) = (s - t)^2⋅(A(s)⋅τ + B(s))`
+//  * * `C`, `D`: `E1(s, τ) = C(s)⋅τ + D(s)`
+ * * `H`: the reduced eliminant in `s`
+ *
+ * Here `τ` is the ray parameter in `q(τ) = p + τ⋅v`.
+ * For a candidate root `s` of `H`, recover `τ` from
+ * `A(s)⋅τ + B(s) = 0`, i.e. `τ = -B(s)/A(s)` when `A(s) ≠ 0`.
+ *
+ * In this same-curve setup, the full eliminant has the form
+ * `H_full(s) = (s - t)^4*(l1*s + l0)`. This function returns only the reduced
+ * factor `H(s) = l1*s + l0` (the repeated `(s - t)^4` factor is omitted).
+ *
+ * @param t parameter on `ps` where `p = b(t)`
+ * @param v ray direction from `p`; assumed normal to `ps` at `t`
+ * @param ps order 2 bezier control points, e.g. `[[0,0],[1,1],[2,1]]`
+ */
+function ddGetMedialPointCoeffsBez2_SameCurve(t, v, ps) {
+    // -----------------------------------------------------
+    // See get-medial-points.md for implementation details.
+    // -----------------------------------------------------
+    // Quadratic bezier in power basis: b(s) = a⋅s² + b⋅s + c
+    const [[ax, bx], [ay, by]] = toPowerBasisDd(ps);
+    const [vx, vy] = v;
+    // Same-curve assumption with explicit parameter `t`:
+    // p = b(t) => u0 = p - c = t*(a*t + b).
+    // const g0x = ax*t + bx;
+    // const g0y = ay*t + by;
+    const g0x = qaq(qmd(t, ax), bx);
+    const g0y = qaq(qmd(t, ay), by);
+    // Reuse core quadratic-form terms to reduce repeated multiplications.
+    // const br2 = ax*ax + ay*ay;
+    // const ab = ax*bx + ay*by;
+    // const bb = bx*bx + by*by;
+    // const ag = ax*g0x + ay*g0y;
+    const br2 = qaq(qmq(ax, ax), qmq(ay, ay));
+    const ab = qaq(qmq(ax, bx), qmq(ay, by));
+    const ag = qaq(qmq(ax, g0x), qmq(ay, g0y));
+    // -----------------------------------------------------
+    // E1(s,t): (u(s) + t⋅v) ⋅ b'(s) = 0
+    // => C(s)⋅t + D(s) = 0
+    // const c1 = 2*(vx*ax + vy*ay);
+    const c1 = qm2(qaq(qmq(vx, ax), qmq(vy, ay)));
+    // v is normal at parameter t, so C(t) = v⋅w(t) = 0.
+    // Since C(s) = c1*s + c0, enforce c0 = -c1*t exactly.
+    // const c0 = -c1*t;
+    // const d3 = -2*br2;
+    // const d2 = -3*ab;
+    // const d1 = 2*t*ag - bb;
+    // p = b(t) => D(t) = u(t)⋅w(t) = 0.
+    // const d0 = -t*(t*(d3*t + d2) + d1);
+    // -----------------------------------------------------
+    // -----------------------------------------------------
+    // E2(s,t): |t⋅v|² - |u(s) + t⋅v|² = 0
+    //         => 2⋅(v⋅u(s))⋅t + |u(s)|² = 0
+    //         => (s - t)^2⋅(A(s)⋅t + B(s)) = 0
+    // Return only the reduced constant A(s) = ar0.
+    // const ar0 = -c1;
+    const ar0 = qno(c1);
+    // Full B(s) factorizes as: B_full(s) = (s - t)^2*B_reduced(s).
+    // Return only B_reduced(s) = br2*s^2 + br1*s + br0.
+    // const br1 = 2*ag;
+    // const br0 = g0x*g0x + g0y*g0y;
+    const br1 = qm2(ag);
+    const br0 = qaq(qmq(g0x, g0x), qmq(g0y, g0y));
+    // -----------------------------------------------------
+    // Eliminate t from:
+    //   A(s)⋅t + B(s) = 0
+    //   C(s)⋅t + D(s) = 0
+    // by taking A(s)⋅D(s) - B(s)⋅C(s) = 0 (degree ≤ 5 in s)
+    // For this same-curve quadratic case:
+    //   d3 = -2*b4, d2 = -(3/2)*b3, a2 = -c1, a1 = -2*c0.
+    // Therefore H5 and H4 reduce to:
+    //   H5 = b4*c1
+    //   H4 = 3*b4*c0 + (1/2)*b3*c1
+    // const H5 = br2*c1;
+    const H5 = qmq(br2, c1);
+    // In the same-curve quadratic case the full eliminant has the form
+    // H_full(s) = (s - t)^4*(l1*s + l0).
+    // Match the leading coefficients of H_full(s):
+    //   H5 = l1
+    //   H4 = l0 - 4*t*l1
+    // const l1 = H5;
+    // const l0 = H4 + 4*t*l1;
+    // Using c0 = -c1*t and b3 = 2*ab, this simplifies to l0 = c1*(br2*t + ab).
+    // const l0 = c1*(br2*t + ab);
+    const l0 = qmq(c1, qaq(qmd(t, br2), ab));
+    return {
+        A: [ar0],
+        B: [br2, br1, br0],
+        // C: [c1, c0],
+        // D: [d3, d2, d1, d0],
+        H: [H5, l0]
+    };
+}
+export { ddGetMedialPointCoeffsBez2_SameCurve };
+//# sourceMappingURL=dd-get-medial-point-coeffs-bez2-same-curve.js.map

package/node/get-medial-points/double-double/dd-get-medial-point-coeffs-bez2-same-curve.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"dd-get-medial-point-coeffs-bez2-same-curve.js","sourceRoot":"","sources":["../../../src/get-medial-points/double-double/dd-get-medial-point-coeffs-bez2-same-curve.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,cAAc,EAAE,MAAM,aAAa,CAAC;AAC7C,OAAO,EAAE,aAAa,EAAE,OAAO,EAAE,QAAQ,EAAE,YAAY,EAAE,SAAS,EAAE,MAAM,eAAe,CAAC;AAE1F,MAAM,GAAG,GAAG,aAAa,CAAC;AAC1B,MAAM,GAAG,GAAG,OAAO,CAAC;AACpB,MAAM,GAAG,GAAG,QAAQ,CAAC;AACrB,MAAM,GAAG,GAAG,YAAY,CAAC;AACzB,MAAM,GAAG,GAAG,SAAS,CAAC;AAGtB;;;;;;;;;;;;;;;;;;;;;GAqBG;AACH,SAAS,oCAAoC,CACrC,CAAS,EACT,CAAa,EACb,EAAc;IAElB,wDAAwD;IACxD,uDAAuD;IACvD,wDAAwD;IACxD,yDAAyD;IACzD,MAAM,CAAC,CAAC,EAAE,EAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC,GAAG,cAAc,CAAC,EAAE,CAAC,CAAC;IAE9C,MAAM,CAAC,EAAE,EAAC,EAAE,CAAC,GAAG,CAAC,CAAC;IAElB,qDAAqD;IACrD,wCAAwC;IACxC,yBAAyB;IACzB,yBAAyB;IACzB,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,EAAC,EAAE,CAAC,EAAC,EAAE,CAAC,CAAC;IAC9B,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,EAAC,EAAE,CAAC,EAAC,EAAE,CAAC,CAAC;IAE9B,sEAAsE;IACtE,6BAA6B;IAC7B,4BAA4B;IAC5B,4BAA4B;IAC5B,8BAA8B;IAC9B,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,EAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC,CAAC;IACvC,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,EAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC,CAAC;IACtC,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,GAAG,CAAC,EAAC,GAAG,CAAC,EAAE,EAAC,GAAG,CAAC,CAAC,CAAC;IAExC,wDAAwD;IACxD,oCAAoC;IACpC,uBAAuB;IACvB,gCAAgC;IAChC,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,EAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAC3C,oDAAoD;IACpD,sDAAsD;IACtD,oBAAoB;IAEpB,qBAAqB;IACrB,oBAAoB;IACpB,0BAA0B;IAC1B,oCAAoC;IACpC,sCAAsC;IACtC,wDAAwD;IAExD,wDAAwD;IACxD,sCAAsC;IACtC,wCAAwC;IACxC,2CAA2C;IAC3C,+CAA+C;IAC/C,mBAAmB;IACnB,MAAM,GAAG,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC;IAEpB,+DAA+D;IAC/D,oDAAoD;IACpD,oBAAoB;IACpB,iCAAiC;IACjC,MAAM,GAAG,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC;IACpB,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,EAAC,GAAG,CAAC,EAAC,GAAG,CAAC,GAAG,EAAC,GAAG,CAAC,CAAC,CAAC;IAC3C,wDAAwD;IAGxD,oBAAoB;IACpB,sBAAsB;IACtB,sBAAsB;IACtB,wDAAwD;IAExD,sCAAsC;IACtC,sDAAsD;IACtD,iCAAiC;IACjC,eAAe;IACf,+BAA+B;IAC/B,qBAAqB;IACrB,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,EAAC,EAAE,CAAC,CAAC;IAEvB,mEAAmE;IACnE,qCAAqC;IACrC,+CAA+C;IAC/C,YAAY;IACZ,qBAAqB;IACrB,iBAAiB;IACjB,0BAA0B;IAC1B,2EAA2E;IAC3E,8BAA8B;IAC9B,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,EAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAC,GAAG,CAAC,EAAC,EAAE,CAAC,CAAC,CAAC;IAEtC,OAAO;QACH,CAAC,EAAE,CAAC,GAAG,CAAC;QACR,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC;QAClB,eAAe;QACf,uBAAuB;QACvB,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC;KACd,CAAC;AACN,CAAC;AAGD,OAAO,EAAE,oCAAoC,EAAE,CAAA"}

package/node/get-medial-points/double-double/dd-get-medial-point-coeffs-bez2.d.ts

@@ -0,0 +1,29 @@
+/**
+ * Returns the polynomial coefficients for the ray parameter `t` and the
+ * curve parameter `s` that encode the medial condition for `q(t) = p + t⋅v`
+ * and a quadratic bezier curve `ps`.
+ *
+ * The returned coefficients describe the equations whose common solutions
+ * satisfy:
+ * * `q(t)` is equidistant from `p` and the nearest point on `ps`
+ * * that common distance is locally minimal among such candidates
+ *
+ * More specifically, this function returns:
+ * * `A` and `B`: the coefficients of `E2(s,t) = A(s)⋅t + B(s)`
+ * * `C` and `D`: the coefficients of `E1(s,t) = C(s)⋅t + D(s)`
+ * * `H`: the eliminated polynomial `A(s)⋅D(s) - B(s)⋅C(s)` whose roots are
+ *   candidate `s` values for medial points
+ *
+ * @param p base point
+ * @param v ray direction from `p`
+ * @param ps quadratic bezier control points, i.e. an order 2 bezier curve
+ * given as an array of control points, e.g. `[[0,0],[1,1],[2,1]]`
+ */
+declare function ddGetMedialPointCoeffsBez2(p: number[][], v: number[][], ps: number[][]): {
+    A: number[][];
+    B: number[][];
+    C: number[][];
+    D: number[][];
+    H: number[][];
+};
+export { ddGetMedialPointCoeffsBez2 };

package/node/get-medial-points/double-double/dd-get-medial-point-coeffs-bez2.js

@@ -0,0 +1,123 @@
+import { toPowerBasisDd } from 'flo-bezier3';
+import { ddMultDouble2, ddAddDd, ddMultDd, ddNegativeOf, ddMultByNeg2, ddMultBy2, ddDiffDouble } from 'double-double';
+const qmd = ddMultDouble2;
+const qaq = ddAddDd;
+const qmq = ddMultDd;
+const qno = ddNegativeOf;
+const qmn2 = ddMultByNeg2;
+const qm2 = ddMultBy2;
+const qdd = ddDiffDouble;
+/**
+ * Returns the polynomial coefficients for the ray parameter `t` and the
+ * curve parameter `s` that encode the medial condition for `q(t) = p + t⋅v`
+ * and a quadratic bezier curve `ps`.
+ *
+ * The returned coefficients describe the equations whose common solutions
+ * satisfy:
+ * * `q(t)` is equidistant from `p` and the nearest point on `ps`
+ * * that common distance is locally minimal among such candidates
+ *
+ * More specifically, this function returns:
+ * * `A` and `B`: the coefficients of `E2(s,t) = A(s)⋅t + B(s)`
+ * * `C` and `D`: the coefficients of `E1(s,t) = C(s)⋅t + D(s)`
+ * * `H`: the eliminated polynomial `A(s)⋅D(s) - B(s)⋅C(s)` whose roots are
+ *   candidate `s` values for medial points
+ *
+ * @param p base point
+ * @param v ray direction from `p`
+ * @param ps quadratic bezier control points, i.e. an order 2 bezier curve
+ * given as an array of control points, e.g. `[[0,0],[1,1],[2,1]]`
+ */
+function ddGetMedialPointCoeffsBez2(p, v, ps) {
+    // -----------------------------------------------------
+    // See get-medial-points.md for implementation details.
+    // -----------------------------------------------------
+    const [px, py] = p;
+    const [vx, vy] = v;
+    const [[x0, y0]] = ps;
+    // Quadratic bezier in power basis: b(s) = a⋅s² + b⋅s + c
+    const [[ax, bx], [ay, by]] = toPowerBasisDd(ps);
+    // u(s) = p - b(s) = u2⋅s² + u1⋅s + u0
+    // const u0x = px - x0;
+    // const u0y = py - y0;
+    const u0x = qdd(px, x0);
+    const u0y = qdd(py, y0);
+    // Shared dot products.
+    // const c1 = 2*(vx*ax + vy*ay);
+    // const c0 = vx*bx + vy*by;
+    // const vu0 = vx*u0x + vy*u0y;
+    // const b4 = ax*ax + ay*ay;
+    // const ab = ax*bx + ay*by;
+    // const bb = bx*bx + by*by;
+    // const u0a = u0x*ax + u0y*ay;
+    // const d0 = u0x*bx + u0y*by;
+    // const b0 = u0x*u0x + u0y*u0y;
+    const c1 = qm2(qaq(qmq(vx, ax), qmq(vy, ay)));
+    const c0 = qaq(qmq(vx, bx), qmq(vy, by));
+    const vu0 = qaq(qmq(vx, u0x), qmq(vy, u0y));
+    const b4 = qaq(qmq(ax, ax), qmq(ay, ay));
+    const ab = qaq(qmq(bx, ax), qmq(by, ay));
+    const bb = qaq(qmq(bx, bx), qmq(by, by));
+    const u0a = qaq(qmq(u0x, ax), qmq(u0y, ay));
+    const d0 = qaq(qmq(bx, u0x), qmq(by, u0y));
+    const b0 = qaq(qmq(u0x, u0x), qmq(u0y, u0y));
+    // -----------------------------------------------------
+    // E1(s,t): (u(s) + t⋅v) ⋅ b'(s) = 0
+    // => C(s)⋅t + D(s) = 0
+    // -----------------------------------------------------
+    // -----------------------------------------------------
+    // const d3 = -2*b4;
+    // const d2 = -3*ab;
+    // const d1 = 2*u0a - bb;
+    const d3 = qmn2(b4);
+    const d2 = qmd(-3, ab);
+    const d1 = qaq(qm2(u0a), qno(bb));
+    // -----------------------------------------------------
+    // -----------------------------------------------------
+    // E2(s,t): |t⋅v|² - |u(s) + t⋅v|² = 0
+    //         => 2⋅(v⋅u(s))⋅t + |u(s)|² = 0
+    //         => A(s)⋅t + B(s) = 0
+    // const a2 = -c1;
+    // const a1 = -2*c0;
+    // const a0 = 2*vu0;
+    const a2 = qno(c1);
+    const a1 = qmn2(c0);
+    const a0 = qm2(vu0);
+    // -----------------------------------------------------
+    // -----------------------------------------------------
+    // const b3 = 2*ab;
+    // const b2 = -d1;
+    // const b1 = -2*d0;
+    const b3 = qm2(ab);
+    const b2 = qno(d1);
+    const b1 = qmn2(d0);
+    // -----------------------------------------------------
+    // Eliminate t from:
+    //   A(s)⋅t + B(s) = 0
+    //   C(s)⋅t + D(s) = 0
+    // by taking A(s)⋅D(s) - B(s)⋅C(s) = 0 (degree ≤ 5 in s)
+    // Using:
+    //   a2 = -c1, d3 = -2*b4, d1 = -b2, b1 = -2*d0
+    // we can compute H = A⋅D - B⋅C directly.
+    // const H5 = b4*c1;
+    // const H4 = ab*c1 + 3*b4*c0;
+    // const H3 = 4*(ab*c0 - vu0*b4);
+    // const H2 = c1*d0 + c0*b2 - 6*vu0*ab;
+    // const H1 = -2*vu0*b2 - c1*b0;
+    // const H0 = 2*vu0*d0 - c0*b0;
+    const H5 = qmq(b4, c1);
+    const H4 = qaq(qmq(ab, c1), qmd(3, qmq(b4, c0)));
+    const H3 = qm2(qm2(qaq(qmq(ab, c0), qno(qmq(vu0, b4)))));
+    const H2 = qaq(qaq(qmq(c1, d0), qmq(c0, b2)), qno(qmd(6, qmq(vu0, ab))));
+    const H1 = qaq(qmn2(qmq(vu0, b2)), qno(qmq(c1, b0)));
+    const H0 = qaq(qm2(qmq(vu0, d0)), qno(qmq(c0, b0)));
+    return {
+        A: [a2, a1, a0],
+        B: [b4, b3, b2, b1, b0],
+        C: [c1, c0],
+        D: [d3, d2, d1, d0],
+        H: [H5, H4, H3, H2, H1, H0]
+    };
+}
+export { ddGetMedialPointCoeffsBez2 };
+//# sourceMappingURL=dd-get-medial-point-coeffs-bez2.js.map

package/node/get-medial-points/double-double/dd-get-medial-point-coeffs-bez2.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"dd-get-medial-point-coeffs-bez2.js","sourceRoot":"","sources":["../../../src/get-medial-points/double-double/dd-get-medial-point-coeffs-bez2.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,cAAc,EAAE,MAAM,aAAa,CAAC;AAC7C,OAAO,EAAE,aAAa,EAAE,OAAO,EAAE,QAAQ,EAAE,YAAY,EAAE,YAAY,EAAE,SAAS,EAAE,YAAY,EAAE,MAAM,eAAe,CAAC;AAEtH,MAAM,GAAG,GAAG,aAAa,CAAC;AAC1B,MAAM,GAAG,GAAG,OAAO,CAAC;AACpB,MAAM,GAAG,GAAG,QAAQ,CAAC;AACrB,MAAM,GAAG,GAAG,YAAY,CAAC;AACzB,MAAM,IAAI,GAAG,YAAY,CAAC;AAC1B,MAAM,GAAG,GAAG,SAAS,CAAC;AACtB,MAAM,GAAG,GAAG,YAAY,CAAC;AAGzB;;;;;;;;;;;;;;;;;;;;GAoBG;AACH,SAAS,0BAA0B,CAC3B,CAAa,EACb,CAAa,EACb,EAAc;IAElB,wDAAwD;IACxD,uDAAuD;IACvD,wDAAwD;IACxD,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC;IACnB,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC;IACnB,MAAM,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC;IAEtB,yDAAyD;IACzD,MAAM,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,GAAG,cAAc,CAAC,EAAE,CAAC,CAAC;IAEhD,sCAAsC;IACtC,uBAAuB;IACvB,uBAAuB;IACvB,MAAM,GAAG,GAAG,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;IACvB,MAAM,GAAG,GAAG,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;IAEvB,uBAAuB;IACvB,gCAAgC;IAChC,4BAA4B;IAC5B,+BAA+B;IAC/B,4BAA4B;IAC5B,4BAA4B;IAC5B,4BAA4B;IAC5B,+BAA+B;IAC/B,8BAA8B;IAC9B,gCAAgC;IAChC,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,EAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAC3C,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,EAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC,CAAC;IACtC,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,GAAG,CAAC,EAAC,GAAG,CAAC,EAAE,EAAC,GAAG,CAAC,CAAC,CAAC;IACzC,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,EAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC,CAAC;IACtC,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,EAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC,CAAC;IACtC,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,EAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC,CAAC;IACtC,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,EAAC,EAAE,CAAC,EAAC,GAAG,CAAC,GAAG,EAAC,EAAE,CAAC,CAAC,CAAC;IACzC,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,GAAG,CAAC,EAAC,GAAG,CAAC,EAAE,EAAC,GAAG,CAAC,CAAC,CAAC;IACxC,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,EAAC,GAAG,CAAC,EAAC,GAAG,CAAC,GAAG,EAAC,GAAG,CAAC,CAAC,CAAC;IAE1C,wDAAwD;IACxD,oCAAoC;IACpC,uBAAuB;IACvB,wDAAwD;IAExD,wDAAwD;IACxD,oBAAoB;IACpB,oBAAoB;IACpB,yBAAyB;IACzB,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;IACpB,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,EAAC,EAAE,CAAC,CAAC;IACtB,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,EAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;IACjC,wDAAwD;IAExD,wDAAwD;IACxD,sCAAsC;IACtC,wCAAwC;IACxC,+BAA+B;IAC/B,kBAAkB;IAClB,oBAAoB;IACpB,oBAAoB;IACpB,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC;IACnB,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;IACpB,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IACpB,wDAAwD;IAExD,wDAAwD;IACxD,mBAAmB;IACnB,kBAAkB;IAClB,oBAAoB;IACpB,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC;IACnB,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC;IACnB,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;IACpB,wDAAwD;IAGxD,oBAAoB;IACpB,sBAAsB;IACtB,sBAAsB;IACtB,wDAAwD;IAExD,SAAS;IACT,+CAA+C;IAC/C,yCAAyC;IACzC,oBAAoB;IACpB,8BAA8B;IAC9B,iCAAiC;IACjC,uCAAuC;IACvC,gCAAgC;IAChC,+BAA+B;IAC/B,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;IACtB,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,EAAC,GAAG,CAAC,CAAC,EAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAC7C,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,EAAC,GAAG,CAAC,GAAG,CAAC,GAAG,EAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACtD,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,EAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC,EAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAC,GAAG,CAAC,GAAG,EAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACnE,MAAM,EAAE,GAAG,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,EAAC,EAAE,CAAC,CAAC,EAAC,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAClD,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,EAAC,EAAE,CAAC,CAAC,EAAC,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAEjD,OAAO;QACH,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;QACf,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;QACvB,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC;QACX,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;QACnB,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;KAC9B,CAAC;AACN,CAAC;AAGD,OAAO,EAAE,0BAA0B,EAAE,CAAA"}