@johnhenry/math 0.0.0

package/dist/Rotor4.js

@@ -0,0 +1,388 @@
+/**
+ * Rotor4 — a rotation operator in 4D space: `scalar + bivector (6) + pseudoscalar (1)`,
+ * 8 components total. This is the 4D generalization of {@link Quaternion}
+ * (which is really the 3D rotor, `scalar + bivector (3)`, that happens to
+ * have the same component count as a 3D vector — a coincidence unique to 3D
+ * that doesn't hold in 4D or any other dimension).
+ *
+ * **Double rotations.** A general 3D rotor always describes rotation in a
+ * single 2-plane. A general 4D rotor can describe a *double rotation*:
+ * simultaneous, independent rotation in two mutually orthogonal 2-planes at
+ * once — this is the origin of the nonzero pseudoscalar component, and it
+ * has no 3D analogue. See {@link Rotor4.multiply} composing two
+ * single-plane rotors in orthogonal planes for how a double rotation arises.
+ *
+ * **Normalization caveat.** {@link Rotor4.normalize} performs the naive
+ * quaternion-style renormalization (divide by magnitude). This is only
+ * exactly correct for *simple* rotors (single rotation plane) — for a
+ * genuine double rotation it does not lie on the same manifold as a simple
+ * rotor, so naive renormalization after numerical drift can silently produce
+ * an invalid rotor. {@link Rotor4.factor} and {@link Rotor4.renormalize} are
+ * the proper fix (Perwass-style rotor factorization): decompose into two
+ * simple rotors in orthogonal planes, each renormalizable exactly, then
+ * recompose — see the `N-Dimensional Rigid Body Dynamics` (SIGGRAPH 2020)
+ * paper §3.3 for why this matters for repeated physics integration.
+ */
+import { Bivector4 } from "./Bivector4.js";
+import { E1, E2, E3, E4, geometricProductMV, mvZero, PSEUDOSCALAR4 } from "./Clifford4Internal.js";
+import { bivectorFromMV, bivectorToMV, leftContraction, wedgeProduct } from "./GeometricAlgebra4.js";
+import { MatrixMath } from "./MatrixMath.js";
+import { Vec4 } from "./Vec4.js";
+import { Vector } from "./Vector.js";
+export class Rotor4 {
+    scalar;
+    bivector;
+    pseudoscalar;
+    static Identity = new Rotor4(1, Bivector4.Zero, 0);
+    constructor(scalar = 1, bivector = Bivector4.Zero, pseudoscalar = 0) {
+        this.scalar = scalar;
+        this.bivector = bivector;
+        this.pseudoscalar = pseudoscalar;
+    }
+    toMV() {
+        const mv = mvZero();
+        mv[0] = this.scalar;
+        const bmv = bivectorToMV(this.bivector);
+        for (let i = 0; i < bmv.length; i++)
+            mv[i] += bmv[i];
+        mv[PSEUDOSCALAR4] = this.pseudoscalar;
+        return mv;
+    }
+    static fromMV(mv) {
+        return new Rotor4(mv[0], bivectorFromMV(mv), mv[PSEUDOSCALAR4]);
+    }
+    add(r) {
+        return new Rotor4(this.scalar + r.scalar, this.bivector.add(r.bivector), this.pseudoscalar + r.pseudoscalar);
+    }
+    subtract(r) {
+        return new Rotor4(this.scalar - r.scalar, this.bivector.subtract(r.bivector), this.pseudoscalar - r.pseudoscalar);
+    }
+    negate() {
+        return new Rotor4(-this.scalar, this.bivector.negate(), -this.pseudoscalar);
+    }
+    scale(s) {
+        return new Rotor4(this.scalar * s, this.bivector.scale(s), this.pseudoscalar * s);
+    }
+    /** Component-wise inner product over all 8 components (scalar, 6 bivector, pseudoscalar) -- treats a Rotor4 as an 8-vector, exactly {@link Quaternion.dot}'s role in that class's own `slerp`. */
+    dot(r) {
+        return this.scalar * r.scalar + this.bivector.dot(r.bivector) + this.pseudoscalar * r.pseudoscalar;
+    }
+    /** The geometric product `this * r` — composes two rotations (apply `r` first, then `this`). */
+    multiply(r) {
+        return Rotor4.fromMV(geometricProductMV(this.toMV(), r.toMV()));
+    }
+    /**
+     * The reverse `R~`: reverses the order of vector factors, which negates the
+     * bivector (grade 2, picks up (-1)^1) but leaves the scalar (grade 0) and
+     * pseudoscalar (grade 4, picks up (-1)^6 = +1) unchanged. For a unit rotor,
+     * this is also the inverse — the 4D analogue of a quaternion's conjugate
+     * ({@link Quaternion.conjugate}). **Not** the inverse for a non-unit rotor
+     * (see {@link inverse}): a rotor that has drifted off the unit manifold
+     * during repeated composition (e.g. per-timestep physics integration,
+     * before renormalization runs) will silently produce a wrong,
+     * `|R|²`-scaled result from `reverse()` where an inverse was intended.
+     */
+    reverse() {
+        return new Rotor4(this.scalar, this.bivector.negate(), this.pseudoscalar);
+    }
+    /**
+     * Multiplicative inverse `R⁻¹ = R~ / |R|²`, distinct from {@link reverse}
+     * for a non-unit rotor — the 4D analogue of {@link Quaternion.inverse}.
+     * `this.multiply(this.inverse())` is `Identity` regardless of `this`'s
+     * magnitude; `this.multiply(this.reverse())` only is when `this` is
+     * already a unit rotor.
+     */
+    inverse() {
+        const m2 = this.magnitudeSquared;
+        if (m2 === 0)
+            throw new Error("The zero rotor has no inverse.");
+        return this.reverse().scale(1 / m2);
+    }
+    get magnitudeSquared() {
+        return this.scalar * this.scalar + this.bivector.magnitudeSquared + this.pseudoscalar * this.pseudoscalar;
+    }
+    get magnitude() {
+        return Math.sqrt(this.magnitudeSquared);
+    }
+    /**
+     * Naive renormalization (divide by magnitude). Exact for simple rotors;
+     * see the class-level doc comment for why this is *not* generally correct
+     * for double rotations, and where the real fix belongs.
+     */
+    normalize() {
+        const m = this.magnitude;
+        if (m === 0)
+            throw new Error("Cannot normalize the zero rotor.");
+        return this.scale(1 / m);
+    }
+    /** Apply this rotor to a vector via the sandwich product `R v R~`. Rotates `v`. */
+    apply(v) {
+        const vmv = mvZero();
+        vmv[E1] = v.x;
+        vmv[E2] = v.y;
+        vmv[E3] = v.z;
+        vmv[E4] = v.w;
+        const mv = geometricProductMV(geometricProductMV(this.toMV(), vmv), this.reverse().toMV());
+        return new Vec4(mv[E1], mv[E2], mv[E3], mv[E4]);
+    }
+    /** Apply this rotor to a bivector via the sandwich product `R B R~`. Rotates the 2-plane `B`. */
+    applyToBivector(b) {
+        const mv = geometricProductMV(geometricProductMV(this.toMV(), bivectorToMV(b)), this.reverse().toMV());
+        return bivectorFromMV(mv);
+    }
+    /**
+     * Perwass-style rotor factorization: decompose this rotor into two simple
+     * rotors `[R1, R2]` in mutually orthogonal planes such that `R2.multiply(R1)`
+     * has the same *action* as this rotor (equal up to the usual rotor
+     * double-cover sign, like any rotor comparison) — the proper fix for
+     * {@link normalize}'s documented limitation on double rotations. Works
+     * even when this rotor has drifted away from being a valid unit rotor
+     * (the intended use case: repeated physics integration).
+     *
+     * Method: find this rotor's two invariant 2-planes as eigenspaces of the
+     * *symmetric* part `S = M + Mᵀ` of its 4x4 matrix representation `M`
+     * (restricted to an invariant plane with rotation angle `θ`, `M + Mᵀ` is
+     * `2cos(θ)` times the identity on that plane — so same-angle-magnitude
+     * planes share an eigenvalue, and a repeated eigenvalue of `S` signals an
+     * isoclinic-type degeneracy where *any* orthonormal basis of that
+     * eigenspace is a genuinely valid invariant-plane choice, not just an
+     * approximation). `eigenSymmetric`'s eigenvector order isn't guaranteed to
+     * group same-plane pairs together when eigenvalues coincide (or are
+     * numerically indistinguishable), so the 3 possible pairings of the 4
+     * eigenvectors into 2 pairs are explicitly checked for genuine
+     * `M`-invariance (zero "leakage" into the other pair) before use, rather
+     * than assumed from sort order. Once a valid pair `(u, v)` spanning an
+     * invariant plane is known, the angle and correctly-oriented plane come
+     * directly from how `M` actually acts on `u` within that plane —
+     * `atan2(v·(Mu), u·(Mu))` — not from this rotor's own (possibly drifted)
+     * scalar/bivector/pseudoscalar components.
+     */
+    factor() {
+        const basis = [Vec4.Ex, Vec4.Ey, Vec4.Ez, Vec4.Ew];
+        const columns = basis.map((e) => this.apply(e).toArray());
+        const M = (i, j) => columns[j]?.[i];
+        const S = Vector.fromArray([0, 1, 2, 3].map((i) => Vector.fromArray([0, 1, 2, 3].map((j) => M(i, j) + M(j, i)))));
+        const { vectors } = MatrixMath.eigenSymmetric(S);
+        const eig = [0, 1, 2, 3].map((i) => new Vec4(vectors[0]?.[i], vectors[1]?.[i], vectors[2]?.[i], vectors[3]?.[i]));
+        const pairings = [
+            [
+                [0, 1],
+                [2, 3],
+            ],
+            [
+                [0, 2],
+                [1, 3],
+            ],
+            [
+                [0, 3],
+                [1, 2],
+            ],
+        ];
+        let bestLeak = Number.POSITIVE_INFINITY;
+        let best = pairings[0];
+        for (const [[i, j], [k, l]] of pairings) {
+            const ei = eig[i];
+            const ej = eig[j];
+            const ek = eig[k];
+            const el = eig[l];
+            const mi = this.apply(ei);
+            const mj = this.apply(ej);
+            const leak = Math.abs(mi.dot(ek)) + Math.abs(mi.dot(el)) + Math.abs(mj.dot(ek)) + Math.abs(mj.dot(el));
+            if (leak < bestLeak) {
+                bestLeak = leak;
+                best = [
+                    [i, j],
+                    [k, l],
+                ];
+            }
+        }
+        const planeAndAngle = (u, v) => {
+            const mu = this.apply(u);
+            const angle = Math.atan2(v.dot(mu), u.dot(mu));
+            const plane = wedgeProduct(v, u).normalize();
+            return Rotor4.fromBivectorAngle(plane, angle);
+        };
+        const [[a, b], [c, d]] = best;
+        return [planeAndAngle(eig[a], eig[b]), planeAndAngle(eig[c], eig[d])];
+    }
+    /**
+     * The proper (Perwass-style) fix for numerical drift: {@link factor} this
+     * rotor into two simple rotors — each exactly renormalizable, unlike a
+     * compound rotor — and recompose. Prefer this over {@link normalize} once
+     * a rotor may have accumulated drift through repeated integration and
+     * could be a genuine double rotation; `normalize` remains correct (and
+     * cheaper) for rotors known to be simple.
+     */
+    renormalize() {
+        // factor() reads this rotor's action via apply() (the sandwich product
+        // R v R~), which scales vectors by |R|^2 for a non-unit R -- corrupting
+        // the "M is a pure rotation matrix" assumption factor()'s eigenstructure
+        // analysis depends on. A cheap magnitude normalize first removes that
+        // radial drift; the exact plane/angle recovery in factor() then handles
+        // the part naive normalize alone gets wrong (the double-rotation
+        // manifold), on an already magnitude-correct input.
+        const [r1, r2] = this.normalize().factor();
+        return r2.multiply(r1);
+    }
+    /**
+     * Build the rotor for a rotation of `angle` radians in the 2-plane `plane`
+     * (analogous to {@link Quaternion.fromAxisAngle}, but a *plane* rather than
+     * an *axis* — axes only work as rotation generators in 3D). `plane` must be
+     * simple (see {@link Bivector4.isSimple}); to build a double rotation,
+     * construct two simple rotors in orthogonal planes and {@link multiply}
+     * them together.
+     */
+    static fromBivectorAngle(plane, angle) {
+        if (!plane.isSimple()) {
+            throw new Error("Rotor4.fromBivectorAngle requires a simple (single-plane) bivector. " +
+                "To build a double rotation, multiply() two simple rotors in orthogonal planes instead.");
+        }
+        const unitPlane = plane.normalize();
+        const half = angle / 2;
+        return new Rotor4(Math.cos(half), unitPlane.scale(Math.sin(half)), 0);
+    }
+    /**
+     * `exp(B)`: the rotor generated by the bivector `B` -- for a *simple* `B`,
+     * this is a rotation of `|B|` radians in `B`'s plane, matching
+     * {@link fromBivectorAngle}'s convention for the common physics case of an
+     * already-scaled generator (e.g. `angularVelocity * dt`, where the
+     * bivector's magnitude *is* the rotation angle and its direction *is* the
+     * plane). For a *compound* `B` (e.g. an angular velocity spinning
+     * simultaneously in two independent planes, entirely possible physically),
+     * the naive `cos|B| + sin|B|·B̂` formula this reduces to for simple `B` is
+     * wrong: `B² ≠ -|B|²` in general (it has a nonzero pseudoscalar term too),
+     * so this instead decomposes `B` into its two orthogonal simple
+     * components `m1·Π1 + m2·Π2` and returns `exp(m2·Π2)·exp(m1·Π1)` (valid
+     * since orthogonal-plane generators commute). Decomposition method:
+     * `B`'s invariant planes are the eigenspaces of `-A²`, where `A` is `B`'s
+     * contraction-action matrix (`A(v) = v⌋B`, antisymmetric, so `A²` is
+     * symmetric negative-semidefinite) -- the bivector analogue of how
+     * {@link factor} finds a *rotor's* invariant planes via `M + Mᵀ`, right
+     * down to needing the same explicit pairing-invariance check (eigenvalue
+     * sort order isn't reliable when eigenvalues coincide). Returns
+     * {@link Identity} for the zero bivector.
+     */
+    static exp(bivector) {
+        if (bivector.magnitude === 0)
+            return Rotor4.Identity;
+        const basis = [Vec4.Ex, Vec4.Ey, Vec4.Ez, Vec4.Ew];
+        const columns = basis.map((e) => leftContraction(e, bivector).toArray());
+        const A = (i, j) => columns[j]?.[i];
+        const A2 = Vector.fromArray([0, 1, 2, 3].map((i) => Vector.fromArray([0, 1, 2, 3].map((j) => {
+            let sum = 0;
+            for (let k = 0; k < 4; k++)
+                sum += A(i, k) * A(k, j);
+            return sum;
+        }))));
+        const { vectors } = MatrixMath.eigenSymmetric(A2);
+        const eig = [0, 1, 2, 3].map((i) => new Vec4(vectors[0]?.[i], vectors[1]?.[i], vectors[2]?.[i], vectors[3]?.[i]));
+        const pairings = [
+            [
+                [0, 1],
+                [2, 3],
+            ],
+            [
+                [0, 2],
+                [1, 3],
+            ],
+            [
+                [0, 3],
+                [1, 2],
+            ],
+        ];
+        let bestLeak = Number.POSITIVE_INFINITY;
+        let best = pairings[0];
+        for (const [[i, j], [k, l]] of pairings) {
+            const ei = eig[i];
+            const ej = eig[j];
+            const ek = eig[k];
+            const el = eig[l];
+            const ai = leftContraction(ei, bivector);
+            const aj = leftContraction(ej, bivector);
+            const leak = Math.abs(ai.dot(ek)) + Math.abs(ai.dot(el)) + Math.abs(aj.dot(ek)) + Math.abs(aj.dot(el));
+            if (leak < bestLeak) {
+                bestLeak = leak;
+                best = [
+                    [i, j],
+                    [k, l],
+                ];
+            }
+        }
+        const planeAndMagnitude = (u, v) => {
+            const au = leftContraction(u, bivector);
+            return { m: -v.dot(au), plane: wedgeProduct(v, u).normalize() };
+        };
+        const [[a, b], [c, d]] = best;
+        const p1 = planeAndMagnitude(eig[a], eig[b]);
+        const p2 = planeAndMagnitude(eig[c], eig[d]);
+        return Rotor4.fromBivectorAngle(p2.plane, p2.m).multiply(Rotor4.fromBivectorAngle(p1.plane, p1.m));
+    }
+    /**
+     * `log()`: the bivector `B` such that `exp(B) = this`, for a *simple*
+     * rotor -- the inverse of {@link exp}. Throws for a compound
+     * (double-rotation) rotor, same as {@link toBivectorAngle} (which this is
+     * built on): decomposing a general rotor's log is the Perwass-factorization
+     * work deferred to the physics layer.
+     */
+    log(epsilon = 1e-9) {
+        const { plane, angle } = this.toBivectorAngle(epsilon);
+        return plane.scale(angle);
+    }
+    /**
+     * Recover `(plane, angle)` for a *simple* rotor (pseudoscalar ≈ 0). Throws
+     * for a compound (double-rotation) rotor — see the class-level doc comment;
+     * decomposing a general rotor's log is the Perwass-factorization work
+     * deferred to the physics layer.
+     */
+    toBivectorAngle(epsilon = 1e-9) {
+        if (Math.abs(this.pseudoscalar) > epsilon) {
+            throw new Error("Rotor4.toBivectorAngle only supports simple rotors (pseudoscalar ≈ 0); " +
+                "this rotor is a compound double rotation.");
+        }
+        const bMag = this.bivector.magnitude;
+        const angle = 2 * Math.atan2(bMag, this.scalar);
+        if (bMag === 0)
+            return { plane: Bivector4.Zero, angle: 0 };
+        return { plane: this.bivector.normalize(), angle };
+    }
+    /**
+     * Spherical linear interpolation between two rotors, `t` in `[0, 1]` --
+     * the 4D analogue of {@link Quaternion.slerp}, same component-wise
+     * treatment (an 8-vector of scalar/bivector(6)/pseudoscalar, via
+     * {@link dot}) and the same near-parallel fallback to a normalized
+     * linear interpolation to avoid dividing by a near-zero `sin(theta0)`.
+     *
+     * **Same manifold caveat as {@link normalize}**: this is only exactly
+     * correct for interpolating between *simple* rotors (single rotation
+     * plane) — a compound double rotation doesn't lie on the sphere this
+     * construction assumes, so slerping one is an approximation, not a
+     * geodesic. For a double rotation, {@link factor} each endpoint into its
+     * two orthogonal simple rotors and slerp each pair independently instead.
+     */
+    static slerp(a, b, t) {
+        const ra = a.normalize();
+        let rb = b.normalize();
+        let cos = ra.dot(rb);
+        if (cos < 0) {
+            rb = rb.scale(-1);
+            cos = -cos;
+        }
+        if (cos > 0.9995)
+            return ra.add(rb.subtract(ra).scale(t)).normalize();
+        const theta0 = Math.acos(cos);
+        const theta = theta0 * t;
+        const sinTheta0 = Math.sin(theta0);
+        const s0 = Math.sin(theta0 - theta) / sinTheta0;
+        const s1 = Math.sin(theta) / sinTheta0;
+        return ra.scale(s0).add(rb.scale(s1));
+    }
+    equals(r, epsilon = 0) {
+        return ((epsilon === 0 ? this.scalar === r.scalar : Math.abs(this.scalar - r.scalar) <= epsilon) &&
+            this.bivector.equals(r.bivector, epsilon) &&
+            (epsilon === 0 ? this.pseudoscalar === r.pseudoscalar : Math.abs(this.pseudoscalar - r.pseudoscalar) <= epsilon));
+    }
+    toString() {
+        return `${this.scalar} + (${this.bivector.toString()}) + ${this.pseudoscalar}e1234`;
+    }
+}
+//# sourceMappingURL=Rotor4.js.map

package/dist/Rotor4.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Rotor4.js","sourceRoot":"","sources":["../src/Rotor4.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AACH,OAAO,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAC3C,OAAO,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,kBAAkB,EAAqB,MAAM,EAAE,aAAa,EAAE,MAAM,wBAAwB,CAAC;AACtH,OAAO,EAAE,cAAc,EAAE,YAAY,EAAE,eAAe,EAAE,YAAY,EAAE,MAAM,wBAAwB,CAAC;AACrG,OAAO,EAAE,UAAU,EAAE,MAAM,iBAAiB,CAAC;AAC7C,OAAO,EAAE,IAAI,EAAE,MAAM,WAAW,CAAC;AACjC,OAAO,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AAErC,MAAM,OAAO,MAAM;IACR,MAAM,CAAS;IACf,QAAQ,CAAY;IACpB,YAAY,CAAS;IAE9B,MAAM,CAAU,QAAQ,GAAG,IAAI,MAAM,CAAC,CAAC,EAAE,SAAS,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;IAE5D,YAAY,MAAM,GAAG,CAAC,EAAE,WAAsB,SAAS,CAAC,IAAI,EAAE,YAAY,GAAG,CAAC;QAC5E,IAAI,CAAC,MAAM,GAAG,MAAM,CAAC;QACrB,IAAI,CAAC,QAAQ,GAAG,QAAQ,CAAC;QACzB,IAAI,CAAC,YAAY,GAAG,YAAY,CAAC;IACnC,CAAC;IAEO,IAAI;QACV,MAAM,EAAE,GAAG,MAAM,EAAE,CAAC;QACpB,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC;QACpB,MAAM,GAAG,GAAG,YAAY,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;QACxC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,EAAE;YAAE,EAAE,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAW,CAAC;QAC/D,EAAE,CAAC,aAAa,CAAC,GAAG,IAAI,CAAC,YAAY,CAAC;QACtC,OAAO,EAAE,CAAC;IACZ,CAAC;IAEO,MAAM,CAAC,MAAM,CAAC,EAAgB;QACpC,OAAO,IAAI,MAAM,CAAC,EAAE,CAAC,CAAC,CAAW,EAAE,cAAc,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,aAAa,CAAW,CAAC,CAAC;IACtF,CAAC;IAED,GAAG,CAAC,CAAS;QACX,OAAO,IAAI,MAAM,CAAC,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,MAAM,EAAE,IAAI,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,YAAY,GAAG,CAAC,CAAC,YAAY,CAAC,CAAC;IAC/G,CAAC;IAED,QAAQ,CAAC,CAAS;QAChB,OAAO,IAAI,MAAM,CAAC,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,MAAM,EAAE,IAAI,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC,CAAC,QAAQ,CAAC,EAAE,IAAI,CAAC,YAAY,GAAG,CAAC,CAAC,YAAY,CAAC,CAAC;IACpH,CAAC;IAED,MAAM;QACJ,OAAO,IAAI,MAAM,CAAC,CAAC,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC,QAAQ,CAAC,MAAM,EAAE,EAAE,CAAC,IAAI,CAAC,YAAY,CAAC,CAAC;IAC9E,CAAC;IAED,KAAK,CAAC,CAAS;QACb,OAAO,IAAI,MAAM,CAAC,IAAI,CAAC,MAAM,GAAG,CAAC,EAAE,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC;IACpF,CAAC;IAED,kMAAkM;IAClM,GAAG,CAAC,CAAS;QACX,OAAO,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,MAAM,GAAG,IAAI,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,QAAQ,CAAC,GAAG,IAAI,CAAC,YAAY,GAAG,CAAC,CAAC,YAAY,CAAC;IACrG,CAAC;IAED,gGAAgG;IAChG,QAAQ,CAAC,CAAS;QAChB,OAAO,MAAM,CAAC,MAAM,CAAC,kBAAkB,CAAC,IAAI,CAAC,IAAI,EAAE,EAAE,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;IAClE,CAAC;IAED;;;;;;;;;;OAUG;IACH,OAAO;QACL,OAAO,IAAI,MAAM,CAAC,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC,QAAQ,CAAC,MAAM,EAAE,EAAE,IAAI,CAAC,YAAY,CAAC,CAAC;IAC5E,CAAC;IAED;;;;;;OAMG;IACH,OAAO;QACL,MAAM,EAAE,GAAG,IAAI,CAAC,gBAAgB,CAAC;QACjC,IAAI,EAAE,KAAK,CAAC;YAAE,MAAM,IAAI,KAAK,CAAC,gCAAgC,CAAC,CAAC;QAChE,OAAO,IAAI,CAAC,OAAO,EAAE,CAAC,KAAK,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;IACtC,CAAC;IAED,IAAI,gBAAgB;QAClB,OAAO,IAAI,CAAC,MAAM,GAAG,IAAI,CAAC,MAAM,GAAG,IAAI,CAAC,QAAQ,CAAC,gBAAgB,GAAG,IAAI,CAAC,YAAY,GAAG,IAAI,CAAC,YAAY,CAAC;IAC5G,CAAC;IAED,IAAI,SAAS;QACX,OAAO,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,gBAAgB,CAAC,CAAC;IAC1C,CAAC;IAED;;;;OAIG;IACH,SAAS;QACP,MAAM,CAAC,GAAG,IAAI,CAAC,SAAS,CAAC;QACzB,IAAI,CAAC,KAAK,CAAC;YAAE,MAAM,IAAI,KAAK,CAAC,kCAAkC,CAAC,CAAC;QACjE,OAAO,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC3B,CAAC;IAED,mFAAmF;IACnF,KAAK,CAAC,CAAO;QACX,MAAM,GAAG,GAAG,MAAM,EAAE,CAAC;QACrB,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACd,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACd,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACd,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACd,MAAM,EAAE,GAAG,kBAAkB,CAAC,kBAAkB,CAAC,IAAI,CAAC,IAAI,EAAE,EAAE,GAAG,CAAC,EAAE,IAAI,CAAC,OAAO,EAAE,CAAC,IAAI,EAAE,CAAC,CAAC;QAC3F,OAAO,IAAI,IAAI,CAAC,EAAE,CAAC,EAAE,CAAW,EAAE,EAAE,CAAC,EAAE,CAAW,EAAE,EAAE,CAAC,EAAE,CAAW,EAAE,EAAE,CAAC,EAAE,CAAW,CAAC,CAAC;IAC1F,CAAC;IAED,iGAAiG;IACjG,eAAe,CAAC,CAAY;QAC1B,MAAM,EAAE,GAAG,kBAAkB,CAAC,kBAAkB,CAAC,IAAI,CAAC,IAAI,EAAE,EAAE,YAAY,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,OAAO,EAAE,CAAC,IAAI,EAAE,CAAC,CAAC;QACvG,OAAO,cAAc,CAAC,EAAE,CAAC,CAAC;IAC5B,CAAC;IAED;;;;;;;;;;;;;;;;;;;;;;;;;;OA0BG;IACH,MAAM;QACJ,MAAM,KAAK,GAAG,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;QACnD,MAAM,OAAO,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,OAAO,EAAE,CAAC,CAAC;QAC1D,MAAM,CAAC,GAAG,CAAC,CAAS,EAAE,CAAS,EAAU,EAAE,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAW,CAAC;QACtE,MAAM,CAAC,GAAG,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAClH,MAAM,EAAE,OAAO,EAAE,GAAG,UAAU,CAAC,cAAc,CAAC,CAAC,CAAC,CAAC;QACjD,MAAM,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAC1B,CAAC,CAAC,EAAE,EAAE,CACJ,IAAI,IAAI,CACN,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAW,EACzB,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAW,EACzB,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAW,EACzB,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAW,CAC1B,CACJ,CAAC;QAEF,MAAM,QAAQ,GAA2C;YACvD;gBACE,CAAC,CAAC,EAAE,CAAC,CAAC;gBACN,CAAC,CAAC,EAAE,CAAC,CAAC;aACP;YACD;gBACE,CAAC,CAAC,EAAE,CAAC,CAAC;gBACN,CAAC,CAAC,EAAE,CAAC,CAAC;aACP;YACD;gBACE,CAAC,CAAC,EAAE,CAAC,CAAC;gBACN,CAAC,CAAC,EAAE,CAAC,CAAC;aACP;SACF,CAAC;QACF,IAAI,QAAQ,GAAG,MAAM,CAAC,iBAAiB,CAAC;QACxC,IAAI,IAAI,GAAyC,QAAQ,CAAC,CAAC,CAAyC,CAAC;QACrG,KAAK,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,QAAQ,EAAE,CAAC;YACxC,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAS,CAAC;YAC1B,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAS,CAAC;YAC1B,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAS,CAAC;YAC1B,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAS,CAAC;YAC1B,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC;YAC1B,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC;YAC1B,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;YACvG,IAAI,IAAI,GAAG,QAAQ,EAAE,CAAC;gBACpB,QAAQ,GAAG,IAAI,CAAC;gBAChB,IAAI,GAAG;oBACL,CAAC,CAAC,EAAE,CAAC,CAAC;oBACN,CAAC,CAAC,EAAE,CAAC,CAAC;iBACP,CAAC;YACJ,CAAC;QACH,CAAC;QAED,MAAM,aAAa,GAAG,CAAC,CAAO,EAAE,CAAO,EAAU,EAAE;YACjD,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;YACzB,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;YAC/C,MAAM,KAAK,GAAG,YAAY,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,SAAS,EAAE,CAAC;YAC7C,OAAO,MAAM,CAAC,iBAAiB,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;QAChD,CAAC,CAAC;QACF,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;QAC9B,OAAO,CAAC,aAAa,CAAC,GAAG,CAAC,CAAC,CAAS,EAAE,GAAG,CAAC,CAAC,CAAS,CAAC,EAAE,aAAa,CAAC,GAAG,CAAC,CAAC,CAAS,EAAE,GAAG,CAAC,CAAC,CAAS,CAAC,CAAC,CAAC;IACxG,CAAC;IAED;;;;;;;OAOG;IACH,WAAW;QACT,uEAAuE;QACvE,wEAAwE;QACxE,yEAAyE;QACzE,sEAAsE;QACtE,wEAAwE;QACxE,iEAAiE;QACjE,oDAAoD;QACpD,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,IAAI,CAAC,SAAS,EAAE,CAAC,MAAM,EAAE,CAAC;QAC3C,OAAO,EAAE,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC;IACzB,CAAC;IAED;;;;;;;OAOG;IACH,MAAM,CAAC,iBAAiB,CAAC,KAAgB,EAAE,KAAa;QACtD,IAAI,CAAC,KAAK,CAAC,QAAQ,EAAE,EAAE,CAAC;YACtB,MAAM,IAAI,KAAK,CACb,sEAAsE;gBACpE,wFAAwF,CAC3F,CAAC;QACJ,CAAC;QACD,MAAM,SAAS,GAAG,KAAK,CAAC,SAAS,EAAE,CAAC;QACpC,MAAM,IAAI,GAAG,KAAK,GAAG,CAAC,CAAC;QACvB,OAAO,IAAI,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,SAAS,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACxE,CAAC;IAED;;;;;;;;;;;;;;;;;;;;OAoBG;IACH,MAAM,CAAC,GAAG,CAAC,QAAmB;QAC5B,IAAI,QAAQ,CAAC,SAAS,KAAK,CAAC;YAAE,OAAO,MAAM,CAAC,QAAQ,CAAC;QAErD,MAAM,KAAK,GAAG,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;QACnD,MAAM,OAAO,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,eAAe,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,OAAO,EAAE,CAAC,CAAC;QACzE,MAAM,CAAC,GAAG,CAAC,CAAS,EAAE,CAAS,EAAU,EAAE,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAW,CAAC;QACtE,MAAM,EAAE,GAAG,MAAM,CAAC,SAAS,CACzB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CACrB,MAAM,CAAC,SAAS,CACd,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE;YACrB,IAAI,GAAG,GAAG,CAAC,CAAC;YACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,GAAG,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YACrD,OAAO,GAAG,CAAC;QACb,CAAC,CAAC,CACH,CACF,CACF,CAAC;QACF,MAAM,EAAE,OAAO,EAAE,GAAG,UAAU,CAAC,cAAc,CAAC,EAAE,CAAC,CAAC;QAClD,MAAM,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAC1B,CAAC,CAAC,EAAE,EAAE,CACJ,IAAI,IAAI,CACN,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAW,EACzB,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAW,EACzB,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAW,EACzB,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAW,CAC1B,CACJ,CAAC;QAEF,MAAM,QAAQ,GAA2C;YACvD;gBACE,CAAC,CAAC,EAAE,CAAC,CAAC;gBACN,CAAC,CAAC,EAAE,CAAC,CAAC;aACP;YACD;gBACE,CAAC,CAAC,EAAE,CAAC,CAAC;gBACN,CAAC,CAAC,EAAE,CAAC,CAAC;aACP;YACD;gBACE,CAAC,CAAC,EAAE,CAAC,CAAC;gBACN,CAAC,CAAC,EAAE,CAAC,CAAC;aACP;SACF,CAAC;QACF,IAAI,QAAQ,GAAG,MAAM,CAAC,iBAAiB,CAAC;QACxC,IAAI,IAAI,GAAyC,QAAQ,CAAC,CAAC,CAAyC,CAAC;QACrG,KAAK,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,QAAQ,EAAE,CAAC;YACxC,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAS,CAAC;YAC1B,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAS,CAAC;YAC1B,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAS,CAAC;YAC1B,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAS,CAAC;YAC1B,MAAM,EAAE,GAAG,eAAe,CAAC,EAAE,EAAE,QAAQ,CAAC,CAAC;YACzC,MAAM,EAAE,GAAG,eAAe,CAAC,EAAE,EAAE,QAAQ,CAAC,CAAC;YACzC,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;YACvG,IAAI,IAAI,GAAG,QAAQ,EAAE,CAAC;gBACpB,QAAQ,GAAG,IAAI,CAAC;gBAChB,IAAI,GAAG;oBACL,CAAC,CAAC,EAAE,CAAC,CAAC;oBACN,CAAC,CAAC,EAAE,CAAC,CAAC;iBACP,CAAC;YACJ,CAAC;QACH,CAAC;QAED,MAAM,iBAAiB,GAAG,CAAC,CAAO,EAAE,CAAO,EAAmC,EAAE;YAC9E,MAAM,EAAE,GAAG,eAAe,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC;YACxC,OAAO,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,KAAK,EAAE,YAAY,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,SAAS,EAAE,EAAE,CAAC;QAClE,CAAC,CAAC;QACF,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;QAC9B,MAAM,EAAE,GAAG,iBAAiB,CAAC,GAAG,CAAC,CAAC,CAAS,EAAE,GAAG,CAAC,CAAC,CAAS,CAAC,CAAC;QAC7D,MAAM,EAAE,GAAG,iBAAiB,CAAC,GAAG,CAAC,CAAC,CAAS,EAAE,GAAG,CAAC,CAAC,CAAS,CAAC,CAAC;QAC7D,OAAO,MAAM,CAAC,iBAAiB,CAAC,EAAE,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,MAAM,CAAC,iBAAiB,CAAC,EAAE,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACrG,CAAC;IAED;;;;;;OAMG;IACH,GAAG,CAAC,OAAO,GAAG,IAAI;QAChB,MAAM,EAAE,KAAK,EAAE,KAAK,EAAE,GAAG,IAAI,CAAC,eAAe,CAAC,OAAO,CAAC,CAAC;QACvD,OAAO,KAAK,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC;IAC5B,CAAC;IAED;;;;;OAKG;IACH,eAAe,CAAC,OAAO,GAAG,IAAI;QAC5B,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,YAAY,CAAC,GAAG,OAAO,EAAE,CAAC;YAC1C,MAAM,IAAI,KAAK,CACb,yEAAyE;gBACvE,2CAA2C,CAC9C,CAAC;QACJ,CAAC;QACD,MAAM,IAAI,GAAG,IAAI,CAAC,QAAQ,CAAC,SAAS,CAAC;QACrC,MAAM,KAAK,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;QAChD,IAAI,IAAI,KAAK,CAAC;YAAE,OAAO,EAAE,KAAK,EAAE,SAAS,CAAC,IAAI,EAAE,KAAK,EAAE,CAAC,EAAE,CAAC;QAC3D,OAAO,EAAE,KAAK,EAAE,IAAI,CAAC,QAAQ,CAAC,SAAS,EAAE,EAAE,KAAK,EAAE,CAAC;IACrD,CAAC;IAED;;;;;;;;;;;;;OAaG;IACH,MAAM,CAAC,KAAK,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS;QAC1C,MAAM,EAAE,GAAG,CAAC,CAAC,SAAS,EAAE,CAAC;QACzB,IAAI,EAAE,GAAG,CAAC,CAAC,SAAS,EAAE,CAAC;QACvB,IAAI,GAAG,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;QACrB,IAAI,GAAG,GAAG,CAAC,EAAE,CAAC;YACZ,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;YAClB,GAAG,GAAG,CAAC,GAAG,CAAC;QACb,CAAC;QACD,IAAI,GAAG,GAAG,MAAM;YAAE,OAAO,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,SAAS,EAAE,CAAC;QACtE,MAAM,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QAC9B,MAAM,KAAK,GAAG,MAAM,GAAG,CAAC,CAAC;QACzB,MAAM,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC;QACnC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,KAAK,CAAC,GAAG,SAAS,CAAC;QAChD,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,SAAS,CAAC;QACvC,OAAO,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC;IACxC,CAAC;IAED,MAAM,CAAC,CAAS,EAAE,OAAO,GAAG,CAAC;QAC3B,OAAO,CACL,CAAC,OAAO,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,KAAK,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,MAAM,CAAC,IAAI,OAAO,CAAC;YACxF,IAAI,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,QAAQ,EAAE,OAAO,CAAC;YACzC,CAAC,OAAO,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,YAAY,KAAK,CAAC,CAAC,YAAY,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,YAAY,GAAG,CAAC,CAAC,YAAY,CAAC,IAAI,OAAO,CAAC,CACjH,CAAC;IACJ,CAAC;IAED,QAAQ;QACN,OAAO,GAAG,IAAI,CAAC,MAAM,OAAO,IAAI,CAAC,QAAQ,CAAC,QAAQ,EAAE,OAAO,IAAI,CAAC,YAAY,OAAO,CAAC;IACtF,CAAC"}

package/dist/SpecialFunctions.d.ts

@@ -0,0 +1,32 @@
+/**
+ * SpecialFunctions — the higher transcendental functions that statistics and
+ * analysis rely on: the gamma and log-gamma functions (Lanczos), the beta
+ * function, the error function `erf`/`erfc`, and the regularised incomplete
+ * gamma and beta functions (which give the distribution CDFs in
+ * {@link Distributions} exact closed forms rather than numeric integration).
+ */
+export declare class SpecialFunctions {
+    /** The gamma function Γ(x) (Lanczos approximation). */
+    static gamma(x: number): number;
+    /** The natural log of |Γ(x)| — stable for large arguments. */
+    static lnGamma(x: number): number;
+    /** The beta function B(a, b) = Γ(a)Γ(b)/Γ(a+b). */
+    static beta(a: number, b: number): number;
+    /** The error function erf(x). */
+    static erf(x: number): number;
+    /** The complementary error function erfc(x) = 1 − erf(x). */
+    static erfc(x: number): number;
+    /**
+     * The regularised lower incomplete gamma function P(a, x) = γ(a, x)/Γ(a).
+     * (This is the CDF of the gamma/chi-square distributions.)
+     */
+    static regularizedGammaP(a: number, x: number): number;
+    /** The regularised upper incomplete gamma function Q(a, x) = 1 − P(a, x). */
+    static regularizedGammaQ(a: number, x: number): number;
+    /**
+     * The regularised incomplete beta function I_x(a, b). (This is the CDF of the
+     * beta distribution and underlies the Student-t and F distributions.)
+     */
+    static regularizedIncompleteBeta(x: number, a: number, b: number): number;
+}
+//# sourceMappingURL=SpecialFunctions.d.ts.map

package/dist/SpecialFunctions.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"SpecialFunctions.d.ts","sourceRoot":"","sources":["../src/SpecialFunctions.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AAQH,qBAAa,gBAAgB;IAC3B,uDAAuD;IACvD,MAAM,CAAC,KAAK,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM;IAS/B,8DAA8D;IAC9D,MAAM,CAAC,OAAO,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM;IASjC,mDAAmD;IACnD,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM;IAIzC,iCAAiC;IACjC,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM;IAY7B,6DAA6D;IAC7D,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM;IAI9B;;;OAGG;IACH,MAAM,CAAC,iBAAiB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM;IAsCtD,6EAA6E;IAC7E,MAAM,CAAC,iBAAiB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM;IAItD;;;OAGG;IACH,MAAM,CAAC,yBAAyB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM;CAS1E"}

package/dist/SpecialFunctions.js

@@ -0,0 +1,161 @@
+/**
+ * SpecialFunctions — the higher transcendental functions that statistics and
+ * analysis rely on: the gamma and log-gamma functions (Lanczos), the beta
+ * function, the error function `erf`/`erfc`, and the regularised incomplete
+ * gamma and beta functions (which give the distribution CDFs in
+ * {@link Distributions} exact closed forms rather than numeric integration).
+ */
+const LANCZOS_G = 7;
+const LANCZOS = [
+    0.9999999999998099, 676.5203681218851, -1259.1392167224028, 771.3234287776531, -176.6150291621406, 12.507343278686905,
+    -0.13857109526572012, 9.984369578019572e-6, 1.5056327351493116e-7,
+];
+export class SpecialFunctions {
+    /** The gamma function Γ(x) (Lanczos approximation). */
+    static gamma(x) {
+        if (x < 0.5)
+            return Math.PI / (Math.sin(Math.PI * x) * SpecialFunctions.gamma(1 - x));
+        const z = x - 1;
+        let a = LANCZOS[0];
+        const t = z + LANCZOS_G + 0.5;
+        for (let i = 1; i < LANCZOS.length; i++)
+            a += LANCZOS[i] / (z + i);
+        return Math.sqrt(2 * Math.PI) * t ** (z + 0.5) * Math.exp(-t) * a;
+    }
+    /** The natural log of |Γ(x)| — stable for large arguments. */
+    static lnGamma(x) {
+        if (x < 0.5)
+            return Math.log(Math.abs(Math.PI / Math.sin(Math.PI * x))) - SpecialFunctions.lnGamma(1 - x);
+        const z = x - 1;
+        let a = LANCZOS[0];
+        const t = z + LANCZOS_G + 0.5;
+        for (let i = 1; i < LANCZOS.length; i++)
+            a += LANCZOS[i] / (z + i);
+        return 0.5 * Math.log(2 * Math.PI) + (z + 0.5) * Math.log(t) - t + Math.log(a);
+    }
+    /** The beta function B(a, b) = Γ(a)Γ(b)/Γ(a+b). */
+    static beta(a, b) {
+        return Math.exp(SpecialFunctions.lnGamma(a) + SpecialFunctions.lnGamma(b) - SpecialFunctions.lnGamma(a + b));
+    }
+    /** The error function erf(x). */
+    static erf(x) {
+        const sign = x < 0 ? -1 : 1;
+        const ax = Math.abs(x);
+        const t = 1 / (1 + 0.3275911 * ax);
+        const y = 1 -
+            ((((1.061405429 * t - 1.453152027) * t + 1.421413741) * t - 0.284496736) * t + 0.254829592) *
+                t *
+                Math.exp(-ax * ax);
+        return sign * y;
+    }
+    /** The complementary error function erfc(x) = 1 − erf(x). */
+    static erfc(x) {
+        return 1 - SpecialFunctions.erf(x);
+    }
+    /**
+     * The regularised lower incomplete gamma function P(a, x) = γ(a, x)/Γ(a).
+     * (This is the CDF of the gamma/chi-square distributions.)
+     */
+    static regularizedGammaP(a, x) {
+        if (x < 0 || a <= 0)
+            return Number.NaN;
+        if (x === 0)
+            return 0;
+        if (x < a + 1) {
+            // series expansion
+            let ap = a;
+            let sum = 1 / a;
+            let del = sum;
+            for (let n = 0; n < 200; n++) {
+                ap++;
+                del *= x / ap;
+                sum += del;
+                if (Math.abs(del) < Math.abs(sum) * 1e-15)
+                    break;
+            }
+            return sum * Math.exp(-x + a * Math.log(x) - SpecialFunctions.lnGamma(a));
+        }
+        // continued fraction for the upper part Q, then P = 1 - Q
+        const tiny = 1e-300;
+        let b = x + 1 - a;
+        let c = 1 / tiny;
+        let d = 1 / b;
+        let h = d;
+        for (let i = 1; i < 200; i++) {
+            const an = -i * (i - a);
+            b += 2;
+            d = an * d + b;
+            if (Math.abs(d) < tiny)
+                d = tiny;
+            c = b + an / c;
+            if (Math.abs(c) < tiny)
+                c = tiny;
+            d = 1 / d;
+            const del = d * c;
+            h *= del;
+            if (Math.abs(del - 1) < 1e-15)
+                break;
+        }
+        const q = Math.exp(-x + a * Math.log(x) - SpecialFunctions.lnGamma(a)) * h;
+        return 1 - q;
+    }
+    /** The regularised upper incomplete gamma function Q(a, x) = 1 − P(a, x). */
+    static regularizedGammaQ(a, x) {
+        return 1 - SpecialFunctions.regularizedGammaP(a, x);
+    }
+    /**
+     * The regularised incomplete beta function I_x(a, b). (This is the CDF of the
+     * beta distribution and underlies the Student-t and F distributions.)
+     */
+    static regularizedIncompleteBeta(x, a, b) {
+        if (x <= 0)
+            return 0;
+        if (x >= 1)
+            return 1;
+        const lbeta = SpecialFunctions.lnGamma(a + b) - SpecialFunctions.lnGamma(a) - SpecialFunctions.lnGamma(b);
+        const front = Math.exp(lbeta + a * Math.log(x) + b * Math.log(1 - x));
+        // Use the continued fraction that converges fastest, swapping (a,b,x) as needed.
+        if (x < (a + 1) / (a + b + 2))
+            return (front * betaContinuedFraction(a, b, x)) / a;
+        return 1 - (front * betaContinuedFraction(b, a, 1 - x)) / b;
+    }
+}
+/** Lentz's continued fraction for the incomplete beta (Numerical Recipes `betacf`). */
+function betaContinuedFraction(a, b, x) {
+    const tiny = 1e-300;
+    const qab = a + b;
+    const qap = a + 1;
+    const qam = a - 1;
+    let c = 1;
+    let d = 1 - (qab * x) / qap;
+    if (Math.abs(d) < tiny)
+        d = tiny;
+    d = 1 / d;
+    let h = d;
+    for (let m = 1; m < 200; m++) {
+        const m2 = 2 * m;
+        let aa = (m * (b - m) * x) / ((qam + m2) * (a + m2));
+        d = 1 + aa * d;
+        if (Math.abs(d) < tiny)
+            d = tiny;
+        c = 1 + aa / c;
+        if (Math.abs(c) < tiny)
+            c = tiny;
+        d = 1 / d;
+        h *= d * c;
+        aa = (-(a + m) * (qab + m) * x) / ((a + m2) * (qap + m2));
+        d = 1 + aa * d;
+        if (Math.abs(d) < tiny)
+            d = tiny;
+        c = 1 + aa / c;
+        if (Math.abs(c) < tiny)
+            c = tiny;
+        d = 1 / d;
+        const del = d * c;
+        h *= del;
+        if (Math.abs(del - 1) < 1e-15)
+            break;
+    }
+    return h;
+}
+//# sourceMappingURL=SpecialFunctions.js.map

package/dist/SpecialFunctions.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"SpecialFunctions.js","sourceRoot":"","sources":["../src/SpecialFunctions.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AAEH,MAAM,SAAS,GAAG,CAAC,CAAC;AACpB,MAAM,OAAO,GAAG;IACd,kBAAkB,EAAE,iBAAiB,EAAE,CAAC,kBAAkB,EAAE,iBAAiB,EAAE,CAAC,iBAAiB,EAAE,kBAAkB;IACrH,CAAC,mBAAmB,EAAE,oBAAoB,EAAE,qBAAqB;CAClE,CAAC;AAEF,MAAM,OAAO,gBAAgB;IAC3B,uDAAuD;IACvD,MAAM,CAAC,KAAK,CAAC,CAAS;QACpB,IAAI,CAAC,GAAG,GAAG;YAAE,OAAO,IAAI,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,gBAAgB,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACtF,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAChB,IAAI,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;QACnB,MAAM,CAAC,GAAG,CAAC,GAAG,SAAS,GAAG,GAAG,CAAC;QAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,EAAE;YAAE,CAAC,IAAI,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QACnE,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACpE,CAAC;IAED,8DAA8D;IAC9D,MAAM,CAAC,OAAO,CAAC,CAAS;QACtB,IAAI,CAAC,GAAG,GAAG;YAAE,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,gBAAgB,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC1G,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAChB,IAAI,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;QACnB,MAAM,CAAC,GAAG,CAAC,GAAG,SAAS,GAAG,GAAG,CAAC;QAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,EAAE;YAAE,CAAC,IAAI,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QACnE,OAAO,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACjF,CAAC;IAED,mDAAmD;IACnD,MAAM,CAAC,IAAI,CAAC,CAAS,EAAE,CAAS;QAC9B,OAAO,IAAI,CAAC,GAAG,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC,CAAC,GAAG,gBAAgB,CAAC,OAAO,CAAC,CAAC,CAAC,GAAG,gBAAgB,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAC/G,CAAC;IAED,iCAAiC;IACjC,MAAM,CAAC,GAAG,CAAC,CAAS;QAClB,MAAM,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC5B,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACvB,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,SAAS,GAAG,EAAE,CAAC,CAAC;QACnC,MAAM,CAAC,GACL,CAAC;YACD,CAAC,CAAC,CAAC,CAAC,WAAW,GAAG,CAAC,GAAG,WAAW,CAAC,GAAG,CAAC,GAAG,WAAW,CAAC,GAAG,CAAC,GAAG,WAAW,CAAC,GAAG,CAAC,GAAG,WAAW,CAAC;gBACzF,CAAC;gBACD,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;QACvB,OAAO,IAAI,GAAG,CAAC,CAAC;IAClB,CAAC;IAED,6DAA6D;IAC7D,MAAM,CAAC,IAAI,CAAC,CAAS;QACnB,OAAO,CAAC,GAAG,gBAAgB,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACrC,CAAC;IAED;;;OAGG;IACH,MAAM,CAAC,iBAAiB,CAAC,CAAS,EAAE,CAAS;QAC3C,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC;YAAE,OAAO,MAAM,CAAC,GAAG,CAAC;QACvC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,CAAC,CAAC;QACtB,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC;YACd,mBAAmB;YACnB,IAAI,EAAE,GAAG,CAAC,CAAC;YACX,IAAI,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;YAChB,IAAI,GAAG,GAAG,GAAG,CAAC;YACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC7B,EAAE,EAAE,CAAC;gBACL,GAAG,IAAI,CAAC,GAAG,EAAE,CAAC;gBACd,GAAG,IAAI,GAAG,CAAC;gBACX,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,KAAK;oBAAE,MAAM;YACnD,CAAC;YACD,OAAO,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,gBAAgB,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;QAC5E,CAAC;QACD,0DAA0D;QAC1D,MAAM,IAAI,GAAG,MAAM,CAAC;QACpB,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAClB,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC;QACjB,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACd,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;YAC7B,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACxB,CAAC,IAAI,CAAC,CAAC;YACP,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;YACf,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI;gBAAE,CAAC,GAAG,IAAI,CAAC;YACjC,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;YACf,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI;gBAAE,CAAC,GAAG,IAAI,CAAC;YACjC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YACV,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;YAClB,CAAC,IAAI,GAAG,CAAC;YACT,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,KAAK;gBAAE,MAAM;QACvC,CAAC;QACD,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,gBAAgB,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QAC3E,OAAO,CAAC,GAAG,CAAC,CAAC;IACf,CAAC;IAED,6EAA6E;IAC7E,MAAM,CAAC,iBAAiB,CAAC,CAAS,EAAE,CAAS;QAC3C,OAAO,CAAC,GAAG,gBAAgB,CAAC,iBAAiB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACtD,CAAC;IAED;;;OAGG;IACH,MAAM,CAAC,yBAAyB,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS;QAC9D,IAAI,CAAC,IAAI,CAAC;YAAE,OAAO,CAAC,CAAC;QACrB,IAAI,CAAC,IAAI,CAAC;YAAE,OAAO,CAAC,CAAC;QACrB,MAAM,KAAK,GAAG,gBAAgB,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,gBAAgB,CAAC,OAAO,CAAC,CAAC,CAAC,GAAG,gBAAgB,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC;QAC1G,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACtE,iFAAiF;QACjF,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YAAE,OAAO,CAAC,KAAK,GAAG,qBAAqB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACnF,OAAO,CAAC,GAAG,CAAC,KAAK,GAAG,qBAAqB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAC9D,CAAC;CACF;AAED,uFAAuF;AACvF,SAAS,qBAAqB,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS;IAC5D,MAAM,IAAI,GAAG,MAAM,CAAC;IACpB,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAClB,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAClB,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAClB,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;IAC5B,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI;QAAE,CAAC,GAAG,IAAI,CAAC;IACjC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACV,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;QAC7B,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QACjB,IAAI,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACrD,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QACf,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI;YAAE,CAAC,GAAG,IAAI,CAAC;QACjC,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QACf,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI;YAAE,CAAC,GAAG,IAAI,CAAC;QACjC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACV,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACX,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC;QAC1D,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QACf,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI;YAAE,CAAC,GAAG,IAAI,CAAC;QACjC,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QACf,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI;YAAE,CAAC,GAAG,IAAI,CAAC;QACjC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACV,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;QAClB,CAAC,IAAI,GAAG,CAAC;QACT,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,KAAK;YAAE,MAAM;IACvC,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC"}

package/dist/SpecialOperator.d.ts

@@ -0,0 +1,21 @@
+/**
+ * SpecialOperator — a named symbolic operator (`+`, `-`, `*`, …) pairing a
+ * textual representation (`rep`) with the name of the function that implements it
+ * (`funct`). Ported from Mallory's ActionScript `SpecialOperator`; used by the
+ * expression evaluator.
+ */
+export declare class SpecialOperator {
+    readonly rep: string;
+    readonly funct: string | null;
+    constructor(rep: string, funct: string | null);
+    toString(): string;
+    static readonly Plus: SpecialOperator;
+    static readonly Minus: SpecialOperator;
+    static readonly Times: SpecialOperator;
+    static readonly Divided: SpecialOperator;
+    static readonly Power: SpecialOperator;
+    static readonly Modulous: SpecialOperator;
+    static readonly Elipses: SpecialOperator;
+    static readonly System: SpecialOperator;
+}
+//# sourceMappingURL=SpecialOperator.d.ts.map

package/dist/SpecialOperator.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"SpecialOperator.d.ts","sourceRoot":"","sources":["../src/SpecialOperator.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AACH,qBAAa,eAAe;IAC1B,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,QAAQ,CAAC,KAAK,EAAE,MAAM,GAAG,IAAI,CAAC;gBAElB,GAAG,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,GAAG,IAAI;IAK7C,QAAQ,IAAI,MAAM;IAIlB,MAAM,CAAC,QAAQ,CAAC,IAAI,kBAAmC;IACvD,MAAM,CAAC,QAAQ,CAAC,KAAK,kBAAwC;IAC7D,MAAM,CAAC,QAAQ,CAAC,KAAK,kBAAwC;IAC7D,MAAM,CAAC,QAAQ,CAAC,OAAO,kBAAsC;IAC7D,MAAM,CAAC,QAAQ,CAAC,KAAK,kBAAqC;IAC1D,MAAM,CAAC,QAAQ,CAAC,QAAQ,kBAAmC;IAC3D,MAAM,CAAC,QAAQ,CAAC,OAAO,kBAAoC;IAC3D,MAAM,CAAC,QAAQ,CAAC,MAAM,kBAAkC;CACzD"}