@hexpro/qr-code-studio 0.0.0-stage → 0.1.0

package/dist/core/math.js

@@ -0,0 +1,609 @@
+/**
+ * Deterministic replacements for the transcendental members of Math.
+ *
+ * ECMAScript specifies `+ - * /` and `Math.sqrt` to IEEE 754 exactly and
+ * leaves sin, cos, atan2, exp, log, pow and hypot "implementation-approximated".
+ * V8, JavaScriptCore and SpiderMonkey really do differ in the last place, and
+ * one disagreement in a curve flattening step moves an anti-aliased edge by a
+ * level, the strict reader thresholds that edge, and the scan check gives a
+ * different verdict in Safari than in Chrome. So everything here is built from
+ * the specified operations alone, plus `Math.abs`, `Math.round` and integer
+ * BigInt arithmetic, which are exact. The results are the same in every
+ * engine, and accurate to better than 1e-12 relative (in practice within a
+ * few units in the last place).
+ *
+ * The algorithms are the classical ones, after fdlibm: Cody-Waite range
+ * reduction with pi/2 split into three 33 bit pieces, exp and log reduced by
+ * ln 2, atan reduced by a table of atan(i/8). The polynomial kernels use
+ * Taylor coefficients, exact rationals, with enough terms that truncation sits
+ * far below rounding; that costs a few multiplications over a minimax fit and
+ * means no coefficient here rests on anybody's memory of a table. The
+ * constants that are not rationals are derived from series in the tests.
+ *
+ * Domains: sin and cos take any finite argument (beyond about 1.6 million the
+ * reduction is done exactly in BigInt against 2/pi to 1216 bits, the way
+ * fdlibm's Payne-Hanek path does it in floating point); atan2, hypot, exp and
+ * ln follow Math for every input, special values included; pow takes a base
+ * of zero or more (see its comment).
+ */
+export const PI = 3.141592653589793;
+export const TAU = 6.283185307179586;
+/* pi/2 in three 33 bit pieces and a tail (fdlibm's pio2_1, pio2_2, pio2_3, pio2_3t). */
+const PIO2_1 = 1.5707963267341256;
+const PIO2_2 = 6.077100506303966e-11;
+const PIO2_3 = 2.0222662487111665e-21;
+const PIO2_3T = 8.4784276603689e-32;
+const INV_PIO2 = 0.6366197723675814;
+const PIO4 = 0.7853981633974483;
+/** Past this, n * PIO2_1 would need more than 53 bits: n < 2^20. */
+const MEDIUM_LIMIT = 1.6e6;
+const PI_LO = 1.2246467991473532e-16;
+const PIO2_HI = 1.5707963267948966;
+const PIO2_LO = 6.123233995736766e-17;
+/** ln 2 as 32 significant bits and the rest, so k * LN2_HI is exact for any exponent k. */
+const LN2_HI = 0.6931471803691238;
+const LN2_LO = 1.9082149292705877e-10;
+const INV_LN2 = 1.4426950408889634;
+const SQRT2 = 1.4142135623730951;
+/** ln of the largest double, and ln of half the smallest subnormal: past these exp is Infinity or 0. */
+const EXP_MAX = 709.782712893384;
+const EXP_MIN = -745.1332191019411;
+const TWO_M27 = 1 / 134217728;
+const TWO_54 = 18014398509481984;
+/** 2^27 + 1, Veltkamp's splitter for exact products without a fused multiply-add. */
+const SPLITTER = 134217729;
+/* Exact helpers. */
+let scratch;
+function view() {
+    scratch ??= new DataView(new ArrayBuffer(8));
+    return scratch;
+}
+/** The unbiased binary exponent of a positive finite x, subnormals included. */
+function exponentOf(x) {
+    const v = view();
+    v.setFloat64(0, x);
+    const biased = (v.getUint32(0) >>> 20) & 0x7ff;
+    if (biased === 0)
+        return exponentOf(x * TWO_54) - 54;
+    return biased - 1023;
+}
+/** 2^k exactly, for an integer k from -1022 to 1023. */
+function powerOfTwo(k) {
+    const v = view();
+    v.setUint32(0, (k + 1023) << 20);
+    v.setUint32(4, 0);
+    return v.getFloat64(0);
+}
+/** x * 2^k, in steps the exponent range allows, so only the last step can round. */
+function scale(x, k) {
+    let value = x;
+    let rest = k;
+    while (rest > 1000) {
+        value *= powerOfTwo(1000);
+        rest -= 1000;
+    }
+    while (rest < -1000) {
+        value *= powerOfTwo(-1000);
+        rest += 1000;
+    }
+    return value * powerOfTwo(rest);
+}
+let sumError = 0;
+/** a + b, leaving the exact rounding error in `sumError` (Knuth's TwoSum). */
+function twoSum(a, b) {
+    const s = a + b;
+    const bb = s - a;
+    sumError = a - (s - bb) + (b - bb);
+    return s;
+}
+/** The exact error of the product p = a * b, for |a|, |b| below 2^996 (Dekker). */
+function productError(a, b, p) {
+    let c = SPLITTER * a;
+    const ah = c - (c - a);
+    const al = a - ah;
+    c = SPLITTER * b;
+    const bh = c - (c - b);
+    const bl = b - bh;
+    return ah * bh - p + ah * bl + al * bh + al * bl;
+}
+/* sin and cos. */
+/* Taylor coefficients of sin: (-1)^k / (2k + 1)!. */
+const S1 = -1 / 6;
+const S2 = 1 / 120;
+const S3 = -1 / 5040;
+const S4 = 1 / 362880;
+const S5 = -1 / 39916800;
+const S6 = 1 / 6227020800;
+const S7 = -1 / 1307674368000;
+const S8 = 1 / 355687428096000;
+/* Taylor coefficients of cos from x^4: (-1)^k / (2k)!. */
+const C1 = 1 / 24;
+const C2 = -1 / 720;
+const C3 = 1 / 40320;
+const C4 = -1 / 3628800;
+const C5 = 1 / 479001600;
+const C6 = -1 / 87178291200;
+const C7 = 1 / 20922789888000;
+const C8 = -1 / 6402373705728000;
+/**
+ * sin(x + y) for |x| <= pi/4 and y a tail below an ulp of x. The series runs
+ * to x^17, whose successor is below 1e-19 of the result.
+ */
+function kernelSin(x, y) {
+    const z = x * x;
+    const v = z * x;
+    const r = S2 + z * (S3 + z * (S4 + z * (S5 + z * (S6 + z * (S7 + z * S8)))));
+    return x - (z * (0.5 * y - v * r) - y - v * S1);
+}
+/** cos(x + y), same conditions. 1 - x^2/2 is formed with its rounding error recovered. */
+function kernelCos(x, y) {
+    const z = x * x;
+    const r = z * (C1 + z * (C2 + z * (C3 + z * (C4 + z * (C5 + z * (C6 + z * (C7 + z * C8)))))));
+    const hz = 0.5 * z;
+    const w = 1 - hz;
+    return w + (1 - w - hz + (z * r - x * y));
+}
+let reducedHi = 0;
+let reducedLo = 0;
+/**
+ * x - n pi/2 for the nearest integer n, as reducedHi + reducedLo, returning n.
+ *
+ * Each piece of pi/2 has 33 bits and n fewer than 21, so every product is
+ * exact, and the subtractions are TwoSums, so the reduced argument carries
+ * about 150 bits of pi/2: enough for the closest any double in range comes to
+ * a multiple of pi/2.
+ */
+function reduceMedium(x) {
+    const n = Math.round(x * INV_PIO2);
+    const r1 = x - n * PIO2_1;
+    const r2 = twoSum(r1, -(n * PIO2_2));
+    let tail = sumError;
+    const r3 = twoSum(r2, -(n * PIO2_3));
+    tail += sumError - n * PIO2_3T;
+    reducedHi = r3 + tail;
+    reducedLo = tail - (reducedHi - r3);
+    return n;
+}
+/** Fraction bits of 2/pi and of pi/2 for the exact reduction. */
+const TWO_OVER_PI_BITS = 1216;
+const PIO2_BITS = 128;
+let exact;
+/**
+ * floor(2/pi * 2^1216) and floor(pi/2 * 2^128), from Machin's formula in
+ * BigInt, on first use. Integer arithmetic throughout, so the constants are
+ * the same everywhere and nothing has to be copied from a table.
+ */
+function exactConstants() {
+    if (exact)
+        return exact;
+    const working = TWO_OVER_PI_BITS + PIO2_BITS + 64;
+    const one = 1n << BigInt(working);
+    const arctanInverse = (n) => {
+        const n2 = n * n;
+        let term = one / n;
+        let sum = term;
+        for (let k = 1n; term !== 0n; k += 1n) {
+            term /= n2;
+            sum += (k & 1n ? -term : term) / (2n * k + 1n);
+        }
+        return sum;
+    };
+    // pi = 16 atan(1/5) - 4 atan(1/239), scaled by 2^working.
+    const pi = 16n * arctanInverse(5n) - 4n * arctanInverse(239n);
+    exact = {
+        twoOverPi: (1n << BigInt(working + TWO_OVER_PI_BITS + 1)) / pi,
+        pio2: pi >> BigInt(working - PIO2_BITS + 1),
+    };
+    return exact;
+}
+function bitLength(value) {
+    const hex = value.toString(16);
+    return (hex.length - 1) * 4 + (32 - Math.clz32(parseInt(hex[0], 16)));
+}
+/**
+ * The same reduction for |x| >= 1.6e6, done exactly: x is an integer times a
+ * power of two, so x * 2/pi mod 4 is a product of integers, and only the
+ * final reduced argument is rounded.
+ */
+function reduceLarge(x) {
+    const { twoOverPi, pio2 } = exactConstants();
+    const ax = Math.abs(x);
+    const e = exponentOf(ax) - 52;
+    const mantissa = BigInt(scale(ax, -e));
+    // x * 2/pi * 2^shift, and its value mod 4.
+    const shift = TWO_OVER_PI_BITS - e;
+    const big = BigInt(shift);
+    const wrapped = (mantissa * twoOverPi) & ((1n << (big + 2n)) - 1n);
+    const nearest = (wrapped + (1n << (big - 1n))) >> big;
+    const fraction = wrapped - (nearest << big);
+    // fraction / 2^shift * pi/2, kept to 110 bits as a double-double.
+    const product = fraction * pio2;
+    const magnitude = product < 0n ? -product : product;
+    const drop = Math.max(0, bitLength(magnitude) - 110);
+    const top = magnitude >> BigInt(drop);
+    const hi = Number(top);
+    const lo = Number(top - BigInt(hi));
+    const exponent = drop - shift - PIO2_BITS;
+    const sign = product < 0n !== x < 0 ? -1 : 1;
+    const yHi = sign * scale(hi, exponent);
+    const yLo = sign * scale(lo, exponent);
+    reducedHi = yHi + yLo;
+    reducedLo = yLo - (reducedHi - yHi);
+    const n = Number(nearest & 3n);
+    return x < 0 ? (4 - n) & 3 : n;
+}
+function reduce(x) {
+    return Math.abs(x) < MEDIUM_LIMIT ? reduceMedium(x) & 3 : reduceLarge(x);
+}
+export function sin(x) {
+    const ax = Math.abs(x);
+    if (!(ax < Infinity))
+        return Number.NaN;
+    if (ax < TWO_M27)
+        return x;
+    if (ax <= PIO4)
+        return kernelSin(x, 0);
+    switch (reduce(x)) {
+        case 0:
+            return kernelSin(reducedHi, reducedLo);
+        case 1:
+            return kernelCos(reducedHi, reducedLo);
+        case 2:
+            return -kernelSin(reducedHi, reducedLo);
+        default:
+            return -kernelCos(reducedHi, reducedLo);
+    }
+}
+export function cos(x) {
+    const ax = Math.abs(x);
+    if (!(ax < Infinity))
+        return Number.NaN;
+    if (ax < TWO_M27)
+        return 1;
+    if (ax <= PIO4)
+        return kernelCos(x, 0);
+    switch (reduce(x)) {
+        case 0:
+            return kernelCos(reducedHi, reducedLo);
+        case 1:
+            return -kernelSin(reducedHi, reducedLo);
+        case 2:
+            return -kernelCos(reducedHi, reducedLo);
+        default:
+            return kernelSin(reducedHi, reducedLo);
+    }
+}
+/* atan2. */
+/** atan(i/8) for i = 0..8, as the nearest double and the rest. */
+const ATAN_HI = [
+    0, 0.12435499454676144, 0.24497866312686414, 0.35877067027057225, 0.4636476090008061,
+    0.5585993153435624, 0.6435011087932844, 0.7188299996216245, 0.7853981633974483,
+];
+const ATAN_LO = [
+    0, -3.1253241424539383e-18, 1.0698755618734451e-17, -2.4623815582638635e-17,
+    2.2698777452961687e-17, -5.4556305485916264e-18, 1.5834785051444286e-17, -2.1478388444456983e-17,
+    3.061616997868383e-17,
+];
+/** atan(u) - u for |u| <= 1/16, by its series to u^21: the next term is below 1e-28 of u. */
+function atanTail(u) {
+    const z = u * u;
+    let p = 1 / 21;
+    p = -1 / 19 + z * p;
+    p = 1 / 17 + z * p;
+    p = -1 / 15 + z * p;
+    p = 1 / 13 + z * p;
+    p = -1 / 11 + z * p;
+    p = 1 / 9 + z * p;
+    p = -1 / 7 + z * p;
+    p = 1 / 5 + z * p;
+    p = -1 / 3 + z * p;
+    return u * z * p;
+}
+/**
+ * atan(t) for t in [0, 1]: t is moved to the nearest c = i/8, and
+ * atan(t) = atan(c) + atan((t - c) / (1 + t c)), where t - c is exact and the
+ * quotient is at most 1/16.
+ */
+function atanUnit(t) {
+    const i = Math.round(t * 8);
+    if (i === 0)
+        return t + atanTail(t);
+    const c = i / 8;
+    const u = (t - c) / (1 + t * c);
+    return ATAN_HI[i] + (ATAN_LO[i] + (u + atanTail(u)));
+}
+/** atan(y / x) in the right quadrant, matching Math.atan2 on every special value, signed zeros included. */
+export function atan2(y, x) {
+    if (y !== y || x !== x)
+        return Number.NaN;
+    const ax = Math.abs(x);
+    const ay = Math.abs(y);
+    const negative = y < 0 || Object.is(y, -0);
+    let angle;
+    if (ay === 0) {
+        // On the x axis: +0 or pi, with the sign of y.
+        angle = x > 0 || Object.is(x, 0) ? 0 : PI;
+    }
+    else if (ax === 0) {
+        angle = PIO2_HI;
+    }
+    else if (ax === Infinity) {
+        angle = ay === Infinity ? (x > 0 ? PIO4 : 3 * PIO4) : x > 0 ? 0 : PI;
+    }
+    else if (ay === Infinity) {
+        angle = PIO2_HI;
+    }
+    else if (ay <= ax) {
+        // Divide the smaller by the larger, so the quotient is in [0, 1], and
+        // finish each quadrant in one rounding: pi - a, pi/2 - a and pi/2 + a all
+        // take a's own value with the tail of the constant, rather than
+        // reflecting an angle that has already been rounded.
+        const a = atanUnit(ay / ax);
+        angle = x > 0 ? a : PI - (a - PI_LO);
+    }
+    else {
+        const a = atanUnit(ax / ay);
+        angle = x > 0 ? PIO2_HI - (a - PIO2_LO) : PIO2_HI + (a + PIO2_LO);
+    }
+    return negative ? -angle : angle;
+}
+/* hypot. */
+const BIG = 2.5822498780869086e120; // 2^400
+const SMALL = 3.8725919148493183e-121; // 2^-400
+const DOWN = 2.409919865102884e-181; // 2^-600
+const UP = 4.149515568880993e180; // 2^600
+/**
+ * sqrt(a^2 + b^2) for a >= b with a in [2^-474, 2^424], where both squares and
+ * their rounding errors are normal numbers. The sum of squares is formed
+ * exactly as a pair, and one Newton step on the square root folds the low
+ * part in, so the result is almost always the correctly rounded one (the
+ * plain formula, and V8's Math.hypot, are off by an ulp about a sixth and a
+ * third of the time respectively).
+ */
+function hypotKernel(a, b) {
+    const aa = a * a;
+    const bb = b * b;
+    const high = twoSum(aa, bb);
+    const low = sumError + productError(a, a, aa) + productError(b, b, bb);
+    const root = Math.sqrt(high);
+    const square = root * root;
+    // high - square is exact: the two are within an ulp of each other.
+    return root + (high - square - productError(root, root, square) + low) / (2 * root);
+}
+/**
+ * sqrt(x^2 + y^2) without overflow or underflow on the way: both are scaled
+ * by a power of two (exactly) when the larger is far from 1. Matches
+ * Math.hypot on special values: Infinity wins over NaN.
+ */
+export function hypot(x, y) {
+    let a = Math.abs(x);
+    let b = Math.abs(y);
+    if (a === Infinity || b === Infinity)
+        return Infinity;
+    if (a !== a || b !== b)
+        return Number.NaN;
+    if (a < b) {
+        const swap = a;
+        a = b;
+        b = swap;
+    }
+    if (b === 0)
+        return a;
+    if (a > BIG)
+        return hypotKernel(a * DOWN, b * DOWN) * UP;
+    if (a < SMALL)
+        return hypotKernel(a * UP, b * UP) * DOWN;
+    return hypotKernel(a, b);
+}
+/* exp, ln and pow. */
+/* 1 / k! for k = 2..15. */
+const E2 = 1 / 2;
+const E3 = 1 / 6;
+const E4 = 1 / 24;
+const E5 = 1 / 120;
+const E6 = 1 / 720;
+const E7 = 1 / 5040;
+const E8 = 1 / 40320;
+const E9 = 1 / 362880;
+const E10 = 1 / 3628800;
+const E11 = 1 / 39916800;
+const E12 = 1 / 479001600;
+const E13 = 1 / 6227020800;
+const E14 = 1 / 87178291200;
+const E15 = 1 / 1307674368000;
+/**
+ * e^(x + tail), x already checked to be in range.
+ *
+ * x = k ln 2 + r with |r| <= ln 2 / 2, r carried with its rounding error, then
+ * e^r by its series to r^15 (the next term is below 3e-21) and a scaling by
+ * 2^k, which is exact unless the result is subnormal.
+ */
+function expKernel(x, tail) {
+    const k = Math.round(x * INV_LN2);
+    const hi = x - k * LN2_HI;
+    const r = twoSum(hi, tail - k * LN2_LO);
+    const c = sumError;
+    let p = E15;
+    p = E14 + r * p;
+    p = E13 + r * p;
+    p = E12 + r * p;
+    p = E11 + r * p;
+    p = E10 + r * p;
+    p = E9 + r * p;
+    p = E8 + r * p;
+    p = E7 + r * p;
+    p = E6 + r * p;
+    p = E5 + r * p;
+    p = E4 + r * p;
+    p = E3 + r * p;
+    p = E2 + r * p;
+    const series = r * r * p;
+    const em1 = r + series;
+    return scale(1 + (em1 + c * (1 + em1)), k);
+}
+export function exp(x) {
+    if (x !== x)
+        return Number.NaN;
+    if (x > EXP_MAX)
+        return Infinity;
+    if (x < EXP_MIN)
+        return 0;
+    return expKernel(x, 0);
+}
+/* 2 / (2j + 1) for j = 1..12: the series of log((1 + s) / (1 - s)) / s - 2 in s^2. */
+const L1 = 2 / 3;
+const L2 = 2 / 5;
+const L3 = 2 / 7;
+const L4 = 2 / 9;
+const L5 = 2 / 11;
+const L6 = 2 / 13;
+const L7 = 2 / 15;
+const L8 = 2 / 17;
+const L9 = 2 / 19;
+const L10 = 2 / 21;
+const L11 = 2 / 23;
+const L12 = 2 / 25;
+let lnLo = 0;
+/**
+ * ln(x) for positive finite x, as a double, with a second-order tail left in
+ * `lnLo` for pow.
+ *
+ * x = 2^k m with m in [sqrt(1/2), sqrt(2)), f = m - 1 (exact), and
+ * ln(1 + f) = f - f^2/2 + s (f^2/2 + R) with s = f / (2 + f), as fdlibm
+ * arranges it so the large terms are exact and only small ones round. R is
+ * the series in s^2 to s^24, past which the terms are below 1e-21.
+ */
+function lnKernel(x) {
+    let k = exponentOf(x);
+    let m = scale(x, -k);
+    if (m > SQRT2) {
+        m *= 0.5;
+        k += 1;
+    }
+    const f = m - 1;
+    const s = f / (2 + f);
+    const z = s * s;
+    let R = L12;
+    R = L11 + z * R;
+    R = L10 + z * R;
+    R = L9 + z * R;
+    R = L8 + z * R;
+    R = L7 + z * R;
+    R = L6 + z * R;
+    R = L5 + z * R;
+    R = L4 + z * R;
+    R = L3 + z * R;
+    R = L2 + z * R;
+    R = L1 + z * R;
+    R *= z;
+    // f^2 / 2 exactly, as a pair, so the one term of order f^2 does not round.
+    const square = f * f;
+    const hfsq = 0.5 * square;
+    const hfsqLo = 0.5 * productError(f, f, square);
+    const small = s * (hfsq + R) + k * LN2_LO - hfsqLo;
+    // k ln2_hi + f - hfsq + small, summed so the rounding errors are kept.
+    const a = twoSum(k * LN2_HI, f);
+    let tail = sumError;
+    const b = twoSum(a, -hfsq);
+    tail += sumError;
+    const c = twoSum(b, small);
+    tail += sumError;
+    const result = c + tail;
+    lnLo = tail - (result - c);
+    return result;
+}
+export function ln(x) {
+    if (x !== x || x < 0)
+        return Number.NaN;
+    if (x === 0)
+        return -Infinity;
+    if (x === Infinity)
+        return Infinity;
+    return lnKernel(x);
+}
+/**
+ * b^n exactly, when every step of binary exponentiation is exact and stays
+ * where Dekker's product check is valid; otherwise undefined. So pow(10, 3) is
+ * 1000 and pow(2, -10) is 1/1024, rather than the nearest thing exp and ln
+ * can reach.
+ */
+function exactPower(base, exponent) {
+    // 2^-969 and 2^996: Dekker's split overflows above the one, and below the
+    // other the error term of a product can underflow and look exact.
+    const inRange = (v) => v >= 2.004168360008973e-292 && v < 6.696928794914171e299;
+    let result = 1;
+    let square = base;
+    for (let n = Math.abs(exponent);;) {
+        if (n & 1) {
+            const product = result * square;
+            if (!inRange(product) || productError(result, square, product) !== 0)
+                return undefined;
+            result = product;
+        }
+        n >>>= 1;
+        if (n === 0)
+            break;
+        const next = square * square;
+        if (!inRange(next) || productError(square, square, next) !== 0)
+            return undefined;
+        square = next;
+    }
+    if (exponent > 0)
+        return result;
+    const inverse = 1 / result;
+    return inRange(inverse) && inverse * result === 1 && productError(inverse, result, 1) === 0
+        ? inverse
+        : undefined;
+}
+/**
+ * base^exponent for base >= 0.
+ *
+ * A negative base throws RangeError: every use in this package raises a length
+ * or a ratio, a negative one is a bug, and a NaN would travel silently into an
+ * export. A zero base, of either sign, is treated as +0. Everything else
+ * matches Math.pow, including pow(1, Infinity) being NaN and pow(NaN, 0) being
+ * 1. Integer exponents up to 64 in size give exact results whenever the exact
+ * result is a double; otherwise the result is exp(exponent ln base) carried in
+ * double-double, which keeps it within about 1e-14 relative even at the edges
+ * of the range.
+ */
+export function pow(base, exponent) {
+    if (exponent === 0)
+        return 1;
+    if (base !== base || exponent !== exponent)
+        return Number.NaN;
+    if (base < 0)
+        throw new RangeError('pow takes a base of zero or more');
+    if (base === 0)
+        return exponent > 0 ? 0 : Infinity;
+    if (exponent === Infinity)
+        return base > 1 ? Infinity : base < 1 ? 0 : Number.NaN;
+    if (exponent === -Infinity)
+        return base > 1 ? 0 : base < 1 ? Infinity : Number.NaN;
+    if (base === Infinity)
+        return exponent > 0 ? Infinity : 0;
+    if (base === 1 || exponent === 1)
+        return base;
+    if (exponent === 0.5)
+        return Math.sqrt(base);
+    if (Number.isInteger(exponent) && Math.abs(exponent) <= 64) {
+        const found = exactPower(base, exponent);
+        if (found !== undefined)
+            return found;
+    }
+    const hi = lnKernel(base);
+    const lo = lnLo;
+    const product = exponent * hi;
+    if (product > EXP_MAX + 1)
+        return Infinity;
+    if (product < EXP_MIN - 1)
+        return 0;
+    const tail = productError(exponent, hi, product) + exponent * lo;
+    const total = product + tail;
+    if (total > EXP_MAX)
+        return Infinity;
+    if (total < EXP_MIN)
+        return 0;
+    return expKernel(product, tail);
+}
+//# sourceMappingURL=math.js.map

package/dist/core/math.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"math.js","sourceRoot":"","sources":["../../src/core/math.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;GA2BG;AAEH,MAAM,CAAC,MAAM,EAAE,GAAG,iBAAiB,CAAC;AACpC,MAAM,CAAC,MAAM,GAAG,GAAG,iBAAiB,CAAC;AAErC,wFAAwF;AACxF,MAAM,MAAM,GAAG,kBAAkB,CAAC;AAClC,MAAM,MAAM,GAAG,qBAAqB,CAAC;AACrC,MAAM,MAAM,GAAG,sBAAsB,CAAC;AACtC,MAAM,OAAO,GAAG,mBAAmB,CAAC;AACpC,MAAM,QAAQ,GAAG,kBAAkB,CAAC;AACpC,MAAM,IAAI,GAAG,kBAAkB,CAAC;AAChC,oEAAoE;AACpE,MAAM,YAAY,GAAG,KAAK,CAAC;AAE3B,MAAM,KAAK,GAAG,sBAAsB,CAAC;AACrC,MAAM,OAAO,GAAG,kBAAkB,CAAC;AACnC,MAAM,OAAO,GAAG,qBAAqB,CAAC;AAEtC,2FAA2F;AAC3F,MAAM,MAAM,GAAG,kBAAkB,CAAC;AAClC,MAAM,MAAM,GAAG,sBAAsB,CAAC;AACtC,MAAM,OAAO,GAAG,kBAAkB,CAAC;AACnC,MAAM,KAAK,GAAG,kBAAkB,CAAC;AACjC,wGAAwG;AACxG,MAAM,OAAO,GAAG,gBAAgB,CAAC;AACjC,MAAM,OAAO,GAAG,CAAC,iBAAiB,CAAC;AAEnC,MAAM,OAAO,GAAG,CAAC,GAAG,SAAS,CAAC;AAC9B,MAAM,MAAM,GAAG,iBAAiB,CAAC;AACjC,qFAAqF;AACrF,MAAM,QAAQ,GAAG,SAAS,CAAC;AAE3B,oBAAoB;AAEpB,IAAI,OAA6B,CAAC;AAElC,SAAS,IAAI;IACZ,OAAO,KAAK,IAAI,QAAQ,CAAC,IAAI,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7C,OAAO,OAAO,CAAC;AAChB,CAAC;AAED,gFAAgF;AAChF,SAAS,UAAU,CAAC,CAAS;IAC5B,MAAM,CAAC,GAAG,IAAI,EAAE,CAAC;IACjB,CAAC,CAAC,UAAU,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACnB,MAAM,MAAM,GAAG,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,GAAG,KAAK,CAAC;IAC/C,IAAI,MAAM,KAAK,CAAC;QAAE,OAAO,UAAU,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,EAAE,CAAC;IACrD,OAAO,MAAM,GAAG,IAAI,CAAC;AACtB,CAAC;AAED,wDAAwD;AACxD,SAAS,UAAU,CAAC,CAAS;IAC5B,MAAM,CAAC,GAAG,IAAI,EAAE,CAAC;IACjB,CAAC,CAAC,SAAS,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,EAAE,CAAC,CAAC;IACjC,CAAC,CAAC,SAAS,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAClB,OAAO,CAAC,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC;AACxB,CAAC;AAED,oFAAoF;AACpF,SAAS,KAAK,CAAC,CAAS,EAAE,CAAS;IAClC,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,OAAO,IAAI,GAAG,IAAI,EAAE,CAAC;QACpB,KAAK,IAAI,UAAU,CAAC,IAAI,CAAC,CAAC;QAC1B,IAAI,IAAI,IAAI,CAAC;IACd,CAAC;IACD,OAAO,IAAI,GAAG,CAAC,IAAI,EAAE,CAAC;QACrB,KAAK,IAAI,UAAU,CAAC,CAAC,IAAI,CAAC,CAAC;QAC3B,IAAI,IAAI,IAAI,CAAC;IACd,CAAC;IACD,OAAO,KAAK,GAAG,UAAU,CAAC,IAAI,CAAC,CAAC;AACjC,CAAC;AAED,IAAI,QAAQ,GAAG,CAAC,CAAC;AAEjB,8EAA8E;AAC9E,SAAS,MAAM,CAAC,CAAS,EAAE,CAAS;IACnC,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAChB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACjB,QAAQ,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;IACnC,OAAO,CAAC,CAAC;AACV,CAAC;AAED,mFAAmF;AACnF,SAAS,YAAY,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS;IACpD,IAAI,CAAC,GAAG,QAAQ,GAAG,CAAC,CAAC;IACrB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACvB,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;IAClB,CAAC,GAAG,QAAQ,GAAG,CAAC,CAAC;IACjB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACvB,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;IAClB,OAAO,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;AAClD,CAAC;AAED,kBAAkB;AAElB,qDAAqD;AACrD,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;AAClB,MAAM,EAAE,GAAG,CAAC,GAAG,GAAG,CAAC;AACnB,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC;AACrB,MAAM,EAAE,GAAG,CAAC,GAAG,MAAM,CAAC;AACtB,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,QAAQ,CAAC;AACzB,MAAM,EAAE,GAAG,CAAC,GAAG,UAAU,CAAC;AAC1B,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,aAAa,CAAC;AAC9B,MAAM,EAAE,GAAG,CAAC,GAAG,eAAe,CAAC;AAC/B,0DAA0D;AAC1D,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;AAClB,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;AACpB,MAAM,EAAE,GAAG,CAAC,GAAG,KAAK,CAAC;AACrB,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC;AACxB,MAAM,EAAE,GAAG,CAAC,GAAG,SAAS,CAAC;AACzB,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,WAAW,CAAC;AAC5B,MAAM,EAAE,GAAG,CAAC,GAAG,cAAc,CAAC;AAC9B,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,gBAAgB,CAAC;AAEjC;;;GAGG;AACH,SAAS,SAAS,CAAC,CAAS,EAAE,CAAS;IACtC,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAChB,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAChB,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7E,OAAO,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;AACjD,CAAC;AAED,0FAA0F;AAC1F,SAAS,SAAS,CAAC,CAAS,EAAE,CAAS;IACtC,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAChB,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC9F,MAAM,EAAE,GAAG,GAAG,GAAG,CAAC,CAAC;IACnB,MAAM,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC;IACjB,OAAO,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AAC3C,CAAC;AAED,IAAI,SAAS,GAAG,CAAC,CAAC;AAClB,IAAI,SAAS,GAAG,CAAC,CAAC;AAElB;;;;;;;GAOG;AACH,SAAS,YAAY,CAAC,CAAS;IAC9B,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC;IACnC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC;IAC1B,MAAM,EAAE,GAAG,MAAM,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;IACrC,IAAI,IAAI,GAAG,QAAQ,CAAC;IACpB,MAAM,EAAE,GAAG,MAAM,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;IACrC,IAAI,IAAI,QAAQ,GAAG,CAAC,GAAG,OAAO,CAAC;IAC/B,SAAS,GAAG,EAAE,GAAG,IAAI,CAAC;IACtB,SAAS,GAAG,IAAI,GAAG,CAAC,SAAS,GAAG,EAAE,CAAC,CAAC;IACpC,OAAO,CAAC,CAAC;AACV,CAAC;AAED,iEAAiE;AACjE,MAAM,gBAAgB,GAAG,IAAI,CAAC;AAC9B,MAAM,SAAS,GAAG,GAAG,CAAC;AACtB,IAAI,KAAwE,CAAC;AAE7E;;;;GAIG;AACH,SAAS,cAAc;IACtB,IAAI,KAAK;QAAE,OAAO,KAAK,CAAC;IACxB,MAAM,OAAO,GAAG,gBAAgB,GAAG,SAAS,GAAG,EAAE,CAAC;IAClD,MAAM,GAAG,GAAG,EAAE,IAAI,MAAM,CAAC,OAAO,CAAC,CAAC;IAClC,MAAM,aAAa,GAAG,CAAC,CAAS,EAAU,EAAE;QAC3C,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QACjB,IAAI,IAAI,GAAG,GAAG,GAAG,CAAC,CAAC;QACnB,IAAI,GAAG,GAAG,IAAI,CAAC;QACf,KAAK,IAAI,CAAC,GAAG,EAAE,EAAE,IAAI,KAAK,EAAE,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC;YACvC,IAAI,IAAI,EAAE,CAAC;YACX,GAAG,IAAI,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;QAChD,CAAC;QACD,OAAO,GAAG,CAAC;IACZ,CAAC,CAAC;IACF,0DAA0D;IAC1D,MAAM,EAAE,GAAG,GAAG,GAAG,aAAa,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,aAAa,CAAC,IAAI,CAAC,CAAC;IAC9D,KAAK,GAAG;QACP,SAAS,EAAE,CAAC,EAAE,IAAI,MAAM,CAAC,OAAO,GAAG,gBAAgB,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE;QAC9D,IAAI,EAAE,EAAE,IAAI,MAAM,CAAC,OAAO,GAAG,SAAS,GAAG,CAAC,CAAC;KAC3C,CAAC;IACF,OAAO,KAAK,CAAC;AACd,CAAC;AAED,SAAS,SAAS,CAAC,KAAa;IAC/B,MAAM,GAAG,GAAG,KAAK,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC;IAC/B,OAAO,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAW,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;AACjF,CAAC;AAED;;;;GAIG;AACH,SAAS,WAAW,CAAC,CAAS;IAC7B,MAAM,EAAE,SAAS,EAAE,IAAI,EAAE,GAAG,cAAc,EAAE,CAAC;IAC7C,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACvB,MAAM,CAAC,GAAG,UAAU,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC;IAC9B,MAAM,QAAQ,GAAG,MAAM,CAAC,KAAK,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACvC,2CAA2C;IAC3C,MAAM,KAAK,GAAG,gBAAgB,GAAG,CAAC,CAAC;IACnC,MAAM,GAAG,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC;IAC1B,MAAM,OAAO,GAAG,CAAC,QAAQ,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;IACnE,MAAM,OAAO,GAAG,CAAC,OAAO,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC;IACtD,MAAM,QAAQ,GAAG,OAAO,GAAG,CAAC,OAAO,IAAI,GAAG,CAAC,CAAC;IAC5C,kEAAkE;IAClE,MAAM,OAAO,GAAG,QAAQ,GAAG,IAAI,CAAC;IAChC,MAAM,SAAS,GAAG,OAAO,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,OAAO,CAAC;IACpD,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,SAAS,CAAC,SAAS,CAAC,GAAG,GAAG,CAAC,CAAC;IACrD,MAAM,GAAG,GAAG,SAAS,IAAI,MAAM,CAAC,IAAI,CAAC,CAAC;IACtC,MAAM,EAAE,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC;IACvB,MAAM,EAAE,GAAG,MAAM,CAAC,GAAG,GAAG,MAAM,CAAC,EAAE,CAAC,CAAC,CAAC;IACpC,MAAM,QAAQ,GAAG,IAAI,GAAG,KAAK,GAAG,SAAS,CAAC;IAC1C,MAAM,IAAI,GAAG,OAAO,GAAG,EAAE,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7C,MAAM,GAAG,GAAG,IAAI,GAAG,KAAK,CAAC,EAAE,EAAE,QAAQ,CAAC,CAAC;IACvC,MAAM,GAAG,GAAG,IAAI,GAAG,KAAK,CAAC,EAAE,EAAE,QAAQ,CAAC,CAAC;IACvC,SAAS,GAAG,GAAG,GAAG,GAAG,CAAC;IACtB,SAAS,GAAG,GAAG,GAAG,CAAC,SAAS,GAAG,GAAG,CAAC,CAAC;IACpC,MAAM,CAAC,GAAG,MAAM,CAAC,OAAO,GAAG,EAAE,CAAC,CAAC;IAC/B,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAChC,CAAC;AAED,SAAS,MAAM,CAAC,CAAS;IACxB,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,YAAY,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC;AAC1E,CAAC;AAED,MAAM,UAAU,GAAG,CAAC,CAAS;IAC5B,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACvB,IAAI,CAAC,CAAC,EAAE,GAAG,QAAQ,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IACxC,IAAI,EAAE,GAAG,OAAO;QAAE,OAAO,CAAC,CAAC;IAC3B,IAAI,EAAE,IAAI,IAAI;QAAE,OAAO,SAAS,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACvC,QAAQ,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC;QACnB,KAAK,CAAC;YACL,OAAO,SAAS,CAAC,SAAS,EAAE,SAAS,CAAC,CAAC;QACxC,KAAK,CAAC;YACL,OAAO,SAAS,CAAC,SAAS,EAAE,SAAS,CAAC,CAAC;QACxC,KAAK,CAAC;YACL,OAAO,CAAC,SAAS,CAAC,SAAS,EAAE,SAAS,CAAC,CAAC;QACzC;YACC,OAAO,CAAC,SAAS,CAAC,SAAS,EAAE,SAAS,CAAC,CAAC;IAC1C,CAAC;AACF,CAAC;AAED,MAAM,UAAU,GAAG,CAAC,CAAS;IAC5B,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACvB,IAAI,CAAC,CAAC,EAAE,GAAG,QAAQ,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IACxC,IAAI,EAAE,GAAG,OAAO;QAAE,OAAO,CAAC,CAAC;IAC3B,IAAI,EAAE,IAAI,IAAI;QAAE,OAAO,SAAS,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACvC,QAAQ,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC;QACnB,KAAK,CAAC;YACL,OAAO,SAAS,CAAC,SAAS,EAAE,SAAS,CAAC,CAAC;QACxC,KAAK,CAAC;YACL,OAAO,CAAC,SAAS,CAAC,SAAS,EAAE,SAAS,CAAC,CAAC;QACzC,KAAK,CAAC;YACL,OAAO,CAAC,SAAS,CAAC,SAAS,EAAE,SAAS,CAAC,CAAC;QACzC;YACC,OAAO,SAAS,CAAC,SAAS,EAAE,SAAS,CAAC,CAAC;IACzC,CAAC;AACF,CAAC;AAED,YAAY;AAEZ,kEAAkE;AAClE,MAAM,OAAO,GAAG;IACf,CAAC,EAAE,mBAAmB,EAAE,mBAAmB,EAAE,mBAAmB,EAAE,kBAAkB;IACpF,kBAAkB,EAAE,kBAAkB,EAAE,kBAAkB,EAAE,kBAAkB;CAC9E,CAAC;AACF,MAAM,OAAO,GAAG;IACf,CAAC,EAAE,CAAC,sBAAsB,EAAE,sBAAsB,EAAE,CAAC,sBAAsB;IAC3E,sBAAsB,EAAE,CAAC,sBAAsB,EAAE,sBAAsB,EAAE,CAAC,sBAAsB;IAChG,qBAAqB;CACrB,CAAC;AAEF,6FAA6F;AAC7F,SAAS,QAAQ,CAAC,CAAS;IAC1B,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAChB,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC;IACf,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACpB,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACnB,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACpB,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACnB,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACpB,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAClB,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACnB,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAClB,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACnB,OAAO,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;AAClB,CAAC;AAED;;;;GAIG;AACH,SAAS,QAAQ,CAAC,CAAS;IAC1B,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC5B,IAAI,CAAC,KAAK,CAAC;QAAE,OAAO,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;IACpC,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAChB,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IAChC,OAAQ,OAAO,CAAC,CAAC,CAAY,GAAG,CAAE,OAAO,CAAC,CAAC,CAAY,GAAG,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAC9E,CAAC;AAED,4GAA4G;AAC5G,MAAM,UAAU,KAAK,CAAC,CAAS,EAAE,CAAS;IACzC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IAC1C,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACvB,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACvB,MAAM,QAAQ,GAAG,CAAC,GAAG,CAAC,IAAI,MAAM,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAC3C,IAAI,KAAa,CAAC;IAClB,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;QACd,+CAA+C;QAC/C,KAAK,GAAG,CAAC,GAAG,CAAC,IAAI,MAAM,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;IAC3C,CAAC;SAAM,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;QACrB,KAAK,GAAG,OAAO,CAAC;IACjB,CAAC;SAAM,IAAI,EAAE,KAAK,QAAQ,EAAE,CAAC;QAC5B,KAAK,GAAG,EAAE,KAAK,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;IACtE,CAAC;SAAM,IAAI,EAAE,KAAK,QAAQ,EAAE,CAAC;QAC5B,KAAK,GAAG,OAAO,CAAC;IACjB,CAAC;SAAM,IAAI,EAAE,IAAI,EAAE,EAAE,CAAC;QACrB,sEAAsE;QACtE,0EAA0E;QAC1E,gEAAgE;QAChE,qDAAqD;QACrD,MAAM,CAAC,GAAG,QAAQ,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;QAC5B,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;IACtC,CAAC;SAAM,CAAC;QACP,MAAM,CAAC,GAAG,QAAQ,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;QAC5B,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,OAAO,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,OAAO,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC;IACnE,CAAC;IACD,OAAO,QAAQ,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC;AAClC,CAAC;AAED,YAAY;AAEZ,MAAM,GAAG,GAAG,sBAAsB,CAAC,CAAC,QAAQ;AAC5C,MAAM,KAAK,GAAG,uBAAuB,CAAC,CAAC,SAAS;AAChD,MAAM,IAAI,GAAG,sBAAsB,CAAC,CAAC,SAAS;AAC9C,MAAM,EAAE,GAAG,qBAAqB,CAAC,CAAC,QAAQ;AAE1C;;;;;;;GAOG;AACH,SAAS,WAAW,CAAC,CAAS,EAAE,CAAS;IACxC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACjB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACjB,MAAM,IAAI,GAAG,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IAC5B,MAAM,GAAG,GAAG,QAAQ,GAAG,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACvE,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IAC7B,MAAM,MAAM,GAAG,IAAI,GAAG,IAAI,CAAC;IAC3B,mEAAmE;IACnE,OAAO,IAAI,GAAG,CAAC,IAAI,GAAG,MAAM,GAAG,YAAY,CAAC,IAAI,EAAE,IAAI,EAAE,MAAM,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC;AACrF,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,KAAK,CAAC,CAAS,EAAE,CAAS;IACzC,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACpB,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACpB,IAAI,CAAC,KAAK,QAAQ,IAAI,CAAC,KAAK,QAAQ;QAAE,OAAO,QAAQ,CAAC;IACtD,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IAC1C,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QACX,MAAM,IAAI,GAAG,CAAC,CAAC;QACf,CAAC,GAAG,CAAC,CAAC;QACN,CAAC,GAAG,IAAI,CAAC;IACV,CAAC;IACD,IAAI,CAAC,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC;IACtB,IAAI,CAAC,GAAG,GAAG;QAAE,OAAO,WAAW,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,EAAE,CAAC;IACzD,IAAI,CAAC,GAAG,KAAK;QAAE,OAAO,WAAW,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC;IACzD,OAAO,WAAW,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;AAC1B,CAAC;AAED,sBAAsB;AAEtB,2BAA2B;AAC3B,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;AACjB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;AACjB,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;AAClB,MAAM,EAAE,GAAG,CAAC,GAAG,GAAG,CAAC;AACnB,MAAM,EAAE,GAAG,CAAC,GAAG,GAAG,CAAC;AACnB,MAAM,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC;AACpB,MAAM,EAAE,GAAG,CAAC,GAAG,KAAK,CAAC;AACrB,MAAM,EAAE,GAAG,CAAC,GAAG,MAAM,CAAC;AACtB,MAAM,GAAG,GAAG,CAAC,GAAG,OAAO,CAAC;AACxB,MAAM,GAAG,GAAG,CAAC,GAAG,QAAQ,CAAC;AACzB,MAAM,GAAG,GAAG,CAAC,GAAG,SAAS,CAAC;AAC1B,MAAM,GAAG,GAAG,CAAC,GAAG,UAAU,CAAC;AAC3B,MAAM,GAAG,GAAG,CAAC,GAAG,WAAW,CAAC;AAC5B,MAAM,GAAG,GAAG,CAAC,GAAG,aAAa,CAAC;AAE9B;;;;;;GAMG;AACH,SAAS,SAAS,CAAC,CAAS,EAAE,IAAY;IACzC,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC;IAClC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC;IAC1B,MAAM,CAAC,GAAG,MAAM,CAAC,EAAE,EAAE,IAAI,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC;IACxC,MAAM,CAAC,GAAG,QAAQ,CAAC;IACnB,IAAI,CAAC,GAAG,GAAG,CAAC;IACZ,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAChB,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAChB,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAChB,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAChB,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAChB,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACf,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACf,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACf,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACf,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACf,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACf,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACf,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACf,MAAM,MAAM,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACzB,MAAM,GAAG,GAAG,CAAC,GAAG,MAAM,CAAC;IACvB,OAAO,KAAK,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;AAC5C,CAAC;AAED,MAAM,UAAU,GAAG,CAAC,CAAS;IAC5B,IAAI,CAAC,KAAK,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IAC/B,IAAI,CAAC,GAAG,OAAO;QAAE,OAAO,QAAQ,CAAC;IACjC,IAAI,CAAC,GAAG,OAAO;QAAE,OAAO,CAAC,CAAC;IAC1B,OAAO,SAAS,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;AACxB,CAAC;AAED,sFAAsF;AACtF,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;AACjB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;AACjB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;AACjB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;AACjB,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;AAClB,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;AAClB,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;AAClB,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;AAClB,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;AAClB,MAAM,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC;AACnB,MAAM,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC;AACnB,MAAM,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC;AAEnB,IAAI,IAAI,GAAG,CAAC,CAAC;AAEb;;;;;;;;GAQG;AACH,SAAS,QAAQ,CAAC,CAAS;IAC1B,IAAI,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;IACtB,IAAI,CAAC,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IACrB,IAAI,CAAC,GAAG,KAAK,EAAE,CAAC;QACf,CAAC,IAAI,GAAG,CAAC;QACT,CAAC,IAAI,CAAC,CAAC;IACR,CAAC;IACD,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAChB,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACtB,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAChB,IAAI,CAAC,GAAG,GAAG,CAAC;IACZ,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAChB,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAChB,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACf,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACf,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACf,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACf,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACf,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACf,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACf,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACf,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACf,CAAC,IAAI,CAAC,CAAC;IACP,2EAA2E;IAC3E,MAAM,MAAM,GAAG,CAAC,GAAG,CAAC,CAAC;IACrB,MAAM,IAAI,GAAG,GAAG,GAAG,MAAM,CAAC;IAC1B,MAAM,MAAM,GAAG,GAAG,GAAG,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;IAChD,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,MAAM,GAAG,MAAM,CAAC;IACnD,uEAAuE;IACvE,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,GAAG,MAAM,EAAE,CAAC,CAAC,CAAC;IAChC,IAAI,IAAI,GAAG,QAAQ,CAAC;IACpB,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC;IAC3B,IAAI,IAAI,QAAQ,CAAC;IACjB,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;IAC3B,IAAI,IAAI,QAAQ,CAAC;IACjB,MAAM,MAAM,GAAG,CAAC,GAAG,IAAI,CAAC;IACxB,IAAI,GAAG,IAAI,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IAC3B,OAAO,MAAM,CAAC;AACf,CAAC;AAED,MAAM,UAAU,EAAE,CAAC,CAAS;IAC3B,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IACxC,IAAI,CAAC,KAAK,CAAC;QAAE,OAAO,CAAC,QAAQ,CAAC;IAC9B,IAAI,CAAC,KAAK,QAAQ;QAAE,OAAO,QAAQ,CAAC;IACpC,OAAO,QAAQ,CAAC,CAAC,CAAC,CAAC;AACpB,CAAC;AAED;;;;;GAKG;AACH,SAAS,UAAU,CAAC,IAAY,EAAE,QAAgB;IACjD,0EAA0E;IAC1E,kEAAkE;IAClE,MAAM,OAAO,GAAG,CAAC,CAAS,EAAW,EAAE,CAAC,CAAC,IAAI,sBAAsB,IAAI,CAAC,GAAG,qBAAqB,CAAC;IACjG,IAAI,MAAM,GAAG,CAAC,CAAC;IACf,IAAI,MAAM,GAAG,IAAI,CAAC;IAClB,KAAK,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,IAAK,CAAC;QACpC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;YACX,MAAM,OAAO,GAAG,MAAM,GAAG,MAAM,CAAC;YAChC,IAAI,CAAC,OAAO,CAAC,OAAO,CAAC,IAAI,YAAY,CAAC,MAAM,EAAE,MAAM,EAAE,OAAO,CAAC,KAAK,CAAC;gBAAE,OAAO,SAAS,CAAC;YACvF,MAAM,GAAG,OAAO,CAAC;QAClB,CAAC;QACD,CAAC,MAAM,CAAC,CAAC;QACT,IAAI,CAAC,KAAK,CAAC;YAAE,MAAM;QACnB,MAAM,IAAI,GAAG,MAAM,GAAG,MAAM,CAAC;QAC7B,IAAI,CAAC,OAAO,CAAC,IAAI,CAAC,IAAI,YAAY,CAAC,MAAM,EAAE,MAAM,EAAE,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,SAAS,CAAC;QACjF,MAAM,GAAG,IAAI,CAAC;IACf,CAAC;IACD,IAAI,QAAQ,GAAG,CAAC;QAAE,OAAO,MAAM,CAAC;IAChC,MAAM,OAAO,GAAG,CAAC,GAAG,MAAM,CAAC;IAC3B,OAAO,OAAO,CAAC,OAAO,CAAC,IAAI,OAAO,GAAG,MAAM,KAAK,CAAC,IAAI,YAAY,CAAC,OAAO,EAAE,MAAM,EAAE,CAAC,CAAC,KAAK,CAAC;QAC1F,CAAC,CAAC,OAAO;QACT,CAAC,CAAC,SAAS,CAAC;AACd,CAAC;AAED;;;;;;;;;;;GAWG;AACH,MAAM,UAAU,GAAG,CAAC,IAAY,EAAE,QAAgB;IACjD,IAAI,QAAQ,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC;IAC7B,IAAI,IAAI,KAAK,IAAI,IAAI,QAAQ,KAAK,QAAQ;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IAC9D,IAAI,IAAI,GAAG,CAAC;QAAE,MAAM,IAAI,UAAU,CAAC,kCAAkC,CAAC,CAAC;IACvE,IAAI,IAAI,KAAK,CAAC;QAAE,OAAO,QAAQ,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC;IACnD,IAAI,QAAQ,KAAK,QAAQ;QAAE,OAAO,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC;IAClF,IAAI,QAAQ,KAAK,CAAC,QAAQ;QAAE,OAAO,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC;IACnF,IAAI,IAAI,KAAK,QAAQ;QAAE,OAAO,QAAQ,GAAG,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1D,IAAI,IAAI,KAAK,CAAC,IAAI,QAAQ,KAAK,CAAC;QAAE,OAAO,IAAI,CAAC;IAC9C,IAAI,QAAQ,KAAK,GAAG;QAAE,OAAO,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IAC7C,IAAI,MAAM,CAAC,SAAS,CAAC,QAAQ,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,IAAI,EAAE,EAAE,CAAC;QAC5D,MAAM,KAAK,GAAG,UAAU,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QACzC,IAAI,KAAK,KAAK,SAAS;YAAE,OAAO,KAAK,CAAC;IACvC,CAAC;IACD,MAAM,EAAE,GAAG,QAAQ,CAAC,IAAI,CAAC,CAAC;IAC1B,MAAM,EAAE,GAAG,IAAI,CAAC;IAChB,MAAM,OAAO,GAAG,QAAQ,GAAG,EAAE,CAAC;IAC9B,IAAI,OAAO,GAAG,OAAO,GAAG,CAAC;QAAE,OAAO,QAAQ,CAAC;IAC3C,IAAI,OAAO,GAAG,OAAO,GAAG,CAAC;QAAE,OAAO,CAAC,CAAC;IACpC,MAAM,IAAI,GAAG,YAAY,CAAC,QAAQ,EAAE,EAAE,EAAE,OAAO,CAAC,GAAG,QAAQ,GAAG,EAAE,CAAC;IACjE,MAAM,KAAK,GAAG,OAAO,GAAG,IAAI,CAAC;IAC7B,IAAI,KAAK,GAAG,OAAO;QAAE,OAAO,QAAQ,CAAC;IACrC,IAAI,KAAK,GAAG,OAAO;QAAE,OAAO,CAAC,CAAC;IAC9B,OAAO,SAAS,CAAC,OAAO,EAAE,IAAI,CAAC,CAAC;AACjC,CAAC"}

package/dist/core/prng.d.ts

@@ -0,0 +1,27 @@
+/**
+ * Seeded randomness, for everything in the determinism zone that needs it:
+ * the scan check's simulated damage, dithering, a random-looking design that
+ * must look the same tomorrow.
+ *
+ * xoshiro128** 1.1 (Blackman and Vigna), seeded by splitmix32. Both are
+ * integer arithmetic over 32 bits with `Math.imul` and shifts, which
+ * ECMAScript fixes exactly, so a seed gives the same stream in every engine.
+ * `Math.random` is banned here for the opposite reason: it is unseeded and
+ * differs everywhere, and a verdict that depends on it cannot be reproduced.
+ *
+ * `fork` gives an independent stream per label, derived from the stream's
+ * starting state rather than its current position. So a job can hand
+ * `fork(1)` to one stage and `fork(2)` to another, and adding a draw to one
+ * stage never shifts the numbers the other sees.
+ */
+export interface Prng {
+    /** A number in [0, 1) with 32 bits of randomness. */
+    next(): number;
+    /** An integer in [0, maxExclusive), unbiased. maxExclusive is 1 to 2^32. */
+    int(maxExclusive: number): number;
+    /** An independent stream for `label` (0 to 2^32 - 1), the same however much of this one is used. */
+    fork(label: number): Prng;
+}
+/** A stream seeded by `seed`, an integer from 0 to 2^32 - 1. */
+export declare function createPrng(seed: number): Prng;
+//# sourceMappingURL=prng.d.ts.map

package/dist/core/prng.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"prng.d.ts","sourceRoot":"","sources":["../../src/core/prng.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,MAAM,WAAW,IAAI;IACpB,qDAAqD;IACrD,IAAI,IAAI,MAAM,CAAC;IACf,4EAA4E;IAC5E,GAAG,CAAC,YAAY,EAAE,MAAM,GAAG,MAAM,CAAC;IAClC,oGAAoG;IACpG,IAAI,CAAC,KAAK,EAAE,MAAM,GAAG,IAAI,CAAC;CAC1B;AAoGD,gEAAgE;AAChE,wBAAgB,UAAU,CAAC,IAAI,EAAE,MAAM,GAAG,IAAI,CAG7C"}