@johnhenry/math 0.0.0

package/dist/Quaternion.js

@@ -0,0 +1,106 @@
+/**
+ * Quaternion — a number of the form `w + x·i + y·j + z·k`, the four-dimensional
+ * division algebra. Extends the idea of {@link ComplexNumber} and is the
+ * standard tool for representing and composing 3D rotations.
+ */
+export class Quaternion {
+    w;
+    x;
+    y;
+    z;
+    static Identity = new Quaternion(1, 0, 0, 0);
+    static Zero = new Quaternion(0, 0, 0, 0);
+    constructor(w = 0, x = 0, y = 0, z = 0) {
+        this.w = w;
+        this.x = x;
+        this.y = y;
+        this.z = z;
+    }
+    add(q) {
+        return new Quaternion(this.w + q.w, this.x + q.x, this.y + q.y, this.z + q.z);
+    }
+    subtract(q) {
+        return new Quaternion(this.w - q.w, this.x - q.x, this.y - q.y, this.z - q.z);
+    }
+    scale(s) {
+        return new Quaternion(this.w * s, this.x * s, this.y * s, this.z * s);
+    }
+    negate() {
+        return new Quaternion(-this.w, -this.x, -this.y, -this.z);
+    }
+    /** The Hamilton product (non-commutative). */
+    multiply(q) {
+        return new Quaternion(this.w * q.w - this.x * q.x - this.y * q.y - this.z * q.z, this.w * q.x + this.x * q.w + this.y * q.z - this.z * q.y, this.w * q.y - this.x * q.z + this.y * q.w + this.z * q.x, this.w * q.z + this.x * q.y - this.y * q.x + this.z * q.w);
+    }
+    conjugate() {
+        return new Quaternion(this.w, -this.x, -this.y, -this.z);
+    }
+    dot(q) {
+        return this.w * q.w + this.x * q.x + this.y * q.y + this.z * q.z;
+    }
+    get norm() {
+        return Math.hypot(this.w, this.x, this.y, this.z);
+    }
+    normalize() {
+        const n = this.norm;
+        if (n === 0)
+            throw new Error("Cannot normalize the zero quaternion.");
+        return this.scale(1 / n);
+    }
+    /** Multiplicative inverse `q⁻¹ = conj(q) / |q|²`. */
+    inverse() {
+        const n2 = this.w * this.w + this.x * this.x + this.y * this.y + this.z * this.z;
+        if (n2 === 0)
+            throw new Error("The zero quaternion has no inverse.");
+        return this.conjugate().scale(1 / n2);
+    }
+    equals(q) {
+        return this.w === q.w && this.x === q.x && this.y === q.y && this.z === q.z;
+    }
+    /** A unit quaternion representing a rotation of `angle` radians about `axis`. */
+    static fromAxisAngle(axis, angle) {
+        const [ax, ay, az] = axis;
+        const len = Math.hypot(ax, ay, az) || 1;
+        const half = angle / 2;
+        const s = Math.sin(half) / len;
+        return new Quaternion(Math.cos(half), ax * s, ay * s, az * s);
+    }
+    /** Rotate a 3D vector by this (unit) quaternion, returning `[x, y, z]`. */
+    rotateVector(v) {
+        const q = this.normalize();
+        const p = new Quaternion(0, v[0], v[1], v[2]);
+        const r = q.multiply(p).multiply(q.conjugate());
+        return [r.x, r.y, r.z];
+    }
+    /** The 3×3 rotation matrix (as `number[][]`) for this quaternion. */
+    toRotationMatrix() {
+        const { w, x, y, z } = this.normalize();
+        return [
+            [1 - 2 * (y * y + z * z), 2 * (x * y - z * w), 2 * (x * z + y * w)],
+            [2 * (x * y + z * w), 1 - 2 * (x * x + z * z), 2 * (y * z - x * w)],
+            [2 * (x * z - y * w), 2 * (y * z + x * w), 1 - 2 * (x * x + y * y)],
+        ];
+    }
+    /** Spherical linear interpolation between two unit quaternions. */
+    static slerp(a, b, t) {
+        const qa = a.normalize();
+        let qb = b.normalize();
+        let cos = qa.dot(qb);
+        if (cos < 0) {
+            qb = qb.scale(-1);
+            cos = -cos;
+        }
+        if (cos > 0.9995)
+            return qa.add(qb.subtract(qa).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 qa.scale(s0).add(qb.scale(s1));
+    }
+    toString() {
+        return `${this.w} + ${this.x}i + ${this.y}j + ${this.z}k`;
+    }
+}
+//# sourceMappingURL=Quaternion.js.map

package/dist/Quaternion.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Quaternion.js","sourceRoot":"","sources":["../src/Quaternion.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AACH,MAAM,OAAO,UAAU;IACZ,CAAC,CAAS;IACV,CAAC,CAAS;IACV,CAAC,CAAS;IACV,CAAC,CAAS;IAEnB,MAAM,CAAU,QAAQ,GAAG,IAAI,UAAU,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACtD,MAAM,CAAU,IAAI,GAAG,IAAI,UAAU,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAElD,YAAY,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC;QACpC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC;QACX,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC;QACX,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC;QACX,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC;IACb,CAAC;IAED,GAAG,CAAC,CAAa;QACf,OAAO,IAAI,UAAU,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAChF,CAAC;IAED,QAAQ,CAAC,CAAa;QACpB,OAAO,IAAI,UAAU,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAChF,CAAC;IAED,KAAK,CAAC,CAAS;QACb,OAAO,IAAI,UAAU,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACxE,CAAC;IAED,MAAM;QACJ,OAAO,IAAI,UAAU,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC5D,CAAC;IAED,8CAA8C;IAC9C,QAAQ,CAAC,CAAa;QACpB,OAAO,IAAI,UAAU,CACnB,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EACzD,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EACzD,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EACzD,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAC1D,CAAC;IACJ,CAAC;IAED,SAAS;QACP,OAAO,IAAI,UAAU,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC3D,CAAC;IAED,GAAG,CAAC,CAAa;QACf,OAAO,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACnE,CAAC;IAED,IAAI,IAAI;QACN,OAAO,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;IACpD,CAAC;IAED,SAAS;QACP,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC;QACpB,IAAI,CAAC,KAAK,CAAC;YAAE,MAAM,IAAI,KAAK,CAAC,uCAAuC,CAAC,CAAC;QACtE,OAAO,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC3B,CAAC;IAED,qDAAqD;IACrD,OAAO;QACL,MAAM,EAAE,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC;QACjF,IAAI,EAAE,KAAK,CAAC;YAAE,MAAM,IAAI,KAAK,CAAC,qCAAqC,CAAC,CAAC;QACrE,OAAO,IAAI,CAAC,SAAS,EAAE,CAAC,KAAK,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;IACxC,CAAC;IAED,MAAM,CAAC,CAAa;QAClB,OAAO,IAAI,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IAC9E,CAAC;IAED,iFAAiF;IACjF,MAAM,CAAC,aAAa,CAAC,IAA8B,EAAE,KAAa;QAChE,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,IAAI,CAAC;QAC1B,MAAM,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,IAAI,CAAC,CAAC;QACxC,MAAM,IAAI,GAAG,KAAK,GAAG,CAAC,CAAC;QACvB,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC;QAC/B,OAAO,IAAI,UAAU,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC;IAChE,CAAC;IAED,2EAA2E;IAC3E,YAAY,CAAC,CAA2B;QACtC,MAAM,CAAC,GAAG,IAAI,CAAC,SAAS,EAAE,CAAC;QAC3B,MAAM,CAAC,GAAG,IAAI,UAAU,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC9C,MAAM,CAAC,GAAG,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,SAAS,EAAE,CAAC,CAAC;QAChD,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACzB,CAAC;IAED,qEAAqE;IACrE,gBAAgB;QACd,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,IAAI,CAAC,SAAS,EAAE,CAAC;QACxC,OAAO;YACL,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;YACnE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;YACnE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;SACpE,CAAC;IACJ,CAAC;IAED,mEAAmE;IACnE,MAAM,CAAC,KAAK,CAAC,CAAa,EAAE,CAAa,EAAE,CAAS;QAClD,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,QAAQ;QACN,OAAO,GAAG,IAAI,CAAC,CAAC,MAAM,IAAI,CAAC,CAAC,OAAO,IAAI,CAAC,CAAC,OAAO,IAAI,CAAC,CAAC,GAAG,CAAC;IAC5D,CAAC"}

package/dist/Rational.d.ts

@@ -0,0 +1,34 @@
+/**
+ * Rational — an exact rational number `numerator / denominator` backed by
+ * `bigint`, always stored in lowest terms with a positive denominator. Gives the
+ * library exact (error-free) arithmetic, e.g. for exact linear algebra or
+ * polynomial GCDs.
+ */
+export declare class Rational {
+    readonly numerator: bigint;
+    readonly denominator: bigint;
+    static readonly Zero: Rational;
+    static readonly One: Rational;
+    constructor(numerator: bigint | number, denominator?: bigint | number);
+    static from(value: bigint | number | Rational): Rational;
+    /** Best rational approximation of a float via continued fractions. */
+    static fromNumber(value: number, maxDenominator?: number): Rational;
+    add(other: Rational): Rational;
+    subtract(other: Rational): Rational;
+    multiply(other: Rational): Rational;
+    divide(other: Rational): Rational;
+    negate(): Rational;
+    reciprocal(): Rational;
+    abs(): Rational;
+    /** Raise to an integer power (negative powers invert). */
+    pow(exponent: number): Rational;
+    equals(other: Rational): boolean;
+    /** -1, 0 or 1 depending on ordering against `other`. */
+    compare(other: Rational): number;
+    lessThan(other: Rational): boolean;
+    get sign(): number;
+    isZero(): boolean;
+    toNumber(): number;
+    toString(): string;
+}
+//# sourceMappingURL=Rational.d.ts.map

package/dist/Rational.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"Rational.d.ts","sourceRoot":"","sources":["../src/Rational.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AAeH,qBAAa,QAAQ;IACnB,QAAQ,CAAC,SAAS,EAAE,MAAM,CAAC;IAC3B,QAAQ,CAAC,WAAW,EAAE,MAAM,CAAC;IAE7B,MAAM,CAAC,QAAQ,CAAC,IAAI,WAAwB;IAC5C,MAAM,CAAC,QAAQ,CAAC,GAAG,WAAwB;gBAE/B,SAAS,EAAE,MAAM,GAAG,MAAM,EAAE,WAAW,GAAE,MAAM,GAAG,MAAW;IAazE,MAAM,CAAC,IAAI,CAAC,KAAK,EAAE,MAAM,GAAG,MAAM,GAAG,QAAQ,GAAG,QAAQ;IAMxD,sEAAsE;IACtE,MAAM,CAAC,UAAU,CAAC,KAAK,EAAE,MAAM,EAAE,cAAc,SAAY,GAAG,QAAQ;IAoBtE,GAAG,CAAC,KAAK,EAAE,QAAQ,GAAG,QAAQ;IAO9B,QAAQ,CAAC,KAAK,EAAE,QAAQ,GAAG,QAAQ;IAOnC,QAAQ,CAAC,KAAK,EAAE,QAAQ,GAAG,QAAQ;IAInC,MAAM,CAAC,KAAK,EAAE,QAAQ,GAAG,QAAQ;IAKjC,MAAM,IAAI,QAAQ;IAIlB,UAAU,IAAI,QAAQ;IAKtB,GAAG,IAAI,QAAQ;IAIf,0DAA0D;IAC1D,GAAG,CAAC,QAAQ,EAAE,MAAM,GAAG,QAAQ;IAc/B,MAAM,CAAC,KAAK,EAAE,QAAQ,GAAG,OAAO;IAIhC,wDAAwD;IACxD,OAAO,CAAC,KAAK,EAAE,QAAQ,GAAG,MAAM;IAKhC,QAAQ,CAAC,KAAK,EAAE,QAAQ,GAAG,OAAO;IAIlC,IAAI,IAAI,IAAI,MAAM,CAEjB;IAED,MAAM,IAAI,OAAO;IAIjB,QAAQ,IAAI,MAAM;IAIlB,QAAQ,IAAI,MAAM;CAGnB"}

package/dist/Rational.js

@@ -0,0 +1,136 @@
+/**
+ * Rational — an exact rational number `numerator / denominator` backed by
+ * `bigint`, always stored in lowest terms with a positive denominator. Gives the
+ * library exact (error-free) arithmetic, e.g. for exact linear algebra or
+ * polynomial GCDs.
+ */
+function bgcd(a, b) {
+    let x = a < 0n ? -a : a;
+    let y = b < 0n ? -b : b;
+    while (y)
+        [x, y] = [y, x % y];
+    return x;
+}
+const toBig = (x) => {
+    if (typeof x === "bigint")
+        return x;
+    if (!Number.isInteger(x))
+        throw new Error(`Rational requires integer inputs; got ${x}`);
+    return BigInt(x);
+};
+export class Rational {
+    numerator;
+    denominator;
+    static Zero = new Rational(0n, 1n);
+    static One = new Rational(1n, 1n);
+    constructor(numerator, denominator = 1n) {
+        let n = toBig(numerator);
+        let d = toBig(denominator);
+        if (d === 0n)
+            throw new Error("Rational denominator must be non-zero.");
+        if (d < 0n) {
+            n = -n;
+            d = -d;
+        }
+        const g = bgcd(n, d) || 1n;
+        this.numerator = n / g;
+        this.denominator = d / g;
+    }
+    static from(value) {
+        if (value instanceof Rational)
+            return value;
+        if (typeof value === "bigint")
+            return new Rational(value, 1n);
+        return Rational.fromNumber(value);
+    }
+    /** Best rational approximation of a float via continued fractions. */
+    static fromNumber(value, maxDenominator = 1_000_000) {
+        if (Number.isInteger(value))
+            return new Rational(BigInt(value), 1n);
+        if (!Number.isFinite(value))
+            throw new Error("Cannot make a Rational from a non-finite number.");
+        let [h0, h1] = [0n, 1n];
+        let [k0, k1] = [1n, 0n];
+        let x = value;
+        const maxD = BigInt(maxDenominator);
+        for (let i = 0; i < 64; i++) {
+            const a = Math.floor(x);
+            const aB = BigInt(a);
+            [h0, h1] = [h1, aB * h1 + h0];
+            [k0, k1] = [k1, aB * k1 + k0];
+            if (k1 > maxD)
+                return new Rational(h0, k0);
+            const frac = x - a;
+            if (frac < 1e-15)
+                break;
+            x = 1 / frac;
+        }
+        return new Rational(h1, k1);
+    }
+    add(other) {
+        return new Rational(this.numerator * other.denominator + other.numerator * this.denominator, this.denominator * other.denominator);
+    }
+    subtract(other) {
+        return new Rational(this.numerator * other.denominator - other.numerator * this.denominator, this.denominator * other.denominator);
+    }
+    multiply(other) {
+        return new Rational(this.numerator * other.numerator, this.denominator * other.denominator);
+    }
+    divide(other) {
+        if (other.numerator === 0n)
+            throw new Error("Division by zero Rational.");
+        return new Rational(this.numerator * other.denominator, this.denominator * other.numerator);
+    }
+    negate() {
+        return new Rational(-this.numerator, this.denominator);
+    }
+    reciprocal() {
+        if (this.numerator === 0n)
+            throw new Error("Zero has no reciprocal.");
+        return new Rational(this.denominator, this.numerator);
+    }
+    abs() {
+        return this.numerator < 0n ? this.negate() : this;
+    }
+    /** Raise to an integer power (negative powers invert). */
+    pow(exponent) {
+        if (!Number.isInteger(exponent))
+            throw new Error("Rational.pow requires an integer exponent.");
+        if (exponent < 0)
+            return this.reciprocal().pow(-exponent);
+        let result = Rational.One;
+        let base = this;
+        let e = exponent;
+        while (e > 0) {
+            if (e & 1)
+                result = result.multiply(base);
+            base = base.multiply(base);
+            e >>= 1;
+        }
+        return result;
+    }
+    equals(other) {
+        return this.numerator === other.numerator && this.denominator === other.denominator;
+    }
+    /** -1, 0 or 1 depending on ordering against `other`. */
+    compare(other) {
+        const diff = this.numerator * other.denominator - other.numerator * this.denominator;
+        return diff < 0n ? -1 : diff > 0n ? 1 : 0;
+    }
+    lessThan(other) {
+        return this.compare(other) < 0;
+    }
+    get sign() {
+        return this.numerator < 0n ? -1 : this.numerator > 0n ? 1 : 0;
+    }
+    isZero() {
+        return this.numerator === 0n;
+    }
+    toNumber() {
+        return Number(this.numerator) / Number(this.denominator);
+    }
+    toString() {
+        return this.denominator === 1n ? `${this.numerator}` : `${this.numerator}/${this.denominator}`;
+    }
+}
+//# sourceMappingURL=Rational.js.map

package/dist/Rational.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Rational.js","sourceRoot":"","sources":["../src/Rational.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AAEH,SAAS,IAAI,CAAC,CAAS,EAAE,CAAS;IAChC,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACxB,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACxB,OAAO,CAAC;QAAE,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;IAC9B,OAAO,CAAC,CAAC;AACX,CAAC;AAED,MAAM,KAAK,GAAG,CAAC,CAAkB,EAAU,EAAE;IAC3C,IAAI,OAAO,CAAC,KAAK,QAAQ;QAAE,OAAO,CAAC,CAAC;IACpC,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,yCAAyC,CAAC,EAAE,CAAC,CAAC;IACxF,OAAO,MAAM,CAAC,CAAC,CAAC,CAAC;AACnB,CAAC,CAAC;AAEF,MAAM,OAAO,QAAQ;IACV,SAAS,CAAS;IAClB,WAAW,CAAS;IAE7B,MAAM,CAAU,IAAI,GAAG,IAAI,QAAQ,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IAC5C,MAAM,CAAU,GAAG,GAAG,IAAI,QAAQ,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IAE3C,YAAY,SAA0B,EAAE,cAA+B,EAAE;QACvE,IAAI,CAAC,GAAG,KAAK,CAAC,SAAS,CAAC,CAAC;QACzB,IAAI,CAAC,GAAG,KAAK,CAAC,WAAW,CAAC,CAAC;QAC3B,IAAI,CAAC,KAAK,EAAE;YAAE,MAAM,IAAI,KAAK,CAAC,wCAAwC,CAAC,CAAC;QACxE,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC;YACX,CAAC,GAAG,CAAC,CAAC,CAAC;YACP,CAAC,GAAG,CAAC,CAAC,CAAC;QACT,CAAC;QACD,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,EAAE,CAAC;QAC3B,IAAI,CAAC,SAAS,GAAG,CAAC,GAAG,CAAC,CAAC;QACvB,IAAI,CAAC,WAAW,GAAG,CAAC,GAAG,CAAC,CAAC;IAC3B,CAAC;IAED,MAAM,CAAC,IAAI,CAAC,KAAiC;QAC3C,IAAI,KAAK,YAAY,QAAQ;YAAE,OAAO,KAAK,CAAC;QAC5C,IAAI,OAAO,KAAK,KAAK,QAAQ;YAAE,OAAO,IAAI,QAAQ,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC;QAC9D,OAAO,QAAQ,CAAC,UAAU,CAAC,KAAK,CAAC,CAAC;IACpC,CAAC;IAED,sEAAsE;IACtE,MAAM,CAAC,UAAU,CAAC,KAAa,EAAE,cAAc,GAAG,SAAS;QACzD,IAAI,MAAM,CAAC,SAAS,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,QAAQ,CAAC,MAAM,CAAC,KAAK,CAAC,EAAE,EAAE,CAAC,CAAC;QACpE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;YAAE,MAAM,IAAI,KAAK,CAAC,kDAAkD,CAAC,CAAC;QACjG,IAAI,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;QACxB,IAAI,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;QACxB,IAAI,CAAC,GAAG,KAAK,CAAC;QACd,MAAM,IAAI,GAAG,MAAM,CAAC,cAAc,CAAC,CAAC;QACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5B,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;YACxB,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;YACrB,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;YAC9B,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;YAC9B,IAAI,EAAE,GAAG,IAAI;gBAAE,OAAO,IAAI,QAAQ,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;YAC3C,MAAM,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC;YACnB,IAAI,IAAI,GAAG,KAAK;gBAAE,MAAM;YACxB,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC;QACf,CAAC;QACD,OAAO,IAAI,QAAQ,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IAC9B,CAAC;IAED,GAAG,CAAC,KAAe;QACjB,OAAO,IAAI,QAAQ,CACjB,IAAI,CAAC,SAAS,GAAG,KAAK,CAAC,WAAW,GAAG,KAAK,CAAC,SAAS,GAAG,IAAI,CAAC,WAAW,EACvE,IAAI,CAAC,WAAW,GAAG,KAAK,CAAC,WAAW,CACrC,CAAC;IACJ,CAAC;IAED,QAAQ,CAAC,KAAe;QACtB,OAAO,IAAI,QAAQ,CACjB,IAAI,CAAC,SAAS,GAAG,KAAK,CAAC,WAAW,GAAG,KAAK,CAAC,SAAS,GAAG,IAAI,CAAC,WAAW,EACvE,IAAI,CAAC,WAAW,GAAG,KAAK,CAAC,WAAW,CACrC,CAAC;IACJ,CAAC;IAED,QAAQ,CAAC,KAAe;QACtB,OAAO,IAAI,QAAQ,CAAC,IAAI,CAAC,SAAS,GAAG,KAAK,CAAC,SAAS,EAAE,IAAI,CAAC,WAAW,GAAG,KAAK,CAAC,WAAW,CAAC,CAAC;IAC9F,CAAC;IAED,MAAM,CAAC,KAAe;QACpB,IAAI,KAAK,CAAC,SAAS,KAAK,EAAE;YAAE,MAAM,IAAI,KAAK,CAAC,4BAA4B,CAAC,CAAC;QAC1E,OAAO,IAAI,QAAQ,CAAC,IAAI,CAAC,SAAS,GAAG,KAAK,CAAC,WAAW,EAAE,IAAI,CAAC,WAAW,GAAG,KAAK,CAAC,SAAS,CAAC,CAAC;IAC9F,CAAC;IAED,MAAM;QACJ,OAAO,IAAI,QAAQ,CAAC,CAAC,IAAI,CAAC,SAAS,EAAE,IAAI,CAAC,WAAW,CAAC,CAAC;IACzD,CAAC;IAED,UAAU;QACR,IAAI,IAAI,CAAC,SAAS,KAAK,EAAE;YAAE,MAAM,IAAI,KAAK,CAAC,yBAAyB,CAAC,CAAC;QACtE,OAAO,IAAI,QAAQ,CAAC,IAAI,CAAC,WAAW,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC;IACxD,CAAC;IAED,GAAG;QACD,OAAO,IAAI,CAAC,SAAS,GAAG,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC;IACpD,CAAC;IAED,0DAA0D;IAC1D,GAAG,CAAC,QAAgB;QAClB,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,QAAQ,CAAC;YAAE,MAAM,IAAI,KAAK,CAAC,4CAA4C,CAAC,CAAC;QAC/F,IAAI,QAAQ,GAAG,CAAC;YAAE,OAAO,IAAI,CAAC,UAAU,EAAE,CAAC,GAAG,CAAC,CAAC,QAAQ,CAAC,CAAC;QAC1D,IAAI,MAAM,GAAG,QAAQ,CAAC,GAAG,CAAC;QAC1B,IAAI,IAAI,GAAa,IAAI,CAAC;QAC1B,IAAI,CAAC,GAAG,QAAQ,CAAC;QACjB,OAAO,CAAC,GAAG,CAAC,EAAE,CAAC;YACb,IAAI,CAAC,GAAG,CAAC;gBAAE,MAAM,GAAG,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;YAC1C,IAAI,GAAG,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;YAC3B,CAAC,KAAK,CAAC,CAAC;QACV,CAAC;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;IAED,MAAM,CAAC,KAAe;QACpB,OAAO,IAAI,CAAC,SAAS,KAAK,KAAK,CAAC,SAAS,IAAI,IAAI,CAAC,WAAW,KAAK,KAAK,CAAC,WAAW,CAAC;IACtF,CAAC;IAED,wDAAwD;IACxD,OAAO,CAAC,KAAe;QACrB,MAAM,IAAI,GAAG,IAAI,CAAC,SAAS,GAAG,KAAK,CAAC,WAAW,GAAG,KAAK,CAAC,SAAS,GAAG,IAAI,CAAC,WAAW,CAAC;QACrF,OAAO,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC5C,CAAC;IAED,QAAQ,CAAC,KAAe;QACtB,OAAO,IAAI,CAAC,OAAO,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;IACjC,CAAC;IAED,IAAI,IAAI;QACN,OAAO,IAAI,CAAC,SAAS,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,SAAS,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAChE,CAAC;IAED,MAAM;QACJ,OAAO,IAAI,CAAC,SAAS,KAAK,EAAE,CAAC;IAC/B,CAAC;IAED,QAAQ;QACN,OAAO,MAAM,CAAC,IAAI,CAAC,SAAS,CAAC,GAAG,MAAM,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC;IAC3D,CAAC;IAED,QAAQ;QACN,OAAO,IAAI,CAAC,WAAW,KAAK,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,SAAS,IAAI,IAAI,CAAC,WAAW,EAAE,CAAC;IACjG,CAAC"}

package/dist/RealMath.d.ts

@@ -0,0 +1,112 @@
+import { Vector } from "./Vector.ts";
+/**
+ * RealMath — arithmetic, trigonometry, numeric calculus and probability over
+ * real numbers, as plain exported functions (no wrapping class — a `static`
+ * method bag buys nothing over top-level functions for pure scalar math).
+ * Ported from Mallory's ActionScript `RealMath`.
+ *
+ * Matrix operations were dropped: `Structure.realField()` already provides
+ * generic linear algebra, and {@link MatrixMath} covers heavier real-only
+ * needs (LU/QR/Cholesky/SVD). Sample statistics (`mean`/`variance`/percentiles/
+ * regression/etc.) moved to {@link Statistics}, which is real-only by nature
+ * (order-dependent stats don't generalize to complex numbers).
+ *
+ * Bugs fixed from the AS3 original:
+ *  - `subtract` returned `a*b` (!).
+ *  - `equal` ignored its tolerance argument (`Math.abs(a-b) <= 0`).
+ *  - `integrateN` used a nonsensical sample point and had an unreachable return.
+ *  - `differentiateN` was a forward difference despite claiming to be symmetric.
+ *  - `crossProduct` read index 3 for the z-component and used truthiness tests.
+ */
+type RVec = Vector<number>;
+type Unary = (x: number) => number;
+export declare const Zero = 0;
+export declare const One = 1;
+export declare const E: number;
+export declare const PI: number;
+export declare const PHI = 1.61803399;
+export declare const PositiveInfinity: number;
+export declare const NegativeInfinity: number;
+export declare function radiansToDegrees(r: number): number;
+export declare function degreesToRadians(d: number): number;
+/** Equality within `distance` (bug fix: AS3 ignored the tolerance). */
+export declare function equal(a: number, b: number, distance?: number): boolean;
+export declare function compare(a: number, b: number): number;
+export declare function lessThan(a: number, b: number): boolean;
+export declare function greaterThan(a: number, b: number): boolean;
+export declare function negative(a?: number): number;
+/** `1/a` (matches `divide`'s directed-infinity convention rather than always returning NaN). */
+export declare function reciprocal(a?: number): number;
+/** Round away tiny floating-point noise (rounds to 10 decimal places). */
+export declare function identity(a?: number): number;
+export declare function add(a?: number, b?: number): number;
+/** Subtraction (bug fix: AS3 returned `a*b`). */
+export declare function subtract(a?: number, b?: number): number;
+export declare function multiply(a?: number, b?: number): number;
+export declare function divide(a?: number, b?: number): number;
+export declare function power(a?: number, exponent?: number): number;
+export declare function square(a?: number): number;
+export declare function squareRoot(a?: number): number;
+export declare function logarithm(a?: number, base?: number): number;
+export declare function sine(a?: number): number;
+export declare function cosine(a?: number): number;
+export declare function tangent(a?: number): number;
+export declare function cosecant(a?: number): number;
+export declare function secant(a?: number): number;
+export declare function cotangent(a?: number): number;
+export declare function hyperbolicSine(a: number): number;
+export declare function hyperbolicCosine(a: number): number;
+export declare function hyperbolicTangent(a: number): number;
+export declare function hyperbolicCosecant(a: number): number;
+export declare function hyperbolicSecant(a: number): number;
+export declare function hyperbolicCotangent(a: number): number;
+export declare function sumSeries(unaryOperation: Unary, start?: number, end?: number): number;
+export declare function productSeries(unaryOperation: Unary, start?: number, end?: number): number;
+/**
+ * Newton–Raphson root finder for `unaryOperation(x) == target`.
+ *
+ * Bug fix: the AS3 version updated `guess -= (f(guess) - target)`, i.e. Newton
+ * with an assumed derivative of 1 — it diverges for almost any nonlinear
+ * function. Here the derivative is estimated numerically ({@link differentiateN})
+ * so the iteration actually converges, and unbounded recursion is replaced by a
+ * bounded loop.
+ */
+export declare function solveN(unaryOperation: Unary, target: number, desiredAccuracy?: number, guess?: number, maxIterations?: number): number;
+/**
+ * Numeric integral via the midpoint rule (bug fix: the AS3 code sampled at
+ * `2*lowerLimit+interval` and had an unreachable return).
+ */
+export declare function integrateN(unaryOperation: Unary, lowerLimit: number, upperLimit: number, interval?: number): number;
+/**
+ * Numeric derivative via the symmetric difference quotient (bug fix: the AS3
+ * code was a forward difference despite the "symmetric" comment).
+ */
+export declare function differentiateN(unaryOperation: Unary, point: number, limit?: number): number;
+export declare function normalDistribution(mean: number, sDev: number): Unary;
+export declare function standardNormalDistribution(x: number): number;
+export declare function normalProbability(x: number, mean: number, sDev: number, negativeInfinityApproximation?: number, integrationInterval?: number): number;
+export declare function standardNormalProbability(z: number, negativeInfinityApproximation?: number, interval?: number): number;
+/** Round to `precision` decimal places (replaces the buggy Flex NumberFormatter). */
+export declare function roundTo(num: number, precision?: number): number;
+/** Uniform random in `[min(r1,r2), max(r1,r2)]`. */
+export declare function random(r1?: number, r2?: number, _inclusive?: boolean): number;
+/** Placeholder for the random.org fetch; always falls back to {@link random}. */
+export declare function randomOrg(lower?: number, upper?: number, inclusive?: boolean, fallback?: boolean): number;
+export declare function jitter(a: number, magnitude?: number): number;
+export declare function addVector(alpha: RVec, beta: RVec): RVec;
+export declare function scaleVector(alpha: RVec, scalar: number): RVec;
+export declare function negativeVector(alpha: RVec): RVec;
+export declare function subtractVector(alpha: RVec, beta: RVec): RVec;
+export declare function dotProduct(alpha: RVec, beta: RVec): number;
+/**
+ * 3D cross product (bug fix: the AS3 code read index 3 for the z-component
+ * and used truthiness tests that dropped legitimate zero components).
+ */
+export declare function crossProduct(alpha: RVec, beta: RVec): RVec;
+export declare function kroneckerProduct(A: RVec, B: RVec): RVec;
+export declare function distanceVector(alpha: RVec, beta: RVec, norm?: number): number;
+export declare function pNorm(alpha: RVec, norm?: number): number;
+export declare function magnitudeVector(alpha: RVec): number;
+export declare function angleBetween(alpha: RVec, beta: RVec): number;
+export {};
+//# sourceMappingURL=RealMath.d.ts.map

package/dist/RealMath.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"RealMath.d.ts","sourceRoot":"","sources":["../src/RealMath.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AAGrC;;;;;;;;;;;;;;;;;;GAkBG;AAEH,KAAK,IAAI,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC;AAC3B,KAAK,KAAK,GAAG,CAAC,CAAC,EAAE,MAAM,KAAK,MAAM,CAAC;AAInC,eAAO,MAAM,IAAI,IAAI,CAAC;AACtB,eAAO,MAAM,GAAG,IAAI,CAAC;AACrB,eAAO,MAAM,CAAC,QAAS,CAAC;AACxB,eAAO,MAAM,EAAE,QAAU,CAAC;AAC1B,eAAO,MAAM,GAAG,aAAa,CAAC;AAC9B,eAAO,MAAM,gBAAgB,QAAW,CAAC;AACzC,eAAO,MAAM,gBAAgB,QAAY,CAAC;AAI1C,wBAAgB,gBAAgB,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAElD;AAED,wBAAgB,gBAAgB,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAElD;AAID,uEAAuE;AACvE,wBAAgB,KAAK,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,QAAQ,SAAI,GAAG,OAAO,CAEjE;AAED,wBAAgB,OAAO,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAIpD;AAED,wBAAgB,QAAQ,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,OAAO,CAEtD;AAED,wBAAgB,WAAW,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,OAAO,CAEzD;AAID,wBAAgB,QAAQ,CAAC,CAAC,SAAI,GAAG,MAAM,CAEtC;AAED,gGAAgG;AAChG,wBAAgB,UAAU,CAAC,CAAC,SAAI,GAAG,MAAM,CAExC;AAED,0EAA0E;AAC1E,wBAAgB,QAAQ,CAAC,CAAC,SAAI,GAAG,MAAM,CAEtC;AAID,wBAAgB,GAAG,CAAC,CAAC,SAAI,EAAE,CAAC,SAAI,GAAG,MAAM,CAExC;AAED,iDAAiD;AACjD,wBAAgB,QAAQ,CAAC,CAAC,SAAI,EAAE,CAAC,SAAI,GAAG,MAAM,CAE7C;AAED,wBAAgB,QAAQ,CAAC,CAAC,SAAI,EAAE,CAAC,SAAI,GAAG,MAAM,CAE7C;AAED,wBAAgB,MAAM,CAAC,CAAC,SAAI,EAAE,CAAC,SAAI,GAAG,MAAM,CAO3C;AAID,wBAAgB,KAAK,CAAC,CAAC,SAAI,EAAE,QAAQ,SAAI,GAAG,MAAM,CAEjD;AAED,wBAAgB,MAAM,CAAC,CAAC,SAAI,GAAG,MAAM,CAEpC;AAED,wBAAgB,UAAU,CAAC,CAAC,SAAI,GAAG,MAAM,CAExC;AAED,wBAAgB,SAAS,CAAC,CAAC,SAAI,EAAE,IAAI,SAAS,GAAG,MAAM,CAEtD;AAID,wBAAgB,IAAI,CAAC,CAAC,SAAI,GAAG,MAAM,CAElC;AACD,wBAAgB,MAAM,CAAC,CAAC,SAAI,GAAG,MAAM,CAEpC;AACD,wBAAgB,OAAO,CAAC,CAAC,SAAI,GAAG,MAAM,CAErC;AACD,wBAAgB,QAAQ,CAAC,CAAC,SAAI,GAAG,MAAM,CAEtC;AACD,wBAAgB,MAAM,CAAC,CAAC,SAAI,GAAG,MAAM,CAEpC;AACD,wBAAgB,SAAS,CAAC,CAAC,SAAI,GAAG,MAAM,CAEvC;AAED,wBAAgB,cAAc,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAEhD;AACD,wBAAgB,gBAAgB,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAElD;AACD,wBAAgB,iBAAiB,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAEnD;AACD,wBAAgB,kBAAkB,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAEpD;AACD,wBAAgB,gBAAgB,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAElD;AACD,wBAAgB,mBAAmB,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAErD;AAID,wBAAgB,SAAS,CAAC,cAAc,EAAE,KAAK,EAAE,KAAK,SAAI,EAAE,GAAG,SAAI,GAAG,MAAM,CAM3E;AAED,wBAAgB,aAAa,CAAC,cAAc,EAAE,KAAK,EAAE,KAAK,SAAI,EAAE,GAAG,SAAI,GAAG,MAAM,CAM/E;AAED;;;;;;;;GAQG;AACH,wBAAgB,MAAM,CACpB,cAAc,EAAE,KAAK,EACrB,MAAM,EAAE,MAAM,EACd,eAAe,SAAO,EACtB,KAAK,SAAI,EACT,aAAa,SAAQ,GACpB,MAAM,CAUR;AAED;;;GAGG;AACH,wBAAgB,UAAU,CAAC,cAAc,EAAE,KAAK,EAAE,UAAU,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,QAAQ,SAAO,GAAG,MAAM,CAOjH;AAED;;;GAGG;AACH,wBAAgB,cAAc,CAAC,cAAc,EAAE,KAAK,EAAE,KAAK,EAAE,MAAM,EAAE,KAAK,SAAO,GAAG,MAAM,CAGzF;AAID,wBAAgB,kBAAkB,CAAC,IAAI,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,GAAG,KAAK,CAEpE;AAED,wBAAgB,0BAA0B,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAE5D;AAED,wBAAgB,iBAAiB,CAC/B,CAAC,EAAE,MAAM,EACT,IAAI,EAAE,MAAM,EACZ,IAAI,EAAE,MAAM,EACZ,6BAA6B,SAAK,EAClC,mBAAmB,SAAS,GAC3B,MAAM,CAER;AAED,wBAAgB,yBAAyB,CAAC,CAAC,EAAE,MAAM,EAAE,6BAA6B,SAAK,EAAE,QAAQ,SAAS,GAAG,MAAM,CAElH;AAID,qFAAqF;AACrF,wBAAgB,OAAO,CAAC,GAAG,EAAE,MAAM,EAAE,SAAS,SAAK,GAAG,MAAM,CAI3D;AAED,oDAAoD;AACpD,wBAAgB,MAAM,CAAC,EAAE,SAAI,EAAE,EAAE,SAAI,EAAE,UAAU,UAAO,GAAG,MAAM,CAIhE;AAED,iFAAiF;AACjF,wBAAgB,SAAS,CAAC,KAAK,SAAI,EAAE,KAAK,SAAI,EAAE,SAAS,UAAO,EAAE,QAAQ,UAAO,GAAG,MAAM,CAGzF;AAED,wBAAgB,MAAM,CAAC,CAAC,EAAE,MAAM,EAAE,SAAS,SAAM,GAAG,MAAM,CAEzD;AAID,wBAAgB,SAAS,CAAC,KAAK,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,GAAG,IAAI,CAEvD;AAED,wBAAgB,WAAW,CAAC,KAAK,EAAE,IAAI,EAAE,MAAM,EAAE,MAAM,GAAG,IAAI,CAE7D;AAED,wBAAgB,cAAc,CAAC,KAAK,EAAE,IAAI,GAAG,IAAI,CAEhD;AAED,wBAAgB,cAAc,CAAC,KAAK,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,GAAG,IAAI,CAE5D;AAED,wBAAgB,UAAU,CAAC,KAAK,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,GAAG,MAAM,CAG1D;AAED;;;GAGG;AACH,wBAAgB,YAAY,CAAC,KAAK,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,GAAG,IAAI,CAW1D;AAED,wBAAgB,gBAAgB,CAAC,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,IAAI,GAAG,IAAI,CAIvD;AAED,wBAAgB,cAAc,CAAC,KAAK,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,SAAI,GAAG,MAAM,CAExE;AAED,wBAAgB,KAAK,CAAC,KAAK,EAAE,IAAI,EAAE,IAAI,SAAI,GAAG,MAAM,CAKnD;AAED,wBAAgB,eAAe,CAAC,KAAK,EAAE,IAAI,GAAG,MAAM,CAEnD;AAED,wBAAgB,YAAY,CAAC,KAAK,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,GAAG,MAAM,CAE5D"}