flo-mat 3.0.1 → 3.0.3

package/node/closest-boundary-point/get-best-distance-squared.js.map

@@ -1 +1 @@
-{"version":3,"file":"get-best-distance-squared.js","sourceRoot":"","sources":["../../src/closest-boundary-point/get-best-distance-squared.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,eAAe,EAAE,MAAM,aAAa,CAAC;AAC9C,OAAO,EAAE,sBAAsB,EAAE,MAAM,cAAc,CAAC;AAItD;;;;GAIG;AACH,SAAS,sBAAsB,CACvB,YAA2B,EAC3B,CAAW;IAEf,IAAI,mBAAmB,GAAG,MAAM,CAAC,iBAAiB,CAAC;IACnD,KAAK,IAAI,CAAC,GAAC,CAAC,EAAE,CAAC,GAAC,YAAY,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QACtC,MAAM,WAAW,GAAG,YAAY,CAAC,CAAC,CAAC,CAAC;QAEpC,MAAM,EAAE,GAAG,WAAW,CAAC,KAAK,CAAC,EAAE,CAAC;QAEhC,MAAM,EAAE,GAAG,eAAe,CAAC,EAAE,EAAE,WAAW,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAClD,MAAM,EAAE,GAAG,eAAe,CAAC,EAAE,EAAE,WAAW,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAElD,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CACd,sBAAsB,CAAC,CAAC,EAAE,EAAE,CAAC,EAC7B,sBAAsB,CAAC,CAAC,EAAE,EAAE,CAAC,CAChC,CAAC;QAEF,IAAI,CAAC,GAAG,mBAAmB,EAAE;YACzB,mBAAmB,GAAG,CAAC,CAAC;SAC3B;KACJ;IAED,mEAAmE;IACnE,OAAO,CAAC,CAAC,GAAG,CAAC,IAAE,CAAC,EAAE,CAAC,GAAC,mBAAmB,CAAC;AAC5C,CAAC;AAGD,OAAO,EAAE,sBAAsB,EAAE,CAAA"}
+{"version":3,"file":"get-best-distance-squared.js","sourceRoot":"","sources":["../../src/closest-boundary-point/get-best-distance-squared.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,eAAe,EAAE,MAAM,aAAa,CAAC;AAC9C,OAAO,EAAE,sBAAsB,EAAE,MAAM,cAAc,CAAC;AAItD;;;;GAIG;AACH,SAAS,sBAAsB,CACvB,YAA0B,EAC1B,CAAW;IAEf,IAAI,mBAAmB,GAAG,MAAM,CAAC,iBAAiB,CAAC;IACnD,KAAK,IAAI,CAAC,GAAC,CAAC,EAAE,CAAC,GAAC,YAAY,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACvC,MAAM,WAAW,GAAG,YAAY,CAAC,CAAC,CAAC,CAAC;QAEpC,MAAM,EAAE,GAAG,WAAW,CAAC,KAAK,CAAC,EAAE,CAAC;QAEhC,MAAM,EAAE,GAAG,eAAe,CAAC,EAAE,EAAE,WAAW,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAClD,MAAM,EAAE,GAAG,eAAe,CAAC,EAAE,EAAE,WAAW,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAElD,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CACd,sBAAsB,CAAC,CAAC,EAAE,EAAE,CAAC,EAC7B,sBAAsB,CAAC,CAAC,EAAE,EAAE,CAAC,CAChC,CAAC;QAEF,IAAI,CAAC,GAAG,mBAAmB,EAAE,CAAC;YAC1B,mBAAmB,GAAG,CAAC,CAAC;QAC5B,CAAC;IACL,CAAC;IAED,mEAAmE;IACnE,OAAO,CAAC,CAAC,GAAG,CAAC,IAAE,CAAC,EAAE,CAAC,GAAC,mBAAmB,CAAC;AAC5C,CAAC;AAGD,OAAO,EAAE,sBAAsB,EAAE,CAAA"}

package/node/closest-boundary-point/get-close-boundary-points-certified.d.ts

@@ -1,15 +1,18 @@
-import { BezierPiece } from '../mat/bezier-piece.js';
+import { CurvePiece } from '../mat/curve-piece.js';
 import { Curve } from '../curve/curve.js';
 import { PointOnShape } from '../point-on-shape/point-on-shape.js';
 /**
  * @internal
  * Returns the closest boundary point to the given point, limited to the given
  * bezier pieces, including the beziers actually checked after culling.
+ *
+ * @param pow
  * @param bezierPieces
  * @param x
  * @param touchedCurve
  * @param t
- * @param extreme
+ * @param for1Prong defaults to `false`;
+ * @param angle defaults to `0`
  */
-declare function getCloseBoundaryPointsCertified(bezierPieces: BezierPiece[], x: number[], touchedCurve?: Curve | undefined, t?: number | undefined, for1Prong?: boolean, angle?: number): PointOnShape[];
+declare function getCloseBoundaryPointsCertified(pow: number, bezierPieces: CurvePiece[], x: number[], touchedCurve?: Curve | undefined, t?: number | undefined, for1Prong?: boolean, angle?: number): PointOnShape[];
 export { getCloseBoundaryPointsCertified };

package/node/closest-boundary-point/get-close-boundary-points-certified.js

@@ -5,18 +5,21 @@ import { createPos } from '../point-on-shape/create-pos.js';
  * @internal
  * Returns the closest boundary point to the given point, limited to the given
  * bezier pieces, including the beziers actually checked after culling.
+ *
+ * @param pow
  * @param bezierPieces
  * @param x
  * @param touchedCurve
  * @param t
- * @param extreme
+ * @param for1Prong defaults to `false`;
+ * @param angle defaults to `0`
  */
-function getCloseBoundaryPointsCertified(bezierPieces, x, touchedCurve = undefined, t = undefined, for1Prong = false, angle = 0) {
+function getCloseBoundaryPointsCertified(pow, bezierPieces, x, touchedCurve = undefined, t = undefined, for1Prong = false, angle = 0) {
     bezierPieces = cullBezierPieces1(bezierPieces, x);
     const pInfoss = [];
     for (let i = 0; i < bezierPieces.length; i++) {
         const bezierPiece = bezierPieces[i];
-        const pInfos = getPotentialClosestPointsOnCurveCertified(bezierPiece.curve, x, bezierPiece.ts, touchedCurve, t, for1Prong, angle);
+        const pInfos = getPotentialClosestPointsOnCurveCertified(pow, bezierPiece.curve, x, bezierPiece.ts, touchedCurve, t, for1Prong, angle);
         pInfoss.push(...pInfos);
     }
     /** the minimum max interval value */

package/node/closest-boundary-point/get-close-boundary-points-certified.js.map

@@ -1 +1 @@
-{"version":3,"file":"get-close-boundary-points-certified.js","sourceRoot":"","sources":["../../src/closest-boundary-point/get-close-boundary-points-certified.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,yCAAyC,EAAE,MAAM,sDAAsD,CAAC;AACjH,OAAO,EAAE,iBAAiB,EAAE,MAAM,yBAAyB,CAAC;AAI5D,OAAO,EAAE,SAAS,EAAE,MAAM,iCAAiC,CAAC;AAG5D;;;;;;;;;GASG;AACH,SAAS,+BAA+B,CAChC,YAA2B,EAC3B,CAAW,EACjB,eAAkC,SAAS,EAC3C,IAAwB,SAAS,EACjC,SAAS,GAAG,KAAK,EACjB,KAAK,GAAG,CAAC;IAEV,YAAY,GAAG,iBAAiB,CAAC,YAAY,EAAE,CAAC,CAAC,CAAC;IAElD,MAAM,OAAO,GAA0B,EAAE,CAAC;IAC1C,KAAK,IAAI,CAAC,GAAC,CAAC,EAAE,CAAC,GAAC,YAAY,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QACzC,MAAM,WAAW,GAAG,YAAY,CAAC,CAAC,CAAC,CAAC;QAEpC,MAAM,MAAM,GAAG,yCAAyC,CACvD,WAAW,CAAC,KAAK,EACjB,CAAC,EACD,WAAW,CAAC,EAAE,EACd,YAAY,EACZ,CAAC,EACD,SAAS,EACT,KAAK,CACL,CAAC;QAEF,OAAO,CAAC,IAAI,CAAC,GAAG,MAAM,CAAC,CAAC;KACxB;IAED,qCAAqC;IACrC,IAAI,MAAM,GAAG,MAAM,CAAC,iBAAiB,CAAC;IACtC,KAAK,IAAI,CAAC,GAAC,CAAC,EAAE,CAAC,GAAC,OAAO,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QACpC,MAAM,KAAK,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;QACtC,IAAI,KAAK,GAAG,MAAM,EAAE;YACnB,MAAM,GAAG,KAAK,CAAC;SACf;KACD;IAED,MAAM,iBAAiB,GAA0B,EAAE,CAAC;IAEpD,KAAK,IAAI,CAAC,GAAC,CAAC,EAAE,CAAC,GAAC,OAAO,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QACpC,MAAM,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;QACxB,IAAI,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,IAAI,MAAM,EAAE;YAChC,iBAAiB,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;SAC7B;KACD;IAED,OAAO,iBAAiB,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,SAAS,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC;AAC5E,CAAC;AAGD,OAAO,EAAE,+BAA+B,EAAE,CAAC"}
+{"version":3,"file":"get-close-boundary-points-certified.js","sourceRoot":"","sources":["../../src/closest-boundary-point/get-close-boundary-points-certified.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,yCAAyC,EAAE,MAAM,sDAAsD,CAAC;AACjH,OAAO,EAAE,iBAAiB,EAAE,MAAM,yBAAyB,CAAC;AAI5D,OAAO,EAAE,SAAS,EAAE,MAAM,iCAAiC,CAAC;AAG5D;;;;;;;;;;;;GAYG;AACH,SAAS,+BAA+B,CACtC,GAAW,EACL,YAA0B,EAC1B,CAAW,EACjB,eAAkC,SAAS,EAC3C,IAAwB,SAAS,EACjC,SAAS,GAAG,KAAK,EACjB,KAAK,GAAG,CAAC;IAEV,YAAY,GAAG,iBAAiB,CAAC,YAAY,EAAE,CAAC,CAAC,CAAC;IAElD,MAAM,OAAO,GAA0B,EAAE,CAAC;IAC1C,KAAK,IAAI,CAAC,GAAC,CAAC,EAAE,CAAC,GAAC,YAAY,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAC1C,MAAM,WAAW,GAAG,YAAY,CAAC,CAAC,CAAC,CAAC;QAEpC,MAAM,MAAM,GAAG,yCAAyC,CACvD,GAAG,EACH,WAAW,CAAC,KAAK,EACjB,CAAC,EACD,WAAW,CAAC,EAAE,EACd,YAAY,EACZ,CAAC,EACD,SAAS,EACT,KAAK,CACL,CAAC;QAEF,OAAO,CAAC,IAAI,CAAC,GAAG,MAAM,CAAC,CAAC;IACzB,CAAC;IAED,qCAAqC;IACrC,IAAI,MAAM,GAAG,MAAM,CAAC,iBAAiB,CAAC;IACtC,KAAK,IAAI,CAAC,GAAC,CAAC,EAAE,CAAC,GAAC,OAAO,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACrC,MAAM,KAAK,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;QACtC,IAAI,KAAK,GAAG,MAAM,EAAE,CAAC;YACpB,MAAM,GAAG,KAAK,CAAC;QAChB,CAAC;IACF,CAAC;IAED,MAAM,iBAAiB,GAA0B,EAAE,CAAC;IAEpD,KAAK,IAAI,CAAC,GAAC,CAAC,EAAE,CAAC,GAAC,OAAO,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACrC,MAAM,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;QACxB,IAAI,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,IAAI,MAAM,EAAE,CAAC;YACjC,iBAAiB,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QAC9B,CAAC;IACF,CAAC;IAED,OAAO,iBAAiB,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,SAAS,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC;AAC5E,CAAC;AAGD,OAAO,EAAE,+BAA+B,EAAE,CAAC"}

package/node/closest-boundary-point/get-potential-closest-points-on-curve-certified.d.ts

@@ -2,11 +2,15 @@ import { Curve } from "../curve/curve.js";
 import { FootAndEndpointInfo } from './foot-and-endpoint-info.js';
 /**
  * @internal
+ *
+ * @param pow
  * @param curve The curve
  * @param x The point from which to check
  * @param tRange The allowed t range
  * @param touchedCurve The bezier on which p is located
  * @param t The t value of the bezier that locates p
+ * @param for1Prong defaults to `false;
+ * @param angle defaults to `0`;
  */
-declare function getPotentialClosestPointsOnCurveCertified(curve: Curve, x: number[], [tS, tE]?: number[], touchedCurve?: Curve | undefined, t?: number | undefined, for1Prong?: boolean, angle?: number): FootAndEndpointInfo[];
+declare function getPotentialClosestPointsOnCurveCertified(pow: number, curve: Curve, x: number[], tRange?: number[], touchedCurve?: Curve | undefined, t?: number | undefined, for1Prong?: boolean, angle?: number): FootAndEndpointInfo[];
 export { getPotentialClosestPointsOnCurveCertified };

package/node/closest-boundary-point/get-potential-closest-points-on-curve-certified.js

@@ -9,13 +9,18 @@ const γγ3 = γγ(3);
 const { sqrt } = Math;
 /**
  * @internal
+ *
+ * @param pow
  * @param curve The curve
  * @param x The point from which to check
  * @param tRange The allowed t range
  * @param touchedCurve The bezier on which p is located
  * @param t The t value of the bezier that locates p
+ * @param for1Prong defaults to `false;
+ * @param angle defaults to `0`;
  */
-function getPotentialClosestPointsOnCurveCertified(curve, x, [tS, tE] = [0, 1], touchedCurve = undefined, t = undefined, for1Prong = false, angle = 0) {
+function getPotentialClosestPointsOnCurveCertified(pow, curve, x, tRange = [0, 1], touchedCurve = undefined, t = undefined, for1Prong = false, angle = 0) {
+    const [tS, tE] = tRange;
     const ps = curve.ps;
     const shouldDeflate = angle === 0 && curve === touchedCurve;
     let { polyDd: polyDdO, polyE: polyEO, getPolyExact: getPolyExactO } = getFootPointsOnBezierPolysCertified(ps, x);
@@ -29,10 +34,22 @@ function getPotentialClosestPointsOnCurveCertified(curve, x, [tS, tE] = [0, 1],
         ? ddDeflateWithRunningError(def2.coeffs, def2.errBound.map(e => e / γγ3), t)
         : undefined;
     const { polyDd, polyE, getPolyExact } = def3 !== undefined
-        ? { polyDd: def3.coeffs, polyE: def3.errBound.map(e => e / γγ3), getPolyExact: () => eDeflate(eDeflate(eDeflate(getPolyExactO(), t), t), t) }
+        ? {
+            polyDd: def3.coeffs,
+            polyE: def3.errBound,
+            getPolyExact: () => eDeflate(eDeflate(eDeflate(getPolyExactO(), t), t), t)
+        }
         : def !== undefined
-            ? { polyDd: def.coeffs, polyE: def.errBound.map(e => e / γγ3), getPolyExact: () => eDeflate(getPolyExactO(), t) }
-            : { polyDd: polyDdO, polyE: polyEO, getPolyExact: getPolyExactO };
+            ? {
+                polyDd: def.coeffs,
+                polyE: def.errBound,
+                getPolyExact: () => eDeflate(getPolyExactO(), t)
+            }
+            : {
+                polyDd: polyDdO,
+                polyE: polyEO,
+                getPolyExact: getPolyExactO
+            };
     const ris = allRootsCertified(polyDd, tS, tE, polyE, getPolyExact);
     const dontPush0 = ((t === 1 && curve === touchedCurve.next) ||
         (t === 0 && curve === touchedCurve));
@@ -74,7 +91,7 @@ function getPotentialClosestPointsOnCurveCertified(curve, x, [tS, tE] = [0, 1],
         const t = _t < 0 ? 0 : _t > 1 ? 1 : _t;
         const box = getIntervalBox(ps, [ts, te]);
         const p_ = getPFromBox(box);
-        const dSquaredI = rootIntervalToDistanceSquaredInterval(box, x);
+        const dSquaredI = rootIntervalToDistanceSquaredInterval(pow, box, x);
         const t_ = t === 0 ? 1 : t;
         const curve_ = t === 0 ? curve.prev : curve;
         return {

package/node/closest-boundary-point/get-potential-closest-points-on-curve-certified.js.map

@@ -1 +1 @@
-{"version":3,"file":"get-potential-closest-points-on-curve-certified.js","sourceRoot":"","sources":["../../src/closest-boundary-point/get-potential-closest-points-on-curve-certified.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,iBAAiB,EAAE,MAAM,UAAU,CAAC;AAC7C,OAAO,EAAE,mCAAmC,EAAE,cAAc,EAAE,MAAM,aAAa,CAAC;AAElF,OAAO,EAAE,yBAAyB,EAAE,MAAM,oCAAoC,CAAC;AAC/E,OAAO,EAAE,QAAQ,EAAE,MAAM,UAAU,CAAC;AACpC,OAAO,EAAE,EAAE,EAAE,MAAM,4BAA4B,CAAC;AAChD,OAAO,EAAE,qCAAqC,EAAE,MAAM,iDAAiD,CAAC;AACxG,OAAO,EAAE,WAAW,EAAE,MAAM,qBAAqB,CAAC;AAIlD,MAAM,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;AAElB,MAAM,EAAE,IAAI,EAAE,GAAG,IAAI,CAAC;AAGtB;;;;;;;GAOG;AACH,SAAS,yCAAyC,CAC1C,KAAY,EACZ,CAAW,EACX,CAAC,EAAE,EAAC,EAAE,IAAc,CAAC,CAAC,EAAC,CAAC,CAAC,EACzB,eAAkC,SAAS,EAC3C,IAAwB,SAAS,EACjC,SAAS,GAAG,KAAK,EACjB,KAAK,GAAG,CAAC;IAEb,MAAM,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC;IAEpB,MAAM,aAAa,GAAG,KAAK,KAAK,CAAC,IAAI,KAAK,KAAK,YAAY,CAAC;IAE5D,IAAI,EAAE,MAAM,EAAE,OAAO,EAAE,KAAK,EAAE,MAAM,EAAE,YAAY,EAAE,aAAa,EAAE,GAC/D,mCAAmC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;IAE/C,MAAM,GAAG,GAAG,aAAa;QACjB,CAAC,CAAC,yBAAyB,CAAC,OAAO,EAAE,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAC,GAAG,CAAC,EAAE,CAAE,CAAC;QAChE,CAAC,CAAC,SAAS,CAAC;IAEpB,MAAM,IAAI,GAAG,SAAS,IAAI,aAAa;QAC/B,CAAC,CAAC,yBAAyB,CAAC,GAAI,CAAC,MAAM,EAAE,GAAI,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAC,GAAG,CAAC,EAAE,CAAE,CAAC;QAC3E,CAAC,CAAC,SAAS,CAAC;IAEpB,MAAM,IAAI,GAAG,SAAS,IAAI,aAAa;QAC/B,CAAC,CAAC,yBAAyB,CAAC,IAAK,CAAC,MAAM,EAAE,IAAK,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAC,GAAG,CAAC,EAAE,CAAE,CAAC;QAC7E,CAAC,CAAC,SAAS,CAAC;IAEpB,MAAM,EAAE,MAAM,EAAE,KAAK,EAAE,YAAY,EAAE,GAC3B,IAAI,KAAK,SAAS;QACpB,CAAC,CAAC,EAAE,MAAM,EAAE,IAAI,CAAC,MAAM,EAAE,KAAK,EAAE,IAAI,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAC,GAAG,CAAC,EAAE,YAAY,EAAG,GAAG,EAAE,CAAC,QAAQ,CAAC,QAAQ,CAAC,QAAQ,CAAC,aAAa,EAAE,EAAE,CAAE,CAAC,EAAE,CAAE,CAAC,EAAE,CAAE,CAAC,EAAE;QAC/I,CAAC,CAAC,GAAG,KAAK,SAAS;YACnB,CAAC,CAAC,EAAE,MAAM,EAAE,GAAG,CAAC,MAAM,EAAE,KAAK,EAAE,GAAG,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAC,GAAG,CAAC,EAAE,YAAY,EAAE,GAAG,EAAE,CAAC,QAAQ,CAAC,aAAa,EAAE,EAAE,CAAE,CAAC,EAAE;YAChH,CAAC,CAAC,EAAE,MAAM,EAAE,OAAO,EAAE,KAAK,EAAE,MAAM,EAAE,YAAY,EAAE,aAAa,EAAE,CAAC;IAE1E,MAAM,GAAG,GAAG,iBAAiB,CAAC,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,KAAK,EAAE,YAAY,CAAC,CAAC;IAEnE,MAAM,SAAS,GAAG,CACd,CAAC,CAAC,KAAK,CAAC,IAAI,KAAK,KAAK,YAAa,CAAC,IAAI,CAAC;QACzC,CAAC,CAAC,KAAK,CAAC,IAAI,KAAK,KAAK,YAAY,CAAC,CACtC,CAAC;IACF,MAAM,SAAS,GAAG,CACd,CAAC,CAAC,KAAK,CAAC,IAAI,KAAK,KAAK,YAAa,CAAC,IAAI,CAAC;QACzC,CAAC,CAAC,KAAK,CAAC,IAAI,KAAK,KAAK,YAAY,CAAC,CACtC,CAAC;IAEF,IAAI,EAAE,KAAK,CAAC,EAAE;QACV,IAAI,CAAC,SAAS,EAAE;YAAE,GAAG,CAAC,IAAI,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,CAAC,CAAC;SAAE;KACnE;SAAM,IAAI,EAAE,KAAK,CAAC,EAAE;QACjB,IAAI,CAAC,SAAS,EAAE;YAAE,GAAG,CAAC,IAAI,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,CAAC,CAAC;SAAE;KACnE;SAAM;QACH,GAAG,CAAC,IAAI,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,YAAY,EAAE,CAAC,EAAE,CAAC,CAAC;KACjD;IAED,IAAI,EAAE,KAAK,CAAC,EAAE;QACV,IAAI,CAAC,SAAS,EAAE;YAAE,GAAG,CAAC,IAAI,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,CAAC,CAAC;SAAE;KACnE;SAAM,IAAI,EAAE,KAAK,CAAC,EAAE;QACjB,IAAI,CAAC,SAAS,EAAE;YAAE,GAAG,CAAC,IAAI,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,CAAC,CAAC;SAAE;KACnE;SAAM;QACH,GAAG,CAAC,IAAI,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,YAAY,EAAE,CAAC,EAAE,CAAC,CAAC;KACjD;IAED,MAAM,KAAK,GAAG,GAAG;SAChB,GAAG,CAAC,EAAE,CAAC,EAAE;QACN,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,EAAE,CAAC;QAE9B,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,EAAE;YACpB,OAAO,SAAS,CAAC;SACpB;QAED,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAC,CAAC,CAAC;QACvB,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;QAEvC,MAAM,GAAG,GAAG,cAAc,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QACzC,MAAM,EAAE,GAAG,WAAW,CAAC,GAAG,CAAC,CAAC;QAE5B,MAAM,SAAS,GAAG,qCAAqC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QAEhE,MAAM,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC3B,MAAM,MAAM,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC;QAE5C,OAAO;YACH,KAAK,EAAE,MAAM;YACb,CAAC,EAAE,EAAE;YACL,CAAC,EAAE,EAAE;YACL,CAAC,EAAE,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,GAAC,CAAC;YAC9C,SAAS;YACT,GAAG;SACN,CAAA;IACL,CAAC,CAAC,CAAC;IAEH,MAAM,MAAM,GAAG,KAAK,CAAC,MAAM,CAAC,IAAI,CAAC,EAAE,CAAC,IAAI,KAAK,SAAS,CAA0B,CAAC;IAEjF,OAAO,MAAM,CAAC;AAClB,CAAC;AAGD,OAAO,EAAE,yCAAyC,EAAE,CAAA"}
+{"version":3,"file":"get-potential-closest-points-on-curve-certified.js","sourceRoot":"","sources":["../../src/closest-boundary-point/get-potential-closest-points-on-curve-certified.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,iBAAiB,EAAE,MAAM,UAAU,CAAC;AAC7C,OAAO,EAAE,mCAAmC,EAAE,cAAc,EAAE,MAAM,aAAa,CAAC;AAElF,OAAO,EAAE,yBAAyB,EAAE,MAAM,oCAAoC,CAAC;AAC/E,OAAO,EAAE,QAAQ,EAAE,MAAM,UAAU,CAAC;AACpC,OAAO,EAAE,EAAE,EAAE,MAAM,4BAA4B,CAAC;AAChD,OAAO,EAAE,qCAAqC,EAAE,MAAM,iDAAiD,CAAC;AACxG,OAAO,EAAE,WAAW,EAAE,MAAM,qBAAqB,CAAC;AAIlD,MAAM,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;AAElB,MAAM,EAAE,IAAI,EAAE,GAAG,IAAI,CAAC;AAGtB;;;;;;;;;;;GAWG;AACH,SAAS,yCAAyC,CAC1C,GAAW,EACX,KAAY,EACZ,CAAW,EACX,SAAmB,CAAC,CAAC,EAAC,CAAC,CAAC,EACxB,eAAkC,SAAS,EAC3C,IAAwB,SAAS,EACjC,SAAS,GAAG,KAAK,EACjB,KAAK,GAAG,CAAC;IAEb,MAAM,CAAC,EAAE,EAAC,EAAE,CAAC,GAAG,MAAM,CAAC;IACvB,MAAM,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC;IAEpB,MAAM,aAAa,GAAG,KAAK,KAAK,CAAC,IAAI,KAAK,KAAK,YAAY,CAAC;IAE5D,IAAI,EAAE,MAAM,EAAE,OAAO,EAAE,KAAK,EAAE,MAAM,EAAE,YAAY,EAAE,aAAa,EAAE,GAC/D,mCAAmC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;IAE/C,MAAM,GAAG,GAAG,aAAa;QACjB,CAAC,CAAC,yBAAyB,CAAC,OAAO,EAAE,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAC,GAAG,CAAC,EAAE,CAAE,CAAC;QAChE,CAAC,CAAC,SAAS,CAAC;IAEpB,MAAM,IAAI,GAAG,SAAS,IAAI,aAAa;QAC/B,CAAC,CAAC,yBAAyB,CAAC,GAAI,CAAC,MAAM,EAAE,GAAI,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAC,GAAG,CAAC,EAAE,CAAE,CAAC;QAC3E,CAAC,CAAC,SAAS,CAAC;IAEpB,MAAM,IAAI,GAAG,SAAS,IAAI,aAAa;QAC/B,CAAC,CAAC,yBAAyB,CAAC,IAAK,CAAC,MAAM,EAAE,IAAK,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAC,GAAG,CAAC,EAAE,CAAE,CAAC;QAC7E,CAAC,CAAC,SAAS,CAAC;IAEpB,MAAM,EAAE,MAAM,EAAE,KAAK,EAAE,YAAY,EAAE,GAC3B,IAAI,KAAK,SAAS;QACpB,CAAC,CAAC;YACE,MAAM,EAAE,IAAI,CAAC,MAAM;YACnB,KAAK,EAAE,IAAI,CAAC,QAAQ;YACpB,YAAY,EAAG,GAAG,EAAE,CAAC,QAAQ,CAAC,QAAQ,CAAC,QAAQ,CAAC,aAAa,EAAE,EAAE,CAAE,CAAC,EAAE,CAAE,CAAC,EAAE,CAAE,CAAC;SACjF;QACD,CAAC,CAAC,GAAG,KAAK,SAAS;YACnB,CAAC,CAAC;gBACE,MAAM,EAAE,GAAG,CAAC,MAAM;gBAClB,KAAK,EAAE,GAAG,CAAC,QAAQ;gBACnB,YAAY,EAAE,GAAG,EAAE,CAAC,QAAQ,CAAC,aAAa,EAAE,EAAE,CAAE,CAAC;aACpD;YACD,CAAC,CAAC;gBACE,MAAM,EAAE,OAAO;gBACf,KAAK,EAAE,MAAM;gBACb,YAAY,EAAE,aAAa;aAC9B,CAAC;IAEV,MAAM,GAAG,GAAG,iBAAiB,CAAC,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,KAAK,EAAE,YAAY,CAAC,CAAC;IAEnE,MAAM,SAAS,GAAG,CACd,CAAC,CAAC,KAAK,CAAC,IAAI,KAAK,KAAK,YAAa,CAAC,IAAI,CAAC;QACzC,CAAC,CAAC,KAAK,CAAC,IAAI,KAAK,KAAK,YAAY,CAAC,CACtC,CAAC;IACF,MAAM,SAAS,GAAG,CACd,CAAC,CAAC,KAAK,CAAC,IAAI,KAAK,KAAK,YAAa,CAAC,IAAI,CAAC;QACzC,CAAC,CAAC,KAAK,CAAC,IAAI,KAAK,KAAK,YAAY,CAAC,CACtC,CAAC;IAEF,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;QACX,IAAI,CAAC,SAAS,EAAE,CAAC;YAAC,GAAG,CAAC,IAAI,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,CAAC,CAAC;QAAC,CAAC;IACpE,CAAC;SAAM,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;QAClB,IAAI,CAAC,SAAS,EAAE,CAAC;YAAC,GAAG,CAAC,IAAI,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,CAAC,CAAC;QAAC,CAAC;IACpE,CAAC;SAAM,CAAC;QACJ,GAAG,CAAC,IAAI,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,YAAY,EAAE,CAAC,EAAE,CAAC,CAAC;IAClD,CAAC;IAED,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;QACX,IAAI,CAAC,SAAS,EAAE,CAAC;YAAC,GAAG,CAAC,IAAI,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,CAAC,CAAC;QAAC,CAAC;IACpE,CAAC;SAAM,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;QAClB,IAAI,CAAC,SAAS,EAAE,CAAC;YAAC,GAAG,CAAC,IAAI,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,CAAC,CAAC;QAAC,CAAC;IACpE,CAAC;SAAM,CAAC;QACJ,GAAG,CAAC,IAAI,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,YAAY,EAAE,CAAC,EAAE,CAAC,CAAC;IAClD,CAAC;IAED,MAAM,KAAK,GAAG,GAAG;SAChB,GAAG,CAAC,EAAE,CAAC,EAAE;QACN,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,EAAE,CAAC;QAE9B,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,EAAE,CAAC;YACrB,OAAO,SAAS,CAAC;QACrB,CAAC;QAED,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAC,CAAC,CAAC;QACvB,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;QAEvC,MAAM,GAAG,GAAG,cAAc,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QACzC,MAAM,EAAE,GAAG,WAAW,CAAC,GAAG,CAAC,CAAC;QAE5B,MAAM,SAAS,GAAG,qCAAqC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC;QAErE,MAAM,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC3B,MAAM,MAAM,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC;QAE5C,OAAO;YACH,KAAK,EAAE,MAAM;YACb,CAAC,EAAE,EAAE;YACL,CAAC,EAAE,EAAE;YACL,CAAC,EAAE,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,GAAC,CAAC;YAC9C,SAAS;YACT,GAAG;SACN,CAAA;IACL,CAAC,CAAC,CAAC;IAEH,MAAM,MAAM,GAAG,KAAK,CAAC,MAAM,CAAC,IAAI,CAAC,EAAE,CAAC,IAAI,KAAK,SAAS,CAA0B,CAAC;IAEjF,OAAO,MAAM,CAAC;AAClB,CAAC;AAGD,OAAO,EAAE,yCAAyC,EAAE,CAAA"}

package/node/closest-boundary-point/root-interval-to-distance-squared-interval.d.ts

@@ -9,5 +9,5 @@
  *
  * @internal
  */
-declare function rootIntervalToDistanceSquaredInterval(box: number[][], p: number[]): number[];
+declare function rootIntervalToDistanceSquaredInterval(pow: number, box: number[][], p: number[]): number[];
 export { rootIntervalToDistanceSquaredInterval };

package/node/closest-boundary-point/root-interval-to-distance-squared-interval.js

@@ -1,9 +1,12 @@
-import { twoDiff, eEstimate, eMult, eAdd } from 'big-float-ts';
+import { twoDiff, eEstimate } from 'big-float-ts';
+import { ddMultDd, ddAddDd } from 'double-double';
 // We *have* to do the below to improve performance with bundlers❗ The assignee is a getter❗ The assigned is a pure function❗
 const estimate = eEstimate;
 const td = twoDiff;
-const emult = eMult;
-const eadd = eAdd;
+// const emult = eMult;
+// const eadd = eAdd;
+const qmq = ddMultDd;
+const qaq = ddAddDd;
 const eps = Number.EPSILON;
 /**
  * Returns the distance interval squared given the root interval (currently
@@ -16,7 +19,7 @@ const eps = Number.EPSILON;
  *
  * @internal
  */
-function rootIntervalToDistanceSquaredInterval(box, p) {
+function rootIntervalToDistanceSquaredInterval(pow, box, p) {
     const bl = box[0];
     const tr = box[1];
     const minX = bl[0];
@@ -41,10 +44,12 @@ function rootIntervalToDistanceSquaredInterval(box, p) {
         /** distance to 1st corner of interval box - `distance² = x² + y²` */
         const ax = td(a, x); // a - x
         const by = td(b, y); // b - y
-        const dc1Exact = eadd(emult(ax, ax), emult(by, by)); // ax**2 + bx**2
-        const dc1 = estimate(dc1Exact);
-        const dc1Min = dc1 * (1 - eps); // distance minus max error
-        const dc1Max = dc1 * (1 + eps); // distance plus max error
+        // const dc1Exact = eadd(emult(ax,ax),emult(by,by));  // ax**2 + bx**2
+        // const dc1 = estimate(dc1Exact);
+        const dc1Dd = qaq(qmq(ax, ax), qmq(by, by)); // ax**2 + bx**2
+        const dc1 = dc1Dd[1];
+        const dc1Min = dc1 * (1 - 2 ** pow * eps); // distance minus max error
+        const dc1Max = dc1 * (1 + 2 ** pow * eps); // distance plus max error
         if (dc1Min <= minDSquared) {
             minDSquared = dc1Min;
         }

package/node/closest-boundary-point/root-interval-to-distance-squared-interval.js.map

@@ -1 +1 @@
-{"version":3,"file":"root-interval-to-distance-squared-interval.js","sourceRoot":"","sources":["../../src/closest-boundary-point/root-interval-to-distance-squared-interval.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,EAAE,SAAS,EAAE,KAAK,EAAE,IAAI,EAAE,MAAM,cAAc,CAAC;AAG/D,6HAA6H;AAC7H,MAAM,QAAQ,GAAG,SAAS,CAAC;AAC3B,MAAM,EAAE,GAAG,OAAO,CAAC;AACnB,MAAM,KAAK,GAAG,KAAK,CAAC;AACpB,MAAM,IAAI,GAAG,IAAI,CAAC;AAElB,MAAM,GAAG,GAAG,MAAM,CAAC,OAAO,CAAC;AAG3B;;;;;;;;;;GAUG;AACH,SAAS,qCAAqC,CACtC,GAAe,EACf,CAAW;IAEf,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IAClB,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IAClB,MAAM,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IACnB,MAAM,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IACnB,MAAM,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IACnB,MAAM,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IAEnB,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAE,MAAM;IACvB,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAE,MAAM;IAEvB,IAAI,WAAW,GAAG,MAAM,CAAC,iBAAiB,CAAC;IAC3C,IAAI,WAAW,GAAG,MAAM,CAAC,iBAAiB,CAAC;IAE3C,sCAAsC;IACtC,KAAK,MAAM,CAAC,CAAC,EAAC,CAAC,CAAC,IAAI,CAAC,CAAC,IAAI,EAAC,IAAI,CAAC,EAAC,CAAC,IAAI,EAAC,IAAI,CAAC,EAAC,CAAC,IAAI,EAAC,IAAI,CAAC,EAAC,CAAC,IAAI,EAAC,IAAI,CAAC,CAAC,EAAE;QACnE;;;;;;;;UAQE;QAEF,qEAAqE;QACrE,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,EAAC,CAAC,CAAC,CAAC,CAAE,QAAQ;QAC7B,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,EAAC,CAAC,CAAC,CAAC,CAAE,QAAQ;QAC7B,MAAM,QAAQ,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAC,EAAE,CAAC,EAAC,KAAK,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC,CAAC,CAAE,gBAAgB;QACnE,MAAM,GAAG,GAAG,QAAQ,CAAC,QAAQ,CAAC,CAAC;QAE/B,MAAM,MAAM,GAAG,GAAG,GAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAE,2BAA2B;QAC1D,MAAM,MAAM,GAAG,GAAG,GAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAE,0BAA0B;QAEzD,IAAI,MAAM,IAAI,WAAW,EAAE;YACvB,WAAW,GAAG,MAAM,CAAC;SACxB;QAED,IAAI,MAAM,IAAI,WAAW,EAAE;YACvB,WAAW,GAAG,MAAM,CAAC;SACxB;KACJ;IAED,OAAO,CAAC,WAAW,EAAC,WAAW,CAAC,CAAC;AACrC,CAAC;AAGD,OAAO,EAAE,qCAAqC,EAAE,CAAA"}
+{"version":3,"file":"root-interval-to-distance-squared-interval.js","sourceRoot":"","sources":["../../src/closest-boundary-point/root-interval-to-distance-squared-interval.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,EAAE,SAAS,EAAe,MAAM,cAAc,CAAC;AAC/D,OAAO,EAAE,QAAQ,EAAE,OAAO,EAAE,MAAM,eAAe,CAAC;AAGlD,6HAA6H;AAC7H,MAAM,QAAQ,GAAG,SAAS,CAAC;AAC3B,MAAM,EAAE,GAAG,OAAO,CAAC;AACnB,uBAAuB;AACvB,qBAAqB;AACrB,MAAM,GAAG,GAAG,QAAQ,CAAC;AACrB,MAAM,GAAG,GAAG,OAAO,CAAC;AAEpB,MAAM,GAAG,GAAG,MAAM,CAAC,OAAO,CAAC;AAG3B;;;;;;;;;;GAUG;AACH,SAAS,qCAAqC,CACtC,GAAW,EACX,GAAe,EACf,CAAW;IAEf,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IAClB,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IAClB,MAAM,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IACnB,MAAM,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IACnB,MAAM,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IACnB,MAAM,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IAEnB,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAE,MAAM;IACvB,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAE,MAAM;IAEvB,IAAI,WAAW,GAAG,MAAM,CAAC,iBAAiB,CAAC;IAC3C,IAAI,WAAW,GAAG,MAAM,CAAC,iBAAiB,CAAC;IAE3C,sCAAsC;IACtC,KAAK,MAAM,CAAC,CAAC,EAAC,CAAC,CAAC,IAAI,CAAC,CAAC,IAAI,EAAC,IAAI,CAAC,EAAC,CAAC,IAAI,EAAC,IAAI,CAAC,EAAC,CAAC,IAAI,EAAC,IAAI,CAAC,EAAC,CAAC,IAAI,EAAC,IAAI,CAAC,CAAC,EAAE,CAAC;QACpE;;;;;;;;UAQE;QAEF,qEAAqE;QACrE,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,EAAC,CAAC,CAAC,CAAC,CAAE,QAAQ;QAC7B,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,EAAC,CAAC,CAAC,CAAC,CAAE,QAAQ;QAC7B,sEAAsE;QACtE,kCAAkC;QAClC,MAAM,KAAK,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,EAAC,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC,CAAC,CAAE,gBAAgB;QAC3D,MAAM,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QAErB,MAAM,MAAM,GAAG,GAAG,GAAC,CAAC,CAAC,GAAG,CAAC,IAAE,GAAG,GAAC,GAAG,CAAC,CAAC,CAAE,2BAA2B;QACjE,MAAM,MAAM,GAAG,GAAG,GAAC,CAAC,CAAC,GAAG,CAAC,IAAE,GAAG,GAAC,GAAG,CAAC,CAAC,CAAE,0BAA0B;QAEhE,IAAI,MAAM,IAAI,WAAW,EAAE,CAAC;YACxB,WAAW,GAAG,MAAM,CAAC;QACzB,CAAC;QAED,IAAI,MAAM,IAAI,WAAW,EAAE,CAAC;YACxB,WAAW,GAAG,MAAM,CAAC;QACzB,CAAC;IACL,CAAC;IAED,OAAO,CAAC,WAAW,EAAC,WAAW,CAAC,CAAC;AACrC,CAAC;AAGD,OAAO,EAAE,qCAAqC,EAAE,CAAA"}

package/node/contact-point/contact-point.js.map

@@ -1 +1 @@
-{"version":3,"file":"contact-point.js","sourceRoot":"","sources":["../../src/contact-point/contact-point.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,WAAW,EAAE,MAAM,mCAAmC,CAAC;AAoBhE;;;;;;;GAOG;AACH,SAAS,UAAU,CAAC,CAAe,EAAE,CAAe;IACnD,IAAI,GAAG,GAAG,WAAW,CAAC,CAAC,CAAC,YAAY,EAAE,CAAC,CAAC,YAAY,CAAC,CAAC;IAEtD,IAAI,GAAG,KAAK,CAAC,EAAE;QAAE,OAAO,GAAG,CAAC;KAAE;IAE9B,GAAG,GAAG,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,KAAK,CAAC;IACxB,IAAI,GAAG,KAAK,CAAC,EAAE;QAAE,OAAO,GAAG,CAAC;KAAE;IAE9B,OAAO,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,MAAM,CAAC;AAC5B,CAAC;AAGD,OAAO,EAAgB,UAAU,EAAE,CAAA"}
+{"version":3,"file":"contact-point.js","sourceRoot":"","sources":["../../src/contact-point/contact-point.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,WAAW,EAAE,MAAM,mCAAmC,CAAC;AAoBhE;;;;;;;GAOG;AACH,SAAS,UAAU,CAAC,CAAe,EAAE,CAAe;IACnD,IAAI,GAAG,GAAG,WAAW,CAAC,CAAC,CAAC,YAAY,EAAE,CAAC,CAAC,YAAY,CAAC,CAAC;IAEtD,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;QAAC,OAAO,GAAG,CAAC;IAAC,CAAC;IAE9B,GAAG,GAAG,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,KAAK,CAAC;IACxB,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;QAAC,OAAO,GAAG,CAAC;IAAC,CAAC;IAE9B,OAAO,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,MAAM,CAAC;AAC5B,CAAC;AAGD,OAAO,EAAgB,UAAU,EAAE,CAAA"}

package/node/corner/compare-curvatures-at-interface.js

@@ -1,5 +1,6 @@
-import { eSign, twoProduct, eDiff, expansionProduct, fastExpansionSum, scaleExpansion } from "big-float-ts";
-import { tangentAt0, evaluate2ndDerivativeAt0, toPowerBasis_3rdDerivative } from "flo-bezier3";
+import { eSign, eDiff, fastExpansionSum, scaleExpansion, eMult } from "big-float-ts";
+import { tangentAt0Exact, evaluate2ndDerivativeAt0Exact, toPowerBasis_3rdDerivativeExact } from "flo-bezier3";
+const { sign } = Math;
 /**
  * Compare the curvature, κ, between two curves at t === 0.
  *
@@ -19,16 +20,14 @@ import { tangentAt0, evaluate2ndDerivativeAt0, toPowerBasis_3rdDerivative } from
  */
 function compareCurvaturesAtInterface(psI, psO) {
     // Get x' and y' for incoming curve evaluated at 0
-    const [dxI, dyI] = tangentAt0(psI); // max bitlength increase / max shift === 3
+    const [dxI, dyI] = tangentAt0Exact(psI);
     // Get x'' and y'' for incoming curve evaluated at 0
-    const [ddxI, ddyI] = evaluate2ndDerivativeAt0(psI); // max bitlength increase / max shift === 5
+    const [ddxI, ddyI] = evaluate2ndDerivativeAt0Exact(psI);
     // Get x' and y' for outgoing curve evaluated at 0
-    const [dxO, dyO] = tangentAt0(psO); // max bitlength increase / max shift === 3
+    const [dxO, dyO] = tangentAt0Exact(psO);
     // Get x'' and y'' for outgoing curve evaluated at 0
-    const [ddxO, ddyO] = evaluate2ndDerivativeAt0(psO); // max bitlength increase / max shift === 5
-    //console.log('κI: ', κ(psI, 0));
-    //console.log('κO: ', κ(psO, 0));
-    // Remember the formula for the signed curvature of a parametric curve:
+    const [ddxO, ddyO] = evaluate2ndDerivativeAt0Exact(psO);
+    // Recall the formula for the signed curvature of a parametric curve:
     // κ = x′y′′ - y′x′′ / sqrt(x′² + y′²)³
     // κ² = (x′y′′ - y′x′′)² / (x′² + y′²)³
     // This allows us to do an exact comparison of curvatures
@@ -38,44 +37,28 @@ function compareCurvaturesAtInterface(psI, psO) {
     // <=>  (x₁′y₁′′ - y₁′x₁′′)²(x₂′² + y₂′²)³ > (x₂′y₂′′ - y₂′x₂′′)²(x₁′² + y₁′²)³
     // <=>  a²b³ > c²d³
     // Note b³ > 0 and d³ > 0
-    // max aggregate bitlength increase (let original bitlength === p):
-    // a -> 2 x ((p+3)+(p+5) + 1) === 4p + 18 -> max p in double precision === 8 -> too low
-    //let a = (dxI*ddyI - dyI*ddxI)**2;
-    // b -> 3 x ((p+3) + 1) === 3p + 12
-    //let b = (dxO*dxO  + dyO*dyO )**3;
-    // c -> 2 x ((p+3)+(p+5) + 1) === 4p + 18
-    //let c = (dxO*ddyO - dyO*ddxO)**2;
-    // d -> 3 x ((p+3) + 1) === 3p + 12
-    //let d = (dxI*dxI  + dyI*dyI )**3;
     // We need to resort to exact floating point arithmetic at this point
-    const a = eDiff(twoProduct(dxI, ddyI), twoProduct(dyI, ddxI));
-    const c = eDiff(twoProduct(dxO, ddyO), twoProduct(dyO, ddxO));
+    const a = eDiff(eMult(dxI, ddyI), eMult(dyI, ddxI));
+    const c = eDiff(eMult(dxO, ddyO), eMult(dyO, ddxO));
     const signA = eSign(a);
     const signC = eSign(c);
-    if (signA !== signC) {
-        //console.log('branch 3');
+    if (sign(signA) !== sign(signC)) {
+        // return a negative value if κI < κO
         return signA - signC;
     }
-    const b = fastExpansionSum(twoProduct(dxO, dxO), twoProduct(dyO, dyO));
-    const d = fastExpansionSum(twoProduct(dxI, dxI), twoProduct(dyI, dyI));
-    const b2 = expansionProduct(b, b);
-    const b3 = expansionProduct(b2, b);
-    const d2 = expansionProduct(d, d);
-    const d3 = expansionProduct(d2, d);
+    const b = fastExpansionSum(eMult(dxO, dxO), eMult(dyO, dyO));
+    const d = fastExpansionSum(eMult(dxI, dxI), eMult(dyI, dyI));
+    const b2 = eMult(b, b);
+    const b3 = eMult(b2, b);
+    const d2 = eMult(d, d);
+    const d3 = eMult(d2, d);
     if (signA !== 0 || signC !== 0) {
-        //console.log('branch 4');
-        const a2 = expansionProduct(a, a);
-        const c2 = expansionProduct(c, c);
-        // max aggregate bitlength increase (let original bitlength === p):
-        // κ -> (2 x ((p+3)+(p+5) + 1)) + (3 x ((p+3) + 1)) === 7p + 30
-        // e.g. for bit-aligned input bitlength p of 10 we get output bitlength 
-        // of 100, or for p === 3 (the max exact bitlength allowed to have exact
-        // results without resorting to infinite precision) we get 51 bits.
-        const κI = expansionProduct(a2, b3);
-        const κO = expansionProduct(c2, d3);
+        const a2 = eMult(a, a);
+        const c2 = eMult(c, c);
+        const κI = eMult(a2, b3);
+        const κO = eMult(c2, d3);
         const δκ = eSign(eDiff(κI, κO));
         if (δκ !== 0) {
-            //console.log('branch 5');
             // At this point signA === signC, both +tive or -tive
             return signA > 0 ? δκ : -δκ;
         }
@@ -91,31 +74,31 @@ function compareCurvaturesAtInterface(psI, psO) {
     // <=> (de - 3af)²b⁵ > (bg - 3ch)²d⁵
     // <=> i²b⁵ > j²d⁵
     // Get x′′′ and y′′′ for incoming curve evaluated at 1
-    const [[dddxI], [dddyI]] = toPowerBasis_3rdDerivative(psI); // max bitlength increase === max shift === 6
-    const [[dddxO], [dddyO]] = toPowerBasis_3rdDerivative(psO); // max bitlength increase === max shift === 6
-    const e = eDiff(twoProduct(dxI, dddyI), twoProduct(dyI, dddxI));
-    const f = fastExpansionSum(twoProduct(dxI, ddxI), twoProduct(dyI, ddyI));
-    const g = eDiff(twoProduct(dxO, dddyO), twoProduct(dyO, dddxO));
-    const h = fastExpansionSum(twoProduct(dxO, ddxO), twoProduct(dyO, ddyO));
+    const [[dddxI], [dddyI]] = toPowerBasis_3rdDerivativeExact(psI); // max bitlength increase === max shift === 6
+    const [[dddxO], [dddyO]] = toPowerBasis_3rdDerivativeExact(psO); // max bitlength increase === max shift === 6
+    const e = eDiff(eMult(dxI, dddyI), eMult(dyI, dddxI));
+    const f = fastExpansionSum(eMult(dxI, ddxI), eMult(dyI, ddyI));
+    const g = eDiff(eMult(dxO, dddyO), eMult(dyO, dddxO));
+    const h = fastExpansionSum(eMult(dxO, ddxO), eMult(dyO, ddyO));
     // (de - 3af)²b⁵ > (bg - 3ch)²d⁵
     // i²b⁵ > j²d⁵
-    const i = eDiff(expansionProduct(d, e), scaleExpansion(expansionProduct(a, f), 3));
-    const j = eDiff(expansionProduct(b, g), scaleExpansion(expansionProduct(c, h), 3));
+    const i = eDiff(eMult(d, e), scaleExpansion(eMult(a, f), 3));
+    const j = eDiff(eMult(b, g), scaleExpansion(eMult(c, h), 3));
     const signI = eSign(i);
     const signJ = eSign(j);
-    if (signA !== signC) {
+    if (signI !== signJ) {
         return signI - signJ;
     }
     if (signI === 0 && signJ === 0) {
         // Both curve extensions are identical, i.e. in the same K-family
         return 0;
     }
-    const i2 = expansionProduct(i, i);
-    const b5 = expansionProduct(b2, b3);
-    const j2 = expansionProduct(j, j);
-    const d5 = expansionProduct(d2, d3);
-    const dκI = expansionProduct(i2, b5);
-    const dκO = expansionProduct(j2, d5);
+    const i2 = eMult(i, i);
+    const b5 = eMult(b2, b3);
+    const j2 = eMult(j, j);
+    const d5 = eMult(d2, d3);
+    const dκI = eMult(i2, b5);
+    const dκO = eMult(j2, d5);
     const sgn = eSign(eDiff(dκI, dκO));
     return signI > 0 ? sgn : -sgn;
     // If the above returned value is still zero then the two curve extensions 

package/node/corner/compare-curvatures-at-interface.js.map

@@ -1 +1 @@
-{"version":3,"file":"compare-curvatures-at-interface.js","sourceRoot":"","sources":["../../src/corner/compare-curvatures-at-interface.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,KAAK,EAAE,UAAU,EAAE,KAAK,EAAE,gBAAgB,EAAE,gBAAgB,EAAE,cAAc,EAAE,MAAM,cAAc,CAAC;AAC5G,OAAO,EAAE,UAAU,EAAE,wBAAwB,EAAE,0BAA0B,EAAE,MAAM,aAAa,CAAC;AAG/F;;;;;;;;;;;;;;;;GAgBG;AACH,SAAS,4BAA4B,CAC7B,GAAe,EACf,GAAe;IAEnB,kDAAkD;IAClD,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,UAAU,CAAC,GAAG,CAAC,CAAC,CAAC,2CAA2C;IAC/E,oDAAoD;IACpD,MAAM,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,wBAAwB,CAAC,GAAG,CAAC,CAAC,CAAC,2CAA2C;IAE/F,kDAAkD;IAClD,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,UAAU,CAAC,GAAG,CAAC,CAAC,CAAC,2CAA2C;IAC/E,oDAAoD;IACpD,MAAM,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,wBAAwB,CAAC,GAAG,CAAC,CAAC,CAAC,2CAA2C;IAE/F,iCAAiC;IACjC,iCAAiC;IAEjC,uEAAuE;IACvE,uCAAuC;IACvC,uCAAuC;IAEvC,yDAAyD;IACzD,4EAA4E;IAC5E,mCAAmC;IACnC,6BAA6B;IAC7B,+EAA+E;IAC/E,mBAAmB;IACnB,yBAAyB;IAEzB,mEAAmE;IACnE,uFAAuF;IACvF,mCAAmC;IACnC,mCAAmC;IACnC,mCAAmC;IACnC,yCAAyC;IACzC,mCAAmC;IACnC,mCAAmC;IACnC,mCAAmC;IAEnC,qEAAqE;IACrE,MAAM,CAAC,GAAG,KAAK,CACX,UAAU,CAAC,GAAG,EAAE,IAAI,CAAC,EACrB,UAAU,CAAC,GAAG,EAAE,IAAI,CAAC,CACxB,CAAC;IACF,MAAM,CAAC,GAAG,KAAK,CACX,UAAU,CAAC,GAAG,EAAE,IAAI,CAAC,EACrB,UAAU,CAAC,GAAG,EAAE,IAAI,CAAC,CACxB,CAAC;IAEF,MAAM,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IACvB,MAAM,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IACvB,IAAI,KAAK,KAAK,KAAK,EAAE;QACjB,0BAA0B;QAC1B,OAAO,KAAK,GAAG,KAAK,CAAC;KACxB;IAED,MAAM,CAAC,GAAG,gBAAgB,CACtB,UAAU,CAAC,GAAG,EAAE,GAAG,CAAC,EACpB,UAAU,CAAC,GAAG,EAAE,GAAG,CAAC,CACvB,CAAC;IACF,MAAM,CAAC,GAAG,gBAAgB,CACtB,UAAU,CAAC,GAAG,EAAE,GAAG,CAAC,EACpB,UAAU,CAAC,GAAG,EAAE,GAAG,CAAC,CACvB,CAAC;IAEF,MAAM,EAAE,GAAG,gBAAgB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAClC,MAAM,EAAE,GAAG,gBAAgB,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;IACnC,MAAM,EAAE,GAAG,gBAAgB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAClC,MAAM,EAAE,GAAG,gBAAgB,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;IAEnC,IAAI,KAAK,KAAK,CAAC,IAAI,KAAK,KAAK,CAAC,EAAE;QAC5B,0BAA0B;QAC1B,MAAM,EAAE,GAAG,gBAAgB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAClC,MAAM,EAAE,GAAG,gBAAgB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAElC,mEAAmE;QACnE,+DAA+D;QAC/D,wEAAwE;QACxE,wEAAwE;QACxE,mEAAmE;QAGnE,MAAM,EAAE,GAAG,gBAAgB,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;QACnC,MAAM,EAAE,GAAG,gBAAgB,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;QACnC,MAAM,EAAE,GAAG,KAAK,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QAEhC,IAAI,EAAE,KAAK,CAAC,EAAE;YACV,0BAA0B;YAC1B,qDAAqD;YACrD,OAAO,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;SAC/B;KACJ;IAED,0DAA0D;IAE1D,yEAAyE;IACzE,QAAQ;IACR,kFAAkF;IAClF,2FAA2F;IAC3F,0BAA0B;IAC1B,4FAA4F;IAC5F,0FAA0F;IAC1F,oCAAoC;IACpC,kBAAkB;IAElB,sDAAsD;IACtD,MAAM,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,GAAG,0BAA0B,CAAC,GAAG,CAAC,CAAC,CAAC,6CAA6C;IACzG,MAAM,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,GAAG,0BAA0B,CAAC,GAAG,CAAC,CAAC,CAAC,6CAA6C;IAEzG,MAAM,CAAC,GAAG,KAAK,CACX,UAAU,CAAC,GAAG,EAAE,KAAK,CAAC,EACtB,UAAU,CAAC,GAAG,EAAE,KAAK,CAAC,CACzB,CAAC;IAEF,MAAM,CAAC,GAAG,gBAAgB,CACtB,UAAU,CAAC,GAAG,EAAE,IAAI,CAAC,EACrB,UAAU,CAAC,GAAG,EAAE,IAAI,CAAC,CACxB,CAAC;IAEF,MAAM,CAAC,GAAG,KAAK,CACX,UAAU,CAAC,GAAG,EAAE,KAAK,CAAC,EACtB,UAAU,CAAC,GAAG,EAAE,KAAK,CAAC,CACzB,CAAC;IAEF,MAAM,CAAC,GAAG,gBAAgB,CACtB,UAAU,CAAC,GAAG,EAAE,IAAI,CAAC,EACrB,UAAU,CAAC,GAAG,EAAE,IAAI,CAAC,CACxB,CAAC;IAEF,gCAAgC;IAChC,cAAc;IACd,MAAM,CAAC,GAAG,KAAK,CACX,gBAAgB,CAAC,CAAC,EAAE,CAAC,CAAC,EACtB,cAAc,CACV,gBAAgB,CAAC,CAAC,EAAE,CAAC,CAAC,EACtB,CAAC,CACJ,CACJ,CAAC;IAEF,MAAM,CAAC,GAAG,KAAK,CACX,gBAAgB,CAAC,CAAC,EAAE,CAAC,CAAC,EACtB,cAAc,CACV,gBAAgB,CAAC,CAAC,EAAE,CAAC,CAAC,EACtB,CAAC,CACJ,CACJ,CAAC;IAEF,MAAM,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IACvB,MAAM,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IACvB,IAAI,KAAK,KAAK,KAAK,EAAE;QACjB,OAAO,KAAK,GAAG,KAAK,CAAC;KACxB;IAED,IAAI,KAAK,KAAK,CAAC,IAAI,KAAK,KAAK,CAAC,EAAE;QAC5B,iEAAiE;QACjE,OAAO,CAAC,CAAC;KACZ;IAED,MAAM,EAAE,GAAG,gBAAgB,CAAC,CAAC,EAAC,CAAC,CAAC,CAAC;IACjC,MAAM,EAAE,GAAG,gBAAgB,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;IACnC,MAAM,EAAE,GAAG,gBAAgB,CAAC,CAAC,EAAC,CAAC,CAAC,CAAC;IACjC,MAAM,EAAE,GAAG,gBAAgB,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;IAEnC,MAAM,GAAG,GAAG,gBAAgB,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;IACpC,MAAM,GAAG,GAAG,gBAAgB,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;IAEpC,MAAM,GAAG,GAAG,KAAK,CAAC,KAAK,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC;IAEnC,OAAO,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;IAE9B,2EAA2E;IAC3E,2CAA2C;AAC/C,CAAC;AAGD,OAAO,EAAE,4BAA4B,EAAE,CAAA"}
+{"version":3,"file":"compare-curvatures-at-interface.js","sourceRoot":"","sources":["../../src/corner/compare-curvatures-at-interface.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,KAAK,EAAE,KAAK,EAAE,gBAAgB,EAAE,cAAc,EAAE,KAAK,EAAE,MAAM,cAAc,CAAC;AACrF,OAAO,EAAE,eAAe,EAAE,6BAA6B,EAAE,+BAA+B,EAAE,MAAM,aAAa,CAAC;AAE9G,MAAM,EAAE,IAAI,EAAE,GAAG,IAAI,CAAC;AAGtB;;;;;;;;;;;;;;;;GAgBG;AACH,SAAS,4BAA4B,CAC7B,GAAe,EACf,GAAe;IAEnB,kDAAkD;IAClD,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,eAAe,CAAC,GAAG,CAAC,CAAC;IACxC,oDAAoD;IACpD,MAAM,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,6BAA6B,CAAC,GAAG,CAAC,CAAC;IAExD,kDAAkD;IAClD,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,eAAe,CAAC,GAAG,CAAC,CAAC;IACxC,oDAAoD;IACpD,MAAM,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,6BAA6B,CAAC,GAAG,CAAC,CAAC;IAExD,qEAAqE;IACrE,uCAAuC;IACvC,uCAAuC;IAEvC,yDAAyD;IACzD,4EAA4E;IAC5E,mCAAmC;IACnC,6BAA6B;IAC7B,+EAA+E;IAC/E,mBAAmB;IACnB,yBAAyB;IAEzB,qEAAqE;IACrE,MAAM,CAAC,GAAG,KAAK,CACX,KAAK,CAAC,GAAG,EAAE,IAAI,CAAC,EAChB,KAAK,CAAC,GAAG,EAAE,IAAI,CAAC,CACnB,CAAC;IACF,MAAM,CAAC,GAAG,KAAK,CACX,KAAK,CAAC,GAAG,EAAE,IAAI,CAAC,EAChB,KAAK,CAAC,GAAG,EAAE,IAAI,CAAC,CACnB,CAAC;IAEF,MAAM,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IACvB,MAAM,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IACvB,IAAI,IAAI,CAAC,KAAK,CAAC,KAAK,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QAC9B,qCAAqC;QACrC,OAAO,KAAK,GAAG,KAAK,CAAC;IACzB,CAAC;IAED,MAAM,CAAC,GAAG,gBAAgB,CACtB,KAAK,CAAC,GAAG,EAAE,GAAG,CAAC,EACf,KAAK,CAAC,GAAG,EAAE,GAAG,CAAC,CAClB,CAAC;IACF,MAAM,CAAC,GAAG,gBAAgB,CACtB,KAAK,CAAC,GAAG,EAAE,GAAG,CAAC,EACf,KAAK,CAAC,GAAG,EAAE,GAAG,CAAC,CAClB,CAAC;IAEF,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACvB,MAAM,EAAE,GAAG,KAAK,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;IACxB,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACvB,MAAM,EAAE,GAAG,KAAK,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;IAExB,IAAI,KAAK,KAAK,CAAC,IAAI,KAAK,KAAK,CAAC,EAAE,CAAC;QAC7B,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACvB,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAEvB,MAAM,EAAE,GAAG,KAAK,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;QACxB,MAAM,EAAE,GAAG,KAAK,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;QACxB,MAAM,EAAE,GAAG,KAAK,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QAEhC,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;YACX,qDAAqD;YACrD,OAAO,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;QAChC,CAAC;IACL,CAAC;IAED,0DAA0D;IAE1D,yEAAyE;IACzE,QAAQ;IACR,kFAAkF;IAClF,2FAA2F;IAC3F,0BAA0B;IAC1B,4FAA4F;IAC5F,0FAA0F;IAC1F,oCAAoC;IACpC,kBAAkB;IAElB,sDAAsD;IACtD,MAAM,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,GAAG,+BAA+B,CAAC,GAAG,CAAC,CAAC,CAAC,6CAA6C;IAC9G,MAAM,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,GAAG,+BAA+B,CAAC,GAAG,CAAC,CAAC,CAAC,6CAA6C;IAE9G,MAAM,CAAC,GAAG,KAAK,CACX,KAAK,CAAC,GAAG,EAAE,KAAK,CAAC,EACjB,KAAK,CAAC,GAAG,EAAE,KAAK,CAAC,CACpB,CAAC;IAEF,MAAM,CAAC,GAAG,gBAAgB,CACtB,KAAK,CAAC,GAAG,EAAE,IAAI,CAAC,EAChB,KAAK,CAAC,GAAG,EAAE,IAAI,CAAC,CACnB,CAAC;IAEF,MAAM,CAAC,GAAG,KAAK,CACX,KAAK,CAAC,GAAG,EAAE,KAAK,CAAC,EACjB,KAAK,CAAC,GAAG,EAAE,KAAK,CAAC,CACpB,CAAC;IAEF,MAAM,CAAC,GAAG,gBAAgB,CACtB,KAAK,CAAC,GAAG,EAAE,IAAI,CAAC,EAChB,KAAK,CAAC,GAAG,EAAE,IAAI,CAAC,CACnB,CAAC;IAEF,gCAAgC;IAChC,cAAc;IACd,MAAM,CAAC,GAAG,KAAK,CACX,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EACX,cAAc,CACV,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EACX,CAAC,CACJ,CACJ,CAAC;IAEF,MAAM,CAAC,GAAG,KAAK,CACX,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EACX,cAAc,CACV,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EACX,CAAC,CACJ,CACJ,CAAC;IAEF,MAAM,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IACvB,MAAM,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IACvB,IAAI,KAAK,KAAK,KAAK,EAAE,CAAC;QAClB,OAAO,KAAK,GAAG,KAAK,CAAC;IACzB,CAAC;IAED,IAAI,KAAK,KAAK,CAAC,IAAI,KAAK,KAAK,CAAC,EAAE,CAAC;QAC7B,iEAAiE;QACjE,OAAO,CAAC,CAAC;IACb,CAAC;IAED,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,EAAC,CAAC,CAAC,CAAC;IACtB,MAAM,EAAE,GAAG,KAAK,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;IACxB,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,EAAC,CAAC,CAAC,CAAC;IACtB,MAAM,EAAE,GAAG,KAAK,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;IAExB,MAAM,GAAG,GAAG,KAAK,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;IACzB,MAAM,GAAG,GAAG,KAAK,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;IAEzB,MAAM,GAAG,GAAG,KAAK,CAAC,KAAK,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC;IAEnC,OAAO,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;IAE9B,2EAA2E;IAC3E,2CAA2C;AAC/C,CAAC;AAGD,OAAO,EAAE,4BAA4B,EAAE,CAAA"}

package/node/corner/corner.d.ts

@@ -10,7 +10,6 @@ interface Corner {
      */
     tangents: number[][];
     /** The cross of the unit tangents. */
-    crossTangents: number;
     /**
      * True if the corner is sharp, i.e. > 0 radians as one goes in a positive
      * direction around the shape boundary.

package/node/corner/get-corner.d.ts

@@ -1,4 +1,4 @@
-import { Corner } from "./corner";
+import { Corner } from "./corner.js";
 /**
  * @internal
  * Returns a new corner with properties.

package/node/corner/get-corner.js

@@ -1,4 +1,4 @@
-import { cross, dot, toUnitVector } from "flo-vector2d";
+import { cross, toUnitVector } from "flo-vector2d";
 import { getInterfaceCcw } from './get-interface-ccw.js';
 /**
  * @internal
@@ -30,32 +30,20 @@ const DEGREE_LIMIT = DEGREES[4];
 function getCorner(psI, psO) {
     // getInterfaceCcw must return a number !== 0 if psI and psO are not the
     // same as seen as a curve extension with t ∈ [-∞,+∞]
-    const ccw = getInterfaceCcw(psI, psO);
+    const { ccw, tangentI, tangentO, dotTangents } = getInterfaceCcw(psI, psO);
     const isSharp = ccw < 0;
     const isDull = ccw > 0;
-    // Find (non-normalized) tangent of curve.ps at t === 1
-    const p0E = psI[psI.length - 2];
-    const p1E = psI[psI.length - 1];
-    const xE = p1E[0] - p0E[0];
-    const yE = p1E[1] - p0E[1];
-    const tangentAtEnd = [xE, yE];
-    // Find (non-normalized) tangent of curve.next.ps at t === 0
-    const p0S = psO[0];
-    const p1S = psO[1];
-    const xS = p1S[0] - p0S[0];
-    const yS = p1S[1] - p0S[1];
-    const tangentAtStart = [xS, yS];
     // These use square root and are thus not exact
-    const tangents_ = [
-        toUnitVector(tangentAtEnd),
-        toUnitVector(tangentAtStart),
+    const unitTangents = [
+        toUnitVector(tangentI.map(v => v[1])),
+        toUnitVector(tangentO.map(v => v[1])),
     ];
     // The cross calculated below should be exact due to beziers having been
     // normalized!
-    const crossTangents = cross(tangents_[0], tangents_[1]);
+    const crossTangents = cross(unitTangents[0], unitTangents[1]);
     let isQuiteSharp;
     let isQuiteDull;
-    const dotTangents = dot(tangentAtEnd, tangentAtStart);
+    // const dotTangents = dot(tangentAtIncoming, tangentAtOutgoing);
     if (dotTangents > 0) {
         // Curves go in same direction
         isQuiteSharp = crossTangents < -DEGREE_LIMIT;
@@ -66,8 +54,8 @@ function getCorner(psI, psO) {
         isQuiteDull = isDull;
     }
     return {
-        tangents: tangents_,
-        crossTangents,
+        tangents: unitTangents,
+        // crossTangents, 
         isSharp,
         isDull,
         isQuiteSharp,

package/node/corner/get-corner.js.map

@@ -1 +1 @@
-{"version":3,"file":"get-corner.js","sourceRoot":"","sources":["../../src/corner/get-corner.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,KAAK,EAAE,GAAG,EAAE,YAAY,EAAE,MAAM,cAAc,CAAC;AAExD,OAAO,EAAE,eAAe,EAAE,MAAM,wBAAwB,CAAC;AAGzD;;;GAGG;AACH,MAAM,OAAO,GAAG;IACf,kBAAkB;IAClB,IAAI,EAAG,MAAM;IACb,CAAC,EAAM,MAAM;IACb,CAAC,EAAM,MAAM;IACb,EAAE,EAAK,MAAM;CACb,CAAC;AAGF,gBAAgB;AAChB,kCAAkC;AAClC,MAAM,YAAY,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;AAChC,oCAAoC;AAGpC;;;;;;;;;;;GAWG;AACH,SAAS,SAAS,CAChB,GAAe,EACf,GAAe;IAEhB,wEAAwE;IACxE,qDAAqD;IACrD,MAAM,GAAG,GAAG,eAAe,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;IACtC,MAAM,OAAO,GAAG,GAAG,GAAG,CAAC,CAAC;IACxB,MAAM,MAAM,GAAI,GAAG,GAAG,CAAC,CAAC;IAExB,uDAAuD;IACvD,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;IAC9B,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;IAC9B,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IAC3B,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IAC3B,MAAM,YAAY,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;IAE7B,4DAA4D;IAC5D,MAAM,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IACnB,MAAM,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IACnB,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IAC3B,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IAC3B,MAAM,cAAc,GAAG,CAAC,EAAE,EAAC,EAAE,CAAC,CAAC;IAE/B,+CAA+C;IAC/C,MAAM,SAAS,GAAG;QACjB,YAAY,CAAC,YAAY,CAAC;QAC1B,YAAY,CAAC,cAAc,CAAC;KAC5B,CAAC;IAEF,wEAAwE;IACxE,cAAc;IACd,MAAM,aAAa,GAAG,KAAK,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC;IAExD,IAAI,YAAqB,CAAC;IAC1B,IAAI,WAAoB,CAAC;IACzB,MAAM,WAAW,GAAG,GAAG,CAAC,YAAY,EAAE,cAAc,CAAC,CAAC;IACtD,IAAI,WAAW,GAAG,CAAC,EAAE;QACpB,8BAA8B;QAC9B,YAAY,GAAG,aAAa,GAAG,CAAC,YAAY,CAAC;QAC7C,WAAW,GAAG,aAAa,GAAG,CAAC,YAAY,CAAC;KAC5C;SAAM;QACN,YAAY,GAAG,OAAO,CAAC;QACvB,WAAW,GAAI,MAAM,CAAC;KACtB;IAED,OAAO;QACN,QAAQ,EAAE,SAAS;QACnB,aAAa;QACb,OAAO;QACP,MAAM;QACN,YAAY;QACZ,WAAW;KACX,CAAC;AACH,CAAC;AAGD,OAAO,EAAE,SAAS,EAAE,CAAA"}
+{"version":3,"file":"get-corner.js","sourceRoot":"","sources":["../../src/corner/get-corner.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,KAAK,EAAO,YAAY,EAAE,MAAM,cAAc,CAAC;AAExD,OAAO,EAAE,eAAe,EAAE,MAAM,wBAAwB,CAAC;AAGzD;;;GAGG;AACH,MAAM,OAAO,GAAG;IACf,kBAAkB;IAClB,IAAI,EAAG,MAAM;IACb,CAAC,EAAM,MAAM;IACb,CAAC,EAAM,MAAM;IACb,EAAE,EAAK,MAAM;CACb,CAAC;AAGF,gBAAgB;AAChB,kCAAkC;AAClC,MAAM,YAAY,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;AAChC,oCAAoC;AAGpC;;;;;;;;;;;GAWG;AACH,SAAS,SAAS,CAChB,GAAe,EACf,GAAe;IAEhB,wEAAwE;IACxE,qDAAqD;IACrD,MAAM,EAAE,GAAG,EAAE,QAAQ,EAAE,QAAQ,EAAE,WAAW,EAAE,GAAG,eAAe,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;IAC3E,MAAM,OAAO,GAAG,GAAG,GAAG,CAAC,CAAC;IACxB,MAAM,MAAM,GAAI,GAAG,GAAG,CAAC,CAAC;IAExB,+CAA+C;IAC/C,MAAM,YAAY,GAAG;QACpB,YAAY,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACrC,YAAY,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;KACrC,CAAC;IAEF,wEAAwE;IACxE,cAAc;IACd,MAAM,aAAa,GAAG,KAAK,CAAC,YAAY,CAAC,CAAC,CAAC,EAAE,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC;IAE9D,IAAI,YAAqB,CAAC;IAC1B,IAAI,WAAoB,CAAC;IACzB,iEAAiE;IACjE,IAAI,WAAW,GAAG,CAAC,EAAE,CAAC;QACrB,8BAA8B;QAC9B,YAAY,GAAG,aAAa,GAAG,CAAC,YAAY,CAAC;QAC7C,WAAW,GAAG,aAAa,GAAG,CAAC,YAAY,CAAC;IAC7C,CAAC;SAAM,CAAC;QACP,YAAY,GAAG,OAAO,CAAC;QACvB,WAAW,GAAI,MAAM,CAAC;IACvB,CAAC;IAED,OAAO;QACN,QAAQ,EAAE,YAAY;QACtB,kBAAkB;QAClB,OAAO;QACP,MAAM;QACN,YAAY;QACZ,WAAW;KACX,CAAC;AACH,CAAC;AAGD,OAAO,EAAE,SAAS,EAAE,CAAA"}

package/node/corner/get-interface-ccw.d.ts

@@ -9,16 +9,15 @@
  * The above obviously necessitates that their endpoints coincide as described.
  *
  * Preconditions (for robustness):
- * * The beziers has control points with max bit-length of 25 and bit-aligned.
  * * The bezier does not have infinite curvature at either endpoint
  *
- * This is so the vectors between control points can be
- * calculated exactly without resorting to adaptive infinite precision floating
- * point operations. Note: aligned to 'grid' here means if you bitwise-and all
- * values together the resulting bitlength === the max bithlength of any value.
- *
  * @param psI The incoming bezier that ends at the interface
  * @param psO The outgoing bezier that starts at the interface
  */
-declare function getInterfaceCcw(psI: number[][], psO: number[][]): number;
+declare function getInterfaceCcw(psI: number[][], psO: number[][]): {
+    ccw: number;
+    tangentI: number[][];
+    tangentO: number[][];
+    dotTangents: number;
+};
 export { getInterfaceCcw };