flo-mat 4.0.1 → 4.0.2
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- package/README.md +1 -1
- package/__tests__/find-2-prong/get-medial/get-closest-points.spec.ts +49 -0
- package/browser/index.js +1153 -1205
- package/browser/index.min.js +1 -1
- package/node/closest-boundary-point/foot-and-endpoint-info.d.ts +0 -3
- package/node/closest-boundary-point/get-close-boundary-points-certified.d.ts +1 -1
- package/node/closest-boundary-point/get-close-boundary-points-certified.js +2 -2
- package/node/closest-boundary-point/get-close-boundary-points-certified.js.map +1 -1
- package/node/closest-boundary-point/get-potential-closest-points-on-curve-certified.d.ts +1 -2
- package/node/closest-boundary-point/get-potential-closest-points-on-curve-certified.js +34 -24
- package/node/closest-boundary-point/get-potential-closest-points-on-curve-certified.js.map +1 -1
- package/node/closest-boundary-point/root-interval-to-distance-squared-interval.d.ts +1 -1
- package/node/closest-boundary-point/root-interval-to-distance-squared-interval.js +6 -9
- package/node/closest-boundary-point/root-interval-to-distance-squared-interval.js.map +1 -1
- package/node/cp-node/fs/insert-cp-node.d.ts +2 -2
- package/node/cp-node/fs/insert-cp-node.js.map +1 -1
- package/node/cp-node/fs/is-order-correct.d.ts +2 -2
- package/node/cp-node/fs/is-order-correct.js.map +1 -1
- package/node/cp-node/fs/traverse-cp.js.map +1 -1
- package/node/cp-node/fs/traverse-vertices.d.ts +1 -1
- package/node/cp-node/fs/traverse-vertices.js.map +1 -1
- package/node/debug/functions/draw-elem/draw-branch.d.ts +1 -1
- package/node/debug/functions/draw-elem/draw-branch.js.map +1 -1
- package/node/debug/functions/draw-elem/draw-cp-node.d.ts +1 -1
- package/node/debug/functions/draw-elem/draw-cp-node.js +2 -2
- package/node/debug/functions/draw-elem/draw-cp-node.js.map +1 -1
- package/node/debug/functions/draw-elem/draw-leaves.d.ts +1 -1
- package/node/debug/functions/draw-elem/draw-leaves.js +2 -2
- package/node/debug/functions/draw-elem/draw-leaves.js.map +1 -1
- package/node/debug/functions/draw-elem/draw-mat.d.ts +1 -1
- package/node/debug/functions/draw-elem/draw-mat.js +1 -2
- package/node/debug/functions/draw-elem/draw-mat.js.map +1 -1
- package/node/debug/functions/draw-elem/draw-min-y.d.ts +1 -1
- package/node/debug/functions/draw-elem/draw-min-y.js.map +1 -1
- package/node/debug/functions/draw-elem/draw-three-prong.d.ts +1 -1
- package/node/debug/functions/draw-elem/draw-three-prong.js +2 -2
- package/node/debug/functions/draw-elem/draw-three-prong.js.map +1 -1
- package/node/debug/functions/draw-elem/draw-vertex.d.ts +1 -1
- package/node/debug/functions/draw-elem/draw-vertex.js +0 -2
- package/node/debug/functions/draw-elem/draw-vertex.js.map +1 -1
- package/node/error-analysis/gamma.d.ts +32 -0
- package/node/error-analysis/gamma.js +39 -0
- package/node/error-analysis/gamma.js.map +1 -0
- package/node/find-2-prong/find-2-prong.js +5 -5
- package/node/find-2-prong/find-2-prong.js.map +1 -1
- package/node/find-2-prong/get-medial/get-closest-points.js +27 -23
- package/node/find-2-prong/get-medial/get-closest-points.js.map +1 -1
- package/node/find-3-prong/calc-better-x.d.ts +1 -1
- package/node/find-3-prong/calc-better-x.js +6 -4
- package/node/find-3-prong/calc-better-x.js.map +1 -1
- package/node/find-3-prong/calc-initial-3-prong-center.d.ts +1 -1
- package/node/find-3-prong/calc-initial-3-prong-center.js +4 -3
- package/node/find-3-prong/calc-initial-3-prong-center.js.map +1 -1
- package/node/find-3-prong/find-3-prong-for-delta3s.js +22 -5
- package/node/find-3-prong/find-3-prong-for-delta3s.js.map +1 -1
- package/node/find-3-prong/find-and-add-3-prong.d.ts +2 -2
- package/node/find-3-prong/find-and-add-3-prong.js +15 -0
- package/node/find-3-prong/find-and-add-3-prong.js.map +1 -1
- package/node/find-3-prong/find-and-add-3-prongs-on-loop.js +8 -21
- package/node/find-3-prong/find-and-add-3-prongs-on-loop.js.map +1 -1
- package/node/find-3-prong/find-and-add-3-prongs.d.ts +1 -4
- package/node/find-3-prong/find-and-add-3-prongs.js +1 -2
- package/node/find-3-prong/find-and-add-3-prongs.js.map +1 -1
- package/node/find-3-prong/get-closest-points.d.ts +1 -1
- package/node/find-3-prong/get-closest-points.js +2 -2
- package/node/find-3-prong/get-closest-points.js.map +1 -1
- package/node/find-3-prong/timed-cc.d.ts +3 -0
- package/node/find-3-prong/timed-cc.js +5 -0
- package/node/find-3-prong/timed-cc.js.map +1 -0
- package/node/find-mat/create-initial-cp-tree.d.ts +2 -2
- package/node/find-mat/create-initial-cp-tree.js +2 -2
- package/node/find-mat/create-initial-cp-tree.js.map +1 -1
- package/node/find-mat/find-mat.d.ts +4 -4
- package/node/find-mat/find-mat.js.map +1 -1
- package/node/find-mat/find-mats.js +2 -2
- package/node/find-mat/find-mats.js.map +1 -1
- package/node/find-mat/get-meta.d.ts +2 -2
- package/node/find-mat/get-meta.js +1 -0
- package/node/find-mat/get-meta.js.map +1 -1
- package/node/index.d.ts +1 -0
- package/node/index.js +1 -0
- package/node/index.js.map +1 -1
- package/node/mat/create-new-cp-trees.d.ts +2 -2
- package/node/mat/create-new-cp-trees.js +2 -2
- package/node/mat/create-new-cp-trees.js.map +1 -1
- package/node/mat/get-cp-node-to-left-or-same.d.ts +2 -2
- package/node/mat/get-cp-node-to-left-or-same.js +2 -1
- package/node/mat/get-cp-node-to-left-or-same.js.map +1 -1
- package/node/mat/mat-meta.d.ts +2 -2
- package/node/point-on-shape/get-osculating-circle.js +2 -1
- package/node/point-on-shape/get-osculating-circle.js.map +1 -1
- package/package.json +5 -6
- package/src/closest-boundary-point/foot-and-endpoint-info.ts +0 -3
- package/src/closest-boundary-point/get-close-boundary-points-certified.ts +0 -2
- package/src/closest-boundary-point/get-potential-closest-points-on-curve-certified.ts +41 -27
- package/src/closest-boundary-point/root-interval-to-distance-squared-interval.ts +5 -10
- package/src/cp-node/fs/insert-cp-node.ts +2 -2
- package/src/cp-node/fs/is-order-correct.ts +2 -2
- package/src/cp-node/fs/traverse-cp.ts +1 -1
- package/src/cp-node/fs/traverse-vertices.ts +1 -2
- package/src/debug/functions/check-ordering.ts +0 -1
- package/src/debug/functions/draw-elem/draw-branch.ts +1 -3
- package/src/debug/functions/draw-elem/draw-cp-node.ts +3 -4
- package/src/debug/functions/draw-elem/draw-leaves.ts +3 -3
- package/src/debug/functions/draw-elem/draw-mat.ts +2 -5
- package/src/debug/functions/draw-elem/draw-min-y.ts +2 -2
- package/src/debug/functions/draw-elem/draw-three-prong.ts +3 -3
- package/src/debug/functions/draw-elem/draw-vertex.ts +1 -4
- package/src/error-analysis/gamma.ts +9 -2
- package/src/find-2-prong/find-2-prong.ts +7 -6
- package/src/find-2-prong/get-medial/get-closest-points.ts +24 -26
- package/src/find-3-prong/calc-better-x.ts +7 -4
- package/src/find-3-prong/calc-initial-3-prong-center.ts +2 -3
- package/src/find-3-prong/find-3-prong-for-delta3s.ts +23 -5
- package/src/find-3-prong/find-and-add-3-prong.ts +18 -4
- package/src/find-3-prong/find-and-add-3-prongs-on-loop.ts +7 -23
- package/src/find-3-prong/find-and-add-3-prongs.ts +5 -3
- package/src/find-3-prong/get-closest-points.ts +0 -2
- package/src/find-3-prong/timed-cc.ts +9 -0
- package/src/find-mat/create-initial-cp-tree.ts +3 -3
- package/src/find-mat/find-mat.ts +5 -6
- package/src/find-mat/find-mats.ts +4 -2
- package/src/find-mat/get-meta.ts +3 -3
- package/src/index.ts +2 -1
- package/src/mat/create-new-cp-trees.ts +3 -3
- package/src/mat/get-cp-node-to-left-or-same.ts +4 -3
- package/src/mat/mat-meta.ts +2 -2
- package/src/point-on-shape/get-osculating-circle.ts +2 -1
- package/wallaby.mjs +1 -0
- package/node/closest-boundary-point/get-p-from-box.d.ts +0 -3
- package/node/closest-boundary-point/get-p-from-box.js +0 -11
- package/node/closest-boundary-point/get-p-from-box.js.map +0 -1
|
@@ -1,10 +1,7 @@
|
|
|
1
1
|
import type { Curve } from "flo-boolean";
|
|
2
2
|
interface FootAndEndpointInfo {
|
|
3
3
|
readonly curve: Curve;
|
|
4
|
-
readonly p: number[];
|
|
5
4
|
readonly t: number;
|
|
6
|
-
readonly d: number;
|
|
7
|
-
readonly box: number[][];
|
|
8
5
|
readonly dSquaredI: number[];
|
|
9
6
|
}
|
|
10
7
|
export { FootAndEndpointInfo };
|
|
@@ -10,5 +10,5 @@ import type { PrePointOnShape } from '../point-on-shape/point-on-shape.js';
|
|
|
10
10
|
*
|
|
11
11
|
* @internal
|
|
12
12
|
*/
|
|
13
|
-
declare function getCloseBoundaryPointsCertified(
|
|
13
|
+
declare function getCloseBoundaryPointsCertified(curvePieces: CurvePiece[], x: number[]): PrePointOnShape[];
|
|
14
14
|
export { getCloseBoundaryPointsCertified };
|
|
@@ -11,12 +11,12 @@ import { toP } from '../point-on-shape/to-p.js';
|
|
|
11
11
|
*
|
|
12
12
|
* @internal
|
|
13
13
|
*/
|
|
14
|
-
function getCloseBoundaryPointsCertified(
|
|
14
|
+
function getCloseBoundaryPointsCertified(curvePieces, x) {
|
|
15
15
|
curvePieces = cullCurvePieces1(curvePieces, x);
|
|
16
16
|
const pInfos = [];
|
|
17
17
|
for (let i = 0; i < curvePieces.length; i++) {
|
|
18
18
|
const curvePiece = curvePieces[i];
|
|
19
|
-
const _pInfos = getPotentialClosestPointsOnCurveCertified(
|
|
19
|
+
const _pInfos = getPotentialClosestPointsOnCurveCertified(curvePiece.curve, x, curvePiece.ts);
|
|
20
20
|
pInfos.push(..._pInfos);
|
|
21
21
|
}
|
|
22
22
|
/** the minimum max interval value */
|
|
@@ -1 +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":"AAGA,OAAO,EAAE,yCAAyC,EAAE,MAAM,sDAAsD,CAAC;AACjH,OAAO,EAAE,gBAAgB,EAAE,MAAM,2BAA2B,CAAC;AAC7D,OAAO,EAAE,GAAG,EAAE,MAAM,2BAA2B,CAAC;AAGhD;;;;;;;;;GASG;AACH,SAAS,+BAA+B,CAChC,
|
|
1
|
+
{"version":3,"file":"get-close-boundary-points-certified.js","sourceRoot":"","sources":["../../src/closest-boundary-point/get-close-boundary-points-certified.ts"],"names":[],"mappings":"AAGA,OAAO,EAAE,yCAAyC,EAAE,MAAM,sDAAsD,CAAC;AACjH,OAAO,EAAE,gBAAgB,EAAE,MAAM,2BAA2B,CAAC;AAC7D,OAAO,EAAE,GAAG,EAAE,MAAM,2BAA2B,CAAC;AAGhD;;;;;;;;;GASG;AACH,SAAS,+BAA+B,CAChC,WAAyB,EACzB,CAAW;IAEf,WAAW,GAAG,gBAAgB,CAAC,WAAW,EAAE,CAAC,CAAC,CAAC;IAE/C,MAAM,MAAM,GAA0B,EAAE,CAAC;IACzC,KAAK,IAAI,CAAC,GAAC,CAAC,EAAE,CAAC,GAAC,WAAW,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACtC,MAAM,UAAU,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC;QAElC,MAAM,OAAO,GAAG,yCAAyC,CACrD,UAAU,CAAC,KAAK,EAChB,CAAC,EACD,UAAU,CAAC,EAAqB,CACnC,CAAC;QAEF,MAAM,CAAC,IAAI,CAAC,GAAG,OAAO,CAAC,CAAC;IAC5B,CAAC;IAED,qCAAqC;IACrC,IAAI,MAAM,GAAG,QAAQ,CAAC;IACtB,KAAK,IAAI,CAAC,GAAC,CAAC,EAAE,CAAC,GAAC,MAAM,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACjC,MAAM,KAAK,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;QACrC,IAAI,KAAK,GAAG,MAAM,EAAE,CAAC;YACjB,MAAM,GAAG,KAAK,CAAC;QACnB,CAAC;IACL,CAAC;IAED,MAAM,iBAAiB,GAA0B,EAAE,CAAC;IAEpD,KAAK,IAAI,CAAC,GAAC,CAAC,EAAE,CAAC,GAAC,MAAM,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACjC,MAAM,IAAI,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QACvB,IAAI,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,IAAI,MAAM,EAAE,CAAC;YAC9B,iBAAiB,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QACjC,CAAC;IACL,CAAC;IAED,OAAO,iBAAiB,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE;QAChC,MAAM,EAAE,KAAK,EAAE,CAAC,EAAE,GAAG,IAAI,CAAC;QAC1B,MAAM,EAAE,EAAE,EAAE,GAAG,KAAK,CAAC;QACrB,MAAM,CAAC,GAAG,GAAG,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;QAErB,OAAO;YACH,CAAC,EAAE,CAAC,EAAE,KAAK,EAAE,QAAQ,EAAE,KAAK;SAC/B,CAAC;IACN,CAAC,CAAC,CAAC;AACP,CAAC;AAGD,OAAO,EAAE,+BAA+B,EAAE,CAAC"}
|
|
@@ -3,10 +3,9 @@ import type { FootAndEndpointInfo } from './foot-and-endpoint-info.js';
|
|
|
3
3
|
/**
|
|
4
4
|
* @internal
|
|
5
5
|
*
|
|
6
|
-
* @param maxCoordPowerOf2
|
|
7
6
|
* @param curve the curve
|
|
8
7
|
* @param x the point from which to check
|
|
9
8
|
* @param tRange The allowed t range
|
|
10
9
|
*/
|
|
11
|
-
declare function getPotentialClosestPointsOnCurveCertified(
|
|
10
|
+
declare function getPotentialClosestPointsOnCurveCertified(curve: Curve, x: number[], tRange: [number, number]): FootAndEndpointInfo[];
|
|
12
11
|
export { getPotentialClosestPointsOnCurveCertified };
|
|
@@ -1,42 +1,52 @@
|
|
|
1
|
-
import {
|
|
2
|
-
import { getFootPointsOnBezierPolysCertified, getIntervalBox } from 'flo-bezier3';
|
|
1
|
+
import { roots } from 'flo-poly';
|
|
2
|
+
import { evalDeCasteljauDd, getFootPointsOnBezierPolysCertified, getIntervalBox } from 'flo-bezier3';
|
|
3
3
|
import { rootIntervalToDistanceSquaredInterval } from './root-interval-to-distance-squared-interval.js';
|
|
4
|
-
import {
|
|
5
|
-
|
|
4
|
+
import { squaredDistanceBetweenDd } from "../find-2-prong/squared-distance-between-dd.js";
|
|
5
|
+
import { eps } from '../error-analysis/gamma.js';
|
|
6
6
|
/**
|
|
7
7
|
* @internal
|
|
8
8
|
*
|
|
9
|
-
* @param maxCoordPowerOf2
|
|
10
9
|
* @param curve the curve
|
|
11
10
|
* @param x the point from which to check
|
|
12
11
|
* @param tRange The allowed t range
|
|
13
12
|
*/
|
|
14
|
-
function getPotentialClosestPointsOnCurveCertified(
|
|
13
|
+
function getPotentialClosestPointsOnCurveCertified(curve, x, tRange) {
|
|
15
14
|
const [tS, tE] = tRange;
|
|
16
|
-
const ps = curve
|
|
17
|
-
|
|
18
|
-
|
|
19
|
-
|
|
20
|
-
|
|
21
|
-
|
|
22
|
-
|
|
23
|
-
|
|
24
|
-
|
|
25
|
-
|
|
26
|
-
|
|
27
|
-
|
|
28
|
-
|
|
15
|
+
const { ps } = curve;
|
|
16
|
+
const infos = [];
|
|
17
|
+
if (tS !== tE) {
|
|
18
|
+
const { polyDd: pDd, polyE: pDd_, getPolyExact } = getFootPointsOnBezierPolysCertified(ps, x);
|
|
19
|
+
const ris = roots(pDd, tS, tE, pDd_, getPolyExact) || [];
|
|
20
|
+
infos.push(...ris.map(ri => {
|
|
21
|
+
const { tS: ts, tE: te, t: _t } = ri;
|
|
22
|
+
const t = _t < 0 ? 0 : _t > 1 ? 1 : _t;
|
|
23
|
+
const box = getIntervalBox(ps, [ts, te]);
|
|
24
|
+
const dSquaredI = rootIntervalToDistanceSquaredInterval(box, x);
|
|
25
|
+
const t_ = t === 0 ? 1 : t;
|
|
26
|
+
const curve_ = t === 0 ? curve.prev : curve;
|
|
27
|
+
return {
|
|
28
|
+
curve: curve_,
|
|
29
|
+
t: t_,
|
|
30
|
+
dSquaredI
|
|
31
|
+
};
|
|
32
|
+
}));
|
|
33
|
+
}
|
|
34
|
+
const ts = tRange.slice();
|
|
35
|
+
if (tS === tE) {
|
|
36
|
+
ts.push(tS);
|
|
37
|
+
}
|
|
38
|
+
infos.push(...ts.map(t => {
|
|
39
|
+
const p = evalDeCasteljauDd(ps, [0, t]).map(v => v[0] + v[1]);
|
|
40
|
+
const d = squaredDistanceBetweenDd(p, x);
|
|
41
|
+
const dSquaredI = [d * (1 - eps), d * (1 + eps)];
|
|
29
42
|
const t_ = t === 0 ? 1 : t;
|
|
30
43
|
const curve_ = t === 0 ? curve.prev : curve;
|
|
31
44
|
return {
|
|
32
45
|
curve: curve_,
|
|
33
|
-
p: p_,
|
|
34
46
|
t: t_,
|
|
35
|
-
|
|
36
|
-
dSquaredI,
|
|
37
|
-
box
|
|
47
|
+
dSquaredI
|
|
38
48
|
};
|
|
39
|
-
});
|
|
49
|
+
}));
|
|
40
50
|
return infos;
|
|
41
51
|
}
|
|
42
52
|
export { getPotentialClosestPointsOnCurveCertified };
|
|
@@ -1 +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":"AAEA,OAAO,EAAE,
|
|
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":"AAEA,OAAO,EAAE,KAAK,EAAE,MAAM,UAAU,CAAC;AACjC,OAAO,EAAE,iBAAiB,EAAE,mCAAmC,EAAE,cAAc,EAAE,MAAM,aAAa,CAAC;AACrG,OAAO,EAAE,qCAAqC,EAAE,MAAM,iDAAiD,CAAC;AACxG,OAAO,EAAE,wBAAwB,EAAE,MAAM,gDAAgD,CAAC;AAC1F,OAAO,EAAE,GAAG,EAAE,MAAM,4BAA4B,CAAC;AAGjD;;;;;;GAMG;AACH,SAAS,yCAAyC,CAC1C,KAAY,EACZ,CAAW,EACX,MAAwB;IAE5B,MAAM,CAAC,EAAE,EAAC,EAAE,CAAC,GAAG,MAAM,CAAC;IACvB,MAAM,EAAE,EAAE,EAAE,GAAG,KAAK,CAAC;IAErB,MAAM,KAAK,GAA0B,EAAE,CAAC;IAExC,IAAI,EAAE,KAAK,EAAE,EAAE,CAAC;QACZ,MAAM,EAAE,MAAM,EAAE,GAAG,EAAE,KAAK,EAAE,IAAI,EAAE,YAAY,EAAE,GAC5C,mCAAmC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;QAE/C,MAAM,GAAG,GAAG,KAAK,CAAC,GAAG,EAAE,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,YAAY,CAAC,IAAI,EAAE,CAAC;QAEzD,KAAK,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YACvB,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,GAAG,EAAE,CAAC;YAErC,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;YAEvC,MAAM,GAAG,GAAG,cAAc,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;YAEzC,MAAM,SAAS,GAAG,qCAAqC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;YAEhE,MAAM,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC3B,MAAM,MAAM,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC;YAE5C,OAAO;gBACH,KAAK,EAAE,MAAM;gBACb,CAAC,EAAE,EAAE;gBACL,SAAS;aACZ,CAAA;QACL,CAAC,CAAC,CAAC,CAAC;IACR,CAAC;IAED,MAAM,EAAE,GAAG,MAAM,CAAC,KAAK,EAAE,CAAC;IAC1B,IAAI,EAAE,KAAK,EAAE,EAAE,CAAC;QACZ,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IAChB,CAAC;IAED,KAAK,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE;QACrB,MAAM,CAAC,GAAG,iBAAiB,CAAC,EAAE,EAAE,CAAC,CAAC,EAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAE7D,MAAM,CAAC,GAAG,wBAAwB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAEzC,MAAM,SAAS,GAAG,CAAC,CAAC,GAAC,CAAC,CAAC,GAAG,GAAG,CAAC,EAAE,CAAC,GAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;QAE7C,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,SAAS;SACZ,CAAA;IACL,CAAC,CAAC,CAAC,CAAC;IAEJ,OAAO,KAAK,CAAC;AACjB,CAAC;AAGD,OAAO,EAAE,yCAAyC,EAAE,CAAA"}
|
|
@@ -9,5 +9,5 @@
|
|
|
9
9
|
*
|
|
10
10
|
* @internal
|
|
11
11
|
*/
|
|
12
|
-
declare function rootIntervalToDistanceSquaredInterval(
|
|
12
|
+
declare function rootIntervalToDistanceSquaredInterval(box: number[][], p: number[]): number[];
|
|
13
13
|
export { rootIntervalToDistanceSquaredInterval };
|
|
@@ -1,10 +1,5 @@
|
|
|
1
|
-
import { twoDiff,
|
|
2
|
-
import { ddMultDd, ddAddDd } from 'double-double';
|
|
3
|
-
// We *have* to do the below to improve performance with bundlers❗ The assignee is a getter❗ The assigned is a pure function❗
|
|
4
|
-
const estimate = eEstimate;
|
|
1
|
+
import { twoDiff, ddMultDd, ddAddDd } from 'double-double';
|
|
5
2
|
const td = twoDiff;
|
|
6
|
-
// const emult = eMult;
|
|
7
|
-
// const eadd = eAdd;
|
|
8
3
|
const qmq = ddMultDd;
|
|
9
4
|
const qaq = ddAddDd;
|
|
10
5
|
const eps = Number.EPSILON;
|
|
@@ -19,7 +14,7 @@ const eps = Number.EPSILON;
|
|
|
19
14
|
*
|
|
20
15
|
* @internal
|
|
21
16
|
*/
|
|
22
|
-
function rootIntervalToDistanceSquaredInterval(
|
|
17
|
+
function rootIntervalToDistanceSquaredInterval(box, p) {
|
|
23
18
|
const bl = box[0];
|
|
24
19
|
const tr = box[1];
|
|
25
20
|
const minX = bl[0];
|
|
@@ -48,8 +43,10 @@ function rootIntervalToDistanceSquaredInterval(maxCoordPowerOf2, box, p) {
|
|
|
48
43
|
// const dc1 = estimate(dc1Exact);
|
|
49
44
|
const dc1Dd = qaq(qmq(ax, ax), qmq(by, by)); // ax**2 + bx**2
|
|
50
45
|
const dc1 = dc1Dd[1];
|
|
51
|
-
const dc1Min = dc1
|
|
52
|
-
const dc1Max = dc1
|
|
46
|
+
// const dc1Min = dc1*(1 - 2**maxCoordPowerOf2*eps); // distance minus max error
|
|
47
|
+
// const dc1Max = dc1*(1 + 2**maxCoordPowerOf2*eps); // distance plus max error
|
|
48
|
+
const dc1Min = dc1 * (1 - eps); // distance minus max error
|
|
49
|
+
const dc1Max = dc1 * (1 + eps); // distance plus max error
|
|
53
50
|
if (dc1Min <= minDSquared) {
|
|
54
51
|
minDSquared = dc1Min;
|
|
55
52
|
}
|
|
@@ -1 +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,
|
|
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,QAAQ,EAAE,OAAO,EAAE,MAAM,eAAe,CAAC;AAE3D,MAAM,EAAE,GAAG,OAAO,CAAC;AACnB,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,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,CAAC,QAAQ,CAAC;IAC5B,IAAI,WAAW,GAAG,CAAC,QAAQ,CAAC;IAE5B,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,iFAAiF;QACjF,gFAAgF;QAChF,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,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"}
|
|
@@ -1,7 +1,7 @@
|
|
|
1
1
|
import type { CpNode } from "../cp-node.js";
|
|
2
2
|
import type { PointOnShape } from "../../point-on-shape/point-on-shape.js";
|
|
3
|
-
import {
|
|
4
|
-
declare function insertCpNode(insertIfOrderWrong: boolean, isHoleClosing: boolean, isIntersection: boolean, cpTree:
|
|
3
|
+
import { RbTree } from "flo-ll-rb-tree";
|
|
4
|
+
declare function insertCpNode(insertIfOrderWrong: boolean, isHoleClosing: boolean, isIntersection: boolean, cpTree: RbTree<CpNode>, pos: PointOnShape, _prev: CpNode | undefined, lastInsertId: {
|
|
5
5
|
id: number;
|
|
6
6
|
}): CpNode | undefined;
|
|
7
7
|
export { insertCpNode };
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"insert-cp-node.js","sourceRoot":"","sources":["../../../src/cp-node/fs/insert-cp-node.ts"],"names":[],"mappings":"AAIA,OAAO,EAAE,cAAc,EAAE,MAAM,uBAAuB,CAAC;AAGvD,SAAS,YAAY,CACb,kBAA2B,EAC3B,aAAsB,EACtB,cAAuB,EACvB,
|
|
1
|
+
{"version":3,"file":"insert-cp-node.js","sourceRoot":"","sources":["../../../src/cp-node/fs/insert-cp-node.ts"],"names":[],"mappings":"AAIA,OAAO,EAAE,cAAc,EAAE,MAAM,uBAAuB,CAAC;AAGvD,SAAS,YAAY,CACb,kBAA2B,EAC3B,aAAsB,EACtB,cAAuB,EACvB,MAAsB,EACtB,GAAiB,EACjB,KAAyB,EACzB,YAA4B;IAEhC,IAAI,KAAK,KAAK,SAAS;QACnB,CAAC,cAAc,CAAC,MAAM,EAAE,GAAG,EAAE,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC;QAE3C,IAAI,CAAC,kBAAkB,EAAE,CAAC;YACtB,OAAO,SAAS,CAAC;QACrB,CAAC;IACL,CAAC;IAED,MAAM,MAAM,GAAoB;QAC5B,YAAY,EAAE,GAAG;QACjB,aAAa;QACb,cAAc;QACd,EAAE,EAAE,YAAY,CAAC,EAAE;QACnB,IAAI,EAAE,SAAU;QAChB,IAAI,EAAE,SAAU;QAChB,YAAY,EAAE,SAAU;QACxB,YAAY,EAAE,SAAU;KAC3B,CAAC;IACF,YAAY,CAAC,EAAE,EAAE,CAAC;IAElB,MAAM,IAAI,GAAG,KAAK,KAAK,SAAS,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,KAAK,CAAC;IAClD,MAAM,IAAI,GAAG,KAAK,KAAK,SAAS,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC;IAErD,IAAwB,CAAC,IAAI,GAAG,MAAM,CAAC;IACvC,IAAwB,CAAC,IAAI,GAAG,MAAM,CAAC;IACxC,MAAM,CAAC,IAAI,GAAG,IAAI,CAAC;IACnB,MAAM,CAAC,IAAI,GAAG,IAAI,CAAC;IAEnB,MAAM,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC;IAEtB,OAAO,MAAgB,CAAC;AAC5B,CAAC;AAGD,OAAO,EAAE,YAAY,EAAE,CAAA"}
|
|
@@ -1,5 +1,5 @@
|
|
|
1
1
|
import type { CpNode } from "../cp-node.js";
|
|
2
|
-
import {
|
|
2
|
+
import { RbTree } from "flo-ll-rb-tree";
|
|
3
3
|
import { PointOnShape } from "../../point-on-shape/point-on-shape.js";
|
|
4
|
-
declare function isOrderCorrect(cpTree:
|
|
4
|
+
declare function isOrderCorrect(cpTree: RbTree<CpNode>, pos: PointOnShape, next: CpNode): boolean;
|
|
5
5
|
export { isOrderCorrect };
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"is-order-correct.js","sourceRoot":"","sources":["../../../src/cp-node/fs/is-order-correct.ts"],"names":[],"mappings":"AAGA,OAAO,EAAE,WAAW,EAAE,MAAM,sCAAsC,CAAC;AAGnE,SAAS,cAAc,CACf,
|
|
1
|
+
{"version":3,"file":"is-order-correct.js","sourceRoot":"","sources":["../../../src/cp-node/fs/is-order-correct.ts"],"names":[],"mappings":"AAGA,OAAO,EAAE,WAAW,EAAE,MAAM,sCAAsC,CAAC;AAGnE,SAAS,cAAc,CACf,MAAsB,EACtB,GAAiB,EACjB,IAAY;IAEhB,MAAM,CAAC,GAAG,WAAW,CAAC,GAAG,EAAE,IAAI,CAAC,YAAY,CAAC,CAAC;IAC9C,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QAAC,OAAO,IAAI,CAAC;IAAC,CAAC;IAC3B,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QACV,OAAO,KAAK,CAAC;IACjB,CAAC;IAED,MAAM,OAAO,GAAG,MAAM,CAAC,UAAU,EAAE,CAAC;IACpC,IAAI,OAAO,KAAK,SAAS,EAAE,CAAC,CAAE,oBAAoB;QAC9C,OAAO,IAAI,CAAC;IAChB,CAAC;IAED,MAAM,GAAG,GAAG,OAAO,CAAC,KAAK,CAAC;IAC1B,MAAM,GAAG,GAAG,MAAM,CAAC,UAAU,EAAG,CAAC,KAAK,CAAC;IAEvC,2CAA2C;IAC3C,OAAO,CAAC,WAAW,CAAC,GAAG,CAAC,YAAa,EAAE,GAAG,CAAC,GAAG,CAAC,IAAI,GAAG,KAAK,IAAI,CAAC,CAAC;AACrE,CAAC;AAGD,OAAO,EAAE,cAAc,EAAE,CAAA"}
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"traverse-cp.js","sourceRoot":"","sources":["../../../src/cp-node/fs/traverse-cp.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,aAAa,EAAE,MAAM,qBAAqB,CAAC;AAGpD;;;;;;;;GAQG;AACH,SAAS,UAAU,CACX,OAAe;IAEnB,IAAI,MAAM,GAAG,OAAO,CAAC;IAErB,IAAI,aAAa,CAAC,MAAM,CAAC,EAAE,CAAC;QACxB,OAAO,CAAC,MAAM,CAAC,CAAC,CAAE,mBAAmB;IACzC,CAAC;IAED,MAAM,UAAU,GAAa,EAAE,CAAC;IAChC,GAAG,CAAC;QACA,UAAU,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QAExB,MAAM,IAAI,GAAG,MAAM,CAAC,IAAI,CAAC,YAAY,CAAC;QACtC,MAAM,GAAG,MAAM,KAAK,IAAI;YAChB,CAAC,CAAC,MAAM,GAAG,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,kBAAkB;YAC9C,CAAC,CAAC,MAAM,GAAG,IAAI,CAAA,
|
|
1
|
+
{"version":3,"file":"traverse-cp.js","sourceRoot":"","sources":["../../../src/cp-node/fs/traverse-cp.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,aAAa,EAAE,MAAM,qBAAqB,CAAC;AAGpD;;;;;;;;GAQG;AACH,SAAS,UAAU,CACX,OAAe;IAEnB,IAAI,MAAM,GAAG,OAAO,CAAC;IAErB,IAAI,aAAa,CAAC,MAAM,CAAC,EAAE,CAAC;QACxB,OAAO,CAAC,MAAM,CAAC,CAAC,CAAE,mBAAmB;IACzC,CAAC;IAED,MAAM,UAAU,GAAa,EAAE,CAAC;IAChC,GAAG,CAAC;QACA,UAAU,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QAExB,MAAM,IAAI,GAAG,MAAM,CAAC,IAAI,CAAC,YAAY,CAAC;QACtC,MAAM,GAAG,MAAM,KAAK,IAAI;YAChB,CAAC,CAAC,MAAM,GAAG,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,kBAAkB;YAC9C,CAAC,CAAC,MAAM,GAAG,IAAI,CAAA,CAAa,iBAAiB;IAEzD,CAAC,QAAQ,MAAM,KAAK,OAAO,EAAE;IAE7B,OAAO,UAAU,CAAC;AACtB,CAAC;AAGD,OAAO,EAAE,UAAU,EAAE,CAAA"}
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"traverse-vertices.js","sourceRoot":"","sources":["../../../src/cp-node/fs/traverse-vertices.ts"],"names":[],"mappings":"
|
|
1
|
+
{"version":3,"file":"traverse-vertices.js","sourceRoot":"","sources":["../../../src/cp-node/fs/traverse-vertices.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,cAAc,EAAE,MAAM,wBAAwB,CAAC;AACxD,OAAO,EAAE,aAAa,EAAE,MAAM,qBAAqB,CAAC;AACpD,OAAO,EAAE,wBAAwB,EAAE,MAAM,kCAAkC,CAAC;AAG5E;;;;;;;;GAQG;AACH,SAAS,gBAAgB,CACjB,OAAe,EACf,wBAAkD;IAEtD,wBAAwB,CAAC,OAAO,CAAC,CAAC;IAElC,wEAAwE;IACxE,gBAAgB;IAChB,MAAM,OAAO,GACT,cAAc,CAAC,OAAO,CAAC;SACtB,MAAM,CAAC,MAAM,CAAC,EAAE,CAAC,CAAC,aAAa,CAAC,MAAM,CAAC,CAAC;SACxC,GAAG,CAAC,MAAM,CAAC,EAAE,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC;IAEhC,KAAK,MAAM,MAAM,IAAI,OAAO,EAAE,CAAC;QAC3B,MAAM,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC;QACrB,OAAO,GAAG,CAAC,MAAM,EAAE,CAAC;YAChB,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,EAAG,CAAC;YACtB,wBAAwB,CAAC,EAAE,CAAC,CAAC;YAC7B,GAAG,CAAC,IAAI,CAAC,GAAG,wBAAwB,CAAC,EAAE,CAAC,CAAC,CAAC;QAC9C,CAAC;IACL,CAAC;AACL,CAAC;AAGD,OAAO,EAAE,gBAAgB,EAAE,CAAA"}
|
|
@@ -1,4 +1,4 @@
|
|
|
1
|
-
import { CpNode } from '../../../cp-node/cp-node.js';
|
|
1
|
+
import type { CpNode } from '../../../cp-node/cp-node.js';
|
|
2
2
|
/** @internal */
|
|
3
3
|
declare function drawBranch(g: SVGGElement, branch: CpNode[], classes?: string, delay?: number, scaleFactor?: number): SVGElement[];
|
|
4
4
|
export { drawBranch };
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"draw-branch.js","sourceRoot":"","sources":["../../../../src/debug/functions/draw-elem/draw-branch.ts"],"names":[],"mappings":"
|
|
1
|
+
{"version":3,"file":"draw-branch.js","sourceRoot":"","sources":["../../../../src/debug/functions/draw-elem/draw-branch.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,iBAAiB,EAAE,MAAM,8CAA8C,CAAC;AACjF,OAAO,EAAE,aAAa,EAAE,MAAM,uCAAuC,CAAC;AACtE,OAAO,EAAE,uBAAuB,EAAE,MAAM,wCAAwC,CAAC;AAGjF,gBAAgB;AAChB,SAAS,UAAU,CACX,CAAc,EACd,MAAgB,EAChB,OAAO,GAAG,qBAAqB,EAC/B,KAAK,GAAG,CAAC,EACT,WAAW,GAAG,CAAC;IAGnB,MAAM,KAAK,GAAiB,EAAE,CAAC;IAC/B,IAAI,CAAC,GAAG,CAAC,CAAC;IAEV,MAAM,GAAG,GAAiB,EAAE,CAAC;IAC7B,KAAK,MAAM,MAAM,IAAI,MAAM,EAAE,CAAC;QAC1B,IAAI,aAAa,CAAC,MAAM,CAAC,EAAE,CAAC;YAAC,SAAS;QAAC,CAAC;QACxC,MAAM,MAAM,GAAG,iBAAiB,CAAC,MAAM,CAAC,CAAC;QACzC,IAAI,CAAC,MAAM,EAAE,CAAC;YAAC,SAAS;QAAC,CAAC;QAC1B,CAAC,EAAE,CAAC;QACJ,4DAA4D;QAC5D,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;IACrB,CAAC;IAED,gBAAgB;IAEhB,KAAK,CAAC,IAAI,CAAC,GAAG,uBAAuB,CAAC,CAAC,EAAE,GAAG,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC,CAAC;IAE/D,OAAO,KAAK,CAAC;AACjB,CAAC;AAGD,OAAO,EAAE,UAAU,EAAE,CAAA"}
|
|
@@ -1,4 +1,4 @@
|
|
|
1
|
-
import { CpNode } from '../../../cp-node/cp-node.js';
|
|
1
|
+
import type { CpNode } from '../../../cp-node/cp-node.js';
|
|
2
2
|
/** @internal */
|
|
3
3
|
declare function drawCpNode(g: SVGGElement, cpNode: CpNode, classes?: string, delay?: number, scaleFactor?: number): SVGElement[];
|
|
4
4
|
export { drawCpNode };
|
|
@@ -1,11 +1,11 @@
|
|
|
1
1
|
import { drawFs } from 'flo-draw';
|
|
2
2
|
import { scaleCircle } from '../../../geometry/circle.js';
|
|
3
|
-
import {
|
|
3
|
+
import { getAllOnCircle } from '../../../cp-node/fs/get-all-on-circle.js';
|
|
4
4
|
/** @internal */
|
|
5
5
|
function drawCpNode(g, cpNode, classes = 'blue thin2 nofill', delay = 0, scaleFactor = 1) {
|
|
6
6
|
const circle = scaleCircle(cpNode.pointOnShape.circle, 1);
|
|
7
7
|
const { center: c, radius: r } = circle;
|
|
8
|
-
const poss =
|
|
8
|
+
const poss = getAllOnCircle(cpNode)
|
|
9
9
|
.map(cpNode => cpNode.pointOnShape);
|
|
10
10
|
const { dot, circle: drawCircle, crossHair } = drawFs;
|
|
11
11
|
const $circle = drawCircle(g, circle, classes, delay);
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"draw-cp-node.js","sourceRoot":"","sources":["../../../../src/debug/functions/draw-elem/draw-cp-node.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,MAAM,EAAE,MAAM,UAAU,CAAC;AAClC,OAAO,EAAE,WAAW,EAAE,MAAM,6BAA6B,CAAC;
|
|
1
|
+
{"version":3,"file":"draw-cp-node.js","sourceRoot":"","sources":["../../../../src/debug/functions/draw-elem/draw-cp-node.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,MAAM,EAAE,MAAM,UAAU,CAAC;AAClC,OAAO,EAAE,WAAW,EAAE,MAAM,6BAA6B,CAAC;AAC1D,OAAO,EAAE,cAAc,EAAE,MAAM,0CAA0C,CAAC;AAG1E,gBAAgB;AAChB,SAAS,UAAU,CACX,CAAc,EACd,MAAc,EACd,OAAO,GAAG,mBAAmB,EAC7B,KAAK,GAAG,CAAC,EACT,WAAW,GAAG,CAAC;IAEnB,MAAM,MAAM,GAAG,WAAW,CACtB,MAAM,CAAC,YAAY,CAAC,MAAM,EAC1B,CAAC,CACJ,CAAC;IACF,MAAM,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,GAAG,MAAM,CAAC;IAExC,MAAM,IAAI,GAAG,cAAc,CAAC,MAAM,CAAC;SAC9B,GAAG,CAAC,MAAM,CAAC,EAAE,CAAC,MAAM,CAAC,YAAY,CAAC,CAAA;IAEvC,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,MAAM,CAAC;IAEtD,MAAM,OAAO,GAAG,UAAU,CAAC,CAAC,EAAE,MAAM,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;IAEtD,MAAM,IAAI,GAAiB,EAAE,CAAC;IAC9B,MAAM,GAAG,GAAiB,EAAE,CAAC;IAC7B,KAAK,IAAI,CAAC,GAAC,CAAC,EAAE,CAAC,GAAC,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACpB,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,GAAC,CAAC,CAAC,GAAC,CAAC,CAAC,GAAC,WAAW,EAAE,MAAM,EAAE,KAAK,CAAC,CAAC,CAAC;QAC/D,GAAG,CAAC,IAAI,CAAC,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,WAAW,EAAE,KAAK,CAAC,CAAC,CAAC;IAC5D,CAAC;IAED,MAAM,MAAM,GAAG,SAAS,CAAC,CAAC,EAAE,CAAC,EAAE,sBAAsB,EAAE,CAAC,GAAC,WAAW,EAAE,KAAK,CAAC,CAAC;IAE7E,iBAAiB;IACjB,+FAA+F;IAC/F,IAAI;IAEJ,OAAO,CAAC,GAAG,IAAI,EAAE,GAAG,GAAG,EAAE,GAAG,OAAO,EAAE,GAAG,MAAM,CAAiB,CAAC;AACpE,CAAC;AAGD,OAAO,EAAE,UAAU,EAAE,CAAA"}
|
|
@@ -1,4 +1,4 @@
|
|
|
1
|
-
import { CpNode } from '../../../cp-node/cp-node.js';
|
|
1
|
+
import type { CpNode } from '../../../cp-node/cp-node.js';
|
|
2
2
|
/** @internal */
|
|
3
3
|
declare function drawLeaves(g: SVGGElement, leaves: CpNode[], classes?: string, delay?: number, scaleFactor?: number): SVGElement[];
|
|
4
4
|
export { drawLeaves };
|
|
@@ -1,6 +1,6 @@
|
|
|
1
1
|
import { drawFs } from 'flo-draw';
|
|
2
|
-
import { CpNodeFs } from '../../../cp-node/cp-node-fs.js';
|
|
3
2
|
import { drawCirclePercent } from './draw-circle-percent.js';
|
|
3
|
+
import { getAllOnCircle } from '../../../cp-node/fs/get-all-on-circle.js';
|
|
4
4
|
/** @internal */
|
|
5
5
|
function drawLeaves(g, leaves, classes = 'thin10 deeppink nofill', delay = 0, scaleFactor = 1) {
|
|
6
6
|
const $elems = [];
|
|
@@ -10,7 +10,7 @@ function drawLeaves(g, leaves, classes = 'thin10 deeppink nofill', delay = 0, sc
|
|
|
10
10
|
const { circle } = cp;
|
|
11
11
|
const { center: c, radius: r } = circle;
|
|
12
12
|
const $center = [drawCirclePercent(g, c, 0.5, 'pinker thin5 nofill')];
|
|
13
|
-
const poss =
|
|
13
|
+
const poss = getAllOnCircle(cpNode)
|
|
14
14
|
.map(cpNode => cpNode.pointOnShape);
|
|
15
15
|
const $ls = [];
|
|
16
16
|
const $cps = [];
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"draw-leaves.js","sourceRoot":"","sources":["../../../../src/debug/functions/draw-elem/draw-leaves.ts"],"names":[],"mappings":"
|
|
1
|
+
{"version":3,"file":"draw-leaves.js","sourceRoot":"","sources":["../../../../src/debug/functions/draw-elem/draw-leaves.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,MAAM,EAAE,MAAM,UAAU,CAAC;AAClC,OAAO,EAAE,iBAAiB,EAAE,MAAM,0BAA0B,CAAC;AAC7D,OAAO,EAAE,cAAc,EAAE,MAAM,0CAA0C,CAAC;AAG1E,gBAAgB;AAChB,SAAS,UAAU,CACX,CAAc,EACd,MAAgB,EAChB,OAAO,GAAG,wBAAwB,EAClC,KAAK,GAAG,CAAC,EACT,WAAW,GAAG,CAAC;IAEnB,MAAM,MAAM,GAAiB,EAAE,CAAC;IAChC,MAAM,EAAE,IAAI,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,GAAG,MAAM,CAAC;IAEjD,KAAK,MAAM,MAAM,IAAI,MAAM,EAAE,CAAC;QAC1B,MAAM,EAAE,GAAG,MAAM,CAAC,YAAY,CAAC;QAC/B,MAAM,EAAE,MAAM,EAAE,GAAG,EAAE,CAAC;QACtB,MAAM,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,GAAG,MAAM,CAAC;QAExC,MAAM,OAAO,GAAG,CAAC,iBAAiB,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,qBAAqB,CAAC,CAAC,CAAC;QAEtE,MAAM,IAAI,GAAG,cAAc,CAAC,MAAM,CAAC;aAC9B,GAAG,CAAC,MAAM,CAAC,EAAE,CAAC,MAAM,CAAC,YAAY,CAAC,CAAC;QAExC,MAAM,GAAG,GAAiB,EAAE,CAAC;QAC7B,MAAM,IAAI,GAAiB,EAAE,CAAC;QAC9B,KAAK,IAAI,CAAC,GAAC,CAAC,EAAE,CAAC,GAAC,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACpB,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,GAAC,CAAC,CAAC,GAAC,CAAC,CAAC,GAAC,WAAW,EAAE,QAAQ,EAAE,KAAK,CAAC,CAAC,CAAC;YACjE,GAAG,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,WAAW,EAAE,KAAK,CAAC,CAAC,CAAC;QACrD,CAAC;QAED,MAAM,OAAO,GAAG,UAAU,CAAC,CAAC,EAAE,MAAM,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;QACtD,6DAA6D;QAE7D,MAAM,CAAC,IAAI,CAAC,GAAG;YACX,GAAG,OAAO;YACV,GAAG,IAAI;YACP,GAAG,GAAG;YACN,GAAG,OAAO;SACb,CAAC,CAAC;IACP,CAAC;IAED,OAAO,MAAM,CAAC;AAClB,CAAC;AAGD,OAAO,EAAE,UAAU,EAAE,CAAA"}
|
|
@@ -1,8 +1,7 @@
|
|
|
1
1
|
import { drawFs } from 'flo-draw';
|
|
2
2
|
import { traverseEdges } from '../../../cp-node/fs/traverse-edges.js';
|
|
3
3
|
import { getMatCurveToNext } from '../../../cp-node/fs/get-mat-curve-to-next.js';
|
|
4
|
-
import {
|
|
5
|
-
const { isTerminating } = CpNodeFs;
|
|
4
|
+
import { isTerminating } from '../../../cp-node/fs/is-terminating.js';
|
|
6
5
|
const altClasses = [
|
|
7
6
|
'thin10 blue nofill',
|
|
8
7
|
'thin10 red nofill'
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"draw-mat.js","sourceRoot":"","sources":["../../../../src/debug/functions/draw-elem/draw-mat.ts"],"names":[],"mappings":"
|
|
1
|
+
{"version":3,"file":"draw-mat.js","sourceRoot":"","sources":["../../../../src/debug/functions/draw-elem/draw-mat.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,MAAM,EAAE,MAAM,UAAU,CAAC;AAClC,OAAO,EAAE,aAAa,EAAE,MAAM,uCAAuC,CAAC;AACtE,OAAO,EAAE,iBAAiB,EAAE,MAAM,8CAA8C,CAAC;AACjF,OAAO,EAAE,aAAa,EAAE,MAAM,uCAAuC,CAAC;AAGtE,MAAM,UAAU,GAAG;IACf,oBAAoB;IACpB,mBAAmB;CACtB,CAAC;AAEF,gBAAgB;AAChB,SAAS,OAAO,CACR,CAAc,EACd,GAAQ,EACR,QAAiB,EACjB,KAAK,GAAG,CAAC,EACT,WAAW,GAAG,CAAC;IAEnB,IAAI,MAAM,GAAG,GAAG,CAAC,MAAM,CAAC;IACxB,IAAI,CAAC,MAAM,EAAE,CAAC;QAAC,OAAO,EAAE,CAAC;IAAC,CAAC;IAE3B,OAAO,CAAC,aAAa,CAAC,MAAM,CAAC,EAAE,CAAC;QAAC,MAAM,GAAG,MAAM,CAAC,IAAI,CAAC;IAAC,CAAC;IAExD,MAAM,KAAK,GAAiB,EAAE,CAAC;IAE/B,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,aAAa,CAAC,MAAM,EAAE,MAAM,CAAC,EAAE;QAC3B,IAAI,aAAa,CAAC,MAAM,CAAC,EAAE,CAAC;YAAC,OAAO;QAAC,CAAC;QAEtC,MAAM,MAAM,GAAG,iBAAiB,CAAC,MAAM,CAAC,CAAC;QACzC,IAAI,CAAC,MAAM,EAAE,CAAC;YAAC,OAAO;QAAC,CAAC;QAExB,KAAK,CAAC,IAAI,CAAE,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,EAAE,MAAM,EAAE,UAAU,CAAC,CAAC,GAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC;QACjE,CAAC,EAAE,CAAC;IACR,CAAC,CAAC,CAAC;IAEH,OAAO,KAAK,CAAC;AACjB,CAAC;AAGD,OAAO,EAAE,OAAO,EAAE,CAAA"}
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"draw-min-y.js","sourceRoot":"","sources":["../../../../src/debug/functions/draw-elem/draw-min-y.ts"],"names":[],"mappings":"
|
|
1
|
+
{"version":3,"file":"draw-min-y.js","sourceRoot":"","sources":["../../../../src/debug/functions/draw-elem/draw-min-y.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,eAAe,EAAE,MAAM,aAAa,CAAC;AAC9C,OAAO,EAAE,MAAM,EAAE,MAAM,UAAU,CAAC;AAGlC,gBAAgB;AAChB,SAAS,IAAI,CAAC,CAAc,EAAE,GAAiB;IAC3C,MAAM,CAAC,GAAG,eAAe,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;IAE/C,MAAM,MAAM,GAAG,MAAM,CAAC,SAAS,CACvB,CAAC,EAAE,CAAC,EAAE,mBAAmB,CAChC,CAAC;IAEF,OAAO,MAAM,CAAC;AAClB,CAAC;AAGD,OAAO,EAAE,IAAI,EAAE,CAAA"}
|
|
@@ -1,4 +1,4 @@
|
|
|
1
|
-
import { CpNode } from '../../../cp-node/cp-node.js';
|
|
1
|
+
import type { CpNode } from '../../../cp-node/cp-node.js';
|
|
2
2
|
/** @internal */
|
|
3
3
|
declare function drawThreeProng(g: SVGGElement, threeProng: CpNode, classes?: string, delay?: number, scaleFactor?: number): SVGElement[];
|
|
4
4
|
export { drawThreeProng };
|
|
@@ -1,6 +1,6 @@
|
|
|
1
1
|
import { drawFs } from 'flo-draw';
|
|
2
2
|
import { scaleCircle } from '../../../geometry/circle.js';
|
|
3
|
-
import {
|
|
3
|
+
import { getAllOnCircle } from '../../../cp-node/fs/get-all-on-circle.js';
|
|
4
4
|
/** @internal */
|
|
5
5
|
function drawThreeProng(g,
|
|
6
6
|
// threeProng: ThreeProngForDebugging,
|
|
@@ -8,7 +8,7 @@ threeProng, classes, delay = 0, scaleFactor = 1) {
|
|
|
8
8
|
const circle = scaleCircle(threeProng.pointOnShape.circle, 1);
|
|
9
9
|
const { center: c, radius: r } = circle;
|
|
10
10
|
// const poss = threeProng.poss;
|
|
11
|
-
const poss =
|
|
11
|
+
const poss = getAllOnCircle(threeProng)
|
|
12
12
|
.map(cpNode => cpNode.pointOnShape);
|
|
13
13
|
const { dot, circle: drawCircle, crossHair } = drawFs;
|
|
14
14
|
const $circle = drawCircle(g, circle, 'blue thin2 nofill', delay);
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"draw-three-prong.js","sourceRoot":"","sources":["../../../../src/debug/functions/draw-elem/draw-three-prong.ts"],"names":[],"mappings":"
|
|
1
|
+
{"version":3,"file":"draw-three-prong.js","sourceRoot":"","sources":["../../../../src/debug/functions/draw-elem/draw-three-prong.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,MAAM,EAAE,MAAM,UAAU,CAAC;AAClC,OAAO,EAAE,WAAW,EAAE,MAAM,6BAA6B,CAAC;AAC1D,OAAO,EAAE,cAAc,EAAE,MAAM,0CAA0C,CAAC;AAG1E,gBAAgB;AAChB,SAAS,cAAc,CACf,CAAc;AACd,sCAAsC;AACtC,UAAkB,EAClB,OAAgB,EAChB,KAAK,GAAG,CAAC,EACT,WAAW,GAAG,CAAC;IAEnB,MAAM,MAAM,GAAG,WAAW,CACtB,UAAU,CAAC,YAAY,CAAC,MAAM,EAC9B,CAAC,CACJ,CAAC;IACF,MAAM,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,GAAG,MAAM,CAAC;IAExC,gCAAgC;IAChC,MAAM,IAAI,GAAG,cAAc,CAAC,UAAU,CAAC;SAClC,GAAG,CAAC,MAAM,CAAC,EAAE,CAAC,MAAM,CAAC,YAAY,CAAC,CAAA;IAEvC,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,MAAM,CAAC;IAEtD,MAAM,OAAO,GAAG,UAAU,CAAC,CAAC,EAAE,MAAM,EAAE,mBAAmB,EAAE,KAAK,CAAC,CAAC;IAElE,MAAM,IAAI,GAAiB,EAAE,CAAC;IAC9B,MAAM,GAAG,GAAiB,EAAE,CAAC;IAC7B,KAAK,IAAI,CAAC,GAAC,CAAC,EAAE,CAAC,GAAC,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACpB,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,GAAC,CAAC,CAAC,GAAC,CAAC,CAAC,GAAC,WAAW,EAAE,MAAM,EAAE,KAAK,CAAC,CAAC,CAAC;QAC/D,GAAG,CAAC,IAAI,CAAC,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,WAAW,EAAE,KAAK,CAAC,CAAC,CAAC;IAC5D,CAAC;IAED,MAAM,MAAM,GAAG,SAAS,CAAC,CAAC,EAAE,MAAM,CAAC,MAAM,EAAE,kBAAkB,EAAE,IAAI,GAAC,WAAW,EAAE,KAAK,CAAC,CAAC;IAExF,OAAO,CAAC,GAAG,IAAI,EAAE,GAAG,GAAG,EAAE,GAAG,OAAO,EAAE,GAAG,MAAM,CAAiB,CAAC;AACpE,CAAC;AAGD,OAAO,EAAE,cAAc,EAAE,CAAA"}
|
|
@@ -1,4 +1,4 @@
|
|
|
1
|
-
import { CpNode } from '../../../cp-node/cp-node.js';
|
|
1
|
+
import type { CpNode } from '../../../cp-node/cp-node.js';
|
|
2
2
|
/** @internal */
|
|
3
3
|
declare function vertex(g: SVGGElement, cpNode: CpNode, classes?: string, delay?: number, scaleFactor?: number): SVGElement[];
|
|
4
4
|
export { vertex };
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"draw-vertex.js","sourceRoot":"","sources":["../../../../src/debug/functions/draw-elem/draw-vertex.ts"],"names":[],"mappings":"
|
|
1
|
+
{"version":3,"file":"draw-vertex.js","sourceRoot":"","sources":["../../../../src/debug/functions/draw-elem/draw-vertex.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,MAAM,EAAE,MAAM,UAAU,CAAC;AAGlC,gBAAgB;AAChB,SAAS,MAAM,CACP,CAAc,EACd,MAAc,EACd,OAAgB,EAChB,KAAK,GAAG,CAAC,EACT,WAAW,GAAG,CAAC;IAEnB,MAAM,MAAM,GAAG,MAAM,CAAC,YAAY,CAAC,MAAM,CAAC;IAE1C,IAAI,KAAK,GAAiB,EAAE,CAAC;IAE7B,+EAA+E;IAC/E,+FAA+F;IAC/F,MAAM,IAAI,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,CAAC,MAAM,EAAE,IAAI,GAAC,WAAW,EAAE,WAAW,EAAE,KAAK,CAAC,CAAC;IAEhF,KAAK,GAAG,CAAC,GAAG,IAAI,CAAC,CAAC;IAElB,OAAO,KAAK,CAAC;AACjB,CAAC;AAGD,OAAO,EAAE,MAAM,EAAE,CAAA"}
|
|
@@ -0,0 +1,32 @@
|
|
|
1
|
+
/** `2 * 2^-53` -> 2x the standard round-of unit `=== Number.EPSILON` */
|
|
2
|
+
declare const eps: number;
|
|
3
|
+
/** `2^-53` -> the standard round-of unit `=== eps/2` */
|
|
4
|
+
declare const u: number;
|
|
5
|
+
/** `2^-106` -> the standard round-of unit for double-double precision `=== (eps/2)**2` */
|
|
6
|
+
declare const uu: number;
|
|
7
|
+
/**
|
|
8
|
+
* The canonical floating point error function, γ.
|
|
9
|
+
*
|
|
10
|
+
* * roughly `=== n * (Number.EPSILON / 2)`
|
|
11
|
+
* * see e.g. [Algorithms for Accurate, Validated and Fast Polynomial Evaluation](https://hal.archives-ouvertes.fr/hal-00285603/document)
|
|
12
|
+
* @param n the parameter - typically a small positive integer, e.g. for
|
|
13
|
+
* polynomial evaluation this === 2*d + 1, where d is the degree of the
|
|
14
|
+
* polynomial
|
|
15
|
+
*
|
|
16
|
+
* @doc
|
|
17
|
+
*/
|
|
18
|
+
declare function γ(n: number): number;
|
|
19
|
+
/**
|
|
20
|
+
* The canonical, once compensated (implying double-double precision),
|
|
21
|
+
* floating point error function.
|
|
22
|
+
*
|
|
23
|
+
* * roughly `=== n * (Number.EPSILON / 2)**2`
|
|
24
|
+
* * see e.g. [Algorithms for Accurate, Validated and Fast Polynomial Evaluation](https://hal.archives-ouvertes.fr/hal-00285603/document)
|
|
25
|
+
* @param n the parameter - typically a small positive integer, e.g. for
|
|
26
|
+
* polynomial evaluation this === 2*d + 1, where d is the degree of the
|
|
27
|
+
* polynomial
|
|
28
|
+
*
|
|
29
|
+
* @doc
|
|
30
|
+
*/
|
|
31
|
+
declare function γγ(n: number): number;
|
|
32
|
+
export { γ, γγ, u, uu, eps };
|
|
@@ -0,0 +1,39 @@
|
|
|
1
|
+
/** `2 * 2^-53` -> 2x the standard round-of unit `=== Number.EPSILON` */
|
|
2
|
+
const eps = Number.EPSILON;
|
|
3
|
+
/** `2^-53` -> the standard round-of unit `=== eps/2` */
|
|
4
|
+
const u = eps / 2;
|
|
5
|
+
/** `2^-106` -> the standard round-of unit for double-double precision `=== (eps/2)**2` */
|
|
6
|
+
const uu = u * u;
|
|
7
|
+
/**
|
|
8
|
+
* The canonical floating point error function, γ.
|
|
9
|
+
*
|
|
10
|
+
* * roughly `=== n * (Number.EPSILON / 2)`
|
|
11
|
+
* * see e.g. [Algorithms for Accurate, Validated and Fast Polynomial Evaluation](https://hal.archives-ouvertes.fr/hal-00285603/document)
|
|
12
|
+
* @param n the parameter - typically a small positive integer, e.g. for
|
|
13
|
+
* polynomial evaluation this === 2*d + 1, where d is the degree of the
|
|
14
|
+
* polynomial
|
|
15
|
+
*
|
|
16
|
+
* @doc
|
|
17
|
+
*/
|
|
18
|
+
function γ(n) {
|
|
19
|
+
const nu = n * u;
|
|
20
|
+
return nu / (1 - nu);
|
|
21
|
+
}
|
|
22
|
+
/**
|
|
23
|
+
* The canonical, once compensated (implying double-double precision),
|
|
24
|
+
* floating point error function.
|
|
25
|
+
*
|
|
26
|
+
* * roughly `=== n * (Number.EPSILON / 2)**2`
|
|
27
|
+
* * see e.g. [Algorithms for Accurate, Validated and Fast Polynomial Evaluation](https://hal.archives-ouvertes.fr/hal-00285603/document)
|
|
28
|
+
* @param n the parameter - typically a small positive integer, e.g. for
|
|
29
|
+
* polynomial evaluation this === 2*d + 1, where d is the degree of the
|
|
30
|
+
* polynomial
|
|
31
|
+
*
|
|
32
|
+
* @doc
|
|
33
|
+
*/
|
|
34
|
+
function γγ(n) {
|
|
35
|
+
const nuu = n * uu;
|
|
36
|
+
return nuu / (1 - nuu);
|
|
37
|
+
}
|
|
38
|
+
export { γ, γγ, u, uu, eps };
|
|
39
|
+
//# sourceMappingURL=gamma.js.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"gamma.js","sourceRoot":"","sources":["../../src/error-analysis/gamma.ts"],"names":[],"mappings":"AACA,wEAAwE;AACxE,MAAM,GAAG,GAAG,MAAM,CAAC,OAAO,CAAC;AAE3B,wDAAwD;AACxD,MAAM,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;AAElB,0FAA0F;AAC1F,MAAM,EAAE,GAAG,CAAC,GAAC,CAAC,CAAC;AAGf;;;;;;;;;;GAUG;AACH,SAAS,CAAC,CAAC,CAAS;IAChB,MAAM,EAAE,GAAG,CAAC,GAAC,CAAC,CAAC;IAEf,OAAO,EAAE,GAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC;AACrB,CAAC;AAED;;;;;;;;;;;GAWG;AACH,SAAS,EAAE,CAAC,CAAS;IACjB,MAAM,GAAG,GAAG,CAAC,GAAC,EAAE,CAAC;IAEjB,OAAO,GAAG,GAAC,CAAC,CAAC,GAAC,GAAG,CAAC,CAAC;AACvB,CAAC;AAGD,OAAO,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,GAAG,EAAE,CAAA"}
|
|
@@ -57,10 +57,6 @@ function find2Prong(meta) {
|
|
|
57
57
|
// The lines below is an optimization.
|
|
58
58
|
const xy2 = reduceRadius(maxCoordPowerOf2, curvePieces, y, nnorm);
|
|
59
59
|
const xy = sqrt(xy2);
|
|
60
|
-
if (for1Prong && rO < (1 - oneProngTolerance) * xy) {
|
|
61
|
-
add1Prong(meta, rO, xO, yPos);
|
|
62
|
-
return undefined;
|
|
63
|
-
}
|
|
64
60
|
const l = sqrt(nnorm[0] ** 2 + nnorm[1] ** 2);
|
|
65
61
|
xO = [y[0] + nnorm[0] / l * xy, y[1] + nnorm[1] / l * xy];
|
|
66
62
|
cullCurvePieces(curvePieces, xO, xy2);
|
|
@@ -77,6 +73,10 @@ function find2Prong(meta) {
|
|
|
77
73
|
}
|
|
78
74
|
}
|
|
79
75
|
const { x, d, curve, s } = antipodalPoints[0];
|
|
76
|
+
if (for1Prong && rO < (1 - oneProngTolerance) * d) {
|
|
77
|
+
add1Prong(meta, rO, xO, yPos);
|
|
78
|
+
return undefined;
|
|
79
|
+
}
|
|
80
80
|
const z = {
|
|
81
81
|
curve, t: s, isSource: false,
|
|
82
82
|
p: toP(curve.ps, s)
|
|
@@ -90,7 +90,7 @@ function find2Prong(meta) {
|
|
|
90
90
|
}
|
|
91
91
|
// The lines below eliminate cases where the found 2-prong is too
|
|
92
92
|
// close to a one-prong
|
|
93
|
-
const squaredSeperationTolerance = max(calcSeperationTolerance(rO, sqrt(xz), 2 ** 2 * errorTolerance), minSquaredSeperationTolerance);
|
|
93
|
+
const squaredSeperationTolerance = max(calcSeperationTolerance(rO, sqrt(xz), (2 ** 2) * errorTolerance), minSquaredSeperationTolerance);
|
|
94
94
|
if (yz <= squaredSeperationTolerance) {
|
|
95
95
|
return undefined;
|
|
96
96
|
}
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"find-2-prong.js","sourceRoot":"","sources":["../../src/find-2-prong/find-2-prong.ts"],"names":[],"mappings":"AAIA,OAAO,EAAE,qBAAqB,EAAE,MAAM,gCAAgC,CAAC;AACvE,OAAO,EAAE,kBAAkB,EAAE,MAAM,sCAAsC,CAAC;AAC1E,OAAO,EAAE,YAAY,EAAE,MAAM,oBAAoB,CAAC;AAClD,OAAO,EAAE,wBAAwB,EAAE,MAAM,kCAAkC,CAAC;AAC5E,OAAO,EAAE,eAAe,EAAE,MAAM,yBAAyB,CAAC;AAC1D,OAAO,EAAE,SAAS,EAAE,MAAM,kBAAkB,CAAC;AAC7C,OAAO,EAAE,GAAG,EAAE,MAAM,2BAA2B,CAAC;AAChD,OAAO,EAAE,uBAAuB,EAAE,MAAM,gCAAgC,CAAC;AACzE,OAAO,EAAE,WAAW,EAAE,MAAM,oBAAoB,CAAC;AACjD,OAAO,EAAE,QAAQ,EAAE,MAAM,aAAa,CAAC;AACvC,OAAO,EAAE,MAAM,EAAE,MAAM,cAAc,CAAC;AAGtC,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,IAAI,CAAC;AAGrC;;;;;;;;;;;;;;;;;;;;;;GAsBG;AACH,SAAS,UAAU,CACX,IAAa;IAEjB,MAAM,EAAE,KAAK,EAAE,gBAAgB,EAAE,qBAAqB,EAAE,GAAG,IAAI,CAAC;IAEhE,8BAA8B;IAC9B,6BAA6B;IAC7B,8BAA8B;IAC9B,MAAM,6BAA6B,GAAG,CAAC,CAAC,IAAE,CAAC,gBAAgB,GAAG,EAAE,CAAC,CAAC,IAAE,CAAC,CAAC;IACtE,MAAM,cAAc,GAAG,CAAC,IAAE,CAAC,gBAAgB,GAAG,EAAE,CAAC,CAAC;IAClD,MAAM,yBAAyB,GAAG,IAAI,CAAC,qBAAqB,CAAC,CAAC;IAC9D,MAAM,iBAAiB,GAAG,CAAC,IAAE,CAAC,gBAAgB,GAAG,EAAE,CAAC,CAAC;IAErD,OAAO,SAAS,eAAe,CACvB,aAAsB,EACtB,SAAkB,EAClB,KAAa,EACb,IAAqB;QAEzB,MAAM,EAAE,KAAK,EAAE,EAAE,IAAI,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,IAAI,CAAC;QAC9C,MAAM,MAAM,GAAG,QAAQ,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;QAC/B,MAAM,MAAM,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAE,6BAA6B;QAC7E,MAAM,KAAK,GAAG,aAAa;YACvB,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,yBAAyB,CAAC;YACjC,CAAC,CAAC,CAAC,KAAK,KAAK,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;QAEtE,mCAAmC;QACnC,IAAI,CAAC,EAAE,EAAC,EAAE,CAAC,GAAG,WAAW,CACrB,KAAK,EAAE,aAAa,EAAE,yBAAyB,EAAE,IAAI,EAAE,KAAK,CAC/D,CAAC;QAGF,6DAA6D;QAC7D,yDAAyD;QACzD,MAAM,WAAW,GAAG,qBAAqB,CACrC,KAAK,EAAE,aAAa,EAAE,IAAI,EAAE,KAAK,EAAE,IAAI,EAAE,IAAI,EAAE,EAAE,CACpD,CAAC;QAEF,sCAAsC;QACtC,MAAM,GAAG,GAAG,YAAY,CAAC,gBAAgB,EAAE,WAAW,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;QAClE,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;QAErB,
|
|
1
|
+
{"version":3,"file":"find-2-prong.js","sourceRoot":"","sources":["../../src/find-2-prong/find-2-prong.ts"],"names":[],"mappings":"AAIA,OAAO,EAAE,qBAAqB,EAAE,MAAM,gCAAgC,CAAC;AACvE,OAAO,EAAE,kBAAkB,EAAE,MAAM,sCAAsC,CAAC;AAC1E,OAAO,EAAE,YAAY,EAAE,MAAM,oBAAoB,CAAC;AAClD,OAAO,EAAE,wBAAwB,EAAE,MAAM,kCAAkC,CAAC;AAC5E,OAAO,EAAE,eAAe,EAAE,MAAM,yBAAyB,CAAC;AAC1D,OAAO,EAAE,SAAS,EAAE,MAAM,kBAAkB,CAAC;AAC7C,OAAO,EAAE,GAAG,EAAE,MAAM,2BAA2B,CAAC;AAChD,OAAO,EAAE,uBAAuB,EAAE,MAAM,gCAAgC,CAAC;AACzE,OAAO,EAAE,WAAW,EAAE,MAAM,oBAAoB,CAAC;AACjD,OAAO,EAAE,QAAQ,EAAE,MAAM,aAAa,CAAC;AACvC,OAAO,EAAE,MAAM,EAAE,MAAM,cAAc,CAAC;AAGtC,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,IAAI,CAAC;AAGrC;;;;;;;;;;;;;;;;;;;;;;GAsBG;AACH,SAAS,UAAU,CACX,IAAa;IAEjB,MAAM,EAAE,KAAK,EAAE,gBAAgB,EAAE,qBAAqB,EAAE,GAAG,IAAI,CAAC;IAEhE,8BAA8B;IAC9B,6BAA6B;IAC7B,8BAA8B;IAC9B,MAAM,6BAA6B,GAAG,CAAC,CAAC,IAAE,CAAC,gBAAgB,GAAG,EAAE,CAAC,CAAC,IAAE,CAAC,CAAC;IACtE,MAAM,cAAc,GAAG,CAAC,IAAE,CAAC,gBAAgB,GAAG,EAAE,CAAC,CAAC;IAClD,MAAM,yBAAyB,GAAG,IAAI,CAAC,qBAAqB,CAAC,CAAC;IAC9D,MAAM,iBAAiB,GAAG,CAAC,IAAE,CAAC,gBAAgB,GAAG,EAAE,CAAC,CAAC;IAErD,OAAO,SAAS,eAAe,CACvB,aAAsB,EACtB,SAAkB,EAClB,KAAa,EACb,IAAqB;QAEzB,MAAM,EAAE,KAAK,EAAE,EAAE,IAAI,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,IAAI,CAAC;QAC9C,MAAM,MAAM,GAAG,QAAQ,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;QAC/B,MAAM,MAAM,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAE,6BAA6B;QAC7E,MAAM,KAAK,GAAG,aAAa;YACvB,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,yBAAyB,CAAC;YACjC,CAAC,CAAC,CAAC,KAAK,KAAK,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;QAEtE,mCAAmC;QACnC,IAAI,CAAC,EAAE,EAAC,EAAE,CAAC,GAAG,WAAW,CACrB,KAAK,EAAE,aAAa,EAAE,yBAAyB,EAAE,IAAI,EAAE,KAAK,CAC/D,CAAC;QAGF,6DAA6D;QAC7D,yDAAyD;QACzD,MAAM,WAAW,GAAG,qBAAqB,CACrC,KAAK,EAAE,aAAa,EAAE,IAAI,EAAE,KAAK,EAAE,IAAI,EAAE,IAAI,EAAE,EAAE,CACpD,CAAC;QAEF,sCAAsC;QACtC,MAAM,GAAG,GAAG,YAAY,CAAC,gBAAgB,EAAE,WAAW,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;QAClE,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;QAErB,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,IAAE,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,IAAE,CAAC,CAAC,CAAC;QAC1C,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAC,CAAC,GAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAC,CAAC,GAAC,EAAE,CAAC,CAAC;QAElD,eAAe,CAAC,WAAW,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC;QAEtC,MAAM,eAAe,GAAG,kBAAkB,CACtC,gBAAgB,EAAE,KAAK,EAAE,IAAI,EAC7B,SAAS,EAAE,KAAK,EAChB,WAAW,CACd,CAAC;QAEF,IAAI,eAAe,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;YAC/B,IAAI,eAAe,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;gBAAC,OAAO,SAAS,CAAC;YAAC,CAAC;YACvD,IAAI,CAAC,mBAAmB,CAAC,eAAe,CAAC,EAAE,CAAC;gBACxC,sEAAsE;gBACtE,kEAAkE;gBAClE,oBAAoB;gBACpB,OAAO,SAAS,CAAC;YACrB,CAAC;QACL,CAAC;QAED,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC,EAAE,GAAG,eAAe,CAAC,CAAC,CAAC,CAAC;QAE9C,IAAI,SAAS,IAAI,EAAE,GAAG,CAAC,CAAC,GAAG,iBAAiB,CAAC,GAAC,CAAC,EAAE,CAAC;YAC9C,SAAS,CAAC,IAAI,EAAE,EAAE,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;YAC9B,OAAO,SAAS,CAAC;QACrB,CAAC;QAED,MAAM,CAAC,GAAoB;YACvB,KAAK,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,EAAE,KAAK;YAC5B,CAAC,EAAE,GAAG,CAAC,KAAK,CAAC,EAAE,EAAE,CAAC,CAAC;SACtB,CAAC;QAEF,MAAM,MAAM,GAAG,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,CAAC;QAExC,MAAM,EAAE,GAAG,wBAAwB,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QACjD,MAAM,EAAE,GAAG,wBAAwB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAE5C,IAAI,CAAC,aAAa,EAAE,CAAC;YACjB,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;gBAAC,OAAO,SAAS,CAAC;YAAC,CAAC;YAEnC,iEAAiE;YACjE,uBAAuB;YACvB,MAAM,0BAA0B,GAAG,GAAG,CAClC,uBAAuB,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,IAAE,CAAC,CAAC,GAAC,cAAc,CAAC,EAC5D,6BAA6B,CAChC,CAAC;YACF,IAAI,EAAE,IAAI,0BAA0B,EAAE,CAAC;gBACnC,OAAO,SAAS,CAAC;YACrB,CAAC;QACL,CAAC;QAED,OAAO,EAAE,MAAM,EAAE,CAAC,EAAE,CAAC;IACzB,CAAC,CAAA;AACL,CAAC;AAGD,gBAAgB;AAChB,SAAS,mBAAmB,CACpB,GAAqB;IAEzB,MAAM,EAAE,KAAK,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IAC/B,KAAK,IAAI,CAAC,GAAC,CAAC,EAAE,CAAC,GAAC,GAAG,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAC9B,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QAExC,IAAI,KAAK,KAAK,MAAM,IAAI,CAAC,KAAK,EAAE,EAAE,CAAC;YAC/B,OAAO,KAAK,CAAC;QACjB,CAAC;IACL,CAAC;IAED,OAAO,IAAI,CAAC;AAChB,CAAC;AAGD,OAAO,EAAE,UAAU,EAAE,CAAA"}
|