flo-mat 3.0.14 → 4.0.1

package/src/cp-node/fs/get-smoothed-speed.ts

@@ -1,97 +0,0 @@
-import { memoize } from "flo-memoize";
-import { totalCurvature$ } from "../../bezier/total-curvature.js";
-import { totalLength$ } from "../../bezier/total-length.js";
-import { memoizePrimitive } from "../../memoize-primitive.js";
-import { CpNode } from "../cp-node.js";
-import { getMatCurveToNext$ } from "./get-mat-curve-to-next.js";
-import { getRealProngCount } from "./get-real-prong-count.js";
-import { isFullyTerminating } from "./is-fully-terminating.js";
-
-
-const RLen = 40;
-
-const getSmoothedSpeed$ = memoizePrimitive(
-    function (L: number) {
-        return memoize(
-            function getSmoothedSpeed(
-                    cpNode: CpNode): number | undefined {
-
-                let cpStart = cpNode;
-
-                let cumLB = 0;
-                let cumLF = 0;
-                let curB = cpStart.prev;
-                let curF = cpStart;
-                let cumKB = 0;
-                let cumKF = 0;
-                let cumLRB = 0;
-                let cumLRF = 0;
-                while (cumLRF < L/2/RLen &&
-                       curF !== cpStart.prev &&
-                       !isFullyTerminating(curF) &&
-                       getRealProngCount(curF) <= 2) {
-                    const { radius: r } = curF.cp.circle;
-                    const ps = getMatCurveToNext$(curF);
-                    const turn = totalCurvature$(ps);
-
-                    let rot = 0;
-                    // if (getProngCount(curF) >= 3) {
-                    //     rot = -getInterfaceRotation(
-                    //         tangentAt0(ps),
-                    //         tangentAt1(curF.prev.cp.pointOnShape.curve.ps)
-                    //     );
-                    // }
-
-                    const l = totalLength$(getMatCurveToNext$(curF));
-                    cumKF += r*(turn + rot);
-                    cumLF += l;
-                    cumLRF += l/r;
-
-                    curF = curF.next;
-                }
-                while (cumLRB < L/2/RLen &&
-                       curB !== cpStart &&
-                       !isFullyTerminating(curB) &&
-                       getRealProngCount(curB) <= 2) {
-                    const { radius: r } = curB.cp.circle;
-                    const ps = getMatCurveToNext$(curB);
-                    const turn = totalCurvature$(ps);
-
-                    let rot = 0;
-                    // if (getProngCount(curB) >= 3) {
-                    //     rot = -getInterfaceRotation(
-                    //         tangentAt0(ps),
-                    //         tangentAt1(curB.prev.cp.pointOnShape.curve.ps)
-                    //     )
-                    // }
-
-                    const l = totalLength$(getMatCurveToNext$(curB));
-                    cumKB += r*(turn + rot);
-                    cumLB += l;
-                    cumLRB += l/r;
-
-                    curB = curB.prev;
-                }
-                const cumL = cumLB + cumLF;
-                const cumK = cumKB + cumKF;
-
-                if (cumL === 0) {
-                    return undefined;
-                }
-
-                const x = -cumK/cumL;
-
-                // f(x) = 1 - exp(-x)
-                // f(x) = 1 - 1/(1 + x)
-                // f(x) = x/(a + x)
-                // f(x) = tanh(x)
-
-                // return 1 - Math.E**(-x)
-                return x;
-            }
-        )
-    }
-);
-
-
-export { getSmoothedSpeed$ }

package/src/cp-node/fs/get-smoothed.ts

@@ -1,62 +0,0 @@
-import { totalAbsoluteCurvature, totalCurvature, totalLength } from "flo-bezier3";
-import { memoize } from "flo-memoize";
-import { totalCurvature$ } from "../../bezier/total-curvature.js";
-import { totalLength$ } from "../../bezier/total-length.js";
-import { memoizePrimitive } from "../../memoize-primitive.js";
-import { CpNode } from "../cp-node.js";
-import { getAvgSpeed, getAvgSpeed$ } from "./get-avg-speed.js";
-import { getMatCurveToNext, getMatCurveToNext$ } from "./get-mat-curve-to-next.js";
-import { getSpeed } from "./get-speed.js";
-
-
-const getSmoothed$ = memoizePrimitive(
-    function (L: number) {
-        return memoize(
-            function getSmoothed(
-                    cpNode: CpNode) {
-
-                // return getAvg$(cpNode)
-                // return getAvgCurvature(cpNode)
-                let cpStart = cpNode;
-
-                let cumLB = 0;
-                let cumLF = 0;
-                let curB = cpStart.prev;
-                let curF = cpStart;
-                let cumRTurnB = 0;
-                let cumRTurnF = 0;
-                while (cumLB < L/2 && curB !== cpStart) {
-                    const { radius: r } = curB.cp.circle;
-                    const ps = getMatCurveToNext$(curB);
-                    const absK = totalCurvature$(ps);
-                    const l = totalLength$(getMatCurveToNext$(curB));
-                    cumRTurnB += r*absK;
-                    cumLB += l;
-
-                    curB = curB.prev;
-                }
-                while (cumLF < L/2 && curF !== cpStart.prev) {
-                    const { radius: r } = curF.cp.circle;
-                    const ps = getMatCurveToNext$(curF);
-                    const absK = totalCurvature$(ps);
-                    const l = totalLength$(getMatCurveToNext$(curF));
-                    cumRTurnF += r*absK;
-                    cumLF += l;
-
-                    curF = curF.next;
-                }
-                const cumL = cumLB + cumLF;
-                const cumRTurn = cumRTurnB + cumRTurnF;
-
-                if (cumL === 0) {
-                    return undefined;
-                }
-
-                return { cumRTurn, cumL };
-            }
-        )
-    }
-);
-
-
-export { getSmoothed$ }

package/src/cp-node/fs/get-speed.ts

@@ -1,56 +0,0 @@
-import { totalAbsoluteCurvature, totalLength } from "flo-bezier3";
-import { memoize } from "flo-memoize";
-import { CpNode } from "../cp-node.js";
-import { getMatCurveToNext } from "./get-mat-curve-to-next.js";
-import { getRealProngCount } from "./get-real-prong-count.js";
-
-const getSpeed$ = memoize(getSpeed);
-
-
-function getSpeed(cpNode: CpNode): number | undefined {
-    if (getRealProngCount(cpNode) !== 2) {
-        return 0;
-    }
-
-    const { radius: r } = cpNode.cp.circle;
-    const matPs = getMatCurveToNext(cpNode);
-    const absK = totalAbsoluteCurvature(matPs, [0,1]);
-    const l = totalLength(matPs);
-    if (l === 0) {
-        return undefined;  // TODO2
-    }
-    // const avgK = absK/l;
-    const avgK = absK;
-
-    return avgK*r;
-}
-
-
-
-
-// function getSpeed() {
-    // const cpNodeF = cpNode;
-    // const cpNodeB = cpNode.next.prevOnCircle;
-    // const matPs = getMatCurveToNext(cpNodeF);
-    // const boundaryPsF = getBoundaryBeziersToNext(cpNodeF);
-    // const boundaryPsB = getBoundaryBeziersToNext(cpNodeB);
-
-    // const lF = sum(boundaryPsF.map(ps => totalLength(ps)));
-    // const lB = sum(boundaryPsB.map(ps => totalLength(ps)));
-
-    // const l = totalLength(matPs);
-    // const ll = l === 0 ? 0.001 : l;
-
-    // const _speedF = lF/ll;
-    // const _speedB = lB/ll;
-
-    // const speedF = _speedF < 1 ? _speedF : 1;
-    // const speedB = _speedB < 1 ? _speedB : 1;
-
-
-    // const speed = min(speedF, speedB);
-
-    // return speed;
-// }
-
-export { getSpeed, getSpeed$ }

package/src/curve/curve.ts

@@ -1,37 +0,0 @@
-import { Loop } from 'flo-boolean';
-
-
-/**
- * Represents a bezier curve on the shape boundary / loop.
- */
-interface Curve {
-    /**
-     * The closed loop of bezier curves representing the shape 
-     * boundary this curve belongs to.
-     */
-    loop: Loop;
-    /**
-     * The bezier control points.
-     */
-    ps: number[][];
-    /**
-     * The previous curve (when going in a negative direction around
-     * the shape boundary, i.e. clockwise for the outer shape and anti-clockwise
-     * for the holes (if any))
-     */
-    prev: Curve; 
-    /**
-     * The next curve (when going in a positive direction around
-     * the shape boundary, i.e. anti-clockwise for the outer shape and clockwise
-     * for the holes (if any))
-     */
-    next: Curve;
-    /**
-     * The curve's ordered index in the loop. This imposes a cycling
-     * ordering of the curves in the loop
-     */
-    idx: number;
-}
-
-
-export { Curve }

package/src/dbscan/dbscan.ts

@@ -1,91 +0,0 @@
-
-/**
- * DBSCAN Clustering Algorithm
- * 
- * @param DB 
- * @param distFunc 
- * @param eps 
- * @param minPts 
- */
-function DBSCAN<T>(
-        DB: T[],
-        distFunc: (a: T, b: T) => number,
-        eps: number,
-        minPts: number) {
-
-    const labelMap = new Map<T,number>();
-    let C = 0;                                                  /* Cluster counter */
-    for (const P of DB) {
-        if (labelMap.has(P)) { continue; }                      /* Previously processed in inner loop */
-
-        const N = RangeQuery(DB, distFunc, P, eps);             /* Find neighbors */
-        if (N.size < minPts) {                                  /* Density check */
-            labelMap.set(P,-1);                                 /* Label as Noise */
-            continue;
-        }
-
-        C++;                                                    /* next cluster label */
-        labelMap.set(P,C);                                      /* Label initial point */
-        const S: Set<T> = new Set(N); S.delete(P);              /* Neighbors to expand */
-        const Ss = Array.from(S);
-
-        for (let i=0; i<Ss.length; i++) {                       /* Process every seed point Q */
-            const Q = Ss[i];
-            if (labelMap.get(Q) === -1) { labelMap.set(Q,C); }  /* Change Noise to border point */
-            if (labelMap.get(Q) !== undefined) { continue; }    /* Previously processed (e.g., border point) */
-            labelMap.set(Q,C);                                  /* Label neighbor */
-            const N = RangeQuery(DB, distFunc, Q, eps);         /* Find neighbors */
-            if (N.size >= minPts) {                             /* Density check (if Q is a core point) */
-                Ss.push(...N);                                  /* Add new neighbors to seed set */
-            }
-        }
-    }
-
-    return labelMap;
-}
-
-
-function RangeQuery<T>(
-        DB: T[],
-        distFunc: (a: T, b: T) => number,
-        Q: T,
-        eps: number) {
-
-    let N = new Set<T>();
-    for (const P of DB) {              /* Scan all points in the database */
-        if (distFunc(Q, P) <= eps) {   /* Compute distance and check epsilon */
-            N.add(P);                  /* Add to result */
-        }
-    }
-
-    return N;
-}
-
-
-export { DBSCAN }
-
-
-// function makeMore(p: number[]) {
-//     const ps: number[][] = [];
-//     for (let i=0; i<5; i++) {
-//         const [x,y] = p;
-//         ps.push([
-//             x + 5*(Math.random() - 0.5),
-//             y + 5*(Math.random() - 0.5),
-//         ]);
-//     }
-
-//     return ps;
-// }
-// const ps = [
-//     [10,10],[10,30],[20,50],[10,40],[70,50],[30,30]
-// ]
-// .flatMap(makeMore);
-// function d(a: number[], b: number[]) {
-//     const x = a[0] - b[0];
-//     const y = a[1] - b[1];
-
-//     return Math.sqrt(x**2 + y**2);
-// }
-// DBSCAN(ps, d, 5, 3);
-

package/src/debug/cp-node-for-debugging.ts

@@ -1,15 +0,0 @@
-import { CpNode } from '../cp-node/cp-node.js';
-import { Generated } from './debug.js';
-
-
-/**
- * @internal
- * Used for debugging only. 
- */
-interface ICpNodeForDebugging {
-    generated: Generated,
-    cpNode: CpNode
-}
-
-
-export { ICpNodeForDebugging }

package/src/debug/functions/draw-elem/draw-speed.ts

@@ -1,35 +0,0 @@
-import { normal } from 'flo-bezier3';
-import { drawFs } from 'flo-draw';
-import { toLength, translate } from 'flo-vector2d';
-import { CpNode } from '../../../cp-node/cp-node.js';
-import { getSmoothedSpeed$ } from '../../../cp-node/fs/get-smoothed-speed.js';
-
-
-/** @internal */
-function drawSpeed(
-        g: SVGGElement,
-        cpNode: CpNode,
-        classes_?: string,
-        delay = 0,
-        scaleFactor = 1) {    
-
-    const $elems: SVGElement[] = [];
-
-    const s = getSmoothedSpeed$(0)(cpNode);
-
-    if (s === undefined) { return []; }
-
-    const pos = cpNode.cp.pointOnShape;
-    const { p, t } = pos;
-    const ps = pos.curve.ps;
-
-    const ss = s*scaleFactor/8;
-
-    const n = toLength(normal(ps,t), ss);
-    $elems.push(...drawFs.line(g, [p,translate(p,n)], 'red thin5', delay));
-    
-    return $elems;
-}
-
-
-export { drawSpeed }

package/src/debug/three-prong-for-debugging.ts

@@ -1,58 +0,0 @@
-import { CpNode } from '../cp-node/cp-node.js';
-import { PointOnShape } from '../point-on-shape/point-on-shape.js';
-import { posToHumanString } from '../point-on-shape/pos-to-human-string.js';
-import { Circle } from '../geometry/circle.js';
-import { Generated } from './debug.js';
-
-
-/**
- * @internal
- * Used for debugging only. A three-prong is a maximally inscribed circle that
- * touches the shape boundary (tangentially) at 3 points.
- */
-interface ThreeProngForDebugging {
-    generated  : Generated;
-    /** The best fit circle found for the 3-prong. */
-    circle     : Circle;
-    /** The best fit 3 points found for the 3-prong. */
-    poss       : PointOnShape[];
-    /** The 3 boundary pieces on which the three prong points were found. */
-    cp3ss      : CpNode[][];
-    /** The boundary pieces that were used to search the three prong on. */
-    cpss       : CpNode[][];                 
-    bestIndx   : number; 
-    /** An array of 3-prongs, each of which may be a best fit 3-prong. */
-    candidateThreeProngs : { 
-        circle: Circle; 
-        ps: PointOnShape[]; 
-    }[];
-    visitedCps : CpNode[];
-    boundaries : number[][][][];
-    traces     : number[][][];
-}
-
-
-function createEmptyThreeProngForDebugging(): ThreeProngForDebugging {
-    return {
-        generated  : undefined!,
-        circle     : undefined!,
-        poss       : undefined!,
-        cp3ss      : undefined!,
-        cpss       : undefined!,
-        bestIndx   : undefined!,
-        candidateThreeProngs : undefined!,
-        visitedCps : undefined!,
-        boundaries : undefined!,
-        traces     : undefined!,
-    }
-}
-
-function cpsSimple(threeProngForDebugging: ThreeProngForDebugging) {
-    return threeProngForDebugging.cpss.map(δ =>
-        [posToHumanString( δ[0].cp.pointOnShape ),
-         posToHumanString( δ[1].cp.pointOnShape )] 
-    );
-}
-
-
-export { ThreeProngForDebugging, createEmptyThreeProngForDebugging }

package/src/debug/two-prong-for-debugging.ts

@@ -1,30 +0,0 @@
-// import { CpNode } from '../cp-node/cp-node.js';
-// import { PointOnShape } from '../point-on-shape/point-on-shape.js';
-// import { Circle } from '../geometry/circle.js';
-// import { CurvePiece } from '../mat/bezier-piece.js';
-// import { Generated } from './debug.js';
-
-
-// /** @internal */
-// interface TwoProngForDebugging {
-//     generated: Generated;
-//     bezierPieces: CurvePiece[];
-//     pos: PointOnShape;
-//     δ: CpNode[]; 
-//     z: number[];
-//     circle: Circle; 
-//     xs: { 
-//         x: number[]; 
-//         y: PointOnShape; 
-//         z: PointOnShape; 
-//         t: number; 
-//     }[]; 
-//     failed: boolean; 
-//     holeClosing: boolean;
-//     notAdded: boolean;
-//     deleted: boolean;
-//     cpNode: CpNode;
-// }
-
-
-// export { TwoProngForDebugging }