@ifc-lite/renderer 6.0.0 → 6.3.0

package/dist/deviation/deviation-readback.js

@@ -1,22 +1,33 @@
 /* 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/. */
+import { summarizeDeviationAssetsAsync, } from './deviation-statistics.js';
 /**
- * Read signed distances only for an explicit export. One staging buffer is
- * alive at a time, so a large scan never needs a second full copy in RAM or
- * VRAM. The normal deviation compute and render path does no readback.
+ * Read every computed point's signed distance back from the GPU, on demand.
+ * The normal compute and render path does no readback.
+ *
+ * The result is ONE `Float32Array` (4 bytes per computed point, the size of
+ * the GPU deviation buffers themselves) with each scan asset as a contiguous
+ * range, so per-asset and whole-run statistics share it without a second
+ * copy. One staging buffer is alive at a time, so VRAM grows by one chunk.
  */
-export async function readDeviationAssetStats(device, nodes, wasComputed) {
-    const results = [];
+export async function readDeviationDistances(device, nodes, wasComputed) {
+    const assets = [];
+    const plan = [];
+    let total = 0;
     for (const node of nodes) {
-        let count = 0;
-        let finiteCount = 0;
-        let sum = 0;
-        let min = Infinity;
-        let max = -Infinity;
-        for (const chunk of node.chunks) {
-            if (!wasComputed(chunk.deviationBuffer) || chunk.pointCount === 0)
-                continue;
+        const chunks = node.chunks.filter((chunk) => chunk.pointCount > 0 && wasComputed(chunk.deviationBuffer));
+        const count = chunks.reduce((sum, chunk) => sum + chunk.pointCount, 0);
+        if (count === 0)
+            continue;
+        assets.push({ expressId: node.meta.expressId, modelIndex: node.meta.modelIndex ?? 0, offset: total, count });
+        plan.push(chunks);
+        total += count;
+    }
+    const values = new Float32Array(total);
+    let cursor = 0;
+    for (const chunks of plan) {
+        for (const chunk of chunks) {
             const size = chunk.pointCount * Float32Array.BYTES_PER_ELEMENT;
             const staging = device.createBuffer({
                 size,
@@ -24,21 +35,13 @@ export async function readDeviationAssetStats(device, nodes, wasComputed) {
             });
             let mapped = false;
             try {
-                const encoder = device.createCommandEncoder({ label: 'deviation-csv-readback' });
+                const encoder = device.createCommandEncoder({ label: 'deviation-readback' });
                 encoder.copyBufferToBuffer(chunk.deviationBuffer, 0, staging, 0, size);
                 device.queue.submit([encoder.finish()]);
                 await staging.mapAsync(GPUMapMode.READ);
                 mapped = true;
-                const values = new Float32Array(staging.getMappedRange());
-                count += values.length;
-                for (const value of values) {
-                    if (!Number.isFinite(value))
-                        continue;
-                    finiteCount++;
-                    sum += value;
-                    min = Math.min(min, value);
-                    max = Math.max(max, value);
-                }
+                values.set(new Float32Array(staging.getMappedRange(), 0, chunk.pointCount), cursor);
+                cursor += chunk.pointCount;
             }
             finally {
                 try {
@@ -50,18 +53,25 @@ export async function readDeviationAssetStats(device, nodes, wasComputed) {
                 }
             }
         }
-        if (count === 0)
-            continue;
-        results.push({
-            expressId: node.meta.expressId,
-            modelIndex: node.meta.modelIndex ?? 0,
-            pointsProcessed: count,
-            finitePoints: finiteCount,
-            minimumDeviation: finiteCount > 0 ? min : null,
-            maximumDeviation: finiteCount > 0 ? max : null,
-            meanDeviation: finiteCount > 0 ? sum / finiteCount : null,
-        });
     }
-    return results;
+    return { values, assets };
+}
+/**
+ * Per-scan-asset statistics over {@link readDeviationDistances}. Holds the
+ * readback (4 bytes per point) for the duration of the call; the statistics
+ * themselves add a fixed ~1 MiB and run in slices that yield to the event loop.
+ */
+export async function readDeviationAssetStats(device, nodes, wasComputed) {
+    const distances = await readDeviationDistances(device, nodes, wasComputed);
+    return (await summarizeDeviationAssetsAsync(distances)).map(({ expressId, modelIndex, statistics }) => ({
+        expressId,
+        modelIndex,
+        pointsProcessed: statistics.count,
+        finitePoints: statistics.validCount,
+        minimumDeviation: statistics.min,
+        maximumDeviation: statistics.max,
+        meanDeviation: statistics.mean,
+        statistics,
+    }));
 }
 //# sourceMappingURL=deviation-readback.js.map

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

@@ -1 +1 @@
-{"version":3,"file":"deviation-readback.js","sourceRoot":"","sources":["../../src/deviation/deviation-readback.ts"],"names":[],"mappings":"AAAA;;+DAE+D;AAe/D;;;;GAIG;AACH,MAAM,CAAC,KAAK,UAAU,uBAAuB,CACzC,MAAiB,EACjB,KAA+B,EAC/B,WAA2C;IAE3C,MAAM,OAAO,GAA0B,EAAE,CAAC;IAC1C,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACvB,IAAI,KAAK,GAAG,CAAC,CAAC;QACd,IAAI,WAAW,GAAG,CAAC,CAAC;QACpB,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,IAAI,GAAG,GAAG,QAAQ,CAAC;QACnB,IAAI,GAAG,GAAG,CAAC,QAAQ,CAAC;QACpB,KAAK,MAAM,KAAK,IAAI,IAAI,CAAC,MAAM,EAAE,CAAC;YAC9B,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,eAAe,CAAC,IAAI,KAAK,CAAC,UAAU,KAAK,CAAC;gBAAE,SAAS;YAC5E,MAAM,IAAI,GAAG,KAAK,CAAC,UAAU,GAAG,YAAY,CAAC,iBAAiB,CAAC;YAC/D,MAAM,OAAO,GAAG,MAAM,CAAC,YAAY,CAAC;gBAChC,IAAI;gBACJ,KAAK,EAAE,cAAc,CAAC,QAAQ,GAAG,cAAc,CAAC,QAAQ;aAC3D,CAAC,CAAC;YACH,IAAI,MAAM,GAAG,KAAK,CAAC;YACnB,IAAI,CAAC;gBACD,MAAM,OAAO,GAAG,MAAM,CAAC,oBAAoB,CAAC,EAAE,KAAK,EAAE,wBAAwB,EAAE,CAAC,CAAC;gBACjF,OAAO,CAAC,kBAAkB,CAAC,KAAK,CAAC,eAAe,EAAE,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC;gBACvE,MAAM,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,OAAO,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;gBACxC,MAAM,OAAO,CAAC,QAAQ,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC;gBACxC,MAAM,GAAG,IAAI,CAAC;gBACd,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,OAAO,CAAC,cAAc,EAAE,CAAC,CAAC;gBAC1D,KAAK,IAAI,MAAM,CAAC,MAAM,CAAC;gBACvB,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;oBACzB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;wBAAE,SAAS;oBACtC,WAAW,EAAE,CAAC;oBACd,GAAG,IAAI,KAAK,CAAC;oBACb,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,KAAK,CAAC,CAAC;oBAC3B,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,KAAK,CAAC,CAAC;gBAC/B,CAAC;YACL,CAAC;oBAAS,CAAC;gBACP,IAAI,CAAC;oBACD,IAAI,MAAM;wBAAE,OAAO,CAAC,KAAK,EAAE,CAAC;gBAChC,CAAC;wBAAS,CAAC;oBACP,OAAO,CAAC,OAAO,EAAE,CAAC;gBACtB,CAAC;YACL,CAAC;QACL,CAAC;QACD,IAAI,KAAK,KAAK,CAAC;YAAE,SAAS;QAC1B,OAAO,CAAC,IAAI,CAAC;YACT,SAAS,EAAE,IAAI,CAAC,IAAI,CAAC,SAAS;YAC9B,UAAU,EAAE,IAAI,CAAC,IAAI,CAAC,UAAU,IAAI,CAAC;YACrC,eAAe,EAAE,KAAK;YACtB,YAAY,EAAE,WAAW;YACzB,gBAAgB,EAAE,WAAW,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI;YAC9C,gBAAgB,EAAE,WAAW,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI;YAC9C,aAAa,EAAE,WAAW,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,WAAW,CAAC,CAAC,CAAC,IAAI;SAC5D,CAAC,CAAC;IACP,CAAC;IACD,OAAO,OAAO,CAAC;AACnB,CAAC"}
+{"version":3,"file":"deviation-readback.js","sourceRoot":"","sources":["../../src/deviation/deviation-readback.ts"],"names":[],"mappings":"AAAA;;+DAE+D;AAG/D,OAAO,EACH,6BAA6B,GAIhC,MAAM,2BAA2B,CAAC;AAenC;;;;;;;;GAQG;AACH,MAAM,CAAC,KAAK,UAAU,sBAAsB,CACxC,MAAiB,EACjB,KAA+B,EAC/B,WAA2C;IAE3C,MAAM,MAAM,GAA0B,EAAE,CAAC;IACzC,MAAM,IAAI,GAAoC,EAAE,CAAC;IACjD,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACvB,MAAM,MAAM,GAAG,IAAI,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,KAAK,CAAC,UAAU,GAAG,CAAC,IAAI,WAAW,CAAC,KAAK,CAAC,eAAe,CAAC,CAAC,CAAC;QACzG,MAAM,KAAK,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,KAAK,EAAE,EAAE,CAAC,GAAG,GAAG,KAAK,CAAC,UAAU,EAAE,CAAC,CAAC,CAAC;QACvE,IAAI,KAAK,KAAK,CAAC;YAAE,SAAS;QAC1B,MAAM,CAAC,IAAI,CAAC,EAAE,SAAS,EAAE,IAAI,CAAC,IAAI,CAAC,SAAS,EAAE,UAAU,EAAE,IAAI,CAAC,IAAI,CAAC,UAAU,IAAI,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,KAAK,EAAE,CAAC,CAAC;QAC7G,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QAClB,KAAK,IAAI,KAAK,CAAC;IACnB,CAAC;IACD,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,KAAK,CAAC,CAAC;IACvC,IAAI,MAAM,GAAG,CAAC,CAAC;IACf,KAAK,MAAM,MAAM,IAAI,IAAI,EAAE,CAAC;QACxB,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;YACzB,MAAM,IAAI,GAAG,KAAK,CAAC,UAAU,GAAG,YAAY,CAAC,iBAAiB,CAAC;YAC/D,MAAM,OAAO,GAAG,MAAM,CAAC,YAAY,CAAC;gBAChC,IAAI;gBACJ,KAAK,EAAE,cAAc,CAAC,QAAQ,GAAG,cAAc,CAAC,QAAQ;aAC3D,CAAC,CAAC;YACH,IAAI,MAAM,GAAG,KAAK,CAAC;YACnB,IAAI,CAAC;gBACD,MAAM,OAAO,GAAG,MAAM,CAAC,oBAAoB,CAAC,EAAE,KAAK,EAAE,oBAAoB,EAAE,CAAC,CAAC;gBAC7E,OAAO,CAAC,kBAAkB,CAAC,KAAK,CAAC,eAAe,EAAE,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC;gBACvE,MAAM,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,OAAO,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;gBACxC,MAAM,OAAO,CAAC,QAAQ,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC;gBACxC,MAAM,GAAG,IAAI,CAAC;gBACd,MAAM,CAAC,GAAG,CAAC,IAAI,YAAY,CAAC,OAAO,CAAC,cAAc,EAAE,EAAE,CAAC,EAAE,KAAK,CAAC,UAAU,CAAC,EAAE,MAAM,CAAC,CAAC;gBACpF,MAAM,IAAI,KAAK,CAAC,UAAU,CAAC;YAC/B,CAAC;oBAAS,CAAC;gBACP,IAAI,CAAC;oBACD,IAAI,MAAM;wBAAE,OAAO,CAAC,KAAK,EAAE,CAAC;gBAChC,CAAC;wBAAS,CAAC;oBACP,OAAO,CAAC,OAAO,EAAE,CAAC;gBACtB,CAAC;YACL,CAAC;QACL,CAAC;IACL,CAAC;IACD,OAAO,EAAE,MAAM,EAAE,MAAM,EAAE,CAAC;AAC9B,CAAC;AAED;;;;GAIG;AACH,MAAM,CAAC,KAAK,UAAU,uBAAuB,CACzC,MAAiB,EACjB,KAA+B,EAC/B,WAA2C;IAE3C,MAAM,SAAS,GAAG,MAAM,sBAAsB,CAAC,MAAM,EAAE,KAAK,EAAE,WAAW,CAAC,CAAC;IAC3E,OAAO,CAAC,MAAM,6BAA6B,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,SAAS,EAAE,UAAU,EAAE,UAAU,EAAE,EAAE,EAAE,CAAC,CAAC;QACpG,SAAS;QACT,UAAU;QACV,eAAe,EAAE,UAAU,CAAC,KAAK;QACjC,YAAY,EAAE,UAAU,CAAC,UAAU;QACnC,gBAAgB,EAAE,UAAU,CAAC,GAAG;QAChC,gBAAgB,EAAE,UAAU,CAAC,GAAG;QAChC,aAAa,EAAE,UAAU,CAAC,IAAI;QAC9B,UAAU;KACb,CAAC,CAAC,CAAC;AACR,CAAC"}

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