pca-web 0.1.0

package/dist/webgpu/gpu-pca.d.ts

@@ -0,0 +1,90 @@
+import { Matrix, type MatrixInput } from '../matrix.js';
+import type { PCAModel } from '../model.js';
+import { type PCAOptions, type SvdSolver } from '../pca.js';
+import { type FitAsyncOptions } from '../progress.js';
+import type { FloatArray } from '../types.js';
+import { type WebGPUDeviceOptions } from './engine.js';
+export interface WebGPUPCAOptions extends PCAOptions, WebGPUDeviceOptions {
+    /**
+     * Minimum element count (rows × cols) before the GPU path is used;
+     * smaller problems run on the CPU, where they are faster anyway.
+     * Default 1<<18 (≈262k elements).
+     */
+    minGpuElements?: number;
+}
+export declare class WebGPUPCA {
+    private cpu;
+    private readonly deviceOptions;
+    private readonly minGpuElements;
+    private enginePromise;
+    private engine;
+    private backend_;
+    private gpuFitting;
+    constructor(options?: WebGPUPCAOptions);
+    /** True when a WebGPU implementation is present in this environment. */
+    static isSupported(): boolean;
+    /** 'webgpu' when the last fit/transform used the GPU, else 'cpu'. */
+    get backend(): 'webgpu' | 'cpu';
+    /** Adapter description of the engine in use, or null on the CPU fallback. */
+    get gpuAdapterInfo(): string | null;
+    /** Measured df64 GEMM relative error from the device self-check (NaN on CPU). */
+    get gpuMeasuredRelError(): number;
+    /** Releases GPU resources (does not destroy an injected device). */
+    dispose(): void;
+    fit(X: MatrixInput, options?: FitAsyncOptions): Promise<this>;
+    fitTransform(X: MatrixInput, options?: FitAsyncOptions): Promise<Matrix>;
+    /** Concurrent fits on one instance would interleave GPU state; throw early. */
+    private guardConcurrentFit;
+    /**
+     * Validates the input (before any GPU spend or in-place mutation) and
+     * decides which backend the fit runs on. Mutates nothing.
+     */
+    private route;
+    /** GPU Gram product + shared CPU covariance_eigh tail. Mutates nothing on failure. */
+    private fitGramGpu;
+    /**
+     * GPU-driven randomized fit: centers on the CPU (respecting copy
+     * semantics), keeps X resident on the device, and drives the shared
+     * algorithm generator with GPU GEMMs. Returns the (sign-flipped)
+     * decomposition for the fitTransform fast path.
+     */
+    private fitRandomizedGpu;
+    /**
+     * CPU refit after a GPU failure — unless the input was already mutated,
+     * or the "failure" is a cancellation: an abort must reject with the abort
+     * error, never silently refit on the CPU.
+     */
+    private recoverOnGpuError;
+    transform(X: MatrixInput): Promise<Matrix>;
+    /** CPU delegate (small k×p GEMM; not worth a device round-trip). */
+    inverseTransform(X: MatrixInput): Matrix;
+    get components(): Matrix;
+    get explainedVariance(): FloatArray;
+    get explainedVarianceRatio(): FloatArray;
+    get singularValues(): FloatArray;
+    get mean(): FloatArray;
+    get nComponents(): number;
+    get nSamples(): number;
+    get nFeaturesIn(): number;
+    get noiseVariance(): number;
+    get whiten(): boolean;
+    get resolvedSvdSolver(): Exclude<SvdSolver, 'auto'>;
+    getCovariance(): Matrix;
+    getPrecision(): Matrix;
+    scoreSamples(X: MatrixInput): FloatArray;
+    score(X: MatrixInput): number;
+    getFeatureNamesOut(): string[];
+    /** The fitted model as a plain object — identical to the CPU class's. */
+    toModel(): PCAModel;
+    /**
+     * Rehydrates a fitted WebGPUPCA. `options` supplies the device knobs
+     * (powerPreference, injected device, minGpuElements); the PCA options
+     * come from the model itself.
+     */
+    static fromModel(model: PCAModel, options?: WebGPUDeviceOptions & {
+        minGpuElements?: number;
+    }): WebGPUPCA;
+    private ensureEngine;
+    private gpuEligible;
+}
+//# sourceMappingURL=gpu-pca.d.ts.map

package/dist/webgpu/gpu-pca.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"gpu-pca.d.ts","sourceRoot":"","sources":["../../src/webgpu/gpu-pca.ts"],"names":[],"mappings":"AAiBA,OAAO,EAAY,MAAM,EAAE,KAAK,WAAW,EAAE,MAAM,cAAc,CAAC;AAClE,OAAO,KAAK,EAAE,QAAQ,EAAE,MAAM,aAAa,CAAC;AAK5C,OAAO,EAEL,KAAK,UAAU,EAEf,KAAK,SAAS,EAEf,MAAM,WAAW,CAAC;AACnB,OAAO,EAAE,KAAK,eAAe,EAAgB,MAAM,gBAAgB,CAAC;AAEpE,OAAO,KAAK,EAAE,UAAU,EAAE,MAAM,aAAa,CAAC;AAE9C,OAAO,EAAgC,KAAK,mBAAmB,EAAE,MAAM,aAAa,CAAC;AAErF,MAAM,WAAW,gBAAiB,SAAQ,UAAU,EAAE,mBAAmB;IACvE;;;;OAIG;IACH,cAAc,CAAC,EAAE,MAAM,CAAC;CACzB;AAWD,qBAAa,SAAS;IACpB,OAAO,CAAC,GAAG,CAAM;IACjB,OAAO,CAAC,QAAQ,CAAC,aAAa,CAAsB;IACpD,OAAO,CAAC,QAAQ,CAAC,cAAc,CAAS;IACxC,OAAO,CAAC,aAAa,CAA0C;IAC/D,OAAO,CAAC,MAAM,CAA0B;IACxC,OAAO,CAAC,QAAQ,CAA2B;IAC3C,OAAO,CAAC,UAAU,CAAS;gBAEf,OAAO,GAAE,gBAAqB;IAO1C,wEAAwE;IACxE,MAAM,CAAC,WAAW,IAAI,OAAO;IAI7B,qEAAqE;IACrE,IAAI,OAAO,IAAI,QAAQ,GAAG,KAAK,CAE9B;IAED,6EAA6E;IAC7E,IAAI,cAAc,IAAI,MAAM,GAAG,IAAI,CAElC;IAED,iFAAiF;IACjF,IAAI,mBAAmB,IAAI,MAAM,CAEhC;IAED,oEAAoE;IACpE,OAAO,IAAI,IAAI;IAUT,GAAG,CAAC,CAAC,EAAE,WAAW,EAAE,OAAO,GAAE,eAAoB,GAAG,OAAO,CAAC,IAAI,CAAC;IA2BjE,YAAY,CAAC,CAAC,EAAE,WAAW,EAAE,OAAO,GAAE,eAAoB,GAAG,OAAO,CAAC,MAAM,CAAC;IA4ClF,+EAA+E;IAC/E,OAAO,CAAC,kBAAkB;IAS1B;;;OAGG;YACW,KAAK;IAoBnB,sFAAsF;YACxE,UAAU;IA0BxB;;;;;OAKG;YACW,gBAAgB;IA6D9B;;;;OAIG;YACW,iBAAiB;IAsBzB,SAAS,CAAC,CAAC,EAAE,WAAW,GAAG,OAAO,CAAC,MAAM,CAAC;IA2BhD,oEAAoE;IACpE,gBAAgB,CAAC,CAAC,EAAE,WAAW,GAAG,MAAM;IAQxC,IAAI,UAAU,IAAI,MAAM,CAEvB;IACD,IAAI,iBAAiB,IAAI,UAAU,CAElC;IACD,IAAI,sBAAsB,IAAI,UAAU,CAEvC;IACD,IAAI,cAAc,IAAI,UAAU,CAE/B;IACD,IAAI,IAAI,IAAI,UAAU,CAErB;IACD,IAAI,WAAW,IAAI,MAAM,CAExB;IACD,IAAI,QAAQ,IAAI,MAAM,CAErB;IACD,IAAI,WAAW,IAAI,MAAM,CAExB;IACD,IAAI,aAAa,IAAI,MAAM,CAE1B;IACD,IAAI,MAAM,IAAI,OAAO,CAEpB;IACD,IAAI,iBAAiB,IAAI,OAAO,CAAC,SAAS,EAAE,MAAM,CAAC,CAElD;IAED,aAAa,IAAI,MAAM;IAGvB,YAAY,IAAI,MAAM;IAGtB,YAAY,CAAC,CAAC,EAAE,WAAW,GAAG,UAAU;IAGxC,KAAK,CAAC,CAAC,EAAE,WAAW,GAAG,MAAM;IAG7B,kBAAkB,IAAI,MAAM,EAAE;IAQ9B,yEAAyE;IACzE,OAAO,IAAI,QAAQ;IAInB;;;;OAIG;IACH,MAAM,CAAC,SAAS,CACd,KAAK,EAAE,QAAQ,EACf,OAAO,GAAE,mBAAmB,GAAG;QAAE,cAAc,CAAC,EAAE,MAAM,CAAA;KAAO,GAC9D,SAAS;YAQE,YAAY;IAU1B,OAAO,CAAC,WAAW;CASpB"}

package/dist/webgpu/gpu-pca.js

@@ -0,0 +1,407 @@
+/**
+ * WebGPU-accelerated PCA with a numerically-equivalent CPU fallback.
+ *
+ * The GPU takes over exactly the large-GEMM hotspots — the XᵀX Gram product
+ * of covariance_eigh, the power-iteration panel products of the randomized
+ * solver (X stays resident on the device across iterations), and large
+ * transform projections. Everything else (eigh, panel QR/LU, the small SVD,
+ * sign flips, variance bookkeeping, validation) runs through the same code
+ * as the CPU `PCA` class via internal bridges, so semantics cannot drift.
+ *
+ * When WebGPU is unavailable — or the fit is too small to benefit, or the
+ * solver is inherently sequential (full, arpack) — fit() delegates to the
+ * CPU implementation and produces bit-identical results. When the GPU is
+ * used, results match the CPU path within the df64 GEMM accuracy measured
+ * by the equivalence test suite (~1e-12 relative on components).
+ */
+import { castTo } from '../base.js';
+import { asMatrix, Matrix } from '../matrix.js';
+import { randomizedSvdSteps } from '../numeric/randomized.js';
+import { checkRandomState } from '../numeric/rng.js';
+import { centerInPlace, colMeans } from '../numeric/stats.js';
+import { PCA, resolveSvdSolver, validateNcForSolver, } from '../pca.js';
+import { makeReporter } from '../progress.js';
+import { throwIfAborted } from '../scheduling.js';
+import { assertAllFinite, assertMinSamplesForFit, checkFeatureCount } from '../validation.js';
+import { GpuEngine, isWebGPUSupported } from './engine.js';
+/** Wraps a GPU failure that must not trigger a CPU refit (input was mutated). */
+class GpuFitUnrecoverableError extends Error {
+    cause;
+    constructor(cause) {
+        super('WebGPU fit failed after in-place mutation');
+        this.cause = cause;
+    }
+}
+export class WebGPUPCA {
+    cpu;
+    deviceOptions;
+    minGpuElements;
+    enginePromise = null;
+    engine = null;
+    backend_ = 'cpu';
+    gpuFitting = false;
+    constructor(options = {}) {
+        const { device, powerPreference, minGpuElements, ...pcaOptions } = options;
+        this.cpu = new PCA(pcaOptions);
+        this.deviceOptions = { device, powerPreference };
+        this.minGpuElements = minGpuElements ?? 1 << 18;
+    }
+    /** True when a WebGPU implementation is present in this environment. */
+    static isSupported() {
+        return isWebGPUSupported();
+    }
+    /** 'webgpu' when the last fit/transform used the GPU, else 'cpu'. */
+    get backend() {
+        return this.backend_;
+    }
+    /** Adapter description of the engine in use, or null on the CPU fallback. */
+    get gpuAdapterInfo() {
+        return this.engine?.adapterInfo ?? null;
+    }
+    /** Measured df64 GEMM relative error from the device self-check (NaN on CPU). */
+    get gpuMeasuredRelError() {
+        return this.engine?.measuredRelError ?? Number.NaN;
+    }
+    /** Releases GPU resources (does not destroy an injected device). */
+    dispose() {
+        this.engine?.dispose();
+        this.engine = null;
+        this.enginePromise = null;
+    }
+    // ------------------------------------------------------------------
+    // Fitting
+    // ------------------------------------------------------------------
+    async fit(X, options = {}) {
+        const xm = asMatrix(X);
+        this.guardConcurrentFit();
+        try {
+            const route = await this.route(xm);
+            if (route === 'cpu') {
+                // The fallback is non-blocking, observable, and abortable too —
+                // and still bit-identical to the plain PCA class.
+                await this.cpu.fitAsync(xm, options);
+                this.backend_ = 'cpu';
+                return this;
+            }
+            try {
+                if (route === 'covariance_eigh') {
+                    await this.fitGramGpu(xm, options);
+                }
+                else {
+                    await this.fitRandomizedGpu(xm, options);
+                }
+            }
+            catch (err) {
+                await this.recoverOnGpuError(xm, err, options);
+            }
+            return this;
+        }
+        finally {
+            this.gpuFitting = false;
+        }
+    }
+    async fitTransform(X, options = {}) {
+        const xm = asMatrix(X);
+        this.guardConcurrentFit();
+        try {
+            const route = await this.route(xm);
+            if (route === 'cpu') {
+                this.backend_ = 'cpu';
+                return await this.cpu.fitTransformAsync(xm, options);
+            }
+            if (route === 'covariance_eigh') {
+                try {
+                    await this.fitGramGpu(xm, options);
+                }
+                catch (err) {
+                    await this.recoverOnGpuError(xm, err, options);
+                }
+                // sklearn's fit_transform for covariance_eigh IS transform-after-fit
+                // (no U is computed at fit time), so this matches exactly.
+                return await this.transform(xm);
+            }
+            let dec;
+            try {
+                dec = await this.fitRandomizedGpu(xm, options);
+            }
+            catch (err) {
+                await this.recoverOnGpuError(xm, err, options, 'skipFit');
+                return await this.cpu.fitTransformAsync(xm, options);
+            }
+            // Randomized fast path: X_new = U·S (or U·√(n−1) when whitening),
+            // exactly like PCA.fitTransform.
+            const n = xm.rows;
+            const k = this.cpu.nComponents;
+            const out = new Float64Array(n * k);
+            const whiten = this.cpu.whiten;
+            const f = Math.sqrt(n - 1);
+            for (let i = 0; i < n; i++) {
+                for (let c = 0; c < k; c++) {
+                    out[i * k + c] = whiten ? dec.u[i * k + c] * f : dec.u[i * k + c] * dec.s[c];
+                }
+            }
+            return new Matrix(castTo(out, xm.dtype), n, k);
+        }
+        finally {
+            this.gpuFitting = false;
+        }
+    }
+    /** Concurrent fits on one instance would interleave GPU state; throw early. */
+    guardConcurrentFit() {
+        if (this.gpuFitting) {
+            throw new Error('This WebGPUPCA instance is already fitting; concurrent fits on one instance are not supported');
+        }
+        this.gpuFitting = true;
+    }
+    /**
+     * Validates the input (before any GPU spend or in-place mutation) and
+     * decides which backend the fit runs on. Mutates nothing.
+     */
+    async route(xm) {
+        if (xm.rows < 1 || xm.cols < 1) {
+            throw new Error(`Found array with shape (${xm.rows}, ${xm.cols}); PCA requires at least 1 sample and 1 feature`);
+        }
+        assertMinSamplesForFit(xm, 'PCA');
+        assertAllFinite(xm, 'PCA.fit');
+        const o = this.cpu._resolvedOptions();
+        const solver = resolveSvdSolver(xm.rows, xm.cols, o.nComponents, o.svdSolver);
+        const engine = await this.ensureEngine();
+        if (!this.gpuEligible(engine, xm, solver)) {
+            return 'cpu';
+        }
+        // Fail doomed configurations before uploading anything.
+        const nc = o.nComponents === null ? Math.min(xm.rows, xm.cols) : o.nComponents;
+        validateNcForSolver(nc, xm.rows, xm.cols, solver);
+        return solver;
+    }
+    /** GPU Gram product + shared CPU covariance_eigh tail. Mutates nothing on failure. */
+    async fitGramGpu(xm, options) {
+        const eng = this.engine;
+        const reporter = makeReporter(options, {
+            estimator: 'PCA',
+            solver: 'covariance_eigh',
+            nRows: xm.rows,
+            whiten: this.cpu.whiten,
+        });
+        throwIfAborted(options.signal);
+        const gx = eng.upload(xm.data, xm.rows, xm.cols);
+        try {
+            reporter?.emit({ phase: 'gram', step: 0, totalSteps: 1 });
+            const gram = await eng.syrk(gx);
+            reporter?.emit({ phase: 'gram', step: 1, totalSteps: 1 });
+            throwIfAborted(options.signal);
+            reporter?.emit({ phase: 'decompose', step: 0, totalSteps: null });
+            this.cpu._fitGram(xm, gram);
+            this.backend_ = 'webgpu';
+            if (reporter) {
+                this.cpu._emitFinalize(reporter, xm, null, 0, null);
+            }
+        }
+        finally {
+            gx.destroy();
+        }
+    }
+    /**
+     * GPU-driven randomized fit: centers on the CPU (respecting copy
+     * semantics), keeps X resident on the device, and drives the shared
+     * algorithm generator with GPU GEMMs. Returns the (sign-flipped)
+     * decomposition for the fitTransform fast path.
+     */
+    async fitRandomizedGpu(xm, options) {
+        const eng = this.engine;
+        const o = this.cpu._resolvedOptions();
+        const n = xm.rows;
+        const p = xm.cols;
+        const nc = (o.nComponents === null ? Math.min(n, p) : o.nComponents);
+        const reporter = makeReporter(options, {
+            estimator: 'PCA',
+            solver: 'randomized',
+            nRows: n,
+            whiten: this.cpu.whiten,
+        });
+        throwIfAborted(options.signal);
+        const meanF64 = colMeans(xm.data, n, p);
+        const xc = o.copy ? xm.data.slice() : xm.data;
+        centerInPlace(xc, n, p, meanF64);
+        const gx = eng.upload(xc, n, p);
+        try {
+            const rng = checkRandomState(o.randomState);
+            // The reporter rides the solver's hooks: power-iteration events and
+            // snapshot decompositions surface as ordinary GEMM requests below, so
+            // X stays device-resident throughout.
+            const gen = randomizedSvdSteps(n, p, nc, {
+                nOversamples: o.nOversamples,
+                nIter: o.iteratedPower,
+                powerIterationNormalizer: o.powerIterationNormalizer,
+                rng,
+                float32Stream: xm.dtype === 'float32',
+                hooks: reporter,
+            });
+            let step = gen.next();
+            while (!step.done) {
+                throwIfAborted(options.signal);
+                const req = step.value;
+                const result = req.op === 'mulA'
+                    ? await eng.mulA(gx, req.b, req.w)
+                    : req.op === 'mulAT'
+                        ? await eng.mulAT(gx, req.b, req.w)
+                        : await eng.mulTA(gx, req.b, req.w);
+                step = gen.next(result);
+            }
+            const dec = step.value;
+            this.cpu._fitDecomposed(xm, xc, meanF64, dec);
+            this.backend_ = 'webgpu';
+            if (reporter) {
+                this.cpu._emitFinalize(reporter, xm, dec.u, nc, dec.s);
+            }
+            return dec;
+        }
+        catch (err) {
+            if (!o.copy) {
+                // X was centered in place; a silent CPU retry would re-center the
+                // mutated buffer and corrupt the fit. Surface the GPU failure.
+                throw new GpuFitUnrecoverableError(err);
+            }
+            throw err;
+        }
+        finally {
+            gx.destroy();
+        }
+    }
+    /**
+     * CPU refit after a GPU failure — unless the input was already mutated,
+     * or the "failure" is a cancellation: an abort must reject with the abort
+     * error, never silently refit on the CPU.
+     */
+    async recoverOnGpuError(xm, err, options, mode) {
+        if (err instanceof GpuFitUnrecoverableError) {
+            throw err.cause;
+        }
+        if (options.signal?.aborted || err?.name === 'AbortError') {
+            throw err;
+        }
+        this.backend_ = 'cpu';
+        if (mode !== 'skipFit') {
+            await this.cpu.fitAsync(xm, options);
+        }
+    }
+    // ------------------------------------------------------------------
+    // Transforms (GPU projection for large inputs, CPU otherwise)
+    // ------------------------------------------------------------------
+    async transform(X) {
+        const xm = asMatrix(X);
+        const engine = this.engine;
+        const comp = this.cpu.components; // asserts fitted
+        if (engine === null ||
+            xm.rows * xm.cols < this.minGpuElements ||
+            !engine.canFit(xm.rows, xm.cols) ||
+            !engine.canFit(comp.rows, comp.cols)) {
+            return this.cpu.transform(xm);
+        }
+        assertAllFinite(xm, 'PCA.transform');
+        checkFeatureCount(xm, this.cpu.nFeaturesIn, 'PCA');
+        const gx = engine.upload(xm.data, xm.rows, xm.cols);
+        const gc = engine.upload(comp.data, comp.rows, comp.cols);
+        try {
+            const proj = await engine.project(gx, gc);
+            return this.cpu._transformFromProjection(proj, xm.rows, xm.dtype, false);
+        }
+        catch {
+            return this.cpu.transform(xm);
+        }
+        finally {
+            gx.destroy();
+            gc.destroy();
+        }
+    }
+    /** CPU delegate (small k×p GEMM; not worth a device round-trip). */
+    inverseTransform(X) {
+        return this.cpu.inverseTransform(X);
+    }
+    // ------------------------------------------------------------------
+    // Fitted state and CPU-delegated methods
+    // ------------------------------------------------------------------
+    get components() {
+        return this.cpu.components;
+    }
+    get explainedVariance() {
+        return this.cpu.explainedVariance;
+    }
+    get explainedVarianceRatio() {
+        return this.cpu.explainedVarianceRatio;
+    }
+    get singularValues() {
+        return this.cpu.singularValues;
+    }
+    get mean() {
+        return this.cpu.mean;
+    }
+    get nComponents() {
+        return this.cpu.nComponents;
+    }
+    get nSamples() {
+        return this.cpu.nSamples;
+    }
+    get nFeaturesIn() {
+        return this.cpu.nFeaturesIn;
+    }
+    get noiseVariance() {
+        return this.cpu.noiseVariance;
+    }
+    get whiten() {
+        return this.cpu.whiten;
+    }
+    get resolvedSvdSolver() {
+        return this.cpu.resolvedSvdSolver;
+    }
+    getCovariance() {
+        return this.cpu.getCovariance();
+    }
+    getPrecision() {
+        return this.cpu.getPrecision();
+    }
+    scoreSamples(X) {
+        return this.cpu.scoreSamples(X);
+    }
+    score(X) {
+        return this.cpu.score(X);
+    }
+    getFeatureNamesOut() {
+        return this.cpu.getFeatureNamesOut();
+    }
+    // ------------------------------------------------------------------
+    // Model serialization
+    // ------------------------------------------------------------------
+    /** The fitted model as a plain object — identical to the CPU class's. */
+    toModel() {
+        return this.cpu.toModel();
+    }
+    /**
+     * Rehydrates a fitted WebGPUPCA. `options` supplies the device knobs
+     * (powerPreference, injected device, minGpuElements); the PCA options
+     * come from the model itself.
+     */
+    static fromModel(model, options = {}) {
+        const gpu = new WebGPUPCA({ ...model.options, ...options });
+        gpu.cpu = PCA.fromModel(model);
+        return gpu;
+    }
+    // ------------------------------------------------------------------
+    async ensureEngine() {
+        if (this.enginePromise === null) {
+            this.enginePromise = GpuEngine.create(this.deviceOptions).then((e) => {
+                this.engine = e;
+                return e;
+            });
+        }
+        return this.enginePromise;
+    }
+    gpuEligible(engine, xm, solver) {
+        return (engine !== null &&
+            (solver === 'covariance_eigh' || solver === 'randomized') &&
+            xm.rows * xm.cols >= this.minGpuElements &&
+            engine.canFit(xm.rows, xm.cols) &&
+            engine.canFit(xm.cols, xm.cols));
+    }
+}
+//# sourceMappingURL=gpu-pca.js.map

package/dist/webgpu/gpu-pca.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"gpu-pca.js","sourceRoot":"","sources":["../../src/webgpu/gpu-pca.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AACH,OAAO,EAAE,MAAM,EAAE,MAAM,YAAY,CAAC;AACpC,OAAO,EAAE,QAAQ,EAAE,MAAM,EAAoB,MAAM,cAAc,CAAC;AAElE,OAAO,EAAE,kBAAkB,EAAE,MAAM,0BAA0B,CAAC;AAC9D,OAAO,EAAE,gBAAgB,EAAE,MAAM,mBAAmB,CAAC;AACrD,OAAO,EAAE,aAAa,EAAE,QAAQ,EAAE,MAAM,qBAAqB,CAAC;AAE9D,OAAO,EACL,GAAG,EAEH,gBAAgB,EAEhB,mBAAmB,GACpB,MAAM,WAAW,CAAC;AACnB,OAAO,EAAwB,YAAY,EAAE,MAAM,gBAAgB,CAAC;AACpE,OAAO,EAAE,cAAc,EAAE,MAAM,kBAAkB,CAAC;AAElD,OAAO,EAAE,eAAe,EAAE,sBAAsB,EAAE,iBAAiB,EAAE,MAAM,kBAAkB,CAAC;AAC9F,OAAO,EAAE,SAAS,EAAE,iBAAiB,EAA4B,MAAM,aAAa,CAAC;AAWrF,iFAAiF;AACjF,MAAM,wBAAyB,SAAQ,KAAK;IACxB,KAAK,CAAU;IACjC,YAAY,KAAc;QACxB,KAAK,CAAC,2CAA2C,CAAC,CAAC;QACnD,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC;IACrB,CAAC;CACF;AAED,MAAM,OAAO,SAAS;IACZ,GAAG,CAAM;IACA,aAAa,CAAsB;IACnC,cAAc,CAAS;IAChC,aAAa,GAAqC,IAAI,CAAC;IACvD,MAAM,GAAqB,IAAI,CAAC;IAChC,QAAQ,GAAqB,KAAK,CAAC;IACnC,UAAU,GAAG,KAAK,CAAC;IAE3B,YAAY,UAA4B,EAAE;QACxC,MAAM,EAAE,MAAM,EAAE,eAAe,EAAE,cAAc,EAAE,GAAG,UAAU,EAAE,GAAG,OAAO,CAAC;QAC3E,IAAI,CAAC,GAAG,GAAG,IAAI,GAAG,CAAC,UAAU,CAAC,CAAC;QAC/B,IAAI,CAAC,aAAa,GAAG,EAAE,MAAM,EAAE,eAAe,EAAE,CAAC;QACjD,IAAI,CAAC,cAAc,GAAG,cAAc,IAAI,CAAC,IAAI,EAAE,CAAC;IAClD,CAAC;IAED,wEAAwE;IACxE,MAAM,CAAC,WAAW;QAChB,OAAO,iBAAiB,EAAE,CAAC;IAC7B,CAAC;IAED,qEAAqE;IACrE,IAAI,OAAO;QACT,OAAO,IAAI,CAAC,QAAQ,CAAC;IACvB,CAAC;IAED,6EAA6E;IAC7E,IAAI,cAAc;QAChB,OAAO,IAAI,CAAC,MAAM,EAAE,WAAW,IAAI,IAAI,CAAC;IAC1C,CAAC;IAED,iFAAiF;IACjF,IAAI,mBAAmB;QACrB,OAAO,IAAI,CAAC,MAAM,EAAE,gBAAgB,IAAI,MAAM,CAAC,GAAG,CAAC;IACrD,CAAC;IAED,oEAAoE;IACpE,OAAO;QACL,IAAI,CAAC,MAAM,EAAE,OAAO,EAAE,CAAC;QACvB,IAAI,CAAC,MAAM,GAAG,IAAI,CAAC;QACnB,IAAI,CAAC,aAAa,GAAG,IAAI,CAAC;IAC5B,CAAC;IAED,qEAAqE;IACrE,UAAU;IACV,qEAAqE;IAErE,KAAK,CAAC,GAAG,CAAC,CAAc,EAAE,UAA2B,EAAE;QACrD,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;QACvB,IAAI,CAAC,kBAAkB,EAAE,CAAC;QAC1B,IAAI,CAAC;YACH,MAAM,KAAK,GAAG,MAAM,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC;YACnC,IAAI,KAAK,KAAK,KAAK,EAAE,CAAC;gBACpB,gEAAgE;gBAChE,kDAAkD;gBAClD,MAAM,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;gBACrC,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC;gBACtB,OAAO,IAAI,CAAC;YACd,CAAC;YACD,IAAI,CAAC;gBACH,IAAI,KAAK,KAAK,iBAAiB,EAAE,CAAC;oBAChC,MAAM,IAAI,CAAC,UAAU,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;gBACrC,CAAC;qBAAM,CAAC;oBACN,MAAM,IAAI,CAAC,gBAAgB,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;gBAC3C,CAAC;YACH,CAAC;YAAC,OAAO,GAAG,EAAE,CAAC;gBACb,MAAM,IAAI,CAAC,iBAAiB,CAAC,EAAE,EAAE,GAAG,EAAE,OAAO,CAAC,CAAC;YACjD,CAAC;YACD,OAAO,IAAI,CAAC;QACd,CAAC;gBAAS,CAAC;YACT,IAAI,CAAC,UAAU,GAAG,KAAK,CAAC;QAC1B,CAAC;IACH,CAAC;IAED,KAAK,CAAC,YAAY,CAAC,CAAc,EAAE,UAA2B,EAAE;QAC9D,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;QACvB,IAAI,CAAC,kBAAkB,EAAE,CAAC;QAC1B,IAAI,CAAC;YACH,MAAM,KAAK,GAAG,MAAM,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC;YACnC,IAAI,KAAK,KAAK,KAAK,EAAE,CAAC;gBACpB,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC;gBACtB,OAAO,MAAM,IAAI,CAAC,GAAG,CAAC,iBAAiB,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;YACvD,CAAC;YACD,IAAI,KAAK,KAAK,iBAAiB,EAAE,CAAC;gBAChC,IAAI,CAAC;oBACH,MAAM,IAAI,CAAC,UAAU,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;gBACrC,CAAC;gBAAC,OAAO,GAAG,EAAE,CAAC;oBACb,MAAM,IAAI,CAAC,iBAAiB,CAAC,EAAE,EAAE,GAAG,EAAE,OAAO,CAAC,CAAC;gBACjD,CAAC;gBACD,qEAAqE;gBACrE,2DAA2D;gBAC3D,OAAO,MAAM,IAAI,CAAC,SAAS,CAAC,EAAE,CAAC,CAAC;YAClC,CAAC;YACD,IAAI,GAAc,CAAC;YACnB,IAAI,CAAC;gBACH,GAAG,GAAG,MAAM,IAAI,CAAC,gBAAgB,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;YACjD,CAAC;YAAC,OAAO,GAAG,EAAE,CAAC;gBACb,MAAM,IAAI,CAAC,iBAAiB,CAAC,EAAE,EAAE,GAAG,EAAE,OAAO,EAAE,SAAS,CAAC,CAAC;gBAC1D,OAAO,MAAM,IAAI,CAAC,GAAG,CAAC,iBAAiB,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;YACvD,CAAC;YACD,kEAAkE;YAClE,iCAAiC;YACjC,MAAM,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC;YAClB,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,WAAW,CAAC;YAC/B,MAAM,GAAG,GAAG,IAAI,YAAY,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACpC,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC;YAC/B,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC3B,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC/E,CAAC;YACH,CAAC;YACD,OAAO,IAAI,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QACjD,CAAC;gBAAS,CAAC;YACT,IAAI,CAAC,UAAU,GAAG,KAAK,CAAC;QAC1B,CAAC;IACH,CAAC;IAED,+EAA+E;IACvE,kBAAkB;QACxB,IAAI,IAAI,CAAC,UAAU,EAAE,CAAC;YACpB,MAAM,IAAI,KAAK,CACb,+FAA+F,CAChG,CAAC;QACJ,CAAC;QACD,IAAI,CAAC,UAAU,GAAG,IAAI,CAAC;IACzB,CAAC;IAED;;;OAGG;IACK,KAAK,CAAC,KAAK,CAAC,EAAU;QAC5B,IAAI,EAAE,CAAC,IAAI,GAAG,CAAC,IAAI,EAAE,CAAC,IAAI,GAAG,CAAC,EAAE,CAAC;YAC/B,MAAM,IAAI,KAAK,CACb,2BAA2B,EAAE,CAAC,IAAI,KAAK,EAAE,CAAC,IAAI,iDAAiD,CAChG,CAAC;QACJ,CAAC;QACD,sBAAsB,CAAC,EAAE,EAAE,KAAK,CAAC,CAAC;QAClC,eAAe,CAAC,EAAE,EAAE,SAAS,CAAC,CAAC;QAC/B,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,gBAAgB,EAAE,CAAC;QACtC,MAAM,MAAM,GAAG,gBAAgB,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,IAAI,EAAE,CAAC,CAAC,WAAW,EAAE,CAAC,CAAC,SAAS,CAAC,CAAC;QAC9E,MAAM,MAAM,GAAG,MAAM,IAAI,CAAC,YAAY,EAAE,CAAC;QACzC,IAAI,CAAC,IAAI,CAAC,WAAW,CAAC,MAAM,EAAE,EAAE,EAAE,MAAM,CAAC,EAAE,CAAC;YAC1C,OAAO,KAAK,CAAC;QACf,CAAC;QACD,wDAAwD;QACxD,MAAM,EAAE,GAAG,CAAC,CAAC,WAAW,KAAK,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,WAAW,CAAC;QAC/E,mBAAmB,CAAC,EAAE,EAAE,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,IAAI,EAAE,MAA0C,CAAC,CAAC;QACtF,OAAO,MAA0C,CAAC;IACpD,CAAC;IAED,sFAAsF;IAC9E,KAAK,CAAC,UAAU,CAAC,EAAU,EAAE,OAAwB;QAC3D,MAAM,GAAG,GAAG,IAAI,CAAC,MAAmB,CAAC;QACrC,MAAM,QAAQ,GAAG,YAAY,CAAC,OAAO,EAAE;YACrC,SAAS,EAAE,KAAK;YAChB,MAAM,EAAE,iBAAiB;YACzB,KAAK,EAAE,EAAE,CAAC,IAAI;YACd,MAAM,EAAE,IAAI,CAAC,GAAG,CAAC,MAAM;SACxB,CAAC,CAAC;QACH,cAAc,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;QAC/B,MAAM,EAAE,GAAG,GAAG,CAAC,MAAM,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,IAAI,CAAC,CAAC;QACjD,IAAI,CAAC;YACH,QAAQ,EAAE,IAAI,CAAC,EAAE,KAAK,EAAE,MAAM,EAAE,IAAI,EAAE,CAAC,EAAE,UAAU,EAAE,CAAC,EAAE,CAAC,CAAC;YAC1D,MAAM,IAAI,GAAG,MAAM,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;YAChC,QAAQ,EAAE,IAAI,CAAC,EAAE,KAAK,EAAE,MAAM,EAAE,IAAI,EAAE,CAAC,EAAE,UAAU,EAAE,CAAC,EAAE,CAAC,CAAC;YAC1D,cAAc,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;YAC/B,QAAQ,EAAE,IAAI,CAAC,EAAE,KAAK,EAAE,WAAW,EAAE,IAAI,EAAE,CAAC,EAAE,UAAU,EAAE,IAAI,EAAE,CAAC,CAAC;YAClE,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,EAAE,IAAI,CAAC,CAAC;YAC5B,IAAI,CAAC,QAAQ,GAAG,QAAQ,CAAC;YACzB,IAAI,QAAQ,EAAE,CAAC;gBACb,IAAI,CAAC,GAAG,CAAC,aAAa,CAAC,QAAQ,EAAE,EAAE,EAAE,IAAI,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC;YACtD,CAAC;QACH,CAAC;gBAAS,CAAC;YACT,EAAE,CAAC,OAAO,EAAE,CAAC;QACf,CAAC;IACH,CAAC;IAED;;;;;OAKG;IACK,KAAK,CAAC,gBAAgB,CAAC,EAAU,EAAE,OAAwB;QACjE,MAAM,GAAG,GAAG,IAAI,CAAC,MAAmB,CAAC;QACrC,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,gBAAgB,EAAE,CAAC;QACtC,MAAM,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC;QAClB,MAAM,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC;QAClB,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,WAAW,KAAK,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,WAAW,CAAW,CAAC;QAC/E,MAAM,QAAQ,GAAG,YAAY,CAAC,OAAO,EAAE;YACrC,SAAS,EAAE,KAAK;YAChB,MAAM,EAAE,YAAY;YACpB,KAAK,EAAE,CAAC;YACR,MAAM,EAAE,IAAI,CAAC,GAAG,CAAC,MAAM;SACxB,CAAC,CAAC;QACH,cAAc,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;QAC/B,MAAM,OAAO,GAAG,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QACxC,MAAM,EAAE,GAAe,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC;QAC1D,aAAa,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC;QACjC,MAAM,EAAE,GAAG,GAAG,CAAC,MAAM,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAChC,IAAI,CAAC;YACH,MAAM,GAAG,GAAG,gBAAgB,CAAC,CAAC,CAAC,WAAW,CAAC,CAAC;YAC5C,oEAAoE;YACpE,sEAAsE;YACtE,sCAAsC;YACtC,MAAM,GAAG,GAAG,kBAAkB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE;gBACvC,YAAY,EAAE,CAAC,CAAC,YAAY;gBAC5B,KAAK,EAAE,CAAC,CAAC,aAAa;gBACtB,wBAAwB,EAAE,CAAC,CAAC,wBAAwB;gBACpD,GAAG;gBACH,aAAa,EAAE,EAAE,CAAC,KAAK,KAAK,SAAS;gBACrC,KAAK,EAAE,QAAQ;aAChB,CAAC,CAAC;YACH,IAAI,IAAI,GAAG,GAAG,CAAC,IAAI,EAAE,CAAC;YACtB,OAAO,CAAC,IAAI,CAAC,IAAI,EAAE,CAAC;gBAClB,cAAc,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;gBAC/B,MAAM,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC;gBACvB,MAAM,MAAM,GACV,GAAG,CAAC,EAAE,KAAK,MAAM;oBACf,CAAC,CAAC,MAAM,GAAG,CAAC,IAAI,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;oBAClC,CAAC,CAAC,GAAG,CAAC,EAAE,KAAK,OAAO;wBAClB,CAAC,CAAC,MAAM,GAAG,CAAC,KAAK,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;wBACnC,CAAC,CAAC,MAAM,GAAG,CAAC,KAAK,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;gBAC1C,IAAI,GAAG,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;YAC1B,CAAC;YACD,MAAM,GAAG,GAAc,IAAI,CAAC,KAAK,CAAC;YAClC,IAAI,CAAC,GAAG,CAAC,cAAc,CAAC,EAAE,EAAE,EAAE,EAAE,OAAO,EAAE,GAAG,CAAC,CAAC;YAC9C,IAAI,CAAC,QAAQ,GAAG,QAAQ,CAAC;YACzB,IAAI,QAAQ,EAAE,CAAC;gBACb,IAAI,CAAC,GAAG,CAAC,aAAa,CAAC,QAAQ,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;YACzD,CAAC;YACD,OAAO,GAAG,CAAC;QACb,CAAC;QAAC,OAAO,GAAG,EAAE,CAAC;YACb,IAAI,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC;gBACZ,kEAAkE;gBAClE,+DAA+D;gBAC/D,MAAM,IAAI,wBAAwB,CAAC,GAAG,CAAC,CAAC;YAC1C,CAAC;YACD,MAAM,GAAG,CAAC;QACZ,CAAC;gBAAS,CAAC;YACT,EAAE,CAAC,OAAO,EAAE,CAAC;QACf,CAAC;IACH,CAAC;IAED;;;;OAIG;IACK,KAAK,CAAC,iBAAiB,CAC7B,EAAU,EACV,GAAY,EACZ,OAAwB,EACxB,IAAgB;QAEhB,IAAI,GAAG,YAAY,wBAAwB,EAAE,CAAC;YAC5C,MAAM,GAAG,CAAC,KAAK,CAAC;QAClB,CAAC;QACD,IAAI,OAAO,CAAC,MAAM,EAAE,OAAO,IAAK,GAAiC,EAAE,IAAI,KAAK,YAAY,EAAE,CAAC;YACzF,MAAM,GAAG,CAAC;QACZ,CAAC;QACD,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC;QACtB,IAAI,IAAI,KAAK,SAAS,EAAE,CAAC;YACvB,MAAM,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;QACvC,CAAC;IACH,CAAC;IAED,qEAAqE;IACrE,8DAA8D;IAC9D,qEAAqE;IAErE,KAAK,CAAC,SAAS,CAAC,CAAc;QAC5B,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;QACvB,MAAM,MAAM,GAAG,IAAI,CAAC,MAAM,CAAC;QAC3B,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC,iBAAiB;QACnD,IACE,MAAM,KAAK,IAAI;YACf,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,IAAI,GAAG,IAAI,CAAC,cAAc;YACvC,CAAC,MAAM,CAAC,MAAM,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,IAAI,CAAC;YAChC,CAAC,MAAM,CAAC,MAAM,CAAC,IAAI,CAAC,IAAI,EAAE,IAAI,CAAC,IAAI,CAAC,EACpC,CAAC;YACD,OAAO,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,EAAE,CAAC,CAAC;QAChC,CAAC;QACD,eAAe,CAAC,EAAE,EAAE,eAAe,CAAC,CAAC;QACrC,iBAAiB,CAAC,EAAE,EAAE,IAAI,CAAC,GAAG,CAAC,WAAW,EAAE,KAAK,CAAC,CAAC;QACnD,MAAM,EAAE,GAAG,MAAM,CAAC,MAAM,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,IAAI,CAAC,CAAC;QACpD,MAAM,EAAE,GAAG,MAAM,CAAC,MAAM,CAAC,IAAI,CAAC,IAAI,EAAE,IAAI,CAAC,IAAI,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC;QAC1D,IAAI,CAAC;YACH,MAAM,IAAI,GAAG,MAAM,MAAM,CAAC,OAAO,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;YAC1C,OAAO,IAAI,CAAC,GAAG,CAAC,wBAAwB,CAAC,IAAI,EAAE,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;QAC3E,CAAC;QAAC,MAAM,CAAC;YACP,OAAO,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,EAAE,CAAC,CAAC;QAChC,CAAC;gBAAS,CAAC;YACT,EAAE,CAAC,OAAO,EAAE,CAAC;YACb,EAAE,CAAC,OAAO,EAAE,CAAC;QACf,CAAC;IACH,CAAC;IAED,oEAAoE;IACpE,gBAAgB,CAAC,CAAc;QAC7B,OAAO,IAAI,CAAC,GAAG,CAAC,gBAAgB,CAAC,CAAC,CAAC,CAAC;IACtC,CAAC;IAED,qEAAqE;IACrE,yCAAyC;IACzC,qEAAqE;IAErE,IAAI,UAAU;QACZ,OAAO,IAAI,CAAC,GAAG,CAAC,UAAU,CAAC;IAC7B,CAAC;IACD,IAAI,iBAAiB;QACnB,OAAO,IAAI,CAAC,GAAG,CAAC,iBAAiB,CAAC;IACpC,CAAC;IACD,IAAI,sBAAsB;QACxB,OAAO,IAAI,CAAC,GAAG,CAAC,sBAAsB,CAAC;IACzC,CAAC;IACD,IAAI,cAAc;QAChB,OAAO,IAAI,CAAC,GAAG,CAAC,cAAc,CAAC;IACjC,CAAC;IACD,IAAI,IAAI;QACN,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC;IACvB,CAAC;IACD,IAAI,WAAW;QACb,OAAO,IAAI,CAAC,GAAG,CAAC,WAAW,CAAC;IAC9B,CAAC;IACD,IAAI,QAAQ;QACV,OAAO,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC;IAC3B,CAAC;IACD,IAAI,WAAW;QACb,OAAO,IAAI,CAAC,GAAG,CAAC,WAAW,CAAC;IAC9B,CAAC;IACD,IAAI,aAAa;QACf,OAAO,IAAI,CAAC,GAAG,CAAC,aAAa,CAAC;IAChC,CAAC;IACD,IAAI,MAAM;QACR,OAAO,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC;IACzB,CAAC;IACD,IAAI,iBAAiB;QACnB,OAAO,IAAI,CAAC,GAAG,CAAC,iBAAiB,CAAC;IACpC,CAAC;IAED,aAAa;QACX,OAAO,IAAI,CAAC,GAAG,CAAC,aAAa,EAAE,CAAC;IAClC,CAAC;IACD,YAAY;QACV,OAAO,IAAI,CAAC,GAAG,CAAC,YAAY,EAAE,CAAC;IACjC,CAAC;IACD,YAAY,CAAC,CAAc;QACzB,OAAO,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC,CAAC,CAAC;IAClC,CAAC;IACD,KAAK,CAAC,CAAc;QAClB,OAAO,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IAC3B,CAAC;IACD,kBAAkB;QAChB,OAAO,IAAI,CAAC,GAAG,CAAC,kBAAkB,EAAE,CAAC;IACvC,CAAC;IAED,qEAAqE;IACrE,sBAAsB;IACtB,qEAAqE;IAErE,yEAAyE;IACzE,OAAO;QACL,OAAO,IAAI,CAAC,GAAG,CAAC,OAAO,EAAE,CAAC;IAC5B,CAAC;IAED;;;;OAIG;IACH,MAAM,CAAC,SAAS,CACd,KAAe,EACf,UAA6D,EAAE;QAE/D,MAAM,GAAG,GAAG,IAAI,SAAS,CAAC,EAAE,GAAG,KAAK,CAAC,OAAO,EAAE,GAAG,OAAO,EAAE,CAAC,CAAC;QAC5D,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC;QAC/B,OAAO,GAAG,CAAC;IACb,CAAC;IAED,qEAAqE;IAE7D,KAAK,CAAC,YAAY;QACxB,IAAI,IAAI,CAAC,aAAa,KAAK,IAAI,EAAE,CAAC;YAChC,IAAI,CAAC,aAAa,GAAG,SAAS,CAAC,MAAM,CAAC,IAAI,CAAC,aAAa,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE;gBACnE,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC;gBAChB,OAAO,CAAC,CAAC;YACX,CAAC,CAAC,CAAC;QACL,CAAC;QACD,OAAO,IAAI,CAAC,aAAa,CAAC;IAC5B,CAAC;IAEO,WAAW,CAAC,MAAwB,EAAE,EAAU,EAAE,MAAc;QACtE,OAAO,CACL,MAAM,KAAK,IAAI;YACf,CAAC,MAAM,KAAK,iBAAiB,IAAI,MAAM,KAAK,YAAY,CAAC;YACzD,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,IAAI,IAAI,IAAI,CAAC,cAAc;YACxC,MAAM,CAAC,MAAM,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,IAAI,CAAC;YAC/B,MAAM,CAAC,MAAM,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,IAAI,CAAC,CAChC,CAAC;IACJ,CAAC;CACF"}

package/dist/webgpu/index.d.ts

@@ -0,0 +1,19 @@
+/**
+ * `pca-web/webgpu` — WebGPU-accelerated PCA.
+ *
+ * import { WebGPUPCA, isWebGPUSupported } from 'pca-web/webgpu';
+ *
+ * TypeScript: the declarations reference the standard WebGPU globals.
+ * With `lib: ["dom"]` (TS ≥ 5.9) they resolve natively; Node-only
+ * consumers add the types-only package @webgpu/types instead (see README).
+ */
+export { asMatrix, Matrix, type MatrixInput } from '../matrix.js';
+export type { PCAModel, PCAModelOptions } from '../model.js';
+export { assertValidModel, modelFromJSON, modelToJSON } from '../model.js';
+export { PCA, type PCAOptions, type PowerIterationNormalizer, type SvdSolver } from '../pca.js';
+export type { FitAsyncOptions, FitObserver, FitPhase, FitSolverId, PCAFitProgress, PCAFitSnapshot, SnapshotOptions, } from '../progress.js';
+export type { AbortSignalLike } from '../scheduling.js';
+export { FitAbortError } from '../scheduling.js';
+export { GpuEngine, isWebGPUSupported, type WebGPUDeviceOptions } from './engine.js';
+export { WebGPUPCA, type WebGPUPCAOptions } from './gpu-pca.js';
+//# sourceMappingURL=index.d.ts.map

package/dist/webgpu/index.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../src/webgpu/index.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AAEH,OAAO,EAAE,QAAQ,EAAE,MAAM,EAAE,KAAK,WAAW,EAAE,MAAM,cAAc,CAAC;AAClE,YAAY,EAAE,QAAQ,EAAE,eAAe,EAAE,MAAM,aAAa,CAAC;AAC7D,OAAO,EAAE,gBAAgB,EAAE,aAAa,EAAE,WAAW,EAAE,MAAM,aAAa,CAAC;AAC3E,OAAO,EAAE,GAAG,EAAE,KAAK,UAAU,EAAE,KAAK,wBAAwB,EAAE,KAAK,SAAS,EAAE,MAAM,WAAW,CAAC;AAChG,YAAY,EACV,eAAe,EACf,WAAW,EACX,QAAQ,EACR,WAAW,EACX,cAAc,EACd,cAAc,EACd,eAAe,GAChB,MAAM,gBAAgB,CAAC;AACxB,YAAY,EAAE,eAAe,EAAE,MAAM,kBAAkB,CAAC;AACxD,OAAO,EAAE,aAAa,EAAE,MAAM,kBAAkB,CAAC;AACjD,OAAO,EAAE,SAAS,EAAE,iBAAiB,EAAE,KAAK,mBAAmB,EAAE,MAAM,aAAa,CAAC;AACrF,OAAO,EAAE,SAAS,EAAE,KAAK,gBAAgB,EAAE,MAAM,cAAc,CAAC"}

package/dist/webgpu/index.js

@@ -0,0 +1,16 @@
+/**
+ * `pca-web/webgpu` — WebGPU-accelerated PCA.
+ *
+ * import { WebGPUPCA, isWebGPUSupported } from 'pca-web/webgpu';
+ *
+ * TypeScript: the declarations reference the standard WebGPU globals.
+ * With `lib: ["dom"]` (TS ≥ 5.9) they resolve natively; Node-only
+ * consumers add the types-only package @webgpu/types instead (see README).
+ */
+export { asMatrix, Matrix } from '../matrix.js';
+export { assertValidModel, modelFromJSON, modelToJSON } from '../model.js';
+export { PCA } from '../pca.js';
+export { FitAbortError } from '../scheduling.js';
+export { GpuEngine, isWebGPUSupported } from './engine.js';
+export { WebGPUPCA } from './gpu-pca.js';
+//# sourceMappingURL=index.js.map

package/dist/webgpu/index.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/webgpu/index.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AAEH,OAAO,EAAE,QAAQ,EAAE,MAAM,EAAoB,MAAM,cAAc,CAAC;AAElE,OAAO,EAAE,gBAAgB,EAAE,aAAa,EAAE,WAAW,EAAE,MAAM,aAAa,CAAC;AAC3E,OAAO,EAAE,GAAG,EAAkE,MAAM,WAAW,CAAC;AAWhG,OAAO,EAAE,aAAa,EAAE,MAAM,kBAAkB,CAAC;AACjD,OAAO,EAAE,SAAS,EAAE,iBAAiB,EAA4B,MAAM,aAAa,CAAC;AACrF,OAAO,EAAE,SAAS,EAAyB,MAAM,cAAc,CAAC"}

package/dist/webgpu/kernels.d.ts

@@ -0,0 +1,47 @@
+/**
+ * WGSL GEMM kernels with near-f64 accuracy on f32-only hardware.
+ *
+ * Design constraint discovered empirically (Chrome/Dawn on Metal): the MSL
+ * compiler applies fast-math — classic Dekker/Knuth error-free
+ * transformations (two_sum, split) are algebraically folded to zero, and
+ * multiply-add pairs are contracted into fma, so double-single (df64)
+ * accumulation silently degrades to f32. What DOES survive, verified by
+ * on-device probes:
+ *   - fma(a, b, -a*b) computes the exact f32 multiplication error
+ *   - single float ops are correctly rounded (no cross-op guarantees)
+ *   - i32 conversions and integer adds are exact and untouched by fast-math
+ *
+ * So products use fma-exact two_prod, and ACCUMULATION is exact integer
+ * binning: every product term is decomposed — exactly, via power-of-two
+ * scaling and fma residuals — into 5 levels of 13-bit integers on a fixed
+ * per-dispatch grid (exponent E bounds every |a·b| term), accumulated in
+ * i32 (no rounding at all), and recombined in float64 on the CPU. The only
+ * error is the truncated tail below the last grid level:
+ *   |error| ≤ K·7·2^(E-65) ≤ maxProduct·2^-45
+ * comparable to true f64 accumulation. Inputs are (hi, lo) f32 pairs of the
+ * f64 values; the three significant cross products (+ the tiny lo·lo term)
+ * feed the binner, preserving ~48 significant bits end to end.
+ *
+ * One tiled kernel template covers all four op(A)·op(B) transpose variants;
+ * only the injected index expressions differ.
+ */
+export declare const TILE = 16;
+/**
+ * Per-level bin capacity is 2^13; with ≤7 terms per MAC each contributing
+ * |n| ≤ 2^13 per level, i32 overflows at K ≈ 2^31/(7·2^13) ≈ 37k. The
+ * engine splits the K dimension at this bound and sums chunks in f64.
+ */
+export declare const MAX_K_PER_DISPATCH = 30000;
+/** i32 bin levels per output element (13 bits each → 65 bits of grid). */
+export declare const BIN_LEVELS = 5;
+/** i32 slots per output element in the C buffer (BIN_LEVELS padded to 8). */
+export declare const OUT_SLOTS = 8;
+/**
+ * Builds the GEMM shader for C(M×N) = op(A)·op(B) over K = [kStart, kStart+kCount).
+ *
+ * Physical shapes: A is M×K_total when !ta, K_total×M when ta; B is
+ * K_total×N when !tb, N×K_total when tb. All row-major, ds-encoded
+ * (array<vec2f> of hi/lo pairs). C is array<i32>, OUT_SLOTS per element.
+ */
+export declare function gemmShader(ta: boolean, tb: boolean): string;
+//# sourceMappingURL=kernels.d.ts.map

package/dist/webgpu/kernels.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"kernels.d.ts","sourceRoot":"","sources":["../../src/webgpu/kernels.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;GA0BG;AAEH,eAAO,MAAM,IAAI,KAAK,CAAC;AAEvB;;;;GAIG;AACH,eAAO,MAAM,kBAAkB,QAAQ,CAAC;AAExC,0EAA0E;AAC1E,eAAO,MAAM,UAAU,IAAI,CAAC;AAE5B,6EAA6E;AAC7E,eAAO,MAAM,SAAS,IAAI,CAAC;AAE3B;;;;;;GAMG;AACH,wBAAgB,UAAU,CAAC,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,GAAG,MAAM,CAoG3D"}