@danielsimonjr/mathts-functions 0.37.0 → 0.38.0

package/dist/index.js

@@ -6876,60 +6876,73 @@ function digammaScalar(x) {
   }
   return result;
 }
-function besselKAsym(nu, x) {
-  const mu = 4 * nu * nu;
-  let a = 1;
-  let xk = 1;
-  let sum3 = 1;
-  let prevMag = Infinity;
-  for (let k = 1; k <= 40; k++) {
-    const t1 = 2 * k - 1;
-    a = a * (mu - t1 * t1) / (8 * k);
-    xk *= x;
-    const t = a / xk;
-    if (Math.abs(t) > prevMag) break;
-    prevMag = Math.abs(t);
-    sum3 += t;
-  }
-  return Math.sqrt(Math.PI / (2 * x)) * Math.exp(-x) * sum3;
-}
-function besselK0Series(x) {
-  const z = 0.25 * x * x;
-  let t = 1;
-  let I0 = 1;
-  let h = 0;
-  let s = 0;
-  for (let k = 1; k <= 80; k++) {
-    t *= z / (k * k);
-    I0 += t;
-    h += 1 / k;
-    s += h * t;
-    if (k > 2 && t <= I0 * 1e-18) break;
-  }
-  return -(Math.log(0.5 * x) + 0.5772156649015329) * I0 + s;
-}
-function besselK1Series(x) {
-  const z = 0.25 * x * x;
-  let w = 1;
-  let S1 = 1;
-  let hk = 0;
-  let hk1 = 1;
-  let S2 = (hk + hk1) * 1;
-  for (let k = 1; k <= 80; k++) {
-    w *= z / (k * (k + 1));
-    S1 += w;
-    hk += 1 / k;
-    hk1 += 1 / (k + 1);
-    S2 += (hk + hk1) * w;
-    if (k > 2 && w <= S1 * 1e-18) break;
-  }
-  return 1 / x + 0.5 * x * (Math.log(0.5 * x) + 0.5772156649015329) * S1 - 0.25 * x * S2;
+function besselK01(x) {
+  const EPS4 = 1e-16;
+  const xi = 1 / x;
+  const xi2 = 2 * xi;
+  let rkmu;
+  let rk1;
+  if (x < 2) {
+    const x2 = 0.5 * x;
+    const d = -Math.log(x2);
+    let ff = -0.5772156649015329 + d;
+    let sum3 = ff;
+    let p = 0.5;
+    let q = 0.5;
+    let c = 1;
+    const dd = x2 * x2;
+    let sum1 = p;
+    for (let i = 1; i <= 1e4; i++) {
+      ff = (i * ff + p + q) / (i * i);
+      c *= dd / i;
+      p /= i;
+      q /= i;
+      const del = c * ff;
+      sum3 += del;
+      sum1 += c * (p - i * ff);
+      if (Math.abs(del) < Math.abs(sum3) * EPS4) break;
+    }
+    rkmu = sum3;
+    rk1 = sum1 * xi2;
+  } else {
+    const a1 = 0.25;
+    let b = 2 * (1 + x);
+    let d = 1 / b;
+    let h = d;
+    let delh = d;
+    let q1 = 0;
+    let q2 = 1;
+    let q = a1;
+    let cc = a1;
+    let a = -a1;
+    let s = 1 + q * delh;
+    for (let i = 2; i <= 1e4; i++) {
+      a -= 2 * (i - 1);
+      cc = -a * cc / i;
+      const qnew = (q1 - b * q2) / a;
+      q1 = q2;
+      q2 = qnew;
+      q += cc * qnew;
+      b += 2;
+      d = 1 / (b + a * d);
+      delh = (b * d - 1) * delh;
+      h += delh;
+      const dels = q * delh;
+      s += dels;
+      if (Math.abs(dels / s) < EPS4) break;
+    }
+    h = a1 * h;
+    rkmu = Math.sqrt(Math.PI / (2 * x)) * Math.exp(-x) / s;
+    rk1 = rkmu * (x + 0.5 - h) * xi;
+  }
+  return [rkmu, rk1];
 }
 function besselKScalar(n, x) {
   const ni = Math.round(Math.abs(n));
   if (x <= 0) return NaN;
-  const k0 = x <= 8 ? besselK0Series(x) : besselKAsym(0, x);
-  const k1 = x <= 8 ? besselK1Series(x) : besselKAsym(1, x);
+  const k01 = besselK01(x);
+  const k0 = k01[0];
+  const k1 = k01[1];
   if (ni === 0) return k0;
   if (ni === 1) return k1;
   let kPrev = k0;
@@ -7349,10 +7362,7 @@ var besselK = mathTyped10("besselK", {
   "number, Float64Array": (n, x) => mapArray(
     x,
     (v) => besselKScalar(n, v),
-    () => kernelSource(
-      [besselKAsym, besselK0Series, besselK1Series, besselKScalar],
-      `(x) => besselKScalar(${n}, x)`
-    )
+    () => kernelSource([besselK01, besselKScalar], `(x) => besselKScalar(${n}, x)`)
   )
 });
 var ellipticK = mathTyped10("ellipticK", {
@@ -36212,6 +36222,21 @@ var createVariance = /* @__PURE__ */ factory(
 );
 
 // src/statistics/quantileSeq.ts
+var QUANTILE_MODES = [
+  "linear",
+  "lower",
+  "higher",
+  "nearest",
+  "midpoint"
+];
+function _assertMode(mode2) {
+  if (!QUANTILE_MODES.includes(mode2)) {
+    throw new Error(
+      `Unknown quantile interpolation mode "${mode2}". Expected one of: ${QUANTILE_MODES.join(", ")}`
+    );
+  }
+  return mode2;
+}
 var name222 = "quantileSeq";
 var dependencies222 = [
   "typed",
@@ -36249,24 +36274,29 @@ var createQuantileSeq = /* @__PURE__ */ factory(
     return typed3(name222, {
       "Array | Matrix, number | BigNumber": (data, p) => _quantileSeqProbNumber(data, p, false),
       "Array | Matrix, number | BigNumber, number": (data, prob, dim) => _quantileSeqDim(data, prob, false, dim, _quantileSeqProbNumber),
+      // Third arg as a string selects the interpolation mode (linear default).
+      "Array | Matrix, number | BigNumber, string": (data, prob, mode2) => _quantileSeqProbNumber(data, prob, false, _assertMode(mode2)),
       "Array | Matrix, number | BigNumber, boolean": _quantileSeqProbNumber,
       "Array | Matrix, number | BigNumber, boolean, number": (data, prob, sorted, dim) => _quantileSeqDim(data, prob, sorted, dim, _quantileSeqProbNumber),
+      "Array | Matrix, number | BigNumber, boolean, string": (data, prob, sorted, mode2) => _quantileSeqProbNumber(data, prob, sorted, _assertMode(mode2)),
       "Array | Matrix, Array | Matrix": (data, p) => _quantileSeqProbCollection(data, p, false),
       "Array | Matrix, Array | Matrix, number": (data, prob, dim) => _quantileSeqDim(data, prob, false, dim, _quantileSeqProbCollection),
+      "Array | Matrix, Array | Matrix, string": (data, p, mode2) => _quantileSeqProbCollection(data, p, false, _assertMode(mode2)),
       "Array | Matrix, Array | Matrix, boolean": _quantileSeqProbCollection,
-      "Array | Matrix, Array | Matrix, boolean, number": (data, prob, sorted, dim) => _quantileSeqDim(data, prob, sorted, dim, _quantileSeqProbCollection)
+      "Array | Matrix, Array | Matrix, boolean, number": (data, prob, sorted, dim) => _quantileSeqDim(data, prob, sorted, dim, _quantileSeqProbCollection),
+      "Array | Matrix, Array | Matrix, boolean, string": (data, p, sorted, mode2) => _quantileSeqProbCollection(data, p, sorted, _assertMode(mode2))
     });
     function _quantileSeqDim(data, prob, sorted, dim, fn) {
       return mapSlices2(data, dim, (x) => fn(x, prob, sorted));
     }
-    function _quantileSeqProbNumber(data, probOrN, sorted) {
+    function _quantileSeqProbNumber(data, probOrN, sorted, mode2 = "linear") {
       let probArr;
       const dataArr = data.valueOf();
       if (smaller2(probOrN, 0)) {
         throw new Error("N/prob must be non-negative");
       }
       if (smallerEq2(probOrN, 1)) {
-        return isNumber(probOrN) ? _quantileSeq2(dataArr, probOrN, sorted) : bignumber2(_quantileSeq2(dataArr, probOrN, sorted));
+        return isNumber(probOrN) ? _quantileSeq2(dataArr, probOrN, sorted, mode2) : bignumber2(_quantileSeq2(dataArr, probOrN, sorted, mode2));
       }
       if (larger2(probOrN, 1)) {
         if (!isInteger3(probOrN)) {
@@ -36281,22 +36311,22 @@ var createQuantileSeq = /* @__PURE__ */ factory(
         probArr = [];
         for (let i = 0; smaller2(i, probOrN); i++) {
           const prob = divide2(i + 1, nPlusOne);
-          probArr.push(_quantileSeq2(dataArr, prob, sorted));
+          probArr.push(_quantileSeq2(dataArr, prob, sorted, mode2));
         }
         return isNumber(probOrN) ? probArr : bignumber2(probArr);
       }
       return void 0;
     }
-    function _quantileSeqProbCollection(data, probOrN, sorted) {
+    function _quantileSeqProbCollection(data, probOrN, sorted, mode2 = "linear") {
       const dataArr = data.valueOf();
       const probOrNArr = probOrN.valueOf();
       const probArr = [];
       for (let i = 0; i < probOrNArr.length; ++i) {
-        probArr.push(_quantileSeq2(dataArr, probOrNArr[i], sorted));
+        probArr.push(_quantileSeq2(dataArr, probOrNArr[i], sorted, mode2));
       }
       return probArr;
     }
-    function _quantileSeq2(array, prob, sorted) {
+    function _quantileSeq2(array, prob, sorted, mode2 = "linear") {
       const flat = flatten2(array);
       const len = flat.length;
       if (len === 0) {
@@ -36322,11 +36352,28 @@ var createQuantileSeq = /* @__PURE__ */ factory(
           }
         }
       }
-      const fracPartConverted = isBigNumber(left) && isNumber(fracPart) ? bignumber2(fracPart) : fracPart;
-      return add3(
-        multiply2(left, subtract3(1, fracPartConverted)),
-        multiply2(right, fracPartConverted)
-      );
+      switch (mode2) {
+        case "lower":
+          return left;
+        case "higher":
+          return right;
+        case "nearest": {
+          const c = compare2(fracPart, 0.5);
+          if (c < 0) return left;
+          if (c > 0) return right;
+          return integerPart % 2 === 0 ? left : right;
+        }
+        case "midpoint":
+          return divide2(add3(left, right), 2);
+        case "linear":
+        default: {
+          const fracPartConverted = isBigNumber(left) && isNumber(fracPart) ? bignumber2(fracPart) : fracPart;
+          return add3(
+            multiply2(left, subtract3(1, fracPartConverted)),
+            multiply2(right, fracPartConverted)
+          );
+        }
+      }
     }
   }
 );

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

@@ -1 +1 @@
-{"version":3,"file":"quantileSeq.d.ts","sourceRoot":"","sources":["../../src/statistics/quantileSeq.ts"],"names":[],"mappings":"AAGA,OAAO,KAAK,EAAE,aAAa,EAAE,MAAM,2BAA2B,CAAC;AAO/D,UAAU,aAAa;IACrB,KAAK,CAAC,CAAC,EAAE,MAAM,GAAG,aAAa,CAAC;IAChC,KAAK,IAAI,aAAa,CAAC;IACvB,QAAQ,IAAI,MAAM,CAAC;IACnB,GAAG,CAAC,CAAC,EAAE,MAAM,GAAG,aAAa,GAAG,aAAa,CAAC;IAC9C,OAAO,IAAI,MAAM,CAAC;CACnB;AAED,UAAU,uBAAuB;IAC/B,KAAK,EAAE,aAAa,CAAC;IACrB,SAAS,CAAC,EAAE,CAAC,KAAK,EAAE,OAAO,KAAK,aAAa,CAAC;IAC9C,GAAG,EAAE,aAAa,CAAC;IACnB,QAAQ,EAAE,aAAa,CAAC;IACxB,MAAM,EAAE,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,KAAK,MAAM,GAAG,aAAa,CAAC;IAC3D,QAAQ,EAAE,aAAa,CAAC;IACxB,eAAe,EAAE,aAAa,CAAC;IAC/B,OAAO,EAAE,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,KAAK,MAAM,CAAC;IAC5C,SAAS,EAAE,aAAa,CAAC;IACzB,OAAO,EAAE,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,KAAK,OAAO,CAAC;IAC7C,SAAS,EAAE,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,KAAK,OAAO,CAAC;IAC/C,MAAM,EAAE,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,KAAK,OAAO,CAAC;IAC5C,SAAS,EAAE,aAAa,CAAC;CAC1B;AAmBD,eAAO,MAAM,iBAAiB,iFAkO7B,CAAC"}
+{"version":3,"file":"quantileSeq.d.ts","sourceRoot":"","sources":["../../src/statistics/quantileSeq.ts"],"names":[],"mappings":"AAGA,OAAO,KAAK,EAAE,aAAa,EAAE,MAAM,2BAA2B,CAAC;AAO/D,UAAU,aAAa;IACrB,KAAK,CAAC,CAAC,EAAE,MAAM,GAAG,aAAa,CAAC;IAChC,KAAK,IAAI,aAAa,CAAC;IACvB,QAAQ,IAAI,MAAM,CAAC;IACnB,GAAG,CAAC,CAAC,EAAE,MAAM,GAAG,aAAa,GAAG,aAAa,CAAC;IAC9C,OAAO,IAAI,MAAM,CAAC;CACnB;AAED,UAAU,uBAAuB;IAC/B,KAAK,EAAE,aAAa,CAAC;IACrB,SAAS,CAAC,EAAE,CAAC,KAAK,EAAE,OAAO,KAAK,aAAa,CAAC;IAC9C,GAAG,EAAE,aAAa,CAAC;IACnB,QAAQ,EAAE,aAAa,CAAC;IACxB,MAAM,EAAE,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,KAAK,MAAM,GAAG,aAAa,CAAC;IAC3D,QAAQ,EAAE,aAAa,CAAC;IACxB,eAAe,EAAE,aAAa,CAAC;IAC/B,OAAO,EAAE,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,KAAK,MAAM,CAAC;IAC5C,SAAS,EAAE,aAAa,CAAC;IACzB,OAAO,EAAE,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,KAAK,OAAO,CAAC;IAC7C,SAAS,EAAE,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,KAAK,OAAO,CAAC;IAC/C,MAAM,EAAE,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,KAAK,OAAO,CAAC;IAC5C,SAAS,EAAE,aAAa,CAAC;CAC1B;AA0CD,eAAO,MAAM,iBAAiB,iFAiR7B,CAAC"}

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

@@ -1 +1 @@
-{"version":3,"file":"special.d.ts","sourceRoot":"","sources":["../../src/typed/special.ts"],"names":[],"mappings":"AAEA;;;;;;;;;;;;;GAaG;AAwEH,0CAA0C;AAC1C,KAAK,GAAG,GAAG,MAAM,CAAC;AA0DlB,yDAAyD;AACzD,wBAAgB,UAAU,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAKtC;AA6ZD;;;;;;;;;GASG;AACH,eAAO,MAAM,IAAI,wCASf,CAAC;AAEH;;;;;GAKG;AACH,eAAO,MAAM,IAAI,wCAIf,CAAC;AAqDH,eAAO,MAAM,MAAM,wCASjB,CAAC;AAEH;;;;;;;;;GASG;AACH,eAAO,MAAM,IAAI,wCAQf,CAAC;AAEH;;;;;;;;;GASG;AACH,eAAO,MAAM,QAAQ,wCAQnB,CAAC;AAEH;;;;;;GAMG;AACH,eAAO,MAAM,SAAS,wCASpB,CAAC;AAEH;;;;;;;GAOG;AACH,eAAO,MAAM,OAAO,wCAQlB,CAAC;AAEH;;;;;;;;GAQG;AACH,eAAO,MAAM,OAAO,wCAIlB,CAAC;AAMH;;;;;;;;;GASG;AACH,eAAO,MAAM,QAAQ,wCAUnB,CAAC;AAEH;;;;;;;;GAQG;AACH,eAAO,MAAM,QAAQ,wCAUnB,CAAC;AAEH;;;;;;;;GAQG;AACH,eAAO,MAAM,QAAQ,wCAanB,CAAC;AAEH;;;;;;;;GAQG;AACH,eAAO,MAAM,QAAQ,wCAanB,CAAC;AAEH;;;;;;;;;GASG;AACH,eAAO,MAAM,OAAO,wCAuBlB,CAAC;AAEH;;;;;;;;;GASG;AACH,eAAO,MAAM,OAAO,wCA2BlB,CAAC;AAEH;;;;;;GAMG;AACH,eAAO,MAAM,OAAO,wCAQlB,CAAC;AAEH;;;;;;GAMG;AACH,eAAO,MAAM,OAAO,wCAYlB,CAAC;AAMH;;;;;;;;;;;;;;;;;;GAkBG;AACH,eAAO,MAAM,SAAS,wCAUpB,CAAC;AAEH;;;;;;;;;;;;;;;;;;;GAmBG;AACH,eAAO,MAAM,SAAS,wCAWpB,CAAC;AAMH;;;;;;GAMG;AACH,eAAO,MAAM,UAAU,wCAQrB,CAAC;AAEH;;;;;;GAMG;AACH,eAAO,MAAM,QAAQ,wCAQnB,CAAC;AAEH;;;;;;GAMG;AACH,eAAO,MAAM,SAAS,wCAQpB,CAAC;AAEH;;;;;;GAMG;AACH,eAAO,MAAM,SAAS,wCAQpB,CAAC;AAMH;;;;;;;;;;;GAWG;AACH,eAAO,MAAM,QAAQ,wCASnB,CAAC;AAMH;;;;;GAKG;AACH,eAAO,MAAM,WAAW,wCAMtB,CAAC;AAEH;;;;;GAKG;AACH,eAAO,MAAM,WAAW,wCAMtB,CAAC;AAEH;;;;;GAKG;AACH,eAAO,MAAM,WAAW,wCAMtB,CAAC;AAEH;;;;;GAKG;AACH,eAAO,MAAM,aAAa,wCAMxB,CAAC;AAMH;;;;;GAKG;AACH,eAAO,MAAM,QAAQ,wCAInB,CAAC;AAEH;;;;;GAKG;AACH,eAAO,MAAM,QAAQ,wCAInB,CAAC;AAMH;;;;;;;;;;;;;;;;;;GAkBG;AACH,eAAO,MAAM,MAAM,wCAUjB,CAAC;AAEH;;;;;;;;;;;;;GAaG;AACH,eAAO,MAAM,MAAM,wCAUjB,CAAC;AAUH;;;;;;;;;;;GAWG;AACH,eAAO,MAAM,SAAS,wCAMpB,CAAC;AAEH;;;;;;;;;GASG;AACH,eAAO,MAAM,SAAS,wCAUpB,CAAC;AAEH;;;;;;;;GAQG;AACH,eAAO,MAAM,SAAS,wCAUpB,CAAC;AAEH;;;;;;GAMG;AACH,eAAO,MAAM,SAAS,wCAWpB,CAAC;AAYH;;;;;;;;;;;;;GAaG;AACH,eAAO,MAAM,SAAS,wCAMpB,CAAC;AAEH;;;;;;;;;;;;;GAaG;AACH,eAAO,MAAM,mBAAmB,wCAM9B,CAAC;AAEH;;;;;;;;;;;GAWG;AACH,eAAO,MAAM,UAAU,wCAUrB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,YAAY;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;CAuCxB,CAAC"}
+{"version":3,"file":"special.d.ts","sourceRoot":"","sources":["../../src/typed/special.ts"],"names":[],"mappings":"AAEA;;;;;;;;;;;;;GAaG;AAwEH,0CAA0C;AAC1C,KAAK,GAAG,GAAG,MAAM,CAAC;AA0DlB,yDAAyD;AACzD,wBAAgB,UAAU,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAKtC;AAuaD;;;;;;;;;GASG;AACH,eAAO,MAAM,IAAI,wCASf,CAAC;AAEH;;;;;GAKG;AACH,eAAO,MAAM,IAAI,wCAIf,CAAC;AAqDH,eAAO,MAAM,MAAM,wCASjB,CAAC;AAEH;;;;;;;;;GASG;AACH,eAAO,MAAM,IAAI,wCAQf,CAAC;AAEH;;;;;;;;;GASG;AACH,eAAO,MAAM,QAAQ,wCAQnB,CAAC;AAEH;;;;;;GAMG;AACH,eAAO,MAAM,SAAS,wCASpB,CAAC;AAEH;;;;;;;GAOG;AACH,eAAO,MAAM,OAAO,wCAQlB,CAAC;AAEH;;;;;;;;GAQG;AACH,eAAO,MAAM,OAAO,wCAIlB,CAAC;AAMH;;;;;;;;;GASG;AACH,eAAO,MAAM,QAAQ,wCAUnB,CAAC;AAEH;;;;;;;;GAQG;AACH,eAAO,MAAM,QAAQ,wCAUnB,CAAC;AAEH;;;;;;;;GAQG;AACH,eAAO,MAAM,QAAQ,wCAanB,CAAC;AAEH;;;;;;;;GAQG;AACH,eAAO,MAAM,QAAQ,wCAanB,CAAC;AAEH;;;;;;;;;GASG;AACH,eAAO,MAAM,OAAO,wCAuBlB,CAAC;AAEH;;;;;;;;;GASG;AACH,eAAO,MAAM,OAAO,wCA2BlB,CAAC;AAEH;;;;;;GAMG;AACH,eAAO,MAAM,OAAO,wCAQlB,CAAC;AAEH;;;;;;GAMG;AACH,eAAO,MAAM,OAAO,wCAQlB,CAAC;AAMH;;;;;;;;;;;;;;;;;;GAkBG;AACH,eAAO,MAAM,SAAS,wCAUpB,CAAC;AAEH;;;;;;;;;;;;;;;;;;;GAmBG;AACH,eAAO,MAAM,SAAS,wCAWpB,CAAC;AAMH;;;;;;GAMG;AACH,eAAO,MAAM,UAAU,wCAQrB,CAAC;AAEH;;;;;;GAMG;AACH,eAAO,MAAM,QAAQ,wCAQnB,CAAC;AAEH;;;;;;GAMG;AACH,eAAO,MAAM,SAAS,wCAQpB,CAAC;AAEH;;;;;;GAMG;AACH,eAAO,MAAM,SAAS,wCAQpB,CAAC;AAMH;;;;;;;;;;;GAWG;AACH,eAAO,MAAM,QAAQ,wCASnB,CAAC;AAMH;;;;;GAKG;AACH,eAAO,MAAM,WAAW,wCAMtB,CAAC;AAEH;;;;;GAKG;AACH,eAAO,MAAM,WAAW,wCAMtB,CAAC;AAEH;;;;;GAKG;AACH,eAAO,MAAM,WAAW,wCAMtB,CAAC;AAEH;;;;;GAKG;AACH,eAAO,MAAM,aAAa,wCAMxB,CAAC;AAMH;;;;;GAKG;AACH,eAAO,MAAM,QAAQ,wCAInB,CAAC;AAEH;;;;;GAKG;AACH,eAAO,MAAM,QAAQ,wCAInB,CAAC;AAMH;;;;;;;;;;;;;;;;;;GAkBG;AACH,eAAO,MAAM,MAAM,wCAUjB,CAAC;AAEH;;;;;;;;;;;;;GAaG;AACH,eAAO,MAAM,MAAM,wCAUjB,CAAC;AAUH;;;;;;;;;;;GAWG;AACH,eAAO,MAAM,SAAS,wCAMpB,CAAC;AAEH;;;;;;;;;GASG;AACH,eAAO,MAAM,SAAS,wCAUpB,CAAC;AAEH;;;;;;;;GAQG;AACH,eAAO,MAAM,SAAS,wCAUpB,CAAC;AAEH;;;;;;GAMG;AACH,eAAO,MAAM,SAAS,wCAWpB,CAAC;AAYH;;;;;;;;;;;;;GAaG;AACH,eAAO,MAAM,SAAS,wCAMpB,CAAC;AAEH;;;;;;;;;;;;;GAaG;AACH,eAAO,MAAM,mBAAmB,wCAM9B,CAAC;AAEH;;;;;;;;;;;GAWG;AACH,eAAO,MAAM,UAAU,wCAUrB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,YAAY;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;CAuCxB,CAAC"}

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@danielsimonjr/mathts-functions",
-  "version": "0.37.0",
+  "version": "0.38.0",
   "description": "Mathematical functions for MathTS - arithmetic, algebra, trigonometry, statistics, and more",
   "author": "Daniel Simon Jr.",
   "license": "MIT",
@@ -37,7 +37,7 @@
     "@danielsimonjr/mathts-core": "^0.10.0",
     "@danielsimonjr/mathts-expression": "^0.6.4",
     "@danielsimonjr/mathts-gpu": "^0.2.0",
-    "@danielsimonjr/mathts-matrix": "^0.4.6",
+    "@danielsimonjr/mathts-matrix": "^0.5.0",
     "@danielsimonjr/mathts-parallel": "^0.6.0",
     "bignumber.js": "^9.1.2",
     "complex.js": "^2.2.5",