@danielsimonjr/mathts-functions 0.68.2 → 0.68.4

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
@@ -9,7 +9,7 @@ interface SimplifyOptions {
9
9
  }
10
10
  interface SymbolicEqualDependencies {
11
11
  parse: (expr: string) => MathNode;
12
- simplify: (node: MathNode, rules: unknown[], options: SimplifyOptions) => MathNode;
12
+ simplify: (node: MathNode, scope: Map<string, unknown>, options: SimplifyOptions) => MathNode;
13
13
  typed: TypedFunction;
14
14
  OperatorNode: OperatorNodeConstructor;
15
15
  }
@@ -1 +1 @@
1
- {"version":3,"file":"symbolicEqual.d.ts","sourceRoot":"","sources":["../../src/algebra/symbolicEqual.ts"],"names":[],"mappings":"AAEA,OAAO,KAAK,EAAE,QAAQ,EAAE,MAAM,kBAAkB,CAAC;AACjD,OAAO,KAAK,EAAE,aAAa,EAAE,MAAM,2BAA2B,CAAC;AAO/D,UAAU,uBAAuB;IAC/B,KAAK,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,IAAI,EAAE,QAAQ,EAAE,GAAG,QAAQ,CAAC;CAC1D;AAED,UAAU,eAAe;IACvB,OAAO,CAAC,EAAE,MAAM,CAAC,MAAM,EAAE,OAAO,CAAC,CAAC;IAClC,CAAC,GAAG,EAAE,MAAM,GAAG,OAAO,CAAC;CACxB;AAED,UAAU,yBAAyB;IACjC,KAAK,EAAE,CAAC,IAAI,EAAE,MAAM,KAAK,QAAQ,CAAC;IAClC,QAAQ,EAAE,CAAC,IAAI,EAAE,QAAQ,EAAE,KAAK,EAAE,OAAO,EAAE,EAAE,OAAO,EAAE,eAAe,KAAK,QAAQ,CAAC;IACnF,KAAK,EAAE,aAAa,CAAC;IACrB,YAAY,EAAE,uBAAuB,CAAC;CACvC;AAKD,eAAO,MAAM,mBAAmB,yFAwD/B,CAAC"}
1
+ {"version":3,"file":"symbolicEqual.d.ts","sourceRoot":"","sources":["../../src/algebra/symbolicEqual.ts"],"names":[],"mappings":"AAEA,OAAO,KAAK,EAAE,QAAQ,EAAE,MAAM,kBAAkB,CAAC;AACjD,OAAO,KAAK,EAAE,aAAa,EAAE,MAAM,2BAA2B,CAAC;AAO/D,UAAU,uBAAuB;IAC/B,KAAK,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,IAAI,EAAE,QAAQ,EAAE,GAAG,QAAQ,CAAC;CAC1D;AAED,UAAU,eAAe;IACvB,OAAO,CAAC,EAAE,MAAM,CAAC,MAAM,EAAE,OAAO,CAAC,CAAC;IAClC,CAAC,GAAG,EAAE,MAAM,GAAG,OAAO,CAAC;CACxB;AAED,UAAU,yBAAyB;IACjC,KAAK,EAAE,CAAC,IAAI,EAAE,MAAM,KAAK,QAAQ,CAAC;IAClC,QAAQ,EAAE,CAAC,IAAI,EAAE,QAAQ,EAAE,KAAK,EAAE,GAAG,CAAC,MAAM,EAAE,OAAO,CAAC,EAAE,OAAO,EAAE,eAAe,KAAK,QAAQ,CAAC;IAC9F,KAAK,EAAE,aAAa,CAAC;IACrB,YAAY,EAAE,uBAAuB,CAAC;CACvC;AAKD,eAAO,MAAM,mBAAmB,yFA0D/B,CAAC"}
@@ -27,4 +27,11 @@ export declare const laplaceQuantile: (p: number, mu?: number, b?: number) => nu
27
27
  export declare const logisticPDF: (x: number, mu?: number, s?: number) => number;
28
28
  export declare const logisticCDF: (x: number, mu?: number, s?: number) => number;
29
29
  export declare const logisticQuantile: (p: number, mu?: number, s?: number) => number;
30
+ /**
31
+ * Two-sided p-value of a Student-t statistic `t` with `df` degrees of freedom,
32
+ * computed DIRECTLY in the tail as the regularized incomplete beta
33
+ * `I_{df/(df+t²)}(df/2, 1/2)`. `2·(1 − studentTCDF(|t|, df))` rounds a far
34
+ * tail to 0, and a steep regression's slope test IS that tail.
35
+ */
36
+ export declare function studentTTwoSidedP(t: number, df: number): number;
30
37
  //# sourceMappingURL=distribution-functions.d.ts.map
@@ -1 +1 @@
1
- {"version":3,"file":"distribution-functions.d.ts","sourceRoot":"","sources":["../src/distribution-functions.ts"],"names":[],"mappings":"AAmBA,2DAA2D;AAC3D,eAAO,MAAM,cAAc,MAAO,MAAM,kCAAsB,MAC/B,CAAC;AAEhC,+DAA+D;AAC/D,eAAO,MAAM,WAAW,MAAO,MAAM,MAAM,MAAM,KAAG,MAA6B,CAAC;AAClF,eAAO,MAAM,gBAAgB,MAAO,MAAM,MAAM,MAAM,KAAG,MAAkC,CAAC;AAE5F,qEAAqE;AACrE,eAAO,MAAM,aAAa,MAAO,MAAM,MAAM,MAAM,KAAG,MAA+B,CAAC;AACtF,eAAO,MAAM,kBAAkB,MAAO,MAAM,MAAM,MAAM,KAAG,MAAoC,CAAC;AAEhG,wEAAwE;AACxE,eAAO,MAAM,IAAI,MAAO,MAAM,MAAM,MAAM,MAAM,MAAM,KAAG,MAAgC,CAAC;AAC1F,eAAO,MAAM,SAAS,MAAO,MAAM,MAAM,MAAM,MAAM,MAAM,KAAG,MAAqC,CAAC;AAEpG,oDAAoD;AACpD,eAAO,MAAM,QAAQ,MAAO,MAAM,SAAS,MAAM,oBAAa,MAAoC,CAAC;AACnG,eAAO,MAAM,aAAa,MAAO,MAAM,SAAS,MAAM,oBAAa,MAClC,CAAC;AAElC,0DAA0D;AAC1D,eAAO,MAAM,OAAO,MAAO,MAAM,KAAK,MAAM,KAAK,MAAM,KAAG,MAA4B,CAAC;AACvF,eAAO,MAAM,YAAY,MAAO,MAAM,KAAK,MAAM,KAAK,MAAM,KAAG,MAAiC,CAAC;AAIjG,4EAA4E;AAC5E,eAAO,MAAM,SAAS,MAAO,MAAM,kCAAsB,MACF,CAAC;AACxD,eAAO,MAAM,SAAS,MAAO,MAAM,kCAAsB,MACZ,CAAC;AAC9C,eAAO,MAAM,cAAc,MAAO,MAAM,kCAAsB,MAClB,CAAC;AAE7C,oFAAoF;AACpF,eAAO,MAAM,UAAU,MAAO,MAAM,8BAAkB,MACX,CAAC;AAC5C,eAAO,MAAM,UAAU,MAAO,MAAM,8BAAkB,MACqB,CAAC;AAC5E,eAAO,MAAM,eAAe,MAAO,MAAM,8BAAkB,MACQ,CAAC;AAEpE,gEAAgE;AAChE,eAAO,MAAM,WAAW,MAAO,MAAM,8BAAkB,MAGtD,CAAC;AACF,eAAO,MAAM,WAAW,MAAO,MAAM,8BAAkB,MAA2C,CAAC;AACnG,eAAO,MAAM,gBAAgB,MAAO,MAAM,8BAAkB,MAC5B,CAAC"}
1
+ {"version":3,"file":"distribution-functions.d.ts","sourceRoot":"","sources":["../src/distribution-functions.ts"],"names":[],"mappings":"AAoBA,2DAA2D;AAC3D,eAAO,MAAM,cAAc,MAAO,MAAM,kCAAsB,MAC/B,CAAC;AAEhC,+DAA+D;AAC/D,eAAO,MAAM,WAAW,MAAO,MAAM,MAAM,MAAM,KAAG,MAA6B,CAAC;AAClF,eAAO,MAAM,gBAAgB,MAAO,MAAM,MAAM,MAAM,KAAG,MAAkC,CAAC;AAE5F,qEAAqE;AACrE,eAAO,MAAM,aAAa,MAAO,MAAM,MAAM,MAAM,KAAG,MAA+B,CAAC;AACtF,eAAO,MAAM,kBAAkB,MAAO,MAAM,MAAM,MAAM,KAAG,MAAoC,CAAC;AAEhG,wEAAwE;AACxE,eAAO,MAAM,IAAI,MAAO,MAAM,MAAM,MAAM,MAAM,MAAM,KAAG,MAAgC,CAAC;AAC1F,eAAO,MAAM,SAAS,MAAO,MAAM,MAAM,MAAM,MAAM,MAAM,KAAG,MAAqC,CAAC;AAEpG,oDAAoD;AACpD,eAAO,MAAM,QAAQ,MAAO,MAAM,SAAS,MAAM,oBAAa,MAAoC,CAAC;AACnG,eAAO,MAAM,aAAa,MAAO,MAAM,SAAS,MAAM,oBAAa,MAClC,CAAC;AAElC,0DAA0D;AAC1D,eAAO,MAAM,OAAO,MAAO,MAAM,KAAK,MAAM,KAAK,MAAM,KAAG,MAA4B,CAAC;AACvF,eAAO,MAAM,YAAY,MAAO,MAAM,KAAK,MAAM,KAAK,MAAM,KAAG,MAAiC,CAAC;AAIjG,4EAA4E;AAC5E,eAAO,MAAM,SAAS,MAAO,MAAM,kCAAsB,MACF,CAAC;AACxD,eAAO,MAAM,SAAS,MAAO,MAAM,kCAAsB,MACZ,CAAC;AAC9C,eAAO,MAAM,cAAc,MAAO,MAAM,kCAAsB,MAClB,CAAC;AAE7C,oFAAoF;AACpF,eAAO,MAAM,UAAU,MAAO,MAAM,8BAAkB,MACX,CAAC;AAC5C,eAAO,MAAM,UAAU,MAAO,MAAM,8BAAkB,MACqB,CAAC;AAC5E,eAAO,MAAM,eAAe,MAAO,MAAM,8BAAkB,MACQ,CAAC;AAEpE,gEAAgE;AAChE,eAAO,MAAM,WAAW,MAAO,MAAM,8BAAkB,MAGtD,CAAC;AACF,eAAO,MAAM,WAAW,MAAO,MAAM,8BAAkB,MAA2C,CAAC;AACnG,eAAO,MAAM,gBAAgB,MAAO,MAAM,8BAAkB,MAC5B,CAAC;AAEjC;;;;;GAKG;AACH,wBAAgB,iBAAiB,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,GAAG,MAAM,CAM/D"}
package/dist/index.js CHANGED
@@ -6363,7 +6363,27 @@ function gammaincScalar(a, x) {
6363
6363
  }
6364
6364
  }
6365
6365
  function gammaincpScalar(a, x) {
6366
- return 1 - gammaincScalar(a, x);
6366
+ if (x < 0) return NaN;
6367
+ if (x === 0) return 1;
6368
+ if (x < a + 1) return 1 - gammaincScalar(a, x);
6369
+ const TINY = 1e-30;
6370
+ const b0 = x + 1 - a;
6371
+ let c = 1 / TINY;
6372
+ let d = 1 / b0;
6373
+ let f = d;
6374
+ for (let n = 1; n < 300; n++) {
6375
+ const an = n * (a - n);
6376
+ const bn = x + 2 * n + 1 - a;
6377
+ d = bn + an * d;
6378
+ if (Math.abs(d) < TINY) d = TINY;
6379
+ c = bn + an / c;
6380
+ if (Math.abs(c) < TINY) c = TINY;
6381
+ d = 1 / d;
6382
+ const delta = d * c;
6383
+ f *= delta;
6384
+ if (Math.abs(delta - 1) < 1e-15) break;
6385
+ }
6386
+ return f * Math.exp(-x + a * Math.log(x) - _lgamma(a));
6367
6387
  }
6368
6388
  function betaincScalar(a, b, x) {
6369
6389
  if (x < 0 || x > 1) return NaN;
@@ -14109,11 +14129,12 @@ function gaussKronrod15(f, a, b) {
14109
14129
  var ABS_FLOOR = 1e-14;
14110
14130
  function quad(f, a, b, opts) {
14111
14131
  const tol = opts?.tol ?? 1e-10;
14132
+ const atol = opts?.atol ?? ABS_FLOOR;
14112
14133
  const maxDepth = opts?.maxDepth ?? 50;
14113
14134
  function panel(a2, b2, depth) {
14114
14135
  const { k, g } = gaussKronrod15(f, a2, b2);
14115
14136
  const panelError = Math.abs(k - g);
14116
- if (depth >= maxDepth || panelError <= Math.max(tol * Math.abs(k), ABS_FLOOR)) {
14137
+ if (depth >= maxDepth || panelError <= Math.max(tol * Math.abs(k), atol)) {
14117
14138
  return { value: k, error: panelError };
14118
14139
  }
14119
14140
  const mid = (a2 + b2) / 2;
@@ -14935,9 +14956,13 @@ function gaussLegendre4(f, y0, tspan, options) {
14935
14956
  // src/typed/numeric.ts
14936
14957
  function findRoot(f, a, b, opts) {
14937
14958
  const tol = opts?.tol ?? 1e-12;
14959
+ const xtol = opts?.xtol ?? tol;
14938
14960
  const maxIter = opts?.maxIter ?? 100;
14939
14961
  let fa = f(a);
14940
14962
  let fb = f(b);
14963
+ if (Number.isNaN(fa) || Number.isNaN(fb)) {
14964
+ throw new Error("findRoot: f(a) or f(b) is not a number (NaN)");
14965
+ }
14941
14966
  if (fa * fb > 0) {
14942
14967
  throw new Error("findRoot: f(a) and f(b) must have opposite signs");
14943
14968
  }
@@ -14951,7 +14976,7 @@ function findRoot(f, a, b, opts) {
14951
14976
  let d = 0;
14952
14977
  let s;
14953
14978
  for (let iter = 0; iter < maxIter; iter++) {
14954
- if (Math.abs(fb) < tol || Math.abs(b - a) < tol) {
14979
+ if (Math.abs(fb) < tol || Math.abs(b - a) < xtol) {
14955
14980
  return b;
14956
14981
  }
14957
14982
  if (fa !== fc && fb !== fc) {
@@ -18376,6 +18401,13 @@ var logisticPDF = (x, mu = 0, s = 1) => {
18376
18401
  };
18377
18402
  var logisticCDF = (x, mu = 0, s = 1) => 1 / (1 + Math.exp(-(x - mu) / s));
18378
18403
  var logisticQuantile = (p, mu = 0, s = 1) => mu + s * Math.log(p / (1 - p));
18404
+ function studentTTwoSidedP(t, df) {
18405
+ if (!(df > 0)) return NaN;
18406
+ if (Number.isNaN(t)) return NaN;
18407
+ if (!Number.isFinite(t)) return 0;
18408
+ if (t === 0) return 1;
18409
+ return betaincScalar(df / 2, 0.5, df / (df + t * t));
18410
+ }
18379
18411
 
18380
18412
  // src/stats/inference-extra.ts
18381
18413
  function chi2Contingency(table, opts = {}) {
@@ -42265,7 +42297,7 @@ var createSymbolicEqual = /* @__PURE__ */ factory(
42265
42297
  ({ parse: _parse2, simplify: simplify2, typed: typed3, OperatorNode }) => {
42266
42298
  function _symbolicEqual(e1, e2, options = {}) {
42267
42299
  const diff2 = new OperatorNode("-", "subtract", [e1, e2]);
42268
- const simplified = simplify2(diff2, [], options);
42300
+ const simplified = simplify2(diff2, /* @__PURE__ */ new Map(), options);
42269
42301
  return isConstantNode(simplified) && !simplified.value;
42270
42302
  }
42271
42303
  return typed3(name248, {
@@ -48728,6 +48760,15 @@ function matVec(a, v) {
48728
48760
  return a.map((row2) => row2.reduce((sum3, val, j) => sum3 + val * v[j], 0));
48729
48761
  }
48730
48762
  function invert(a) {
48763
+ const scale4 = a.map((row2, i) => Math.sqrt(Math.abs(row2[i])));
48764
+ if (scale4.some((s) => !(s > 0) || !Number.isFinite(s))) {
48765
+ throw new Error("ols: X^T X is singular (predictors are collinear)");
48766
+ }
48767
+ const equilibrated = a.map((row2, i) => row2.map((v, j) => v / (scale4[i] * scale4[j])));
48768
+ const inv3 = invertUnscaled(equilibrated);
48769
+ return inv3.map((row2, i) => row2.map((v, j) => v / (scale4[i] * scale4[j])));
48770
+ }
48771
+ function invertUnscaled(a) {
48731
48772
  const n = a.length;
48732
48773
  const M = a.map((row2, i) => [...row2, ...Array.from({ length: n }, (_, j) => i === j ? 1 : 0)]);
48733
48774
  for (let col = 0; col < n; col++) {
@@ -48739,7 +48780,7 @@ function invert(a) {
48739
48780
  pivotRow = row2;
48740
48781
  }
48741
48782
  }
48742
- if (maxVal < 1e-15) throw new Error("ols: X^T X is singular (predictors are collinear)");
48783
+ if (maxVal < 1e-12) throw new Error("ols: X^T X is singular (predictors are collinear)");
48743
48784
  if (pivotRow !== col) [M[col], M[pivotRow]] = [M[pivotRow], M[col]];
48744
48785
  const pivot = M[col][col];
48745
48786
  for (let j = 0; j < 2 * n; j++) M[col][j] /= pivot;
@@ -48776,7 +48817,7 @@ function ols(X, y, opts) {
48776
48817
  const sigma2 = rss / dfResid;
48777
48818
  const stderr = XtXinv.map((row2, j) => Math.sqrt(sigma2 * row2[j]));
48778
48819
  const tValues = coefficients.map((b, j) => b / stderr[j]);
48779
- const pValues = tValues.map((t) => 2 * (1 - studentTCDF(Math.abs(t), dfResid)));
48820
+ const pValues = tValues.map((t) => studentTTwoSidedP(t, dfResid));
48780
48821
  const r2 = tss === 0 ? 1 : 1 - rss / tss;
48781
48822
  const adjR2 = 1 - (1 - r2) * (n - 1) / dfResid;
48782
48823
  const fStat = p > 1 ? (tss - rss) / (p - 1) / (rss / dfResid) : NaN;
@@ -48928,7 +48969,7 @@ function linearRegression(x, y) {
48928
48969
  const df = n - 2;
48929
48970
  const stdErr = Math.sqrt((syy - slope * sxy) / df / sxx);
48930
48971
  const t = slope / stdErr;
48931
- const pValue = 2 * (1 - studentTCDF(Math.abs(t), df));
48972
+ const pValue = studentTTwoSidedP(t, df);
48932
48973
  return { slope, intercept, rValue, rSquared: rValue * rValue, pValue, stdErr };
48933
48974
  }
48934
48975
 
@@ -1 +1 @@
1
- {"version":3,"file":"ols.d.ts","sourceRoot":"","sources":["../../src/ml/ols.ts"],"names":[],"mappings":"AAUA,yBAAyB;AACzB,MAAM,WAAW,UAAU;IACzB,+CAA+C;IAC/C,SAAS,CAAC,EAAE,OAAO,CAAC;CACrB;AAED,wEAAwE;AACxE,MAAM,WAAW,SAAS;IACxB,0DAA0D;IAC1D,YAAY,EAAE,MAAM,EAAE,CAAC;IACvB,0CAA0C;IAC1C,MAAM,EAAE,MAAM,EAAE,CAAC;IACjB,8DAA8D;IAC9D,OAAO,EAAE,MAAM,EAAE,CAAC;IAClB,4DAA4D;IAC5D,OAAO,EAAE,MAAM,EAAE,CAAC;IAClB,oCAAoC;IACpC,EAAE,EAAE,MAAM,CAAC;IACX,gDAAgD;IAChD,KAAK,EAAE,MAAM,CAAC;IACd,kEAAkE;IAClE,KAAK,EAAE,MAAM,CAAC;IACd,wBAAwB;IACxB,SAAS,EAAE,MAAM,EAAE,CAAC;CACrB;AAqED;;;;;;;;;;;;GAYG;AACH,wBAAgB,GAAG,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,IAAI,CAAC,EAAE,UAAU,GAAG,SAAS,CAmC5E"}
1
+ {"version":3,"file":"ols.d.ts","sourceRoot":"","sources":["../../src/ml/ols.ts"],"names":[],"mappings":"AAUA,yBAAyB;AACzB,MAAM,WAAW,UAAU;IACzB,+CAA+C;IAC/C,SAAS,CAAC,EAAE,OAAO,CAAC;CACrB;AAED,wEAAwE;AACxE,MAAM,WAAW,SAAS;IACxB,0DAA0D;IAC1D,YAAY,EAAE,MAAM,EAAE,CAAC;IACvB,0CAA0C;IAC1C,MAAM,EAAE,MAAM,EAAE,CAAC;IACjB,8DAA8D;IAC9D,OAAO,EAAE,MAAM,EAAE,CAAC;IAClB,4DAA4D;IAC5D,OAAO,EAAE,MAAM,EAAE,CAAC;IAClB,oCAAoC;IACpC,EAAE,EAAE,MAAM,CAAC;IACX,gDAAgD;IAChD,KAAK,EAAE,MAAM,CAAC;IACd,kEAAkE;IAClE,KAAK,EAAE,MAAM,CAAC;IACd,wBAAwB;IACxB,SAAS,EAAE,MAAM,EAAE,CAAC;CACrB;AAsFD;;;;;;;;;;;;GAYG;AACH,wBAAgB,GAAG,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,IAAI,CAAC,EAAE,UAAU,GAAG,SAAS,CAmC5E"}
@@ -19,6 +19,12 @@ type i32 = number;
19
19
  export interface QuadOptions {
20
20
  /** Relative tolerance on each panel's `|K - G|` vs `|K|` (default 1e-10). */
21
21
  tol?: f64;
22
+ /**
23
+ * Absolute floor on a panel's `|K - G|` below which it is accepted whatever
24
+ * `tol · |K|` is (default 1e-14). An integrand whose integral is itself below
25
+ * the floor (an SI-scaled quantity) needs a smaller floor, or 0.
26
+ */
27
+ atol?: f64;
22
28
  /** Maximum bisection recursion depth per panel (default 50). */
23
29
  maxDepth?: i32;
24
30
  }
@@ -1 +1 @@
1
- {"version":3,"file":"adaptive-quad.d.ts","sourceRoot":"","sources":["../../src/numeric/adaptive-quad.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;GAcG;AAEH,KAAK,GAAG,GAAG,MAAM,CAAC;AAClB,KAAK,GAAG,GAAG,MAAM,CAAC;AAElB,gCAAgC;AAChC,MAAM,WAAW,WAAW;IAC1B,6EAA6E;IAC7E,GAAG,CAAC,EAAE,GAAG,CAAC;IACV,gEAAgE;IAChE,QAAQ,CAAC,EAAE,GAAG,CAAC;CAChB;AAED,8BAA8B;AAC9B,MAAM,WAAW,UAAU;IACzB,8BAA8B;IAC9B,KAAK,EAAE,GAAG,CAAC;IACX,+DAA+D;IAC/D,KAAK,EAAE,GAAG,CAAC;CACZ;AAgED;;;;;;;;;;;;;;;;;;;;;GAqBG;AACH,wBAAgB,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,KAAK,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,EAAE,WAAW,GAAG,UAAU,CAmBvF"}
1
+ {"version":3,"file":"adaptive-quad.d.ts","sourceRoot":"","sources":["../../src/numeric/adaptive-quad.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;GAcG;AAEH,KAAK,GAAG,GAAG,MAAM,CAAC;AAClB,KAAK,GAAG,GAAG,MAAM,CAAC;AAElB,gCAAgC;AAChC,MAAM,WAAW,WAAW;IAC1B,6EAA6E;IAC7E,GAAG,CAAC,EAAE,GAAG,CAAC;IACV;;;;OAIG;IACH,IAAI,CAAC,EAAE,GAAG,CAAC;IACX,gEAAgE;IAChE,QAAQ,CAAC,EAAE,GAAG,CAAC;CAChB;AAED,8BAA8B;AAC9B,MAAM,WAAW,UAAU;IACzB,8BAA8B;IAC9B,KAAK,EAAE,GAAG,CAAC;IACX,+DAA+D;IAC/D,KAAK,EAAE,GAAG,CAAC;CACZ;AAgED;;;;;;;;;;;;;;;;;;;;;GAqBG;AACH,wBAAgB,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,KAAK,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,EAAE,WAAW,GAAG,UAAU,CAoBvF"}
@@ -21,8 +21,14 @@ type i32 = number;
21
21
  * Options for root-finding algorithms.
22
22
  */
23
23
  export interface FindRootOptions {
24
- /** Absolute tolerance for convergence (default 1e-12) */
24
+ /** Absolute tolerance on the residual |f(x)| (default 1e-12). */
25
25
  tol?: f64;
26
+ /**
27
+ * Absolute tolerance on the bracket width |b − a| (default: `tol`). Pass 0
28
+ * to stop on the residual alone, which honours a tolerance stated in the
29
+ * units of f when x is in other units.
30
+ */
31
+ xtol?: f64;
26
32
  /** Maximum iterations (default 100) */
27
33
  maxIter?: i32;
28
34
  }
@@ -1 +1 @@
1
- {"version":3,"file":"numeric.d.ts","sourceRoot":"","sources":["../../src/typed/numeric.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAIH,OAAO,EACL,iBAAiB,EACjB,KAAK,wBAAwB,EAC7B,KAAK,uBAAuB,GAC7B,MAAM,kCAAkC,CAAC;AAI1C,OAAO,EACL,iBAAiB,EACjB,KAAK,wBAAwB,EAC7B,KAAK,oBAAoB,EACzB,KAAK,WAAW,EAChB,KAAK,gBAAgB,EACrB,KAAK,aAAa,EAClB,KAAK,eAAe,GACrB,MAAM,iCAAiC,CAAC;AAGzC,OAAO,EACL,QAAQ,EACR,KAAK,eAAe,EACpB,KAAK,WAAW,EAChB,KAAK,eAAe,EACpB,KAAK,aAAa,EAClB,KAAK,iBAAiB,GACvB,MAAM,wBAAwB,CAAC;AAGhC,OAAO,EACL,QAAQ,EACR,KAAK,eAAe,EACpB,KAAK,WAAW,EAChB,KAAK,UAAU,EACf,KAAK,UAAU,GAChB,MAAM,wBAAwB,CAAC;AAGhC,OAAO,EACL,cAAc,EACd,KAAK,EACL,KAAK,EACL,KAAK,EACL,mBAAmB,EACnB,KAAK,qBAAqB,EAC1B,KAAK,sBAAsB,GAC5B,MAAM,+BAA+B,CAAC;AAMvC,KAAK,GAAG,GAAG,MAAM,CAAC;AAClB,KAAK,GAAG,GAAG,MAAM,CAAC;AAMlB;;GAEG;AACH,MAAM,WAAW,eAAe;IAC9B,yDAAyD;IACzD,GAAG,CAAC,EAAE,GAAG,CAAC;IACV,uCAAuC;IACvC,OAAO,CAAC,EAAE,GAAG,CAAC;CACf;AAED;;;;;;;;;;;;;;GAcG;AACH,wBAAgB,QAAQ,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,KAAK,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,EAAE,eAAe,GAAG,GAAG,CAwExF;AAED;;;;;;;;;GASG;AACH,wBAAgB,QAAQ,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE,CAyC7D;AAMD;;GAEG;AACH,MAAM,WAAW,eAAe;IAC9B,wCAAwC;IACxC,GAAG,CAAC,EAAE,GAAG,CAAC;IACV,wCAAwC;IACxC,OAAO,CAAC,EAAE,GAAG,CAAC;IACd,8CAA8C;IAC9C,IAAI,CAAC,EAAE,GAAG,CAAC;CACZ;AAED;;;;;;;;;;GAUG;AACH,wBAAgB,QAAQ,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,KAAK,GAAG,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,IAAI,CAAC,EAAE,eAAe,GAAG,MAAM,EAAE,CA+FhG;AAED;;;;;;;GAOG;AACH,wBAAgB,QAAQ,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,KAAK,GAAG,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,IAAI,CAAC,EAAE,eAAe,GAAG,MAAM,EAAE,CAEhG;AAED;;;;;;;GAOG;AACH,wBAAgB,cAAc,CAC5B,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,KAAK,GAAG,EACvB,MAAM,EAAE,CAAC,MAAM,EAAE,MAAM,CAAC,EAAE,EAC1B,IAAI,CAAC,EAAE,eAAe,GAAG;IAAE,KAAK,CAAC,EAAE,GAAG,CAAA;CAAE,GACvC,MAAM,EAAE,CAsBV;AAED;;;;;;GAMG;AACH,wBAAgB,YAAY,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE,CAuBjE;AAMD;;;;;;;;;;;;;GAaG;AACH,wBAAgB,UAAU,CACxB,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,KAAK,GAAG,EAClB,CAAC,EAAE,GAAG,EACN,CAAC,EAAE,GAAG,EACN,IAAI,CAAC,EAAE;IAAE,GAAG,CAAC,EAAE,GAAG,CAAC;IAAC,QAAQ,CAAC,EAAE,GAAG,CAAA;CAAE,GACnC,GAAG,CAEL;AAED;;;;;;;;GAQG;AACH,wBAAgB,QAAQ,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,KAAK,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,GAAE,GAAS,GAAG,GAAG,CAQ9E;AAMD;;;;;;;GAOG;AACH,wBAAgB,WAAW,CACzB,EAAE,EAAE,MAAM,EAAE,EACZ,EAAE,EAAE,MAAM,EAAE,EACZ,MAAM,GAAE,QAAQ,GAAG,UAAU,GAAG,QAAmB,GAClD,CAAC,CAAC,EAAE,GAAG,KAAK,GAAG,CA4EjB;AAED;;GAEG;AACH,wBAAgB,OAAO,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,KAAK,GAAG,CAEnE;AAED;;GAEG;AACH,wBAAgB,KAAK,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,GAAG,GAAG,GAAG,CA+C7D;AAED;;;;;;GAMG;AACH,wBAAgB,WAAW,CAAC,aAAa,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,GAAG,GAAG,MAAM,EAAE,CAiBvE;AAED;;;;;;;GAOG;AACH,wBAAgB,OAAO,CAAC,aAAa,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,MAAM,EAAE,CA2ChF;AAED;;;;;;;;GAQG;AACH,wBAAgB,KAAK,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,GAAG,EAAE,SAAS,GAAE,GAAS,GAAG,GAAG,CAuCnF;AAED;;;;;;;;GAQG;AACH,wBAAgB,QAAQ,CACtB,MAAM,EAAE,MAAM,EAAE,EAAE,EAClB,MAAM,EAAE,MAAM,EAAE,EAChB,EAAE,EAAE,MAAM,EAAE,EACZ,EAAE,EAAE,MAAM,EAAE,GACX,MAAM,EAAE,EAAE,CAgCZ;AAED;;;;;;;;GAQG;AACH,wBAAgB,cAAc,CAC5B,MAAM,EAAE,MAAM,EAAE,EAAE,EAClB,MAAM,EAAE,MAAM,EAAE,EAChB,EAAE,EAAE,MAAM,EAAE,EAAE,EACd,MAAM,GAAE,UAAU,GAAG,cAAc,GAAG,WAAwB,GAC7D,MAAM,EAAE,CA8DV;AAMD;;;;;;;;GAQG;AACH,wBAAgB,QAAQ,CACtB,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,EAAE,MAAM,EAAE,MAAM,EAAE,KAAK,GAAG,EACpC,EAAE,EAAE,MAAM,EAAE,EACZ,EAAE,EAAE,MAAM,EAAE,EACZ,EAAE,EAAE,MAAM,EAAE,GACX,MAAM,EAAE,CA+EV;AAED;;;;;;GAMG;AACH,wBAAgB,MAAM,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAiB7D;AAED;;;;;;GAMG;AACH,wBAAgB,MAAM,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAgB7D;AAED;;;;;;GAMG;AACH,wBAAgB,QAAQ,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAkB/D;AAMD;;GAEG;AACH,MAAM,WAAW,WAAW;IAC1B,CAAC,EAAE,MAAM,EAAE,CAAC;IACZ,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC;CACf;AAuGD;;;;;;;;;;;;;;GAcG;AACH,wBAAgB,cAAc,CAC5B,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,MAAM,EAAE,EACpC,EAAE,EAAE,MAAM,EAAE,EACZ,KAAK,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,EACjB,IAAI,CAAC,EAAE;IAAE,GAAG,CAAC,EAAE,GAAG,CAAC;IAAC,QAAQ,CAAC,EAAE,GAAG,CAAC;IAAC,EAAE,CAAC,EAAE,GAAG,CAAA;CAAE,GAC7C,WAAW,CAwCb;AAED;;;;;;;;;;;;;GAaG;AACH,wBAAgB,cAAc,CAC5B,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,MAAM,EAAE,EACpC,EAAE,EAAE,MAAM,EAAE,EACZ,KAAK,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,GAChB,WAAW,CAYb;AAED;;;;;;;;;;;;;;;;;;;;;;GAsBG;AACH,wBAAgB,QAAQ,CACtB,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,MAAM,EAAE,EACpC,EAAE,EAAE,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE,MAAM,EAAE,KAAK,MAAM,EAAE,EAC5C,IAAI,EAAE,MAAM,EAAE,EACd,OAAO,GAAE,MAAM,EAAW,GACzB,WAAW,CA8Cb;AAED;;;;;;;;;GASG;AACH,wBAAgB,eAAe,CAC7B,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,MAAM,EAAE,EACpC,EAAE,EAAE,MAAM,EAAE,EACZ,EAAE,EAAE,GAAG,EACP,CAAC,EAAE,GAAG,EACN,GAAG,GAAE,GAAU,GACd;IAAE,CAAC,EAAE,MAAM,EAAE,CAAC;IAAC,CAAC,EAAE,GAAG,CAAC;IAAC,CAAC,EAAE,GAAG,CAAA;CAAE,CAgDjC;AAED;;;;;;;;;GASG;AACH,wBAAgB,cAAc,CAC5B,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,MAAM,EAAE,EACpC,EAAE,EAAE,MAAM,EAAE,EACZ,KAAK,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,EACjB,KAAK,EAAE,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,GAAG,GAClC,WAAW,GAAG;IAAE,SAAS,CAAC,EAAE,GAAG,CAAA;CAAE,CA6EnC;AAMD;;;;;;GAMG;AACH,wBAAgB,IAAI,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,GAAG,GAAE,GAAW,GAAG,GAAG,CA6BzD;AAED;;;;;GAKG;AACH,wBAAgB,SAAS,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,GAAG,MAAM,EAAE,EAAE,CAqDnD;AAED;;;;;;;GAOG;AACH,wBAAgB,OAAO,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG;IAAE,QAAQ,EAAE,MAAM,EAAE,CAAC;IAAC,KAAK,EAAE,MAAM,EAAE,CAAA;CAAE,CAwBzF;AAwDD;;;;;;;;GAQG;AACH,wBAAgB,eAAe,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,KAAK,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,GAAE,GAAQ,GAAG,CAAC,CAAC,EAAE,GAAG,KAAK,GAAG,CAiChG;AAED;;;;;;;GAOG;AACH,wBAAgB,eAAe,CAC7B,MAAM,EAAE,MAAM,EAAE,EAChB,CAAC,EAAE,GAAG,EACN,CAAC,EAAE,GAAG,GACL;IAAE,GAAG,EAAE,MAAM,EAAE,CAAC;IAAC,GAAG,EAAE,MAAM,EAAE,CAAA;CAAE,CA2ClC;AAED;;;;;;;;;GASG;AACH,wBAAgB,QAAQ,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE,CAuCzF;AAED;;;GAGG;AACH,MAAM,WAAW,cAAc;IAC7B,IAAI,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC;IAClB,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC;IAChB,IAAI,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC;IAClB,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC;IAChB,mGAAmG;IACnG,MAAM,CAAC,EAAE,SAAS,CAAC,SAAS,CAAC,MAAM,GAAG,IAAI,EAAE,MAAM,GAAG,IAAI,CAAC,CAAC,EAAE,CAAC;CAC/D;AAED,qDAAqD;AACrD,MAAM,WAAW,aAAa;IAC5B,CAAC,EAAE,MAAM,EAAE,CAAC;IACZ,GAAG,EAAE,MAAM,CAAC;IACZ,OAAO,EAAE,OAAO,CAAC;IACjB,MAAM,EAAE,SAAS,GAAG,YAAY,GAAG,WAAW,CAAC;CAChD;AA4bD;;;;;;;;;GASG;AACH,wBAAgB,OAAO,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,IAAI,EAAE,MAAM,EAAE,EAAE,EAAE,IAAI,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE,GAAG,IAAI,CAAC;AACxF,wBAAgB,OAAO,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,IAAI,EAAE,cAAc,GAAG,aAAa,CAAC;AAkB1E;;;;;;;;;GASG;AACH,wBAAgB,QAAQ,CACtB,GAAG,EAAE;IAAE,KAAK,EAAE,GAAG,CAAA;CAAE,EACnB,MAAM,EAAE;IAAE,CAAC,EAAE,GAAG,CAAC;IAAC,EAAE,EAAE,GAAG,CAAC;IAAC,EAAE,EAAE,GAAG,CAAC;IAAC,CAAC,EAAE,GAAG,CAAA;CAAE,EAC5C,EAAE,EAAE;IAAE,IAAI,EAAE,GAAG,CAAC;IAAC,KAAK,EAAE,GAAG,CAAC;IAAC,OAAO,EAAE,CAAC,CAAC,EAAE,GAAG,KAAK,GAAG,CAAA;CAAE,GACtD;IAAE,CAAC,EAAE,MAAM,EAAE,CAAC;IAAC,CAAC,EAAE,MAAM,EAAE,CAAA;CAAE,CA+B9B"}
1
+ {"version":3,"file":"numeric.d.ts","sourceRoot":"","sources":["../../src/typed/numeric.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAIH,OAAO,EACL,iBAAiB,EACjB,KAAK,wBAAwB,EAC7B,KAAK,uBAAuB,GAC7B,MAAM,kCAAkC,CAAC;AAI1C,OAAO,EACL,iBAAiB,EACjB,KAAK,wBAAwB,EAC7B,KAAK,oBAAoB,EACzB,KAAK,WAAW,EAChB,KAAK,gBAAgB,EACrB,KAAK,aAAa,EAClB,KAAK,eAAe,GACrB,MAAM,iCAAiC,CAAC;AAGzC,OAAO,EACL,QAAQ,EACR,KAAK,eAAe,EACpB,KAAK,WAAW,EAChB,KAAK,eAAe,EACpB,KAAK,aAAa,EAClB,KAAK,iBAAiB,GACvB,MAAM,wBAAwB,CAAC;AAGhC,OAAO,EACL,QAAQ,EACR,KAAK,eAAe,EACpB,KAAK,WAAW,EAChB,KAAK,UAAU,EACf,KAAK,UAAU,GAChB,MAAM,wBAAwB,CAAC;AAGhC,OAAO,EACL,cAAc,EACd,KAAK,EACL,KAAK,EACL,KAAK,EACL,mBAAmB,EACnB,KAAK,qBAAqB,EAC1B,KAAK,sBAAsB,GAC5B,MAAM,+BAA+B,CAAC;AAMvC,KAAK,GAAG,GAAG,MAAM,CAAC;AAClB,KAAK,GAAG,GAAG,MAAM,CAAC;AAMlB;;GAEG;AACH,MAAM,WAAW,eAAe;IAC9B,iEAAiE;IACjE,GAAG,CAAC,EAAE,GAAG,CAAC;IACV;;;;OAIG;IACH,IAAI,CAAC,EAAE,GAAG,CAAC;IACX,uCAAuC;IACvC,OAAO,CAAC,EAAE,GAAG,CAAC;CACf;AAED;;;;;;;;;;;;;;GAcG;AACH,wBAAgB,QAAQ,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,KAAK,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,EAAE,eAAe,GAAG,GAAG,CA4ExF;AAED;;;;;;;;;GASG;AACH,wBAAgB,QAAQ,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE,CAyC7D;AAMD;;GAEG;AACH,MAAM,WAAW,eAAe;IAC9B,wCAAwC;IACxC,GAAG,CAAC,EAAE,GAAG,CAAC;IACV,wCAAwC;IACxC,OAAO,CAAC,EAAE,GAAG,CAAC;IACd,8CAA8C;IAC9C,IAAI,CAAC,EAAE,GAAG,CAAC;CACZ;AAED;;;;;;;;;;GAUG;AACH,wBAAgB,QAAQ,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,KAAK,GAAG,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,IAAI,CAAC,EAAE,eAAe,GAAG,MAAM,EAAE,CA+FhG;AAED;;;;;;;GAOG;AACH,wBAAgB,QAAQ,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,KAAK,GAAG,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,IAAI,CAAC,EAAE,eAAe,GAAG,MAAM,EAAE,CAEhG;AAED;;;;;;;GAOG;AACH,wBAAgB,cAAc,CAC5B,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,KAAK,GAAG,EACvB,MAAM,EAAE,CAAC,MAAM,EAAE,MAAM,CAAC,EAAE,EAC1B,IAAI,CAAC,EAAE,eAAe,GAAG;IAAE,KAAK,CAAC,EAAE,GAAG,CAAA;CAAE,GACvC,MAAM,EAAE,CAsBV;AAED;;;;;;GAMG;AACH,wBAAgB,YAAY,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE,CAuBjE;AAMD;;;;;;;;;;;;;GAaG;AACH,wBAAgB,UAAU,CACxB,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,KAAK,GAAG,EAClB,CAAC,EAAE,GAAG,EACN,CAAC,EAAE,GAAG,EACN,IAAI,CAAC,EAAE;IAAE,GAAG,CAAC,EAAE,GAAG,CAAC;IAAC,QAAQ,CAAC,EAAE,GAAG,CAAA;CAAE,GACnC,GAAG,CAEL;AAED;;;;;;;;GAQG;AACH,wBAAgB,QAAQ,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,KAAK,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,GAAE,GAAS,GAAG,GAAG,CAQ9E;AAMD;;;;;;;GAOG;AACH,wBAAgB,WAAW,CACzB,EAAE,EAAE,MAAM,EAAE,EACZ,EAAE,EAAE,MAAM,EAAE,EACZ,MAAM,GAAE,QAAQ,GAAG,UAAU,GAAG,QAAmB,GAClD,CAAC,CAAC,EAAE,GAAG,KAAK,GAAG,CA4EjB;AAED;;GAEG;AACH,wBAAgB,OAAO,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,KAAK,GAAG,CAEnE;AAED;;GAEG;AACH,wBAAgB,KAAK,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,GAAG,GAAG,GAAG,CA+C7D;AAED;;;;;;GAMG;AACH,wBAAgB,WAAW,CAAC,aAAa,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,GAAG,GAAG,MAAM,EAAE,CAiBvE;AAED;;;;;;;GAOG;AACH,wBAAgB,OAAO,CAAC,aAAa,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,MAAM,EAAE,CA2ChF;AAED;;;;;;;;GAQG;AACH,wBAAgB,KAAK,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,GAAG,EAAE,SAAS,GAAE,GAAS,GAAG,GAAG,CAuCnF;AAED;;;;;;;;GAQG;AACH,wBAAgB,QAAQ,CACtB,MAAM,EAAE,MAAM,EAAE,EAAE,EAClB,MAAM,EAAE,MAAM,EAAE,EAChB,EAAE,EAAE,MAAM,EAAE,EACZ,EAAE,EAAE,MAAM,EAAE,GACX,MAAM,EAAE,EAAE,CAgCZ;AAED;;;;;;;;GAQG;AACH,wBAAgB,cAAc,CAC5B,MAAM,EAAE,MAAM,EAAE,EAAE,EAClB,MAAM,EAAE,MAAM,EAAE,EAChB,EAAE,EAAE,MAAM,EAAE,EAAE,EACd,MAAM,GAAE,UAAU,GAAG,cAAc,GAAG,WAAwB,GAC7D,MAAM,EAAE,CA8DV;AAMD;;;;;;;;GAQG;AACH,wBAAgB,QAAQ,CACtB,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,EAAE,MAAM,EAAE,MAAM,EAAE,KAAK,GAAG,EACpC,EAAE,EAAE,MAAM,EAAE,EACZ,EAAE,EAAE,MAAM,EAAE,EACZ,EAAE,EAAE,MAAM,EAAE,GACX,MAAM,EAAE,CA+EV;AAED;;;;;;GAMG;AACH,wBAAgB,MAAM,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAiB7D;AAED;;;;;;GAMG;AACH,wBAAgB,MAAM,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAgB7D;AAED;;;;;;GAMG;AACH,wBAAgB,QAAQ,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAkB/D;AAMD;;GAEG;AACH,MAAM,WAAW,WAAW;IAC1B,CAAC,EAAE,MAAM,EAAE,CAAC;IACZ,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC;CACf;AAuGD;;;;;;;;;;;;;;GAcG;AACH,wBAAgB,cAAc,CAC5B,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,MAAM,EAAE,EACpC,EAAE,EAAE,MAAM,EAAE,EACZ,KAAK,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,EACjB,IAAI,CAAC,EAAE;IAAE,GAAG,CAAC,EAAE,GAAG,CAAC;IAAC,QAAQ,CAAC,EAAE,GAAG,CAAC;IAAC,EAAE,CAAC,EAAE,GAAG,CAAA;CAAE,GAC7C,WAAW,CAwCb;AAED;;;;;;;;;;;;;GAaG;AACH,wBAAgB,cAAc,CAC5B,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,MAAM,EAAE,EACpC,EAAE,EAAE,MAAM,EAAE,EACZ,KAAK,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,GAChB,WAAW,CAYb;AAED;;;;;;;;;;;;;;;;;;;;;;GAsBG;AACH,wBAAgB,QAAQ,CACtB,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,MAAM,EAAE,EACpC,EAAE,EAAE,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE,MAAM,EAAE,KAAK,MAAM,EAAE,EAC5C,IAAI,EAAE,MAAM,EAAE,EACd,OAAO,GAAE,MAAM,EAAW,GACzB,WAAW,CA8Cb;AAED;;;;;;;;;GASG;AACH,wBAAgB,eAAe,CAC7B,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,MAAM,EAAE,EACpC,EAAE,EAAE,MAAM,EAAE,EACZ,EAAE,EAAE,GAAG,EACP,CAAC,EAAE,GAAG,EACN,GAAG,GAAE,GAAU,GACd;IAAE,CAAC,EAAE,MAAM,EAAE,CAAC;IAAC,CAAC,EAAE,GAAG,CAAC;IAAC,CAAC,EAAE,GAAG,CAAA;CAAE,CAgDjC;AAED;;;;;;;;;GASG;AACH,wBAAgB,cAAc,CAC5B,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,MAAM,EAAE,EACpC,EAAE,EAAE,MAAM,EAAE,EACZ,KAAK,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,EACjB,KAAK,EAAE,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,GAAG,GAClC,WAAW,GAAG;IAAE,SAAS,CAAC,EAAE,GAAG,CAAA;CAAE,CA6EnC;AAMD;;;;;;GAMG;AACH,wBAAgB,IAAI,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,GAAG,GAAE,GAAW,GAAG,GAAG,CA6BzD;AAED;;;;;GAKG;AACH,wBAAgB,SAAS,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,GAAG,MAAM,EAAE,EAAE,CAqDnD;AAED;;;;;;;GAOG;AACH,wBAAgB,OAAO,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG;IAAE,QAAQ,EAAE,MAAM,EAAE,CAAC;IAAC,KAAK,EAAE,MAAM,EAAE,CAAA;CAAE,CAwBzF;AAwDD;;;;;;;;GAQG;AACH,wBAAgB,eAAe,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,KAAK,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,GAAE,GAAQ,GAAG,CAAC,CAAC,EAAE,GAAG,KAAK,GAAG,CAiChG;AAED;;;;;;;GAOG;AACH,wBAAgB,eAAe,CAC7B,MAAM,EAAE,MAAM,EAAE,EAChB,CAAC,EAAE,GAAG,EACN,CAAC,EAAE,GAAG,GACL;IAAE,GAAG,EAAE,MAAM,EAAE,CAAC;IAAC,GAAG,EAAE,MAAM,EAAE,CAAA;CAAE,CA2ClC;AAED;;;;;;;;;GASG;AACH,wBAAgB,QAAQ,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE,CAuCzF;AAED;;;GAGG;AACH,MAAM,WAAW,cAAc;IAC7B,IAAI,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC;IAClB,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC;IAChB,IAAI,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC;IAClB,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC;IAChB,mGAAmG;IACnG,MAAM,CAAC,EAAE,SAAS,CAAC,SAAS,CAAC,MAAM,GAAG,IAAI,EAAE,MAAM,GAAG,IAAI,CAAC,CAAC,EAAE,CAAC;CAC/D;AAED,qDAAqD;AACrD,MAAM,WAAW,aAAa;IAC5B,CAAC,EAAE,MAAM,EAAE,CAAC;IACZ,GAAG,EAAE,MAAM,CAAC;IACZ,OAAO,EAAE,OAAO,CAAC;IACjB,MAAM,EAAE,SAAS,GAAG,YAAY,GAAG,WAAW,CAAC;CAChD;AA4bD;;;;;;;;;GASG;AACH,wBAAgB,OAAO,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,IAAI,EAAE,MAAM,EAAE,EAAE,EAAE,IAAI,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE,GAAG,IAAI,CAAC;AACxF,wBAAgB,OAAO,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,IAAI,EAAE,cAAc,GAAG,aAAa,CAAC;AAkB1E;;;;;;;;;GASG;AACH,wBAAgB,QAAQ,CACtB,GAAG,EAAE;IAAE,KAAK,EAAE,GAAG,CAAA;CAAE,EACnB,MAAM,EAAE;IAAE,CAAC,EAAE,GAAG,CAAC;IAAC,EAAE,EAAE,GAAG,CAAC;IAAC,EAAE,EAAE,GAAG,CAAC;IAAC,CAAC,EAAE,GAAG,CAAA;CAAE,EAC5C,EAAE,EAAE;IAAE,IAAI,EAAE,GAAG,CAAC;IAAC,KAAK,EAAE,GAAG,CAAC;IAAC,OAAO,EAAE,CAAC,CAAC,EAAE,GAAG,KAAK,GAAG,CAAA;CAAE,GACtD;IAAE,CAAC,EAAE,MAAM,EAAE,CAAC;IAAC,CAAC,EAAE,MAAM,EAAE,CAAA;CAAE,CA+B9B"}
@@ -106,7 +106,14 @@ export declare function factorial(n: number): number;
106
106
  export declare function betaScalar(a: f64, b: f64): f64;
107
107
  /** Regularized lower incomplete gamma function P(a, x). */
108
108
  export declare function gammaincScalar(a: f64, x: f64): f64;
109
- /** Upper regularized incomplete gamma Q(a, x) = 1 - P(a, x). */
109
+ /**
110
+ * Upper regularized incomplete gamma Q(a, x).
111
+ *
112
+ * Computed DIRECTLY in the upper tail (x >= a + 1) by the modified Lentz
113
+ * continued fraction, and as 1 − P below it where P is the series. `1 − P`
114
+ * everywhere rounded a far tail to 0 (Q(5, 80) is 2.5e-28, and 1 − P gave 0);
115
+ * a p-value of a bad χ² fit is that tail.
116
+ */
110
117
  export declare function gammaincpScalar(a: f64, x: f64): f64;
111
118
  /** Regularized incomplete beta function I_x(a, b). */
112
119
  export declare function betaincScalar(a: f64, b: f64, x: f64): f64;
@@ -1 +1 @@
1
- {"version":3,"file":"scalars.d.ts","sourceRoot":"","sources":["../../../src/wasm/special/scalars.ts"],"names":[],"mappings":"AAEA;;;;;;;;;;;;;;;;;;;;;;GAsBG;AAEH,0CAA0C;AAC1C,KAAK,GAAG,GAAG,MAAM,CAAC;AAMlB;;;;GAIG;AACH,wBAAgB,OAAO,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAwBnC;AAWD;;;;;;;;;GASG;AACH,wBAAgB,YAAY,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,KAAK,EAAE,OAAO,GAAG,GAAG,CAyBjE;AAED;;;;GAIG;AACH,wBAAgB,cAAc,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAU1C;AAED;;;;GAIG;AACH,wBAAgB,cAAc,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAU1C;AAGD,2FAA2F;AAC3F,wBAAgB,cAAc,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAiB1C;AAID,2FAA2F;AAC3F,wBAAgB,cAAc,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAmB1C;AAED,yDAAyD;AACzD,wBAAgB,cAAc,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAG1C;AAED,yDAAyD;AACzD,wBAAgB,cAAc,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAK1C;AAED,0DAA0D;AAC1D,wBAAgB,cAAc,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAG1C;AAED,0DAA0D;AAC1D,wBAAgB,cAAc,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAG1C;AAED,0EAA0E;AAC1E,wBAAgB,aAAa,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,GAAG,CAoCjD;AAED;;;;;;;GAOG;AACH,wBAAgB,aAAa,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,GAAG,CAejD;AAMD,iEAAiE;AACjE,wBAAgB,eAAe,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAmB3C;AAED,mFAAmF;AACnF,wBAAgB,uBAAuB,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAyBnD;AASD;;;;;GAKG;AACH,wBAAgB,WAAW,IAAI,MAAM,EAAE,CAMtC;AAED;;;;;GAKG;AACH,wBAAgB,UAAU,CAAC,IAAI,EAAE,GAAG,GAAG;IAAE,CAAC,EAAE,GAAG,CAAC;IAAC,CAAC,EAAE,GAAG,CAAA;CAAE,CAsBxD;AAED,6CAA6C;AAC7C,wBAAgB,YAAY,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CA0DxC;AAED,8CAA8C;AAC9C,wBAAgB,YAAY,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAqDxC;AAYD,2CAA2C;AAC3C,wBAAgB,SAAS,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAK3C;AAED,gEAAgE;AAChE,wBAAgB,UAAU,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,GAAG,CAE9C;AAED,2DAA2D;AAC3D,wBAAgB,cAAc,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,GAAG,CAkClD;AAED,gEAAgE;AAChE,wBAAgB,eAAe,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,GAAG,CAEnD;AAED,sDAAsD;AACtD,wBAAgB,aAAa,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,GAAG,CAyCzD;AAED,0DAA0D;AAC1D,wBAAgB,aAAa,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,GAAG,CAYjD"}
1
+ {"version":3,"file":"scalars.d.ts","sourceRoot":"","sources":["../../../src/wasm/special/scalars.ts"],"names":[],"mappings":"AAEA;;;;;;;;;;;;;;;;;;;;;;GAsBG;AAEH,0CAA0C;AAC1C,KAAK,GAAG,GAAG,MAAM,CAAC;AAMlB;;;;GAIG;AACH,wBAAgB,OAAO,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAwBnC;AAWD;;;;;;;;;GASG;AACH,wBAAgB,YAAY,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,KAAK,EAAE,OAAO,GAAG,GAAG,CAyBjE;AAED;;;;GAIG;AACH,wBAAgB,cAAc,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAU1C;AAED;;;;GAIG;AACH,wBAAgB,cAAc,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAU1C;AAGD,2FAA2F;AAC3F,wBAAgB,cAAc,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAiB1C;AAID,2FAA2F;AAC3F,wBAAgB,cAAc,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAmB1C;AAED,yDAAyD;AACzD,wBAAgB,cAAc,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAG1C;AAED,yDAAyD;AACzD,wBAAgB,cAAc,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAK1C;AAED,0DAA0D;AAC1D,wBAAgB,cAAc,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAG1C;AAED,0DAA0D;AAC1D,wBAAgB,cAAc,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAG1C;AAED,0EAA0E;AAC1E,wBAAgB,aAAa,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,GAAG,CAoCjD;AAED;;;;;;;GAOG;AACH,wBAAgB,aAAa,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,GAAG,CAejD;AAMD,iEAAiE;AACjE,wBAAgB,eAAe,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAmB3C;AAED,mFAAmF;AACnF,wBAAgB,uBAAuB,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAyBnD;AASD;;;;;GAKG;AACH,wBAAgB,WAAW,IAAI,MAAM,EAAE,CAMtC;AAED;;;;;GAKG;AACH,wBAAgB,UAAU,CAAC,IAAI,EAAE,GAAG,GAAG;IAAE,CAAC,EAAE,GAAG,CAAC;IAAC,CAAC,EAAE,GAAG,CAAA;CAAE,CAsBxD;AAED,6CAA6C;AAC7C,wBAAgB,YAAY,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CA0DxC;AAED,8CAA8C;AAC9C,wBAAgB,YAAY,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAqDxC;AAYD,2CAA2C;AAC3C,wBAAgB,SAAS,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAK3C;AAED,gEAAgE;AAChE,wBAAgB,UAAU,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,GAAG,CAE9C;AAED,2DAA2D;AAC3D,wBAAgB,cAAc,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,GAAG,CAkClD;AAED;;;;;;;GAOG;AACH,wBAAgB,eAAe,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,GAAG,CAsBnD;AAED,sDAAsD;AACtD,wBAAgB,aAAa,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,GAAG,CAyCzD;AAED,0DAA0D;AAC1D,wBAAgB,aAAa,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,GAAG,CAYjD"}
package/package.json CHANGED
@@ -1,6 +1,6 @@
1
1
  {
2
2
  "name": "@danielsimonjr/mathts-functions",
3
- "version": "0.68.2",
3
+ "version": "0.68.4",
4
4
  "description": "Mathematical functions for MathTS - arithmetic, algebra, trigonometry, statistics, and more",
5
5
  "author": "Daniel Simon Jr.",
6
6
  "license": "MIT",
@@ -34,7 +34,7 @@
34
34
  "build:prod": "tsup src/index.ts --format esm --clean --minify --treeshake && tsc -p tsconfig.dts.json && node scripts/copy-wasm.mjs"
35
35
  },
36
36
  "dependencies": {
37
- "@danielsimonjr/mathts-core": "^0.18.0",
37
+ "@danielsimonjr/mathts-core": "^0.18.1",
38
38
  "@danielsimonjr/mathts-expression": "^0.10.2",
39
39
  "@danielsimonjr/mathts-gpu": "^0.2.1",
40
40
  "@danielsimonjr/mathts-matrix": "^0.7.8",