@danielsimonjr/mathts-functions 0.2.4 → 0.2.5

package/README.md

@@ -0,0 +1,21 @@
+# @danielsimonjr/mathts-functions
+
+Math functions for [MathTS](https://github.com/danielsimonjr/mathts).
+
+The mathematical function library for MathTS — arithmetic, trigonometry, algebra, statistics, signal processing, and more — dispatched polymorphically via typed-function across number / Complex / Fraction / BigNumber / Matrix.
+
+## Install
+
+```sh
+npm install @danielsimonjr/mathts-functions
+```
+
+## What it provides
+
+- Hundreds of typed functions (`add`, `sin`, `svd`-aware ops, `parallelFFT`, …).
+- A wired `parse` / `evaluate` against the full activated math scope.
+- Per-domain focused entry points are also published: `@danielsimonjr/mathts-{arithmetic,trigonometry,statistics,signal}`.
+
+## License
+
+MIT (c) Daniel Simon Jr.

package/dist/index.js

@@ -5463,8 +5463,8 @@ function besselY1Scalar(x) {
   if (x <= 0) return NaN;
   if (x < 8) {
     const y = x * x;
-    const r1 = -4900604943e3 + y * (127527439e4 + y * (-51534866840 + y * (6227854327e-1 + y * (-3130827714e-3 + y * 7.374753505))));
-    const r2 = 249958057e5 + y * (424441966400 + y * (3733650367 + y * (2245976615e-2 + y * (103832.3184 + y * (365.2510261 + y * 1)))));
+    const r1 = x * (-4900604943e3 + y * (127527439e4 + y * (-51534381390 + y * (7349264551e-1 + y * (-4237922726e-3 + y * 8511.937935)))));
+    const r2 = 249958057e5 + y * (424441966400 + y * (3733650367 + y * (2245904002e-2 + y * (102042.605 + y * (354.9632885 + y)))));
     return r1 / r2 + 0.636619772 * (besselJ1Scalar(x) * Math.log(x) - 1 / x);
   } else {
     const z = 8 / x;
@@ -5732,26 +5732,40 @@ function lambertWScalar(x) {
 function cosIntegralScalar(x) {
   if (x <= 0) return NaN;
   const euler = 0.5772156649015329;
-  if (x <= 4) {
+  if (x <= 35) {
     let sum2 = 0;
     let term = 1;
-    for (let n = 1; n < 100; n++) {
+    for (let n = 1; n < 200; n++) {
       term *= -x * x / ((2 * n - 1) * (2 * n));
       sum2 += term / (2 * n);
-      if (Math.abs(term / (2 * n)) < 1e-15) break;
+      if (Math.abs(term / (2 * n)) < Math.abs(sum2) * 1e-16) break;
     }
     return euler + Math.log(x) + sum2;
   }
-  let f = 0, g = 0;
+  let f = 1, g = 1;
   let fn = 1, gn = 1;
-  for (let n = 1; n <= 30; n++) {
-    fn *= -(2 * n - 1) * (2 * n - 2) / (x * x);
-    gn *= -(2 * n) * (2 * n - 1) / (x * x);
-    f += fn;
-    g += gn;
+  let fDone = false, gDone = false;
+  for (let n = 1; n <= 60; n++) {
+    const fnNext = fn * (-(2 * n) * (2 * n - 1)) / (x * x);
+    const gnNext = gn * (-(2 * n + 1) * (2 * n)) / (x * x);
+    if (!fDone) {
+      if (Math.abs(fnNext) > Math.abs(fn)) fDone = true;
+      else {
+        fn = fnNext;
+        f += fn;
+      }
+    }
+    if (!gDone) {
+      if (Math.abs(gnNext) > Math.abs(gn)) gDone = true;
+      else {
+        gn = gnNext;
+        g += gn;
+      }
+    }
+    if (fDone && gDone) break;
   }
-  f = (1 + f) / x;
-  g = (1 + g) / x;
+  f = f / x;
+  g = g / (x * x);
   return f * Math.sin(x) - g * Math.cos(x);
 }
 function sinIntegralScalar(x) {
@@ -5807,6 +5821,33 @@ function fresnelCScalar(x) {
   }
   return sign3 * sum2;
 }
+var AIRY_U = (() => {
+  const u = [1];
+  for (let k = 1; k <= 10; k++) {
+    u.push(u[k - 1] * ((6 * k - 5) * (6 * k - 3) * (6 * k - 1)) / (216 * k));
+  }
+  return u;
+})();
+function airyAsymPQ(zeta2) {
+  let p = 0, q = 0;
+  let pTerm = Infinity, qTerm = Infinity;
+  const half = Math.floor((AIRY_U.length - 1) / 2);
+  for (let k = 0; k <= half; k++) {
+    const sign3 = k % 2 === 0 ? 1 : -1;
+    const pt = sign3 * AIRY_U[2 * k] / Math.pow(zeta2, 2 * k);
+    if (Math.abs(pt) > Math.abs(pTerm)) break;
+    p += pt;
+    pTerm = pt;
+  }
+  for (let k = 0; k <= half - 1; k++) {
+    const sign3 = k % 2 === 0 ? 1 : -1;
+    const qt = sign3 * AIRY_U[2 * k + 1] / Math.pow(zeta2, 2 * k + 1);
+    if (Math.abs(qt) > Math.abs(qTerm)) break;
+    q += qt;
+    qTerm = qt;
+  }
+  return { p, q };
+}
 function airyAiScalar(x) {
   const XBIG = 4.5;
   const AI0 = 0.3550280538878172;
@@ -5836,14 +5877,8 @@ function airyAiScalar(x) {
     const axp = Math.pow(ax, 0.25);
     const zeta2 = 2 / 3 * ax * Math.sqrt(ax);
     const theta = zeta2 - Math.PI / 4;
-    let p = 0, q = 0, zk = 1, sign3 = 1;
-    for (let k = 0; k < C.length; k++) {
-      if (k % 2 === 0) p += sign3 * C[k] / zk;
-      else q += sign3 * C[k] / zk;
-      zk *= zeta2;
-      sign3 = -sign3;
-    }
-    return (Math.sin(theta) * p + Math.cos(theta) * q) / (Math.sqrt(Math.PI) * axp);
+    const { p, q } = airyAsymPQ(zeta2);
+    return (Math.cos(theta) * p + Math.sin(theta) * q) / (Math.sqrt(Math.PI) * axp);
   }
   const x3 = x * x * x;
   let f = 1, g = x, fTerm = 1, gTerm = x;
@@ -5889,15 +5924,9 @@ function airyBiScalar(x) {
     const ax = Math.abs(x);
     const axp = Math.pow(ax, 0.25);
     const zeta2 = 2 / 3 * ax * Math.sqrt(ax);
-    const theta = zeta2 + Math.PI / 4;
-    let p = 0, q = 0, zk = 1, sign3 = 1;
-    for (let k = 0; k < C.length; k++) {
-      if (k % 2 === 0) p += sign3 * C[k] / zk;
-      else q += sign3 * C[k] / zk;
-      zk *= zeta2;
-      sign3 = -sign3;
-    }
-    return (Math.cos(theta) * p + Math.sin(theta) * q) / (Math.sqrt(Math.PI) * axp);
+    const theta = zeta2 - Math.PI / 4;
+    const { p, q } = airyAsymPQ(zeta2);
+    return (-Math.sin(theta) * p + Math.cos(theta) * q) / (Math.sqrt(Math.PI) * axp);
   }
   const x3 = x * x * x;
   let f = 1, g = x, fTerm = 1, gTerm = x;

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;AAk9BH;;;;;;;;;GASG;AACH,eAAO,MAAM,IAAI,wCAIf,CAAC;AAEH;;;;;GAKG;AACH,eAAO,MAAM,IAAI,wCAIf,CAAC;AAMH;;;;;;;;;;;;;;;;;;GAkBG;AACH,eAAO,MAAM,MAAM,wCAQjB,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,wCAUnB,CAAC;AAEH;;;;;;;;GAQG;AACH,eAAO,MAAM,QAAQ,wCAUnB,CAAC;AAEH;;;;;;;;;GASG;AACH,eAAO,MAAM,OAAO,wCAgBlB,CAAC;AAEH;;;;;;;;;GASG;AACH,eAAO,MAAM,OAAO,wCAgBlB,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;;;;;;;;GAQG;AACH,eAAO,MAAM,QAAQ,wCAInB,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,wCAQjB,CAAC;AAEH;;;;;;;;;;;;;GAaG;AACH,eAAO,MAAM,MAAM,wCAQjB,CAAC;AAUH;;;;;;;;;;;GAWG;AACH,eAAO,MAAM,SAAS,wCAMpB,CAAC;AAEH;;;;;;;;;GASG;AACH,eAAO,MAAM,SAAS,wCAMpB,CAAC;AAEH;;;;;;;;GAQG;AACH,eAAO,MAAM,SAAS,wCAMpB,CAAC;AAEH;;;;;;GAMG;AACH,eAAO,MAAM,SAAS,wCAMpB,CAAC;AAYH;;;;;;;;;;;;;GAaG;AACH,eAAO,MAAM,SAAS,wCAMpB,CAAC;AAEH;;;;;;;;;;;;;GAaG;AACH,eAAO,MAAM,mBAAmB,wCAM9B,CAAC;AAEH;;;;;;;;;;;GAWG;AACH,eAAO,MAAM,UAAU,wCAMrB,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;AAohCH;;;;;;;;;GASG;AACH,eAAO,MAAM,IAAI,wCAIf,CAAC;AAEH;;;;;GAKG;AACH,eAAO,MAAM,IAAI,wCAIf,CAAC;AAMH;;;;;;;;;;;;;;;;;;GAkBG;AACH,eAAO,MAAM,MAAM,wCAQjB,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,wCAUnB,CAAC;AAEH;;;;;;;;GAQG;AACH,eAAO,MAAM,QAAQ,wCAUnB,CAAC;AAEH;;;;;;;;;GASG;AACH,eAAO,MAAM,OAAO,wCAgBlB,CAAC;AAEH;;;;;;;;;GASG;AACH,eAAO,MAAM,OAAO,wCAgBlB,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;;;;;;;;GAQG;AACH,eAAO,MAAM,QAAQ,wCAInB,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,wCAQjB,CAAC;AAEH;;;;;;;;;;;;;GAaG;AACH,eAAO,MAAM,MAAM,wCAQjB,CAAC;AAUH;;;;;;;;;;;GAWG;AACH,eAAO,MAAM,SAAS,wCAMpB,CAAC;AAEH;;;;;;;;;GASG;AACH,eAAO,MAAM,SAAS,wCAMpB,CAAC;AAEH;;;;;;;;GAQG;AACH,eAAO,MAAM,SAAS,wCAMpB,CAAC;AAEH;;;;;;GAMG;AACH,eAAO,MAAM,SAAS,wCAMpB,CAAC;AAYH;;;;;;;;;;;;;GAaG;AACH,eAAO,MAAM,SAAS,wCAMpB,CAAC;AAEH;;;;;;;;;;;;;GAaG;AACH,eAAO,MAAM,mBAAmB,wCAM9B,CAAC;AAEH;;;;;;;;;;;GAWG;AACH,eAAO,MAAM,UAAU,wCAMrB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,YAAY;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;CAuCxB,CAAC"}

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@danielsimonjr/mathts-functions",
-  "version": "0.2.4",
+  "version": "0.2.5",
   "description": "Mathematical functions for MathTS - arithmetic, algebra, trigonometry, statistics, and more",
   "author": "Daniel Simon Jr.",
   "license": "MIT",