@ifc-lite/renderer 6.0.1 → 6.3.1

package/dist/deviation/deviation-statistics-kernel.d.ts

@@ -0,0 +1,111 @@
+/**
+ * Linear passes behind the deviation summary statistics (#6872).
+ *
+ * Every figure is built from at most two O(n) passes over the readback, each
+ * of which can run in one call or in time-boxed slices that yield to the
+ * event loop, so a 25M-point cloud never blocks the main thread in one task.
+ *
+ * Exact percentiles without a copy: for a non-negative float32, the IEEE bit
+ * pattern read as an unsigned integer is monotone in the value. So |d|'s
+ * order statistics are the order statistics of `bits & 0x7fffffff`, found by
+ * a two-level radix select: pass 1 counts the high 16 bits (65,536 bins), the
+ * cumulative counts name the bucket holding each wanted rank, and pass 2
+ * counts the low 16 bits of the values in just those (at most three)
+ * buckets. Memory is four 256 KiB count tables, independent of n, and the
+ * input is never reordered or copied.
+ */
+import type { DeviationHistogram, DeviationHistogramRange } from './deviation-statistics.js';
+/** A pass over `[start, end)`, then an optional step once the pass is done. */
+export interface Phase {
+    scan(start: number, end: number): void;
+    done?(): void;
+}
+export declare function runPhases(length: number, phases: readonly Phase[]): void;
+/**
+ * Run the phases in slices of at most ~{@link SLICE_BUDGET_MS} ms. Always
+ * yields once before starting, so the caller's task (a click, a render)
+ * finishes first and an abort raised in it is seen before any work.
+ */
+export declare function runPhasesSliced(length: number, phases: readonly Phase[], signal?: AbortSignal): Promise<void>;
+export declare function checkMask(values: Float32Array, valid: ArrayLike<number> | undefined): void;
+export declare function checkTolerance(tolerance: number): void;
+/** Valid points with |d| ≤ tolerance. */
+export declare class ToleranceCount implements Phase {
+    private readonly values;
+    private readonly tolerance;
+    private readonly valid?;
+    count: number;
+    constructor(values: Float32Array, tolerance: number, valid?: ArrayLike<number> | undefined);
+    scan(start: number, end: number): void;
+}
+/** Fixed bins over [center − h, center + h], the colour ramp's own range. */
+export declare class HistogramBinner {
+    readonly histogram: DeviationHistogram;
+    private readonly last;
+    constructor(range: DeviationHistogramRange);
+    /** `v` must be finite. */
+    add(v: number): void;
+}
+/** Histogram as a standalone pass. */
+export declare class HistogramPass implements Phase {
+    private readonly values;
+    private readonly valid?;
+    readonly binner: HistogramBinner;
+    constructor(values: Float32Array, range: DeviationHistogramRange, valid?: ArrayLike<number> | undefined);
+    scan(start: number, end: number): void;
+}
+/** Neumaier compensated sum: exact-to-rounding for mixed magnitudes. */
+declare class CompensatedSum {
+    private sum;
+    private compensation;
+    add(x: number): void;
+    value(): number;
+}
+/** Nearest-rank index (0-based) of the `percent`-th percentile of n values. */
+export declare function nearestRank(percent: number, n: number): number;
+/** Percentiles of |d| reported, in percent. */
+export declare const PERCENTILES: readonly [50, 95, 99];
+export interface MomentOptions {
+    valid?: ArrayLike<number>;
+    tolerance?: number;
+    /** Compared against `Math.fround(clipRange)`: the shader pegs in float32. */
+    clipRange?: number;
+    histogram?: DeviationHistogramRange;
+}
+/**
+ * The two passes of a full summary: moments, extremes, counts and the
+ * high-bits table first; then σ about the now-known mean, and the low bits
+ * of the percentile buckets.
+ */
+export declare class StatisticsPasses {
+    private readonly values;
+    private readonly options;
+    readonly phases: readonly Phase[];
+    n: number;
+    within: number;
+    clipped: number;
+    min: number;
+    max: number;
+    maxAbs: number;
+    mean: number;
+    readonly sum: CompensatedSum;
+    readonly sumAbs: CompensatedSum;
+    readonly sumSq: CompensatedSum;
+    readonly residual: CompensatedSum;
+    readonly binner: HistogramBinner | null;
+    /** |d| for each of {@link PERCENTILES}, set after the second pass. */
+    readonly percentiles: number[];
+    private readonly bits;
+    private readonly high;
+    private readonly clipAt;
+    /** Per wanted rank: its high-bits bucket and the rank within that bucket. */
+    private targets;
+    constructor(values: Float32Array, options: MomentOptions);
+    private scanMoments;
+    /** Locate each wanted rank's high-bits bucket from the cumulative counts. */
+    private plan;
+    private scanResiduals;
+    private resolve;
+}
+export {};
+//# sourceMappingURL=deviation-statistics-kernel.d.ts.map

package/dist/deviation/deviation-statistics-kernel.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"deviation-statistics-kernel.d.ts","sourceRoot":"","sources":["../../src/deviation/deviation-statistics-kernel.ts"],"names":[],"mappings":"AAIA;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,KAAK,EAAE,kBAAkB,EAAE,uBAAuB,EAAE,MAAM,2BAA2B,CAAC;AAY7F,+EAA+E;AAC/E,MAAM,WAAW,KAAK;IAClB,IAAI,CAAC,KAAK,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,IAAI,CAAC;IACvC,IAAI,CAAC,IAAI,IAAI,CAAC;CACjB;AAED,wBAAgB,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,SAAS,KAAK,EAAE,GAAG,IAAI,CAKxE;AAqBD;;;;GAIG;AACH,wBAAsB,eAAe,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,SAAS,KAAK,EAAE,EAAE,MAAM,CAAC,EAAE,WAAW,GAAG,OAAO,CAAC,IAAI,CAAC,CAiBnH;AAED,wBAAgB,SAAS,CAAC,MAAM,EAAE,YAAY,EAAE,KAAK,EAAE,SAAS,CAAC,MAAM,CAAC,GAAG,SAAS,GAAG,IAAI,CAI1F;AAED,wBAAgB,cAAc,CAAC,SAAS,EAAE,MAAM,GAAG,IAAI,CAItD;AAaD,yCAAyC;AACzC,qBAAa,cAAe,YAAW,KAAK;IAGpC,OAAO,CAAC,QAAQ,CAAC,MAAM;IACvB,OAAO,CAAC,QAAQ,CAAC,SAAS;IAC1B,OAAO,CAAC,QAAQ,CAAC,KAAK,CAAC;IAJ3B,KAAK,SAAK;gBAEW,MAAM,EAAE,YAAY,EACpB,SAAS,EAAE,MAAM,EACjB,KAAK,CAAC,EAAE,SAAS,CAAC,MAAM,CAAC,YAAA;IAK9C,IAAI,CAAC,KAAK,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,IAAI;CASzC;AAED,6EAA6E;AAC7E,qBAAa,eAAe;IACxB,QAAQ,CAAC,SAAS,EAAE,kBAAkB,CAAC;IACvC,OAAO,CAAC,QAAQ,CAAC,IAAI,CAAS;gBAClB,KAAK,EAAE,uBAAuB;IAc1C,0BAA0B;IAC1B,GAAG,CAAC,CAAC,EAAE,MAAM,GAAG,IAAI;CAMvB;AAED,sCAAsC;AACtC,qBAAa,aAAc,YAAW,KAAK;IAE3B,OAAO,CAAC,QAAQ,CAAC,MAAM;IAAgD,OAAO,CAAC,QAAQ,CAAC,KAAK,CAAC;IAD1G,QAAQ,CAAC,MAAM,EAAE,eAAe,CAAC;gBACJ,MAAM,EAAE,YAAY,EAAE,KAAK,EAAE,uBAAuB,EAAmB,KAAK,CAAC,EAAE,SAAS,CAAC,MAAM,CAAC,YAAA;IAI7H,IAAI,CAAC,KAAK,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,IAAI;CAOzC;AAED,wEAAwE;AACxE,cAAM,cAAc;IAChB,OAAO,CAAC,GAAG,CAAK;IAChB,OAAO,CAAC,YAAY,CAAK;IACzB,GAAG,CAAC,CAAC,EAAE,MAAM,GAAG,IAAI;IAMpB,KAAK,IAAI,MAAM;CAGlB;AAED,+EAA+E;AAC/E,wBAAgB,WAAW,CAAC,OAAO,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAE9D;AAED,+CAA+C;AAC/C,eAAO,MAAM,WAAW,uBAAwB,CAAC;AAEjD,MAAM,WAAW,aAAa;IAC1B,KAAK,CAAC,EAAE,SAAS,CAAC,MAAM,CAAC,CAAC;IAC1B,SAAS,CAAC,EAAE,MAAM,CAAC;IACnB,6EAA6E;IAC7E,SAAS,CAAC,EAAE,MAAM,CAAC;IACnB,SAAS,CAAC,EAAE,uBAAuB,CAAC;CACvC;AAED;;;;GAIG;AACH,qBAAa,gBAAgB;IAsBb,OAAO,CAAC,QAAQ,CAAC,MAAM;IAAgB,OAAO,CAAC,QAAQ,CAAC,OAAO;IArB3E,QAAQ,CAAC,MAAM,EAAE,SAAS,KAAK,EAAE,CAAC;IAClC,CAAC,SAAK;IACN,MAAM,SAAK;IACX,OAAO,SAAK;IACZ,GAAG,SAAY;IACf,GAAG,SAAa;IAChB,MAAM,SAAK;IACX,IAAI,SAAK;IACT,QAAQ,CAAC,GAAG,iBAAwB;IACpC,QAAQ,CAAC,MAAM,iBAAwB;IACvC,QAAQ,CAAC,KAAK,iBAAwB;IACtC,QAAQ,CAAC,QAAQ,iBAAwB;IACzC,QAAQ,CAAC,MAAM,EAAE,eAAe,GAAG,IAAI,CAAC;IACxC,sEAAsE;IACtE,QAAQ,CAAC,WAAW,EAAE,MAAM,EAAE,CAAM;IACpC,OAAO,CAAC,QAAQ,CAAC,IAAI,CAAc;IACnC,OAAO,CAAC,QAAQ,CAAC,IAAI,CAA0B;IAC/C,OAAO,CAAC,QAAQ,CAAC,MAAM,CAAqB;IAC5C,6EAA6E;IAC7E,OAAO,CAAC,OAAO,CAAiE;gBAEnD,MAAM,EAAE,YAAY,EAAmB,OAAO,EAAE,aAAa;IAY1F,OAAO,CAAC,WAAW;IA2BnB,6EAA6E;IAC7E,OAAO,CAAC,IAAI;IAiBZ,OAAO,CAAC,aAAa;IAqBrB,OAAO,CAAC,OAAO;CAQlB"}

package/dist/deviation/deviation-statistics-kernel.js

@@ -0,0 +1,303 @@
+/* This Source Code Form is subject to the terms of the Mozilla Public
+ * License, v. 2.0. If a copy of the MPL was not distributed with this
+ * file, You can obtain one at https://mozilla.org/MPL/2.0/. */
+const ABS_MASK = 0x7fffffff;
+/** |d| bit patterns at or above this are ±Infinity or NaN. */
+const NON_FINITE = 0x7f800000;
+const RADIX = 1 << 16;
+const LOW_MASK = RADIX - 1;
+/** Elements per slice between clock checks. */
+const SLICE = 1 << 16;
+/** Longest stretch of work before yielding, ms. */
+const SLICE_BUDGET_MS = 8;
+export function runPhases(length, phases) {
+    for (const phase of phases) {
+        phase.scan(0, length);
+        phase.done?.();
+    }
+}
+/**
+ * Yield through a MessageChannel task, not `scheduler.yield()`: that
+ * continuation runs ahead of already-queued normal tasks, React renders
+ * included, so a loop yielding through it still starves the UI (the same
+ * reasoning as the IDS validator's yield).
+ */
+function yieldToEventLoop() {
+    return new Promise((resolve) => {
+        const channel = new MessageChannel();
+        channel.port1.onmessage = () => {
+            // An open port keeps a Node process alive.
+            channel.port1.close();
+            channel.port2.close();
+            resolve();
+        };
+        channel.port2.postMessage(null);
+    });
+}
+/**
+ * Run the phases in slices of at most ~{@link SLICE_BUDGET_MS} ms. Always
+ * yields once before starting, so the caller's task (a click, a render)
+ * finishes first and an abort raised in it is seen before any work.
+ */
+export async function runPhasesSliced(length, phases, signal) {
+    await yieldToEventLoop();
+    signal?.throwIfAborted();
+    let sliceStart = performance.now();
+    for (const phase of phases) {
+        for (let start = 0; start < length;) {
+            const end = Math.min(length, start + SLICE);
+            phase.scan(start, end);
+            start = end;
+            if (performance.now() - sliceStart >= SLICE_BUDGET_MS) {
+                await yieldToEventLoop();
+                signal?.throwIfAborted();
+                sliceStart = performance.now();
+            }
+        }
+        phase.done?.();
+    }
+}
+export function checkMask(values, valid) {
+    if (valid && valid.length !== values.length) {
+        throw new RangeError(`deviation statistics: mask length ${valid.length} != values length ${values.length}`);
+    }
+}
+export function checkTolerance(tolerance) {
+    if (!(tolerance >= 0) || !Number.isFinite(tolerance)) {
+        throw new RangeError(`deviation statistics: tolerance must be a finite value ≥ 0, got ${tolerance}`);
+    }
+}
+function bitsOf(values) {
+    return new Uint32Array(values.buffer, values.byteOffset, values.length);
+}
+const bitCast = new Uint32Array(1);
+const floatCast = new Float32Array(bitCast.buffer);
+function floatFromBits(bits) {
+    bitCast[0] = bits;
+    return floatCast[0];
+}
+/** Valid points with |d| ≤ tolerance. */
+export class ToleranceCount {
+    values;
+    tolerance;
+    valid;
+    count = 0;
+    constructor(values, tolerance, valid) {
+        this.values = values;
+        this.tolerance = tolerance;
+        this.valid = valid;
+        checkMask(values, valid);
+        checkTolerance(tolerance);
+    }
+    scan(start, end) {
+        const { values, tolerance, valid } = this;
+        let count = 0;
+        for (let i = start; i < end; i++) {
+            // NaN and ±Infinity fail the comparison on their own.
+            if (Math.abs(values[i]) <= tolerance && (!valid || valid[i] !== 0))
+                count++;
+        }
+        this.count += count;
+    }
+}
+/** Fixed bins over [center − h, center + h], the colour ramp's own range. */
+export class HistogramBinner {
+    histogram;
+    last;
+    constructor(range) {
+        const { center, halfRange, bins } = range;
+        if (!(halfRange > 0) || !Number.isFinite(halfRange) || !Number.isFinite(center)) {
+            throw new RangeError(`deviation histogram: halfRange must be finite and > 0, got ${halfRange}`);
+        }
+        if (!Number.isInteger(bins) || bins < 1) {
+            throw new RangeError(`deviation histogram: bins must be a positive integer, got ${bins}`);
+        }
+        this.histogram = {
+            min: center - halfRange, max: center + halfRange, binWidth: (2 * halfRange) / bins,
+            counts: new Array(bins).fill(0), below: 0, above: 0,
+        };
+        this.last = bins - 1;
+    }
+    /** `v` must be finite. */
+    add(v) {
+        const h = this.histogram;
+        if (v < h.min)
+            h.below++;
+        else if (v > h.max)
+            h.above++;
+        else
+            h.counts[Math.min(this.last, Math.floor((v - h.min) / h.binWidth))]++;
+    }
+}
+/** Histogram as a standalone pass. */
+export class HistogramPass {
+    values;
+    valid;
+    binner;
+    constructor(values, range, valid) {
+        this.values = values;
+        this.valid = valid;
+        checkMask(values, valid);
+        this.binner = new HistogramBinner(range);
+    }
+    scan(start, end) {
+        const { values, valid, binner } = this;
+        for (let i = start; i < end; i++) {
+            const v = values[i];
+            if (Number.isFinite(v) && (!valid || valid[i] !== 0))
+                binner.add(v);
+        }
+    }
+}
+/** Neumaier compensated sum: exact-to-rounding for mixed magnitudes. */
+class CompensatedSum {
+    sum = 0;
+    compensation = 0;
+    add(x) {
+        const t = this.sum + x;
+        if (Math.abs(this.sum) >= Math.abs(x))
+            this.compensation += (this.sum - t) + x;
+        else
+            this.compensation += (x - t) + this.sum;
+        this.sum = t;
+    }
+    value() {
+        return this.sum + this.compensation;
+    }
+}
+/** Nearest-rank index (0-based) of the `percent`-th percentile of n values. */
+export function nearestRank(percent, n) {
+    return Math.min(n - 1, Math.max(0, Math.ceil((percent * n) / 100) - 1));
+}
+/** Percentiles of |d| reported, in percent. */
+export const PERCENTILES = [50, 95, 99];
+/**
+ * The two passes of a full summary: moments, extremes, counts and the
+ * high-bits table first; then σ about the now-known mean, and the low bits
+ * of the percentile buckets.
+ */
+export class StatisticsPasses {
+    values;
+    options;
+    phases;
+    n = 0;
+    within = 0;
+    clipped = 0;
+    min = Infinity;
+    max = -Infinity;
+    maxAbs = 0;
+    mean = 0;
+    sum = new CompensatedSum();
+    sumAbs = new CompensatedSum();
+    sumSq = new CompensatedSum();
+    residual = new CompensatedSum();
+    binner;
+    /** |d| for each of {@link PERCENTILES}, set after the second pass. */
+    percentiles = [];
+    bits;
+    high = new Uint32Array(RADIX);
+    clipAt;
+    /** Per wanted rank: its high-bits bucket and the rank within that bucket. */
+    targets = [];
+    constructor(values, options) {
+        this.values = values;
+        this.options = options;
+        checkMask(values, options.valid);
+        if (options.tolerance !== undefined)
+            checkTolerance(options.tolerance);
+        this.binner = options.histogram ? new HistogramBinner(options.histogram) : null;
+        this.clipAt = options.clipRange === undefined ? undefined : Math.fround(options.clipRange);
+        this.bits = bitsOf(values);
+        this.phases = [
+            { scan: (s, e) => this.scanMoments(s, e), done: () => this.plan() },
+            { scan: (s, e) => this.scanResiduals(s, e), done: () => this.resolve() },
+        ];
+    }
+    scanMoments(start, end) {
+        const { values, bits, high, binner, sum, sumAbs, sumSq, clipAt } = this;
+        const { valid } = this.options;
+        // Absent options compare against NaN / +Infinity, so they never count.
+        const tolerance = this.options.tolerance ?? Number.NaN;
+        const clip = clipAt ?? Number.POSITIVE_INFINITY;
+        let { n, min, max, maxAbs, within, clipped } = this;
+        for (let i = start; i < end; i++) {
+            const b = bits[i] & ABS_MASK;
+            if (b >= NON_FINITE || (valid && valid[i] === 0))
+                continue;
+            const v = values[i];
+            const abs = Math.abs(v);
+            n++;
+            high[b >>> 16]++;
+            sum.add(v);
+            sumAbs.add(abs);
+            sumSq.add(v * v);
+            if (v < min)
+                min = v;
+            if (v > max)
+                max = v;
+            if (abs > maxAbs)
+                maxAbs = abs;
+            if (abs <= tolerance)
+                within++;
+            if (abs >= clip)
+                clipped++;
+            binner?.add(v);
+        }
+        Object.assign(this, { n, min, max, maxAbs, within, clipped });
+    }
+    /** Locate each wanted rank's high-bits bucket from the cumulative counts. */
+    plan() {
+        const n = this.n;
+        if (n === 0)
+            return;
+        this.mean = this.sum.value() / n;
+        const lows = new Map();
+        let bucket = 0;
+        let below = 0;
+        // Ranks ascend with the percentiles, so one forward walk serves all.
+        for (const percent of PERCENTILES) {
+            const rank = nearestRank(percent, n);
+            while (below + this.high[bucket] <= rank)
+                below += this.high[bucket++];
+            let low = lows.get(bucket);
+            if (!low)
+                lows.set(bucket, low = new Uint32Array(RADIX));
+            this.targets.push({ bucket, rank: rank - below, low });
+        }
+    }
+    scanResiduals(start, end) {
+        if (this.n === 0)
+            return;
+        const { values, bits, residual, mean } = this;
+        const { valid } = this.options;
+        // At most three distinct buckets; unused slots never match (> 0xffff).
+        const [t0, t1, t2] = [0, 1, 2].map((k) => this.targets[k]);
+        const b0 = t0.bucket, b1 = t1.bucket === b0 ? -1 : t1.bucket;
+        const b2 = t2.bucket === b0 || t2.bucket === t1.bucket ? -1 : t2.bucket;
+        const l0 = t0.low, l1 = t1.low, l2 = t2.low;
+        for (let i = start; i < end; i++) {
+            const b = bits[i] & ABS_MASK;
+            if (b >= NON_FINITE || (valid && valid[i] === 0))
+                continue;
+            const d = values[i] - mean;
+            residual.add(d * d);
+            const h = b >>> 16;
+            if (h === b0)
+                l0[b & LOW_MASK]++;
+            else if (h === b1)
+                l1[b & LOW_MASK]++;
+            else if (h === b2)
+                l2[b & LOW_MASK]++;
+        }
+    }
+    resolve() {
+        for (const { bucket, rank, low } of this.targets) {
+            let below = 0;
+            let l = 0;
+            while (below + low[l] <= rank)
+                below += low[l++];
+            this.percentiles.push(floatFromBits((bucket << 16) | l));
+        }
+    }
+}
+//# sourceMappingURL=deviation-statistics-kernel.js.map

package/dist/deviation/deviation-statistics-kernel.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"deviation-statistics-kernel.js","sourceRoot":"","sources":["../../src/deviation/deviation-statistics-kernel.ts"],"names":[],"mappings":"AAAA;;+DAE+D;AAqB/D,MAAM,QAAQ,GAAG,UAAU,CAAC;AAC5B,8DAA8D;AAC9D,MAAM,UAAU,GAAG,UAAU,CAAC;AAC9B,MAAM,KAAK,GAAG,CAAC,IAAI,EAAE,CAAC;AACtB,MAAM,QAAQ,GAAG,KAAK,GAAG,CAAC,CAAC;AAC3B,+CAA+C;AAC/C,MAAM,KAAK,GAAG,CAAC,IAAI,EAAE,CAAC;AACtB,mDAAmD;AACnD,MAAM,eAAe,GAAG,CAAC,CAAC;AAQ1B,MAAM,UAAU,SAAS,CAAC,MAAc,EAAE,MAAwB;IAC9D,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;QACzB,KAAK,CAAC,IAAI,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;QACtB,KAAK,CAAC,IAAI,EAAE,EAAE,CAAC;IACnB,CAAC;AACL,CAAC;AAED;;;;;GAKG;AACH,SAAS,gBAAgB;IACrB,OAAO,IAAI,OAAO,CAAO,CAAC,OAAO,EAAE,EAAE;QACjC,MAAM,OAAO,GAAG,IAAI,cAAc,EAAE,CAAC;QACrC,OAAO,CAAC,KAAK,CAAC,SAAS,GAAG,GAAG,EAAE;YAC3B,2CAA2C;YAC3C,OAAO,CAAC,KAAK,CAAC,KAAK,EAAE,CAAC;YACtB,OAAO,CAAC,KAAK,CAAC,KAAK,EAAE,CAAC;YACtB,OAAO,EAAE,CAAC;QACd,CAAC,CAAC;QACF,OAAO,CAAC,KAAK,CAAC,WAAW,CAAC,IAAI,CAAC,CAAC;IACpC,CAAC,CAAC,CAAC;AACP,CAAC;AAED;;;;GAIG;AACH,MAAM,CAAC,KAAK,UAAU,eAAe,CAAC,MAAc,EAAE,MAAwB,EAAE,MAAoB;IAChG,MAAM,gBAAgB,EAAE,CAAC;IACzB,MAAM,EAAE,cAAc,EAAE,CAAC;IACzB,IAAI,UAAU,GAAG,WAAW,CAAC,GAAG,EAAE,CAAC;IACnC,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;QACzB,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,MAAM,GAAG,CAAC;YAClC,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,EAAE,KAAK,GAAG,KAAK,CAAC,CAAC;YAC5C,KAAK,CAAC,IAAI,CAAC,KAAK,EAAE,GAAG,CAAC,CAAC;YACvB,KAAK,GAAG,GAAG,CAAC;YACZ,IAAI,WAAW,CAAC,GAAG,EAAE,GAAG,UAAU,IAAI,eAAe,EAAE,CAAC;gBACpD,MAAM,gBAAgB,EAAE,CAAC;gBACzB,MAAM,EAAE,cAAc,EAAE,CAAC;gBACzB,UAAU,GAAG,WAAW,CAAC,GAAG,EAAE,CAAC;YACnC,CAAC;QACL,CAAC;QACD,KAAK,CAAC,IAAI,EAAE,EAAE,CAAC;IACnB,CAAC;AACL,CAAC;AAED,MAAM,UAAU,SAAS,CAAC,MAAoB,EAAE,KAAoC;IAChF,IAAI,KAAK,IAAI,KAAK,CAAC,MAAM,KAAK,MAAM,CAAC,MAAM,EAAE,CAAC;QAC1C,MAAM,IAAI,UAAU,CAAC,qCAAqC,KAAK,CAAC,MAAM,qBAAqB,MAAM,CAAC,MAAM,EAAE,CAAC,CAAC;IAChH,CAAC;AACL,CAAC;AAED,MAAM,UAAU,cAAc,CAAC,SAAiB;IAC5C,IAAI,CAAC,CAAC,SAAS,IAAI,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,SAAS,CAAC,EAAE,CAAC;QACnD,MAAM,IAAI,UAAU,CAAC,mEAAmE,SAAS,EAAE,CAAC,CAAC;IACzG,CAAC;AACL,CAAC;AAED,SAAS,MAAM,CAAC,MAAoB;IAChC,OAAO,IAAI,WAAW,CAAC,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,UAAU,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC;AAC5E,CAAC;AAED,MAAM,OAAO,GAAG,IAAI,WAAW,CAAC,CAAC,CAAC,CAAC;AACnC,MAAM,SAAS,GAAG,IAAI,YAAY,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;AACnD,SAAS,aAAa,CAAC,IAAY;IAC/B,OAAO,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;IAClB,OAAO,SAAS,CAAC,CAAC,CAAC,CAAC;AACxB,CAAC;AAED,yCAAyC;AACzC,MAAM,OAAO,cAAc;IAGF;IACA;IACA;IAJrB,KAAK,GAAG,CAAC,CAAC;IACV,YACqB,MAAoB,EACpB,SAAiB,EACjB,KAAyB;QAFzB,WAAM,GAAN,MAAM,CAAc;QACpB,cAAS,GAAT,SAAS,CAAQ;QACjB,UAAK,GAAL,KAAK,CAAoB;QAE1C,SAAS,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC;QACzB,cAAc,CAAC,SAAS,CAAC,CAAC;IAC9B,CAAC;IACD,IAAI,CAAC,KAAa,EAAE,GAAW;QAC3B,MAAM,EAAE,MAAM,EAAE,SAAS,EAAE,KAAK,EAAE,GAAG,IAAI,CAAC;QAC1C,IAAI,KAAK,GAAG,CAAC,CAAC;QACd,KAAK,IAAI,CAAC,GAAG,KAAK,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,sDAAsD;YACtD,IAAI,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,IAAI,SAAS,IAAI,CAAC,CAAC,KAAK,IAAI,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;gBAAE,KAAK,EAAE,CAAC;QAChF,CAAC;QACD,IAAI,CAAC,KAAK,IAAI,KAAK,CAAC;IACxB,CAAC;CACJ;AAED,6EAA6E;AAC7E,MAAM,OAAO,eAAe;IACf,SAAS,CAAqB;IACtB,IAAI,CAAS;IAC9B,YAAY,KAA8B;QACtC,MAAM,EAAE,MAAM,EAAE,SAAS,EAAE,IAAI,EAAE,GAAG,KAAK,CAAC;QAC1C,IAAI,CAAC,CAAC,SAAS,GAAG,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,SAAS,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,EAAE,CAAC;YAC9E,MAAM,IAAI,UAAU,CAAC,8DAA8D,SAAS,EAAE,CAAC,CAAC;QACpG,CAAC;QACD,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,IAAI,CAAC,IAAI,IAAI,GAAG,CAAC,EAAE,CAAC;YACtC,MAAM,IAAI,UAAU,CAAC,6DAA6D,IAAI,EAAE,CAAC,CAAC;QAC9F,CAAC;QACD,IAAI,CAAC,SAAS,GAAG;YACb,GAAG,EAAE,MAAM,GAAG,SAAS,EAAE,GAAG,EAAE,MAAM,GAAG,SAAS,EAAE,QAAQ,EAAE,CAAC,CAAC,GAAG,SAAS,CAAC,GAAG,IAAI;YAClF,MAAM,EAAE,IAAI,KAAK,CAAS,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC;SAC9D,CAAC;QACF,IAAI,CAAC,IAAI,GAAG,IAAI,GAAG,CAAC,CAAC;IACzB,CAAC;IACD,0BAA0B;IAC1B,GAAG,CAAC,CAAS;QACT,MAAM,CAAC,GAAG,IAAI,CAAC,SAAS,CAAC;QACzB,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG;YAAE,CAAC,CAAC,KAAK,EAAE,CAAC;aACpB,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG;YAAE,CAAC,CAAC,KAAK,EAAE,CAAC;;YACzB,CAAC,CAAC,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;IAC/E,CAAC;CACJ;AAED,sCAAsC;AACtC,MAAM,OAAO,aAAa;IAEO;IAAuE;IAD3F,MAAM,CAAkB;IACjC,YAA6B,MAAoB,EAAE,KAA8B,EAAmB,KAAyB;QAAhG,WAAM,GAAN,MAAM,CAAc;QAAmD,UAAK,GAAL,KAAK,CAAoB;QACzH,SAAS,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC;QACzB,IAAI,CAAC,MAAM,GAAG,IAAI,eAAe,CAAC,KAAK,CAAC,CAAC;IAC7C,CAAC;IACD,IAAI,CAAC,KAAa,EAAE,GAAW;QAC3B,MAAM,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,GAAG,IAAI,CAAC;QACvC,KAAK,IAAI,CAAC,GAAG,KAAK,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;YACpB,IAAI,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,KAAK,IAAI,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;gBAAE,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACxE,CAAC;IACL,CAAC;CACJ;AAED,wEAAwE;AACxE,MAAM,cAAc;IACR,GAAG,GAAG,CAAC,CAAC;IACR,YAAY,GAAG,CAAC,CAAC;IACzB,GAAG,CAAC,CAAS;QACT,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC;QACvB,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;YAAE,IAAI,CAAC,YAAY,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;;YAC1E,IAAI,CAAC,YAAY,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC;QAC7C,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC;IACjB,CAAC;IACD,KAAK;QACD,OAAO,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,YAAY,CAAC;IACxC,CAAC;CACJ;AAED,+EAA+E;AAC/E,MAAM,UAAU,WAAW,CAAC,OAAe,EAAE,CAAS;IAClD,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,OAAO,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AAC5E,CAAC;AAED,+CAA+C;AAC/C,MAAM,CAAC,MAAM,WAAW,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAU,CAAC;AAUjD;;;;GAIG;AACH,MAAM,OAAO,gBAAgB;IAsBI;IAAuC;IArB3D,MAAM,CAAmB;IAClC,CAAC,GAAG,CAAC,CAAC;IACN,MAAM,GAAG,CAAC,CAAC;IACX,OAAO,GAAG,CAAC,CAAC;IACZ,GAAG,GAAG,QAAQ,CAAC;IACf,GAAG,GAAG,CAAC,QAAQ,CAAC;IAChB,MAAM,GAAG,CAAC,CAAC;IACX,IAAI,GAAG,CAAC,CAAC;IACA,GAAG,GAAG,IAAI,cAAc,EAAE,CAAC;IAC3B,MAAM,GAAG,IAAI,cAAc,EAAE,CAAC;IAC9B,KAAK,GAAG,IAAI,cAAc,EAAE,CAAC;IAC7B,QAAQ,GAAG,IAAI,cAAc,EAAE,CAAC;IAChC,MAAM,CAAyB;IACxC,sEAAsE;IAC7D,WAAW,GAAa,EAAE,CAAC;IACnB,IAAI,CAAc;IAClB,IAAI,GAAG,IAAI,WAAW,CAAC,KAAK,CAAC,CAAC;IAC9B,MAAM,CAAqB;IAC5C,6EAA6E;IACrE,OAAO,GAA8D,EAAE,CAAC;IAEhF,YAA6B,MAAoB,EAAmB,OAAsB;QAA7D,WAAM,GAAN,MAAM,CAAc;QAAmB,YAAO,GAAP,OAAO,CAAe;QACtF,SAAS,CAAC,MAAM,EAAE,OAAO,CAAC,KAAK,CAAC,CAAC;QACjC,IAAI,OAAO,CAAC,SAAS,KAAK,SAAS;YAAE,cAAc,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC;QACvE,IAAI,CAAC,MAAM,GAAG,OAAO,CAAC,SAAS,CAAC,CAAC,CAAC,IAAI,eAAe,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC;QAChF,IAAI,CAAC,MAAM,GAAG,OAAO,CAAC,SAAS,KAAK,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC;QAC3F,IAAI,CAAC,IAAI,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC;QAC3B,IAAI,CAAC,MAAM,GAAG;YACV,EAAE,IAAI,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC,IAAI,CAAC,IAAI,EAAE,EAAE;YACnE,EAAE,IAAI,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,aAAa,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC,IAAI,CAAC,OAAO,EAAE,EAAE;SAC3E,CAAC;IACN,CAAC;IAEO,WAAW,CAAC,KAAa,EAAE,GAAW;QAC1C,MAAM,EAAE,MAAM,EAAE,IAAI,EAAE,IAAI,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,GAAG,IAAI,CAAC;QACxE,MAAM,EAAE,KAAK,EAAE,GAAG,IAAI,CAAC,OAAO,CAAC;QAC/B,uEAAuE;QACvE,MAAM,SAAS,GAAG,IAAI,CAAC,OAAO,CAAC,SAAS,IAAI,MAAM,CAAC,GAAG,CAAC;QACvD,MAAM,IAAI,GAAG,MAAM,IAAI,MAAM,CAAC,iBAAiB,CAAC;QAChD,IAAI,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,MAAM,EAAE,MAAM,EAAE,OAAO,EAAE,GAAG,IAAI,CAAC;QACpD,KAAK,IAAI,CAAC,GAAG,KAAK,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC;YAC7B,IAAI,CAAC,IAAI,UAAU,IAAI,CAAC,KAAK,IAAI,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;gBAAE,SAAS;YAC3D,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;YACpB,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YACxB,CAAC,EAAE,CAAC;YACJ,IAAI,CAAC,CAAC,KAAK,EAAE,CAAC,EAAE,CAAC;YACjB,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YACX,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YAChB,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACjB,IAAI,CAAC,GAAG,GAAG;gBAAE,GAAG,GAAG,CAAC,CAAC;YACrB,IAAI,CAAC,GAAG,GAAG;gBAAE,GAAG,GAAG,CAAC,CAAC;YACrB,IAAI,GAAG,GAAG,MAAM;gBAAE,MAAM,GAAG,GAAG,CAAC;YAC/B,IAAI,GAAG,IAAI,SAAS;gBAAE,MAAM,EAAE,CAAC;YAC/B,IAAI,GAAG,IAAI,IAAI;gBAAE,OAAO,EAAE,CAAC;YAC3B,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;QACnB,CAAC;QACD,MAAM,CAAC,MAAM,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,MAAM,EAAE,MAAM,EAAE,OAAO,EAAE,CAAC,CAAC;IAClE,CAAC;IAED,6EAA6E;IACrE,IAAI;QACR,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC;QACjB,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO;QACpB,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,GAAG,CAAC,CAAC;QACjC,MAAM,IAAI,GAAG,IAAI,GAAG,EAAuB,CAAC;QAC5C,IAAI,MAAM,GAAG,CAAC,CAAC;QACf,IAAI,KAAK,GAAG,CAAC,CAAC;QACd,qEAAqE;QACrE,KAAK,MAAM,OAAO,IAAI,WAAW,EAAE,CAAC;YAChC,MAAM,IAAI,GAAG,WAAW,CAAC,OAAO,EAAE,CAAC,CAAC,CAAC;YACrC,OAAO,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,IAAI,IAAI;gBAAE,KAAK,IAAI,IAAI,CAAC,IAAI,CAAC,MAAM,EAAE,CAAC,CAAC;YACvE,IAAI,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC;YAC3B,IAAI,CAAC,GAAG;gBAAE,IAAI,CAAC,GAAG,CAAC,MAAM,EAAE,GAAG,GAAG,IAAI,WAAW,CAAC,KAAK,CAAC,CAAC,CAAC;YACzD,IAAI,CAAC,OAAO,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,IAAI,EAAE,IAAI,GAAG,KAAK,EAAE,GAAG,EAAE,CAAC,CAAC;QAC3D,CAAC;IACL,CAAC;IAEO,aAAa,CAAC,KAAa,EAAE,GAAW;QAC5C,IAAI,IAAI,CAAC,CAAC,KAAK,CAAC;YAAE,OAAO;QACzB,MAAM,EAAE,MAAM,EAAE,IAAI,EAAE,QAAQ,EAAE,IAAI,EAAE,GAAG,IAAI,CAAC;QAC9C,MAAM,EAAE,KAAK,EAAE,GAAG,IAAI,CAAC,OAAO,CAAC;QAC/B,uEAAuE;QACvE,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;QAC3D,MAAM,EAAE,GAAG,EAAE,CAAC,MAAM,EAAE,EAAE,GAAG,EAAE,CAAC,MAAM,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,MAAM,CAAC;QAC7D,MAAM,EAAE,GAAG,EAAE,CAAC,MAAM,KAAK,EAAE,IAAI,EAAE,CAAC,MAAM,KAAK,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,MAAM,CAAC;QACxE,MAAM,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC;QAC5C,KAAK,IAAI,CAAC,GAAG,KAAK,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC;YAC7B,IAAI,CAAC,IAAI,UAAU,IAAI,CAAC,KAAK,IAAI,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;gBAAE,SAAS;YAC3D,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;YAC3B,QAAQ,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACpB,MAAM,CAAC,GAAG,CAAC,KAAK,EAAE,CAAC;YACnB,IAAI,CAAC,KAAK,EAAE;gBAAE,EAAE,CAAC,CAAC,GAAG,QAAQ,CAAC,EAAE,CAAC;iBAC5B,IAAI,CAAC,KAAK,EAAE;gBAAE,EAAE,CAAC,CAAC,GAAG,QAAQ,CAAC,EAAE,CAAC;iBACjC,IAAI,CAAC,KAAK,EAAE;gBAAE,EAAE,CAAC,CAAC,GAAG,QAAQ,CAAC,EAAE,CAAC;QAC1C,CAAC;IACL,CAAC;IAEO,OAAO;QACX,KAAK,MAAM,EAAE,MAAM,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,IAAI,CAAC,OAAO,EAAE,CAAC;YAC/C,IAAI,KAAK,GAAG,CAAC,CAAC;YACd,IAAI,CAAC,GAAG,CAAC,CAAC;YACV,OAAO,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC,IAAI,IAAI;gBAAE,KAAK,IAAI,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;YACjD,IAAI,CAAC,WAAW,CAAC,IAAI,CAAC,aAAa,CAAC,CAAC,MAAM,IAAI,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC7D,CAAC;IACL,CAAC;CACJ"}

package/dist/deviation/deviation-statistics.d.ts

@@ -0,0 +1,100 @@
+/** A bin layout that lines up with the deviation colour ramp. */
+export interface DeviationHistogramRange {
+    /** Ramp centre in metres (`deviationRange.centerOffset`). */
+    center: number;
+    /** Ramp half-width in metres; bins span [center − h, center + h]. */
+    halfRange: number;
+    /** Number of equal-width bins. */
+    bins: number;
+}
+export interface DeviationHistogram {
+    /** Lower edge of the first bin, metres. */
+    min: number;
+    /** Upper edge of the last bin (inclusive), metres. */
+    max: number;
+    binWidth: number;
+    /** Valid points per bin, low → high. */
+    counts: number[];
+    /** Valid points below `min` / above `max` (the ramp's saturated ends). */
+    below: number;
+    above: number;
+}
+export interface DeviationToleranceShare {
+    tolerance: number;
+    /** Valid points with |d| ≤ tolerance. */
+    count: number;
+    /** `count / validCount`; null when no point is valid. */
+    share: number | null;
+}
+export interface DeviationStatistics {
+    /** Every input point, valid or not. */
+    count: number;
+    /** Points that are unmasked and finite; every other figure is over these. */
+    validCount: number;
+    /** Valid points with |d| ≥ `clipRange` (pegged by the compute clip). */
+    clippedCount: number;
+    min: number | null;
+    max: number | null;
+    /** Signed mean. */
+    mean: number | null;
+    meanAbs: number | null;
+    rms: number | null;
+    /** Population standard deviation of the signed values. */
+    stdDev: number | null;
+    p50Abs: number | null;
+    p95Abs: number | null;
+    p99Abs: number | null;
+    maxAbs: number | null;
+    withinTolerance: DeviationToleranceShare | null;
+    histogram: DeviationHistogram | null;
+}
+export interface DeviationStatisticsOptions {
+    /** Per-point validity; zero excludes the point. Must match `values.length`. */
+    valid?: ArrayLike<number>;
+    /** Report the share of points with |d| ≤ tolerance (metres, ≥ 0). */
+    tolerance?: number;
+    /**
+     * The `maxRange` the compute pass clamped to, to count pegged points. The
+     * shader pegs at the float32 value, so the comparison is against
+     * `Math.fround(clipRange)`: 0.7 m pegs at 0.699999988 m.
+     */
+    clipRange?: number;
+    histogram?: DeviationHistogramRange;
+}
+/** One scan asset's slice of a shared readback array. */
+export interface DeviationAssetRange {
+    expressId: number;
+    modelIndex: number;
+    offset: number;
+    count: number;
+}
+/** Every computed point's signed distance, grouped by scan asset. */
+export interface DeviationDistances {
+    values: Float32Array;
+    assets: DeviationAssetRange[];
+}
+/** Cancellation for the sliced `…Async` variants. */
+export interface DeviationAsyncOptions {
+    /** Aborting rejects with the signal's reason at the next slice boundary. */
+    signal?: AbortSignal;
+}
+export interface DeviationAssetSummary {
+    expressId: number;
+    modelIndex: number;
+    statistics: DeviationStatistics;
+}
+/** Full summary in one synchronous call: two O(n) passes, no copy. */
+export declare function computeDeviationStatistics(values: Float32Array, options?: DeviationStatisticsOptions): DeviationStatistics;
+/** {@link computeDeviationStatistics} in time-boxed slices that yield to the event loop. */
+export declare function computeDeviationStatisticsAsync(values: Float32Array, options?: DeviationStatisticsOptions & DeviationAsyncOptions): Promise<DeviationStatistics>;
+/** Valid points with |d| ≤ tolerance: one O(n) pass in yielding slices, no allocation. */
+export declare function countWithinToleranceAsync(values: Float32Array, tolerance: number, options?: {
+    valid?: ArrayLike<number>;
+} & DeviationAsyncOptions): Promise<number>;
+/** Fixed-bin histogram over the ramp range: one O(n) pass in yielding slices. */
+export declare function deviationHistogramAsync(values: Float32Array, range: DeviationHistogramRange, options?: {
+    valid?: ArrayLike<number>;
+} & DeviationAsyncOptions): Promise<DeviationHistogram>;
+/** Per-asset statistics over each asset's slice of one readback array, in yielding slices. */
+export declare function summarizeDeviationAssetsAsync(distances: DeviationDistances, options?: Omit<DeviationStatisticsOptions, 'valid'> & DeviationAsyncOptions): Promise<DeviationAssetSummary[]>;
+//# sourceMappingURL=deviation-statistics.d.ts.map

package/dist/deviation/deviation-statistics.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"deviation-statistics.d.ts","sourceRoot":"","sources":["../../src/deviation/deviation-statistics.ts"],"names":[],"mappings":"AAsCA,iEAAiE;AACjE,MAAM,WAAW,uBAAuB;IACpC,6DAA6D;IAC7D,MAAM,EAAE,MAAM,CAAC;IACf,qEAAqE;IACrE,SAAS,EAAE,MAAM,CAAC;IAClB,kCAAkC;IAClC,IAAI,EAAE,MAAM,CAAC;CAChB;AAED,MAAM,WAAW,kBAAkB;IAC/B,2CAA2C;IAC3C,GAAG,EAAE,MAAM,CAAC;IACZ,sDAAsD;IACtD,GAAG,EAAE,MAAM,CAAC;IACZ,QAAQ,EAAE,MAAM,CAAC;IACjB,wCAAwC;IACxC,MAAM,EAAE,MAAM,EAAE,CAAC;IACjB,0EAA0E;IAC1E,KAAK,EAAE,MAAM,CAAC;IACd,KAAK,EAAE,MAAM,CAAC;CACjB;AAED,MAAM,WAAW,uBAAuB;IACpC,SAAS,EAAE,MAAM,CAAC;IAClB,yCAAyC;IACzC,KAAK,EAAE,MAAM,CAAC;IACd,yDAAyD;IACzD,KAAK,EAAE,MAAM,GAAG,IAAI,CAAC;CACxB;AAED,MAAM,WAAW,mBAAmB;IAChC,uCAAuC;IACvC,KAAK,EAAE,MAAM,CAAC;IACd,6EAA6E;IAC7E,UAAU,EAAE,MAAM,CAAC;IACnB,wEAAwE;IACxE,YAAY,EAAE,MAAM,CAAC;IACrB,GAAG,EAAE,MAAM,GAAG,IAAI,CAAC;IACnB,GAAG,EAAE,MAAM,GAAG,IAAI,CAAC;IACnB,mBAAmB;IACnB,IAAI,EAAE,MAAM,GAAG,IAAI,CAAC;IACpB,OAAO,EAAE,MAAM,GAAG,IAAI,CAAC;IACvB,GAAG,EAAE,MAAM,GAAG,IAAI,CAAC;IACnB,0DAA0D;IAC1D,MAAM,EAAE,MAAM,GAAG,IAAI,CAAC;IACtB,MAAM,EAAE,MAAM,GAAG,IAAI,CAAC;IACtB,MAAM,EAAE,MAAM,GAAG,IAAI,CAAC;IACtB,MAAM,EAAE,MAAM,GAAG,IAAI,CAAC;IACtB,MAAM,EAAE,MAAM,GAAG,IAAI,CAAC;IACtB,eAAe,EAAE,uBAAuB,GAAG,IAAI,CAAC;IAChD,SAAS,EAAE,kBAAkB,GAAG,IAAI,CAAC;CACxC;AAED,MAAM,WAAW,0BAA0B;IACvC,+EAA+E;IAC/E,KAAK,CAAC,EAAE,SAAS,CAAC,MAAM,CAAC,CAAC;IAC1B,qEAAqE;IACrE,SAAS,CAAC,EAAE,MAAM,CAAC;IACnB;;;;OAIG;IACH,SAAS,CAAC,EAAE,MAAM,CAAC;IACnB,SAAS,CAAC,EAAE,uBAAuB,CAAC;CACvC;AAED,yDAAyD;AACzD,MAAM,WAAW,mBAAmB;IAChC,SAAS,EAAE,MAAM,CAAC;IAClB,UAAU,EAAE,MAAM,CAAC;IACnB,MAAM,EAAE,MAAM,CAAC;IACf,KAAK,EAAE,MAAM,CAAC;CACjB;AAED,qEAAqE;AACrE,MAAM,WAAW,kBAAkB;IAC/B,MAAM,EAAE,YAAY,CAAC;IACrB,MAAM,EAAE,mBAAmB,EAAE,CAAC;CACjC;AAED,qDAAqD;AACrD,MAAM,WAAW,qBAAqB;IAClC,4EAA4E;IAC5E,MAAM,CAAC,EAAE,WAAW,CAAC;CACxB;AAED,MAAM,WAAW,qBAAqB;IAClC,SAAS,EAAE,MAAM,CAAC;IAClB,UAAU,EAAE,MAAM,CAAC;IACnB,UAAU,EAAE,mBAAmB,CAAC;CACnC;AAqCD,sEAAsE;AACtE,wBAAgB,0BAA0B,CACtC,MAAM,EAAE,YAAY,EACpB,OAAO,GAAE,0BAA+B,GACzC,mBAAmB,CAIrB;AAED,4FAA4F;AAC5F,wBAAsB,+BAA+B,CACjD,MAAM,EAAE,YAAY,EACpB,OAAO,GAAE,0BAA0B,GAAG,qBAA0B,GACjE,OAAO,CAAC,mBAAmB,CAAC,CAI9B;AAED,0FAA0F;AAC1F,wBAAsB,yBAAyB,CAC3C,MAAM,EAAE,YAAY,EACpB,SAAS,EAAE,MAAM,EACjB,OAAO,GAAE;IAAE,KAAK,CAAC,EAAE,SAAS,CAAC,MAAM,CAAC,CAAA;CAAE,GAAG,qBAA0B,GACpE,OAAO,CAAC,MAAM,CAAC,CAIjB;AAED,iFAAiF;AACjF,wBAAsB,uBAAuB,CACzC,MAAM,EAAE,YAAY,EACpB,KAAK,EAAE,uBAAuB,EAC9B,OAAO,GAAE;IAAE,KAAK,CAAC,EAAE,SAAS,CAAC,MAAM,CAAC,CAAA;CAAE,GAAG,qBAA0B,GACpE,OAAO,CAAC,kBAAkB,CAAC,CAI7B;AAED,8FAA8F;AAC9F,wBAAsB,6BAA6B,CAC/C,SAAS,EAAE,kBAAkB,EAC7B,OAAO,GAAE,IAAI,CAAC,0BAA0B,EAAE,OAAO,CAAC,GAAG,qBAA0B,GAChF,OAAO,CAAC,qBAAqB,EAAE,CAAC,CAalC"}

package/dist/deviation/deviation-statistics.js

@@ -0,0 +1,100 @@
+/* This Source Code Form is subject to the terms of the Mozilla Public
+ * License, v. 2.0. If a copy of the MPL was not distributed with this
+ * file, You can obtain one at https://mozilla.org/MPL/2.0/. */
+/**
+ * Summary statistics over a BIM ↔ scan deviation readback (#6872).
+ *
+ * Pure functions of a `Float32Array` of signed distances in metres (one per
+ * scan point, positive = outside the surface), plus an optional validity
+ * mask. A value is VALID when the mask (if any) is non-zero AND the value is
+ * finite: NaN and ±Infinity never enter a moment, a percentile or a count.
+ *
+ * Accuracy:
+ * - Sums use Neumaier (improved Kahan) compensated summation in float64, and
+ *   the standard deviation is a second pass over (d − mean)², so neither the
+ *   mean nor σ loses digits to cancellation on large clouds.
+ * - Percentiles of |d| are EXACT nearest-rank order statistics
+ *   (rank ⌈p·n⌉, the convention `frame-timing-stats.ts` uses), found by a
+ *   two-pass radix select on the float32 bit patterns
+ *   (`deviation-statistics-kernel.ts`): O(n) worst case, a fixed 1 MiB of
+ *   count tables whatever n is, and the input is never copied or reordered.
+ *
+ * Main thread: every figure is at most two linear passes. The `…Async`
+ * variants run those passes in slices of a few milliseconds that yield to
+ * the event loop, so a 25M-point summary never blocks one task for long, and
+ * take an `AbortSignal` so a superseded request stops early. The viewer only
+ * uses those; the synchronous {@link computeDeviationStatistics} suits a
+ * worker, a CLI or a test.
+ */
+import { HistogramPass, StatisticsPasses, ToleranceCount, runPhases, runPhasesSliced, } from './deviation-statistics-kernel.js';
+function statisticsOf(passes, count, tolerance) {
+    const { n } = passes;
+    const withinTolerance = tolerance === undefined
+        ? null
+        : { tolerance, count: passes.within, share: n > 0 ? passes.within / n : null };
+    const histogram = passes.binner?.histogram ?? null;
+    if (n === 0) {
+        return {
+            count, validCount: 0, clippedCount: 0,
+            min: null, max: null, mean: null, meanAbs: null, rms: null, stdDev: null,
+            p50Abs: null, p95Abs: null, p99Abs: null, maxAbs: null,
+            withinTolerance, histogram,
+        };
+    }
+    const [p50Abs, p95Abs, p99Abs] = passes.percentiles;
+    return {
+        count,
+        validCount: n,
+        clippedCount: passes.clipped,
+        min: passes.min,
+        max: passes.max,
+        mean: passes.mean,
+        meanAbs: passes.sumAbs.value() / n,
+        rms: Math.sqrt(passes.sumSq.value() / n),
+        // σ from squared residuals about the mean, not E[d²] − mean², which cancels.
+        stdDev: Math.sqrt(passes.residual.value() / n),
+        p50Abs,
+        p95Abs,
+        p99Abs,
+        maxAbs: passes.maxAbs,
+        withinTolerance,
+        histogram,
+    };
+}
+/** Full summary in one synchronous call: two O(n) passes, no copy. */
+export function computeDeviationStatistics(values, options = {}) {
+    const passes = new StatisticsPasses(values, options);
+    runPhases(values.length, passes.phases);
+    return statisticsOf(passes, values.length, options.tolerance);
+}
+/** {@link computeDeviationStatistics} in time-boxed slices that yield to the event loop. */
+export async function computeDeviationStatisticsAsync(values, options = {}) {
+    const passes = new StatisticsPasses(values, options);
+    await runPhasesSliced(values.length, passes.phases, options.signal);
+    return statisticsOf(passes, values.length, options.tolerance);
+}
+/** Valid points with |d| ≤ tolerance: one O(n) pass in yielding slices, no allocation. */
+export async function countWithinToleranceAsync(values, tolerance, options = {}) {
+    const pass = new ToleranceCount(values, tolerance, options.valid);
+    await runPhasesSliced(values.length, [pass], options.signal);
+    return pass.count;
+}
+/** Fixed-bin histogram over the ramp range: one O(n) pass in yielding slices. */
+export async function deviationHistogramAsync(values, range, options = {}) {
+    const pass = new HistogramPass(values, range, options.valid);
+    await runPhasesSliced(values.length, [pass], options.signal);
+    return pass.binner.histogram;
+}
+/** Per-asset statistics over each asset's slice of one readback array, in yielding slices. */
+export async function summarizeDeviationAssetsAsync(distances, options = {}) {
+    const summaries = [];
+    for (const asset of distances.assets) {
+        summaries.push({
+            expressId: asset.expressId,
+            modelIndex: asset.modelIndex,
+            statistics: await computeDeviationStatisticsAsync(distances.values.subarray(asset.offset, asset.offset + asset.count), options),
+        });
+    }
+    return summaries;
+}
+//# sourceMappingURL=deviation-statistics.js.map

package/dist/deviation/deviation-statistics.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"deviation-statistics.js","sourceRoot":"","sources":["../../src/deviation/deviation-statistics.ts"],"names":[],"mappings":"AAAA;;+DAE+D;AAE/D;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AAEH,OAAO,EACH,aAAa,EACb,gBAAgB,EAChB,cAAc,EACd,SAAS,EACT,eAAe,GAClB,MAAM,kCAAkC,CAAC;AAgG1C,SAAS,YAAY,CAAC,MAAwB,EAAE,KAAa,EAAE,SAA6B;IACxF,MAAM,EAAE,CAAC,EAAE,GAAG,MAAM,CAAC;IACrB,MAAM,eAAe,GAAG,SAAS,KAAK,SAAS;QAC3C,CAAC,CAAC,IAAI;QACN,CAAC,CAAC,EAAE,SAAS,EAAE,KAAK,EAAE,MAAM,CAAC,MAAM,EAAE,KAAK,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC;IACnF,MAAM,SAAS,GAAG,MAAM,CAAC,MAAM,EAAE,SAAS,IAAI,IAAI,CAAC;IACnD,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QACV,OAAO;YACH,KAAK,EAAE,UAAU,EAAE,CAAC,EAAE,YAAY,EAAE,CAAC;YACrC,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,OAAO,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,MAAM,EAAE,IAAI;YACxE,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,IAAI;YACtD,eAAe,EAAE,SAAS;SAC7B,CAAC;IACN,CAAC;IACD,MAAM,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,GAAG,MAAM,CAAC,WAAW,CAAC;IACpD,OAAO;QACH,KAAK;QACL,UAAU,EAAE,CAAC;QACb,YAAY,EAAE,MAAM,CAAC,OAAO;QAC5B,GAAG,EAAE,MAAM,CAAC,GAAG;QACf,GAAG,EAAE,MAAM,CAAC,GAAG;QACf,IAAI,EAAE,MAAM,CAAC,IAAI;QACjB,OAAO,EAAE,MAAM,CAAC,MAAM,CAAC,KAAK,EAAE,GAAG,CAAC;QAClC,GAAG,EAAE,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,KAAK,EAAE,GAAG,CAAC,CAAC;QACxC,6EAA6E;QAC7E,MAAM,EAAE,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,EAAE,GAAG,CAAC,CAAC;QAC9C,MAAM;QACN,MAAM;QACN,MAAM;QACN,MAAM,EAAE,MAAM,CAAC,MAAM;QACrB,eAAe;QACf,SAAS;KACZ,CAAC;AACN,CAAC;AAED,sEAAsE;AACtE,MAAM,UAAU,0BAA0B,CACtC,MAAoB,EACpB,UAAsC,EAAE;IAExC,MAAM,MAAM,GAAG,IAAI,gBAAgB,CAAC,MAAM,EAAE,OAAO,CAAC,CAAC;IACrD,SAAS,CAAC,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC;IACxC,OAAO,YAAY,CAAC,MAAM,EAAE,MAAM,CAAC,MAAM,EAAE,OAAO,CAAC,SAAS,CAAC,CAAC;AAClE,CAAC;AAED,4FAA4F;AAC5F,MAAM,CAAC,KAAK,UAAU,+BAA+B,CACjD,MAAoB,EACpB,UAA8D,EAAE;IAEhE,MAAM,MAAM,GAAG,IAAI,gBAAgB,CAAC,MAAM,EAAE,OAAO,CAAC,CAAC;IACrD,MAAM,eAAe,CAAC,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,MAAM,EAAE,OAAO,CAAC,MAAM,CAAC,CAAC;IACpE,OAAO,YAAY,CAAC,MAAM,EAAE,MAAM,CAAC,MAAM,EAAE,OAAO,CAAC,SAAS,CAAC,CAAC;AAClE,CAAC;AAED,0FAA0F;AAC1F,MAAM,CAAC,KAAK,UAAU,yBAAyB,CAC3C,MAAoB,EACpB,SAAiB,EACjB,UAAiE,EAAE;IAEnE,MAAM,IAAI,GAAG,IAAI,cAAc,CAAC,MAAM,EAAE,SAAS,EAAE,OAAO,CAAC,KAAK,CAAC,CAAC;IAClE,MAAM,eAAe,CAAC,MAAM,CAAC,MAAM,EAAE,CAAC,IAAI,CAAC,EAAE,OAAO,CAAC,MAAM,CAAC,CAAC;IAC7D,OAAO,IAAI,CAAC,KAAK,CAAC;AACtB,CAAC;AAED,iFAAiF;AACjF,MAAM,CAAC,KAAK,UAAU,uBAAuB,CACzC,MAAoB,EACpB,KAA8B,EAC9B,UAAiE,EAAE;IAEnE,MAAM,IAAI,GAAG,IAAI,aAAa,CAAC,MAAM,EAAE,KAAK,EAAE,OAAO,CAAC,KAAK,CAAC,CAAC;IAC7D,MAAM,eAAe,CAAC,MAAM,CAAC,MAAM,EAAE,CAAC,IAAI,CAAC,EAAE,OAAO,CAAC,MAAM,CAAC,CAAC;IAC7D,OAAO,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC;AACjC,CAAC;AAED,8FAA8F;AAC9F,MAAM,CAAC,KAAK,UAAU,6BAA6B,CAC/C,SAA6B,EAC7B,UAA6E,EAAE;IAE/E,MAAM,SAAS,GAA4B,EAAE,CAAC;IAC9C,KAAK,MAAM,KAAK,IAAI,SAAS,CAAC,MAAM,EAAE,CAAC;QACnC,SAAS,CAAC,IAAI,CAAC;YACX,SAAS,EAAE,KAAK,CAAC,SAAS;YAC1B,UAAU,EAAE,KAAK,CAAC,UAAU;YAC5B,UAAU,EAAE,MAAM,+BAA+B,CAC7C,SAAS,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,MAAM,EAAE,KAAK,CAAC,MAAM,GAAG,KAAK,CAAC,KAAK,CAAC,EACnE,OAAO,CACV;SACJ,CAAC,CAAC;IACP,CAAC;IACD,OAAO,SAAS,CAAC;AACrB,CAAC"}