@danielsimonjr/mathts-functions 0.12.1 → 0.13.0

package/dist/index.js

@@ -468,10 +468,11 @@ __export(typed_exports, {
 });
 
 // src/typed/arithmetic.ts
-import { mathTyped, Complex, Fraction, BigNumber, Dual } from "@danielsimonjr/mathts-core";
+import { mathTyped as mathTyped2, Complex, Fraction, BigNumber, Dual } from "@danielsimonjr/mathts-core";
 import { DenseMatrix as DenseMatrix2, backendManager as backendManager2, singularValues } from "@danielsimonjr/mathts-matrix";
 
 // src/factories/matrix-bridge.ts
+import { mathTyped } from "@danielsimonjr/mathts-core";
 import { DenseMatrix, backendManager } from "@danielsimonjr/mathts-matrix";
 void backendManager.initialize().catch(() => {
 });
@@ -490,6 +491,29 @@ function matrixIndexSelectors(index) {
   }
   return out;
 }
+function _adaptCallbackArity(callback) {
+  const t = mathTyped;
+  if (t.isTypedFunction(callback)) {
+    let arity;
+    return (...args) => {
+      if (arity === void 0) {
+        for (let k = args.length; k >= 1; k--) {
+          if (t.resolve(callback, args.slice(0, k))) {
+            arity = k;
+            break;
+          }
+        }
+        if (arity === void 0) arity = args.length;
+      }
+      return callback(...args.slice(0, arity));
+    };
+  }
+  const declared = callback.length;
+  if (declared >= 1 && declared < 3) {
+    return (...args) => callback(...args.slice(0, declared));
+  }
+  return callback;
+}
 var MathJSDenseMatrix = class _MathJSDenseMatrix {
   _data;
   _size;
@@ -564,27 +588,38 @@ var MathJSDenseMatrix = class _MathJSDenseMatrix {
   }
   /**
    * Map over all elements, producing a new matrix.
+   *
+   * Callbacks are invoked ARITY-APPROPRIATELY (mathjs contract): a typed or
+   * fixed-arity callback declared as `(value)` or `(value, index)` must not be
+   * force-fed 3 arguments — a typed 2-arg expression lambda (`f(x, i) = …`)
+   * previously threw "Too many arguments" out of the expression language.
    */
   map(callback) {
+    const cb = _adaptCallbackArity(callback);
     if (this._size.length === 1) {
-      const flat = this._data.map((val, i) => callback(val, [i], this));
-      const m = new _MathJSDenseMatrix(flat);
-      m._size = [...this._size];
+      const flat = this._data.map(
+        (val, i) => cb(val, [i], this)
+      );
+      const m = new _MathJSDenseMatrix({
+        data: flat,
+        size: [...this._size]
+      });
       return m;
     }
-    const result = this._data.map((row2, i) => row2.map((val, j) => callback(val, [i, j], this)));
+    const result = this._data.map((row2, i) => row2.map((val, j) => cb(val, [i, j], this)));
     return new _MathJSDenseMatrix(result);
   }
   /**
-   * Iterate over all elements.
+   * Iterate over all elements (callback invoked arity-appropriately; see map).
    */
   forEach(callback) {
+    const cb = _adaptCallbackArity(callback);
     if (this._size.length === 1) {
-      this._data.forEach((val, i) => callback(val, [i], this));
+      this._data.forEach((val, i) => cb(val, [i], this));
       return;
     }
     this._data.forEach((row2, i) => {
-      row2.forEach((val, j) => callback(val, [i, j], this));
+      row2.forEach((val, j) => cb(val, [i, j], this));
     });
   }
   /**
@@ -1333,7 +1368,7 @@ var unitAdd = (a, b) => unitAddSub(a, b, addOp);
 var unitSub = (a, b) => unitAddSub(a, b, (x, y) => subtract(x, y));
 var unitMul = (a, b) => asUnit(a).multiply(b);
 var unitDiv = (a, b) => asUnit(a).divide(b);
-var add = mathTyped("add", {
+var add = mathTyped2("add", {
   "number, number": (a, b) => a + b,
   "bigint, bigint": (a, b) => a + b,
   "Complex, Complex": (a, b) => a.add(b),
@@ -1385,7 +1420,7 @@ var add = mathTyped("add", {
     return result;
   }
 });
-var subtract = mathTyped("subtract", {
+var subtract = mathTyped2("subtract", {
   "number, number": (a, b) => a - b,
   "bigint, bigint": (a, b) => a - b,
   "Complex, Complex": (a, b) => a.subtract(b),
@@ -1409,7 +1444,7 @@ var subtract = mathTyped("subtract", {
     return result.result;
   }
 });
-var multiply = mathTyped("multiply", {
+var multiply = mathTyped2("multiply", {
   "number, number": (a, b) => a * b,
   "bigint, bigint": (a, b) => a * b,
   "Complex, Complex": (a, b) => a.multiply(b),
@@ -1488,7 +1523,7 @@ var multiply = mathTyped("multiply", {
     return result;
   }
 });
-var divide = mathTyped("divide", {
+var divide = mathTyped2("divide", {
   "number, number": (a, b) => a / b,
   "bigint, bigint": (a, b) => a / b,
   "Complex, Complex": (a, b) => a.divide(b),
@@ -1513,7 +1548,7 @@ var divide = mathTyped("divide", {
     return result.result;
   }
 });
-var unaryMinus = mathTyped("unaryMinus", {
+var unaryMinus = mathTyped2("unaryMinus", {
   number: (a) => -a,
   bigint: (a) => -a,
   Complex: (a) => a.negate(),
@@ -1525,7 +1560,7 @@ var unaryMinus = mathTyped("unaryMinus", {
     return result.result;
   }
 });
-var unaryPlus = mathTyped("unaryPlus", {
+var unaryPlus = mathTyped2("unaryPlus", {
   number: (a) => +a,
   bigint: (a) => a,
   Complex: (a) => a,
@@ -1534,7 +1569,7 @@ var unaryPlus = mathTyped("unaryPlus", {
   Float64Array: (a) => a
   // Identity, no copy needed
 });
-var abs = mathTyped("abs", {
+var abs = mathTyped2("abs", {
   number: (a) => Math.abs(a),
   bigint: (a) => a < 0n ? -a : a,
   Complex: (a) => a.abs(),
@@ -1550,7 +1585,7 @@ var abs = mathTyped("abs", {
     return result.result;
   }
 });
-var sign = mathTyped("sign", {
+var sign = mathTyped2("sign", {
   number: (a) => Math.sign(a),
   bigint: (a) => a > 0n ? 1n : a < 0n ? -1n : 0n,
   Fraction: (a) => new Fraction(BigInt(a.sign()), 1n),
@@ -1571,7 +1606,7 @@ var sign = mathTyped("sign", {
     return out;
   }
 });
-var pow = mathTyped("pow", {
+var pow = mathTyped2("pow", {
   "number, number": (a, b) => Math.pow(a, b),
   "bigint, bigint": (a, b) => a ** b,
   "Complex, number": (a, b) => a.pow(b),
@@ -1613,7 +1648,7 @@ var pow = mathTyped("pow", {
     return result.toArray();
   }
 });
-var sqrt = mathTyped("sqrt", {
+var sqrt = mathTyped2("sqrt", {
   number: (a) => {
     if (a < 0) {
       return new Complex(0, Math.sqrt(-a));
@@ -1629,7 +1664,7 @@ var sqrt = mathTyped("sqrt", {
     return result.result;
   }
 });
-var square = mathTyped("square", {
+var square = mathTyped2("square", {
   number: (a) => a * a,
   bigint: (a) => a * a,
   Complex: (a) => a.multiply(a),
@@ -1642,7 +1677,7 @@ var square = mathTyped("square", {
     return result.result;
   }
 });
-var cube = mathTyped("cube", {
+var cube = mathTyped2("cube", {
   number: (a) => a * a * a,
   bigint: (a) => a * a * a,
   Complex: (a) => a.multiply(a).multiply(a),
@@ -1660,7 +1695,7 @@ var cube = mathTyped("cube", {
     return out;
   }
 });
-var cbrt = mathTyped("cbrt", {
+var cbrt = mathTyped2("cbrt", {
   number: (a) => Math.cbrt(a),
   Complex: (a) => a.pow(1 / 3),
   BigNumber: (a) => a.cbrt(),
@@ -1676,12 +1711,12 @@ var cbrt = mathTyped("cbrt", {
     return out;
   }
 });
-var nthRoot = mathTyped("nthRoot", {
+var nthRoot = mathTyped2("nthRoot", {
   "number, number": (a, n) => Math.pow(a, 1 / n),
   "Complex, number": (a, n) => a.pow(1 / n),
   "BigNumber, number": (a, n) => a.pow(1 / n)
 });
-var exp = mathTyped("exp", {
+var exp = mathTyped2("exp", {
   number: (a) => Math.exp(a),
   Complex: (a) => a.exp(),
   BigNumber: (a) => a.exp(),
@@ -1694,7 +1729,7 @@ var exp = mathTyped("exp", {
     return result.result;
   }
 });
-var log = mathTyped("log", {
+var log = mathTyped2("log", {
   number: (a) => Math.log(a),
   Complex: (a) => a.log(),
   BigNumber: (a) => a.ln(),
@@ -1708,7 +1743,7 @@ var log = mathTyped("log", {
     return result.result;
   }
 });
-var log10 = mathTyped("log10", {
+var log10 = mathTyped2("log10", {
   number: (a) => Math.log10(a),
   Complex: (a) => a.log().divide(Complex.fromNumber(Math.LN10)),
   BigNumber: (a) => a.log10(),
@@ -1725,7 +1760,7 @@ var log10 = mathTyped("log10", {
     return out;
   }
 });
-var log2 = mathTyped("log2", {
+var log2 = mathTyped2("log2", {
   number: (a) => Math.log2(a),
   Complex: (a) => a.log().divide(Complex.fromNumber(Math.LN2)),
   BigNumber: (a) => a.log2(),
@@ -1742,7 +1777,7 @@ var log2 = mathTyped("log2", {
     return out;
   }
 });
-var log1p = mathTyped("log1p", {
+var log1p = mathTyped2("log1p", {
   number: (a) => Math.log1p(a),
   // Parallel array log1p
   Float64Array: async (a) => {
@@ -1757,7 +1792,7 @@ var log1p = mathTyped("log1p", {
     return out;
   }
 });
-var expm1 = mathTyped("expm1", {
+var expm1 = mathTyped2("expm1", {
   number: (a) => Math.expm1(a),
   // Parallel array expm1
   Float64Array: async (a) => {
@@ -1772,7 +1807,7 @@ var expm1 = mathTyped("expm1", {
     return out;
   }
 });
-var round = mathTyped("round", {
+var round = mathTyped2("round", {
   number: (a) => Math.round(a),
   Fraction: (a) => a.round(),
   BigNumber: (a) => a.round(),
@@ -1792,7 +1827,7 @@ var round = mathTyped("round", {
     return out;
   }
 });
-var floor = mathTyped("floor", {
+var floor = mathTyped2("floor", {
   number: (a) => Math.floor(a),
   Fraction: (a) => a.floor(),
   BigNumber: (a) => a.floor(),
@@ -1808,7 +1843,7 @@ var floor = mathTyped("floor", {
     return out;
   }
 });
-var ceil = mathTyped("ceil", {
+var ceil = mathTyped2("ceil", {
   number: (a) => Math.ceil(a),
   Fraction: (a) => a.ceil(),
   BigNumber: (a) => a.ceil(),
@@ -1824,7 +1859,7 @@ var ceil = mathTyped("ceil", {
     return out;
   }
 });
-var fix = mathTyped("fix", {
+var fix = mathTyped2("fix", {
   number: (a) => Math.trunc(a),
   Fraction: (a) => a.trunc(),
   BigNumber: (a) => a.trunc(),
@@ -1840,13 +1875,13 @@ var fix = mathTyped("fix", {
     return out;
   }
 });
-var mod = mathTyped("mod", {
+var mod = mathTyped2("mod", {
   "number, number": (a, b) => (a % b + b) % b,
   "bigint, bigint": (a, b) => (a % b + b) % b,
   "Fraction, Fraction": (a, b) => a.mod(b),
   "BigNumber, BigNumber": (a, b) => a.mod(b)
 });
-var gcd = mathTyped("gcd", {
+var gcd = mathTyped2("gcd", {
   "number, number": (a, b) => {
     a = Math.abs(Math.floor(a));
     b = Math.abs(Math.floor(b));
@@ -1878,7 +1913,7 @@ var gcd = mathTyped("gcd", {
     return x;
   }
 });
-var lcm = mathTyped("lcm", {
+var lcm = mathTyped2("lcm", {
   "number, number": (a, b) => {
     const g = gcd(a, b);
     return Math.abs(a * b) / g;
@@ -1895,7 +1930,7 @@ var lcm = mathTyped("lcm", {
     return a.abs().multiply(b.abs()).divide(g);
   }
 });
-var xgcd = mathTyped("xgcd", {
+var xgcd = mathTyped2("xgcd", {
   "number, number": (a, b) => {
     a = Math.floor(a);
     b = Math.floor(b);
@@ -1967,7 +2002,7 @@ function matrixNorm(data, p = "fro") {
   throw new Error(`Unsupported matrix norm p=${String(p)} (use 'fro', 1, 2, or Infinity)`);
 }
 var matrixLike = (m) => m.toArray();
-var norm = mathTyped("norm", {
+var norm = mathTyped2("norm", {
   number: (a) => Math.abs(a),
   Complex: (a) => a.abs(),
   BigNumber: (a) => a.abs(),
@@ -1994,7 +2029,7 @@ var norm = mathTyped("norm", {
     return Math.sqrt(s);
   }
 });
-var sinh = mathTyped("sinh", {
+var sinh = mathTyped2("sinh", {
   number: (a) => Math.sinh(a),
   Complex: (a) => a.sinh(),
   BigNumber: (a) => a.sinh(),
@@ -2011,7 +2046,7 @@ var sinh = mathTyped("sinh", {
     return out;
   }
 });
-var cosh = mathTyped("cosh", {
+var cosh = mathTyped2("cosh", {
   number: (a) => Math.cosh(a),
   Complex: (a) => a.cosh(),
   BigNumber: (a) => a.cosh(),
@@ -2026,7 +2061,7 @@ var cosh = mathTyped("cosh", {
     return out;
   }
 });
-var tanh = mathTyped("tanh", {
+var tanh = mathTyped2("tanh", {
   number: (a) => Math.tanh(a),
   Complex: (a) => a.tanh(),
   BigNumber: (a) => a.tanh(),
@@ -2043,7 +2078,7 @@ var tanh = mathTyped("tanh", {
     return out;
   }
 });
-var equal = mathTyped("equal", {
+var equal = mathTyped2("equal", {
   "number, number": (a, b) => a === b,
   "bigint, bigint": (a, b) => a === b,
   "Complex, Complex": (a, b) => a.equals(b),
@@ -2051,42 +2086,42 @@ var equal = mathTyped("equal", {
   "BigNumber, BigNumber": (a, b) => a.equals(b),
   "Unit, Unit": (a, b) => unitSameDimension(a, b) && asUnit(a).value === asUnit(b).value
 });
-var smaller = mathTyped("smaller", {
+var smaller = mathTyped2("smaller", {
   "number, number": (a, b) => a < b,
   "bigint, bigint": (a, b) => a < b,
   "Fraction, Fraction": (a, b) => a.compare(b) < 0,
   "BigNumber, BigNumber": (a, b) => a.compare(b) < 0,
   "Unit, Unit": (a, b) => unitCmp(a, b) < 0
 });
-var larger = mathTyped("larger", {
+var larger = mathTyped2("larger", {
   "number, number": (a, b) => a > b,
   "bigint, bigint": (a, b) => a > b,
   "Fraction, Fraction": (a, b) => a.compare(b) > 0,
   "BigNumber, BigNumber": (a, b) => a.compare(b) > 0,
   "Unit, Unit": (a, b) => unitCmp(a, b) > 0
 });
-var smallerEq = mathTyped("smallerEq", {
+var smallerEq = mathTyped2("smallerEq", {
   "number, number": (a, b) => a <= b,
   "bigint, bigint": (a, b) => a <= b,
   "Fraction, Fraction": (a, b) => a.compare(b) <= 0,
   "BigNumber, BigNumber": (a, b) => a.compare(b) <= 0,
   "Unit, Unit": (a, b) => unitCmp(a, b) <= 0
 });
-var largerEq = mathTyped("largerEq", {
+var largerEq = mathTyped2("largerEq", {
   "number, number": (a, b) => a >= b,
   "bigint, bigint": (a, b) => a >= b,
   "Fraction, Fraction": (a, b) => a.compare(b) >= 0,
   "BigNumber, BigNumber": (a, b) => a.compare(b) >= 0,
   "Unit, Unit": (a, b) => unitCmp(a, b) >= 0
 });
-var compare = mathTyped("compare", {
+var compare = mathTyped2("compare", {
   "number, number": (a, b) => a > b ? 1 : a < b ? -1 : 0,
   "bigint, bigint": (a, b) => a > b ? 1 : a < b ? -1 : 0,
   "Fraction, Fraction": (a, b) => a.compare(b),
   "BigNumber, BigNumber": (a, b) => a.compare(b),
   "Unit, Unit": (a, b) => unitCmp(a, b)
 });
-var min = mathTyped("min", {
+var min = mathTyped2("min", {
   "bigint, bigint": (a, b) => a < b ? a : b,
   "Fraction, Fraction": (a, b) => a.compare(b) < 0 ? a : b,
   "BigNumber, BigNumber": (a, b) => a.compare(b) < 0 ? a : b,
@@ -2105,7 +2140,7 @@ var min = mathTyped("min", {
     return m;
   }
 });
-var max = mathTyped("max", {
+var max = mathTyped2("max", {
   "bigint, bigint": (a, b) => a > b ? a : b,
   "Fraction, Fraction": (a, b) => a.compare(b) > 0 ? a : b,
   "BigNumber, BigNumber": (a, b) => a.compare(b) > 0 ? a : b,
@@ -2123,7 +2158,7 @@ var max = mathTyped("max", {
     return m;
   }
 });
-var sum = mathTyped("sum", {
+var sum = mathTyped2("sum", {
   Array: (arr7) => {
     let total = 0;
     for (let i = 0; i < arr7.length; i++) {
@@ -2137,7 +2172,7 @@ var sum = mathTyped("sum", {
     return result.result;
   }
 });
-var mean = mathTyped("mean", {
+var mean = mathTyped2("mean", {
   Array: (arr7) => {
     if (arr7.length === 0) return NaN;
     let total = 0;
@@ -2175,7 +2210,7 @@ function m2OfArray(arr7) {
   }
   return m2;
 }
-var variance = mathTyped("variance", {
+var variance = mathTyped2("variance", {
   Array: (arr7) => varianceFromM2(m2OfArray(arr7), arr7.length, "unbiased"),
   "Array, string": (arr7, norm4) => varianceFromM2(m2OfArray(arr7), arr7.length, norm4),
   // Parallel Float64Array variance. computePool returns POPULATION variance
@@ -2189,7 +2224,7 @@ var variance = mathTyped("variance", {
     return varianceFromM2(result.result.variance * a.length, a.length, norm4);
   }
 });
-var std = mathTyped("std", {
+var std = mathTyped2("std", {
   Array: (arr7) => Math.sqrt(varianceFromM2(m2OfArray(arr7), arr7.length, "unbiased")),
   "Array, string": (arr7, norm4) => Math.sqrt(varianceFromM2(m2OfArray(arr7), arr7.length, norm4)),
   // Parallel Float64Array std — derive from the (renormalized) variance so the
@@ -2205,7 +2240,7 @@ var std = mathTyped("std", {
     );
   }
 });
-var dot = mathTyped("dot", {
+var dot = mathTyped2("dot", {
   "Array, Array": (a, b) => {
     if (a.length !== b.length) {
       throw new Error(`Vector lengths must match: ${a.length} vs ${b.length}`);
@@ -2320,9 +2355,9 @@ var typedArithmetic = {
 };
 
 // src/typed/trigonometry.ts
-import { mathTyped as mathTyped2, Complex as Complex2, BigNumber as BigNumber2 } from "@danielsimonjr/mathts-core";
+import { mathTyped as mathTyped3, Complex as Complex2, BigNumber as BigNumber2 } from "@danielsimonjr/mathts-core";
 import { computePool as computePool2 } from "@danielsimonjr/mathts-parallel";
-var sin = mathTyped2("sin", {
+var sin = mathTyped3("sin", {
   number: (a) => Math.sin(a),
   Complex: (a) => a.sin(),
   BigNumber: (a) => a.sin(),
@@ -2335,7 +2370,7 @@ var sin = mathTyped2("sin", {
     return result.result;
   }
 });
-var cos = mathTyped2("cos", {
+var cos = mathTyped3("cos", {
   number: (a) => Math.cos(a),
   Complex: (a) => a.cos(),
   BigNumber: (a) => a.cos(),
@@ -2348,7 +2383,7 @@ var cos = mathTyped2("cos", {
     return result.result;
   }
 });
-var tan = mathTyped2("tan", {
+var tan = mathTyped3("tan", {
   number: (a) => Math.tan(a),
   Complex: (a) => a.tan(),
   BigNumber: (a) => a.tan(),
@@ -2361,7 +2396,7 @@ var tan = mathTyped2("tan", {
     return result.result;
   }
 });
-var csc = mathTyped2("csc", {
+var csc = mathTyped3("csc", {
   number: (a) => 1 / Math.sin(a),
   Complex: (a) => a.sin().inverse(),
   BigNumber: (a) => BigNumber2.fromNumber(1).divide(a.sin()),
@@ -2378,7 +2413,7 @@ var csc = mathTyped2("csc", {
     return out;
   }
 });
-var sec = mathTyped2("sec", {
+var sec = mathTyped3("sec", {
   number: (a) => 1 / Math.cos(a),
   Complex: (a) => a.cos().inverse(),
   BigNumber: (a) => BigNumber2.fromNumber(1).divide(a.cos()),
@@ -2395,7 +2430,7 @@ var sec = mathTyped2("sec", {
     return out;
   }
 });
-var cot = mathTyped2("cot", {
+var cot = mathTyped3("cot", {
   number: (a) => 1 / Math.tan(a),
   Complex: (a) => a.tan().inverse(),
   BigNumber: (a) => BigNumber2.fromNumber(1).divide(a.tan()),
@@ -2412,7 +2447,7 @@ var cot = mathTyped2("cot", {
     return out;
   }
 });
-var asin = mathTyped2("asin", {
+var asin = mathTyped3("asin", {
   number: (a) => {
     if (a < -1 || a > 1) {
       return new Complex2(a, 0).asin();
@@ -2432,7 +2467,7 @@ var asin = mathTyped2("asin", {
     return out;
   }
 });
-var acos = mathTyped2("acos", {
+var acos = mathTyped3("acos", {
   number: (a) => {
     if (a < -1 || a > 1) {
       return new Complex2(a, 0).acos();
@@ -2452,7 +2487,7 @@ var acos = mathTyped2("acos", {
     return out;
   }
 });
-var atan = mathTyped2("atan", {
+var atan = mathTyped3("atan", {
   number: (a) => Math.atan(a),
   Complex: (a) => a.atan(),
   BigNumber: (a) => a.atan(),
@@ -2469,29 +2504,29 @@ var atan = mathTyped2("atan", {
     return out;
   }
 });
-var atan2 = mathTyped2("atan2", {
+var atan2 = mathTyped3("atan2", {
   "number, number": (y, x) => Math.atan2(y, x),
   "BigNumber, BigNumber": (y, x) => y.atan2(x)
 });
-var acsc = mathTyped2("acsc", {
+var acsc = mathTyped3("acsc", {
   number: (a) => Math.asin(1 / a),
   Complex: (a) => a.inverse().asin(),
   // BigNumber: acsc(a) = asin(1/a) — matches the forward csc's BigNumber path.
   BigNumber: (a) => BigNumber2.fromNumber(1).divide(a).asin()
 });
-var asec = mathTyped2("asec", {
+var asec = mathTyped3("asec", {
   number: (a) => Math.acos(1 / a),
   Complex: (a) => a.inverse().acos(),
   // BigNumber: asec(a) = acos(1/a).
   BigNumber: (a) => BigNumber2.fromNumber(1).divide(a).acos()
 });
-var acot = mathTyped2("acot", {
+var acot = mathTyped3("acot", {
   number: (a) => Math.atan(1 / a),
   Complex: (a) => a.inverse().atan(),
   // BigNumber: acot(a) = atan(1/a).
   BigNumber: (a) => BigNumber2.fromNumber(1).divide(a).atan()
 });
-var asinh = mathTyped2("asinh", {
+var asinh = mathTyped3("asinh", {
   number: (a) => Math.asinh(a),
   Complex: (a) => a.asinh(),
   BigNumber: (a) => a.asinh(),
@@ -2506,7 +2541,7 @@ var asinh = mathTyped2("asinh", {
     return out;
   }
 });
-var acosh = mathTyped2("acosh", {
+var acosh = mathTyped3("acosh", {
   number: (a) => Math.acosh(a),
   Complex: (a) => a.acosh(),
   BigNumber: (a) => a.acosh(),
@@ -2521,7 +2556,7 @@ var acosh = mathTyped2("acosh", {
     return out;
   }
 });
-var atanh = mathTyped2("atanh", {
+var atanh = mathTyped3("atanh", {
   number: (a) => Math.atanh(a),
   Complex: (a) => a.atanh(),
   BigNumber: (a) => a.atanh(),
@@ -2538,13 +2573,13 @@ var atanh = mathTyped2("atanh", {
     return out;
   }
 });
-var toRadians = mathTyped2("toRadians", {
+var toRadians = mathTyped3("toRadians", {
   number: (deg) => deg * Math.PI / 180
 });
-var toDegrees = mathTyped2("toDegrees", {
+var toDegrees = mathTyped3("toDegrees", {
   number: (rad) => rad * 180 / Math.PI
 });
-var hypot = mathTyped2("hypot", {
+var hypot = mathTyped3("hypot", {
   number: (a) => Math.abs(a),
   "number, number": (a, b) => Math.hypot(a, b),
   "number, number, ...number": (a, b, rest) => Math.hypot(a, b, ...rest),
@@ -2577,12 +2612,12 @@ var typedTrigonometry = {
 };
 
 // src/typed/statistics.ts
-import { mathTyped as mathTyped3 } from "@danielsimonjr/mathts-core";
+import { mathTyped as mathTyped4 } from "@danielsimonjr/mathts-core";
 import { computePool as computePool3 } from "@danielsimonjr/mathts-parallel";
 function toFloat64Array(arr7) {
   return new Float64Array(arr7);
 }
-var parallelStatSum = mathTyped3("parallelStatSum", {
+var parallelStatSum = mathTyped4("parallelStatSum", {
   // Float64Array - parallel execution
   Float64Array: async (data) => {
     const result = await computePool3.sum(data);
@@ -2599,7 +2634,7 @@ var parallelStatSum = mathTyped3("parallelStatSum", {
   "number, number, number": (a, b, c) => a + b + c,
   "number, number, number, number": (a, b, c, d) => a + b + c + d
 });
-var parallelStatMean = mathTyped3("parallelStatMean", {
+var parallelStatMean = mathTyped4("parallelStatMean", {
   // Float64Array - parallel execution
   Float64Array: async (data) => {
     const result = await computePool3.mean(data);
@@ -2617,7 +2652,7 @@ var parallelStatMean = mathTyped3("parallelStatMean", {
   "number, number, number": (a, b, c) => (a + b + c) / 3,
   "number, number, number, number": (a, b, c, d) => (a + b + c + d) / 4
 });
-var parallelStatVariance = mathTyped3("parallelStatVariance", {
+var parallelStatVariance = mathTyped4("parallelStatVariance", {
   // Float64Array with default normalization (unbiased)
   Float64Array: async (data) => {
     const result = await computePool3.variance(data);
@@ -2664,7 +2699,7 @@ var parallelStatVariance = mathTyped3("parallelStatVariance", {
     return ((a - mean6) ** 2 + (b - mean6) ** 2 + (c - mean6) ** 2) / 2;
   }
 });
-var parallelStatStd = mathTyped3("parallelStatStd", {
+var parallelStatStd = mathTyped4("parallelStatStd", {
   // Float64Array with default normalization (unbiased)
   Float64Array: async (data) => {
     const result = await computePool3.std(data);
@@ -2690,7 +2725,7 @@ var parallelStatStd = mathTyped3("parallelStatStd", {
     return parallelStatStd(data, normalization);
   }
 });
-var parallelStatMin = mathTyped3("parallelStatMin", {
+var parallelStatMin = mathTyped4("parallelStatMin", {
   // Float64Array - parallel execution
   Float64Array: async (data) => {
     const result = await computePool3.min(data);
@@ -2710,7 +2745,7 @@ var parallelStatMin = mathTyped3("parallelStatMin", {
   // Four numbers
   "number, number, number, number": (a, b, c, d) => Math.min(a, b, c, d)
 });
-var parallelStatMax = mathTyped3("parallelStatMax", {
+var parallelStatMax = mathTyped4("parallelStatMax", {
   // Float64Array - parallel execution
   Float64Array: async (data) => {
     const result = await computePool3.max(data);
@@ -2730,7 +2765,7 @@ var parallelStatMax = mathTyped3("parallelStatMax", {
   // Four numbers
   "number, number, number, number": (a, b, c, d) => Math.max(a, b, c, d)
 });
-var parallelStatMinMax = mathTyped3("parallelStatMinMax", {
+var parallelStatMinMax = mathTyped4("parallelStatMinMax", {
   // Float64Array - parallel execution
   Float64Array: async (data) => {
     const result = await computePool3.minMax(data);
@@ -2743,7 +2778,7 @@ var parallelStatMinMax = mathTyped3("parallelStatMinMax", {
     return result.result;
   }
 });
-var parallelStatMedian = mathTyped3("parallelStatMedian", {
+var parallelStatMedian = mathTyped4("parallelStatMedian", {
   // Float64Array
   Float64Array: async (data) => {
     const n = data.length;
@@ -2769,7 +2804,7 @@ var parallelStatMedian = mathTyped3("parallelStatMedian", {
     return arr7[1];
   }
 });
-var parallelStatMode = mathTyped3("parallelStatMode", {
+var parallelStatMode = mathTyped4("parallelStatMode", {
   // Number array
   Array: (arr7) => {
     if (arr7.length === 0) return [];
@@ -2792,7 +2827,7 @@ var parallelStatMode = mathTyped3("parallelStatMode", {
     return parallelStatMode(Array.from(data));
   }
 });
-var parallelStatProd = mathTyped3("parallelStatProd", {
+var parallelStatProd = mathTyped4("parallelStatProd", {
   // Number array
   Array: (arr7) => {
     if (arr7.length === 0) return 1;
@@ -2819,7 +2854,7 @@ var parallelStatProd = mathTyped3("parallelStatProd", {
   "number, number, number": (a, b, c) => a * b * c,
   "number, number, number, number": (a, b, c, d) => a * b * c * d
 });
-var parallelStatNorm = mathTyped3("parallelStatNorm", {
+var parallelStatNorm = mathTyped4("parallelStatNorm", {
   // Float64Array - parallel execution
   Float64Array: async (data) => {
     const result = await computePool3.norm(data);
@@ -2872,7 +2907,7 @@ var parallelStatNorm = mathTyped3("parallelStatNorm", {
     return parallelStatNorm(toFloat64Array(arr7), p);
   }
 });
-var parallelStatDistance = mathTyped3("parallelStatDistance", {
+var parallelStatDistance = mathTyped4("parallelStatDistance", {
   // Float64Array - parallel execution
   "Float64Array, Float64Array": async (a, b) => {
     const result = await computePool3.distance(a, b);
@@ -2888,7 +2923,7 @@ var parallelStatDistance = mathTyped3("parallelStatDistance", {
     return Math.sqrt((x2 - x1) ** 2 + (y2 - y1) ** 2);
   }
 });
-var parallelStatCorr = mathTyped3("parallelStatCorr", {
+var parallelStatCorr = mathTyped4("parallelStatCorr", {
   // Float64Array pairs - compute correlation
   "Float64Array, Float64Array": async (x, y) => {
     const n = Math.min(x.length, y.length);
@@ -2912,7 +2947,7 @@ var parallelStatCorr = mathTyped3("parallelStatCorr", {
     return parallelStatCorr(toFloat64Array(x), toFloat64Array(y));
   }
 });
-var parallelStatMAD = mathTyped3("parallelStatMAD", {
+var parallelStatMAD = mathTyped4("parallelStatMAD", {
   // Float64Array
   Float64Array: async (data) => {
     const n = data.length;
@@ -2930,7 +2965,7 @@ var parallelStatMAD = mathTyped3("parallelStatMAD", {
     return parallelStatMAD(toFloat64Array(arr7));
   }
 });
-var parallelStatCumsum = mathTyped3("parallelStatCumsum", {
+var parallelStatCumsum = mathTyped4("parallelStatCumsum", {
   // Float64Array
   Float64Array: (data) => {
     const result = new Float64Array(data.length);
@@ -2952,7 +2987,7 @@ var parallelStatCumsum = mathTyped3("parallelStatCumsum", {
     return result;
   }
 });
-var parallelStatQuantile = mathTyped3("parallelStatQuantile", {
+var parallelStatQuantile = mathTyped4("parallelStatQuantile", {
   // Float64Array with quantile value
   "Float64Array, number": (data, q) => {
     if (q < 0 || q > 1) {
@@ -2982,7 +3017,7 @@ function parallelStatPercentile(data, p) {
   if (p < 0 || p > 100) throw new Error("Percentile must be between 0 and 100");
   return parallelStatQuantile(data, p / 100);
 }
-var parallelStatHistogram = mathTyped3("parallelStatHistogram", {
+var parallelStatHistogram = mathTyped4("parallelStatHistogram", {
   // Float64Array with number of bins
   "Float64Array, number": async (data, bins) => {
     const result = await computePool3.histogram(data, bins);
@@ -3109,7 +3144,7 @@ async function terminateStatistics() {
 }
 
 // src/typed/signal.ts
-import { mathTyped as mathTyped4 } from "@danielsimonjr/mathts-core";
+import { mathTyped as mathTyped5 } from "@danielsimonjr/mathts-core";
 import { computePool as computePool4 } from "@danielsimonjr/mathts-parallel";
 
 // src/wasm/signal/wasm-bridge.ts
@@ -3549,7 +3584,7 @@ async function fourStepFFT(real, imag, inverse) {
   }
   return { real: outReal, imag: outImag };
 }
-var parallelFFT = mathTyped4("parallelFFT", {
+var parallelFFT = mathTyped5("parallelFFT", {
   // FFT of number array
   Array: (signal) => {
     const n = signal.length;
@@ -3575,7 +3610,7 @@ var parallelFFT = mathTyped4("parallelFFT", {
     return { ...result, originalLength: n };
   }
 });
-var parallelIFFT = mathTyped4("parallelIFFT", {
+var parallelIFFT = mathTyped5("parallelIFFT", {
   // IFFT from real/imag arrays
   "Float64Array, Float64Array": async (real, imag) => {
     return computePool4.shouldParallelize(real.length) ? fourStepFFT(real, imag, true) : fftCoreFloat64(real, imag, true);
@@ -3585,7 +3620,7 @@ var parallelIFFT = mathTyped4("parallelIFFT", {
     return computePool4.shouldParallelize(spectrum.real.length) ? fourStepFFT(spectrum.real, spectrum.imag, true) : fftCoreFloat64(spectrum.real, spectrum.imag, true);
   }
 });
-var parallelFFTMagnitude = mathTyped4("parallelFFTMagnitude", {
+var parallelFFTMagnitude = mathTyped5("parallelFFTMagnitude", {
   "Float64Array, Float64Array": async (real, imag) => {
     if (computePool4.shouldParallelize(real.length)) {
       const r = await computePool4.applyKernel2(
@@ -3605,7 +3640,7 @@ var parallelFFTMagnitude = mathTyped4("parallelFFTMagnitude", {
     return parallelFFTMagnitude(spectrum.real, spectrum.imag);
   }
 });
-var parallelFFTPower = mathTyped4("parallelFFTPower", {
+var parallelFFTPower = mathTyped5("parallelFFTPower", {
   "Float64Array, Float64Array": async (real, imag) => {
     if (computePool4.shouldParallelize(real.length)) {
       const r = await computePool4.applyKernel2(real, imag, "(re, im) => re * re + im * im");
@@ -3621,7 +3656,7 @@ var parallelFFTPower = mathTyped4("parallelFFTPower", {
     return parallelFFTPower(spectrum.real, spectrum.imag);
   }
 });
-var parallelConv = mathTyped4("parallelConv", {
+var parallelConv = mathTyped5("parallelConv", {
   // Convolution of two Float64Arrays
   "Float64Array, Float64Array": async (x, h) => {
     const n = x.length;
@@ -3671,7 +3706,7 @@ var parallelConv = mathTyped4("parallelConv", {
     return parallelConv(new Float64Array(x), new Float64Array(h));
   }
 });
-var parallelXCorr = mathTyped4("parallelXCorr", {
+var parallelXCorr = mathTyped5("parallelXCorr", {
   "Float64Array, Float64Array": async (x, h) => {
     const hReversed = new Float64Array(h.length);
     for (let i = 0; i < h.length; i++) {
@@ -3683,7 +3718,7 @@ var parallelXCorr = mathTyped4("parallelXCorr", {
     return parallelXCorr(new Float64Array(x), new Float64Array(h));
   }
 });
-var parallelAutoCorr = mathTyped4("parallelAutoCorr", {
+var parallelAutoCorr = mathTyped5("parallelAutoCorr", {
   Float64Array: async (x) => {
     return parallelXCorr(x, x);
   },
@@ -4444,7 +4479,7 @@ async function terminateSignal() {
 }
 
 // src/typed/bitwise.ts
-import { mathTyped as mathTyped5, BigNumber as BigNumber3 } from "@danielsimonjr/mathts-core";
+import { mathTyped as mathTyped6, BigNumber as BigNumber3 } from "@danielsimonjr/mathts-core";
 import { computePool as computePool5 } from "@danielsimonjr/mathts-parallel";
 
 // src/wasm/bitwise/wasm-bridge.ts
@@ -4523,7 +4558,7 @@ function logShiftBigInt(x, n) {
   const shifted = masked >> n;
   return shifted >= 0x80000000n ? shifted - 0x100000000n : shifted;
 }
-var bitAnd = mathTyped5("bitAnd", {
+var bitAnd = mathTyped6("bitAnd", {
   "number, number": (a, b) => {
     if (!Number.isInteger(a) || !Number.isInteger(b)) {
       throw new Error("Integers expected in function bitAnd");
@@ -4553,7 +4588,7 @@ var bitAnd = mathTyped5("bitAnd", {
     a & b
   )
 });
-var bitOr = mathTyped5("bitOr", {
+var bitOr = mathTyped6("bitOr", {
   "number, number": (a, b) => {
     if (!Number.isInteger(a) || !Number.isInteger(b)) {
       throw new Error("Integers expected in function bitOr");
@@ -4580,7 +4615,7 @@ var bitOr = mathTyped5("bitOr", {
     a | b
   )
 });
-var bitXor = mathTyped5("bitXor", {
+var bitXor = mathTyped6("bitXor", {
   "number, number": (a, b) => {
     if (!Number.isInteger(a) || !Number.isInteger(b)) {
       throw new Error("Integers expected in function bitXor");
@@ -4607,7 +4642,7 @@ var bitXor = mathTyped5("bitXor", {
     a ^ b
   )
 });
-var bitNot = mathTyped5("bitNot", {
+var bitNot = mathTyped6("bitNot", {
   number: (a) => {
     if (!Number.isInteger(a)) {
       throw new Error("Integer expected in function bitNot");
@@ -4628,7 +4663,7 @@ var bitNot = mathTyped5("bitNot", {
     return result.result;
   }
 });
-var leftShift = mathTyped5("leftShift", {
+var leftShift = mathTyped6("leftShift", {
   "number, number": (a, b) => {
     if (!Number.isInteger(a) || !Number.isInteger(b)) {
       throw new Error("Integers expected in function leftShift");
@@ -4664,7 +4699,7 @@ var leftShift = mathTyped5("leftShift", {
     a << b
   )
 });
-var rightArithShift = mathTyped5("rightArithShift", {
+var rightArithShift = mathTyped6("rightArithShift", {
   "number, number": (a, b) => {
     if (!Number.isInteger(a) || !Number.isInteger(b)) {
       throw new Error("Integers expected in function rightArithShift");
@@ -4700,7 +4735,7 @@ var rightArithShift = mathTyped5("rightArithShift", {
     a >> b
   )
 });
-var rightLogShift = mathTyped5("rightLogShift", {
+var rightLogShift = mathTyped6("rightLogShift", {
   "number, number": (a, b) => {
     if (!Number.isInteger(a) || !Number.isInteger(b)) {
       throw new Error("Integers expected in function rightLogShift");
@@ -4747,7 +4782,7 @@ var typedBitwise = {
 };
 
 // src/typed/logical.ts
-import { mathTyped as mathTyped6 } from "@danielsimonjr/mathts-core";
+import { mathTyped as mathTyped7 } from "@danielsimonjr/mathts-core";
 function isTruthyNumber(x) {
   return x !== 0 && !Number.isNaN(x);
 }
@@ -4760,14 +4795,14 @@ function isTruthyComplex(x) {
 function isTruthy(x) {
   return !!x;
 }
-var not = mathTyped6("not", {
+var not = mathTyped7("not", {
   any: (x) => !isTruthy(x),
   number: (x) => !isTruthyNumber(x),
   bigint: (x) => x === 0n,
   BigNumber: (x) => !isTruthyBigNumber(x),
   Complex: (x) => !isTruthyComplex(x)
 });
-var and = mathTyped6("and", {
+var and = mathTyped7("and", {
   "any, any": (a, b) => isTruthy(a) && isTruthy(b),
   "number, number": (a, b) => isTruthyNumber(a) && isTruthyNumber(b),
   "bigint, bigint": (a, b) => a !== 0n && b !== 0n,
@@ -4781,7 +4816,7 @@ var and = mathTyped6("and", {
     return isTruthy(a) && isTruthy(b) && restArr.every(isTruthy);
   }
 });
-var or = mathTyped6("or", {
+var or = mathTyped7("or", {
   "any, any": (a, b) => isTruthy(a) || isTruthy(b),
   "number, number": (a, b) => isTruthyNumber(a) || isTruthyNumber(b),
   "bigint, bigint": (a, b) => a !== 0n || b !== 0n,
@@ -4795,7 +4830,7 @@ var or = mathTyped6("or", {
     return isTruthy(a) || isTruthy(b) || restArr.some(isTruthy);
   }
 });
-var xor = mathTyped6("xor", {
+var xor = mathTyped7("xor", {
   "any, any": (a, b) => isTruthy(a) !== isTruthy(b),
   "number, number": (a, b) => isTruthyNumber(a) !== isTruthyNumber(b),
   "bigint, bigint": (a, b) => a !== 0n !== (b !== 0n),
@@ -4813,7 +4848,7 @@ var xor = mathTyped6("xor", {
     return result;
   }
 });
-var nullish = mathTyped6("nullish", {
+var nullish = mathTyped7("nullish", {
   // Concrete numeric types: always return a (never null/undefined)
   "number, any": (a, _b) => a,
   "bigint, any": (a, _b) => a,
@@ -4839,13 +4874,13 @@ var typedLogical = {
 };
 
 // src/typed/complex.ts
-import { mathTyped as mathTyped7, BigNumber as BigNumber5 } from "@danielsimonjr/mathts-core";
+import { mathTyped as mathTyped8, BigNumber as BigNumber5 } from "@danielsimonjr/mathts-core";
 function deepMapArray(arr7, fn) {
   return arr7.map(
     (item) => Array.isArray(item) ? deepMapArray(item, fn) : fn(item)
   );
 }
-var arg = mathTyped7("arg", {
+var arg = mathTyped8("arg", {
   number: (x) => Math.atan2(0, x),
   bigint: (x) => {
     if (x > 0n) return 0;
@@ -4861,21 +4896,21 @@ var arg = mathTyped7("arg", {
   Complex: (x) => x.arg(),
   Array: (x) => deepMapArray(x, (item) => arg(item))
 });
-var conj = mathTyped7("conj", {
+var conj = mathTyped8("conj", {
   number: (x) => x,
   bigint: (x) => x,
   BigNumber: (x) => x,
   Complex: (x) => x.conjugate(),
   Array: (x) => deepMapArray(x, (item) => conj(item))
 });
-var im = mathTyped7("im", {
+var im = mathTyped8("im", {
   number: (_x) => 0,
   bigint: (_x) => 0,
   BigNumber: (_x) => BigNumber5.fromNumber(0),
   Complex: (x) => x.im,
   Array: (x) => deepMapArray(x, (item) => im(item))
 });
-var re = mathTyped7("re", {
+var re = mathTyped8("re", {
   number: (x) => x,
   bigint: (x) => x,
   BigNumber: (x) => x,
@@ -4890,7 +4925,7 @@ var typedComplex = {
 };
 
 // src/typed/set.ts
-import { mathTyped as mathTyped8 } from "@danielsimonjr/mathts-core";
+import { mathTyped as mathTyped9 } from "@danielsimonjr/mathts-core";
 function flattenArray(arr7) {
   const result = [];
   const stack = [...arr7];
@@ -4946,7 +4981,7 @@ function identify(sorted) {
   }
   return result;
 }
-var setUnion = mathTyped8("setUnion", {
+var setUnion = mathTyped9("setUnion", {
   "Array, Array": (a1, a2) => {
     const b1 = toFlatArray(a1);
     const b2 = toFlatArray(a2);
@@ -4958,7 +4993,7 @@ var setUnion = mathTyped8("setUnion", {
     return [...diff1, ...inter, ...diff2];
   }
 });
-var setIntersect = mathTyped8("setIntersect", {
+var setIntersect = mathTyped9("setIntersect", {
   "Array, Array": (a1, a2) => {
     const b1 = toFlatArray(a1);
     const b2 = toFlatArray(a2);
@@ -4966,7 +5001,7 @@ var setIntersect = mathTyped8("setIntersect", {
     return _intersect(b1, b2);
   }
 });
-var setDifference = mathTyped8("setDifference", {
+var setDifference = mathTyped9("setDifference", {
   "Array, Array": (a1, a2) => {
     const b1 = toFlatArray(a1);
     const b2 = toFlatArray(a2);
@@ -4975,7 +5010,7 @@ var setDifference = mathTyped8("setDifference", {
     return _difference(b1, b2);
   }
 });
-var setSymDifference = mathTyped8("setSymDifference", {
+var setSymDifference = mathTyped9("setSymDifference", {
   "Array, Array": (a1, a2) => {
     const b1 = toFlatArray(a1);
     const b2 = toFlatArray(a2);
@@ -4984,7 +5019,7 @@ var setSymDifference = mathTyped8("setSymDifference", {
     return [..._difference(b1, b2), ..._difference(b2, b1)];
   }
 });
-var setIsSubset = mathTyped8("setIsSubset", {
+var setIsSubset = mathTyped9("setIsSubset", {
   "Array, Array": (a1, a2) => {
     const b1 = toFlatArray(a1);
     const b2 = toFlatArray(a2);
@@ -5001,7 +5036,7 @@ var setIsSubset = mathTyped8("setIsSubset", {
     return true;
   }
 });
-var setMultiplicity = mathTyped8("setMultiplicity", {
+var setMultiplicity = mathTyped9("setMultiplicity", {
   "any, Array": (e, a) => {
     const b = toFlatArray(a);
     if (b.length === 0) return 0;
@@ -5012,7 +5047,7 @@ var setMultiplicity = mathTyped8("setMultiplicity", {
     return count2;
   }
 });
-var setPowerset = mathTyped8("setPowerset", {
+var setPowerset = mathTyped9("setPowerset", {
   Array: (a) => {
     const b = toFlatArray(a).sort(naturalCompare);
     if (b.length === 0) return [];
@@ -5031,7 +5066,7 @@ var setPowerset = mathTyped8("setPowerset", {
     return result;
   }
 });
-var setDistinct = mathTyped8("setDistinct", {
+var setDistinct = mathTyped9("setDistinct", {
   Array: (a) => {
     const b = toFlatArray(a).sort(naturalCompare);
     if (b.length === 0) return [];
@@ -5044,7 +5079,7 @@ var setDistinct = mathTyped8("setDistinct", {
     return result;
   }
 });
-var setSize = mathTyped8("setSize", {
+var setSize = mathTyped9("setSize", {
   Array: (a) => toFlatArray(a).length,
   "Array, boolean": (a, unique) => {
     const b = toFlatArray(a);
@@ -5057,7 +5092,7 @@ var setSize = mathTyped8("setSize", {
     return count2;
   }
 });
-var setCartesian = mathTyped8("setCartesian", {
+var setCartesian = mathTyped9("setCartesian", {
   "Array, Array": (a1, a2) => {
     const b1 = toFlatArray(a1).sort(naturalCompare);
     const b2 = toFlatArray(a2).sort(naturalCompare);
@@ -5109,7 +5144,7 @@ var typedSet = {
 };
 
 // src/typed/special.ts
-import { mathTyped as mathTyped9, Complex as Complex5 } from "@danielsimonjr/mathts-core";
+import { mathTyped as mathTyped10, Complex as Complex5 } from "@danielsimonjr/mathts-core";
 
 // src/utils/number.ts
 import {
@@ -6864,7 +6899,7 @@ function fresnelSScalar(x) {
   }
   return sign3 * sum3;
 }
-var erfc = mathTyped9("erfc", {
+var erfc = mathTyped10("erfc", {
   number: erfcScalar,
   Float64Array: (x) => {
     const wasm = elementwiseUnaryDispatch("erfc", x);
@@ -6872,7 +6907,7 @@ var erfc = mathTyped9("erfc", {
     return mapArray(x, erfcScalar, () => kernelSource([_erfSeries, _erfcCF, erfcScalar], "(x) => erfcScalar(x)"));
   }
 });
-var erfi = mathTyped9("erfi", {
+var erfi = mathTyped10("erfi", {
   number: erfiScalar,
   Float64Array: (x) => mapArray(x, erfiScalar, () => kernelSource([erfiScalar], "(x) => erfiScalar(x)"))
 });
@@ -6892,7 +6927,7 @@ function lgammaComplex(z) {
   const term = zc.add(new Complex5(0.5, 0)).mul(t.log());
   return halfLog2Pi.add(term).sub(t).add(x.log());
 }
-var lgamma = mathTyped9("lgamma", {
+var lgamma = mathTyped10("lgamma", {
   number: _lgamma,
   Complex: lgammaComplex,
   Float64Array: (x) => {
@@ -6902,7 +6937,7 @@ var lgamma = mathTyped9("lgamma", {
     return mapArray(x, _lgamma, () => kernelSource([_lgamma], "(x) => _lgamma(x)"));
   }
 });
-var beta = mathTyped9("beta", {
+var beta = mathTyped10("beta", {
   "number, number": betaScalar,
   "Float64Array, number": (a, b) => mapArray(
     a,
@@ -6910,7 +6945,7 @@ var beta = mathTyped9("beta", {
     () => kernelSource([_lgamma, betaScalar], `(a) => betaScalar(a, ${b})`)
   )
 });
-var gammainc = mathTyped9("gammainc", {
+var gammainc = mathTyped10("gammainc", {
   "number, number": gammaincScalar,
   "number, Float64Array": (a, x) => mapArray(
     x,
@@ -6918,7 +6953,7 @@ var gammainc = mathTyped9("gammainc", {
     () => kernelSource([_lgamma, gammaincScalar], `(x) => gammaincScalar(${a}, x)`)
   )
 });
-var gammaincp = mathTyped9("gammaincp", {
+var gammaincp = mathTyped10("gammaincp", {
   "number, number": gammaincpScalar,
   "number, Float64Array": (a, x) => mapArray(
     x,
@@ -6926,7 +6961,7 @@ var gammaincp = mathTyped9("gammaincp", {
     () => kernelSource([_lgamma, gammaincScalar, gammaincpScalar], `(x) => gammaincpScalar(${a}, x)`)
   )
 });
-var betainc = mathTyped9("betainc", {
+var betainc = mathTyped10("betainc", {
   "number, number, number": betaincScalar,
   "number, number, Float64Array": (a, b, x) => mapArray(
     x,
@@ -6934,11 +6969,11 @@ var betainc = mathTyped9("betainc", {
     () => kernelSource([_lgamma, betaincScalar], `(x) => betaincScalar(${a}, ${b}, x)`)
   )
 });
-var digamma = mathTyped9("digamma", {
+var digamma = mathTyped10("digamma", {
   number: digammaScalar,
   Float64Array: (x) => mapArray(x, digammaScalar, () => kernelSource([digammaScalar], "(x) => digammaScalar(x)"))
 });
-var besselJ0 = mathTyped9("besselJ0", {
+var besselJ0 = mathTyped10("besselJ0", {
   number: besselJ0Scalar,
   Float64Array: (x) => {
     if (x.length >= WASM_SPECIAL_THRESHOLD) {
@@ -6951,7 +6986,7 @@ var besselJ0 = mathTyped9("besselJ0", {
     );
   }
 });
-var besselJ1 = mathTyped9("besselJ1", {
+var besselJ1 = mathTyped10("besselJ1", {
   number: besselJ1Scalar,
   Float64Array: (x) => {
     if (x.length >= WASM_SPECIAL_THRESHOLD) {
@@ -6964,7 +6999,7 @@ var besselJ1 = mathTyped9("besselJ1", {
     );
   }
 });
-var besselY0 = mathTyped9("besselY0", {
+var besselY0 = mathTyped10("besselY0", {
   number: besselY0Scalar,
   Float64Array: (x) => {
     if (x.length >= WASM_SPECIAL_THRESHOLD) {
@@ -6977,7 +7012,7 @@ var besselY0 = mathTyped9("besselY0", {
     );
   }
 });
-var besselY1 = mathTyped9("besselY1", {
+var besselY1 = mathTyped10("besselY1", {
   number: besselY1Scalar,
   Float64Array: (x) => {
     if (x.length >= WASM_SPECIAL_THRESHOLD) {
@@ -6990,7 +7025,7 @@ var besselY1 = mathTyped9("besselY1", {
     );
   }
 });
-var besselJ = mathTyped9("besselJ", {
+var besselJ = mathTyped10("besselJ", {
   "number, number": besselJScalar,
   "number, Float64Array": (n, x) => {
     if (x.length >= WASM_SPECIAL_THRESHOLD) {
@@ -7006,7 +7041,7 @@ var besselJ = mathTyped9("besselJ", {
     );
   }
 });
-var besselY = mathTyped9("besselY", {
+var besselY = mathTyped10("besselY", {
   "number, number": besselYScalar,
   "number, Float64Array": (n, x) => {
     if (x.length >= WASM_SPECIAL_THRESHOLD) {
@@ -7022,7 +7057,7 @@ var besselY = mathTyped9("besselY", {
     );
   }
 });
-var besselI = mathTyped9("besselI", {
+var besselI = mathTyped10("besselI", {
   "number, number": besselIScalar,
   "number, Float64Array": (n, x) => mapArray(
     x,
@@ -7030,7 +7065,7 @@ var besselI = mathTyped9("besselI", {
     () => kernelSource([factorial, _lgamma, besselIScalar], `(x) => besselIScalar(${n}, x)`)
   )
 });
-var besselK = mathTyped9("besselK", {
+var besselK = mathTyped10("besselK", {
   "number, number": besselKScalar,
   "number, Float64Array": (n, x) => mapArray(
     x,
@@ -7041,7 +7076,7 @@ var besselK = mathTyped9("besselK", {
     )
   )
 });
-var ellipticK = mathTyped9("ellipticK", {
+var ellipticK = mathTyped10("ellipticK", {
   number: ellipticKScalar,
   Float64Array: (m) => {
     if (m.length >= WASM_SPECIAL_THRESHOLD) {
@@ -7054,7 +7089,7 @@ var ellipticK = mathTyped9("ellipticK", {
     );
   }
 });
-var ellipticE = mathTyped9("ellipticE", {
+var ellipticE = mathTyped10("ellipticE", {
   "number, number": ellipticEScalar,
   number: ellipticECompleteScalar,
   Float64Array: (m) => {
@@ -7068,7 +7103,7 @@ var ellipticE = mathTyped9("ellipticE", {
     );
   }
 });
-var chebyshevT = mathTyped9("chebyshevT", {
+var chebyshevT = mathTyped10("chebyshevT", {
   "number, number": chebyshevTScalar,
   "number, Float64Array": (n, x) => mapArray(
     x,
@@ -7076,7 +7111,7 @@ var chebyshevT = mathTyped9("chebyshevT", {
     () => kernelSource([chebyshevTScalar], `(x) => chebyshevTScalar(${n}, x)`)
   )
 });
-var hermiteH = mathTyped9("hermiteH", {
+var hermiteH = mathTyped10("hermiteH", {
   "number, number": hermiteHScalar,
   "number, Float64Array": (n, x) => mapArray(
     x,
@@ -7084,7 +7119,7 @@ var hermiteH = mathTyped9("hermiteH", {
     () => kernelSource([hermiteHScalar], `(x) => hermiteHScalar(${n}, x)`)
   )
 });
-var laguerreL = mathTyped9("laguerreL", {
+var laguerreL = mathTyped10("laguerreL", {
   "number, number": laguerreLScalar,
   "number, Float64Array": (n, x) => mapArray(
     x,
@@ -7092,7 +7127,7 @@ var laguerreL = mathTyped9("laguerreL", {
     () => kernelSource([laguerreLScalar], `(x) => laguerreLScalar(${n}, x)`)
   )
 });
-var legendreP = mathTyped9("legendreP", {
+var legendreP = mathTyped10("legendreP", {
   "number, number": legendrePScalar,
   "number, Float64Array": (n, x) => mapArray(
     x,
@@ -7100,11 +7135,11 @@ var legendreP = mathTyped9("legendreP", {
     () => kernelSource([legendrePScalar], `(x) => legendrePScalar(${n}, x)`)
   )
 });
-var lambertW = mathTyped9("lambertW", {
+var lambertW = mathTyped10("lambertW", {
   number: lambertWScalar,
   Float64Array: (x) => mapArray(x, lambertWScalar, () => kernelSource([lambertWScalar], "(x) => lambertWScalar(x)"))
 });
-var cosIntegral = mathTyped9("cosIntegral", {
+var cosIntegral = mathTyped10("cosIntegral", {
   number: cosIntegralScalar,
   Float64Array: (x) => mapArray(
     x,
@@ -7112,7 +7147,7 @@ var cosIntegral = mathTyped9("cosIntegral", {
     () => kernelSource([cosIntegralScalar], "(x) => cosIntegralScalar(x)")
   )
 });
-var sinIntegral = mathTyped9("sinIntegral", {
+var sinIntegral = mathTyped10("sinIntegral", {
   number: sinIntegralScalar,
   Float64Array: (x) => mapArray(
     x,
@@ -7120,7 +7155,7 @@ var sinIntegral = mathTyped9("sinIntegral", {
     () => kernelSource([sinIntegralScalar], "(x) => sinIntegralScalar(x)")
   )
 });
-var logIntegral = mathTyped9("logIntegral", {
+var logIntegral = mathTyped10("logIntegral", {
   number: logIntegralScalar,
   Float64Array: (x) => mapArray(
     x,
@@ -7128,7 +7163,7 @@ var logIntegral = mathTyped9("logIntegral", {
     () => kernelSource([expIntegralEiScalar, logIntegralScalar], "(x) => logIntegralScalar(x)")
   )
 });
-var expIntegralEi = mathTyped9("expIntegralEi", {
+var expIntegralEi = mathTyped10("expIntegralEi", {
   number: expIntegralEiScalar,
   Float64Array: (x) => mapArray(
     x,
@@ -7136,15 +7171,15 @@ var expIntegralEi = mathTyped9("expIntegralEi", {
     () => kernelSource([expIntegralEiScalar], "(x) => expIntegralEiScalar(x)")
   )
 });
-var fresnelC = mathTyped9("fresnelC", {
+var fresnelC = mathTyped10("fresnelC", {
   number: fresnelCScalar,
   Float64Array: (x) => mapArray(x, fresnelCScalar, () => kernelSource([fresnelCScalar], "(x) => fresnelCScalar(x)"))
 });
-var fresnelS = mathTyped9("fresnelS", {
+var fresnelS = mathTyped10("fresnelS", {
   number: fresnelSScalar,
   Float64Array: (x) => mapArray(x, fresnelSScalar, () => kernelSource([fresnelSScalar], "(x) => fresnelSScalar(x)"))
 });
-var airyAi = mathTyped9("airyAi", {
+var airyAi = mathTyped10("airyAi", {
   number: airyAiScalar,
   Float64Array: (x) => {
     if (x.length >= WASM_SPECIAL_THRESHOLD) {
@@ -7153,7 +7188,7 @@ var airyAi = mathTyped9("airyAi", {
     return mapArray(x, airyAiScalar, () => kernelSource([airyUCoeffs, airyAsymPQ, airyAiScalar], "(x) => airyAiScalar(x)"));
   }
 });
-var airyBi = mathTyped9("airyBi", {
+var airyBi = mathTyped10("airyBi", {
   number: airyBiScalar,
   Float64Array: (x) => {
     if (x.length >= WASM_SPECIAL_THRESHOLD) {
@@ -7162,49 +7197,49 @@ var airyBi = mathTyped9("airyBi", {
     return mapArray(x, airyBiScalar, () => kernelSource([airyUCoeffs, airyAsymPQ, airyBiScalar], "(x) => airyBiScalar(x)"));
   }
 });
-var carlsonRC = mathTyped9("carlsonRC", {
+var carlsonRC = mathTyped10("carlsonRC", {
   "number, number": carlsonRCScalar,
   "Float64Array, Float64Array": (xs, ys) => {
     if (xs.length >= WASM_SPECIAL_THRESHOLD) return carlsonRCDispatch(xs, ys);
     return carlsonRCJS(xs, ys);
   }
 });
-var carlsonRF = mathTyped9("carlsonRF", {
+var carlsonRF = mathTyped10("carlsonRF", {
   "number, number, number": carlsonRFScalar,
   "Float64Array, Float64Array, Float64Array": (xs, ys, zs) => {
     if (xs.length >= WASM_SPECIAL_THRESHOLD) return carlsonRFDispatch(xs, ys, zs);
     return carlsonRFJS(xs, ys, zs);
   }
 });
-var carlsonRD = mathTyped9("carlsonRD", {
+var carlsonRD = mathTyped10("carlsonRD", {
   "number, number, number": carlsonRDScalar,
   "Float64Array, Float64Array, Float64Array": (xs, ys, zs) => {
     if (xs.length >= WASM_SPECIAL_THRESHOLD) return carlsonRDDispatch(xs, ys, zs);
     return carlsonRDJS(xs, ys, zs);
   }
 });
-var carlsonRJ = mathTyped9("carlsonRJ", {
+var carlsonRJ = mathTyped10("carlsonRJ", {
   "number, number, number, number": carlsonRJScalar,
   "Float64Array, Float64Array, Float64Array, Float64Array": (xs, ys, zs, ps) => {
     if (xs.length >= WASM_SPECIAL_THRESHOLD) return carlsonRJDispatch(xs, ys, zs, ps);
     return carlsonRJJS(xs, ys, zs, ps);
   }
 });
-var ellipticF = mathTyped9("ellipticF", {
+var ellipticF = mathTyped10("ellipticF", {
   "number, number": ellipticFIncompleteScalar,
   "Float64Array, Float64Array": (phis, ms) => {
     if (phis.length >= WASM_SPECIAL_THRESHOLD) return ellipticFIncompleteDispatch(phis, ms);
     return ellipticFIncompleteJS(phis, ms);
   }
 });
-var ellipticEIncomplete = mathTyped9("ellipticEIncomplete", {
+var ellipticEIncomplete = mathTyped10("ellipticEIncomplete", {
   "number, number": ellipticEIncompleteScalar,
   "Float64Array, Float64Array": (phis, ms) => {
     if (phis.length >= WASM_SPECIAL_THRESHOLD) return ellipticEIncompleteDispatch(phis, ms);
     return ellipticEIncompleteJS(phis, ms);
   }
 });
-var ellipticPi = mathTyped9("ellipticPi", {
+var ellipticPi = mathTyped10("ellipticPi", {
   "number, number, number": ellipticPiIncompleteScalar,
   "Float64Array, Float64Array, Float64Array": (ns, phis, ms) => {
     if (ns.length >= WASM_SPECIAL_THRESHOLD) return ellipticPiIncompleteDispatch(ns, phis, ms);
@@ -7285,7 +7320,7 @@ function fuseUnaryChain(ops, xs) {
 }
 
 // src/typed/distributions.ts
-import { mathTyped as mathTyped10 } from "@danielsimonjr/mathts-core";
+import { mathTyped as mathTyped11 } from "@danielsimonjr/mathts-core";
 var LN2 = Math.log(2);
 function _lgammaD(x) {
   if (x <= 0 && x === Math.floor(x)) return Infinity;
@@ -7401,7 +7436,7 @@ function bernoulliPMFScalar(k, p) {
   if (k === 0) return 1 - p;
   return NaN;
 }
-var normalPDF = mathTyped10("normalPDF", {
+var normalPDF = mathTyped11("normalPDF", {
   number: (x) => normalPDFScalar(x, 0, 1),
   "number, number, number": (x, mu, sigma) => normalPDFScalar(x, mu, sigma),
   Float64Array: (x) => mapArray(
@@ -7415,7 +7450,7 @@ var normalPDF = mathTyped10("normalPDF", {
     () => kernelSource([normalPDFScalar], `(x) => normalPDFScalar(x, ${mu}, ${sigma})`)
   )
 });
-var normalCDF = mathTyped10("normalCDF", {
+var normalCDF = mathTyped11("normalCDF", {
   number: (x) => normalCDFScalar(x, 0, 1),
   "number, number, number": (x, mu, sigma) => normalCDFScalar(x, mu, sigma),
   Float64Array: (x) => mapArray(
@@ -7429,7 +7464,7 @@ var normalCDF = mathTyped10("normalCDF", {
     () => kernelSource([_erfc, normalCDFScalar], `(x) => normalCDFScalar(x, ${mu}, ${sigma})`)
   )
 });
-var exponentialPDF = mathTyped10("exponentialPDF", {
+var exponentialPDF = mathTyped11("exponentialPDF", {
   "number, number": (x, lambda) => exponentialPDFScalar(x, lambda),
   "Float64Array, number": (x, lambda) => mapArray(
     x,
@@ -7437,7 +7472,7 @@ var exponentialPDF = mathTyped10("exponentialPDF", {
     () => kernelSource([exponentialPDFScalar], `(x) => exponentialPDFScalar(x, ${lambda})`)
   )
 });
-var exponentialCDF = mathTyped10("exponentialCDF", {
+var exponentialCDF = mathTyped11("exponentialCDF", {
   "number, number": (x, lambda) => exponentialCDFScalar(x, lambda),
   "Float64Array, number": (x, lambda) => mapArray(
     x,
@@ -7445,7 +7480,7 @@ var exponentialCDF = mathTyped10("exponentialCDF", {
     () => kernelSource([exponentialCDFScalar], `(x) => exponentialCDFScalar(x, ${lambda})`)
   )
 });
-var poissonPMF = mathTyped10("poissonPMF", {
+var poissonPMF = mathTyped11("poissonPMF", {
   "number, number": (k, lambda) => poissonPMFScalar(k, lambda),
   "Float64Array, number": (k, lambda) => mapArray(
     k,
@@ -7453,7 +7488,7 @@ var poissonPMF = mathTyped10("poissonPMF", {
     () => kernelSource([_logFactorial, poissonPMFScalar], `(k) => poissonPMFScalar(k, ${lambda})`)
   )
 });
-var binomialPMF = mathTyped10("binomialPMF", {
+var binomialPMF = mathTyped11("binomialPMF", {
   "number, number, number": (k, n, p) => binomialPMFScalar(k, n, p),
   "Float64Array, number, number": (k, n, p) => mapArray(
     k,
@@ -7464,7 +7499,7 @@ var binomialPMF = mathTyped10("binomialPMF", {
     )
   )
 });
-var geometricPMF = mathTyped10("geometricPMF", {
+var geometricPMF = mathTyped11("geometricPMF", {
   "number, number": (k, p) => geometricPMFScalar(k, p),
   "Float64Array, number": (k, p) => mapArray(
     k,
@@ -7472,7 +7507,7 @@ var geometricPMF = mathTyped10("geometricPMF", {
     () => kernelSource([geometricPMFScalar], `(k) => geometricPMFScalar(k, ${p})`)
   )
 });
-var bernoulliPMF = mathTyped10("bernoulliPMF", {
+var bernoulliPMF = mathTyped11("bernoulliPMF", {
   "number, number": (k, p) => bernoulliPMFScalar(k, p),
   "Float64Array, number": (k, p) => mapArray(
     k,
@@ -7480,7 +7515,7 @@ var bernoulliPMF = mathTyped10("bernoulliPMF", {
     () => kernelSource([bernoulliPMFScalar], `(k) => bernoulliPMFScalar(k, ${p})`)
   )
 });
-var entropy = mathTyped10("entropy", {
+var entropy = mathTyped11("entropy", {
   Array: (probs) => {
     let h = 0;
     for (const p of probs) {
@@ -7492,7 +7527,7 @@ var entropy = mathTyped10("entropy", {
     return h;
   }
 });
-var jsDivergence = mathTyped10("jsDivergence", {
+var jsDivergence = mathTyped11("jsDivergence", {
   "Array, Array": (p, q) => {
     if (p.length !== q.length) return NaN;
     let jsd = 0;
@@ -7513,7 +7548,7 @@ function betaPDFScalar(x, alpha, beta_) {
   const logB = _lgammaD(alpha) + _lgammaD(beta_) - _lgammaD(alpha + beta_);
   return Math.exp((alpha - 1) * Math.log(x) + (beta_ - 1) * Math.log(1 - x) - logB);
 }
-var betaPDF = mathTyped10("betaPDF", {
+var betaPDF = mathTyped11("betaPDF", {
   "number, number, number": (x, alpha, beta_) => betaPDFScalar(x, alpha, beta_),
   "Float64Array, number, number": (x, alpha, beta_) => {
     const out = new Float64Array(x.length);
@@ -7545,7 +7580,7 @@ function gammaPDFScalar(x, shape, scale) {
     (shape - 1) * Math.log(x) - x / scale - _lgammaD(shape) - shape * Math.log(scale)
   );
 }
-var gammaPDF = mathTyped10("gammaPDF", {
+var gammaPDF = mathTyped11("gammaPDF", {
   "number, number, number": (x, shape, scale) => gammaPDFScalar(x, shape, scale),
   "Float64Array, number, number": (x, shape, scale) => {
     const out = new Float64Array(x.length);
@@ -7570,7 +7605,7 @@ function studentTPDFScalar(x, df) {
   const logNorm = lgHalfNu1 - lgHalfNu - 0.5 * Math.log(df * Math.PI);
   return Math.exp(logNorm - (df + 1) / 2 * Math.log(1 + x * x / df));
 }
-var studentTPDF = mathTyped10("studentTPDF", {
+var studentTPDF = mathTyped11("studentTPDF", {
   "number, number": (x, df) => studentTPDFScalar(x, df),
   "Float64Array, number": (x, df) => {
     const out = new Float64Array(x.length);
@@ -7601,7 +7636,7 @@ function noncentralChi2PDFScalar(x, df, ncp) {
   }
   return sum3;
 }
-var noncentralChi2PDF = mathTyped10("noncentralChi2PDF", {
+var noncentralChi2PDF = mathTyped11("noncentralChi2PDF", {
   "number, number, number": (x, df, ncp) => noncentralChi2PDFScalar(x, df, ncp),
   "Float64Array, number, number": (x, df, ncp) => {
     const out = new Float64Array(x.length);
@@ -11750,8 +11785,8 @@ function solvePDE(pde, domain, bc) {
 }
 
 // src/typed/combinatorics.ts
-import { mathTyped as mathTyped11 } from "@danielsimonjr/mathts-core";
-var fibonacci = mathTyped11("fibonacci", {
+import { mathTyped as mathTyped12 } from "@danielsimonjr/mathts-core";
+var fibonacci = mathTyped12("fibonacci", {
   number: (n) => {
     if (n < 0 || !Number.isInteger(n)) {
       throw new Error("fibonacci requires a non-negative integer");
@@ -11780,7 +11815,7 @@ var fibonacci = mathTyped11("fibonacci", {
     return a;
   }
 });
-var lucas = mathTyped11("lucas", {
+var lucas = mathTyped12("lucas", {
   number: (n) => {
     if (n < 0 || !Number.isInteger(n)) {
       throw new Error("lucas requires a non-negative integer");
@@ -11809,7 +11844,7 @@ var lucas = mathTyped11("lucas", {
     return 2 * b - a;
   }
 });
-var doubleFactorial = mathTyped11("doubleFactorial", {
+var doubleFactorial = mathTyped12("doubleFactorial", {
   number: (n) => {
     if (!Number.isInteger(n) || n < -1) {
       throw new Error("doubleFactorial requires an integer >= -1");
@@ -11822,7 +11857,7 @@ var doubleFactorial = mathTyped11("doubleFactorial", {
     return result;
   }
 });
-var risingFactorial = mathTyped11("risingFactorial", {
+var risingFactorial = mathTyped12("risingFactorial", {
   "number, number": (x, n) => {
     if (!Number.isInteger(n) || n < 0) {
       throw new Error("risingFactorial requires a non-negative integer for n");
@@ -11835,7 +11870,7 @@ var risingFactorial = mathTyped11("risingFactorial", {
     return result;
   }
 });
-var fallingFactorial = mathTyped11("fallingFactorial", {
+var fallingFactorial = mathTyped12("fallingFactorial", {
   "number, number": (x, n) => {
     if (!Number.isInteger(n) || n < 0) {
       throw new Error("fallingFactorial requires a non-negative integer for n");
@@ -11848,7 +11883,7 @@ var fallingFactorial = mathTyped11("fallingFactorial", {
     return result;
   }
 });
-var subfactorial = mathTyped11("subfactorial", {
+var subfactorial = mathTyped12("subfactorial", {
   number: (n) => {
     if (n < 0 || !Number.isInteger(n)) {
       throw new Error("subfactorial requires a non-negative integer");
@@ -11865,7 +11900,7 @@ var subfactorial = mathTyped11("subfactorial", {
     return prev1;
   }
 });
-var prime = mathTyped11("prime", {
+var prime = mathTyped12("prime", {
   number: (n) => {
     if (n < 1 || !Number.isInteger(n)) {
       throw new Error("prime requires a positive integer");
@@ -11896,7 +11931,7 @@ function _sieveOfEratosthenes(limit2) {
   }
   return primes;
 }
-var nextPrime = mathTyped11("nextPrime", {
+var nextPrime = mathTyped12("nextPrime", {
   number: (n) => {
     if (n < 0) throw new Error("nextPrime requires a non-negative number");
     let candidate = Math.floor(n) + 1;
@@ -11917,7 +11952,7 @@ function _isPrime(n) {
   }
   return true;
 }
-var primePi = mathTyped11("primePi", {
+var primePi = mathTyped12("primePi", {
   number: (n) => {
     if (n < 2) return 0;
     const limit2 = Math.floor(n);
@@ -11925,7 +11960,7 @@ var primePi = mathTyped11("primePi", {
     return sieve.length;
   }
 });
-var primeFactors = mathTyped11("primeFactors", {
+var primeFactors = mathTyped12("primeFactors", {
   number: (n) => {
     if (n < 2 || !Number.isInteger(n)) {
       throw new Error("primeFactors requires an integer >= 2");
@@ -11942,7 +11977,7 @@ var primeFactors = mathTyped11("primeFactors", {
     return factors;
   }
 });
-var divisors = mathTyped11("divisors", {
+var divisors = mathTyped12("divisors", {
   number: (n) => {
     if (n < 1 || !Number.isInteger(n)) {
       throw new Error("divisors requires a positive integer");
@@ -11958,7 +11993,7 @@ var divisors = mathTyped11("divisors", {
     return small.concat(large.reverse());
   }
 });
-var eulerPhi = mathTyped11("eulerPhi", {
+var eulerPhi = mathTyped12("eulerPhi", {
   number: (n) => {
     if (n < 1 || !Number.isInteger(n)) {
       throw new Error("eulerPhi requires a positive integer");
@@ -11975,7 +12010,7 @@ var eulerPhi = mathTyped11("eulerPhi", {
     return result;
   }
 });
-var divisorSigma = mathTyped11("divisorSigma", {
+var divisorSigma = mathTyped12("divisorSigma", {
   number: (n) => {
     return _divisorSigmaImpl(n, 1);
   },
@@ -11999,7 +12034,7 @@ function _divisorSigmaImpl(n, k) {
   }
   return sum3;
 }
-var carmichaelLambda = mathTyped11("carmichaelLambda", {
+var carmichaelLambda = mathTyped12("carmichaelLambda", {
   number: (n) => {
     if (n < 1 || !Number.isInteger(n)) {
       throw new Error("carmichaelLambda requires a positive integer");
@@ -12043,7 +12078,7 @@ function _gcd(a, b) {
 function _lcm(a, b) {
   return a / _gcd(a, b) * b;
 }
-var moebiusMu = mathTyped11("moebiusMu", {
+var moebiusMu = mathTyped12("moebiusMu", {
   number: (n) => {
     if (n < 1 || !Number.isInteger(n)) {
       throw new Error("moebiusMu requires a positive integer");
@@ -12062,7 +12097,7 @@ var moebiusMu = mathTyped11("moebiusMu", {
     return numFactors % 2 === 0 ? 1 : -1;
   }
 });
-var jacobiSymbol = mathTyped11("jacobiSymbol", {
+var jacobiSymbol = mathTyped12("jacobiSymbol", {
   "number, number": (a, n) => {
     if (n < 1 || !Number.isInteger(n) || n % 2 === 0) {
       throw new Error("jacobiSymbol: n must be a positive odd integer");
@@ -12087,7 +12122,7 @@ var jacobiSymbol = mathTyped11("jacobiSymbol", {
     return nn === 1 ? result : 0;
   }
 });
-var chineseRemainder = mathTyped11("chineseRemainder", {
+var chineseRemainder = mathTyped12("chineseRemainder", {
   "Array, Array": (remainders, moduli) => {
     if (remainders.length !== moduli.length) {
       throw new Error("chineseRemainder: remainders and moduli must have same length");
@@ -12126,7 +12161,7 @@ function _modInverse(a, m) {
   if (old_r !== 1) return -1;
   return (old_s % m + m) % m;
 }
-var lucasL = mathTyped11("lucasL", {
+var lucasL = mathTyped12("lucasL", {
   number: (n) => {
     if (n < 0 || !Number.isInteger(n)) {
       throw new Error("lucasL requires a non-negative integer");
@@ -12141,7 +12176,7 @@ var lucasL = mathTyped11("lucasL", {
     return b;
   }
 });
-var partitions = mathTyped11("partitions", {
+var partitions = mathTyped12("partitions", {
   number: (n) => {
     if (n < 0 || !Number.isInteger(n)) {
       throw new Error("partitions requires a non-negative integer");
@@ -12156,7 +12191,7 @@ var partitions = mathTyped11("partitions", {
     return dp[n];
   }
 });
-var harmonicNumber = mathTyped11("harmonicNumber", {
+var harmonicNumber = mathTyped12("harmonicNumber", {
   number: (n) => {
     if (n < 1 || !Number.isInteger(n)) {
       throw new Error("harmonicNumber requires a positive integer");
@@ -12168,7 +12203,7 @@ var harmonicNumber = mathTyped11("harmonicNumber", {
     return sum3;
   }
 });
-var integerDigits = mathTyped11("integerDigits", {
+var integerDigits = mathTyped12("integerDigits", {
   number: (n) => {
     return _integerDigitsImpl(n, 10);
   },
@@ -14215,7 +14250,7 @@ import {
   matrixSqrtm as matrixSqrtmPrimitive,
   DenseMatrix as DenseMatrix3
 } from "@danielsimonjr/mathts-matrix";
-import { mathTyped as mathTyped12 } from "@danielsimonjr/mathts-core";
+import { mathTyped as mathTyped13 } from "@danielsimonjr/mathts-core";
 import { computePool as computePool11 } from "@danielsimonjr/mathts-parallel";
 function cloneMatrix(A) {
   return A.map((r) => [...r]);
@@ -14845,7 +14880,7 @@ async function solveApprox(A, b, n) {
   }
   return x;
 }
-var pinv = mathTyped12("pinv", {
+var pinv = mathTyped13("pinv", {
   DenseMatrix: (A) => matrixPinv(A),
   "DenseMatrix, Object": (A, opts) => matrixPinv(A, opts),
   // Array in → Array out (mathjs convention). Previously `pinv([[...]])` threw
@@ -14853,22 +14888,22 @@ var pinv = mathTyped12("pinv", {
   Array: (A) => matrixPinv(DenseMatrix3.fromArray(A)).toArray(),
   "Array, Object": (A, opts) => matrixPinv(DenseMatrix3.fromArray(A), opts).toArray()
 });
-var cond = mathTyped12("cond", {
+var cond = mathTyped13("cond", {
   Array: (A) => matrixCond(A)
 });
-var norm2 = mathTyped12("norm2", {
+var norm2 = mathTyped13("norm2", {
   Array: (A) => matrixNorm2(A)
 });
-var normFro = mathTyped12("normFro", {
+var normFro = mathTyped13("normFro", {
   Array: (A) => matrixNormFro(A)
 });
-var lowRankApprox = mathTyped12("lowRankApprox", {
+var lowRankApprox = mathTyped13("lowRankApprox", {
   "Array, number": (A, k) => matrixLowRankApprox(A, k)
 });
-var singularValues2 = mathTyped12("singularValues", {
+var singularValues2 = mathTyped13("singularValues", {
   Array: (A) => matrixSingularValues(A)
 });
-var matrixExpm = mathTyped12("matrixExpm", {
+var matrixExpm = mathTyped13("matrixExpm", {
   DenseMatrix: (A) => matrixExpmPrimitive(A),
   "DenseMatrix, Object": (A, _opts) => matrixExpmPrimitive(A),
   Array: (A) => {
@@ -14876,7 +14911,7 @@ var matrixExpm = mathTyped12("matrixExpm", {
     return matrixExpmPrimitive(D).toArray();
   }
 });
-var matrixLogm = mathTyped12("matrixLogm", {
+var matrixLogm = mathTyped13("matrixLogm", {
   DenseMatrix: (A) => matrixLogmPrimitive(A),
   "DenseMatrix, Object": (A, opts) => matrixLogmPrimitive(A, opts),
   Array: (A) => {
@@ -14884,7 +14919,7 @@ var matrixLogm = mathTyped12("matrixLogm", {
     return matrixLogmPrimitive(D).toArray();
   }
 });
-var matrixSqrtm = mathTyped12("matrixSqrtm", {
+var matrixSqrtm = mathTyped13("matrixSqrtm", {
   DenseMatrix: (A) => matrixSqrtmPrimitive(A),
   "DenseMatrix, Object": (A, opts) => matrixSqrtmPrimitive(A, opts),
   Array: (A) => {
@@ -14919,7 +14954,7 @@ async function gpuScale(a, scalar) {
 
 // src/typed/relational.ts
 import naturalSort from "javascript-natural-sort";
-import { mathTyped as mathTyped13, Complex as Complex6, BigNumber as BigNumber6, Fraction as Fraction2 } from "@danielsimonjr/mathts-core";
+import { mathTyped as mathTyped14, Complex as Complex6, BigNumber as BigNumber6, Fraction as Fraction2 } from "@danielsimonjr/mathts-core";
 var REL_TOL = 1e-12;
 var ABS_TOL = 1e-15;
 function numbersNearlyEqual(x, y) {
@@ -14967,7 +15002,7 @@ function compareObjects(cmp, x, y) {
   }
   return 0;
 }
-var equalScalar = mathTyped13("equalScalar", {
+var equalScalar = mathTyped14("equalScalar", {
   "boolean, boolean": (x, y) => x === y,
   "number, number": (x, y) => numbersNearlyEqual(x, y),
   "bigint, bigint": (x, y) => x === y,
@@ -14983,7 +15018,7 @@ var equalScalar = mathTyped13("equalScalar", {
     return JSON.stringify(x) === JSON.stringify(y);
   }
 });
-var unequal = mathTyped13("unequal", {
+var unequal = mathTyped14("unequal", {
   "any, any": (x, y) => {
     if (x === null) return y !== null;
     if (y === null) return x !== null;
@@ -14992,7 +15027,7 @@ var unequal = mathTyped13("unequal", {
     return !equalScalar(x, y);
   }
 });
-var deepEqual = mathTyped13("deepEqual", {
+var deepEqual = mathTyped14("deepEqual", {
   "any, any": (x, y) => {
     return _deepEqual(
       x.valueOf ? x.valueOf() : x,
@@ -15012,7 +15047,7 @@ function _deepEqual(x, y) {
   if (Array.isArray(y)) return false;
   return equalScalar(x, y);
 }
-var compareText = mathTyped13("compareText", {
+var compareText = mathTyped14("compareText", {
   "string, string": (x, y) => {
     return x < y ? -1 : x > y ? 1 : 0;
   },
@@ -15022,12 +15057,12 @@ var compareText = mathTyped13("compareText", {
     return sx < sy ? -1 : sx > sy ? 1 : 0;
   }
 });
-var equalText = mathTyped13("equalText", {
+var equalText = mathTyped14("equalText", {
   "any, any": (x, y) => {
     return compareText(x, y) === 0;
   }
 });
-var compareNatural = mathTyped13("compareNatural", {
+var compareNatural = mathTyped14("compareNatural", {
   "any, any": (x, y) => _compareNatural(x, y)
 });
 function _compareNatural(x, y) {
@@ -15102,7 +15137,7 @@ function _numericCompare(x, y) {
   if (Number.isNaN(ny)) return -1;
   return nx > ny ? 1 : nx < ny ? -1 : 0;
 }
-var compareUnits = mathTyped13("compareUnits", {
+var compareUnits = mathTyped14("compareUnits", {
   "any, any": (x, y) => {
     const ux = x;
     const uy = y;
@@ -15126,7 +15161,7 @@ var typedRelational = {
 };
 
 // src/typed/string.ts
-import { mathTyped as mathTyped14, BigNumber as BigNumber7 } from "@danielsimonjr/mathts-core";
+import { mathTyped as mathTyped15, BigNumber as BigNumber7 } from "@danielsimonjr/mathts-core";
 
 // src/utils/is.ts
 function isNumber(x) {
@@ -15542,25 +15577,25 @@ function interpolate2(template, values, options) {
     return original;
   });
 }
-var bin = mathTyped14("bin", {
+var bin = mathTyped15("bin", {
   number: (x) => numberToBase(x, 2),
   "number, number": (x, wordSize) => numberToBase(x, 2, wordSize),
   BigNumber: (x) => bigNumberToBase(x, 2),
   "BigNumber, number": (x, wordSize) => bigNumberToBase(x, 2, wordSize)
 });
-var hex = mathTyped14("hex", {
+var hex = mathTyped15("hex", {
   number: (x) => numberToBase(x, 16),
   "number, number": (x, wordSize) => numberToBase(x, 16, wordSize),
   BigNumber: (x) => bigNumberToBase(x, 16),
   "BigNumber, number": (x, wordSize) => bigNumberToBase(x, 16, wordSize)
 });
-var oct = mathTyped14("oct", {
+var oct = mathTyped15("oct", {
   number: (x) => numberToBase(x, 8),
   "number, number": (x, wordSize) => numberToBase(x, 8, wordSize),
   BigNumber: (x) => bigNumberToBase(x, 8),
   "BigNumber, number": (x, wordSize) => bigNumberToBase(x, 8, wordSize)
 });
-var format4 = mathTyped14("format", {
+var format4 = mathTyped15("format", {
   // Single-arg: format with defaults. For BigNumber, utils/string isBigNumber check
   // may fail (core BigNumber uses .isBigNumber property but not on prototype), so
   // we intercept BigNumber explicitly before delegating to formatValue.
@@ -15588,7 +15623,7 @@ var format4 = mathTyped14("format", {
     return format3(x, options);
   }
 });
-var print = mathTyped14("print", {
+var print = mathTyped15("print", {
   "string, Object | Array": (template, values) => interpolate2(template, values),
   "string, Object | Array, number | Object": (template, values, options) => interpolate2(template, values, options)
 });
@@ -15601,7 +15636,7 @@ var typedString = {
 };
 
 // src/typed/probability.ts
-import { mathTyped as mathTyped15 } from "@danielsimonjr/mathts-core";
+import { mathTyped as mathTyped16 } from "@danielsimonjr/mathts-core";
 import seedrandom from "seedrandom";
 function makeRng(seed) {
   return seed === null ? seedrandom() : seedrandom(String(seed));
@@ -15666,13 +15701,13 @@ function combinationsExact(n, k) {
   }
   return Math.round(r);
 }
-var bernoulli = mathTyped15("bernoulli", {
+var bernoulli = mathTyped16("bernoulli", {
   number: (n) => bernoulliNumber(n)
 });
-var combinations = mathTyped15("combinations", {
+var combinations = mathTyped16("combinations", {
   "number, number": (n, k) => combinationsExact(n, k)
 });
-var combinationsWithRep = mathTyped15("combinationsWithRep", {
+var combinationsWithRep = mathTyped16("combinationsWithRep", {
   "number, number": (n, k) => {
     if (!Number.isInteger(n) || n < 1) {
       throw new TypeError(
@@ -15685,7 +15720,7 @@ var combinationsWithRep = mathTyped15("combinationsWithRep", {
     return combinationsExact(n + k - 1, k);
   }
 });
-var multinomial = mathTyped15("multinomial", {
+var multinomial = mathTyped16("multinomial", {
   Array: (a) => {
     let sum3 = 0;
     let denomLog = 0;
@@ -15701,7 +15736,7 @@ var multinomial = mathTyped15("multinomial", {
     return Math.round(Math.exp(numLog - denomLog));
   }
 });
-var permutations = mathTyped15("permutations", {
+var permutations = mathTyped16("permutations", {
   number: (n) => {
     if (!Number.isInteger(n) || n < 0) {
       throw new TypeError("Non-negative integer expected in function permutations");
@@ -15727,7 +15762,7 @@ var _rng = makeRng(null);
 function seedProbabilityRng(seed) {
   _rng = makeRng(seed);
 }
-var random = mathTyped15("random", {
+var random = mathTyped16("random", {
   "": () => _rng(),
   number: (max2) => {
     if (max2 < 0) throw new RangeError("max must be non-negative");
@@ -15738,7 +15773,7 @@ var random = mathTyped15("random", {
     return min2 + _rng() * (max2 - min2);
   }
 });
-var randomInt = mathTyped15("randomInt", {
+var randomInt = mathTyped16("randomInt", {
   "": () => Math.floor(_rng() * 2),
   number: (max2) => {
     if (!Number.isInteger(max2) || max2 < 1) {
@@ -15753,7 +15788,7 @@ var randomInt = mathTyped15("randomInt", {
     return Math.floor(min2 + _rng() * (max2 - min2));
   }
 });
-var pickRandom = mathTyped15("pickRandom", {
+var pickRandom = mathTyped16("pickRandom", {
   Array: (arr7) => _pick(arr7, 1, void 0)[0],
   "Array, number": (arr7, n) => _pick(arr7, n, void 0),
   "Array, Array": (arr7, weights) => _pick(arr7, 1, weights)[0],
@@ -15809,13 +15844,13 @@ var typedProbability = {
 };
 
 // src/typed/unit.ts
-import { mathTyped as mathTyped16 } from "@danielsimonjr/mathts-core";
+import { mathTyped as mathTyped17 } from "@danielsimonjr/mathts-core";
 import { Unit } from "@danielsimonjr/mathts-core";
-var to = mathTyped16("to", {
+var to = mathTyped17("to", {
   "Unit, string": (value, target) => value.to(target),
   "number, string": (value, target) => new Unit(value, target)
 });
-var toBest = mathTyped16("toBest", {
+var toBest = mathTyped17("toBest", {
   Unit: (value) => value.toBest()
 });
 var typedUnit = {
@@ -15875,7 +15910,7 @@ import {
 } from "@danielsimonjr/mathts-core";
 
 // src/factories/scope.ts
-import { mathTyped as mathTyped17, Complex as Complex7, Fraction as Fraction3, BigNumber as BigNumber8 } from "@danielsimonjr/mathts-core";
+import { mathTyped as mathTyped18, Complex as Complex7, Fraction as Fraction3, BigNumber as BigNumber8 } from "@danielsimonjr/mathts-core";
 
 // src/typed/typed-bridge.ts
 import { registerNativeTypes } from "@danielsimonjr/mathts-core";
@@ -15919,7 +15954,7 @@ var DEFAULT_CONFIG = {
 
 // src/factories/scope.ts
 initTypeBridge();
-var typed = mathTyped17;
+var typed = mathTyped18;
 var extraTypes = [
   // mathjs uses lowercase 'function' in signatures; typed-function registers 'Function'
   {
@@ -16013,7 +16048,7 @@ for (const t of extraTypes) {
   }
 }
 var factoryScope = {
-  typed: mathTyped17,
+  typed: mathTyped18,
   config: DEFAULT_CONFIG,
   // Class constructors (uppercase) — factories use `new BigNumber(...)` etc.
   Complex: Complex7,
@@ -39652,6 +39687,33 @@ function numberFactory(name254, value) {
 }
 
 // src/factories/index.ts
+import {
+  createAndTransform,
+  createBitAndTransform,
+  createBitOrTransform,
+  createColumnTransform,
+  createConcatTransform,
+  createCumSumTransform,
+  createDiffTransform,
+  createFilterTransform,
+  createForEachTransform,
+  createIndexTransform,
+  createMapTransform,
+  createMapSlicesTransform,
+  createMaxTransform,
+  createMeanTransform,
+  createMinTransform,
+  createNullishTransform,
+  createOrTransform,
+  createPrintTransform,
+  createQuantileSeqTransform,
+  createRangeTransform,
+  createRowTransform,
+  createStdTransform,
+  createSubsetTransform,
+  createSumTransform,
+  createVarianceTransform
+} from "@danielsimonjr/mathts-expression";
 var addScalar = createAddScalar(factoryScope);
 var multiplyScalar = createMultiplyScalar(
   factoryScope
@@ -40455,6 +40517,38 @@ var cumsum = factory_cumsum;
 var ctranspose = factory_ctranspose;
 var createUnit = factory_createUnit;
 Object.assign(_mathWithTransform, factoryScope);
+{
+  const transformFactories = [
+    createAndTransform,
+    createBitAndTransform,
+    createBitOrTransform,
+    createColumnTransform,
+    createConcatTransform,
+    createCumSumTransform,
+    createDiffTransform,
+    createFilterTransform,
+    createForEachTransform,
+    createIndexTransform,
+    createMapTransform,
+    createMapSlicesTransform,
+    createMaxTransform,
+    createMeanTransform,
+    createMinTransform,
+    createNullishTransform,
+    createOrTransform,
+    createPrintTransform,
+    createQuantileSeqTransform,
+    createRangeTransform,
+    createRowTransform,
+    createStdTransform,
+    createSubsetTransform,
+    createSumTransform,
+    createVarianceTransform
+  ];
+  for (const create of transformFactories) {
+    _mathWithTransform[create.fn] = create(factoryScope);
+  }
+}
 
 // src/factories/evaluate.ts
 var mathScope = {
@@ -40482,6 +40576,38 @@ var mathScope = {
   true: true,
   false: false
 };
+{
+  const mwt = factoryScope.mathWithTransform;
+  for (const key of [
+    "and",
+    "bitAnd",
+    "bitOr",
+    "column",
+    "concat",
+    "cumsum",
+    "diff",
+    "filter",
+    "forEach",
+    "index",
+    "map",
+    "mapSlices",
+    "max",
+    "mean",
+    "min",
+    "nullish",
+    "or",
+    "print",
+    "quantileSeq",
+    "range",
+    "row",
+    "std",
+    "subset",
+    "sum",
+    "variance"
+  ]) {
+    if (typeof mwt[key] === "function") mathScope[key] = mwt[key];
+  }
+}
 var parse = factoryScope.parse;
 var evaluate = createEvaluate(parse, mathScope);
 function compileExpr(expr) {