@danielsimonjr/mathts-functions 0.70.1 → 0.72.0

package/dist/signal/waveforms.d.ts

@@ -0,0 +1,52 @@
+/**
+ * Test waveform generators (scipy.signal.waveforms parity).
+ *
+ * Phase offsets `phi` are in radians (scipy takes degrees for chirp/sweep_poly).
+ *
+ * @packageDocumentation
+ */
+type Vec = readonly number[] | Float64Array;
+/** Frequency sweep method for {@link chirp}. */
+export type ChirpMethod = 'linear' | 'quadratic' | 'logarithmic' | 'hyperbolic';
+/** Options for {@link chirp}. */
+export interface ChirpOptions {
+    /** Sweep law (default `'linear'`). */
+    method?: ChirpMethod;
+    /** Phase offset in radians (default 0). */
+    phi?: number;
+}
+/**
+ * Cosine chirp sweeping from f0 at t = 0 to f1 at t = t1 (frequencies in Hz).
+ *
+ * @example
+ * chirp([0], 1, 1, 10); // [1]
+ */
+export declare function chirp(t: Vec, f0: number, t1: number, f1: number, opts?: ChirpOptions): number[];
+/**
+ * Square wave: +1 while frac(t / 2 pi) < duty, otherwise -1.
+ */
+export declare function squareWave(t: Vec, duty?: number): number[];
+/**
+ * Sawtooth wave: -1 at phase 0 rising to +1 at phase 2 pi * width, then falling
+ * back. `width = 1` is a rising ramp; `width = 0.5` is a triangle.
+ */
+export declare function sawtooth(t: Vec, width?: number): number[];
+/**
+ * Gaussian-modulated cosine pulse.
+ *
+ * @param fc - Centre frequency (Hz)
+ * @param bw - Fractional bandwidth (default 0.5)
+ * @param bwr - Reference level in dB for the bandwidth (default -6)
+ */
+export declare function gaussPulse(t: Vec, fc: number, bw?: number, bwr?: number): number[];
+/**
+ * Discrete unit impulse of length n. `idx` defaults to 0, `'mid'` is floor(n / 2).
+ */
+export declare function unitImpulse(n: number, idx?: number | 'mid'): number[];
+/**
+ * Cosine with a polynomial instantaneous frequency f(t) = sum_k c_k t^(p-k) Hz.
+ * `polyCoeffs` are highest-order first (numpy poly1d order).
+ */
+export declare function sweepPoly(t: Vec, polyCoeffs: readonly number[], phi?: number): number[];
+export {};
+//# sourceMappingURL=waveforms.d.ts.map

package/dist/signal/waveforms.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"waveforms.d.ts","sourceRoot":"","sources":["../../src/signal/waveforms.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AAEH,KAAK,GAAG,GAAG,SAAS,MAAM,EAAE,GAAG,YAAY,CAAC;AAI5C,gDAAgD;AAChD,MAAM,MAAM,WAAW,GAAG,QAAQ,GAAG,WAAW,GAAG,aAAa,GAAG,YAAY,CAAC;AAEhF,iCAAiC;AACjC,MAAM,WAAW,YAAY;IAC3B,sCAAsC;IACtC,MAAM,CAAC,EAAE,WAAW,CAAC;IACrB,2CAA2C;IAC3C,GAAG,CAAC,EAAE,MAAM,CAAC;CACd;AAED;;;;;GAKG;AACH,wBAAgB,KAAK,CAAC,CAAC,EAAE,GAAG,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,IAAI,CAAC,EAAE,YAAY,GAAG,MAAM,EAAE,CAM/F;AAuCD;;GAEG;AACH,wBAAgB,UAAU,CAAC,CAAC,EAAE,GAAG,EAAE,IAAI,GAAE,MAAY,GAAG,MAAM,EAAE,CAM/D;AAED;;;GAGG;AACH,wBAAgB,QAAQ,CAAC,CAAC,EAAE,GAAG,EAAE,KAAK,GAAE,MAAU,GAAG,MAAM,EAAE,CAO5D;AAED;;;;;;GAMG;AACH,wBAAgB,UAAU,CAAC,CAAC,EAAE,GAAG,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,GAAE,MAAY,EAAE,GAAG,GAAE,MAAW,GAAG,MAAM,EAAE,CAI3F;AAED;;GAEG;AACH,wBAAgB,WAAW,CAAC,CAAC,EAAE,MAAM,EAAE,GAAG,GAAE,MAAM,GAAG,KAAS,GAAG,MAAM,EAAE,CASxE;AAED;;;GAGG;AACH,wBAAgB,SAAS,CAAC,CAAC,EAAE,GAAG,EAAE,UAAU,EAAE,SAAS,MAAM,EAAE,EAAE,GAAG,GAAE,MAAU,GAAG,MAAM,EAAE,CAW1F"}

package/dist/signal/windows-extra.d.ts

@@ -0,0 +1,52 @@
+/**
+ * Extra window functions (scipy.signal.windows parity).
+ *
+ * Every window takes `opts.sym` (default `true`, symmetric). `sym: false` builds
+ * the periodic variant (length n+1, last sample dropped), as scipy does.
+ *
+ * @packageDocumentation
+ */
+/** Window names accepted by {@link getWindow}. */
+export type ExtraWindowName = 'kaiser' | 'gaussian' | 'tukey' | 'blackmanharris' | 'nuttall' | 'flattop' | 'bohman' | 'parzen' | 'cosine' | 'exponential';
+/** Options shared by the window constructors. */
+export interface WindowOptions {
+    /** Symmetric window (default `true`). `false` gives the periodic variant. */
+    sym?: boolean;
+}
+/**
+ * Kaiser window with shape parameter beta (beta = 0 gives a rectangular window).
+ *
+ * @example
+ * kaiserWindow(5, 0); // [1, 1, 1, 1, 1]
+ */
+export declare function kaiserWindow(n: number, beta: number, opts?: WindowOptions): number[];
+/** Gaussian window with standard deviation `std` (in samples). */
+export declare function gaussianWindow(n: number, std: number, opts?: WindowOptions): number[];
+/**
+ * Tukey (tapered cosine) window. alpha = 0 gives rectangular, alpha = 1 gives Hann.
+ */
+export declare function tukeyWindow(n: number, alpha?: number, opts?: WindowOptions): number[];
+/** Blackman-Harris 4-term window (scipy coefficients). */
+export declare function blackmanHarrisWindow(n: number, opts?: WindowOptions): number[];
+/** Nuttall 4-term window. */
+export declare function nuttallWindow(n: number, opts?: WindowOptions): number[];
+/** Flat-top 5-term window (scipy coefficients). */
+export declare function flattopWindow(n: number, opts?: WindowOptions): number[];
+/** Bohman window (zero at both ends). */
+export declare function bohmanWindow(n: number, opts?: WindowOptions): number[];
+/** Parzen window (piecewise cubic). */
+export declare function parzenWindow(n: number, opts?: WindowOptions): number[];
+/** Cosine (half-sine) window: sin(pi (i + 0.5) / M). */
+export declare function cosineWindow(n: number, opts?: WindowOptions): number[];
+/**
+ * Exponential (two-sided) window exp(-|i - center| / tau), centred on (M-1)/2.
+ */
+export declare function exponentialWindow(n: number, tau: number, opts?: WindowOptions): number[];
+/**
+ * Build a window by name. `param` is beta (kaiser, default 0), std (gaussian,
+ * default 1), alpha (tukey, default 0.5) or tau (exponential, default 1). The
+ * classic names (hamming, hann, blackman, bartlett, rectangular) delegate to
+ * {@link windowFunction} and are always symmetric.
+ */
+export declare function getWindow(name: ExtraWindowName | 'hamming' | 'hann' | 'blackman' | 'bartlett' | 'rectangular', n: number, param?: number): number[];
+//# sourceMappingURL=windows-extra.d.ts.map

package/dist/signal/windows-extra.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"windows-extra.d.ts","sourceRoot":"","sources":["../../src/signal/windows-extra.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAKH,kDAAkD;AAClD,MAAM,MAAM,eAAe,GACvB,QAAQ,GACR,UAAU,GACV,OAAO,GACP,gBAAgB,GAChB,SAAS,GACT,SAAS,GACT,QAAQ,GACR,QAAQ,GACR,QAAQ,GACR,aAAa,CAAC;AAElB,iDAAiD;AACjD,MAAM,WAAW,aAAa;IAC5B,6EAA6E;IAC7E,GAAG,CAAC,EAAE,OAAO,CAAC;CACf;AAmCD;;;;;GAKG;AACH,wBAAgB,YAAY,CAAC,CAAC,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,IAAI,CAAC,EAAE,aAAa,GAAG,MAAM,EAAE,CAWpF;AAED,kEAAkE;AAClE,wBAAgB,cAAc,CAAC,CAAC,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,IAAI,CAAC,EAAE,aAAa,GAAG,MAAM,EAAE,CAWrF;AAED;;GAEG;AACH,wBAAgB,WAAW,CAAC,CAAC,EAAE,MAAM,EAAE,KAAK,GAAE,MAAY,EAAE,IAAI,CAAC,EAAE,aAAa,GAAG,MAAM,EAAE,CAoB1F;AAED,0DAA0D;AAC1D,wBAAgB,oBAAoB,CAAC,CAAC,EAAE,MAAM,EAAE,IAAI,CAAC,EAAE,aAAa,GAAG,MAAM,EAAE,CAE9E;AAED,6BAA6B;AAC7B,wBAAgB,aAAa,CAAC,CAAC,EAAE,MAAM,EAAE,IAAI,CAAC,EAAE,aAAa,GAAG,MAAM,EAAE,CAEvE;AAED,mDAAmD;AACnD,wBAAgB,aAAa,CAAC,CAAC,EAAE,MAAM,EAAE,IAAI,CAAC,EAAE,aAAa,GAAG,MAAM,EAAE,CAIvE;AAED,yCAAyC;AACzC,wBAAgB,YAAY,CAAC,CAAC,EAAE,MAAM,EAAE,IAAI,CAAC,EAAE,aAAa,GAAG,MAAM,EAAE,CAStE;AAED,uCAAuC;AACvC,wBAAgB,YAAY,CAAC,CAAC,EAAE,MAAM,EAAE,IAAI,CAAC,EAAE,aAAa,GAAG,MAAM,EAAE,CAYtE;AAED,wDAAwD;AACxD,wBAAgB,YAAY,CAAC,CAAC,EAAE,MAAM,EAAE,IAAI,CAAC,EAAE,aAAa,GAAG,MAAM,EAAE,CAItE;AAED;;GAEG;AACH,wBAAgB,iBAAiB,CAAC,CAAC,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,IAAI,CAAC,EAAE,aAAa,GAAG,MAAM,EAAE,CAMxF;AAED;;;;;GAKG;AACH,wBAAgB,SAAS,CACvB,IAAI,EAAE,eAAe,GAAG,SAAS,GAAG,MAAM,GAAG,UAAU,GAAG,UAAU,GAAG,aAAa,EACpF,CAAC,EAAE,MAAM,EACT,KAAK,CAAC,EAAE,MAAM,GACb,MAAM,EAAE,CA+BV"}

package/dist/signal/zpk2tf.d.ts

@@ -1,3 +1,23 @@
-import type { TypedFunction } from '../types.js';
-export declare const createZpk2tf: import("../utils/factory.js").FactoryFunction<any, TypedFunction>;
+import type { Matrix, Complex, TypedFunction } from '../types.js';
+/**
+ * Transfer function representation [numerator, denominator]
+ */
+export type TransferFunction = [Complex[], Complex[]];
+/**
+ * Zero, pole, or gain value (can be number, Complex, or BigNumber)
+ */
+export type ZPKValue = number | Complex | {
+    type: string;
+    re?: number;
+    im?: number;
+};
+/** Call signatures of `zpk2tf` (zeros, poles, optional gain defaulting to 1). */
+interface Zpk2tfOverloads {
+    (z: ZPKValue[], p: ZPKValue[], k?: number): TransferFunction;
+    (z: Matrix, p: Matrix, k?: number): TransferFunction;
+}
+/** `zpk2tf` with its precise call signatures; still a `TypedFunction`. */
+export type Zpk2tfFunction = Zpk2tfOverloads & TypedFunction;
+export declare const createZpk2tf: import("../utils/factory.js").FactoryFunction<any, Zpk2tfFunction>;
+export {};
 //# sourceMappingURL=zpk2tf.d.ts.map

package/dist/signal/zpk2tf.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"zpk2tf.d.ts","sourceRoot":"","sources":["../../src/signal/zpk2tf.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EAAmB,aAAa,EAAE,MAAM,aAAa,CAAC;AAgBlE,eAAO,MAAM,YAAY,mEAuGxB,CAAC"}
+{"version":3,"file":"zpk2tf.d.ts","sourceRoot":"","sources":["../../src/signal/zpk2tf.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EAAE,MAAM,EAAE,OAAO,EAAE,aAAa,EAAE,MAAM,aAAa,CAAC;AAMlE;;GAEG;AACH,MAAM,MAAM,gBAAgB,GAAG,CAAC,OAAO,EAAE,EAAE,OAAO,EAAE,CAAC,CAAC;AAEtD;;GAEG;AACH,MAAM,MAAM,QAAQ,GAAG,MAAM,GAAG,OAAO,GAAG;IAAE,IAAI,EAAE,MAAM,CAAC;IAAC,EAAE,CAAC,EAAE,MAAM,CAAC;IAAC,EAAE,CAAC,EAAE,MAAM,CAAA;CAAE,CAAC;AAErF,iFAAiF;AACjF,UAAU,eAAe;IACvB,CAAC,CAAC,EAAE,QAAQ,EAAE,EAAE,CAAC,EAAE,QAAQ,EAAE,EAAE,CAAC,CAAC,EAAE,MAAM,GAAG,gBAAgB,CAAC;IAC7D,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,CAAC,EAAE,MAAM,GAAG,gBAAgB,CAAC;CACtD;AAED,0EAA0E;AAC1E,MAAM,MAAM,cAAc,GAAG,eAAe,GAAG,aAAa,CAAC;AAE7D,eAAO,MAAM,YAAY,oEAiHxB,CAAC"}

package/dist/special/angular-momentum.d.ts

@@ -0,0 +1,38 @@
+/**
+ * Angular-momentum coupling (Wigner 3j, Clebsch–Gordan, Wigner 6j) via the Racah
+ * formulas, and spherical harmonics with the Condon–Shortley phase.
+ *
+ * All angular momenta may be integers or half-integers. Inputs that are not
+ * multiples of ½ raise `RangeError`; selection-rule violations return 0.
+ */
+import { Complex } from '@danielsimonjr/mathts-core';
+/**
+ * Wigner 3j symbol `(j1 j2 j3; m1 m2 m3)` by the Racah formula.
+ *
+ * @example
+ * wigner3j(1, 1, 0, 0, 0, 0) // -0.5773502691896258
+ */
+export declare function wigner3j(j1: number, j2: number, j3: number, m1: number, m2: number, m3: number): number;
+/**
+ * Clebsch–Gordan coefficient `⟨j1 m1; j2 m2 | J M⟩`.
+ *
+ * @example
+ * clebschGordan(1, 0, 1, 0, 2, 0) // 0.816496580927726
+ */
+export declare function clebschGordan(j1: number, m1: number, j2: number, m2: number, J: number, M: number): number;
+/**
+ * Wigner 6j symbol `{j1 j2 j3; j4 j5 j6}` by the Racah formula.
+ *
+ * @example
+ * wigner6j(1, 1, 1, 1, 1, 1) // 1/6
+ */
+export declare function wigner6j(j1: number, j2: number, j3: number, j4: number, j5: number, j6: number): number;
+/**
+ * Spherical harmonic `Y_l^m(θ, φ)` with the Condon–Shortley phase.
+ * `θ` is the polar angle and `φ` the azimuth. Returns zero when `|m| > l`.
+ *
+ * @example
+ * sphericalHarmonic(0, 0, 0.3, 1.2) // re 0.28209479177387814, im 0
+ */
+export declare function sphericalHarmonic(l: number, m: number, theta: number, phi: number): Complex;
+//# sourceMappingURL=angular-momentum.d.ts.map

package/dist/special/angular-momentum.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"angular-momentum.d.ts","sourceRoot":"","sources":["../../src/special/angular-momentum.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AAEH,OAAO,EAAE,OAAO,EAAE,MAAM,4BAA4B,CAAC;AAiCrD;;;;;GAKG;AACH,wBAAgB,QAAQ,CACtB,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,GACT,MAAM,CAkCR;AAED;;;;;GAKG;AACH,wBAAgB,aAAa,CAC3B,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,GACR,MAAM,CAER;AAED;;;;;GAKG;AACH,wBAAgB,QAAQ,CACtB,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,GACT,MAAM,CAgDR;AAED;;;;;;GAMG;AACH,wBAAgB,iBAAiB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,OAAO,CAe3F"}

package/dist/special/erf-extra.d.ts

@@ -0,0 +1,47 @@
+/**
+ * Inverse and scaled error functions, and Dawson's integral.
+ *
+ * - `erfinv` / `erfcinv`: Newton iteration seeded from `normalQuantile`, with
+ *   the complementary branch used for `|y| > 1/2` to avoid cancellation.
+ * - `erfcx`: scaled complementary error function `e^{x²} erfc(x)`.
+ * - `dawson`: Dawson's integral `F(x) = e^{-x²} ∫₀ˣ e^{t²} dt`.
+ *
+ * Scalar kernels only (number in, number out).
+ */
+/**
+ * Scaled complementary error function `erfcx(x) = e^{x²} erfc(x)`.
+ *
+ * @example
+ * erfcx(0) // 1
+ * erfcx(1) // 0.42758357615580705
+ */
+export declare function erfcx(x: number): number;
+/**
+ * Inverse error function: `erf(erfinv(y)) = y` for `y ∈ [−1, 1]`.
+ *
+ * @example
+ * erfinv(0)          // 0
+ * erfinv(0.5)        // 0.4769362762044699
+ * erfinv(1)          // Infinity
+ */
+export declare function erfinv(y: number): number;
+/**
+ * Inverse complementary error function: `erfc(erfcinv(y)) = y` for `y ∈ [0, 2]`.
+ *
+ * @example
+ * erfcinv(1)  // 0
+ * erfcinv(0)  // Infinity
+ */
+export declare function erfcinv(y: number): number;
+/**
+ * Dawson's integral `F(x) = e^{-x²} ∫₀ˣ e^{t²} dt`.
+ *
+ * Uses the positive series `F(x) = e^{-x²} Σ x^{2n+1}/(n!(2n+1))` for |x| ≤ 25
+ * (no cancellation) and a four-term asymptotic expansion beyond.
+ *
+ * @example
+ * dawson(0)  // 0
+ * dawson(1)  // 0.5380795069127684
+ */
+export declare function dawson(x: number): number;
+//# sourceMappingURL=erf-extra.d.ts.map

package/dist/special/erf-extra.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"erf-extra.d.ts","sourceRoot":"","sources":["../../src/special/erf-extra.ts"],"names":[],"mappings":"AAAA;;;;;;;;;GASG;AAmCH;;;;;;GAMG;AACH,wBAAgB,KAAK,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAKvC;AAiBD;;;;;;;GAOG;AACH,wBAAgB,MAAM,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAmBxC;AAED;;;;;;GAMG;AACH,wBAAgB,OAAO,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAOzC;AAED;;;;;;;;;GASG;AACH,wBAAgB,MAAM,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAoBxC"}

package/dist/special/exp-integrals.d.ts

@@ -0,0 +1,38 @@
+/**
+ * Exponential-type integrals: E1 / E_n, Shi (sinh integral) and Chi (cosh integral).
+ */
+/**
+ * Generalized exponential integral `E_n(x) = ∫₁^∞ e^{-xt} t^{-n} dt`.
+ *
+ * - `n = 0`: `e^{-x}/x`.
+ * - `x ≤ 1`: E1 by its power series, then the (stable for x < 1) upward recurrence.
+ * - `x > 1`: continued fraction.
+ *
+ * @example
+ * expn(2, 1)  // 0.14849550677592205
+ * expn(3, 0)  // 0.5
+ */
+export declare function expn(n: number, x: number): number;
+/**
+ * Exponential integral E1(x) for real x > 0 (`−Ei(−x)`). Returns Infinity at 0, NaN for x < 0.
+ *
+ * @example
+ * expIntegralE1(1) // 0.21938393439552029
+ */
+export declare function expIntegralE1(x: number): number;
+/**
+ * Sinh integral `Shi(x) = ∫₀ˣ sinh(t)/t dt` (odd function).
+ *
+ * @example
+ * sinhIntegral(1) // 1.0572508753757285
+ */
+export declare function sinhIntegral(x: number): number;
+/**
+ * Cosh integral `Chi(x) = γ + ln x + ∫₀ˣ (cosh t − 1)/t dt` for x > 0.
+ * Returns −Infinity at 0 and NaN for x < 0.
+ *
+ * @example
+ * coshIntegral(1) // 0.8378669409802082
+ */
+export declare function coshIntegral(x: number): number;
+//# sourceMappingURL=exp-integrals.d.ts.map

package/dist/special/exp-integrals.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"exp-integrals.d.ts","sourceRoot":"","sources":["../../src/special/exp-integrals.ts"],"names":[],"mappings":"AAAA;;GAEG;AAyBH;;;;;;;;;;GAUG;AACH,wBAAgB,IAAI,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAoBjD;AAED;;;;;GAKG;AACH,wBAAgB,aAAa,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAI/C;AAED;;;;;GAKG;AACH,wBAAgB,YAAY,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAe9C;AAED;;;;;;GAMG;AACH,wBAAgB,YAAY,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAc9C"}

package/dist/special/spherical-bessel.d.ts

@@ -0,0 +1,28 @@
+/**
+ * Spherical Bessel functions j_n, y_n, i_n, k_n, Hankel functions of real order,
+ * and the first positive zeros of J_ν.
+ *
+ * Conventions: `j_n(x) = √(π/2x) J_{n+1/2}(x)`, `y_n` likewise from Y, `i_n = √(π/2x) I_{n+1/2}(x)`,
+ * `k_n = √(π/2x) K_{n+1/2}(x)` (so `k_0(x) = (π/2) e^{-x}/x`).
+ */
+import { Complex } from '@danielsimonjr/mathts-core';
+/** Spherical Bessel function of the first kind `j_n(x)`. */
+export declare function sphericalBesselJ(n: number, x: number): number;
+/** Spherical Bessel function of the second kind `y_n(x)`. */
+export declare function sphericalBesselY(n: number, x: number): number;
+/** Modified spherical Bessel function of the first kind `i_n(x)`. */
+export declare function sphericalBesselI(n: number, x: number): number;
+/** Modified spherical Bessel function of the second kind `k_n(x)`. */
+export declare function sphericalBesselK(n: number, x: number): number;
+/** Hankel function of the first kind `H¹_ν(x) = J_ν(x) + i Y_ν(x)`. */
+export declare function hankelH1(nu: number, x: number): Complex;
+/** Hankel function of the second kind `H²_ν(x) = J_ν(x) − i Y_ν(x)`. */
+export declare function hankelH2(nu: number, x: number): Complex;
+/**
+ * First `k` positive zeros of `J_ν(x)`, via McMahon's asymptotic guess refined by Newton.
+ *
+ * @example
+ * besselJZeros(0, 2) // [2.404825557695773, 5.520078110286311]
+ */
+export declare function besselJZeros(nu: number, k: number): number[];
+//# sourceMappingURL=spherical-bessel.d.ts.map

package/dist/special/spherical-bessel.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"spherical-bessel.d.ts","sourceRoot":"","sources":["../../src/special/spherical-bessel.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AAEH,OAAO,EAAE,OAAO,EAAE,MAAM,4BAA4B,CAAC;AA4FrD,4DAA4D;AAC5D,wBAAgB,gBAAgB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAG7D;AAED,6DAA6D;AAC7D,wBAAgB,gBAAgB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAG7D;AAED,qEAAqE;AACrE,wBAAgB,gBAAgB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAG7D;AAED,sEAAsE;AACtE,wBAAgB,gBAAgB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAG7D;AAED,uEAAuE;AACvE,wBAAgB,QAAQ,CAAC,EAAE,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,OAAO,CAIvD;AAED,wEAAwE;AACxE,wBAAgB,QAAQ,CAAC,EAAE,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,OAAO,CAIvD;AAYD;;;;;GAKG;AACH,wBAAgB,YAAY,CAAC,EAAE,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,EAAE,CAkB5D"}

package/dist/special/zeta-extra.d.ts

@@ -0,0 +1,36 @@
+/**
+ * Hurwitz zeta, Dirichlet eta, and Bernoulli / Euler polynomials.
+ */
+/**
+ * Hurwitz zeta `ζ(s, a) = Σ_{k≥0} (k + a)^{−s}` for `s > 1`, `a > 0`, by Euler–Maclaurin
+ * with `N = 10` leading terms and eight Bernoulli correction terms. NaN outside the domain.
+ *
+ * @example
+ * hurwitzZeta(2, 1) // 1.6449340668482264
+ */
+export declare function hurwitzZeta(s: number, a: number): number;
+/**
+ * Dirichlet eta function `η(s) = (1 − 2^{1−s}) ζ(s)`, all real s, by Borwein's
+ * alternating-series algorithm (n = 40, error below 1e-20 relative).
+ *
+ * @example
+ * dirichletEta(1) // 0.6931471805599453
+ * dirichletEta(2) // 0.8224670334241132
+ */
+export declare function dirichletEta(s: number): number;
+/**
+ * Bernoulli polynomial `B_n(x) = Σ_k C(n,k) B_k x^{n−k}` (with `B_1 = −1/2`).
+ *
+ * @example
+ * bernoulliPolynomial(2, 0.5) // -0.08333333333333333
+ */
+export declare function bernoulliPolynomial(n: number, x: number): number;
+/**
+ * Euler polynomial `E_n(x) = 2/(n+1) · (B_{n+1}(x) − 2^{n+1} B_{n+1}(x/2))`.
+ *
+ * @example
+ * eulerPolynomial(2, 0.5) // -0.25
+ * eulerPolynomial(1, 0)   // -0.5
+ */
+export declare function eulerPolynomial(n: number, x: number): number;
+//# sourceMappingURL=zeta-extra.d.ts.map

package/dist/special/zeta-extra.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"zeta-extra.d.ts","sourceRoot":"","sources":["../../src/special/zeta-extra.ts"],"names":[],"mappings":"AAAA;;GAEG;AAcH;;;;;;GAMG;AACH,wBAAgB,WAAW,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAoBxD;AAED;;;;;;;GAOG;AACH,wBAAgB,YAAY,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAuB9C;AAED;;;;;GAKG;AACH,wBAAgB,mBAAmB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAWhE;AAED;;;;;;GAMG;AACH,wBAAgB,eAAe,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAM5D"}

package/dist/stats/distances.d.ts

@@ -0,0 +1,28 @@
+/**
+ * Distances between samples and discrete information measures (gap-analysis
+ * 2026-11-04, Batch A).
+ */
+type Vec = readonly number[] | Float64Array;
+/**
+ * 1-D Wasserstein (earth mover's) distance: ∫|F_u(t) − F_v(t)| dt, computed
+ * exactly over the merged support. Default weights are uniform.
+ */
+export declare function wassersteinDistance(u: Vec, v: Vec, opts?: {
+    uWeights?: Vec;
+    vWeights?: Vec;
+}): number;
+/**
+ * Energy distance √(2E|U−V| − E|U−U′| − E|V−V′|) using V-statistic means
+ * (diagonal pairs included).
+ */
+export declare function energyDistance(u: Vec, v: Vec): number;
+/** Mutual information (nats) of two discrete label sequences. */
+export declare function mutualInformation(x: readonly number[], y: readonly number[]): number;
+/** Empirical CDF: step points and an `evaluate(t)` function (fraction of values ≤ t). */
+export declare function ecdf(x: Vec): {
+    x: number[];
+    y: number[];
+    evaluate: (t: number) => number;
+};
+export {};
+//# sourceMappingURL=distances.d.ts.map

package/dist/stats/distances.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"distances.d.ts","sourceRoot":"","sources":["../../src/stats/distances.ts"],"names":[],"mappings":"AAAA;;;GAGG;AACH,KAAK,GAAG,GAAG,SAAS,MAAM,EAAE,GAAG,YAAY,CAAC;AAI5C;;;GAGG;AACH,wBAAgB,mBAAmB,CACjC,CAAC,EAAE,GAAG,EACN,CAAC,EAAE,GAAG,EACN,IAAI,CAAC,EAAE;IAAE,QAAQ,CAAC,EAAE,GAAG,CAAC;IAAC,QAAQ,CAAC,EAAE,GAAG,CAAA;CAAE,GACxC,MAAM,CAwBR;AAgBD;;;GAGG;AACH,wBAAgB,cAAc,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,MAAM,CAWrD;AAED,iEAAiE;AACjE,wBAAgB,iBAAiB,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,EAAE,CAAC,EAAE,SAAS,MAAM,EAAE,GAAG,MAAM,CAoBpF;AAED,yFAAyF;AACzF,wBAAgB,IAAI,CAAC,CAAC,EAAE,GAAG,GAAG;IAAE,CAAC,EAAE,MAAM,EAAE,CAAC;IAAC,CAAC,EAAE,MAAM,EAAE,CAAC;IAAC,QAAQ,EAAE,CAAC,CAAC,EAAE,MAAM,KAAK,MAAM,CAAA;CAAE,CAyB1F"}

package/dist/stats/distributions-extra.d.ts

@@ -0,0 +1,36 @@
+/**
+ * Owen's T function T(h, a) = (1/2π)∫₀ᵃ exp(−h²(1+x²)/2)/(1+x²) dx.
+ */
+export declare function owensT(h: number, a: number): number;
+/** Skew-normal PDF: 2φ(z)Φ(αz)/scale with z = (x − loc)/scale. */
+export declare function skewNormalPDF(x: number, alpha: number, loc?: number, scale?: number): number;
+/** Skew-normal CDF: Φ(z) − 2T(z, α). */
+export declare function skewNormalCDF(x: number, alpha: number, loc?: number, scale?: number): number;
+/** Truncated normal PDF on [a, b]. */
+export declare function truncatedNormalPDF(x: number, a: number, b: number, mu?: number, sigma?: number): number;
+/** Truncated normal CDF on [a, b]. */
+export declare function truncatedNormalCDF(x: number, a: number, b: number, mu?: number, sigma?: number): number;
+/** Mean of the truncated normal on [a, b]: μ + σ(φ(α) − φ(β))/(Φ(β) − Φ(α)), tail-stable. */
+export declare function truncatedNormalMean(a: number, b: number, mu?: number, sigma?: number): number;
+/** Inverse gamma PDF: β^α x^(−α−1) e^(−β/x) / Γ(α). */
+export declare function inverseGammaPDF(x: number, alpha: number, beta: number): number;
+/** Inverse gamma CDF: 1 − P(α, β/x) (regularized upper incomplete gamma). */
+export declare function inverseGammaCDF(x: number, alpha: number, beta: number): number;
+/** GEV CDF exp(−t), t = (1 + ξz)^(−1/ξ) (ξ = 0: t = e^(−z)). Outside support: 0 or 1. */
+export declare function gevCDF(x: number, xi: number, mu?: number, sigma?: number): number;
+/** GEV PDF t^(ξ+1) e^(−t) / σ, with t as in {@link gevCDF}. Zero outside support. */
+export declare function gevPDF(x: number, xi: number, mu?: number, sigma?: number): number;
+/**
+ * Skellam PMF: e^(−(μ₁+μ₂)) (μ₁/μ₂)^(k/2) I_|k|(2√(μ₁μ₂)). Requires μ₁, μ₂ > 0.
+ */
+export declare function skellamPMF(k: number, mu1: number, mu2: number): number;
+/** Zipf PMF on {1..n}: k^(−s) / H_{n,s}. Zero for non-integer or out-of-range k. */
+export declare function zipfPMF(k: number, s: number, n: number): number;
+/** Nakagami PDF: 2 m^m x^(2m−1) e^(−m x²/ω) / (Γ(m) ω^m), x > 0. */
+export declare function nakagamiPDF(x: number, m: number, omega: number): number;
+/**
+ * Rice PDF: (x/σ²) e^(−(x²+ν²)/(2σ²)) I₀(xν/σ²), x ≥ 0. Evaluated with the scaled
+ * I₀ so that exp(−(x−ν)²/(2σ²)) · I₀e does not overflow for large xν/σ².
+ */
+export declare function ricePDF(x: number, nu: number, sigma: number): number;
+//# sourceMappingURL=distributions-extra.d.ts.map

package/dist/stats/distributions-extra.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"distributions-extra.d.ts","sourceRoot":"","sources":["../../src/stats/distributions-extra.ts"],"names":[],"mappings":"AA2BA;;GAEG;AACH,wBAAgB,MAAM,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAMnD;AAED,kEAAkE;AAClE,wBAAgB,aAAa,CAAC,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,GAAG,SAAI,EAAE,KAAK,SAAI,GAAG,MAAM,CAGlF;AAED,wCAAwC;AACxC,wBAAgB,aAAa,CAAC,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,GAAG,SAAI,EAAE,KAAK,SAAI,GAAG,MAAM,CAGlF;AAuCD,sCAAsC;AACtC,wBAAgB,kBAAkB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,SAAI,EAAE,KAAK,SAAI,GAAG,MAAM,CAO7F;AAED,sCAAsC;AACtC,wBAAgB,kBAAkB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,SAAI,EAAE,KAAK,SAAI,GAAG,MAAM,CAU7F;AAED,6FAA6F;AAC7F,wBAAgB,mBAAmB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,SAAI,EAAE,KAAK,SAAI,GAAG,MAAM,CASnF;AAED,uDAAuD;AACvD,wBAAgB,eAAe,CAAC,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,GAAG,MAAM,CAG9E;AAED,6EAA6E;AAC7E,wBAAgB,eAAe,CAAC,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,GAAG,MAAM,CAG9E;AAED,yFAAyF;AACzF,wBAAgB,MAAM,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,SAAI,EAAE,KAAK,SAAI,GAAG,MAAM,CAMvE;AAED,qFAAqF;AACrF,wBAAgB,MAAM,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,SAAI,EAAE,KAAK,SAAI,GAAG,MAAM,CAWvE;AAED;;GAEG;AACH,wBAAgB,UAAU,CAAC,CAAC,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,MAAM,CAQtE;AAED,oFAAoF;AACpF,wBAAgB,OAAO,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAK/D;AAED,oEAAoE;AACpE,wBAAgB,WAAW,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,GAAG,MAAM,CAUvE;AA6BD;;;GAGG;AACH,wBAAgB,OAAO,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,GAAG,MAAM,CAMpE"}

package/dist/stats/effect-sizes.d.ts

@@ -0,0 +1,21 @@
+/**
+ * Standardized effect sizes (gap-analysis 2026-11-04, Batch A).
+ */
+type Vec = readonly number[] | Float64Array;
+/** Cohen's d with the pooled (ddof 1) standard deviation. */
+export declare function cohensD(a: Vec, b: Vec): number;
+/** Hedges' g: Cohen's d with the small-sample correction 1 − 3/(4(n₁+n₂) − 9). */
+export declare function hedgesG(a: Vec, b: Vec): number;
+/** Glass's Δ: (mean a − mean b) / sd(b). */
+export declare function glassDelta(a: Vec, b: Vec): number;
+/** η² = SS_between / SS_total across groups. */
+export declare function etaSquared(groups: readonly Vec[]): number;
+/**
+ * Rank-biserial correlation r = 2U_a/(n₁n₂) − 1, where U_a counts pairs with
+ * aᵢ > bⱼ (ties count ½). r is positive when a tends to exceed b, so its sign
+ * is opposite to `cohensD(a, b)`'s convention only in magnitude form; it equals
+ * 1 − 2U_b/(n₁n₂) with U_b = n₁n₂ − U_a.
+ */
+export declare function rankBiserial(a: Vec, b: Vec): number;
+export {};
+//# sourceMappingURL=effect-sizes.d.ts.map

package/dist/stats/effect-sizes.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"effect-sizes.d.ts","sourceRoot":"","sources":["../../src/stats/effect-sizes.ts"],"names":[],"mappings":"AAAA;;GAEG;AACH,KAAK,GAAG,GAAG,SAAS,MAAM,EAAE,GAAG,YAAY,CAAC;AAS5C,6DAA6D;AAC7D,wBAAgB,OAAO,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,MAAM,CAQ9C;AAED,kFAAkF;AAClF,wBAAgB,OAAO,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,MAAM,CAG9C;AAED,4CAA4C;AAC5C,wBAAgB,UAAU,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,MAAM,CAKjD;AAED,gDAAgD;AAChD,wBAAgB,UAAU,CAAC,MAAM,EAAE,SAAS,GAAG,EAAE,GAAG,MAAM,CAQzD;AAED;;;;;GAKG;AACH,wBAAgB,YAAY,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,MAAM,CAYnD"}

package/dist/stats/nonparametric-extra.d.ts

@@ -0,0 +1,52 @@
+type Vec = readonly number[] | Float64Array;
+/** Sign test on the median. Ties with `median` are dropped; two-sided exact binomial p. */
+export declare function signTest(x: Vec, medianValue?: number): {
+    statistic: number;
+    pValue: number;
+    nPos: number;
+    nNeg: number;
+};
+/**
+ * Wald–Wolfowitz runs test. Booleans are used directly; numeric values are
+ * classified against `cutoff` (default: the median), with values equal to the
+ * cutoff dropped. Normal approximation, two-sided p. With n₁ = 0 or n₂ = 0,
+ * z is NaN and p is 1.
+ */
+export declare function runsTest(x: Vec | readonly boolean[], opts?: {
+    cutoff?: number;
+}): {
+    runs: number;
+    z: number;
+    pValue: number;
+};
+/** Mood's median test: chi-square test of the above/at-or-below-grand-median table. */
+export declare function moodMedianTest(groups: readonly Vec[]): {
+    statistic: number;
+    pValue: number;
+    grandMedian: number;
+};
+/**
+ * Cressie–Read power divergence statistic for observed counts against expected
+ * counts (default: uniform). λ = 1 gives Pearson χ², λ = 0 gives G, λ = −1 is
+ * rejected as singular. Zero observed counts contribute 0 when λ = 0.
+ */
+export declare function powerDivergence(observed: Vec, expected?: Vec, lambda?: number): {
+    statistic: number;
+    pValue: number;
+};
+/** Log-likelihood ratio G-test: `powerDivergence(observed, expected, 0)`. */
+export declare function gTest(observed: Vec, expected?: Vec): {
+    statistic: number;
+    pValue: number;
+};
+/**
+ * Delete-one jackknife. `estimate` is the bias-corrected estimate
+ * θ̂ − bias, `bias` = (n−1)(θ̄₋ᵢ − θ̂), and `se` = √((n−1)/n · Σ(θ₋ᵢ − θ̄₋ᵢ)²).
+ */
+export declare function jackknife(x: Vec, stat: (sample: number[]) => number): {
+    estimate: number;
+    bias: number;
+    se: number;
+};
+export {};
+//# sourceMappingURL=nonparametric-extra.d.ts.map

package/dist/stats/nonparametric-extra.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"nonparametric-extra.d.ts","sourceRoot":"","sources":["../../src/stats/nonparametric-extra.ts"],"names":[],"mappings":"AAOA,KAAK,GAAG,GAAG,SAAS,MAAM,EAAE,GAAG,YAAY,CAAC;AAc5C,2FAA2F;AAC3F,wBAAgB,QAAQ,CACtB,CAAC,EAAE,GAAG,EACN,WAAW,SAAI,GACd;IAAE,SAAS,EAAE,MAAM,CAAC;IAAC,MAAM,EAAE,MAAM,CAAC;IAAC,IAAI,EAAE,MAAM,CAAC;IAAC,IAAI,EAAE,MAAM,CAAA;CAAE,CAUnE;AAQD;;;;;GAKG;AACH,wBAAgB,QAAQ,CACtB,CAAC,EAAE,GAAG,GAAG,SAAS,OAAO,EAAE,EAC3B,IAAI,CAAC,EAAE;IAAE,MAAM,CAAC,EAAE,MAAM,CAAA;CAAE,GACzB;IAAE,IAAI,EAAE,MAAM,CAAC;IAAC,CAAC,EAAE,MAAM,CAAC;IAAC,MAAM,EAAE,MAAM,CAAA;CAAE,CAmB7C;AAED,uFAAuF;AACvF,wBAAgB,cAAc,CAAC,MAAM,EAAE,SAAS,GAAG,EAAE,GAAG;IACtD,SAAS,EAAE,MAAM,CAAC;IAClB,MAAM,EAAE,MAAM,CAAC;IACf,WAAW,EAAE,MAAM,CAAC;CACrB,CAkBA;AAED;;;;GAIG;AACH,wBAAgB,eAAe,CAC7B,QAAQ,EAAE,GAAG,EACb,QAAQ,CAAC,EAAE,GAAG,EACd,MAAM,SAAI,GACT;IAAE,SAAS,EAAE,MAAM,CAAC;IAAC,MAAM,EAAE,MAAM,CAAA;CAAE,CAiBvC;AAED,6EAA6E;AAC7E,wBAAgB,KAAK,CAAC,QAAQ,EAAE,GAAG,EAAE,QAAQ,CAAC,EAAE,GAAG,GAAG;IAAE,SAAS,EAAE,MAAM,CAAC;IAAC,MAAM,EAAE,MAAM,CAAA;CAAE,CAE1F;AAED;;;GAGG;AACH,wBAAgB,SAAS,CACvB,CAAC,EAAE,GAAG,EACN,IAAI,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,KAAK,MAAM,GACjC;IAAE,QAAQ,EAAE,MAAM,CAAC;IAAC,IAAI,EAAE,MAAM,CAAC;IAAC,EAAE,EAAE,MAAM,CAAA;CAAE,CAWhD"}

package/dist/stats/robust-stats.d.ts

@@ -0,0 +1,38 @@
+/**
+ * Generalized means, geometric standard deviation, robust regression slopes
+ * and partial correlation (gap-analysis 2026-11-04, Batch A).
+ */
+type Vec = readonly number[] | Float64Array;
+/**
+ * Power mean M_p = (Σxᵖ / n)^(1/p). p = 0 is the geometric mean, p = 1 the
+ * arithmetic mean and p = −1 the harmonic mean. Non-integer p requires x > 0.
+ */
+export declare function pmean(x: Vec, p: number): number;
+/** Geometric standard deviation: exp(std(ln x)). Default ddof 1 (scipy). */
+export declare function gstd(x: Vec, opts?: {
+    ddof?: number;
+}): number;
+/** Geometric z-score: (ln xᵢ − mean(ln x)) / std(ln x). */
+export declare function gzscore(x: Vec, opts?: {
+    ddof?: number;
+}): number[];
+/**
+ * Theil–Sen estimator: slope = median of all pairwise slopes, intercept =
+ * median(y) − slope·median(x).
+ */
+export declare function theilSen(x: Vec, y: Vec): {
+    slope: number;
+    intercept: number;
+};
+/**
+ * Siegel repeated-medians estimator: slope = median over i of the median of
+ * slopes from point i to every other point; intercept = median(y − slope·x).
+ */
+export declare function siegelSlopes(x: Vec, y: Vec): {
+    slope: number;
+    intercept: number;
+};
+/** Partial correlation of x and y controlling for z (first-order formula). */
+export declare function partialCorrelation(x: Vec, y: Vec, z: Vec): number;
+export {};
+//# sourceMappingURL=robust-stats.d.ts.map

package/dist/stats/robust-stats.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"robust-stats.d.ts","sourceRoot":"","sources":["../../src/stats/robust-stats.ts"],"names":[],"mappings":"AAAA;;;GAGG;AACH,KAAK,GAAG,GAAG,SAAS,MAAM,EAAE,GAAG,YAAY,CAAC;AAoB5C;;;GAGG;AACH,wBAAgB,KAAK,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAe/C;AAED,4EAA4E;AAC5E,wBAAgB,IAAI,CAAC,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,EAAE;IAAE,IAAI,CAAC,EAAE,MAAM,CAAA;CAAE,GAAG,MAAM,CAK7D;AAED,2DAA2D;AAC3D,wBAAgB,OAAO,CAAC,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,EAAE;IAAE,IAAI,CAAC,EAAE,MAAM,CAAA;CAAE,GAAG,MAAM,EAAE,CAQlE;AAqBD;;;GAGG;AACH,wBAAgB,QAAQ,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG;IAAE,KAAK,EAAE,MAAM,CAAC;IAAC,SAAS,EAAE,MAAM,CAAA;CAAE,CAM7E;AAED;;;GAGG;AACH,wBAAgB,YAAY,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG;IAAE,KAAK,EAAE,MAAM,CAAC;IAAC,SAAS,EAAE,MAAM,CAAA;CAAE,CAajF;AAgBD,8EAA8E;AAC9E,wBAAgB,kBAAkB,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,MAAM,CAWjE"}

package/dist/stats/weighted-stats.d.ts

@@ -0,0 +1,27 @@
+/**
+ * Weighted descriptive statistics (gap-analysis 2026-11-04, Batch A).
+ *
+ * Weights are non-negative frequency-style weights. Zero-weight points never
+ * contribute to a quantile.
+ */
+type Vec = readonly number[] | Float64Array;
+/** Weighted arithmetic mean: Σwᵢxᵢ / Σwᵢ. */
+export declare function weightedMean(x: Vec, w: Vec): number;
+/**
+ * Weighted variance. `ddof: 0` (default) gives Σw(x−m)²/Σw; `ddof: 1` divides
+ * by (Σw − Σw²/Σw), the reliability-weights unbiased form (equal weights give
+ * n/(n−1) scaling, and it stays finite for normalised weights).
+ */
+export declare function weightedVariance(x: Vec, w: Vec, opts?: {
+    ddof?: 0 | 1;
+}): number;
+/**
+ * Weighted quantile, lower rule: sort by x and return the first value whose
+ * cumulative weight reaches p·Σw. Points with zero weight are never returned
+ * unless every weight is zero (which is rejected).
+ */
+export declare function weightedQuantile(x: Vec, w: Vec, p: number): number;
+/** Weighted median (lower rule): `weightedQuantile(x, w, 0.5)`. */
+export declare function weightedMedian(x: Vec, w: Vec): number;
+export {};
+//# sourceMappingURL=weighted-stats.d.ts.map

package/dist/stats/weighted-stats.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"weighted-stats.d.ts","sourceRoot":"","sources":["../../src/stats/weighted-stats.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AACH,KAAK,GAAG,GAAG,SAAS,MAAM,EAAE,GAAG,YAAY,CAAC;AAqB5C,6CAA6C;AAC7C,wBAAgB,YAAY,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,MAAM,CAWnD;AAED;;;;GAIG;AACH,wBAAgB,gBAAgB,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,EAAE;IAAE,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAA;CAAE,GAAG,MAAM,CAgBhF;AAED;;;;GAIG;AACH,wBAAgB,gBAAgB,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAgBlE;AAED,mEAAmE;AACnE,wBAAgB,cAAc,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,MAAM,CAErD"}

package/dist/typed/algebra.d.ts

@@ -47,9 +47,9 @@ export declare function polyadd(a: number[], b: number[]): number[];
  * @param b - Second polynomial coefficients
  * @returns Product polynomial coefficients
  *
- * Below `WASM_POLY_THRESHOLD` coefficients, multiplies via the shared
- * `polyMulJS` kernel then strips high-end zeros with `trimPoly`. At or
- * above the threshold, uses `polyMulDispatch` (WASM when policy allows).
+ * Multiplies through `polyMulDispatch` (the shared `polyMulJS` kernel, or the
+ * WASM kernel at or above `WASM_POLY_THRESHOLD` coefficients when the policy
+ * allows), then strips high-end zeros with `trimPoly`.
  *
  * @example
  * ```typescript
@@ -88,6 +88,10 @@ export declare function polyder(coeffs: number[], n?: number): number[];
  * @param b - Second polynomial coefficients (index = power)
  * @returns GCD polynomial coefficients (monic); `[0]` when both inputs are zero
  *
+ * Integer coefficients, and rational coefficients whose common denominator is at most 1e6
+ * (checked to 1e-12 relative), are scaled to integers and solved exactly; other input
+ * uses floating-point Euclid.
+ *
  * @example
  * ```typescript
  * // GCD of x^2-1 and x-1 is x-1