@danielsimonjr/mathts-functions 0.21.0 → 0.23.0

package/dist/arithmetic/norm.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"norm.d.ts","sourceRoot":"","sources":["../../src/arithmetic/norm.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EAAE,aAAa,EAAE,MAAM,2BAA2B,CAAC;AAO/D,UAAU,WAAW;IACnB,GAAG,IAAI,MAAM,CAAC;CACf;AAED,UAAU,aAAa;IACrB,GAAG,IAAI,aAAa,CAAC;CACtB;AAED,UAAU,UAAU;IAClB,IAAI,IAAI,MAAM,EAAE,CAAC;IACjB,OAAO,CACL,QAAQ,EAAE,CACR,KAAK,EAAE,MAAM,GAAG,aAAa,GAAG,WAAW,EAC3C,KAAK,EAAE,MAAM,EAAE,EACf,MAAM,EAAE,UAAU,KACf,IAAI,EACT,SAAS,CAAC,EAAE,OAAO,GAClB,IAAI,CAAC;IACR,OAAO,IAAI,CAAC,MAAM,GAAG,aAAa,GAAG,WAAW,CAAC,EAAE,CAAC;CACrD;AAED,UAAU,UAAU;IAClB,MAAM,EAAE,UAAU,CAAC;CACpB;AAED,UAAU,gBAAgB;IACxB,KAAK,EAAE,aAAa,CAAC;IACrB,GAAG,EAAE,aAAa,CAAC;IACnB,GAAG,EAAE,aAAa,CAAC;IACnB,GAAG,EAAE,aAAa,CAAC;IACnB,IAAI,EAAE,aAAa,CAAC;IACpB,IAAI,EAAE,aAAa,CAAC;IACpB,QAAQ,EAAE,aAAa,CAAC;IACxB,WAAW,EAAE,aAAa,CAAC;IAC3B,MAAM,EAAE,CAAC,CAAC,EAAE,MAAM,GAAG,aAAa,EAAE,CAAC,EAAE,MAAM,GAAG,aAAa,KAAK,OAAO,CAAC;IAC1E,OAAO,EAAE,CAAC,CAAC,EAAE,MAAM,GAAG,aAAa,EAAE,CAAC,EAAE,MAAM,GAAG,aAAa,KAAK,OAAO,CAAC;IAC3E,MAAM,EAAE,CAAC,IAAI,EAAE,OAAO,KAAK,UAAU,CAAC;IACtC,UAAU,EAAE,aAAa,CAAC;IAC1B,IAAI,EAAE,CAAC,CAAC,EAAE,UAAU,KAAK,UAAU,CAAC;CACrC;AAmBD,eAAO,MAAM,UAAU,gFA+VtB,CAAC"}
+{"version":3,"file":"norm.d.ts","sourceRoot":"","sources":["../../src/arithmetic/norm.ts"],"names":[],"mappings":"AAEA,OAAO,KAAK,EAAE,aAAa,EAAE,MAAM,2BAA2B,CAAC;AAO/D,UAAU,WAAW;IACnB,GAAG,IAAI,MAAM,CAAC;CACf;AAED,UAAU,aAAa;IACrB,GAAG,IAAI,aAAa,CAAC;CACtB;AAED,UAAU,UAAU;IAClB,IAAI,IAAI,MAAM,EAAE,CAAC;IACjB,OAAO,CACL,QAAQ,EAAE,CACR,KAAK,EAAE,MAAM,GAAG,aAAa,GAAG,WAAW,EAC3C,KAAK,EAAE,MAAM,EAAE,EACf,MAAM,EAAE,UAAU,KACf,IAAI,EACT,SAAS,CAAC,EAAE,OAAO,GAClB,IAAI,CAAC;IACR,OAAO,IAAI,CAAC,MAAM,GAAG,aAAa,GAAG,WAAW,CAAC,EAAE,CAAC;CACrD;AAED,UAAU,UAAU;IAClB,MAAM,EAAE,UAAU,CAAC;CACpB;AAED,UAAU,gBAAgB;IACxB,KAAK,EAAE,aAAa,CAAC;IACrB,GAAG,EAAE,aAAa,CAAC;IACnB,GAAG,EAAE,aAAa,CAAC;IACnB,GAAG,EAAE,aAAa,CAAC;IACnB,IAAI,EAAE,aAAa,CAAC;IACpB,IAAI,EAAE,aAAa,CAAC;IACpB,QAAQ,EAAE,aAAa,CAAC;IACxB,WAAW,EAAE,aAAa,CAAC;IAC3B,MAAM,EAAE,CAAC,CAAC,EAAE,MAAM,GAAG,aAAa,EAAE,CAAC,EAAE,MAAM,GAAG,aAAa,KAAK,OAAO,CAAC;IAC1E,OAAO,EAAE,CAAC,CAAC,EAAE,MAAM,GAAG,aAAa,EAAE,CAAC,EAAE,MAAM,GAAG,aAAa,KAAK,OAAO,CAAC;IAC3E,MAAM,EAAE,CAAC,IAAI,EAAE,OAAO,KAAK,UAAU,CAAC;IACtC,UAAU,EAAE,aAAa,CAAC;IAC1B,IAAI,EAAE,CAAC,CAAC,EAAE,UAAU,KAAK,UAAU,CAAC;CACrC;AAmBD,eAAO,MAAM,UAAU,gFAiXtB,CAAC"}

package/dist/geometry/distance.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"distance.d.ts","sourceRoot":"","sources":["../../src/geometry/distance.ts"],"names":[],"mappings":"AAGA,OAAO,KAAK,EAAE,eAAe,EAAE,MAAM,aAAa,CAAC;AACnD,OAAO,KAAK,EAAE,aAAa,EAAE,MAAM,2BAA2B,CAAC;AAM/D,UAAU,oBAAoB;IAC5B,KAAK,EAAE,aAAa,CAAC;IACrB,SAAS,EAAE,CAAC,CAAC,EAAE,eAAe,EAAE,CAAC,EAAE,eAAe,KAAK,eAAe,CAAC;IACvE,cAAc,EAAE,CAAC,CAAC,EAAE,eAAe,EAAE,CAAC,EAAE,eAAe,KAAK,eAAe,CAAC;IAC5E,cAAc,EAAE,CAAC,CAAC,EAAE,eAAe,EAAE,CAAC,EAAE,eAAe,KAAK,eAAe,CAAC;IAC5E,YAAY,EAAE,CAAC,CAAC,EAAE,eAAe,EAAE,CAAC,EAAE,eAAe,KAAK,eAAe,CAAC;IAC1E,SAAS,EAAE,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,KAAK,OAAO,CAAC;IAC/C,IAAI,EAAE,CAAC,CAAC,EAAE,eAAe,KAAK,eAAe,CAAC;IAC9C,GAAG,EAAE,CAAC,CAAC,EAAE,eAAe,KAAK,eAAe,CAAC;CAC9C;AAcD,eAAO,MAAM,cAAc,8EAshB1B,CAAC"}
+{"version":3,"file":"distance.d.ts","sourceRoot":"","sources":["../../src/geometry/distance.ts"],"names":[],"mappings":"AAIA,OAAO,KAAK,EAAE,eAAe,EAAE,MAAM,aAAa,CAAC;AACnD,OAAO,KAAK,EAAE,aAAa,EAAE,MAAM,2BAA2B,CAAC;AAG/D,UAAU,oBAAoB;IAC5B,KAAK,EAAE,aAAa,CAAC;IACrB,SAAS,EAAE,CAAC,CAAC,EAAE,eAAe,EAAE,CAAC,EAAE,eAAe,KAAK,eAAe,CAAC;IACvE,cAAc,EAAE,CAAC,CAAC,EAAE,eAAe,EAAE,CAAC,EAAE,eAAe,KAAK,eAAe,CAAC;IAC5E,cAAc,EAAE,CAAC,CAAC,EAAE,eAAe,EAAE,CAAC,EAAE,eAAe,KAAK,eAAe,CAAC;IAC5E,YAAY,EAAE,CAAC,CAAC,EAAE,eAAe,EAAE,CAAC,EAAE,eAAe,KAAK,eAAe,CAAC;IAC1E,SAAS,EAAE,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,KAAK,OAAO,CAAC;IAC/C,IAAI,EAAE,CAAC,CAAC,EAAE,eAAe,KAAK,eAAe,CAAC;IAC9C,GAAG,EAAE,CAAC,CAAC,EAAE,eAAe,KAAK,eAAe,CAAC;CAC9C;AAcD,eAAO,MAAM,cAAc,8EAkgB1B,CAAC"}

package/dist/index.js

@@ -194,6 +194,7 @@ __export(typed_exports, {
   fresnelC: () => fresnelC,
   fresnelS: () => fresnelS,
   friedmanTest: () => friedmanTest,
+  fsum: () => fsum,
   fullSimplify: () => fullSimplify,
   functionExpand: () => functionExpand,
   fuseUnaryChain: () => fuseUnaryChain,
@@ -313,7 +314,7 @@ __export(typed_exports, {
   nintegrate: () => nintegrate,
   noncentralChi2PDF: () => noncentralChi2PDF,
   norm: () => norm,
-  norm2: () => norm2,
+  norm2: () => norm22,
   normFro: () => normFro,
   normalCDF: () => normalCDF,
   normalDist: () => normalDist,
@@ -509,6 +510,7 @@ __export(typed_exports, {
 
 // src/typed/arithmetic.ts
 import { mathTyped as mathTyped2, Complex, Fraction, BigNumber, Dual } from "@danielsimonjr/mathts-core";
+import { pairwiseSum, neumaierSum, norm2, pairwiseDot } from "@danielsimonjr/mathts-core";
 import { DenseMatrix as DenseMatrix2, backendManager as backendManager2, singularValues } from "@danielsimonjr/mathts-matrix";
 
 // src/factories/matrix-bridge.ts
@@ -2020,11 +2022,10 @@ function vectorNorm(arr7, p = 2) {
   if (p === Infinity) return Math.max(...arr7.map((x) => Math.abs(x)));
   if (p === -Infinity) return Math.min(...arr7.map((x) => Math.abs(x)));
   if (p === 0) return arr7.filter((x) => x !== 0).length;
-  if (p === 2) return Math.sqrt(arr7.reduce((sum3, x) => sum3 + x * x, 0));
-  return Math.pow(
-    arr7.reduce((sum3, x) => sum3 + Math.pow(Math.abs(x), p), 0),
-    1 / p
-  );
+  if (p === 2) return norm2(arr7);
+  const terms = new Array(arr7.length);
+  for (let i = 0; i < arr7.length; i++) terms[i] = Math.pow(Math.abs(arr7[i]), p);
+  return Math.pow(pairwiseSum(terms), 1 / p);
 }
 function matrixNorm(data, p = "fro") {
   const rows = data.length;
@@ -2081,9 +2082,7 @@ var norm = mathTyped2("norm", {
     if (computePool.shouldParallelize(a.length, "norm")) {
       return (await computePool.norm(a)).result;
     }
-    let s = 0;
-    for (let i = 0; i < a.length; i++) s += a[i] * a[i];
-    return Math.sqrt(s);
+    return norm2(a);
   }
 });
 var sinh = mathTyped2("sinh", {
@@ -2215,15 +2214,28 @@ var max = mathTyped2("max", {
     return m;
   }
 });
+var fsum = mathTyped2("fsum", {
+  Array: (arr7) => neumaierSum(arr7),
+  Float64Array: (a) => neumaierSum(a)
+});
 var sum = mathTyped2("sum", {
-  Array: (arr7) => {
-    let total = 0;
-    for (let i = 0; i < arr7.length; i++) {
-      total += arr7[i];
-    }
-    return total;
-  },
-  // Parallel Float64Array sum
+  // PAIRWISE summation, not `total += arr[i]`.
+  //
+  // Naive accumulation lets the running total grow large while the addends stay small, so each
+  // addition rounds off a bit more of the total: error grows as O(n)*eps. Measured on 1e6 copies
+  // of 0.1 (exact answer 100000):
+  //
+  //     naive (what shipped)  relative error 1.3e-11
+  //     pairwise (this)       relative error 2.9e-16   <- identical to NumPy's np.sum
+  //
+  // We were ~46,000x less accurate than NumPy on a bog-standard `sum`, and `mean`/`std`/`variance`
+  // all inherit it. Pairwise costs the same number of additions — the naive loop was simply worse.
+  // For catastrophic cancellation (e.g. [1e16, 1, -1e16], which pairwise AND np.sum both
+  // annihilate to 0) use `fsum`.
+  Array: (arr7) => pairwiseSum(arr7),
+  // Float64Array: the pool's sequential branch is pairwise too (see ComputePool.sum), and its
+  // parallel branch sums per-chunk then combines — blocked summation, which is pairwise-like.
+  // Either way the O(n) drift is gone.
   Float64Array: async (a) => {
     const result = await computePool.sum(a);
     return result.result;
@@ -2232,11 +2244,7 @@ var sum = mathTyped2("sum", {
 var mean = mathTyped2("mean", {
   Array: (arr7) => {
     if (arr7.length === 0) return NaN;
-    let total = 0;
-    for (let i = 0; i < arr7.length; i++) {
-      total += arr7[i];
-    }
-    return total / arr7.length;
+    return pairwiseSum(arr7) / arr7.length;
   },
   // Parallel Float64Array mean
   Float64Array: async (a) => {
@@ -2302,11 +2310,7 @@ var dot = mathTyped2("dot", {
     if (a.length !== b.length) {
       throw new Error(`Vector lengths must match: ${a.length} vs ${b.length}`);
     }
-    let total = 0;
-    for (let i = 0; i < a.length; i++) {
-      total += a[i] * b[i];
-    }
-    return total;
+    return pairwiseDot(a, b);
   },
   // Parallel Float64Array dot product
   "Float64Array, Float64Array": async (a, b) => {
@@ -2669,7 +2673,7 @@ var typedTrigonometry = {
 };
 
 // src/typed/statistics.ts
-import { mathTyped as mathTyped4 } from "@danielsimonjr/mathts-core";
+import { mathTyped as mathTyped4, neumaierCumsum } from "@danielsimonjr/mathts-core";
 import { computePool as computePool3 } from "@danielsimonjr/mathts-parallel";
 function toFloat64Array(arr7) {
   return new Float64Array(arr7);
@@ -3023,24 +3027,16 @@ var parallelStatMAD = mathTyped4("parallelStatMAD", {
   }
 });
 var parallelStatCumsum = mathTyped4("parallelStatCumsum", {
-  // Float64Array
+  // Float64Array — Neumaier-compensated prefix scan (exact where np.cumsum drifts O(n)·ε).
   Float64Array: (data) => {
     const result = new Float64Array(data.length);
-    let sum3 = 0;
-    for (let i = 0; i < data.length; i++) {
-      sum3 += data[i];
-      result[i] = sum3;
-    }
+    neumaierCumsum(data, result);
     return result;
   },
   // Number array
   Array: (arr7) => {
-    const result = [];
-    let sum3 = 0;
-    for (let i = 0; i < arr7.length; i++) {
-      sum3 += arr7[i];
-      result.push(sum3);
-    }
+    const result = new Array(arr7.length);
+    neumaierCumsum(arr7, result);
     return result;
   }
 });
@@ -16247,7 +16243,7 @@ var pinv = mathTyped13("pinv", {
 var cond = mathTyped13("cond", {
   Array: (A) => matrixCond(A)
 });
-var norm2 = mathTyped13("norm2", {
+var norm22 = mathTyped13("norm2", {
   Array: (A) => matrixNorm2(A)
 });
 var normFro = mathTyped13("normFro", {
@@ -21182,6 +21178,13 @@ var createUnaryPlus = /* @__PURE__ */ factory(
 );
 
 // src/matrix/dot.ts
+import { pairwiseDot as pairwiseDot2 } from "@danielsimonjr/mathts-core";
+function _isFlatNumbers(arr7) {
+  for (let i = 0; i < arr7.length; i++) {
+    if (typeof arr7[i] !== "number") return false;
+  }
+  return true;
+}
 var name72 = "dot";
 var dependencies72 = ["typed", "addScalar", "multiplyScalar", "conj", "size"];
 var createDot = /* @__PURE__ */ factory(
@@ -21234,6 +21237,9 @@ var createDot = /* @__PURE__ */ factory(
         add2 = typed3.find(addScalar2, [dt, dt]);
         mul = typed3.find(multiplyScalar2, [dt, dt]);
       }
+      if (!aIsColumn && !bIsColumn && add2 === addScalar2 && mul === multiplyScalar2 && _isFlatNumbers(adata) && _isFlatNumbers(bdata)) {
+        return pairwiseDot2(adata, bdata);
+      }
       if (!aIsColumn && !bIsColumn) {
         let c = mul(conj3(adata[0]), bdata[0]);
         for (let i = 1; i < N; i++) {
@@ -26235,7 +26241,7 @@ var createZpk2tf = /* @__PURE__ */ factory(
 );
 
 // src/statistics/cumsum.ts
-var WASM_CUMSUM_THRESHOLD = 100;
+import { neumaierCumsum as neumaierCumsum2 } from "@danielsimonjr/mathts-core";
 function isFlatNumberArray2(arr7) {
   for (let i = 0; i < arr7.length; i++) {
     if (typeof arr7[i] !== "number") {
@@ -26283,24 +26289,10 @@ var createCumSum = /* @__PURE__ */ factory(
       if (array.length === 0) {
         return [];
       }
-      if (array.length >= WASM_CUMSUM_THRESHOLD && isFlatNumberArray2(array)) {
-        const wasm = wasmLoader.getModule();
-        if (wasm) {
-          try {
-            const alloc = wasmLoader.allocateFloat64Array(array);
-            try {
-              wasm.statsCumsum(alloc.ptr, array.length);
-              const result = new Array(array.length);
-              for (let i = 0; i < array.length; i++) {
-                result[i] = alloc.array[i];
-              }
-              return result;
-            } finally {
-              wasmLoader.free(alloc.ptr);
-            }
-          } catch {
-          }
-        }
+      if (isFlatNumberArray2(array)) {
+        const result = new Array(array.length);
+        neumaierCumsum2(array, result);
+        return result;
       }
       const sums = [unaryPlus2(array[0])];
       for (let i = 1; i < array.length; ++i) {
@@ -26348,7 +26340,7 @@ var createCumSum = /* @__PURE__ */ factory(
 );
 
 // src/statistics/sum.ts
-var WASM_SUM_THRESHOLD = 100;
+import { pairwiseSum as pairwiseSum2 } from "@danielsimonjr/mathts-core";
 function isFlatNumberArray3(arr7) {
   for (let i = 0; i < arr7.length; i++) {
     if (typeof arr7[i] !== "number") {
@@ -26381,21 +26373,8 @@ var createSum = /* @__PURE__ */ factory(
       }
     });
     function _sum(array) {
-      if (Array.isArray(array) && array.length >= WASM_SUM_THRESHOLD) {
-        if (isFlatNumberArray3(array)) {
-          const wasm = wasmLoader.getModule();
-          if (wasm) {
-            try {
-              const alloc = wasmLoader.allocateFloat64Array(array);
-              try {
-                return wasm.statsSum(alloc.ptr, array.length);
-              } finally {
-                wasmLoader.free(alloc.ptr);
-              }
-            } catch {
-            }
-          }
-        }
+      if (Array.isArray(array) && isFlatNumberArray3(array)) {
+        return pairwiseSum2(array);
       }
       let sum3;
       deepForEach2(array, function(value) {
@@ -30384,7 +30363,7 @@ var createHypot = /* @__PURE__ */ factory(
 );
 
 // src/geometry/distance.ts
-var WASM_DISTANCE_THRESHOLD = 4;
+import { scaledDistance } from "@danielsimonjr/mathts-core";
 var name181 = "distance";
 var dependencies181 = [
   "typed",
@@ -30669,10 +30648,14 @@ var createDistance = /* @__PURE__ */ factory(
     }
     function _euclideanDistance(x, y) {
       const vectorSize = x.length;
-      if (vectorSize >= WASM_DISTANCE_THRESHOLD) {
-        const wasmResult = _tryWasmDistanceND(x, y, vectorSize);
-        if (wasmResult !== null) return wasmResult;
+      let allNumbers = true;
+      for (let i = 0; i < vectorSize; i++) {
+        if (typeof x[i] !== "number" || typeof y[i] !== "number") {
+          allNumbers = false;
+          break;
+        }
       }
+      if (allNumbers) return scaledDistance(x, y);
       let result = 0;
       for (let i = 0; i < vectorSize; i++) {
         const diff2 = subtractScalar2(x[i], y[i]);
@@ -30680,21 +30663,6 @@ var createDistance = /* @__PURE__ */ factory(
       }
       return sqrt2(result);
     }
-    function _tryWasmDistanceND(x, y, n) {
-      const wasm = wasmLoader.getModule();
-      if (!wasm) return null;
-      for (let i = 0; i < n; i++) {
-        if (typeof x[i] !== "number" || typeof y[i] !== "number") return null;
-      }
-      const p1Alloc = wasmLoader.allocateFloat64Array(x);
-      const p2Alloc = wasmLoader.allocateFloat64Array(y);
-      try {
-        return wasm.distanceND(p1Alloc.ptr, p2Alloc.ptr, n);
-      } finally {
-        wasmLoader.free(p1Alloc.ptr);
-        wasmLoader.free(p2Alloc.ptr);
-      }
-    }
     function _distancePairwise(a) {
       const wasmResult = _tryWasmPairwise(a);
       if (wasmResult !== null) return wasmResult;
@@ -39726,6 +39694,7 @@ var createDerivative = /* @__PURE__ */ factory(
 );
 
 // src/arithmetic/norm.ts
+import { norm2 as norm23 } from "@danielsimonjr/mathts-core";
 var WASM_NORM_THRESHOLD = 100;
 var name245 = "norm";
 var dependencies245 = [
@@ -39922,6 +39891,12 @@ var createNorm = /* @__PURE__ */ factory(
       throw new Error("Unsupported parameter value " + p);
     }
     function _norm(x, p) {
+      if ((p === 2 || p === "fro") && x.size().length === 1) {
+        const v = x.valueOf();
+        if (v.every((e) => typeof e === "number")) {
+          return norm23(v);
+        }
+      }
       const sizeX = x.size();
       if (sizeX.length === 1) {
         return _vectorNorm(x, p);
@@ -44854,7 +44829,7 @@ function levenbergMarquardt(residual, x0, opts = {}) {
   let x = Array.from(x0);
   let lambda = 1e-3;
   const h = 1e-7;
-  const norm22 = (v) => v.reduce((s, c) => s + c * c, 0);
+  const norm24 = (v) => v.reduce((s, c) => s + c * c, 0);
   const jacobian2 = (xc, r0) => {
     const m = r0.length;
     const J = Array.from({ length: m }, () => new Array(n).fill(0));
@@ -44867,7 +44842,7 @@ function levenbergMarquardt(residual, x0, opts = {}) {
     return J;
   };
   let r = residual(x);
-  let cost = norm22(r);
+  let cost = norm24(r);
   let iter = 0;
   for (; iter < maxIter; iter++) {
     const J = jacobian2(x, r);
@@ -44898,7 +44873,7 @@ function levenbergMarquardt(residual, x0, opts = {}) {
     }
     const xNew = x.map((v, i) => v + delta[i]);
     const rNew = residual(xNew);
-    const costNew = norm22(rNew);
+    const costNew = norm24(rNew);
     if (costNew < cost) {
       x = xNew;
       r = rNew;
@@ -45661,6 +45636,7 @@ export {
   fresnelC,
   fresnelS,
   friedmanTest,
+  fsum,
   fullSimplify,
   functionExpand,
   fuseUnaryChain,
@@ -45882,7 +45858,7 @@ export {
   nodeOperations,
   noncentralChi2PDF,
   norm,
-  norm2,
+  norm22 as norm2,
   normFro,
   normalCDF,
   normalDist,

package/dist/matrix/dot.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"dot.d.ts","sourceRoot":"","sources":["../../src/matrix/dot.ts"],"names":[],"mappings":"AAEA,OAAO,KAAK,EAAE,aAAa,EAAE,MAAM,2BAA2B,CAAC;AA8B/D,UAAU,YAAY;IACpB,KAAK,EAAE,aAAa,CAAC;IACrB,SAAS,EAAE,aAAa,CAAC;IACzB,cAAc,EAAE,aAAa,CAAC;IAC9B,IAAI,EAAE,aAAa,CAAC;IACpB,IAAI,EAAE,aAAa,CAAC;CACrB;AAKD,eAAO,MAAM,SAAS,sEAmMrB,CAAC"}
+{"version":3,"file":"dot.d.ts","sourceRoot":"","sources":["../../src/matrix/dot.ts"],"names":[],"mappings":"AAGA,OAAO,KAAK,EAAE,aAAa,EAAE,MAAM,2BAA2B,CAAC;AAsC/D,UAAU,YAAY;IACpB,KAAK,EAAE,aAAa,CAAC;IACrB,SAAS,EAAE,aAAa,CAAC;IACzB,cAAc,EAAE,aAAa,CAAC;IAC9B,IAAI,EAAE,aAAa,CAAC;IACpB,IAAI,EAAE,aAAa,CAAC;CACrB;AAKD,eAAO,MAAM,SAAS,sEAiNrB,CAAC"}

package/dist/statistics/cumsum.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"cumsum.d.ts","sourceRoot":"","sources":["../../src/statistics/cumsum.ts"],"names":[],"mappings":"AAOA,OAAO,KAAK,EAAE,aAAa,EAAE,MAAM,2BAA2B,CAAC;AAwB/D,UAAU,kBAAkB;IAC1B,KAAK,EAAE,aAAa,CAAC;IACrB,GAAG,EAAE,aAAa,CAAC;IACnB,SAAS,EAAE,aAAa,CAAC;CAC1B;AAKD,eAAO,MAAM,YAAY,4EAgLxB,CAAC"}
+{"version":3,"file":"cumsum.d.ts","sourceRoot":"","sources":["../../src/statistics/cumsum.ts"],"names":[],"mappings":"AAOA,OAAO,KAAK,EAAE,aAAa,EAAE,MAAM,2BAA2B,CAAC;AAqB/D,UAAU,kBAAkB;IAC1B,KAAK,EAAE,aAAa,CAAC;IACrB,GAAG,EAAE,aAAa,CAAC;IACnB,SAAS,EAAE,aAAa,CAAC;CAC1B;AAKD,eAAO,MAAM,YAAY,4EAiKxB,CAAC"}

package/dist/statistics/sum.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"sum.d.ts","sourceRoot":"","sources":["../../src/statistics/sum.ts"],"names":[],"mappings":"AAIA,OAAO,KAAK,EAAE,aAAa,EAAE,MAAM,2BAA2B,CAAC;AAC/D,OAAO,KAAK,EAAE,aAAa,EAAE,MAAM,mBAAmB,CAAC;AA2BvD,UAAU,eAAe;IACvB,KAAK,EAAE,aAAa,CAAC;IACrB,MAAM,EAAE,aAAa,CAAC;IACtB,GAAG,EAAE,aAAa,CAAC;IACnB,OAAO,EAAE,aAAa,CAAC;IACvB,qBAAqB,EAAE,CAAC,KAAK,EAAE,MAAM,KAAK,OAAO,CAAC;CACnD;AAKD,eAAO,MAAM,SAAS,yEAmHrB,CAAC"}
+{"version":3,"file":"sum.d.ts","sourceRoot":"","sources":["../../src/statistics/sum.ts"],"names":[],"mappings":"AAIA,OAAO,KAAK,EAAE,aAAa,EAAE,MAAM,2BAA2B,CAAC;AAC/D,OAAO,KAAK,EAAE,aAAa,EAAE,MAAM,mBAAmB,CAAC;AA0BvD,UAAU,eAAe;IACvB,KAAK,EAAE,aAAa,CAAC;IACrB,MAAM,EAAE,aAAa,CAAC;IACtB,GAAG,EAAE,aAAa,CAAC;IACnB,OAAO,EAAE,aAAa,CAAC;IACvB,qBAAqB,EAAE,CAAC,KAAK,EAAE,MAAM,KAAK,OAAO,CAAC;CACnD;AAKD,eAAO,MAAM,SAAS,yEAmHrB,CAAC"}

package/dist/typed/arithmetic.d.ts

@@ -186,6 +186,20 @@ export declare const max: import("typed-function").TypedFunction;
 /**
  * Sum with parallel array support
  */
+/**
+ * Exactly-rounded sum — the equivalent of Python's `math.fsum` / NumPy's `math.fsum` idiom.
+ *
+ * `sum` uses pairwise summation, which is accurate to ~machine epsilon and free. But pairwise
+ * cannot recover a value that has already been annihilated by catastrophic cancellation:
+ *
+ *     sum([1e16, 1, -1e16])   ->  0     (so does np.sum — the 1 is lost the moment it meets 1e16)
+ *     fsum([1e16, 1, -1e16])  ->  1     (exact)
+ *
+ * Neumaier compensation tracks the low-order bits each addition discards and folds them back in.
+ * ~2-4x slower than `sum`, so it is opt-in. Reach for it when the result is a small difference of
+ * large terms: conservation checks, residuals, long-running accumulators.
+ */
+export declare const fsum: import("typed-function").TypedFunction;
 export declare const sum: import("typed-function").TypedFunction;
 /**
  * Mean with parallel array support

package/dist/typed/arithmetic.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"arithmetic.d.ts","sourceRoot":"","sources":["../../src/typed/arithmetic.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAQH,OAAO,EAAe,WAAW,EAAE,MAAM,gCAAgC,CAAC;AAU1E,4BAA4B;AAC5B,KAAK,GAAG,GAAG,MAAM,CAAC;AAkIlB;;;;;;;;;;;;;GAaG;AACH,eAAO,MAAM,GAAG,wCAgEd,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,QAAQ,wCAiCnB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,QAAQ,wCA8FnB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,MAAM,wCAkCjB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,UAAU,wCAYrB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,SAAS,wCAOpB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAgBd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAqBf,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCA2Cd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAgBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,MAAM,wCAajB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAkBf,CAAC;AAaH;;GAEG;AACH,eAAO,MAAM,IAAI,wCA4Bf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,OAAO,wCAIlB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAad,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAcd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAiBhB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAiBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAehB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAehB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAoBhB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAgBhB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAgBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAgBd,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAKd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCA+Bd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAgBd,CAAC;AAEH;;;GAGG;AACH,eAAO,MAAM,IAAI,wCA6Bf,CAAC;AAgEH,eAAO,MAAM,IAAI,wCAmCf,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAiBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAiBf,CAAC;AAMH;;GAEG;AAIH,eAAO,MAAM,KAAK,wCAQhB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,OAAO,wCAMlB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,MAAM,wCAMjB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,SAAS,wCAMpB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,QAAQ,wCAMnB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,OAAO,wCAMlB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAmBd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAkBd,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAcd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAEH;;;;;GAKG;AACH,MAAM,MAAM,gBAAgB,GAAG,UAAU,GAAG,aAAa,GAAG,QAAQ,CAAC;AA6BrE;;;;;;GAMG;AACH,eAAO,MAAM,QAAQ,wCAenB,CAAC;AAEH;;;GAGG;AACH,eAAO,MAAM,GAAG,wCAiBd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAiBd,CAAC;AAMH;;;;;;;GAOG;AACH,wBAAsB,MAAM,CAC1B,CAAC,EAAE,YAAY,EACf,KAAK,EAAE,GAAG,EACV,KAAK,EAAE,GAAG,EACV,CAAC,EAAE,YAAY,EACf,KAAK,EAAE,GAAG,GACT,OAAO,CAAC,YAAY,CAAC,CAGvB;AAED;;GAEG;AACH,wBAAsB,SAAS,CAAC,IAAI,EAAE,YAAY,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,GAAG,OAAO,CAAC,YAAY,CAAC,CAG/F;AAED;;GAEG;AACH,wBAAsB,MAAM,CAC1B,MAAM,EAAE,YAAY,EACpB,IAAI,EAAE,GAAG,EACT,IAAI,EAAE,GAAG,EACT,MAAM,EAAE,YAAY,GACnB,OAAO,CAAC,YAAY,CAAC,CAGvB;AAED;;GAEG;AACH,wBAAsB,KAAK,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,GAAG,OAAO,CAAC,YAAY,CAAC,CAGnF;AAMD;;GAEG;AACH,wBAAsB,cAAc,IAAI,OAAO,CAAC,IAAI,CAAC,CAEpD;AAED;;GAEG;AACH,wBAAsB,aAAa,IAAI,OAAO,CAAC,IAAI,CAAC,CAEnD;AAED;;GAEG;AACH,wBAAgB,iBAAiB,CAAC,YAAY,EAAE,GAAG,GAAG,OAAO,CAE5D;AAED;;GAEG;AACH,wBAAgB,cAAc,IAAI,WAAW,CAE5C;AAMD,eAAO,MAAM,eAAe;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;CA+E3B,CAAC"}
+{"version":3,"file":"arithmetic.d.ts","sourceRoot":"","sources":["../../src/typed/arithmetic.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AASH,OAAO,EAAe,WAAW,EAAE,MAAM,gCAAgC,CAAC;AAU1E,4BAA4B;AAC5B,KAAK,GAAG,GAAG,MAAM,CAAC;AAkIlB;;;;;;;;;;;;;GAaG;AACH,eAAO,MAAM,GAAG,wCAgEd,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,QAAQ,wCAiCnB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,QAAQ,wCA8FnB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,MAAM,wCAkCjB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,UAAU,wCAYrB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,SAAS,wCAOpB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAgBd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAqBf,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCA2Cd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAgBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,MAAM,wCAajB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAkBf,CAAC;AAaH;;GAEG;AACH,eAAO,MAAM,IAAI,wCA4Bf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,OAAO,wCAIlB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAad,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAcd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAiBhB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAiBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAehB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAehB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAoBhB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAgBhB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAgBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAgBd,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAKd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCA+Bd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAgBd,CAAC;AAEH;;;GAGG;AACH,eAAO,MAAM,IAAI,wCA6Bf,CAAC;AAwEH,eAAO,MAAM,IAAI,wCAmCf,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAiBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAiBf,CAAC;AAMH;;GAEG;AAIH,eAAO,MAAM,KAAK,wCAQhB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,OAAO,wCAMlB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,MAAM,wCAMjB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,SAAS,wCAMpB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,QAAQ,wCAMnB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,OAAO,wCAMlB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAmBd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAkBd,CAAC;AAMH;;GAEG;AACH;;;;;;;;;;;;GAYG;AACH,eAAO,MAAM,IAAI,wCAGf,CAAC;AAEH,eAAO,MAAM,GAAG,wCAuBd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAWf,CAAC;AAEH;;;;;GAKG;AACH,MAAM,MAAM,gBAAgB,GAAG,UAAU,GAAG,aAAa,GAAG,QAAQ,CAAC;AA6BrE;;;;;;GAMG;AACH,eAAO,MAAM,QAAQ,wCAenB,CAAC;AAEH;;;GAGG;AACH,eAAO,MAAM,GAAG,wCAiBd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAad,CAAC;AAMH;;;;;;;GAOG;AACH,wBAAsB,MAAM,CAC1B,CAAC,EAAE,YAAY,EACf,KAAK,EAAE,GAAG,EACV,KAAK,EAAE,GAAG,EACV,CAAC,EAAE,YAAY,EACf,KAAK,EAAE,GAAG,GACT,OAAO,CAAC,YAAY,CAAC,CAGvB;AAED;;GAEG;AACH,wBAAsB,SAAS,CAAC,IAAI,EAAE,YAAY,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,GAAG,OAAO,CAAC,YAAY,CAAC,CAG/F;AAED;;GAEG;AACH,wBAAsB,MAAM,CAC1B,MAAM,EAAE,YAAY,EACpB,IAAI,EAAE,GAAG,EACT,IAAI,EAAE,GAAG,EACT,MAAM,EAAE,YAAY,GACnB,OAAO,CAAC,YAAY,CAAC,CAGvB;AAED;;GAEG;AACH,wBAAsB,KAAK,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,GAAG,OAAO,CAAC,YAAY,CAAC,CAGnF;AAMD;;GAEG;AACH,wBAAsB,cAAc,IAAI,OAAO,CAAC,IAAI,CAAC,CAEpD;AAED;;GAEG;AACH,wBAAsB,aAAa,IAAI,OAAO,CAAC,IAAI,CAAC,CAEnD;AAED;;GAEG;AACH,wBAAgB,iBAAiB,CAAC,YAAY,EAAE,GAAG,GAAG,OAAO,CAE5D;AAED;;GAEG;AACH,wBAAgB,cAAc,IAAI,WAAW,CAE5C;AAMD,eAAO,MAAM,eAAe;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;CA+E3B,CAAC"}

package/dist/typed/statistics.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"statistics.d.ts","sourceRoot":"","sources":["../../src/typed/statistics.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAUH,0CAA0C;AAC1C,KAAK,GAAG,GAAG,MAAM,CAAC;AAElB,4BAA4B;AAC5B,KAAK,GAAG,GAAG,MAAM,CAAC;AAMlB;;GAEG;AACH,MAAM,MAAM,iBAAiB,GAAG,UAAU,GAAG,aAAa,GAAG,QAAQ,CAAC;AAiBtE;;GAEG;AACH,eAAO,MAAM,eAAe,wCAoB1B,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,gBAAgB,wCAqB3B,CAAC;AAMH;;;;GAIG;AACH,eAAO,MAAM,oBAAoB,wCA2D/B,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,eAAe,wCAgC1B,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,eAAe,wCAuB1B,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,eAAe,wCAuB1B,CAAC;AAEH;;;GAGG;AACH,eAAO,MAAM,kBAAkB,wCAe7B,CAAC;AAMH;;;;GAIG;AACH,eAAO,MAAM,kBAAkB,wCAgC7B,CAAC;AAMH;;;GAGG;AACH,eAAO,MAAM,gBAAgB,wCA2B3B,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,gBAAgB,wCA8B3B,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,gBAAgB,wCA4D3B,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,oBAAoB,wCAiB/B,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,gBAAgB,wCAgC3B,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,eAAe,wCAsB1B,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,kBAAkB,wCAsB7B,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,oBAAoB,wCAgC/B,CAAC;AAMH;;;;;;;;GAQG;AACH,wBAAgB,sBAAsB,CAAC,IAAI,EAAE,YAAY,EAAE,CAAC,EAAE,GAAG,GAAG,GAAG,CAGtE;AAMD;;GAEG;AACH,eAAO,MAAM,qBAAqB,wCAsBhC,CAAC;AAMH;;;;;;;GAOG;AACH,wBAAgB,WAAW,CAAC,GAAG,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,GAAG,GAAG,GAAG,CAMtD;AA8CD;;;;;GAKG;AACH,wBAAgB,YAAY,CAAC,GAAG,EAAE,MAAM,EAAE,GAAG,GAAG,CAc/C;AAED;;;;;;GAMG;AACH,wBAAgB,SAAS,CAAC,GAAG,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,GAAG,GAAG,MAAM,EAAE,CAWzD;AAED;;;;;;GAMG;AACH,wBAAgB,SAAS,CAAC,GAAG,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,GAAG,GAAG,MAAM,EAAE,CAWzD;AAMD;;GAEG;AACH,eAAO,MAAM,eAAe;;;;;;;;;;;;;;;;;;;;;;;CAuB3B,CAAC;AAEF;;GAEG;AACH,wBAAsB,oBAAoB,IAAI,OAAO,CAAC,IAAI,CAAC,CAE1D;AAED;;GAEG;AACH,wBAAsB,mBAAmB,IAAI,OAAO,CAAC,IAAI,CAAC,CAEzD"}
+{"version":3,"file":"statistics.d.ts","sourceRoot":"","sources":["../../src/typed/statistics.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAUH,0CAA0C;AAC1C,KAAK,GAAG,GAAG,MAAM,CAAC;AAElB,4BAA4B;AAC5B,KAAK,GAAG,GAAG,MAAM,CAAC;AAMlB;;GAEG;AACH,MAAM,MAAM,iBAAiB,GAAG,UAAU,GAAG,aAAa,GAAG,QAAQ,CAAC;AAiBtE;;GAEG;AACH,eAAO,MAAM,eAAe,wCAoB1B,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,gBAAgB,wCAqB3B,CAAC;AAMH;;;;GAIG;AACH,eAAO,MAAM,oBAAoB,wCA2D/B,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,eAAe,wCAgC1B,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,eAAe,wCAuB1B,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,eAAe,wCAuB1B,CAAC;AAEH;;;GAGG;AACH,eAAO,MAAM,kBAAkB,wCAe7B,CAAC;AAMH;;;;GAIG;AACH,eAAO,MAAM,kBAAkB,wCAgC7B,CAAC;AAMH;;;GAGG;AACH,eAAO,MAAM,gBAAgB,wCA2B3B,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,gBAAgB,wCA8B3B,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,gBAAgB,wCA4D3B,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,oBAAoB,wCAiB/B,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,gBAAgB,wCAgC3B,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,eAAe,wCAsB1B,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,kBAAkB,wCAc7B,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,oBAAoB,wCAgC/B,CAAC;AAMH;;;;;;;;GAQG;AACH,wBAAgB,sBAAsB,CAAC,IAAI,EAAE,YAAY,EAAE,CAAC,EAAE,GAAG,GAAG,GAAG,CAGtE;AAMD;;GAEG;AACH,eAAO,MAAM,qBAAqB,wCAsBhC,CAAC;AAMH;;;;;;;GAOG;AACH,wBAAgB,WAAW,CAAC,GAAG,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,GAAG,GAAG,GAAG,CAMtD;AA8CD;;;;;GAKG;AACH,wBAAgB,YAAY,CAAC,GAAG,EAAE,MAAM,EAAE,GAAG,GAAG,CAc/C;AAED;;;;;;GAMG;AACH,wBAAgB,SAAS,CAAC,GAAG,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,GAAG,GAAG,MAAM,EAAE,CAWzD;AAED;;;;;;GAMG;AACH,wBAAgB,SAAS,CAAC,GAAG,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,GAAG,GAAG,MAAM,EAAE,CAWzD;AAMD;;GAEG;AACH,eAAO,MAAM,eAAe;;;;;;;;;;;;;;;;;;;;;;;CAuB3B,CAAC;AAEF;;GAEG;AACH,wBAAsB,oBAAoB,IAAI,OAAO,CAAC,IAAI,CAAC,CAE1D;AAED;;GAEG;AACH,wBAAsB,mBAAmB,IAAI,OAAO,CAAC,IAAI,CAAC,CAEzD"}

package/package.json

@@ -1,76 +1,76 @@
-{
-  "name": "@danielsimonjr/mathts-functions",
-  "version": "0.21.0",
-  "description": "Mathematical functions for MathTS - arithmetic, algebra, trigonometry, statistics, and more",
-  "author": "Daniel Simon Jr.",
-  "license": "MIT",
-  "type": "module",
-  "main": "./dist/index.js",
-  "module": "./dist/index.js",
-  "types": "./dist/index.d.ts",
-  "exports": {
-    ".": {
-      "import": "./dist/index.js",
-      "types": "./dist/index.d.ts"
-    }
-  },
-  "files": [
-    "dist",
-    "types",
-    "README.md"
-  ],
-  "scripts": {
-    "build": "tsup src/index.ts --format esm --clean && tsc -p tsconfig.dts.json && node scripts/copy-wasm.mjs",
-    "dev": "tsup src/index.ts --format esm --dts --watch",
-    "test": "vitest run",
-    "test:diff": "node tests/diff-special.test.mjs && node tests/diff-elementwise.test.mjs && node tests/diff-fusion.test.mjs",
-    "golden:gen": "python -X utf8 tests/golden/gen_special_goldens.py",
-    "test:watch": "vitest",
-    "test:coverage": "vitest run --coverage",
-    "typecheck": "tsc --noEmit",
-    "lint": "eslint src --ext .ts",
-    "lint:fix": "eslint src --ext .ts --fix",
-    "clean": "rm -rf dist",
-    "build:prod": "tsup src/index.ts --format esm --clean --minify --treeshake"
-  },
-  "dependencies": {
-    "@danielsimonjr/mathts-core": "^0.6.0",
-    "@danielsimonjr/mathts-expression": "^0.6.0",
-    "@danielsimonjr/mathts-gpu": "^0.2.0",
-    "@danielsimonjr/mathts-matrix": "^0.4.1",
-    "@danielsimonjr/mathts-parallel": "^0.3.4",
-    "bignumber.js": "^9.1.2",
-    "complex.js": "^2.2.5",
-    "decimal.js": "^10.4.3",
-    "escape-latex": "^1.2.0",
-    "fraction.js": "^5.2.1",
-    "javascript-natural-sort": "^0.7.1",
-    "seedrandom": "^3.0.5",
-    "tiny-emitter": "^2.1.0",
-    "typed-function": "github:danielsimonjr/typed-function"
-  },
-  "devDependencies": {
-    "@types/node": "^25.5.2",
-    "@webgpu/types": "^0.1.67",
-    "tsup": "^8.0.0",
-    "typescript": "^5.3.0",
-    "vitest": "^4.1.5"
-  },
-  "publishConfig": {
-    "access": "public"
-  },
-  "repository": {
-    "type": "git",
-    "url": "https://github.com/danielsimonjr/mathts",
-    "directory": "functions"
-  },
-  "keywords": [
-    "math",
-    "typescript",
-    "functions",
-    "arithmetic",
-    "algebra",
-    "trigonometry",
-    "statistics"
-  ]
-}
+{
+  "name": "@danielsimonjr/mathts-functions",
+  "version": "0.23.0",
+  "description": "Mathematical functions for MathTS - arithmetic, algebra, trigonometry, statistics, and more",
+  "author": "Daniel Simon Jr.",
+  "license": "MIT",
+  "type": "module",
+  "main": "./dist/index.js",
+  "module": "./dist/index.js",
+  "types": "./dist/index.d.ts",
+  "exports": {
+    ".": {
+      "import": "./dist/index.js",
+      "types": "./dist/index.d.ts"
+    }
+  },
+  "files": [
+    "dist",
+    "types",
+    "README.md"
+  ],
+  "scripts": {
+    "build": "tsup src/index.ts --format esm --clean && tsc -p tsconfig.dts.json && node scripts/copy-wasm.mjs",
+    "dev": "tsup src/index.ts --format esm --dts --watch",
+    "test": "vitest run",
+    "test:diff": "node tests/diff-special.test.mjs && node tests/diff-elementwise.test.mjs && node tests/diff-fusion.test.mjs",
+    "golden:gen": "python -X utf8 tests/golden/gen_special_goldens.py",
+    "test:watch": "vitest",
+    "test:coverage": "vitest run --coverage",
+    "typecheck": "tsc --noEmit",
+    "lint": "eslint src --ext .ts",
+    "lint:fix": "eslint src --ext .ts --fix",
+    "clean": "rm -rf dist",
+    "build:prod": "tsup src/index.ts --format esm --clean --minify --treeshake"
+  },
+  "dependencies": {
+    "@danielsimonjr/mathts-core": "^0.8.0",
+    "@danielsimonjr/mathts-expression": "^0.6.2",
+    "@danielsimonjr/mathts-gpu": "^0.2.0",
+    "@danielsimonjr/mathts-matrix": "^0.4.3",
+    "@danielsimonjr/mathts-parallel": "^0.5.0",
+    "bignumber.js": "^9.1.2",
+    "complex.js": "^2.2.5",
+    "decimal.js": "^10.4.3",
+    "escape-latex": "^1.2.0",
+    "fraction.js": "^5.2.1",
+    "javascript-natural-sort": "^0.7.1",
+    "seedrandom": "^3.0.5",
+    "tiny-emitter": "^2.1.0",
+    "typed-function": "github:danielsimonjr/typed-function"
+  },
+  "devDependencies": {
+    "@types/node": "^25.5.2",
+    "@webgpu/types": "^0.1.67",
+    "tsup": "^8.0.0",
+    "typescript": "^5.3.0",
+    "vitest": "^4.1.5"
+  },
+  "publishConfig": {
+    "access": "public"
+  },
+  "repository": {
+    "type": "git",
+    "url": "https://github.com/danielsimonjr/mathts",
+    "directory": "functions"
+  },
+  "keywords": [
+    "math",
+    "typescript",
+    "functions",
+    "arithmetic",
+    "algebra",
+    "trigonometry",
+    "statistics"
+  ]
+}