pca-web 0.1.0

package/dist/incremental-pca.js

@@ -0,0 +1,397 @@
+/**
+ * Incremental PCA matching `sklearn.decomposition.IncrementalPCA` (1.9):
+ * constant-memory fitting via `partialFit` batches (Ross et al. 2008), so the
+ * dataset never needs to reside in memory at once.
+ */
+import { BasePCA, castTo } from './base.js';
+import { asMatrix, Matrix } from './matrix.js';
+import { assertValidModel } from './model.js';
+import { colMeans, incrementalMeanAndVar } from './numeric/stats.js';
+import { svd } from './numeric/svd.js';
+import { svdFlipVBased } from './numeric/svdflip.js';
+import { makeReporter, projectForSnapshot, toFloat64Copy, } from './progress.js';
+import { driveAsync, driveSync } from './scheduling.js';
+import { dtypeOf, epsFor } from './types.js';
+import { assertAllFinite, assertMinSamplesForFit, checkFeatureCount } from './validation.js';
+function validateOptions(o) {
+    const nComponents = o.nComponents ?? null;
+    if (nComponents !== null && !(Number.isInteger(nComponents) && nComponents >= 1)) {
+        throw new Error(`nComponents must be an integer >= 1 or null; got ${nComponents}`);
+    }
+    const batchSize = o.batchSize ?? null;
+    if (batchSize !== null && !(Number.isInteger(batchSize) && batchSize >= 1)) {
+        throw new Error(`batchSize must be an integer >= 1 or null; got ${batchSize}`);
+    }
+    return {
+        nComponents,
+        whiten: o.whiten ?? false,
+        copy: o.copy ?? true,
+        batchSize,
+    };
+}
+/** sklearn's `gen_batches`: full slices, with a too-small tail merged into the last one. */
+function* genBatches(n, batchSize, minBatchSize) {
+    let start = 0;
+    const nFull = Math.floor(n / batchSize);
+    for (let i = 0; i < nFull; i++) {
+        const end = start + batchSize;
+        if (end + minBatchSize > n) {
+            continue;
+        }
+        yield [start, end];
+        start = end;
+    }
+    if (start < n) {
+        yield [start, n];
+    }
+}
+export class IncrementalPCA extends BasePCA {
+    estimatorName = 'IncrementalPCA';
+    opts;
+    var_ = null;
+    meanF64_ = null;
+    nSamplesSeen_ = 0;
+    batchSize_ = null;
+    /** Mirrors sklearn's hasattr(self, 'components_') — set by fit() before batching. */
+    hasComponentsAttr_ = false;
+    constructor(options = {}) {
+        super();
+        this.opts = validateOptions(options);
+        this.whitenOpt = this.opts.whiten;
+    }
+    /** Per-feature running variance (always float64, like sklearn) — `var_`. */
+    get variance() {
+        this.assertFitted();
+        return this.var_;
+    }
+    /** Samples processed so far — sklearn's `n_samples_seen_`. */
+    get nSamplesSeen() {
+        this.assertFitted();
+        return this.nSamplesSeen_;
+    }
+    /** Batch size used by fit — sklearn's `batch_size_` (only set by fit()). */
+    get batchSize() {
+        if (this.batchSize_ === null) {
+            throw new Error("batchSize is only available after fit() (sklearn's batch_size_)");
+        }
+        return this.batchSize_;
+    }
+    // ------------------------------------------------------------------
+    // Fitting
+    // ------------------------------------------------------------------
+    /**
+     * Fit in minibatches of `batchSize` rows — sklearn's `IncrementalPCA.fit`.
+     * The first batch must have at least 2 samples.
+     */
+    fit(X, observer) {
+        driveSync(this._fitSteps(asMatrix(X), observer), observer?.signal);
+        return this;
+    }
+    /**
+     * Non-blocking fit: identical batching, time-sliced on the event loop.
+     * Results are bit-identical to the synchronous fit. An abort between
+     * batches keeps the completed-batch model (the estimator stays usable).
+     */
+    async fitAsync(X, options = {}) {
+        await driveAsync(this._fitSteps(asMatrix(X), options), {
+            budgetMs: options.budgetMs,
+            signal: options.signal,
+        });
+        return this;
+    }
+    /**
+     * @internal The batched fit as a step generator — one `yield` per batch
+     * is the suspension/abort checkpoint for the drivers. Mirrors sklearn's
+     * fit(): state is reset before validation, and each batch runs the
+     * unchanged partialFit update.
+     */
+    *_fitSteps(X, observer) {
+        if (this.fitting) {
+            throw new Error('This IncrementalPCA instance is already fitting; concurrent fits on one instance are not supported');
+        }
+        this.fitting = true;
+        try {
+            // Reset all state, mirroring sklearn's fit().
+            this.components_ = null;
+            this.hasComponentsAttr_ = true; // sklearn sets self.components_ = None
+            this.nSamplesSeen_ = 0;
+            this.meanF64_ = null;
+            this.var_ = null;
+            this.singularValues_ = null;
+            this.explainedVariance_ = null;
+            this.explainedVarianceRatio_ = null;
+            this.noiseVariance_ = 0;
+            this.fitted = false;
+            assertAllFinite(X, 'IncrementalPCA.fit');
+            if (X.rows < 1 || X.cols < 1) {
+                throw new Error(`Found array with shape (${X.rows}, ${X.cols}); at least 1 sample and 1 feature required`);
+            }
+            const work = this.opts.copy ? X.copy() : X;
+            const n = work.rows;
+            const p = work.cols;
+            this.nFeaturesIn_ = p;
+            this.batchSize_ = this.opts.batchSize ?? 5 * p;
+            // Materialized so totalSteps is known up front.
+            const batches = [...genBatches(n, this.batchSize_, this.opts.nComponents ?? 0)];
+            const reporter = makeReporter(observer, {
+                estimator: 'IncrementalPCA',
+                solver: 'incremental',
+                nRows: n,
+                whiten: this.opts.whiten,
+            });
+            for (let i = 0; i < batches.length; i++) {
+                const [start, end] = batches[i];
+                const rows = end - start;
+                const batch = new Matrix(work.data.subarray(start * p, end * p), rows, p);
+                this.partialFitCore(batch);
+                if (reporter) {
+                    reporter.emit({
+                        phase: 'batch',
+                        step: i + 1,
+                        totalSteps: batches.length,
+                        snapshot: reporter.wantSnapshot(i + 1)
+                            ? this.batchSnapshot(X, end, reporter.scoresRequested)
+                            : undefined,
+                    });
+                }
+                yield;
+            }
+            if (reporter) {
+                reporter.emit({
+                    phase: 'finalize',
+                    step: batches.length,
+                    totalSteps: batches.length,
+                    snapshot: reporter.snapshotsEnabled
+                        ? this.batchSnapshot(X, n, reporter.scoresRequested)
+                        : undefined,
+                });
+            }
+        }
+        finally {
+            this.fitting = false;
+        }
+    }
+    /**
+     * Snapshot of the just-stored (batch-prefix) model, as fresh float64
+     * copies. Scores project the original input's rows seen so far; with
+     * `copy: false` the first batch's rows were destructively centered by
+     * the update (sklearn's contract), making those scores approximate.
+     */
+    batchSnapshot(X, rowsSeen, wantScores) {
+        const comp = this.components_;
+        const k = comp.rows;
+        const p = comp.cols;
+        const snapshot = {
+            components: new Matrix(toFloat64Copy(comp.data), k, p),
+            singularValues: toFloat64Copy(this.singularValues_),
+            explainedVariance: toFloat64Copy(this.explainedVariance_),
+        };
+        if (wantScores) {
+            const ev = this.explainedVariance_;
+            snapshot.scores = projectForSnapshot(X.data, rowsSeen, p, comp.data, k, this.mean_, ev, this.whitenOpt, epsFor(dtypeOf(ev)));
+        }
+        return snapshot;
+    }
+    /**
+     * Incremental fit on one batch — sklearn's `partial_fit`. All of X is one
+     * batch. The first call must see at least 2 samples; later calls may
+     * stream batches of any size (including single rows).
+     */
+    partialFit(X) {
+        if (this.fitting) {
+            throw new Error('This IncrementalPCA instance is fitting; partialFit cannot run concurrently with fit');
+        }
+        const xm = asMatrix(X);
+        assertAllFinite(xm, 'IncrementalPCA.partialFit');
+        if (xm.rows < 1 || xm.cols < 1) {
+            throw new Error(`Found array with shape (${xm.rows}, ${xm.cols}); at least 1 sample and 1 feature required`);
+        }
+        if (this.nFeaturesIn_ > 0 && this.hasComponentsAttr_) {
+            checkFeatureCount(xm, this.nFeaturesIn_, this.estimatorName);
+        }
+        else {
+            this.nFeaturesIn_ = xm.cols;
+        }
+        const work = this.opts.copy ? xm.copy() : xm;
+        this.partialFitCore(work);
+        return this;
+    }
+    /** Equivalent to fit(X).transform(X) (sklearn's TransformerMixin behavior). */
+    fitTransform(X, observer) {
+        this.fit(X, observer);
+        return this.transform(X);
+    }
+    /** Non-blocking fitTransform — fitAsync followed by transform. */
+    async fitTransformAsync(X, options = {}) {
+        await this.fitAsync(X, options);
+        return this.transform(X);
+    }
+    // ------------------------------------------------------------------
+    // Model serialization
+    // ------------------------------------------------------------------
+    /**
+     * The fitted model as a plain structured-clone-friendly object. Includes
+     * the running statistics (`nSamplesSeen`, float64 mean/variance), so
+     * `partialFit` resumes bit-exactly after rehydration.
+     */
+    toModel() {
+        const base = this.exportBaseModel();
+        return {
+            ...base,
+            estimator: 'ipca',
+            nSamplesSeen: this.nSamplesSeen_,
+            variance: this.var_.slice(),
+            batchSize: this.batchSize_,
+            options: {
+                nComponents: this.opts.nComponents,
+                whiten: this.opts.whiten,
+                copy: this.opts.copy,
+                batchSize: this.opts.batchSize,
+            },
+        };
+    }
+    /**
+     * Rehydrates a fitted IncrementalPCA from `toModel()` output (validated
+     * first; arrays adopted without copying). Further `partialFit` calls
+     * continue the stream exactly where the exported estimator left off.
+     */
+    static fromModel(model) {
+        assertValidModel(model, 'ipca');
+        const ipca = new IncrementalPCA(model.options);
+        ipca.importBaseModel(model);
+        // mean_ and meanF64_ alias one array in live estimators; preserve that.
+        ipca.meanF64_ = model.mean;
+        ipca.var_ = model.variance;
+        ipca.nSamplesSeen_ = model.nSamplesSeen;
+        ipca.batchSize_ = model.batchSize;
+        ipca.hasComponentsAttr_ = true;
+        return ipca;
+    }
+    /** The whitening/update step; mutates `X`'s data (callers pass a copy unless copy=false). */
+    partialFitCore(X) {
+        // Only the very first batch of data needs 2+ samples (it seeds the
+        // variance estimate); later 1-row updates are well-defined.
+        if (this.nSamplesSeen_ === 0) {
+            assertMinSamplesForFit(X, this.estimatorName);
+        }
+        const firstPass = !this.hasComponentsAttr_;
+        if (firstPass) {
+            this.components_ = null;
+            this.hasComponentsAttr_ = true;
+        }
+        const n = X.rows;
+        const p = X.cols;
+        // Resolve n_components exactly like sklearn.
+        const ncOpt = this.opts.nComponents;
+        if (ncOpt === null) {
+            if (this.components_ === null) {
+                this.nComponents_ = Math.min(n, p);
+            }
+            else {
+                this.nComponents_ = this.components_.rows;
+            }
+        }
+        else if (!(ncOpt <= p)) {
+            throw new Error(`nComponents=${ncOpt} invalid for nFeatures=${p}, need more rows than columns for IncrementalPCA processing`);
+        }
+        else if (ncOpt > n && firstPass) {
+            throw new Error(`nComponents=${ncOpt} must be less or equal to the batch number of samples ${n} for the first partialFit call.`);
+        }
+        else {
+            this.nComponents_ = ncOpt;
+        }
+        if (this.components_ !== null && this.components_.rows !== this.nComponents_) {
+            throw new Error(`Number of components has changed from ${this.components_.rows} to ` +
+                `${this.nComponents_} between calls to partialFit! Try setting nComponents to a fixed value.`);
+        }
+        const k = this.nComponents_;
+        // Update running statistics (float64, like sklearn's float64 accumulators).
+        const stats = incrementalMeanAndVar(X.data, n, p, this.meanF64_, this.var_, this.nSamplesSeen_);
+        const nTotal = stats.count;
+        // Build the matrix to decompose.
+        let stacked;
+        let stackedRows;
+        let stackedDtype;
+        if (this.nSamplesSeen_ === 0) {
+            // First step: just center the batch (in its own dtype, like sklearn).
+            stackedDtype = X.dtype;
+            stackedRows = n;
+            stacked = X.data;
+            for (let i = 0; i < n; i++) {
+                const off = i * p;
+                for (let j = 0; j < p; j++) {
+                    stacked[off + j] -= stats.mean[j];
+                }
+            }
+        }
+        else {
+            // Later steps: numpy's vstack promotes to float64 because the
+            // mean-correction row is float64 — replicated (a float32 model's
+            // attributes really do become float64 after the second partial_fit).
+            const comp = this.components_;
+            const sv = this.singularValues_;
+            const colBatchMean = colMeans(X.data, n, p);
+            stackedDtype = 'float64';
+            stackedRows = k + n + 1;
+            const st = new Float64Array(stackedRows * p);
+            for (let c = 0; c < k; c++) {
+                const f = sv[c];
+                const off = c * p;
+                for (let j = 0; j < p; j++) {
+                    st[off + j] = f * comp.data[off + j];
+                }
+            }
+            // X -= col_batch_mean happens in the batch's dtype in sklearn.
+            const f32 = X.dtype === 'float32';
+            for (let i = 0; i < n; i++) {
+                const src = i * p;
+                const dst = (k + i) * p;
+                for (let j = 0; j < p; j++) {
+                    const v = X.data[src + j] - colBatchMean[j];
+                    st[dst + j] = f32 ? Math.fround(v) : v;
+                }
+            }
+            const corrScale = Math.sqrt((this.nSamplesSeen_ / nTotal) * n);
+            const off = (k + n) * p;
+            for (let j = 0; j < p; j++) {
+                st[off + j] = corrScale * (this.meanF64_[j] - colBatchMean[j]);
+            }
+            stacked = st;
+        }
+        const dec = svd(stacked, stackedRows, p);
+        const minDim = Math.min(stackedRows, p);
+        svdFlipVBased(dec.u, stackedRows, dec.vt, minDim, p);
+        const explainedVariance = new Float64Array(minDim);
+        let varTimesTotal = 0;
+        for (let j = 0; j < p; j++) {
+            varTimesTotal += stats.variance[j] * nTotal;
+        }
+        const ratio = new Float64Array(minDim);
+        for (let i = 0; i < minDim; i++) {
+            const s2 = dec.s[i] * dec.s[i];
+            explainedVariance[i] = s2 / (nTotal - 1);
+            ratio[i] = s2 / varTimesTotal;
+        }
+        this.nSamplesSeen_ = nTotal;
+        this.meanF64_ = stats.mean;
+        this.var_ = stats.variance;
+        this.dtype = stackedDtype;
+        this.components_ = new Matrix(castTo(dec.vt.slice(0, k * p), stackedDtype), k, p);
+        this.singularValues_ = castTo(dec.s.slice(0, k), stackedDtype);
+        this.explainedVariance_ = castTo(explainedVariance.slice(0, k), stackedDtype);
+        this.explainedVarianceRatio_ = castTo(ratio.slice(0, k), stackedDtype);
+        this.mean_ = stats.mean; // float64, like sklearn
+        if (k !== n && k !== p) {
+            let acc = 0;
+            for (let i = k; i < minDim; i++) {
+                acc += explainedVariance[i];
+            }
+            this.noiseVariance_ = acc / (minDim - k);
+        }
+        else {
+            this.noiseVariance_ = 0;
+        }
+        this.fitted = true;
+    }
+}
+//# sourceMappingURL=incremental-pca.js.map

package/dist/incremental-pca.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"incremental-pca.js","sourceRoot":"","sources":["../src/incremental-pca.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AACH,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,MAAM,WAAW,CAAC;AAC5C,OAAO,EAAE,QAAQ,EAAE,MAAM,EAAoB,MAAM,aAAa,CAAC;AACjE,OAAO,EAAE,gBAAgB,EAA4B,MAAM,YAAY,CAAC;AACxE,OAAO,EAAE,QAAQ,EAAE,qBAAqB,EAAE,MAAM,oBAAoB,CAAC;AACrE,OAAO,EAAE,GAAG,EAAE,MAAM,kBAAkB,CAAC;AACvC,OAAO,EAAE,aAAa,EAAE,MAAM,sBAAsB,CAAC;AACrD,OAAO,EAGL,YAAY,EAEZ,kBAAkB,EAClB,aAAa,GACd,MAAM,eAAe,CAAC;AACvB,OAAO,EAAE,UAAU,EAAE,SAAS,EAAE,MAAM,iBAAiB,CAAC;AACxD,OAAO,EAAc,OAAO,EAAE,MAAM,EAAmB,MAAM,YAAY,CAAC;AAC1E,OAAO,EAAE,eAAe,EAAE,sBAAsB,EAAE,iBAAiB,EAAE,MAAM,iBAAiB,CAAC;AAoB7F,SAAS,eAAe,CAAC,CAAwB;IAC/C,MAAM,WAAW,GAAG,CAAC,CAAC,WAAW,IAAI,IAAI,CAAC;IAC1C,IAAI,WAAW,KAAK,IAAI,IAAI,CAAC,CAAC,MAAM,CAAC,SAAS,CAAC,WAAW,CAAC,IAAI,WAAW,IAAI,CAAC,CAAC,EAAE,CAAC;QACjF,MAAM,IAAI,KAAK,CAAC,oDAAoD,WAAW,EAAE,CAAC,CAAC;IACrF,CAAC;IACD,MAAM,SAAS,GAAG,CAAC,CAAC,SAAS,IAAI,IAAI,CAAC;IACtC,IAAI,SAAS,KAAK,IAAI,IAAI,CAAC,CAAC,MAAM,CAAC,SAAS,CAAC,SAAS,CAAC,IAAI,SAAS,IAAI,CAAC,CAAC,EAAE,CAAC;QAC3E,MAAM,IAAI,KAAK,CAAC,kDAAkD,SAAS,EAAE,CAAC,CAAC;IACjF,CAAC;IACD,OAAO;QACL,WAAW;QACX,MAAM,EAAE,CAAC,CAAC,MAAM,IAAI,KAAK;QACzB,IAAI,EAAE,CAAC,CAAC,IAAI,IAAI,IAAI;QACpB,SAAS;KACV,CAAC;AACJ,CAAC;AAED,4FAA4F;AAC5F,QAAQ,CAAC,CAAC,UAAU,CAClB,CAAS,EACT,SAAiB,EACjB,YAAoB;IAEpB,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC;IACxC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,MAAM,GAAG,GAAG,KAAK,GAAG,SAAS,CAAC;QAC9B,IAAI,GAAG,GAAG,YAAY,GAAG,CAAC,EAAE,CAAC;YAC3B,SAAS;QACX,CAAC;QACD,MAAM,CAAC,KAAK,EAAE,GAAG,CAAC,CAAC;QACnB,KAAK,GAAG,GAAG,CAAC;IACd,CAAC;IACD,IAAI,KAAK,GAAG,CAAC,EAAE,CAAC;QACd,MAAM,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;IACnB,CAAC;AACH,CAAC;AAED,MAAM,OAAO,cAAe,SAAQ,OAAO;IACb,aAAa,GAAW,gBAAgB,CAAC;IACpD,IAAI,CAAkB;IAE/B,IAAI,GAAwB,IAAI,CAAC;IACjC,QAAQ,GAAwB,IAAI,CAAC;IACrC,aAAa,GAAG,CAAC,CAAC;IAClB,UAAU,GAAkB,IAAI,CAAC;IACzC,qFAAqF;IAC7E,kBAAkB,GAAG,KAAK,CAAC;IAEnC,YAAY,UAAiC,EAAE;QAC7C,KAAK,EAAE,CAAC;QACR,IAAI,CAAC,IAAI,GAAG,eAAe,CAAC,OAAO,CAAC,CAAC;QACrC,IAAI,CAAC,SAAS,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC;IACpC,CAAC;IAED,4EAA4E;IAC5E,IAAI,QAAQ;QACV,IAAI,CAAC,YAAY,EAAE,CAAC;QACpB,OAAO,IAAI,CAAC,IAAoB,CAAC;IACnC,CAAC;IAED,8DAA8D;IAC9D,IAAI,YAAY;QACd,IAAI,CAAC,YAAY,EAAE,CAAC;QACpB,OAAO,IAAI,CAAC,aAAa,CAAC;IAC5B,CAAC;IAED,4EAA4E;IAC5E,IAAI,SAAS;QACX,IAAI,IAAI,CAAC,UAAU,KAAK,IAAI,EAAE,CAAC;YAC7B,MAAM,IAAI,KAAK,CAAC,iEAAiE,CAAC,CAAC;QACrF,CAAC;QACD,OAAO,IAAI,CAAC,UAAU,CAAC;IACzB,CAAC;IAED,qEAAqE;IACrE,UAAU;IACV,qEAAqE;IAErE;;;OAGG;IACH,GAAG,CAAC,CAAc,EAAE,QAAsB;QACxC,SAAS,CAAC,IAAI,CAAC,SAAS,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,EAAE,QAAQ,EAAE,MAAM,CAAC,CAAC;QACnE,OAAO,IAAI,CAAC;IACd,CAAC;IAED;;;;OAIG;IACH,KAAK,CAAC,QAAQ,CAAC,CAAc,EAAE,UAA2B,EAAE;QAC1D,MAAM,UAAU,CAAC,IAAI,CAAC,SAAS,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,EAAE;YACrD,QAAQ,EAAE,OAAO,CAAC,QAAQ;YAC1B,MAAM,EAAE,OAAO,CAAC,MAAM;SACvB,CAAC,CAAC;QACH,OAAO,IAAI,CAAC;IACd,CAAC;IAED;;;;;OAKG;IACH,CAAC,SAAS,CAAC,CAAS,EAAE,QAAsB;QAC1C,IAAI,IAAI,CAAC,OAAO,EAAE,CAAC;YACjB,MAAM,IAAI,KAAK,CACb,oGAAoG,CACrG,CAAC;QACJ,CAAC;QACD,IAAI,CAAC,OAAO,GAAG,IAAI,CAAC;QACpB,IAAI,CAAC;YACH,8CAA8C;YAC9C,IAAI,CAAC,WAAW,GAAG,IAAI,CAAC;YACxB,IAAI,CAAC,kBAAkB,GAAG,IAAI,CAAC,CAAC,uCAAuC;YACvE,IAAI,CAAC,aAAa,GAAG,CAAC,CAAC;YACvB,IAAI,CAAC,QAAQ,GAAG,IAAI,CAAC;YACrB,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC;YACjB,IAAI,CAAC,eAAe,GAAG,IAAI,CAAC;YAC5B,IAAI,CAAC,kBAAkB,GAAG,IAAI,CAAC;YAC/B,IAAI,CAAC,uBAAuB,GAAG,IAAI,CAAC;YACpC,IAAI,CAAC,cAAc,GAAG,CAAC,CAAC;YACxB,IAAI,CAAC,MAAM,GAAG,KAAK,CAAC;YAEpB,eAAe,CAAC,CAAC,EAAE,oBAAoB,CAAC,CAAC;YACzC,IAAI,CAAC,CAAC,IAAI,GAAG,CAAC,IAAI,CAAC,CAAC,IAAI,GAAG,CAAC,EAAE,CAAC;gBAC7B,MAAM,IAAI,KAAK,CACb,2BAA2B,CAAC,CAAC,IAAI,KAAK,CAAC,CAAC,IAAI,6CAA6C,CAC1F,CAAC;YACJ,CAAC;YACD,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;YAC3C,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC;YACpB,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC;YACpB,IAAI,CAAC,YAAY,GAAG,CAAC,CAAC;YACtB,IAAI,CAAC,UAAU,GAAG,IAAI,CAAC,IAAI,CAAC,SAAS,IAAI,CAAC,GAAG,CAAC,CAAC;YAE/C,gDAAgD;YAChD,MAAM,OAAO,GAAG,CAAC,GAAG,UAAU,CAAC,CAAC,EAAE,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,IAAI,CAAC,WAAW,IAAI,CAAC,CAAC,CAAC,CAAC;YAChF,MAAM,QAAQ,GAAG,YAAY,CAAC,QAAQ,EAAE;gBACtC,SAAS,EAAE,gBAAgB;gBAC3B,MAAM,EAAE,aAAa;gBACrB,KAAK,EAAE,CAAC;gBACR,MAAM,EAAE,IAAI,CAAC,IAAI,CAAC,MAAM;aACzB,CAAC,CAAC;YACH,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;gBACxC,MAAM,CAAC,KAAK,EAAE,GAAG,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;gBAChC,MAAM,IAAI,GAAG,GAAG,GAAG,KAAK,CAAC;gBACzB,MAAM,KAAK,GAAG,IAAI,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,KAAK,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC;gBAC1E,IAAI,CAAC,cAAc,CAAC,KAAK,CAAC,CAAC;gBAC3B,IAAI,QAAQ,EAAE,CAAC;oBACb,QAAQ,CAAC,IAAI,CAAC;wBACZ,KAAK,EAAE,OAAO;wBACd,IAAI,EAAE,CAAC,GAAG,CAAC;wBACX,UAAU,EAAE,OAAO,CAAC,MAAM;wBAC1B,QAAQ,EAAE,QAAQ,CAAC,YAAY,CAAC,CAAC,GAAG,CAAC,CAAC;4BACpC,CAAC,CAAC,IAAI,CAAC,aAAa,CAAC,CAAC,EAAE,GAAG,EAAE,QAAQ,CAAC,eAAe,CAAC;4BACtD,CAAC,CAAC,SAAS;qBACd,CAAC,CAAC;gBACL,CAAC;gBACD,KAAK,CAAC;YACR,CAAC;YACD,IAAI,QAAQ,EAAE,CAAC;gBACb,QAAQ,CAAC,IAAI,CAAC;oBACZ,KAAK,EAAE,UAAU;oBACjB,IAAI,EAAE,OAAO,CAAC,MAAM;oBACpB,UAAU,EAAE,OAAO,CAAC,MAAM;oBAC1B,QAAQ,EAAE,QAAQ,CAAC,gBAAgB;wBACjC,CAAC,CAAC,IAAI,CAAC,aAAa,CAAC,CAAC,EAAE,CAAC,EAAE,QAAQ,CAAC,eAAe,CAAC;wBACpD,CAAC,CAAC,SAAS;iBACd,CAAC,CAAC;YACL,CAAC;QACH,CAAC;gBAAS,CAAC;YACT,IAAI,CAAC,OAAO,GAAG,KAAK,CAAC;QACvB,CAAC;IACH,CAAC;IAED;;;;;OAKG;IACK,aAAa,CAAC,CAAS,EAAE,QAAgB,EAAE,UAAmB;QACpE,MAAM,IAAI,GAAG,IAAI,CAAC,WAAqB,CAAC;QACxC,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC;QACpB,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC;QACpB,MAAM,QAAQ,GAAmB;YAC/B,UAAU,EAAE,IAAI,MAAM,CAAC,aAAa,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;YACtD,cAAc,EAAE,aAAa,CAAC,IAAI,CAAC,eAA6B,CAAC;YACjE,iBAAiB,EAAE,aAAa,CAAC,IAAI,CAAC,kBAAgC,CAAC;SACxE,CAAC;QACF,IAAI,UAAU,EAAE,CAAC;YACf,MAAM,EAAE,GAAG,IAAI,CAAC,kBAAgC,CAAC;YACjD,QAAQ,CAAC,MAAM,GAAG,kBAAkB,CAClC,CAAC,CAAC,IAAI,EACN,QAAQ,EACR,CAAC,EACD,IAAI,CAAC,IAAI,EACT,CAAC,EACD,IAAI,CAAC,KAAmB,EACxB,EAAE,EACF,IAAI,CAAC,SAAS,EACd,MAAM,CAAC,OAAO,CAAC,EAAE,CAAC,CAAC,CACpB,CAAC;QACJ,CAAC;QACD,OAAO,QAAQ,CAAC;IAClB,CAAC;IAED;;;;OAIG;IACH,UAAU,CAAC,CAAc;QACvB,IAAI,IAAI,CAAC,OAAO,EAAE,CAAC;YACjB,MAAM,IAAI,KAAK,CACb,sFAAsF,CACvF,CAAC;QACJ,CAAC;QACD,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;QACvB,eAAe,CAAC,EAAE,EAAE,2BAA2B,CAAC,CAAC;QACjD,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,6CAA6C,CAC5F,CAAC;QACJ,CAAC;QACD,IAAI,IAAI,CAAC,YAAY,GAAG,CAAC,IAAI,IAAI,CAAC,kBAAkB,EAAE,CAAC;YACrD,iBAAiB,CAAC,EAAE,EAAE,IAAI,CAAC,YAAY,EAAE,IAAI,CAAC,aAAa,CAAC,CAAC;QAC/D,CAAC;aAAM,CAAC;YACN,IAAI,CAAC,YAAY,GAAG,EAAE,CAAC,IAAI,CAAC;QAC9B,CAAC;QACD,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;QAC7C,IAAI,CAAC,cAAc,CAAC,IAAI,CAAC,CAAC;QAC1B,OAAO,IAAI,CAAC;IACd,CAAC;IAED,+EAA+E;IAC/E,YAAY,CAAC,CAAc,EAAE,QAAsB;QACjD,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC;QACtB,OAAO,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;IAC3B,CAAC;IAED,kEAAkE;IAClE,KAAK,CAAC,iBAAiB,CAAC,CAAc,EAAE,UAA2B,EAAE;QACnE,MAAM,IAAI,CAAC,QAAQ,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;QAChC,OAAO,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;IAC3B,CAAC;IAED,qEAAqE;IACrE,sBAAsB;IACtB,qEAAqE;IAErE;;;;OAIG;IACH,OAAO;QACL,MAAM,IAAI,GAAG,IAAI,CAAC,eAAe,EAAE,CAAC;QACpC,OAAO;YACL,GAAG,IAAI;YACP,SAAS,EAAE,MAAM;YACjB,YAAY,EAAE,IAAI,CAAC,aAAa;YAChC,QAAQ,EAAG,IAAI,CAAC,IAAqB,CAAC,KAAK,EAAE;YAC7C,SAAS,EAAE,IAAI,CAAC,UAAU;YAC1B,OAAO,EAAE;gBACP,WAAW,EAAE,IAAI,CAAC,IAAI,CAAC,WAAW;gBAClC,MAAM,EAAE,IAAI,CAAC,IAAI,CAAC,MAAM;gBACxB,IAAI,EAAE,IAAI,CAAC,IAAI,CAAC,IAAI;gBACpB,SAAS,EAAE,IAAI,CAAC,IAAI,CAAC,SAAS;aAC/B;SACF,CAAC;IACJ,CAAC;IAED;;;;OAIG;IACH,MAAM,CAAC,SAAS,CAAC,KAA0B;QACzC,gBAAgB,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC;QAChC,MAAM,IAAI,GAAG,IAAI,cAAc,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC;QAC/C,IAAI,CAAC,eAAe,CAAC,KAAK,CAAC,CAAC;QAC5B,wEAAwE;QACxE,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC,IAAoB,CAAC;QAC3C,IAAI,CAAC,IAAI,GAAG,KAAK,CAAC,QAAQ,CAAC;QAC3B,IAAI,CAAC,aAAa,GAAG,KAAK,CAAC,YAAY,CAAC;QACxC,IAAI,CAAC,UAAU,GAAG,KAAK,CAAC,SAAS,CAAC;QAClC,IAAI,CAAC,kBAAkB,GAAG,IAAI,CAAC;QAC/B,OAAO,IAAI,CAAC;IACd,CAAC;IAED,6FAA6F;IACrF,cAAc,CAAC,CAAS;QAC9B,mEAAmE;QACnE,4DAA4D;QAC5D,IAAI,IAAI,CAAC,aAAa,KAAK,CAAC,EAAE,CAAC;YAC7B,sBAAsB,CAAC,CAAC,EAAE,IAAI,CAAC,aAAa,CAAC,CAAC;QAChD,CAAC;QACD,MAAM,SAAS,GAAG,CAAC,IAAI,CAAC,kBAAkB,CAAC;QAC3C,IAAI,SAAS,EAAE,CAAC;YACd,IAAI,CAAC,WAAW,GAAG,IAAI,CAAC;YACxB,IAAI,CAAC,kBAAkB,GAAG,IAAI,CAAC;QACjC,CAAC;QACD,MAAM,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC;QACjB,MAAM,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC;QAEjB,6CAA6C;QAC7C,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,WAAW,CAAC;QACpC,IAAI,KAAK,KAAK,IAAI,EAAE,CAAC;YACnB,IAAI,IAAI,CAAC,WAAW,KAAK,IAAI,EAAE,CAAC;gBAC9B,IAAI,CAAC,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YACrC,CAAC;iBAAM,CAAC;gBACN,IAAI,CAAC,YAAY,GAAG,IAAI,CAAC,WAAW,CAAC,IAAI,CAAC;YAC5C,CAAC;QACH,CAAC;aAAM,IAAI,CAAC,CAAC,KAAK,IAAI,CAAC,CAAC,EAAE,CAAC;YACzB,MAAM,IAAI,KAAK,CACb,eAAe,KAAK,0BAA0B,CAAC,6DAA6D,CAC7G,CAAC;QACJ,CAAC;aAAM,IAAI,KAAK,GAAG,CAAC,IAAI,SAAS,EAAE,CAAC;YAClC,MAAM,IAAI,KAAK,CACb,eAAe,KAAK,yDAAyD,CAAC,iCAAiC,CAChH,CAAC;QACJ,CAAC;aAAM,CAAC;YACN,IAAI,CAAC,YAAY,GAAG,KAAK,CAAC;QAC5B,CAAC;QACD,IAAI,IAAI,CAAC,WAAW,KAAK,IAAI,IAAI,IAAI,CAAC,WAAW,CAAC,IAAI,KAAK,IAAI,CAAC,YAAY,EAAE,CAAC;YAC7E,MAAM,IAAI,KAAK,CACb,yCAAyC,IAAI,CAAC,WAAW,CAAC,IAAI,MAAM;gBAClE,GAAG,IAAI,CAAC,YAAY,yEAAyE,CAChG,CAAC;QACJ,CAAC;QACD,MAAM,CAAC,GAAG,IAAI,CAAC,YAAY,CAAC;QAE5B,4EAA4E;QAC5E,MAAM,KAAK,GAAG,qBAAqB,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,IAAI,EAAE,IAAI,CAAC,aAAa,CAAC,CAAC;QAChG,MAAM,MAAM,GAAG,KAAK,CAAC,KAAK,CAAC;QAE3B,iCAAiC;QACjC,IAAI,OAAmB,CAAC;QACxB,IAAI,WAAmB,CAAC;QACxB,IAAI,YAAmB,CAAC;QACxB,IAAI,IAAI,CAAC,aAAa,KAAK,CAAC,EAAE,CAAC;YAC7B,sEAAsE;YACtE,YAAY,GAAG,CAAC,CAAC,KAAK,CAAC;YACvB,WAAW,GAAG,CAAC,CAAC;YAChB,OAAO,GAAG,CAAC,CAAC,IAAI,CAAC;YACjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;gBAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC3B,OAAO,CAAC,GAAG,GAAG,CAAC,CAAC,IAAI,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;gBACpC,CAAC;YACH,CAAC;QACH,CAAC;aAAM,CAAC;YACN,8DAA8D;YAC9D,iEAAiE;YACjE,qEAAqE;YACrE,MAAM,IAAI,GAAG,IAAI,CAAC,WAAqB,CAAC;YACxC,MAAM,EAAE,GAAG,IAAI,CAAC,eAA6B,CAAC;YAC9C,MAAM,YAAY,GAAG,QAAQ,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YAC5C,YAAY,GAAG,SAAS,CAAC;YACzB,WAAW,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YACxB,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,WAAW,GAAG,CAAC,CAAC,CAAC;YAC7C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;gBAChB,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;gBAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC3B,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;gBACvC,CAAC;YACH,CAAC;YACD,+DAA+D;YAC/D,MAAM,GAAG,GAAG,CAAC,CAAC,KAAK,KAAK,SAAS,CAAC;YAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;gBAClB,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;gBACxB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC3B,MAAM,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,CAAC;oBAC5C,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACzC,CAAC;YACH,CAAC;YACD,MAAM,SAAS,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,aAAa,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC;YAC/D,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;YACxB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,SAAS,GAAG,CAAE,IAAI,CAAC,QAAyB,CAAC,CAAC,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC;YACnF,CAAC;YACD,OAAO,GAAG,EAAE,CAAC;QACf,CAAC;QAED,MAAM,GAAG,GAAG,GAAG,CAAC,OAAO,EAAE,WAAW,EAAE,CAAC,CAAC,CAAC;QACzC,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,WAAW,EAAE,CAAC,CAAC,CAAC;QACxC,aAAa,CAAC,GAAG,CAAC,CAAC,EAAE,WAAW,EAAE,GAAG,CAAC,EAAE,EAAE,MAAM,EAAE,CAAC,CAAC,CAAC;QAErD,MAAM,iBAAiB,GAAG,IAAI,YAAY,CAAC,MAAM,CAAC,CAAC;QACnD,IAAI,aAAa,GAAG,CAAC,CAAC;QACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,aAAa,IAAI,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;QAC9C,CAAC;QACD,MAAM,KAAK,GAAG,IAAI,YAAY,CAAC,MAAM,CAAC,CAAC;QACvC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YAChC,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC/B,iBAAiB,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;YACzC,KAAK,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,aAAa,CAAC;QAChC,CAAC;QAED,IAAI,CAAC,aAAa,GAAG,MAAM,CAAC;QAC5B,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC,IAAI,CAAC;QAC3B,IAAI,CAAC,IAAI,GAAG,KAAK,CAAC,QAAQ,CAAC;QAC3B,IAAI,CAAC,KAAK,GAAG,YAAY,CAAC;QAC1B,IAAI,CAAC,WAAW,GAAG,IAAI,MAAM,CAAC,MAAM,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,YAAY,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAClF,IAAI,CAAC,eAAe,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,YAAY,CAAC,CAAC;QAC/D,IAAI,CAAC,kBAAkB,GAAG,MAAM,CAAC,iBAAiB,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,YAAY,CAAC,CAAC;QAC9E,IAAI,CAAC,uBAAuB,GAAG,MAAM,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,YAAY,CAAC,CAAC;QACvE,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,wBAAwB;QACjD,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;YACvB,IAAI,GAAG,GAAG,CAAC,CAAC;YACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;gBAChC,GAAG,IAAI,iBAAiB,CAAC,CAAC,CAAC,CAAC;YAC9B,CAAC;YACD,IAAI,CAAC,cAAc,GAAG,GAAG,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;QAC3C,CAAC;aAAM,CAAC;YACN,IAAI,CAAC,cAAc,GAAG,CAAC,CAAC;QAC1B,CAAC;QACD,IAAI,CAAC,MAAM,GAAG,IAAI,CAAC;IACrB,CAAC;CACF"}

package/dist/index.d.ts

@@ -0,0 +1,21 @@
+/**
+ * pca-web — scikit-learn-compatible PCA and IncrementalPCA for the browser
+ * and Node.js, with zero runtime dependencies.
+ *
+ * WebGPU acceleration lives behind the `pca-web/webgpu` subpath export so
+ * the core stays free of GPU code for bundlers.
+ */
+export type { IncrementalPCAOptions } from './incremental-pca.js';
+export { IncrementalPCA } from './incremental-pca.js';
+export type { MatrixInput } from './matrix.js';
+export { asMatrix, Matrix } from './matrix.js';
+export type { AnyPCAModel, IncrementalPCAModel, IncrementalPCAModelOptions, PCAModel, PCAModelBaseFields, PCAModelOptions, } from './model.js';
+export { assertValidModel, modelFromJSON, modelToJSON, PCA_MODEL_FORMAT_VERSION, } from './model.js';
+export type { PCAOptions, PowerIterationNormalizer, SvdSolver } from './pca.js';
+export { PCA, RandomState } 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 type { Dtype, FloatArray } from './types.js';
+export { NotFittedError } from './validation.js';
+//# sourceMappingURL=index.d.ts.map

package/dist/index.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AAEH,YAAY,EAAE,qBAAqB,EAAE,MAAM,sBAAsB,CAAC;AAClE,OAAO,EAAE,cAAc,EAAE,MAAM,sBAAsB,CAAC;AACtD,YAAY,EAAE,WAAW,EAAE,MAAM,aAAa,CAAC;AAC/C,OAAO,EAAE,QAAQ,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AAC/C,YAAY,EACV,WAAW,EACX,mBAAmB,EACnB,0BAA0B,EAC1B,QAAQ,EACR,kBAAkB,EAClB,eAAe,GAChB,MAAM,YAAY,CAAC;AACpB,OAAO,EACL,gBAAgB,EAChB,aAAa,EACb,WAAW,EACX,wBAAwB,GACzB,MAAM,YAAY,CAAC;AACpB,YAAY,EAAE,UAAU,EAAE,wBAAwB,EAAE,SAAS,EAAE,MAAM,UAAU,CAAC;AAChF,OAAO,EAAE,GAAG,EAAE,WAAW,EAAE,MAAM,UAAU,CAAC;AAC5C,YAAY,EACV,eAAe,EACf,WAAW,EACX,QAAQ,EACR,WAAW,EACX,cAAc,EACd,cAAc,EACd,eAAe,GAChB,MAAM,eAAe,CAAC;AACvB,YAAY,EAAE,eAAe,EAAE,MAAM,iBAAiB,CAAC;AACvD,OAAO,EAAE,aAAa,EAAE,MAAM,iBAAiB,CAAC;AAChD,YAAY,EAAE,KAAK,EAAE,UAAU,EAAE,MAAM,YAAY,CAAC;AACpD,OAAO,EAAE,cAAc,EAAE,MAAM,iBAAiB,CAAC"}

package/dist/index.js

@@ -0,0 +1,14 @@
+/**
+ * pca-web — scikit-learn-compatible PCA and IncrementalPCA for the browser
+ * and Node.js, with zero runtime dependencies.
+ *
+ * WebGPU acceleration lives behind the `pca-web/webgpu` subpath export so
+ * the core stays free of GPU code for bundlers.
+ */
+export { IncrementalPCA } from './incremental-pca.js';
+export { asMatrix, Matrix } from './matrix.js';
+export { assertValidModel, modelFromJSON, modelToJSON, PCA_MODEL_FORMAT_VERSION, } from './model.js';
+export { PCA, RandomState } from './pca.js';
+export { FitAbortError } from './scheduling.js';
+export { NotFittedError } from './validation.js';
+//# sourceMappingURL=index.js.map

package/dist/index.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AAGH,OAAO,EAAE,cAAc,EAAE,MAAM,sBAAsB,CAAC;AAEtD,OAAO,EAAE,QAAQ,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AAS/C,OAAO,EACL,gBAAgB,EAChB,aAAa,EACb,WAAW,EACX,wBAAwB,GACzB,MAAM,YAAY,CAAC;AAEpB,OAAO,EAAE,GAAG,EAAE,WAAW,EAAE,MAAM,UAAU,CAAC;AAW5C,OAAO,EAAE,aAAa,EAAE,MAAM,iBAAiB,CAAC;AAEhD,OAAO,EAAE,cAAc,EAAE,MAAM,iBAAiB,CAAC"}

package/dist/matrix.d.ts

@@ -0,0 +1,44 @@
+import { type Dtype, type FloatArray } from './types.js';
+/**
+ * A dense row-major matrix backed by a Float64Array or Float32Array.
+ *
+ * This is the container returned by all pca-web methods. It is a thin,
+ * zero-copy wrapper: `data[i * cols + j]` is the element at row `i`,
+ * column `j`.
+ */
+export declare class Matrix {
+    readonly data: FloatArray;
+    readonly rows: number;
+    readonly cols: number;
+    /**
+     * Wraps an existing typed array (no copy). `data.length` must equal
+     * `rows * cols`.
+     */
+    constructor(data: FloatArray, rows: number, cols: number);
+    /** Storage dtype: 'float64' or 'float32'. */
+    get dtype(): Dtype;
+    /** Allocates a zero-initialized matrix. */
+    static zeros(rows: number, cols: number, dtype?: Dtype): Matrix;
+    /**
+     * Builds a matrix from an array of rows, e.g. `Matrix.from2d([[1, 2], [3, 4]])`.
+     * All rows must have equal length.
+     */
+    static from2d(rows2d: ArrayLike<number>[], dtype?: Dtype): Matrix;
+    get(i: number, j: number): number;
+    set(i: number, j: number, v: number): void;
+    /** Returns a deep copy (same dtype). */
+    copy(): Matrix;
+    /** Converts to a plain nested array of numbers. */
+    toArray(): number[][];
+    /** Copy of row `i` as a plain typed array (same dtype). */
+    row(i: number): FloatArray;
+}
+/** Accepted input type for training/transform data. */
+export type MatrixInput = Matrix | ArrayLike<number>[];
+/**
+ * Coerces user input to a Matrix without copying when it is already one.
+ * Nested arrays become float64 (mirrors sklearn coercing non-float dtypes
+ * to float64).
+ */
+export declare function asMatrix(X: MatrixInput): Matrix;
+//# sourceMappingURL=matrix.d.ts.map

package/dist/matrix.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"matrix.d.ts","sourceRoot":"","sources":["../src/matrix.ts"],"names":[],"mappings":"AAAA,OAAO,EAAkB,KAAK,KAAK,EAAW,KAAK,UAAU,EAAE,MAAM,YAAY,CAAC;AAElF;;;;;;GAMG;AACH,qBAAa,MAAM;IACjB,QAAQ,CAAC,IAAI,EAAE,UAAU,CAAC;IAC1B,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;IACtB,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;IAEtB;;;OAGG;gBACS,IAAI,EAAE,UAAU,EAAE,IAAI,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM;IAWxD,6CAA6C;IAC7C,IAAI,KAAK,IAAI,KAAK,CAEjB;IAED,2CAA2C;IAC3C,MAAM,CAAC,KAAK,CAAC,IAAI,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,KAAK,GAAE,KAAiB,GAAG,MAAM;IAK1E;;;OAGG;IACH,MAAM,CAAC,MAAM,CAAC,MAAM,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,EAAE,KAAK,GAAE,KAAiB,GAAG,MAAM;IAsB5E,GAAG,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM;IAIjC,GAAG,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,IAAI;IAI1C,wCAAwC;IACxC,IAAI,IAAI,MAAM;IAId,mDAAmD;IACnD,OAAO,IAAI,MAAM,EAAE,EAAE;IAYrB,2DAA2D;IAC3D,GAAG,CAAC,CAAC,EAAE,MAAM,GAAG,UAAU;CAG3B;AAED,uDAAuD;AACvD,MAAM,MAAM,WAAW,GAAG,MAAM,GAAG,SAAS,CAAC,MAAM,CAAC,EAAE,CAAC;AAEvD;;;;GAIG;AACH,wBAAgB,QAAQ,CAAC,CAAC,EAAE,WAAW,GAAG,MAAM,CAW/C"}

package/dist/matrix.js

@@ -0,0 +1,99 @@
+import { constructorFor, dtypeOf } from './types.js';
+/**
+ * A dense row-major matrix backed by a Float64Array or Float32Array.
+ *
+ * This is the container returned by all pca-web methods. It is a thin,
+ * zero-copy wrapper: `data[i * cols + j]` is the element at row `i`,
+ * column `j`.
+ */
+export class Matrix {
+    data;
+    rows;
+    cols;
+    /**
+     * Wraps an existing typed array (no copy). `data.length` must equal
+     * `rows * cols`.
+     */
+    constructor(data, rows, cols) {
+        if (data.length !== rows * cols) {
+            throw new Error(`Matrix: data length ${data.length} does not match shape ${rows}x${cols} (= ${rows * cols})`);
+        }
+        this.data = data;
+        this.rows = rows;
+        this.cols = cols;
+    }
+    /** Storage dtype: 'float64' or 'float32'. */
+    get dtype() {
+        return dtypeOf(this.data);
+    }
+    /** Allocates a zero-initialized matrix. */
+    static zeros(rows, cols, dtype = 'float64') {
+        const Ctor = constructorFor(dtype);
+        return new Matrix(new Ctor(rows * cols), rows, cols);
+    }
+    /**
+     * Builds a matrix from an array of rows, e.g. `Matrix.from2d([[1, 2], [3, 4]])`.
+     * All rows must have equal length.
+     */
+    static from2d(rows2d, dtype = 'float64') {
+        const rows = rows2d.length;
+        if (rows === 0) {
+            throw new Error('Matrix.from2d: at least one row is required');
+        }
+        const cols = rows2d[0].length;
+        const Ctor = constructorFor(dtype);
+        const data = new Ctor(rows * cols);
+        for (let i = 0; i < rows; i++) {
+            const r = rows2d[i];
+            if (r.length !== cols) {
+                throw new Error(`Matrix.from2d: row ${i} has length ${r.length}, expected ${cols} (rows must be equal length)`);
+            }
+            for (let j = 0; j < cols; j++) {
+                data[i * cols + j] = r[j];
+            }
+        }
+        return new Matrix(data, rows, cols);
+    }
+    get(i, j) {
+        return this.data[i * this.cols + j];
+    }
+    set(i, j, v) {
+        this.data[i * this.cols + j] = v;
+    }
+    /** Returns a deep copy (same dtype). */
+    copy() {
+        return new Matrix(this.data.slice(), this.rows, this.cols);
+    }
+    /** Converts to a plain nested array of numbers. */
+    toArray() {
+        const out = new Array(this.rows);
+        for (let i = 0; i < this.rows; i++) {
+            const row = new Array(this.cols);
+            for (let j = 0; j < this.cols; j++) {
+                row[j] = this.data[i * this.cols + j];
+            }
+            out[i] = row;
+        }
+        return out;
+    }
+    /** Copy of row `i` as a plain typed array (same dtype). */
+    row(i) {
+        return this.data.slice(i * this.cols, (i + 1) * this.cols);
+    }
+}
+/**
+ * Coerces user input to a Matrix without copying when it is already one.
+ * Nested arrays become float64 (mirrors sklearn coercing non-float dtypes
+ * to float64).
+ */
+export function asMatrix(X) {
+    if (X instanceof Matrix) {
+        return X;
+    }
+    if (Array.isArray(X)) {
+        return Matrix.from2d(X);
+    }
+    throw new TypeError('Expected a Matrix or an array of rows (number[][]); ' +
+        'wrap typed-array data with new Matrix(data, rows, cols)');
+}
+//# sourceMappingURL=matrix.js.map

package/dist/matrix.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"matrix.js","sourceRoot":"","sources":["../src/matrix.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,cAAc,EAAc,OAAO,EAAmB,MAAM,YAAY,CAAC;AAElF;;;;;;GAMG;AACH,MAAM,OAAO,MAAM;IACR,IAAI,CAAa;IACjB,IAAI,CAAS;IACb,IAAI,CAAS;IAEtB;;;OAGG;IACH,YAAY,IAAgB,EAAE,IAAY,EAAE,IAAY;QACtD,IAAI,IAAI,CAAC,MAAM,KAAK,IAAI,GAAG,IAAI,EAAE,CAAC;YAChC,MAAM,IAAI,KAAK,CACb,uBAAuB,IAAI,CAAC,MAAM,yBAAyB,IAAI,IAAI,IAAI,OAAO,IAAI,GAAG,IAAI,GAAG,CAC7F,CAAC;QACJ,CAAC;QACD,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC;QACjB,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC;QACjB,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC;IACnB,CAAC;IAED,6CAA6C;IAC7C,IAAI,KAAK;QACP,OAAO,OAAO,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IAC5B,CAAC;IAED,2CAA2C;IAC3C,MAAM,CAAC,KAAK,CAAC,IAAY,EAAE,IAAY,EAAE,QAAe,SAAS;QAC/D,MAAM,IAAI,GAAG,cAAc,CAAC,KAAK,CAAC,CAAC;QACnC,OAAO,IAAI,MAAM,CAAC,IAAI,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;IACvD,CAAC;IAED;;;OAGG;IACH,MAAM,CAAC,MAAM,CAAC,MAA2B,EAAE,QAAe,SAAS;QACjE,MAAM,IAAI,GAAG,MAAM,CAAC,MAAM,CAAC;QAC3B,IAAI,IAAI,KAAK,CAAC,EAAE,CAAC;YACf,MAAM,IAAI,KAAK,CAAC,6CAA6C,CAAC,CAAC;QACjE,CAAC;QACD,MAAM,IAAI,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;QAC9B,MAAM,IAAI,GAAG,cAAc,CAAC,KAAK,CAAC,CAAC;QACnC,MAAM,IAAI,GAAG,IAAI,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC;QACnC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;YAC9B,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;YACpB,IAAI,CAAC,CAAC,MAAM,KAAK,IAAI,EAAE,CAAC;gBACtB,MAAM,IAAI,KAAK,CACb,sBAAsB,CAAC,eAAe,CAAC,CAAC,MAAM,cAAc,IAAI,8BAA8B,CAC/F,CAAC;YACJ,CAAC;YACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC9B,IAAI,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YAC5B,CAAC;QACH,CAAC;QACD,OAAO,IAAI,MAAM,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;IACtC,CAAC;IAED,GAAG,CAAC,CAAS,EAAE,CAAS;QACtB,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;IACtC,CAAC;IAED,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS;QACjC,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;IACnC,CAAC;IAED,wCAAwC;IACxC,IAAI;QACF,OAAO,IAAI,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,IAAI,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC;IAC7D,CAAC;IAED,mDAAmD;IACnD,OAAO;QACL,MAAM,GAAG,GAAe,IAAI,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QAC7C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;YACnC,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,IAAI,CAAC,IAAI,CAAC,CAAC;YACzC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;gBACnC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;YACxC,CAAC;YACD,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QACf,CAAC;QACD,OAAO,GAAG,CAAC;IACb,CAAC;IAED,2DAA2D;IAC3D,GAAG,CAAC,CAAS;QACX,OAAO,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC;IAC7D,CAAC;CACF;AAKD;;;;GAIG;AACH,MAAM,UAAU,QAAQ,CAAC,CAAc;IACrC,IAAI,CAAC,YAAY,MAAM,EAAE,CAAC;QACxB,OAAO,CAAC,CAAC;IACX,CAAC;IACD,IAAI,KAAK,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC;QACrB,OAAO,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;IAC1B,CAAC;IACD,MAAM,IAAI,SAAS,CACjB,sDAAsD;QACpD,yDAAyD,CAC5D,CAAC;AACJ,CAAC"}