@twick/2d 0.11.0

package/lib/curves/Polynomial.js

@@ -0,0 +1,259 @@
+import { clamp } from '@twick/core';
+/**
+ * A polynomial in the form ax^3 + bx^2 + cx + d up to a cubic polynomial.
+ *
+ * Source code liberally taken from:
+ * https://github.com/FreyaHolmer/Mathfs/blob/master/Runtime/Curves/Polynomial.cs
+ */
+export class Polynomial {
+    /**
+     * Constructs a constant polynomial
+     *
+     * @param c0 - The constant coefficient
+     */
+    static constant(c0) {
+        return new Polynomial(c0);
+    }
+    /**
+     * Constructs a linear polynomial
+     *
+     * @param c0 - The constant coefficient
+     * @param c1 - The linear coefficient
+     */
+    static linear(c0, c1) {
+        return new Polynomial(c0, c1);
+    }
+    /**
+     * Constructs a quadratic polynomial
+     *
+     * @param c0 - The constant coefficient
+     * @param c1 - The linear coefficient
+     * @param c2 - The quadratic coefficient
+     */
+    static quadratic(c0, c1, c2) {
+        return new Polynomial(c0, c1, c2);
+    }
+    /**
+     * Constructs a cubic polynomial
+     *
+     * @param c0 - The constant coefficient
+     * @param c1 - The linear coefficient
+     * @param c2 - The quadratic coefficient
+     * @param c3 - The cubic coefficient
+     */
+    static cubic(c0, c1, c2, c3) {
+        return new Polynomial(c0, c1, c2, c3);
+    }
+    /**
+     * The degree of the polynomial
+     */
+    get degree() {
+        if (this.c3 !== 0) {
+            return 3;
+        }
+        else if (this.c2 !== 0) {
+            return 2;
+        }
+        else if (this.c1 !== 0) {
+            return 1;
+        }
+        return 0;
+    }
+    constructor(c0, c1, c2, c3) {
+        this.c0 = c0;
+        this.c1 = c1 ?? 0;
+        this.c2 = c2 ?? 0;
+        this.c3 = c3 ?? 0;
+    }
+    /**
+     * Return the nth derivative of the polynomial.
+     *
+     * @param n - The number of times to differentiate the polynomial.
+     */
+    differentiate(n = 1) {
+        switch (n) {
+            case 0:
+                return this;
+            case 1:
+                return new Polynomial(this.c1, 2 * this.c2, 3 * this.c3, 0);
+            case 2:
+                return new Polynomial(2 * this.c2, 6 * this.c3, 0, 0);
+            case 3:
+                return new Polynomial(6 * this.c3, 0, 0, 0);
+            default:
+                throw new Error('Unsupported derivative');
+        }
+    }
+    eval(t, derivative = 0) {
+        if (derivative !== 0) {
+            return this.differentiate(derivative).eval(t);
+        }
+        return this.c3 * (t * t * t) + this.c2 * (t * t) + this.c1 * t + this.c0;
+    }
+    /**
+     * Split the polynomial into two polynomials of the same overall shape.
+     *
+     * @param u - The point at which to split the polynomial.
+     */
+    split(u) {
+        const d = 1 - u;
+        const pre = new Polynomial(this.c0, this.c1 * u, this.c2 * u * u, this.c3 * u * u * u);
+        const post = new Polynomial(this.eval(0), d * this.differentiate(1).eval(u), ((d * d) / 2) * this.differentiate(2).eval(u), ((d * d * d) / 6) * this.differentiate(3).eval(u));
+        return [pre, post];
+    }
+    /**
+     * Calculate the roots (values where this polynomial = 0).
+     *
+     * @remarks
+     * Depending on the degree of the polynomial, returns between 0 and 3 results.
+     */
+    roots() {
+        switch (this.degree) {
+            case 3:
+                return this.solveCubicRoots();
+            case 2:
+                return this.solveQuadraticRoots();
+            case 1:
+                return this.solveLinearRoot();
+            case 0:
+                return [];
+            default:
+                throw new Error(`Unsupported polynomial degree: ${this.degree}`);
+        }
+    }
+    /**
+     * Calculate the local extrema of the polynomial.
+     */
+    localExtrema() {
+        return this.differentiate().roots();
+    }
+    /**
+     * Calculate the local extrema of the polynomial in the unit interval.
+     */
+    localExtrema01() {
+        const all = this.localExtrema();
+        const valids = [];
+        for (let i = 0; i < all.length; i++) {
+            const t = all[i];
+            if (t >= 0 && t <= 1) {
+                valids.push(all[i]);
+            }
+        }
+        return valids;
+    }
+    /**
+     * Return the output value range within the unit interval.
+     */
+    outputRange01() {
+        let range = [this.eval(0), this.eval(1)];
+        // Expands the minimum or maximum value of the range to contain the given
+        // value.
+        const encapsulate = (value) => {
+            if (range[1] > range[0]) {
+                range = [Math.min(range[0], value), Math.max(range[1], value)];
+            }
+            else {
+                range = [Math.min(range[1], value), Math.max(range[0], value)];
+            }
+        };
+        this.localExtrema01().forEach(t => encapsulate(this.eval(t)));
+        return range;
+    }
+    solveCubicRoots() {
+        const a = this.c0;
+        const b = this.c1;
+        const c = this.c2;
+        const d = this.c3;
+        // First, depress the cubic to make it easier to solve
+        const aa = a * a;
+        const ac = a * c;
+        const bb = b * b;
+        const p = (3 * ac - bb) / (3 * aa);
+        const q = (2 * bb * b - 9 * ac * b + 27 * aa * d) / (27 * aa * a);
+        const dpr = this.solveDepressedCubicRoots(p, q);
+        // We now have the roots of the depressed cubic, now convert back to the
+        // normal cubic
+        const undepressRoot = (r) => r - b / (3 * a);
+        switch (dpr.length) {
+            case 1:
+                return [undepressRoot(dpr[0])];
+            case 2:
+                return [undepressRoot(dpr[0]), undepressRoot(dpr[1])];
+            case 3:
+                return [
+                    undepressRoot(dpr[0]),
+                    undepressRoot(dpr[1]),
+                    undepressRoot(dpr[2]),
+                ];
+            default:
+                return [];
+        }
+    }
+    solveDepressedCubicRoots(p, q) {
+        // t³+pt+q = 0
+        // Triple root - one solution. solve x³+q = 0 => x = cr(-q)
+        if (this.almostZero(p)) {
+            return [Math.cbrt(-q)];
+        }
+        const TAU = Math.PI * 2;
+        const discriminant = 4 * p * p * p + 27 * q * q;
+        if (discriminant < 0.00001) {
+            // Two or three roots guaranteed, use trig solution
+            const pre = 2 * Math.sqrt(-p / 3);
+            const acosInner = ((3 * q) / (2 * p)) * Math.sqrt(-3 / p);
+            const getRoot = (k) => pre *
+                Math.cos((1 / 3) * Math.acos(clamp(-1, 1, acosInner)) - (TAU / 3) * k);
+            // If acos hits 0 or TAU/2, the offsets will have the same value,
+            // which means we have a double root plus one regular root on our hands
+            if (acosInner >= 0.9999) {
+                // two roots - one single and one double root
+                return [getRoot(0), getRoot(2)];
+            }
+            if (acosInner <= -0.9999) {
+                // two roots - one single and one double root
+                return [getRoot(1), getRoot(2)];
+            }
+            return [getRoot(0), getRoot(1), getRoot(2)];
+        }
+        if (discriminant > 0 && p < 0) {
+            // one root
+            const coshInner = (1 / 3) *
+                Math.acosh(((-3 * Math.abs(q)) / (2 * p)) * Math.sqrt(-3 / p));
+            const r = -2 * Math.sign(q) * Math.sqrt(-p / 3) * Math.cosh(coshInner);
+            return [r];
+        }
+        if (p > 0) {
+            // one root
+            const sinhInner = (1 / 3) * Math.asinh(((3 * q) / (2 * p)) * Math.sqrt(3 / p));
+            const r = -2 * Math.sqrt(p / 3) * Math.sinh(sinhInner);
+            return [r];
+        }
+        // no roots
+        return [];
+    }
+    solveQuadraticRoots() {
+        const a = this.c2;
+        const b = this.c1;
+        const c = this.c0;
+        const rootContent = b * b - 4 * a * c;
+        if (this.almostZero(rootContent)) {
+            // two equivalent solutions at one point
+            return [-b / (2 * a)];
+        }
+        if (rootContent >= 0) {
+            const root = Math.sqrt(rootContent);
+            // crosses at two points
+            const r0 = (-b - root) / (2 * a);
+            const r1 = (-b + root) / (2 * a);
+            return [Math.min(r0, r1), Math.max(r0, r1)];
+        }
+        return [];
+    }
+    solveLinearRoot() {
+        return [-this.c0 / this.c1];
+    }
+    almostZero(value) {
+        return Math.abs(0 - value) <= Number.EPSILON;
+    }
+}
+//# sourceMappingURL=data:application/json;base64,eyJ2ZXJzaW9uIjozLCJmaWxlIjoiUG9seW5vbWlhbC5qcyIsInNvdXJjZVJvb3QiOiIiLCJzb3VyY2VzIjpbIi4uLy4uL3NyYy9saWIvY3VydmVzL1BvbHlub21pYWwudHMiXSwibmFtZXMiOltdLCJtYXBwaW5ncyI6IkFBQUEsT0FBTyxFQUFDLEtBQUssRUFBQyxNQUFNLGFBQWEsQ0FBQztBQUVsQzs7Ozs7R0FLRztBQUNILE1BQU0sT0FBTyxVQUFVO0lBS3JCOzs7O09BSUc7SUFDSSxNQUFNLENBQUMsUUFBUSxDQUFDLEVBQVU7UUFDL0IsT0FBTyxJQUFJLFVBQVUsQ0FBQyxFQUFFLENBQUMsQ0FBQztJQUM1QixDQUFDO0lBRUQ7Ozs7O09BS0c7SUFDSSxNQUFNLENBQUMsTUFBTSxDQUFDLEVBQVUsRUFBRSxFQUFVO1FBQ3pDLE9BQU8sSUFBSSxVQUFVLENBQUMsRUFBRSxFQUFFLEVBQUUsQ0FBQyxDQUFDO0lBQ2hDLENBQUM7SUFFRDs7Ozs7O09BTUc7SUFDSSxNQUFNLENBQUMsU0FBUyxDQUFDLEVBQVUsRUFBRSxFQUFVLEVBQUUsRUFBVTtRQUN4RCxPQUFPLElBQUksVUFBVSxDQUFDLEVBQUUsRUFBRSxFQUFFLEVBQUUsRUFBRSxDQUFDLENBQUM7SUFDcEMsQ0FBQztJQUVEOzs7Ozs7O09BT0c7SUFDSSxNQUFNLENBQUMsS0FBSyxDQUNqQixFQUFVLEVBQ1YsRUFBVSxFQUNWLEVBQVUsRUFDVixFQUFVO1FBRVYsT0FBTyxJQUFJLFVBQVUsQ0FBQyxFQUFFLEVBQUUsRUFBRSxFQUFFLEVBQUUsRUFBRSxFQUFFLENBQUMsQ0FBQztJQUN4QyxDQUFDO0lBRUQ7O09BRUc7SUFDSCxJQUFXLE1BQU07UUFDZixJQUFJLElBQUksQ0FBQyxFQUFFLEtBQUssQ0FBQyxFQUFFLENBQUM7WUFDbEIsT0FBTyxDQUFDLENBQUM7UUFDWCxDQUFDO2FBQU0sSUFBSSxJQUFJLENBQUMsRUFBRSxLQUFLLENBQUMsRUFBRSxDQUFDO1lBQ3pCLE9BQU8sQ0FBQyxDQUFDO1FBQ1gsQ0FBQzthQUFNLElBQUksSUFBSSxDQUFDLEVBQUUsS0FBSyxDQUFDLEVBQUUsQ0FBQztZQUN6QixPQUFPLENBQUMsQ0FBQztRQUNYLENBQUM7UUFDRCxPQUFPLENBQUMsQ0FBQztJQUNYLENBQUM7SUF3QkQsWUFDa0IsRUFBVSxFQUMxQixFQUFXLEVBQ1gsRUFBVyxFQUNYLEVBQVc7UUFISyxPQUFFLEdBQUYsRUFBRSxDQUFRO1FBSzFCLElBQUksQ0FBQyxFQUFFLEdBQUcsRUFBRSxJQUFJLENBQUMsQ0FBQztRQUNsQixJQUFJLENBQUMsRUFBRSxHQUFHLEVBQUUsSUFBSSxDQUFDLENBQUM7UUFDbEIsSUFBSSxDQUFDLEVBQUUsR0FBRyxFQUFFLElBQUksQ0FBQyxDQUFDO0lBQ3BCLENBQUM7SUFFRDs7OztPQUlHO0lBQ0ksYUFBYSxDQUFDLENBQUMsR0FBRyxDQUFDO1FBQ3hCLFFBQVEsQ0FBQyxFQUFFLENBQUM7WUFDVixLQUFLLENBQUM7Z0JBQ0osT0FBTyxJQUFJLENBQUM7WUFDZCxLQUFLLENBQUM7Z0JBQ0osT0FBTyxJQUFJLFVBQVUsQ0FBQyxJQUFJLENBQUMsRUFBRSxFQUFFLENBQUMsR0FBRyxJQUFJLENBQUMsRUFBRSxFQUFFLENBQUMsR0FBRyxJQUFJLENBQUMsRUFBRSxFQUFFLENBQUMsQ0FBQyxDQUFDO1lBQzlELEtBQUssQ0FBQztnQkFDSixPQUFPLElBQUksVUFBVSxDQUFDLENBQUMsR0FBRyxJQUFJLENBQUMsRUFBRSxFQUFFLENBQUMsR0FBRyxJQUFJLENBQUMsRUFBRSxFQUFFLENBQUMsRUFBRSxDQUFDLENBQUMsQ0FBQztZQUN4RCxLQUFLLENBQUM7Z0JBQ0osT0FBTyxJQUFJLFVBQVUsQ0FBQyxDQUFDLEdBQUcsSUFBSSxDQUFDLEVBQUUsRUFBRSxDQUFDLEVBQUUsQ0FBQyxFQUFFLENBQUMsQ0FBQyxDQUFDO1lBQzlDO2dCQUNFLE1BQU0sSUFBSSxLQUFLLENBQUMsd0JBQXdCLENBQUMsQ0FBQztRQUM5QyxDQUFDO0lBQ0gsQ0FBQztJQWVNLElBQUksQ0FBQyxDQUFTLEVBQUUsVUFBVSxHQUFHLENBQUM7UUFDbkMsSUFBSSxVQUFVLEtBQUssQ0FBQyxFQUFFLENBQUM7WUFDckIsT0FBTyxJQUFJLENBQUMsYUFBYSxDQUFDLFVBQVUsQ0FBQyxDQUFDLElBQUksQ0FBQyxDQUFDLENBQUMsQ0FBQztRQUNoRCxDQUFDO1FBQ0QsT0FBTyxJQUFJLENBQUMsRUFBRSxHQUFHLENBQUMsQ0FBQyxHQUFHLENBQUMsR0FBRyxDQUFDLENBQUMsR0FBRyxJQUFJLENBQUMsRUFBRSxHQUFHLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxHQUFHLElBQUksQ0FBQyxFQUFFLEdBQUcsQ0FBQyxHQUFHLElBQUksQ0FBQyxFQUFFLENBQUM7SUFDM0UsQ0FBQztJQUVEOzs7O09BSUc7SUFDSSxLQUFLLENBQUMsQ0FBUztRQUNwQixNQUFNLENBQUMsR0FBRyxDQUFDLEdBQUcsQ0FBQyxDQUFDO1FBRWhCLE1BQU0sR0FBRyxHQUFHLElBQUksVUFBVSxDQUN4QixJQUFJLENBQUMsRUFBRSxFQUNQLElBQUksQ0FBQyxFQUFFLEdBQUcsQ0FBQyxFQUNYLElBQUksQ0FBQyxFQUFFLEdBQUcsQ0FBQyxHQUFHLENBQUMsRUFDZixJQUFJLENBQUMsRUFBRSxHQUFHLENBQUMsR0FBRyxDQUFDLEdBQUcsQ0FBQyxDQUNwQixDQUFDO1FBQ0YsTUFBTSxJQUFJLEdBQUcsSUFBSSxVQUFVLENBQ3pCLElBQUksQ0FBQyxJQUFJLENBQUMsQ0FBQyxDQUFDLEVBQ1osQ0FBQyxHQUFHLElBQUksQ0FBQyxhQUFhLENBQUMsQ0FBQyxDQUFDLENBQUMsSUFBSSxDQUFDLENBQUMsQ0FBQyxFQUNqQyxDQUFDLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxHQUFHLElBQUksQ0FBQyxhQUFhLENBQUMsQ0FBQyxDQUFDLENBQUMsSUFBSSxDQUFDLENBQUMsQ0FBQyxFQUM3QyxDQUFDLENBQUMsQ0FBQyxHQUFHLENBQUMsR0FBRyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsR0FBRyxJQUFJLENBQUMsYUFBYSxDQUFDLENBQUMsQ0FBQyxDQUFDLElBQUksQ0FBQyxDQUFDLENBQUMsQ0FDbEQsQ0FBQztRQUVGLE9BQU8sQ0FBQyxHQUFHLEVBQUUsSUFBSSxDQUFDLENBQUM7SUFDckIsQ0FBQztJQUVEOzs7OztPQUtHO0lBQ0ksS0FBSztRQUNWLFFBQVEsSUFBSSxDQUFDLE1BQU0sRUFBRSxDQUFDO1lBQ3BCLEtBQUssQ0FBQztnQkFDSixPQUFPLElBQUksQ0FBQyxlQUFlLEVBQUUsQ0FBQztZQUNoQyxLQUFLLENBQUM7Z0JBQ0osT0FBTyxJQUFJLENBQUMsbUJBQW1CLEVBQUUsQ0FBQztZQUNwQyxLQUFLLENBQUM7Z0JBQ0osT0FBTyxJQUFJLENBQUMsZUFBZSxFQUFFLENBQUM7WUFDaEMsS0FBSyxDQUFDO2dCQUNKLE9BQU8sRUFBRSxDQUFDO1lBQ1o7Z0JBQ0UsTUFBTSxJQUFJLEtBQUssQ0FBQyxrQ0FBa0MsSUFBSSxDQUFDLE1BQU0sRUFBRSxDQUFDLENBQUM7UUFDckUsQ0FBQztJQUNILENBQUM7SUFFRDs7T0FFRztJQUNJLFlBQVk7UUFDakIsT0FBTyxJQUFJLENBQUMsYUFBYSxFQUFFLENBQUMsS0FBSyxFQUFFLENBQUM7SUFDdEMsQ0FBQztJQUVEOztPQUVHO0lBQ0ksY0FBYztRQUNuQixNQUFNLEdBQUcsR0FBRyxJQUFJLENBQUMsWUFBWSxFQUFFLENBQUM7UUFDaEMsTUFBTSxNQUFNLEdBQUcsRUFBRSxDQUFDO1FBQ2xCLEtBQUssSUFBSSxDQUFDLEdBQUcsQ0FBQyxFQUFFLENBQUMsR0FBRyxHQUFHLENBQUMsTUFBTSxFQUFFLENBQUMsRUFBRSxFQUFFLENBQUM7WUFDcEMsTUFBTSxDQUFDLEdBQUcsR0FBRyxDQUFDLENBQUMsQ0FBQyxDQUFDO1lBQ2pCLElBQUksQ0FBQyxJQUFJLENBQUMsSUFBSSxDQUFDLElBQUksQ0FBQyxFQUFFLENBQUM7Z0JBQ3JCLE1BQU0sQ0FBQyxJQUFJLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUM7WUFDdEIsQ0FBQztRQUNILENBQUM7UUFDRCxPQUFPLE1BQU0sQ0FBQztJQUNoQixDQUFDO0lBRUQ7O09BRUc7SUFDSSxhQUFhO1FBQ2xCLElBQUksS0FBSyxHQUFHLENBQUMsSUFBSSxDQUFDLElBQUksQ0FBQyxDQUFDLENBQUMsRUFBRSxJQUFJLENBQUMsSUFBSSxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUM7UUFFekMseUVBQXlFO1FBQ3pFLFNBQVM7UUFDVCxNQUFNLFdBQVcsR0FBRyxDQUFDLEtBQWEsRUFBRSxFQUFFO1lBQ3BDLElBQUksS0FBSyxDQUFDLENBQUMsQ0FBQyxHQUFHLEtBQUssQ0FBQyxDQUFDLENBQUMsRUFBRSxDQUFDO2dCQUN4QixLQUFLLEdBQUcsQ0FBQyxJQUFJLENBQUMsR0FBRyxDQUFDLEtBQUssQ0FBQyxDQUFDLENBQUMsRUFBRSxLQUFLLENBQUMsRUFBRSxJQUFJLENBQUMsR0FBRyxDQUFDLEtBQUssQ0FBQyxDQUFDLENBQUMsRUFBRSxLQUFLLENBQUMsQ0FBQyxDQUFDO1lBQ2pFLENBQUM7aUJBQU0sQ0FBQztnQkFDTixLQUFLLEdBQUcsQ0FBQyxJQUFJLENBQUMsR0FBRyxDQUFDLEtBQUssQ0FBQyxDQUFDLENBQUMsRUFBRSxLQUFLLENBQUMsRUFBRSxJQUFJLENBQUMsR0FBRyxDQUFDLEtBQUssQ0FBQyxDQUFDLENBQUMsRUFBRSxLQUFLLENBQUMsQ0FBQyxDQUFDO1lBQ2pFLENBQUM7UUFDSCxDQUFDLENBQUM7UUFFRixJQUFJLENBQUMsY0FBYyxFQUFFLENBQUMsT0FBTyxDQUFDLENBQUMsQ0FBQyxFQUFFLENBQUMsV0FBVyxDQUFDLElBQUksQ0FBQyxJQUFJLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDO1FBRTlELE9BQU8sS0FBSyxDQUFDO0lBQ2YsQ0FBQztJQUVPLGVBQWU7UUFDckIsTUFBTSxDQUFDLEdBQUcsSUFBSSxDQUFDLEVBQUUsQ0FBQztRQUNsQixNQUFNLENBQUMsR0FBRyxJQUFJLENBQUMsRUFBRSxDQUFDO1FBQ2xCLE1BQU0sQ0FBQyxHQUFHLElBQUksQ0FBQyxFQUFFLENBQUM7UUFDbEIsTUFBTSxDQUFDLEdBQUcsSUFBSSxDQUFDLEVBQUUsQ0FBQztRQUVsQixzREFBc0Q7UUFDdEQsTUFBTSxFQUFFLEdBQUcsQ0FBQyxHQUFHLENBQUMsQ0FBQztRQUNqQixNQUFNLEVBQUUsR0FBRyxDQUFDLEdBQUcsQ0FBQyxDQUFDO1FBQ2pCLE1BQU0sRUFBRSxHQUFHLENBQUMsR0FBRyxDQUFDLENBQUM7UUFDakIsTUFBTSxDQUFDLEdBQUcsQ0FBQyxDQUFDLEdBQUcsRUFBRSxHQUFHLEVBQUUsQ0FBQyxHQUFHLENBQUMsQ0FBQyxHQUFHLEVBQUUsQ0FBQyxDQUFDO1FBQ25DLE1BQU0sQ0FBQyxHQUFHLENBQUMsQ0FBQyxHQUFHLEVBQUUsR0FBRyxDQUFDLEdBQUcsQ0FBQyxHQUFHLEVBQUUsR0FBRyxDQUFDLEdBQUcsRUFBRSxHQUFHLEVBQUUsR0FBRyxDQUFDLENBQUMsR0FBRyxDQUFDLEVBQUUsR0FBRyxFQUFFLEdBQUcsQ0FBQyxDQUFDLENBQUM7UUFFbEUsTUFBTSxHQUFHLEdBQUcsSUFBSSxDQUFDLHdCQUF3QixDQUFDLENBQUMsRUFBRSxDQUFDLENBQUMsQ0FBQztRQUVoRCx3RUFBd0U7UUFDeEUsZUFBZTtRQUNmLE1BQU0sYUFBYSxHQUFHLENBQUMsQ0FBUyxFQUFFLEVBQUUsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxHQUFHLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDO1FBQ3JELFFBQVEsR0FBRyxDQUFDLE1BQU0sRUFBRSxDQUFDO1lBQ25CLEtBQUssQ0FBQztnQkFDSixPQUFPLENBQUMsYUFBYSxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUM7WUFDakMsS0FBSyxDQUFDO2dCQUNKLE9BQU8sQ0FBQyxhQUFhLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQyxDQUFDLEVBQUUsYUFBYSxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUM7WUFDeEQsS0FBSyxDQUFDO2dCQUNKLE9BQU87b0JBQ0wsYUFBYSxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsQ0FBQztvQkFDckIsYUFBYSxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsQ0FBQztvQkFDckIsYUFBYSxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsQ0FBQztpQkFDdEIsQ0FBQztZQUNKO2dCQUNFLE9BQU8sRUFBRSxDQUFDO1FBQ2QsQ0FBQztJQUNILENBQUM7SUFFTyx3QkFBd0IsQ0FBQyxDQUFTLEVBQUUsQ0FBUztRQUNuRCxjQUFjO1FBRWQsMkRBQTJEO1FBQzNELElBQUksSUFBSSxDQUFDLFVBQVUsQ0FBQyxDQUFDLENBQUMsRUFBRSxDQUFDO1lBQ3ZCLE9BQU8sQ0FBQyxJQUFJLENBQUMsSUFBSSxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQztRQUN6QixDQUFDO1FBRUQsTUFBTSxHQUFHLEdBQUcsSUFBSSxDQUFDLEVBQUUsR0FBRyxDQUFDLENBQUM7UUFDeEIsTUFBTSxZQUFZLEdBQUcsQ0FBQyxHQUFHLENBQUMsR0FBRyxDQUFDLEdBQUcsQ0FBQyxHQUFHLEVBQUUsR0FBRyxDQUFDLEdBQUcsQ0FBQyxDQUFDO1FBQ2hELElBQUksWUFBWSxHQUFHLE9BQU8sRUFBRSxDQUFDO1lBQzNCLG1EQUFtRDtZQUNuRCxNQUFNLEdBQUcsR0FBRyxDQUFDLEdBQUcsSUFBSSxDQUFDLElBQUksQ0FBQyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQztZQUNsQyxNQUFNLFNBQVMsR0FBRyxDQUFDLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDLEdBQUcsSUFBSSxDQUFDLElBQUksQ0FBQyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQztZQUUxRCxNQUFNLE9BQU8sR0FBRyxDQUFDLENBQVMsRUFBRSxFQUFFLENBQzVCLEdBQUc7Z0JBQ0gsSUFBSSxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsR0FBRyxJQUFJLENBQUMsSUFBSSxDQUFDLEtBQUssQ0FBQyxDQUFDLENBQUMsRUFBRSxDQUFDLEVBQUUsU0FBUyxDQUFDLENBQUMsR0FBRyxDQUFDLEdBQUcsR0FBRyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQztZQUV6RSxpRUFBaUU7WUFDakUsdUVBQXVFO1lBQ3ZFLElBQUksU0FBUyxJQUFJLE1BQU0sRUFBRSxDQUFDO2dCQUN4Qiw2Q0FBNkM7Z0JBQzdDLE9BQU8sQ0FBQyxPQUFPLENBQUMsQ0FBQyxDQUFDLEVBQUUsT0FBTyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUM7WUFDbEMsQ0FBQztZQUVELElBQUksU0FBUyxJQUFJLENBQUMsTUFBTSxFQUFFLENBQUM7Z0JBQ3pCLDZDQUE2QztnQkFDN0MsT0FBTyxDQUFDLE9BQU8sQ0FBQyxDQUFDLENBQUMsRUFBRSxPQUFPLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQztZQUNsQyxDQUFDO1lBRUQsT0FBTyxDQUFDLE9BQU8sQ0FBQyxDQUFDLENBQUMsRUFBRSxPQUFPLENBQUMsQ0FBQyxDQUFDLEVBQUUsT0FBTyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUM7UUFDOUMsQ0FBQztRQUVELElBQUksWUFBWSxHQUFHLENBQUMsSUFBSSxDQUFDLEdBQUcsQ0FBQyxFQUFFLENBQUM7WUFDOUIsV0FBVztZQUNYLE1BQU0sU0FBUyxHQUNiLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQztnQkFDUCxJQUFJLENBQUMsS0FBSyxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUMsR0FBRyxJQUFJLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsR0FBRyxJQUFJLENBQUMsSUFBSSxDQUFDLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDLENBQUM7WUFDakUsTUFBTSxDQUFDLEdBQUcsQ0FBQyxDQUFDLEdBQUcsSUFBSSxDQUFDLElBQUksQ0FBQyxDQUFDLENBQUMsR0FBRyxJQUFJLENBQUMsSUFBSSxDQUFDLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxHQUFHLElBQUksQ0FBQyxJQUFJLENBQUMsU0FBUyxDQUFDLENBQUM7WUFDdkUsT0FBTyxDQUFDLENBQUMsQ0FBQyxDQUFDO1FBQ2IsQ0FBQztRQUVELElBQUksQ0FBQyxHQUFHLENBQUMsRUFBRSxDQUFDO1lBQ1YsV0FBVztZQUNYLE1BQU0sU0FBUyxHQUNiLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxHQUFHLElBQUksQ0FBQyxLQUFLLENBQUMsQ0FBQyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQyxHQUFHLElBQUksQ0FBQyxJQUFJLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDLENBQUM7WUFDL0QsTUFBTSxDQUFDLEdBQUcsQ0FBQyxDQUFDLEdBQUcsSUFBSSxDQUFDLElBQUksQ0FBQyxDQUFDLEdBQUcsQ0FBQyxDQUFDLEdBQUcsSUFBSSxDQUFDLElBQUksQ0FBQyxTQUFTLENBQUMsQ0FBQztZQUN2RCxPQUFPLENBQUMsQ0FBQyxDQUFDLENBQUM7UUFDYixDQUFDO1FBRUQsV0FBVztRQUNYLE9BQU8sRUFBRSxDQUFDO0lBQ1osQ0FBQztJQUVPLG1CQUFtQjtRQUN6QixNQUFNLENBQUMsR0FBRyxJQUFJLENBQUMsRUFBRSxDQUFDO1FBQ2xCLE1BQU0sQ0FBQyxHQUFHLElBQUksQ0FBQyxFQUFFLENBQUM7UUFDbEIsTUFBTSxDQUFDLEdBQUcsSUFBSSxDQUFDLEVBQUUsQ0FBQztRQUNsQixNQUFNLFdBQVcsR0FBRyxDQUFDLEdBQUcsQ0FBQyxHQUFHLENBQUMsR0FBRyxDQUFDLEdBQUcsQ0FBQyxDQUFDO1FBRXRDLElBQUksSUFBSSxDQUFDLFVBQVUsQ0FBQyxXQUFXLENBQUMsRUFBRSxDQUFDO1lBQ2pDLHdDQUF3QztZQUN4QyxPQUFPLENBQUMsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsQ0FBQztRQUN4QixDQUFDO1FBRUQsSUFBSSxXQUFXLElBQUksQ0FBQyxFQUFFLENBQUM7WUFDckIsTUFBTSxJQUFJLEdBQUcsSUFBSSxDQUFDLElBQUksQ0FBQyxXQUFXLENBQUMsQ0FBQztZQUNwQyx3QkFBd0I7WUFDeEIsTUFBTSxFQUFFLEdBQUcsQ0FBQyxDQUFDLENBQUMsR0FBRyxJQUFJLENBQUMsR0FBRyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQztZQUNqQyxNQUFNLEVBQUUsR0FBRyxDQUFDLENBQUMsQ0FBQyxHQUFHLElBQUksQ0FBQyxHQUFHLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDO1lBRWpDLE9BQU8sQ0FBQyxJQUFJLENBQUMsR0FBRyxDQUFDLEVBQUUsRUFBRSxFQUFFLENBQUMsRUFBRSxJQUFJLENBQUMsR0FBRyxDQUFDLEVBQUUsRUFBRSxFQUFFLENBQUMsQ0FBQyxDQUFDO1FBQzlDLENBQUM7UUFFRCxPQUFPLEVBQUUsQ0FBQztJQUNaLENBQUM7SUFFTyxlQUFlO1FBQ3JCLE9BQU8sQ0FBQyxDQUFDLElBQUksQ0FBQyxFQUFFLEdBQUcsSUFBSSxDQUFDLEVBQUUsQ0FBQyxDQUFDO0lBQzlCLENBQUM7SUFFTyxVQUFVLENBQUMsS0FBYTtRQUM5QixPQUFPLElBQUksQ0FBQyxHQUFHLENBQUMsQ0FBQyxHQUFHLEtBQUssQ0FBQyxJQUFJLE1BQU0sQ0FBQyxPQUFPLENBQUM7SUFDL0MsQ0FBQztDQUNGIn0=

package/lib/curves/Polynomial2D.d.ts

@@ -0,0 +1,22 @@
+import { BBox, Vector2 } from '@twick/core';
+import { Polynomial } from './Polynomial';
+export declare class Polynomial2D {
+    readonly c0: Vector2 | Polynomial;
+    readonly c1: Vector2 | Polynomial;
+    readonly c2?: Vector2 | undefined;
+    readonly c3?: Vector2 | undefined;
+    readonly x: Polynomial;
+    readonly y: Polynomial;
+    constructor(c0: Vector2, c1: Vector2, c2: Vector2, c3: Vector2);
+    constructor(c0: Vector2, c1: Vector2, c2: Vector2);
+    constructor(x: Polynomial, y: Polynomial);
+    eval(t: number, derivative?: number): Vector2;
+    split(u: number): [Polynomial2D, Polynomial2D];
+    differentiate(n?: number): Polynomial2D;
+    evalDerivative(t: number): Vector2;
+    /**
+     * Calculate the tight axis-aligned bounds of the curve in the unit interval.
+     */
+    getBounds(): BBox;
+}
+//# sourceMappingURL=Polynomial2D.d.ts.map

package/lib/curves/Polynomial2D.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"Polynomial2D.d.ts","sourceRoot":"","sources":["../../src/lib/curves/Polynomial2D.ts"],"names":[],"mappings":"AAAA,OAAO,EAAC,IAAI,EAAE,OAAO,EAAC,MAAM,aAAa,CAAC;AAE1C,OAAO,EAAC,UAAU,EAAC,MAAM,cAAc,CAAC;AAExC,qBAAa,YAAY;aAQL,EAAE,EAAE,OAAO,GAAG,UAAU;aACxB,EAAE,EAAE,OAAO,GAAG,UAAU;aACxB,EAAE,CAAC,EAAE,OAAO;aACZ,EAAE,CAAC,EAAE,OAAO;IAV9B,SAAgB,CAAC,EAAE,UAAU,CAAC;IAC9B,SAAgB,CAAC,EAAE,UAAU,CAAC;gBAEX,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO;gBAClD,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO;gBACrC,CAAC,EAAE,UAAU,EAAE,CAAC,EAAE,UAAU;IAmBxC,IAAI,CAAC,CAAC,EAAE,MAAM,EAAE,UAAU,SAAI,GAAG,OAAO;IAOxC,KAAK,CAAC,CAAC,EAAE,MAAM,GAAG,CAAC,YAAY,EAAE,YAAY,CAAC;IAM9C,aAAa,CAAC,CAAC,SAAI,GAAG,YAAY;IAIlC,cAAc,CAAC,CAAC,EAAE,MAAM,GAAG,OAAO;IAIzC;;OAEG;IACI,SAAS,IAAI,IAAI;CAQzB"}

package/lib/curves/Polynomial2D.js

@@ -0,0 +1,45 @@
+import { BBox, Vector2 } from '@twick/core';
+import { Polynomial } from './Polynomial';
+export class Polynomial2D {
+    constructor(c0, c1, c2, c3) {
+        this.c0 = c0;
+        this.c1 = c1;
+        this.c2 = c2;
+        this.c3 = c3;
+        if (c0 instanceof Polynomial) {
+            this.x = c0;
+            this.y = c1;
+        }
+        else if (c3 !== undefined) {
+            this.x = new Polynomial(c0.x, c1.x, c2.x, c3.x);
+            this.y = new Polynomial(c0.y, c1.y, c2.y, c3.y);
+        }
+        else {
+            this.x = new Polynomial(c0.x, c1.x, c2.x);
+            this.y = new Polynomial(c0.y, c1.y, c2.y);
+        }
+    }
+    eval(t, derivative = 0) {
+        return new Vector2(this.x.differentiate(derivative).eval(t), this.y.differentiate(derivative).eval(t));
+    }
+    split(u) {
+        const [xPre, xPost] = this.x.split(u);
+        const [yPre, yPost] = this.y.split(u);
+        return [new Polynomial2D(xPre, yPre), new Polynomial2D(xPost, yPost)];
+    }
+    differentiate(n = 1) {
+        return new Polynomial2D(this.x.differentiate(n), this.y.differentiate(n));
+    }
+    evalDerivative(t) {
+        return this.differentiate().eval(t);
+    }
+    /**
+     * Calculate the tight axis-aligned bounds of the curve in the unit interval.
+     */
+    getBounds() {
+        const rangeX = this.x.outputRange01();
+        const rangeY = this.y.outputRange01();
+        return BBox.fromPoints(new Vector2(Math.min(...rangeX), Math.max(...rangeY)), new Vector2(Math.max(...rangeX), Math.min(...rangeY)));
+    }
+}
+//# sourceMappingURL=data:application/json;base64,eyJ2ZXJzaW9uIjozLCJmaWxlIjoiUG9seW5vbWlhbDJELmpzIiwic291cmNlUm9vdCI6IiIsInNvdXJjZXMiOlsiLi4vLi4vc3JjL2xpYi9jdXJ2ZXMvUG9seW5vbWlhbDJELnRzIl0sIm5hbWVzIjpbXSwibWFwcGluZ3MiOiJBQUFBLE9BQU8sRUFBQyxJQUFJLEVBQUUsT0FBTyxFQUFDLE1BQU0sYUFBYSxDQUFDO0FBRTFDLE9BQU8sRUFBQyxVQUFVLEVBQUMsTUFBTSxjQUFjLENBQUM7QUFFeEMsTUFBTSxPQUFPLFlBQVk7SUFPdkIsWUFDa0IsRUFBd0IsRUFDeEIsRUFBd0IsRUFDeEIsRUFBWSxFQUNaLEVBQVk7UUFIWixPQUFFLEdBQUYsRUFBRSxDQUFzQjtRQUN4QixPQUFFLEdBQUYsRUFBRSxDQUFzQjtRQUN4QixPQUFFLEdBQUYsRUFBRSxDQUFVO1FBQ1osT0FBRSxHQUFGLEVBQUUsQ0FBVTtRQUU1QixJQUFJLEVBQUUsWUFBWSxVQUFVLEVBQUUsQ0FBQztZQUM3QixJQUFJLENBQUMsQ0FBQyxHQUFHLEVBQUUsQ0FBQztZQUNaLElBQUksQ0FBQyxDQUFDLEdBQUcsRUFBZ0IsQ0FBQztRQUM1QixDQUFDO2FBQU0sSUFBSSxFQUFFLEtBQUssU0FBUyxFQUFFLENBQUM7WUFDNUIsSUFBSSxDQUFDLENBQUMsR0FBRyxJQUFJLFVBQVUsQ0FBQyxFQUFFLENBQUMsQ0FBQyxFQUFHLEVBQWMsQ0FBQyxDQUFDLEVBQUUsRUFBRyxDQUFDLENBQUMsRUFBRSxFQUFFLENBQUMsQ0FBQyxDQUFDLENBQUM7WUFDOUQsSUFBSSxDQUFDLENBQUMsR0FBRyxJQUFJLFVBQVUsQ0FBQyxFQUFFLENBQUMsQ0FBQyxFQUFHLEVBQWMsQ0FBQyxDQUFDLEVBQUUsRUFBRyxDQUFDLENBQUMsRUFBRSxFQUFFLENBQUMsQ0FBQyxDQUFDLENBQUM7UUFDaEUsQ0FBQzthQUFNLENBQUM7WUFDTixJQUFJLENBQUMsQ0FBQyxHQUFHLElBQUksVUFBVSxDQUFDLEVBQUUsQ0FBQyxDQUFDLEVBQUcsRUFBYyxDQUFDLENBQUMsRUFBRSxFQUFHLENBQUMsQ0FBQyxDQUFDLENBQUM7WUFDeEQsSUFBSSxDQUFDLENBQUMsR0FBRyxJQUFJLFVBQVUsQ0FBQyxFQUFFLENBQUMsQ0FBQyxFQUFHLEVBQWMsQ0FBQyxDQUFDLEVBQUUsRUFBRyxDQUFDLENBQUMsQ0FBQyxDQUFDO1FBQzFELENBQUM7SUFDSCxDQUFDO0lBRU0sSUFBSSxDQUFDLENBQVMsRUFBRSxVQUFVLEdBQUcsQ0FBQztRQUNuQyxPQUFPLElBQUksT0FBTyxDQUNoQixJQUFJLENBQUMsQ0FBQyxDQUFDLGFBQWEsQ0FBQyxVQUFVLENBQUMsQ0FBQyxJQUFJLENBQUMsQ0FBQyxDQUFDLEVBQ3hDLElBQUksQ0FBQyxDQUFDLENBQUMsYUFBYSxDQUFDLFVBQVUsQ0FBQyxDQUFDLElBQUksQ0FBQyxDQUFDLENBQUMsQ0FDekMsQ0FBQztJQUNKLENBQUM7SUFFTSxLQUFLLENBQUMsQ0FBUztRQUNwQixNQUFNLENBQUMsSUFBSSxFQUFFLEtBQUssQ0FBQyxHQUFHLElBQUksQ0FBQyxDQUFDLENBQUMsS0FBSyxDQUFDLENBQUMsQ0FBQyxDQUFDO1FBQ3RDLE1BQU0sQ0FBQyxJQUFJLEVBQUUsS0FBSyxDQUFDLEdBQUcsSUFBSSxDQUFDLENBQUMsQ0FBQyxLQUFLLENBQUMsQ0FBQyxDQUFDLENBQUM7UUFDdEMsT0FBTyxDQUFDLElBQUksWUFBWSxDQUFDLElBQUksRUFBRSxJQUFJLENBQUMsRUFBRSxJQUFJLFlBQVksQ0FBQyxLQUFLLEVBQUUsS0FBSyxDQUFDLENBQUMsQ0FBQztJQUN4RSxDQUFDO0lBRU0sYUFBYSxDQUFDLENBQUMsR0FBRyxDQUFDO1FBQ3hCLE9BQU8sSUFBSSxZQUFZLENBQUMsSUFBSSxDQUFDLENBQUMsQ0FBQyxhQUFhLENBQUMsQ0FBQyxDQUFDLEVBQUUsSUFBSSxDQUFDLENBQUMsQ0FBQyxhQUFhLENBQUMsQ0FBQyxDQUFDLENBQUMsQ0FBQztJQUM1RSxDQUFDO0lBRU0sY0FBYyxDQUFDLENBQVM7UUFDN0IsT0FBTyxJQUFJLENBQUMsYUFBYSxFQUFFLENBQUMsSUFBSSxDQUFDLENBQUMsQ0FBQyxDQUFDO0lBQ3RDLENBQUM7SUFFRDs7T0FFRztJQUNJLFNBQVM7UUFDZCxNQUFNLE1BQU0sR0FBRyxJQUFJLENBQUMsQ0FBQyxDQUFDLGFBQWEsRUFBRSxDQUFDO1FBQ3RDLE1BQU0sTUFBTSxHQUFHLElBQUksQ0FBQyxDQUFDLENBQUMsYUFBYSxFQUFFLENBQUM7UUFDdEMsT0FBTyxJQUFJLENBQUMsVUFBVSxDQUNwQixJQUFJLE9BQU8sQ0FBQyxJQUFJLENBQUMsR0FBRyxDQUFDLEdBQUcsTUFBTSxDQUFDLEVBQUUsSUFBSSxDQUFDLEdBQUcsQ0FBQyxHQUFHLE1BQU0sQ0FBQyxDQUFDLEVBQ3JELElBQUksT0FBTyxDQUFDLElBQUksQ0FBQyxHQUFHLENBQUMsR0FBRyxNQUFNLENBQUMsRUFBRSxJQUFJLENBQUMsR0FBRyxDQUFDLEdBQUcsTUFBTSxDQUFDLENBQUMsQ0FDdEQsQ0FBQztJQUNKLENBQUM7Q0FDRiJ9

package/lib/curves/PolynomialSegment.d.ts

@@ -0,0 +1,39 @@
+import { type BBox, type Vector2 } from '@twick/core';
+import type { CurvePoint } from './CurvePoint';
+import type { Polynomial2D } from './Polynomial2D';
+import { Segment } from './Segment';
+import { UniformPolynomialCurveSampler } from './UniformPolynomialCurveSampler';
+export declare abstract class PolynomialSegment extends Segment {
+    protected readonly curve: Polynomial2D;
+    protected readonly length: number;
+    protected readonly pointSampler: UniformPolynomialCurveSampler;
+    get arcLength(): number;
+    abstract get points(): Vector2[];
+    protected constructor(curve: Polynomial2D, length: number);
+    getBBox(): BBox;
+    /**
+     * Evaluate the polynomial at the given t value.
+     *
+     * @param t - The t value at which to evaluate the curve.
+     */
+    eval(t: number): CurvePoint;
+    /**
+     * Split the curve into two separate polynomials at the given t value. The two
+     * resulting curves form the same overall shape as the original curve.
+     *
+     * @param t - The t value at which to split the curve.
+     */
+    abstract split(t: number): [PolynomialSegment, PolynomialSegment];
+    getPoint(distance: number): CurvePoint;
+    transformPoints(matrix: DOMMatrix): Vector2[];
+    /**
+     * Return the tangent of the point that sits at the provided t value on the
+     * curve.
+     *
+     * @param t - The t value at which to evaluate the curve.
+     */
+    tangent(t: number): Vector2;
+    draw(context: CanvasRenderingContext2D | Path2D, start?: number, end?: number, move?: boolean): [CurvePoint, CurvePoint];
+    protected abstract doDraw(context: CanvasRenderingContext2D | Path2D): void;
+}
+//# sourceMappingURL=PolynomialSegment.d.ts.map

package/lib/curves/PolynomialSegment.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"PolynomialSegment.d.ts","sourceRoot":"","sources":["../../src/lib/curves/PolynomialSegment.ts"],"names":[],"mappings":"AAAA,OAAO,EAAyB,KAAK,IAAI,EAAE,KAAK,OAAO,EAAC,MAAM,aAAa,CAAC;AAE5E,OAAO,KAAK,EAAC,UAAU,EAAC,MAAM,cAAc,CAAC;AAC7C,OAAO,KAAK,EAAC,YAAY,EAAC,MAAM,gBAAgB,CAAC;AACjD,OAAO,EAAC,OAAO,EAAC,MAAM,WAAW,CAAC;AAClC,OAAO,EAAC,6BAA6B,EAAC,MAAM,iCAAiC,CAAC;AAE9E,8BAAsB,iBAAkB,SAAQ,OAAO;IAUnD,SAAS,CAAC,QAAQ,CAAC,KAAK,EAAE,YAAY;IACtC,SAAS,CAAC,QAAQ,CAAC,MAAM,EAAE,MAAM;IAVnC,SAAS,CAAC,QAAQ,CAAC,YAAY,EAAE,6BAA6B,CAAC;IAE/D,IAAW,SAAS,IAAI,MAAM,CAE7B;IAED,aAA6B,MAAM,IAAI,OAAO,EAAE,CAAC;IAEjD,SAAS,aACY,KAAK,EAAE,YAAY,EACnB,MAAM,EAAE,MAAM;IAM5B,OAAO,IAAI,IAAI;IAItB;;;;OAIG;IACI,IAAI,CAAC,CAAC,EAAE,MAAM,GAAG,UAAU;IAUlC;;;;;OAKG;aACa,KAAK,CAAC,CAAC,EAAE,MAAM,GAAG,CAAC,iBAAiB,EAAE,iBAAiB,CAAC;IAEjE,QAAQ,CAAC,QAAQ,EAAE,MAAM,GAAG,UAAU;IAWtC,eAAe,CAAC,MAAM,EAAE,SAAS,GAAG,OAAO,EAAE;IAIpD;;;;;OAKG;IACI,OAAO,CAAC,CAAC,EAAE,MAAM,GAAG,OAAO;IAI3B,IAAI,CACT,OAAO,EAAE,wBAAwB,GAAG,MAAM,EAC1C,KAAK,SAAI,EACT,GAAG,SAAI,EACP,IAAI,UAAO,GACV,CAAC,UAAU,EAAE,UAAU,CAAC;IAyC3B,SAAS,CAAC,QAAQ,CAAC,MAAM,CAAC,OAAO,EAAE,wBAAwB,GAAG,MAAM,GAAG,IAAI;CAC5E"}

package/lib/curves/PolynomialSegment.js

@@ -0,0 +1,86 @@
+import { transformVectorAsPoint } from '@twick/core';
+import { moveTo } from '../utils';
+import { Segment } from './Segment';
+import { UniformPolynomialCurveSampler } from './UniformPolynomialCurveSampler';
+export class PolynomialSegment extends Segment {
+    get arcLength() {
+        return this.length;
+    }
+    constructor(curve, length) {
+        super();
+        this.curve = curve;
+        this.length = length;
+        this.pointSampler = new UniformPolynomialCurveSampler(this);
+    }
+    getBBox() {
+        return this.curve.getBounds();
+    }
+    /**
+     * Evaluate the polynomial at the given t value.
+     *
+     * @param t - The t value at which to evaluate the curve.
+     */
+    eval(t) {
+        const tangent = this.tangent(t);
+        return {
+            position: this.curve.eval(t),
+            tangent,
+            normal: tangent.perpendicular,
+        };
+    }
+    getPoint(distance) {
+        const closestPoint = this.pointSampler.pointAtDistance(this.arcLength * distance);
+        return {
+            position: closestPoint.position,
+            tangent: closestPoint.tangent,
+            normal: closestPoint.tangent.perpendicular,
+        };
+    }
+    transformPoints(matrix) {
+        return this.points.map(point => transformVectorAsPoint(point, matrix));
+    }
+    /**
+     * Return the tangent of the point that sits at the provided t value on the
+     * curve.
+     *
+     * @param t - The t value at which to evaluate the curve.
+     */
+    tangent(t) {
+        return this.curve.evalDerivative(t).normalized;
+    }
+    draw(context, start = 0, end = 1, move = true) {
+        let curve = null;
+        let startT = start;
+        let endT = end;
+        let points = this.points;
+        if (start !== 0 || end !== 1) {
+            const startDistance = this.length * start;
+            const endDistance = this.length * end;
+            startT = this.pointSampler.distanceToT(startDistance);
+            endT = this.pointSampler.distanceToT(endDistance);
+            const relativeEndT = (endT - startT) / (1 - startT);
+            const [, startSegment] = this.split(startT);
+            [curve] = startSegment.split(relativeEndT);
+            points = curve.points;
+        }
+        if (move) {
+            moveTo(context, points[0]);
+        }
+        (curve ?? this).doDraw(context);
+        const startTangent = this.tangent(startT);
+        const endTangent = this.tangent(endT);
+        return [
+            {
+                position: points[0],
+                tangent: startTangent,
+                normal: startTangent.perpendicular,
+            },
+            {
+                position: points.at(-1),
+                tangent: endTangent,
+                normal: endTangent.perpendicular,
+            },
+        ];
+    }
+}
+//# sourceMappingURL=data:application/json;base64,eyJ2ZXJzaW9uIjozLCJmaWxlIjoiUG9seW5vbWlhbFNlZ21lbnQuanMiLCJzb3VyY2VSb290IjoiIiwic291cmNlcyI6WyIuLi8uLi9zcmMvbGliL2N1cnZlcy9Qb2x5bm9taWFsU2VnbWVudC50cyJdLCJuYW1lcyI6W10sIm1hcHBpbmdzIjoiQUFBQSxPQUFPLEVBQUMsc0JBQXNCLEVBQTBCLE1BQU0sYUFBYSxDQUFDO0FBQzVFLE9BQU8sRUFBQyxNQUFNLEVBQUMsTUFBTSxVQUFVLENBQUM7QUFHaEMsT0FBTyxFQUFDLE9BQU8sRUFBQyxNQUFNLFdBQVcsQ0FBQztBQUNsQyxPQUFPLEVBQUMsNkJBQTZCLEVBQUMsTUFBTSxpQ0FBaUMsQ0FBQztBQUU5RSxNQUFNLE9BQWdCLGlCQUFrQixTQUFRLE9BQU87SUFHckQsSUFBVyxTQUFTO1FBQ2xCLE9BQU8sSUFBSSxDQUFDLE1BQU0sQ0FBQztJQUNyQixDQUFDO0lBSUQsWUFDcUIsS0FBbUIsRUFDbkIsTUFBYztRQUVqQyxLQUFLLEVBQUUsQ0FBQztRQUhXLFVBQUssR0FBTCxLQUFLLENBQWM7UUFDbkIsV0FBTSxHQUFOLE1BQU0sQ0FBUTtRQUdqQyxJQUFJLENBQUMsWUFBWSxHQUFHLElBQUksNkJBQTZCLENBQUMsSUFBSSxDQUFDLENBQUM7SUFDOUQsQ0FBQztJQUVNLE9BQU87UUFDWixPQUFPLElBQUksQ0FBQyxLQUFLLENBQUMsU0FBUyxFQUFFLENBQUM7SUFDaEMsQ0FBQztJQUVEOzs7O09BSUc7SUFDSSxJQUFJLENBQUMsQ0FBUztRQUNuQixNQUFNLE9BQU8sR0FBRyxJQUFJLENBQUMsT0FBTyxDQUFDLENBQUMsQ0FBQyxDQUFDO1FBRWhDLE9BQU87WUFDTCxRQUFRLEVBQUUsSUFBSSxDQUFDLEtBQUssQ0FBQyxJQUFJLENBQUMsQ0FBQyxDQUFDO1lBQzVCLE9BQU87WUFDUCxNQUFNLEVBQUUsT0FBTyxDQUFDLGFBQWE7U0FDOUIsQ0FBQztJQUNKLENBQUM7SUFVTSxRQUFRLENBQUMsUUFBZ0I7UUFDOUIsTUFBTSxZQUFZLEdBQUcsSUFBSSxDQUFDLFlBQVksQ0FBQyxlQUFlLENBQ3BELElBQUksQ0FBQyxTQUFTLEdBQUcsUUFBUSxDQUMxQixDQUFDO1FBQ0YsT0FBTztZQUNMLFFBQVEsRUFBRSxZQUFZLENBQUMsUUFBUTtZQUMvQixPQUFPLEVBQUUsWUFBWSxDQUFDLE9BQU87WUFDN0IsTUFBTSxFQUFFLFlBQVksQ0FBQyxPQUFPLENBQUMsYUFBYTtTQUMzQyxDQUFDO0lBQ0osQ0FBQztJQUVNLGVBQWUsQ0FBQyxNQUFpQjtRQUN0QyxPQUFPLElBQUksQ0FBQyxNQUFNLENBQUMsR0FBRyxDQUFDLEtBQUssQ0FBQyxFQUFFLENBQUMsc0JBQXNCLENBQUMsS0FBSyxFQUFFLE1BQU0sQ0FBQyxDQUFDLENBQUM7SUFDekUsQ0FBQztJQUVEOzs7OztPQUtHO0lBQ0ksT0FBTyxDQUFDLENBQVM7UUFDdEIsT0FBTyxJQUFJLENBQUMsS0FBSyxDQUFDLGNBQWMsQ0FBQyxDQUFDLENBQUMsQ0FBQyxVQUFVLENBQUM7SUFDakQsQ0FBQztJQUVNLElBQUksQ0FDVCxPQUEwQyxFQUMxQyxLQUFLLEdBQUcsQ0FBQyxFQUNULEdBQUcsR0FBRyxDQUFDLEVBQ1AsSUFBSSxHQUFHLElBQUk7UUFFWCxJQUFJLEtBQUssR0FBNkIsSUFBSSxDQUFDO1FBQzNDLElBQUksTUFBTSxHQUFHLEtBQUssQ0FBQztRQUNuQixJQUFJLElBQUksR0FBRyxHQUFHLENBQUM7UUFDZixJQUFJLE1BQU0sR0FBRyxJQUFJLENBQUMsTUFBTSxDQUFDO1FBRXpCLElBQUksS0FBSyxLQUFLLENBQUMsSUFBSSxHQUFHLEtBQUssQ0FBQyxFQUFFLENBQUM7WUFDN0IsTUFBTSxhQUFhLEdBQUcsSUFBSSxDQUFDLE1BQU0sR0FBRyxLQUFLLENBQUM7WUFDMUMsTUFBTSxXQUFXLEdBQUcsSUFBSSxDQUFDLE1BQU0sR0FBRyxHQUFHLENBQUM7WUFFdEMsTUFBTSxHQUFHLElBQUksQ0FBQyxZQUFZLENBQUMsV0FBVyxDQUFDLGFBQWEsQ0FBQyxDQUFDO1lBQ3RELElBQUksR0FBRyxJQUFJLENBQUMsWUFBWSxDQUFDLFdBQVcsQ0FBQyxXQUFXLENBQUMsQ0FBQztZQUNsRCxNQUFNLFlBQVksR0FBRyxDQUFDLElBQUksR0FBRyxNQUFNLENBQUMsR0FBRyxDQUFDLENBQUMsR0FBRyxNQUFNLENBQUMsQ0FBQztZQUVwRCxNQUFNLENBQUMsRUFBRSxZQUFZLENBQUMsR0FBRyxJQUFJLENBQUMsS0FBSyxDQUFDLE1BQU0sQ0FBQyxDQUFDO1lBQzVDLENBQUMsS0FBSyxDQUFDLEdBQUcsWUFBWSxDQUFDLEtBQUssQ0FBQyxZQUFZLENBQUMsQ0FBQztZQUMzQyxNQUFNLEdBQUcsS0FBSyxDQUFDLE1BQU0sQ0FBQztRQUN4QixDQUFDO1FBRUQsSUFBSSxJQUFJLEVBQUUsQ0FBQztZQUNULE1BQU0sQ0FBQyxPQUFPLEVBQUUsTUFBTSxDQUFDLENBQUMsQ0FBQyxDQUFDLENBQUM7UUFDN0IsQ0FBQztRQUNELENBQUMsS0FBSyxJQUFJLElBQUksQ0FBQyxDQUFDLE1BQU0sQ0FBQyxPQUFPLENBQUMsQ0FBQztRQUVoQyxNQUFNLFlBQVksR0FBRyxJQUFJLENBQUMsT0FBTyxDQUFDLE1BQU0sQ0FBQyxDQUFDO1FBQzFDLE1BQU0sVUFBVSxHQUFHLElBQUksQ0FBQyxPQUFPLENBQUMsSUFBSSxDQUFDLENBQUM7UUFFdEMsT0FBTztZQUNMO2dCQUNFLFFBQVEsRUFBRSxNQUFNLENBQUMsQ0FBQyxDQUFDO2dCQUNuQixPQUFPLEVBQUUsWUFBWTtnQkFDckIsTUFBTSxFQUFFLFlBQVksQ0FBQyxhQUFhO2FBQ25DO1lBQ0Q7Z0JBQ0UsUUFBUSxFQUFFLE1BQU0sQ0FBQyxFQUFFLENBQUMsQ0FBQyxDQUFDLENBQUU7Z0JBQ3hCLE9BQU8sRUFBRSxVQUFVO2dCQUNuQixNQUFNLEVBQUUsVUFBVSxDQUFDLGFBQWE7YUFDakM7U0FDRixDQUFDO0lBQ0osQ0FBQztDQUdGIn0=

package/lib/curves/QuadBezierSegment.d.ts

@@ -0,0 +1,17 @@
+import { Vector2 } from '@twick/core';
+import { PolynomialSegment } from './PolynomialSegment';
+/**
+ * A spline segment representing a quadratic Bézier curve.
+ */
+export declare class QuadBezierSegment extends PolynomialSegment {
+    readonly p0: Vector2;
+    readonly p1: Vector2;
+    readonly p2: Vector2;
+    private static el;
+    get points(): Vector2[];
+    constructor(p0: Vector2, p1: Vector2, p2: Vector2);
+    split(t: number): [PolynomialSegment, PolynomialSegment];
+    protected static getLength(p0: Vector2, p1: Vector2, p2: Vector2): number;
+    protected doDraw(context: CanvasRenderingContext2D | Path2D): void;
+}
+//# sourceMappingURL=QuadBezierSegment.d.ts.map

package/lib/curves/QuadBezierSegment.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"QuadBezierSegment.d.ts","sourceRoot":"","sources":["../../src/lib/curves/QuadBezierSegment.ts"],"names":[],"mappings":"AAAA,OAAO,EAAC,OAAO,EAAO,MAAM,aAAa,CAAC;AAG1C,OAAO,EAAC,iBAAiB,EAAC,MAAM,qBAAqB,CAAC;AAEtD;;GAEG;AACH,qBAAa,iBAAkB,SAAQ,iBAAiB;aASpC,EAAE,EAAE,OAAO;aACX,EAAE,EAAE,OAAO;aACX,EAAE,EAAE,OAAO;IAT7B,OAAO,CAAC,MAAM,CAAC,EAAE,CAAiB;IAElC,IAAW,MAAM,IAAI,OAAO,EAAE,CAE7B;gBAGiB,EAAE,EAAE,OAAO,EACX,EAAE,EAAE,OAAO,EACX,EAAE,EAAE,OAAO;IActB,KAAK,CAAC,CAAC,EAAE,MAAM,GAAG,CAAC,iBAAiB,EAAE,iBAAiB,CAAC;IAiB/D,SAAS,CAAC,MAAM,CAAC,SAAS,CAAC,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,GAAG,MAAM;cAUtD,MAAM,CAAC,OAAO,EAAE,wBAAwB,GAAG,MAAM;CAGrE"}

package/lib/curves/QuadBezierSegment.js

@@ -0,0 +1,49 @@
+var __decorate = (this && this.__decorate) || function (decorators, target, key, desc) {
+    var c = arguments.length, r = c < 3 ? target : desc === null ? desc = Object.getOwnPropertyDescriptor(target, key) : desc, d;
+    if (typeof Reflect === "object" && typeof Reflect.decorate === "function") r = Reflect.decorate(decorators, target, key, desc);
+    else for (var i = decorators.length - 1; i >= 0; i--) if (d = decorators[i]) r = (c < 3 ? d(r) : c > 3 ? d(target, key, r) : d(target, key)) || r;
+    return c > 3 && r && Object.defineProperty(target, key, r), r;
+};
+import { Vector2, lazy } from '@twick/core';
+import { quadraticCurveTo } from '../utils';
+import { Polynomial2D } from './Polynomial2D';
+import { PolynomialSegment } from './PolynomialSegment';
+/**
+ * A spline segment representing a quadratic Bézier curve.
+ */
+export class QuadBezierSegment extends PolynomialSegment {
+    get points() {
+        return [this.p0, this.p1, this.p2];
+    }
+    constructor(p0, p1, p2) {
+        super(new Polynomial2D(p0, 
+        // 2*(-p0+p1)
+        p0.flipped.add(p1).scale(2), 
+        // p0-2*p1+p2
+        p0.sub(p1.scale(2)).add(p2)), QuadBezierSegment.getLength(p0, p1, p2));
+        this.p0 = p0;
+        this.p1 = p1;
+        this.p2 = p2;
+    }
+    split(t) {
+        const a = new Vector2(this.p0.x + (this.p1.x - this.p0.x) * t, this.p0.y + (this.p1.y - this.p0.y) * t);
+        const b = new Vector2(this.p1.x + (this.p2.x - this.p1.x) * t, this.p1.y + (this.p2.y - this.p1.y) * t);
+        const p = new Vector2(a.x + (b.x - a.x) * t, a.y + (b.y - a.y) * t);
+        const left = new QuadBezierSegment(this.p0, a, p);
+        const right = new QuadBezierSegment(p, b, this.p2);
+        return [left, right];
+    }
+    static getLength(p0, p1, p2) {
+        // Let the browser do the work for us instead of calculating the arclength
+        // manually.
+        QuadBezierSegment.el.setAttribute('d', `M ${p0.x} ${p0.y} Q ${p1.x} ${p1.y} ${p2.x} ${p2.y}`);
+        return QuadBezierSegment.el.getTotalLength();
+    }
+    doDraw(context) {
+        quadraticCurveTo(context, this.p1, this.p2);
+    }
+}
+__decorate([
+    lazy(() => document.createElementNS('http://www.w3.org/2000/svg', 'path'))
+], QuadBezierSegment, "el", void 0);
+//# sourceMappingURL=data:application/json;base64,eyJ2ZXJzaW9uIjozLCJmaWxlIjoiUXVhZEJlemllclNlZ21lbnQuanMiLCJzb3VyY2VSb290IjoiIiwic291cmNlcyI6WyIuLi8uLi9zcmMvbGliL2N1cnZlcy9RdWFkQmV6aWVyU2VnbWVudC50cyJdLCJuYW1lcyI6W10sIm1hcHBpbmdzIjoiOzs7Ozs7QUFBQSxPQUFPLEVBQUMsT0FBTyxFQUFFLElBQUksRUFBQyxNQUFNLGFBQWEsQ0FBQztBQUMxQyxPQUFPLEVBQUMsZ0JBQWdCLEVBQUMsTUFBTSxVQUFVLENBQUM7QUFDMUMsT0FBTyxFQUFDLFlBQVksRUFBQyxNQUFNLGdCQUFnQixDQUFDO0FBQzVDLE9BQU8sRUFBQyxpQkFBaUIsRUFBQyxNQUFNLHFCQUFxQixDQUFDO0FBRXREOztHQUVHO0FBQ0gsTUFBTSxPQUFPLGlCQUFrQixTQUFRLGlCQUFpQjtJQUl0RCxJQUFXLE1BQU07UUFDZixPQUFPLENBQUMsSUFBSSxDQUFDLEVBQUUsRUFBRSxJQUFJLENBQUMsRUFBRSxFQUFFLElBQUksQ0FBQyxFQUFFLENBQUMsQ0FBQztJQUNyQyxDQUFDO0lBRUQsWUFDa0IsRUFBVyxFQUNYLEVBQVcsRUFDWCxFQUFXO1FBRTNCLEtBQUssQ0FDSCxJQUFJLFlBQVksQ0FDZCxFQUFFO1FBQ0YsYUFBYTtRQUNiLEVBQUUsQ0FBQyxPQUFPLENBQUMsR0FBRyxDQUFDLEVBQUUsQ0FBQyxDQUFDLEtBQUssQ0FBQyxDQUFDLENBQUM7UUFDM0IsYUFBYTtRQUNiLEVBQUUsQ0FBQyxHQUFHLENBQUMsRUFBRSxDQUFDLEtBQUssQ0FBQyxDQUFDLENBQUMsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxFQUFFLENBQUMsQ0FDNUIsRUFDRCxpQkFBaUIsQ0FBQyxTQUFTLENBQUMsRUFBRSxFQUFFLEVBQUUsRUFBRSxFQUFFLENBQUMsQ0FDeEMsQ0FBQztRQWJjLE9BQUUsR0FBRixFQUFFLENBQVM7UUFDWCxPQUFFLEdBQUYsRUFBRSxDQUFTO1FBQ1gsT0FBRSxHQUFGLEVBQUUsQ0FBUztJQVk3QixDQUFDO0lBRU0sS0FBSyxDQUFDLENBQVM7UUFDcEIsTUFBTSxDQUFDLEdBQUcsSUFBSSxPQUFPLENBQ25CLElBQUksQ0FBQyxFQUFFLENBQUMsQ0FBQyxHQUFHLENBQUMsSUFBSSxDQUFDLEVBQUUsQ0FBQyxDQUFDLEdBQUcsSUFBSSxDQUFDLEVBQUUsQ0FBQyxDQUFDLENBQUMsR0FBRyxDQUFDLEVBQ3ZDLElBQUksQ0FBQyxFQUFFLENBQUMsQ0FBQyxHQUFHLENBQUMsSUFBSSxDQUFDLEVBQUUsQ0FBQyxDQUFDLEdBQUcsSUFBSSxDQUFDLEVBQUUsQ0FBQyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQ3hDLENBQUM7UUFDRixNQUFNLENBQUMsR0FBRyxJQUFJLE9BQU8sQ0FDbkIsSUFBSSxDQUFDLEVBQUUsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxJQUFJLENBQUMsRUFBRSxDQUFDLENBQUMsR0FBRyxJQUFJLENBQUMsRUFBRSxDQUFDLENBQUMsQ0FBQyxHQUFHLENBQUMsRUFDdkMsSUFBSSxDQUFDLEVBQUUsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxJQUFJLENBQUMsRUFBRSxDQUFDLENBQUMsR0FBRyxJQUFJLENBQUMsRUFBRSxDQUFDLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FDeEMsQ0FBQztRQUNGLE1BQU0sQ0FBQyxHQUFHLElBQUksT0FBTyxDQUFDLENBQUMsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDLENBQUMsR0FBRyxDQUFDLEVBQUUsQ0FBQyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQyxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsQ0FBQyxHQUFHLENBQUMsQ0FBQyxDQUFDO1FBRXBFLE1BQU0sSUFBSSxHQUFHLElBQUksaUJBQWlCLENBQUMsSUFBSSxDQUFDLEVBQUUsRUFBRSxDQUFDLEVBQUUsQ0FBQyxDQUFDLENBQUM7UUFDbEQsTUFBTSxLQUFLLEdBQUcsSUFBSSxpQkFBaUIsQ0FBQyxDQUFDLEVBQUUsQ0FBQyxFQUFFLElBQUksQ0FBQyxFQUFFLENBQUMsQ0FBQztRQUVuRCxPQUFPLENBQUMsSUFBSSxFQUFFLEtBQUssQ0FBQyxDQUFDO0lBQ3ZCLENBQUM7SUFFUyxNQUFNLENBQUMsU0FBUyxDQUFDLEVBQVcsRUFBRSxFQUFXLEVBQUUsRUFBVztRQUM5RCwwRUFBMEU7UUFDMUUsWUFBWTtRQUNaLGlCQUFpQixDQUFDLEVBQUUsQ0FBQyxZQUFZLENBQy9CLEdBQUcsRUFDSCxLQUFLLEVBQUUsQ0FBQyxDQUFDLElBQUksRUFBRSxDQUFDLENBQUMsTUFBTSxFQUFFLENBQUMsQ0FBQyxJQUFJLEVBQUUsQ0FBQyxDQUFDLElBQUksRUFBRSxDQUFDLENBQUMsSUFBSSxFQUFFLENBQUMsQ0FBQyxFQUFFLENBQ3RELENBQUM7UUFDRixPQUFPLGlCQUFpQixDQUFDLEVBQUUsQ0FBQyxjQUFjLEVBQUUsQ0FBQztJQUMvQyxDQUFDO0lBRWtCLE1BQU0sQ0FBQyxPQUEwQztRQUNsRSxnQkFBZ0IsQ0FBQyxPQUFPLEVBQUUsSUFBSSxDQUFDLEVBQUUsRUFBRSxJQUFJLENBQUMsRUFBRSxDQUFDLENBQUM7SUFDOUMsQ0FBQztDQUNGO0FBckRnQjtJQURkLElBQUksQ0FBQyxHQUFHLEVBQUUsQ0FBQyxRQUFRLENBQUMsZUFBZSxDQUFDLDRCQUE0QixFQUFFLE1BQU0sQ0FBQyxDQUFDO21DQUN6QyJ9

package/lib/curves/Segment.d.ts

@@ -0,0 +1,9 @@
+import type { Vector2 } from '@twick/core';
+import type { CurvePoint } from './CurvePoint';
+export declare abstract class Segment {
+    abstract readonly points: Vector2[];
+    abstract draw(context: CanvasRenderingContext2D | Path2D, start: number, end: number, move: boolean): [CurvePoint, CurvePoint];
+    abstract getPoint(distance: number): CurvePoint;
+    abstract get arcLength(): number;
+}
+//# sourceMappingURL=Segment.d.ts.map

package/lib/curves/Segment.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"Segment.d.ts","sourceRoot":"","sources":["../../src/lib/curves/Segment.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAC,OAAO,EAAC,MAAM,aAAa,CAAC;AACzC,OAAO,KAAK,EAAC,UAAU,EAAC,MAAM,cAAc,CAAC;AAE7C,8BAAsB,OAAO;IAC3B,kBAAyB,MAAM,EAAE,OAAO,EAAE,CAAC;aAE3B,IAAI,CAClB,OAAO,EAAE,wBAAwB,GAAG,MAAM,EAC1C,KAAK,EAAE,MAAM,EACb,GAAG,EAAE,MAAM,EACX,IAAI,EAAE,OAAO,GACZ,CAAC,UAAU,EAAE,UAAU,CAAC;aAEX,QAAQ,CAAC,QAAQ,EAAE,MAAM,GAAG,UAAU;IAEtD,aAAoB,SAAS,IAAI,MAAM,CAAC;CACzC"}

package/lib/curves/Segment.js

@@ -0,0 +1,3 @@
+export class Segment {
+}
+//# sourceMappingURL=data:application/json;base64,eyJ2ZXJzaW9uIjozLCJmaWxlIjoiU2VnbWVudC5qcyIsInNvdXJjZVJvb3QiOiIiLCJzb3VyY2VzIjpbIi4uLy4uL3NyYy9saWIvY3VydmVzL1NlZ21lbnQudHMiXSwibmFtZXMiOltdLCJtYXBwaW5ncyI6IkFBR0EsTUFBTSxPQUFnQixPQUFPO0NBYTVCIn0=

package/lib/curves/UniformPolynomialCurveSampler.d.ts

@@ -0,0 +1,43 @@
+import type { CurvePoint } from './CurvePoint';
+import type { PolynomialSegment } from './PolynomialSegment';
+/**
+ * Class to uniformly sample points on a given polynomial curve.
+ *
+ * @remarks
+ * In order to uniformly sample points from non-linear curves, this sampler
+ * re-parameterizes the curve by arclength.
+ */
+export declare class UniformPolynomialCurveSampler {
+    private readonly curve;
+    private sampledDistances;
+    /**
+     * @param curve - The curve to sample
+     * @param samples - How many points to sample from the provided curve. The
+     *                  more points get sampled, the higher the resolution–and
+     *                  therefore precision–of the sampler.
+     */
+    constructor(curve: PolynomialSegment, samples?: number);
+    /**
+     * Discard all previously sampled points and resample the provided number of
+     * points from the curve.
+     *
+     * @param samples - The number of points to sample.
+     */
+    resample(samples: number): void;
+    /**
+     * Return the point at the provided distance along the sampled curve's
+     * arclength.
+     *
+     * @param distance - The distance along the curve's arclength for which to
+     *                   retrieve the point.
+     */
+    pointAtDistance(distance: number): CurvePoint;
+    /**
+     * Return the t value for the point at the provided distance along the sampled
+     * curve's arc length.
+     *
+     * @param distance - The distance along the arclength
+     */
+    distanceToT(distance: number): number;
+}
+//# sourceMappingURL=UniformPolynomialCurveSampler.d.ts.map

package/lib/curves/UniformPolynomialCurveSampler.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"UniformPolynomialCurveSampler.d.ts","sourceRoot":"","sources":["../../src/lib/curves/UniformPolynomialCurveSampler.ts"],"names":[],"mappings":"AAEA,OAAO,KAAK,EAAC,UAAU,EAAC,MAAM,cAAc,CAAC;AAC7C,OAAO,KAAK,EAAC,iBAAiB,EAAC,MAAM,qBAAqB,CAAC;AAE3D;;;;;;GAMG;AACH,qBAAa,6BAA6B;IAUtC,OAAO,CAAC,QAAQ,CAAC,KAAK;IATxB,OAAO,CAAC,gBAAgB,CAAgB;IAExC;;;;;OAKG;gBAEgB,KAAK,EAAE,iBAAiB,EACzC,OAAO,SAAK;IAKd;;;;;OAKG;IACI,QAAQ,CAAC,OAAO,EAAE,MAAM,GAAG,IAAI;IAsBtC;;;;;;OAMG;IACI,eAAe,CAAC,QAAQ,EAAE,MAAM,GAAG,UAAU;IAIpD;;;;;OAKG;IACI,WAAW,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;CAoB7C"}