@milaboratories/miplots4 1.0.176 → 1.0.178

package/dist/node_modules/.pnpm/@stdlib_math-base-special-kernel-betaincinv@0.2.2/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/ibeta_roots.js

@@ -1,53 +0,0 @@
-import { __require as A } from "../../../../../@stdlib_math-base-special-kernel-betainc@0.2.2/node_modules/@stdlib/math-base-special-kernel-betainc/lib/index.js";
-import { __require as I } from "../../../../../@stdlib_math-base-special-abs@0.2.2/node_modules/@stdlib/math-base-special-abs/lib/index.js";
-import { __require as M } from "../../../../../@stdlib_constants-float64-max@0.2.2/node_modules/@stdlib/constants-float64-max/lib/index.js";
-import { __require as O } from "../../../../../@stdlib_constants-float64-smallest-normal@0.2.2/node_modules/@stdlib/constants-float64-smallest-normal/lib/index.js";
-/**
-* @license Apache-2.0
-*
-* Copyright (c) 2018 The Stdlib Authors.
-*
-* Licensed under the Apache License, Version 2.0 (the "License");
-* you may not use this file except in compliance with the License.
-* You may obtain a copy of the License at
-*
-*    http://www.apache.org/licenses/LICENSE-2.0
-*
-* Unless required by applicable law or agreed to in writing, software
-* distributed under the License is distributed on an "AS IS" BASIS,
-* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
-* See the License for the specific language governing permissions and
-* limitations under the License.
-*
-*
-* ## Notice
-*
-* The original C++ code and copyright notice are from the [Boost library]{@link http://www.boost.org/doc/libs/1_64_0/boost/math/special_functions/detail/ibeta_inverse.hpp}. The implementation has been modified for JavaScript.
-*
-* ```text
-* Copyright John Maddock 2006.
-* Copyright Paul A. Bristow 2007.
-*
-* Use, modification and distribution are subject to the
-* Boost Software License, Version 1.0. (See accompanying file
-* LICENSE or copy at http://www.boost.org/LICENSE_1_0.txt)
-* ```
-*/
-var f, n;
-function T() {
-  if (n) return f;
-  n = 1;
-  var b = A().assign, m = I(), v = M(), u = O();
-  function L(_, s, p, o) {
-    return c;
-    function c(r) {
-      var t, e, i, q, a;
-      return a = 1 - r, t = [0, 0], b(r, _, s, !0, o, t, 1, 0), q = t[0] - p, e = t[1], o && (e = -e), a === 0 && (a = u * 64), r === 0 && (r = u * 64), i = e * (-(a * _) + (s - 2) * r + 1), m(i) < a * r * v && (i /= a * r), o && (i = -i), e === 0 && (e = (o ? -1 : 1) * u * 64), [q, e, i];
-    }
-  }
-  return f = L, f;
-}
-export {
-  T as __require
-};
-//# sourceMappingURL=ibeta_roots.js.map

package/dist/node_modules/.pnpm/@stdlib_math-base-special-kernel-betaincinv@0.2.2/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/ibeta_roots.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"ibeta_roots.js","sources":["../../../../../../../../../../node_modules/.pnpm/@stdlib+math-base-special-kernel-betaincinv@0.2.2/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/ibeta_roots.js"],"sourcesContent":["/**\n* @license Apache-2.0\n*\n* Copyright (c) 2018 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*\n*\n* ## Notice\n*\n* The original C++ code and copyright notice are from the [Boost library]{@link http://www.boost.org/doc/libs/1_64_0/boost/math/special_functions/detail/ibeta_inverse.hpp}. The implementation has been modified for JavaScript.\n*\n* ```text\n* Copyright John Maddock 2006.\n* Copyright Paul A. Bristow 2007.\n*\n* Use, modification and distribution are subject to the\n* Boost Software License, Version 1.0. (See accompanying file\n* LICENSE or copy at http://www.boost.org/LICENSE_1_0.txt)\n* ```\n*/\n\n'use strict';\n\n// MODULES //\n\nvar kernelBetainc = require( '@stdlib/math-base-special-kernel-betainc' ).assign;\nvar abs = require( '@stdlib/math-base-special-abs' );\nvar FLOAT64_MAX = require( '@stdlib/constants-float64-max' );\nvar FLOAT64_MIN_NORM = require( '@stdlib/constants-float64-smallest-normal' );\n\n\n// MAIN //\n\n/**\n* Returns a root finding function.\n*\n* @private\n* @param {PositiveNumber} a - function parameter\n* @param {PositiveNumber} b - function parameter\n* @param {Probability} target - probability value\n* @param {boolean} invert - boolean indicating whether to find the roots of the upper or lower incomplete beta function\n* @returns {Function} root finding function\n*/\nfunction ibetaRoots( a, b, target, invert ) {\n\treturn roots;\n\n\t/**\n\t* Calculates roots.\n\t*\n\t* @private\n\t* @param {number} x - input value\n\t* @returns {Array} roots\n\t*/\n\tfunction roots( x ) {\n\t\tvar buf;\n\t\tvar f1;\n\t\tvar f2;\n\t\tvar f;\n\t\tvar y;\n\n\t\ty = 1.0 - x;\n\n\t\tbuf = [ 0.0, 0.0 ];\n\t\tkernelBetainc( x, a, b, true, invert, buf, 1, 0 );\n\t\tf = buf[ 0 ] - target;\n\t\tf1 = buf[ 1 ];\n\t\tif ( invert ) {\n\t\t\tf1 = -f1;\n\t\t}\n\t\tif ( y === 0.0 ) {\n\t\t\ty = FLOAT64_MIN_NORM * 64.0;\n\t\t}\n\t\tif ( x === 0.0 ) {\n\t\t\tx = FLOAT64_MIN_NORM * 64.0;\n\t\t}\n\t\tf2 = f1 * ( -(y*a) + ( ( b-2.0 ) * x ) + 1.0 );\n\t\tif ( abs( f2 ) < y * x * FLOAT64_MAX ) {\n\t\t\tf2 /= (y * x);\n\t\t}\n\t\tif ( invert ) {\n\t\t\tf2 = -f2;\n\t\t}\n\t\t// Make sure we don't have a zero derivative:\n\t\tif ( f1 === 0.0 ) {\n\t\t\tf1 = ( ( invert ) ? -1.0 : 1.0 ) * FLOAT64_MIN_NORM * 64.0;\n\t\t}\n\t\treturn [ f, f1, f2 ];\n\t}\n}\n\n\n// EXPORTS //\n\nmodule.exports = ibetaRoots;\n"],"names":["kernelBetainc","require$$0","abs","require$$1","FLOAT64_MAX","require$$2","FLOAT64_MIN_NORM","require$$3","ibetaRoots","a","b","target","invert","roots","x","buf","f1","f2","f","y","ibeta_roots"],"mappings":";;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;AAoCA,MAAIA,IAAgBC,EAAA,EAAsD,QACtEC,IAAMC,EAAA,GACNC,IAAcC,EAAA,GACdC,IAAmBC,EAAA;AAevB,WAASC,EAAYC,GAAGC,GAAGC,GAAQC,GAAS;AAC3C,WAAOC;AASP,aAASA,EAAOC,GAAI;AACnB,UAAIC,GACAC,GACAC,GACAC,GACAC;AAEJ,aAAAA,IAAI,IAAML,GAEVC,IAAM,CAAE,GAAK,CAAG,GAChBf,EAAec,GAAGL,GAAGC,GAAG,IAAME,GAAQG,GAAK,GAAG,CAAC,GAC/CG,IAAIH,EAAK,CAAC,IAAKJ,GACfK,IAAKD,EAAK,CAAC,GACNH,MACJI,IAAK,CAACA,IAEFG,MAAM,MACVA,IAAIb,IAAmB,KAEnBQ,MAAM,MACVA,IAAIR,IAAmB,KAExBW,IAAKD,KAAO,EAAEG,IAAEV,MAASC,IAAE,KAAQI,IAAM,IACpCZ,EAAKe,CAAE,IAAKE,IAAIL,IAAIV,MACxBa,KAAOE,IAAIL,IAEPF,MACJK,IAAK,CAACA,IAGFD,MAAO,MACXA,KAASJ,IAAW,KAAO,KAAQN,IAAmB,KAEhD,CAAEY,GAAGF,GAAIC,CAAE;AAAA,IACpB;AAAA,EACA;AAKA,SAAAG,IAAiBZ;;","x_google_ignoreList":[0]}

package/dist/node_modules/.pnpm/@stdlib_math-base-special-kernel-betaincinv@0.2.2/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/index.js

@@ -1,29 +0,0 @@
-import { __require as a } from "./main.js";
-/**
-* @license Apache-2.0
-*
-* Copyright (c) 2018 The Stdlib Authors.
-*
-* Licensed under the Apache License, Version 2.0 (the "License");
-* you may not use this file except in compliance with the License.
-* You may obtain a copy of the License at
-*
-*    http://www.apache.org/licenses/LICENSE-2.0
-*
-* Unless required by applicable law or agreed to in writing, software
-* distributed under the License is distributed on an "AS IS" BASIS,
-* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
-* See the License for the specific language governing permissions and
-* limitations under the License.
-*/
-var r, e;
-function n() {
-  if (e) return r;
-  e = 1;
-  var i = a();
-  return r = i, r;
-}
-export {
-  n as __require
-};
-//# sourceMappingURL=index.js.map

package/dist/node_modules/.pnpm/@stdlib_math-base-special-kernel-betaincinv@0.2.2/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/index.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"index.js","sources":["../../../../../../../../../../node_modules/.pnpm/@stdlib+math-base-special-kernel-betaincinv@0.2.2/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/index.js"],"sourcesContent":["/**\n* @license Apache-2.0\n*\n* Copyright (c) 2018 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*/\n\n'use strict';\n\n/**\n* Evaluate the inverse of the incomplete beta function.\n*\n* @module @stdlib/math-base-special-kernel-betaincinv\n*\n* @example\n* var kernelBetaincinv = require( '@stdlib/math-base-special-kernel-betaincinv' );\n*\n* var y = kernelBetaincinv( 3.0, 3.0, 0.2, 0.8 );\n* // returns [ ~0.327, ~0.673 ]\n*\n* y = kernelBetaincinv( 3.0, 3.0, 0.4, 0.6 );\n* // returns [ ~0.446, ~0.554 ]\n*\n* y = kernelBetaincinv( 1.0, 6.0, 0.4, 0.6 );\n* // returns [ ~0.082, ~0.918 ]\n*\n* y = kernelBetaincinv( 1.0, 6.0, 0.8, 0.2 );\n* // returns [ ~0.235, ~0.765 ]\n*/\n\n// MODULES //\n\nvar main = require( './main.js' );\n\n\n// EXPORTS //\n\nmodule.exports = main;\n"],"names":["main","require$$0","lib"],"mappings":";;;;;;;;;;;;;;;;;;;;;;AA2CA,MAAIA,IAAOC,EAAA;AAKX,SAAAC,IAAiBF;;","x_google_ignoreList":[0]}

package/dist/node_modules/.pnpm/@stdlib_math-base-special-kernel-betaincinv@0.2.2/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/inverse_students_t.js

@@ -1,85 +0,0 @@
-import { __require as g } from "../../../../../@stdlib_math-base-special-erfcinv@0.2.3/node_modules/@stdlib/math-base-special-erfcinv/lib/index.js";
-import { __require as j } from "../../../../../@stdlib_math-base-special-floor@0.2.3/node_modules/@stdlib/math-base-special-floor/lib/index.js";
-import { __require as z } from "../../../../../@stdlib_math-base-special-ldexp@0.2.3/node_modules/@stdlib/math-base-special-ldexp/lib/index.js";
-import { __require as G } from "../../../../../@stdlib_math-base-special-round@0.3.0/node_modules/@stdlib/math-base-special-round/lib/index.js";
-import { __require as J } from "../../../../../@stdlib_math-base-special-acos@0.2.3/node_modules/@stdlib/math-base-special-acos/lib/index.js";
-import { __require as K } from "../../../../../@stdlib_math-base-special-sqrt@0.2.2/node_modules/@stdlib/math-base-special-sqrt/lib/index.js";
-import { __require as M } from "../../../../../@stdlib_math-base-special-abs@0.2.2/node_modules/@stdlib/math-base-special-abs/lib/index.js";
-import { __require as U } from "../../../../../@stdlib_math-base-special-cos@0.3.0/node_modules/@stdlib/math-base-special-cos/lib/index.js";
-import { __require as V } from "../../../../../@stdlib_math-base-special-pow@0.3.0/node_modules/@stdlib/math-base-special-pow/lib/index.js";
-import { __require as W } from "../../../../../@stdlib_math-base-special-sin@0.3.0/node_modules/@stdlib/math-base-special-sin/lib/index.js";
-import { __require as Y } from "../../../../../@stdlib_constants-float64-sqrt-two@0.2.2/node_modules/@stdlib/constants-float64-sqrt-two/lib/index.js";
-import { __require as Z } from "../../../../../@stdlib_constants-float64-pi@0.2.2/node_modules/@stdlib/constants-float64-pi/lib/index.js";
-import { __require as d } from "./inverse_students_t_body_series.js";
-import { __require as rr } from "./inverse_students_t_tail_series.js";
-import { __require as er } from "./inverse_students_t_hill.js";
-/**
-* @license Apache-2.0
-*
-* Copyright (c) 2018 The Stdlib Authors.
-*
-* Licensed under the Apache License, Version 2.0 (the "License");
-* you may not use this file except in compliance with the License.
-* You may obtain a copy of the License at
-*
-*    http://www.apache.org/licenses/LICENSE-2.0
-*
-* Unless required by applicable law or agreed to in writing, software
-* distributed under the License is distributed on an "AS IS" BASIS,
-* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
-* See the License for the specific language governing permissions and
-* limitations under the License.
-*
-*
-* ## Notice
-*
-* The original C++ code and copyright notice are from the [Boost library]{@link http://www.boost.org/doc/libs/1_62_0/boost/math/special_functions/detail/t_distribution_inv.hpp}. The implementation has been modified for JavaScript.
-*
-* ```text
-* Copyright John Maddock 2006.
-* Copyright Paul A. Bristow 2007.
-*
-* Use, modification and distribution are subject to the
-* Boost Software License, Version 1.0. (See accompanying file
-* LICENSE or copy at http://www.boost.org/LICENSE_1_0.txt)
-* ```
-*/
-var f, O;
-function Lr() {
-  if (O) return f;
-  O = 1;
-  var L = g(), c = j(), m = z(), $ = G(), P = J(), _ = K(), A = M(), S = U(), B = V(), C = W(), T = Y(), I = Z(), h = d(), o = rr(), l = er(), k = 268435456, F = 1 / 3, N = 106 / 3, Q = 0.8549879733383485;
-  function X(e, r, t) {
-    var s, x, p, v, i, H, R, w, q, u, D, a, E, b, n, y;
-    if (i = 0, r > t ? (R = t, t = r, r = R, v = !0) : v = !1, c(e) === e && e < 20)
-      switch (x = m(1, N), c(e)) {
-        case 1:
-          r === 0.5 ? i = 0 : i = -S(I * r) / C(I * r);
-          break;
-        case 2:
-          i = (2 * r - 1) / _(2 * r * t);
-          break;
-        case 4:
-          H = 4 * r * t, p = _(H), E = 4 * S(P(p) / 3) / p, b = _(E - 4), i = r - 0.5 < 0 ? -b : b;
-          break;
-        case 6:
-          if (r < 1e-150)
-            return (v ? -1 : 1) * l(e, r);
-          n = 4 * (r - r * r), y = B(n, F), a = 6 * (1 + Q * (1 / y - 1));
-          do
-            q = a * a, u = q * q, D = a * u, w = a, a = 2 * (8 * n * D - 270 * q + 2187) / (5 * (4 * n * u - 216 * a - 243));
-          while (A((a - w) / a) > x);
-          a = _(a - e), i = r - 0.5 < 0 ? -a : a;
-          break;
-        default:
-          e > k ? i = L(2 * r) * T : e < 3 ? (s = 0.2742 - e * 0.0242143, r > s ? i = h(e, r) : i = o(e, r)) : (s = m(1, $(e / -0.654)), r > s ? i = l(e, r) : i = o(e, r));
-      }
-    else e > k ? i = -L(2 * r) * T : e < 3 ? (s = 0.2742 - e * 0.0242143, r > s ? i = h(e, r) : i = o(e, r)) : (s = m(1, $(e / -0.654)), r > s ? i = l(e, r) : i = o(e, r));
-    return v ? -i : i;
-  }
-  return f = X, f;
-}
-export {
-  Lr as __require
-};
-//# sourceMappingURL=inverse_students_t.js.map

package/dist/node_modules/.pnpm/@stdlib_math-base-special-kernel-betaincinv@0.2.2/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/inverse_students_t.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"inverse_students_t.js","sources":["../../../../../../../../../../node_modules/.pnpm/@stdlib+math-base-special-kernel-betaincinv@0.2.2/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/inverse_students_t.js"],"sourcesContent":["/**\n* @license Apache-2.0\n*\n* Copyright (c) 2018 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*\n*\n* ## Notice\n*\n* The original C++ code and copyright notice are from the [Boost library]{@link http://www.boost.org/doc/libs/1_62_0/boost/math/special_functions/detail/t_distribution_inv.hpp}. The implementation has been modified for JavaScript.\n*\n* ```text\n* Copyright John Maddock 2006.\n* Copyright Paul A. Bristow 2007.\n*\n* Use, modification and distribution are subject to the\n* Boost Software License, Version 1.0. (See accompanying file\n* LICENSE or copy at http://www.boost.org/LICENSE_1_0.txt)\n* ```\n*/\n\n'use strict';\n\n// MODULES //\n\nvar erfcinv = require( '@stdlib/math-base-special-erfcinv' );\nvar floor = require( '@stdlib/math-base-special-floor' );\nvar ldexp = require( '@stdlib/math-base-special-ldexp' );\nvar round = require( '@stdlib/math-base-special-round' );\nvar acos = require( '@stdlib/math-base-special-acos' );\nvar sqrt = require( '@stdlib/math-base-special-sqrt' );\nvar abs = require( '@stdlib/math-base-special-abs' );\nvar cos = require( '@stdlib/math-base-special-cos' );\nvar pow = require( '@stdlib/math-base-special-pow' );\nvar sin = require( '@stdlib/math-base-special-sin' );\nvar SQRT2 = require( '@stdlib/constants-float64-sqrt-two' );\nvar PI = require( '@stdlib/constants-float64-pi' );\nvar inverseStudentsTBodySeries = require( './inverse_students_t_body_series.js' );\nvar inverseStudentsTTailSeries = require( './inverse_students_t_tail_series.js' );\nvar inverseStudentsTHill = require( './inverse_students_t_hill.js' );\n\n\n// VARIABLES //\n\nvar DF_THRESHOLD = 0x10000000; // 2^28\nvar ONE_THIRD = 1.0 / 3.0;\nvar EXP = ( 2.0 * 53.0 ) / 3.0;\nvar C = 0.85498797333834849467655443627193;\n\n\n// MAIN //\n\n/**\n* Evaluates Student's t quantiles.\n*\n* @private\n* @param {PositiveNumber} df - degrees of freedom\n* @param {Probability} u - input probability\n* @param {Probability} v - probability equal to `1-u`\n* @returns {number} function value\n*/\nfunction inverseStudentsT( df, u, v ) {\n\tvar crossover;\n\tvar tolerance;\n\tvar rootAlpha;\n\tvar invert;\n\tvar result;\n\tvar alpha;\n\tvar tmp;\n\tvar p0;\n\tvar p2;\n\tvar p4;\n\tvar p5;\n\tvar p;\n\tvar r;\n\tvar x;\n\tvar a;\n\tvar b;\n\n\tresult = 0;\n\tif ( u > v ) {\n\t\t// Function is symmetric, so invert it:\n\t\ttmp = v;\n\t\tv = u;\n\t\tu = tmp;\n\t\tinvert = true;\n\t} else {\n\t\tinvert = false;\n\t}\n\tif ( floor(df) === df && df < 20 ) {\n\t\t// We have integer degrees of freedom, try for the special cases first:\n\t\ttolerance = ldexp( 1.0, EXP );\n\n\t\tswitch ( floor( df ) ) {\n\t\tcase 1:\n\t\t\t// `df = 1` is the same as the Cauchy distribution, see Shaw Eq 35:\n\t\t\tif ( u === 0.5 ) {\n\t\t\t\tresult = 0.0;\n\t\t\t} else {\n\t\t\t\tresult = -cos( PI * u ) / sin( PI * u );\n\t\t\t}\n\t\t\tbreak;\n\t\tcase 2:\n\t\t\t// `df = 2` has an exact result, see Shaw Eq 36:\n\t\t\tresult = ( (2.0*u) - 1.0 ) / sqrt( 2.0 * u * v );\n\t\t\tbreak;\n\t\tcase 4:\n\t\t\t// `df = 4` has an exact result, see Shaw Eq 38 & 39:\n\t\t\talpha = 4.0 * u * v;\n\t\t\trootAlpha = sqrt( alpha );\n\t\t\tr = 4 * cos( acos( rootAlpha ) / 3.0 ) / rootAlpha;\n\t\t\tx = sqrt( r - 4.0 );\n\t\t\tresult = ( u - 0.5 < 0.0 ) ? -x : x;\n\t\t\tbreak;\n\t\tcase 6:\n\t\t\t// We get numeric overflow in this area:\n\t\t\tif ( u < 1.0e-150 ) {\n\t\t\t\treturn ( ( invert ) ? -1 : 1 ) * inverseStudentsTHill( df, u );\n\t\t\t}\n\t\t\t// Newton-Raphson iteration of a polynomial case, choice of seed value is taken from Shaw's online supplement:\n\t\t\ta = 4.0 * ( u - (u*u) ); // 1 - 4 * (u - 0.5f) * (u - 0.5f);\n\t\t\tb = pow( a, ONE_THIRD );\n\t\t\tp = 6.0 * ( 1.0 + ( C * ( (1.0/b) - 1.0 ) ) );\n\t\t\tdo {\n\t\t\t\tp2 = p * p;\n\t\t\t\tp4 = p2 * p2;\n\t\t\t\tp5 = p * p4;\n\t\t\t\tp0 = p;\n\n\t\t\t\t// Next term is given by Eq 41:\n\t\t\t\tp = 2.0 * ( (8.0*a*p5) - (270.0*p2) + 2187 ) /\n\t\t\t\t\t( 5.0 * ( (4.0*a*p4) - (216.0*p) - 243.0 ) );\n\t\t\t} while ( abs( (p - p0) / p ) > tolerance );\n\n\t\t\t// Use Eq 45 to extract the result:\n\t\t\tp = sqrt( p - df );\n\t\t\tresult = ( u - 0.5 < 0.0 ) ? -p : p;\n\t\t\tbreak;\n\t\tdefault:\n\t\t\tif ( df > DF_THRESHOLD ) { // 2^28\n\t\t\t\tresult = erfcinv( 2.0 * u ) * SQRT2;\n\t\t\t} else if ( df < 3 ) {\n\t\t\t\t// Use a roughly linear scheme to choose between Shaw's tail series and body series:\n\t\t\t\tcrossover = 0.2742 - ( df * 0.0242143 );\n\t\t\t\tif ( u > crossover ) {\n\t\t\t\t\tresult = inverseStudentsTBodySeries( df, u );\n\t\t\t\t} else {\n\t\t\t\t\tresult = inverseStudentsTTailSeries( df, u );\n\t\t\t\t}\n\t\t\t} else {\n\t\t\t\t// Use Hill's method except in the extreme tails where we use Shaw's tail series. The crossover point is roughly exponential in -df:\n\t\t\t\tcrossover = ldexp( 1.0, round( df / -0.654 ) );\n\t\t\t\tif ( u > crossover ) {\n\t\t\t\t\tresult = inverseStudentsTHill( df, u );\n\t\t\t\t} else {\n\t\t\t\t\tresult = inverseStudentsTTailSeries( df, u );\n\t\t\t\t}\n\t\t\t}\n\t\t}\n\t} else if ( df > DF_THRESHOLD ) {\n\t\tresult = -erfcinv( 2.0 * u ) * SQRT2;\n\t} else if ( df < 3 ) {\n\t\t// Use a roughly linear scheme to choose between Shaw's tail series and body series:\n\t\tcrossover = 0.2742 - ( df * 0.0242143 );\n\t\tif ( u > crossover ) {\n\t\t\tresult = inverseStudentsTBodySeries( df, u );\n\t\t} else {\n\t\t\tresult = inverseStudentsTTailSeries( df, u );\n\t\t}\n\t} else {\n\t\t// Use Hill's method except in the extreme tails where we use Shaw's tail series. The crossover point is roughly exponential in -df:\n\t\tcrossover = ldexp( 1.0, round( df / -0.654 ) );\n\t\tif ( u > crossover ) {\n\t\t\tresult = inverseStudentsTHill( df, u );\n\t\t} else {\n\t\t\tresult = inverseStudentsTTailSeries( df, u );\n\t\t}\n\t}\n\treturn ( invert ) ? -result : result;\n}\n\n\n// EXPORTS //\n\nmodule.exports = inverseStudentsT;\n"],"names":["erfcinv","require$$0","floor","require$$1","ldexp","require$$2","round","require$$3","acos","require$$4","sqrt","require$$5","abs","require$$6","cos","require$$7","pow","require$$8","sin","require$$9","SQRT2","require$$10","PI","require$$11","inverseStudentsTBodySeries","require$$12","inverseStudentsTTailSeries","require$$13","inverseStudentsTHill","require$$14","DF_THRESHOLD","ONE_THIRD","EXP","C","inverseStudentsT","df","u","v","crossover","tolerance","rootAlpha","invert","result","alpha","tmp","p0","p2","p4","p5","p","r","x","a","b","inverse_students_t"],"mappings":";;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;AAoCA,MAAIA,IAAUC,EAAA,GACVC,IAAQC,EAAA,GACRC,IAAQC,EAAA,GACRC,IAAQC,EAAA,GACRC,IAAOC,EAAA,GACPC,IAAOC,EAAA,GACPC,IAAMC,EAAA,GACNC,IAAMC,EAAA,GACNC,IAAMC,EAAA,GACNC,IAAMC,EAAA,GACNC,IAAQC,EAAA,GACRC,IAAKC,EAAA,GACLC,IAA6BC,EAAA,GAC7BC,IAA6BC,GAAA,GAC7BC,IAAuBC,GAAA,GAKvBC,IAAe,WACfC,IAAY,IAAM,GAClBC,IAAQ,MAAe,GACvBC,IAAI;AAcR,WAASC,EAAkBC,GAAIC,GAAGC,GAAI;AACrC,QAAIC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC;AAYJ,QAVAX,IAAS,GACJN,IAAIC,KAERO,IAAMP,GACNA,IAAID,GACJA,IAAIQ,GACJH,IAAS,MAETA,IAAS,IAELvC,EAAMiC,CAAE,MAAMA,KAAMA,IAAK;AAI7B,cAFAI,IAAYnC,EAAO,GAAK4B,CAAG,GAElB9B,EAAOiC,CAAE,GAAE;AAAA,QACpB,KAAK;AAEJ,UAAKC,MAAM,MACVM,IAAS,IAETA,IAAS,CAAC5B,EAAKQ,IAAKc,CAAC,IAAKlB,EAAKI,IAAKc,CAAC;AAEtC;AAAA,QACD,KAAK;AAEJ,UAAAM,KAAY,IAAIN,IAAK,KAAQ1B,EAAM,IAAM0B,IAAIC,CAAC;AAC9C;AAAA,QACD,KAAK;AAEJ,UAAAM,IAAQ,IAAMP,IAAIC,GAClBG,IAAY9B,EAAMiC,CAAK,GACvBO,IAAI,IAAIpC,EAAKN,EAAMgC,CAAS,IAAK,CAAG,IAAKA,GACzCW,IAAIzC,EAAMwC,IAAI,CAAG,GACjBR,IAAWN,IAAI,MAAM,IAAQ,CAACe,IAAIA;AAClC;AAAA,QACD,KAAK;AAEJ,cAAKf,IAAI;AACR,oBAAWK,IAAW,KAAK,KAAMb,EAAsBO,GAAIC,CAAC;AAG7D,UAAAgB,IAAI,KAAQhB,IAAKA,IAAEA,IACnBiB,IAAIrC,EAAKoC,GAAGrB,CAAS,GACrBkB,IAAI,KAAQ,IAAQhB,KAAO,IAAIoB,IAAK;AACpC;AACC,YAAAP,IAAKG,IAAIA,GACTF,IAAKD,IAAKA,GACVE,IAAKC,IAAIF,GACTF,IAAKI,GAGLA,IAAI,KAAS,IAAIG,IAAEJ,IAAO,MAAMF,IAAM,SACnC,KAAS,IAAIM,IAAEL,IAAO,MAAME,IAAK;AAAA,iBAC3BrC,GAAMqC,IAAIJ,KAAMI,CAAC,IAAKV;AAGhC,UAAAU,IAAIvC,EAAMuC,IAAId,CAAE,GAChBO,IAAWN,IAAI,MAAM,IAAQ,CAACa,IAAIA;AAClC;AAAA,QACD;AACC,UAAKd,IAAKL,IACTY,IAAS1C,EAAS,IAAMoC,CAAC,IAAKhB,IACnBe,IAAK,KAEhBG,IAAY,SAAWH,IAAK,WACvBC,IAAIE,IACRI,IAASlB,EAA4BW,GAAIC,CAAC,IAE1CM,IAAShB,EAA4BS,GAAIC,CAAC,MAI3CE,IAAYlC,EAAO,GAAKE,EAAO6B,IAAK,MAAM,CAAE,GACvCC,IAAIE,IACRI,IAASd,EAAsBO,GAAIC,CAAC,IAEpCM,IAAShB,EAA4BS,GAAIC,CAAC;AAAA,MAG/C;AAAA,QACQ,CAAKD,IAAKL,IAChBY,IAAS,CAAC1C,EAAS,IAAMoC,CAAC,IAAKhB,IACpBe,IAAK,KAEhBG,IAAY,SAAWH,IAAK,WACvBC,IAAIE,IACRI,IAASlB,EAA4BW,GAAIC,CAAC,IAE1CM,IAAShB,EAA4BS,GAAIC,CAAC,MAI3CE,IAAYlC,EAAO,GAAKE,EAAO6B,IAAK,MAAM,CAAE,GACvCC,IAAIE,IACRI,IAASd,EAAsBO,GAAIC,CAAC,IAEpCM,IAAShB,EAA4BS,GAAIC,CAAC;AAG5C,WAASK,IAAW,CAACC,IAASA;AAAA,EAC/B;AAKA,SAAAY,IAAiBpB;;","x_google_ignoreList":[0]}

package/dist/node_modules/.pnpm/@stdlib_math-base-special-kernel-betaincinv@0.2.2/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/inverse_students_t_body_series.js

@@ -1,59 +0,0 @@
-import { __require as I } from "../../../../../@stdlib_math-base-special-gamma-delta-ratio@0.2.2/node_modules/@stdlib/math-base-special-gamma-delta-ratio/lib/index.js";
-import { __require as $ } from "../../../../../@stdlib_math-base-tools-evalpoly@0.2.2/node_modules/@stdlib/math-base-tools-evalpoly/lib/index.js";
-import { __require as R } from "../../../../../@stdlib_math-base-special-sqrt@0.2.2/node_modules/@stdlib/math-base-special-sqrt/lib/index.js";
-import { __require as S } from "../../../../../@stdlib_constants-float64-pi@0.2.2/node_modules/@stdlib/constants-float64-pi/lib/index.js";
-import { __require as g } from "./polyval_co14.js";
-import { __require as h } from "./polyval_co15.js";
-import { __require as x } from "./polyval_co16.js";
-import { __require as B } from "./polyval_co17.js";
-import { __require as D } from "./polyval_co18.js";
-import { __require as T } from "./polyval_co19.js";
-import { __require as j } from "./polyval_co20.js";
-import { __require as k } from "./polyval_co21.js";
-import { __require as w } from "./polyval_co22.js";
-/**
-* @license Apache-2.0
-*
-* Copyright (c) 2018 The Stdlib Authors.
-*
-* Licensed under the Apache License, Version 2.0 (the "License");
-* you may not use this file except in compliance with the License.
-* You may obtain a copy of the License at
-*
-*    http://www.apache.org/licenses/LICENSE-2.0
-*
-* Unless required by applicable law or agreed to in writing, software
-* distributed under the License is distributed on an "AS IS" BASIS,
-* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
-* See the License for the specific language governing permissions and
-* limitations under the License.
-*
-*
-* ## Notice
-*
-* The original C++ code and copyright notice are from the [Boost library]{@link http://www.boost.org/doc/libs/1_62_0/boost/math/special_functions/detail/t_distribution_inv.hpp}. The implementation has been modified for JavaScript.
-*
-* ```text
-* Copyright John Maddock 2006.
-* Copyright Paul A. Bristow 2007.
-*
-* Use, modification and distribution are subject to the
-* Boost Software License, Version 1.0. (See accompanying file
-* LICENSE or copy at http://www.boost.org/LICENSE_1_0.txt)
-* ```
-*/
-var i, _;
-function U() {
-  if (_) return i;
-  _ = 1;
-  var v = I(), l = $(), u = R(), s = S(), t = g(), q = h(), m = x(), p = B(), y = D(), f = T(), n = j(), c = k(), P = w(), d = 0, r = [1, 0, 0, 0, 0, 0, 0, 0, 0, 0];
-  function b(a, L) {
-    var e, o;
-    return o = v(a / 2, 0.5) * u(a * s) * (L - 0.5), e = 1 / a, r[1] = t(e), r[2] = q(e), r[3] = m(e), r[4] = p(e), r[5] = y(e), r[6] = f(e), r[7] = n(e), r[8] = c(e), r[9] = P(e), d + o * l(r, o * o);
-  }
-  return i = b, i;
-}
-export {
-  U as __require
-};
-//# sourceMappingURL=inverse_students_t_body_series.js.map

package/dist/node_modules/.pnpm/@stdlib_math-base-special-kernel-betaincinv@0.2.2/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/inverse_students_t_body_series.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"inverse_students_t_body_series.js","sources":["../../../../../../../../../../node_modules/.pnpm/@stdlib+math-base-special-kernel-betaincinv@0.2.2/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/inverse_students_t_body_series.js"],"sourcesContent":["/**\n* @license Apache-2.0\n*\n* Copyright (c) 2018 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*\n*\n* ## Notice\n*\n* The original C++ code and copyright notice are from the [Boost library]{@link http://www.boost.org/doc/libs/1_62_0/boost/math/special_functions/detail/t_distribution_inv.hpp}. The implementation has been modified for JavaScript.\n*\n* ```text\n* Copyright John Maddock 2006.\n* Copyright Paul A. Bristow 2007.\n*\n* Use, modification and distribution are subject to the\n* Boost Software License, Version 1.0. (See accompanying file\n* LICENSE or copy at http://www.boost.org/LICENSE_1_0.txt)\n* ```\n*/\n\n'use strict';\n\n// MODULES //\n\nvar gammaDeltaRatio = require( '@stdlib/math-base-special-gamma-delta-ratio' );\nvar evalpoly = require( '@stdlib/math-base-tools-evalpoly' );\nvar sqrt = require( '@stdlib/math-base-special-sqrt' );\nvar PI = require( '@stdlib/constants-float64-pi' );\nvar polyval1 = require( './polyval_co14.js' );\nvar polyval2 = require( './polyval_co15.js' );\nvar polyval3 = require( './polyval_co16.js' );\nvar polyval4 = require( './polyval_co17.js' );\nvar polyval5 = require( './polyval_co18.js' );\nvar polyval6 = require( './polyval_co19.js' );\nvar polyval7 = require( './polyval_co20.js' );\nvar polyval8 = require( './polyval_co21.js' );\nvar polyval9 = require( './polyval_co22.js' );\n\n\n// VARIABLES //\n\nvar c0 = 0.0;\n\n// Workspace for the polynomial coefficients:\nvar c = [ 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ]; // WARNING: not thread safe\n\n\n// MAIN //\n\n/**\n* Evaluates Student's t quantiles via a body series expansion. Tail and body series are due to Shaw.\n*\n* ## References\n*\n* -   Shaw, William T. 2006. \"Sampling Student's T distribution – use of the inverse cumulative distribution function.\" _Journal of Computational Finance_ 9 (4): 37–73. [www.mth.kcl.ac.uk/~shaww/web\\_page/papers/Tdistribution06.pdf](www.mth.kcl.ac.uk/~shaww/web_page/papers/Tdistribution06.pdf).\n*\n* @private\n* @param {PositiveNumber} df - degrees of freedom\n* @param {Probability} u - input probability\n* @returns {number} function value\n*/\nfunction inverseStudentsTBodySeries( df, u ) {\n\tvar idf;\n\tvar v;\n\n\t// Body series for small N, start with Eq 56 of Shaw:\n\tv = gammaDeltaRatio( df/2, 0.5 ) * sqrt( df*PI ) * ( u-0.5 );\n\n\t// Figure out what the coefficients are. They depend only on the degrees of freedom (Eq 57 of Shaw):\n\tidf = 1.0 / df;\n\tc[ 1 ] = polyval1( idf );\n\tc[ 2 ] = polyval2( idf );\n\tc[ 3 ] = polyval3( idf );\n\tc[ 4 ] = polyval4( idf );\n\tc[ 5 ] = polyval5( idf );\n\tc[ 6 ] = polyval6( idf );\n\tc[ 7 ] = polyval7( idf );\n\tc[ 8 ] = polyval8( idf );\n\tc[ 9 ] = polyval9( idf );\n\n\t// Result is then an odd polynomial in v (see Eq 56 of Shaw)...\n\treturn c0 + ( v*evalpoly( c, v*v ) );\n}\n\n\n// EXPORTS //\n\nmodule.exports = inverseStudentsTBodySeries;\n"],"names":["gammaDeltaRatio","require$$0","evalpoly","require$$1","sqrt","require$$2","PI","require$$3","polyval1","require$$4","polyval2","require$$5","polyval3","require$$6","polyval4","require$$7","polyval5","require$$8","polyval6","require$$9","polyval7","require$$10","polyval8","require$$11","polyval9","require$$12","c0","c","inverseStudentsTBodySeries","df","u","idf","v","inverse_students_t_body_series"],"mappings":";;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;AAoCA,MAAIA,IAAkBC,EAAA,GAClBC,IAAWC,EAAA,GACXC,IAAOC,EAAA,GACPC,IAAKC,EAAA,GACLC,IAAWC,EAAA,GACXC,IAAWC,EAAA,GACXC,IAAWC,EAAA,GACXC,IAAWC,EAAA,GACXC,IAAWC,EAAA,GACXC,IAAWC,EAAA,GACXC,IAAWC,EAAA,GACXC,IAAWC,EAAA,GACXC,IAAWC,EAAA,GAKXC,IAAK,GAGLC,IAAI,CAAE,GAAK,GAAK,GAAK,GAAK,GAAK,GAAK,GAAK,GAAK,GAAK,CAAG;AAiB1D,WAASC,EAA4BC,GAAIC,GAAI;AAC5C,QAAIC,GACAC;AAGJ,WAAAA,IAAIhC,EAAiB6B,IAAG,GAAG,GAAG,IAAKzB,EAAMyB,IAAGvB,MAASwB,IAAE,MAGvDC,IAAM,IAAMF,GACZF,EAAG,CAAC,IAAKnB,EAAUuB,CAAG,GACtBJ,EAAG,CAAC,IAAKjB,EAAUqB,CAAG,GACtBJ,EAAG,CAAC,IAAKf,EAAUmB,CAAG,GACtBJ,EAAG,CAAC,IAAKb,EAAUiB,CAAG,GACtBJ,EAAG,CAAC,IAAKX,EAAUe,CAAG,GACtBJ,EAAG,CAAC,IAAKT,EAAUa,CAAG,GACtBJ,EAAG,CAAC,IAAKP,EAAUW,CAAG,GACtBJ,EAAG,CAAC,IAAKL,EAAUS,CAAG,GACtBJ,EAAG,CAAC,IAAKH,EAAUO,CAAG,GAGfL,IAAOM,IAAE9B,EAAUyB,GAAGK,IAAEA;EAChC;AAKA,SAAAC,IAAiBL;;","x_google_ignoreList":[0]}

package/dist/node_modules/.pnpm/@stdlib_math-base-special-kernel-betaincinv@0.2.2/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/inverse_students_t_hill.js

@@ -1,52 +0,0 @@
-import { __require as c } from "../../../../../@stdlib_math-base-special-erfcinv@0.2.3/node_modules/@stdlib/math-base-special-erfcinv/lib/index.js";
-import { __require as h } from "../../../../../@stdlib_math-base-special-expm1@0.2.3/node_modules/@stdlib/math-base-special-expm1/lib/index.js";
-import { __require as x } from "../../../../../@stdlib_math-base-special-sqrt@0.2.2/node_modules/@stdlib/math-base-special-sqrt/lib/index.js";
-import { __require as I } from "../../../../../@stdlib_math-base-special-pow@0.3.0/node_modules/@stdlib/math-base-special-pow/lib/index.js";
-import { __require as H } from "../../../../../@stdlib_constants-float64-half-pi@0.2.2/node_modules/@stdlib/constants-float64-half-pi/lib/index.js";
-import { __require as R } from "../../../../../@stdlib_constants-float64-sqrt-two@0.2.2/node_modules/@stdlib/constants-float64-sqrt-two/lib/index.js";
-/**
-* @license Apache-2.0
-*
-* Copyright (c) 2018 The Stdlib Authors.
-*
-* Licensed under the Apache License, Version 2.0 (the "License");
-* you may not use this file except in compliance with the License.
-* You may obtain a copy of the License at
-*
-*    http://www.apache.org/licenses/LICENSE-2.0
-*
-* Unless required by applicable law or agreed to in writing, software
-* distributed under the License is distributed on an "AS IS" BASIS,
-* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
-* See the License for the specific language governing permissions and
-* limitations under the License.
-*
-*
-* ## Notice
-*
-* The original C++ code and copyright notice are from the [Boost library]{@link http://www.boost.org/doc/libs/1_61_0/boost/math/special_functions/detail/t_distribution_inv.hpp}. The implementation has been modified for JavaScript.
-*
-* ```text
-* Copyright John Maddock 2006.
-* Copyright Paul A. Bristow 2007.
-*
-* Use, modification and distribution are subject to the
-* Boost Software License, Version 1.0. (See accompanying file
-* LICENSE or copy at http://www.boost.org/LICENSE_1_0.txt)
-* ```
-*/
-var v, p;
-function F() {
-  if (p) return v;
-  p = 1;
-  var q = c(), b = h(), o = x(), L = I(), n = H(), m = R();
-  function $(e, s) {
-    var i, a, _, u, l, t, r;
-    return e > 1e20 ? -q(2 * s) * m : (i = 1 / (e - 0.5), a = 48 / (i * i), _ = ((20700 * i / a - 98) * i - 16) * i + 96.36, u = ((94.5 / (a + _) - 3) / a + 1) * o(i * n) * e, r = L(u * 2 * s, 2 / e), r > 0.05 + i ? (t = -q(2 * s) * m, r = t * t, e < 5 && (_ += 0.3 * (e - 4.5) * (t + 0.6)), _ += (((0.05 * u * t - 5) * t - 7) * t - 2) * t + a, r = (((((0.4 * r + 6.3) * r + 36) * r + 94.5) / _ - r - 3) / a + 1) * t, r = b(i * r * r)) : r = ((1 / (((e + 6) / (e * r) - 0.089 * u - 0.822) * (e + 2) * 3) + 0.5 / (e + 4)) * r - 1) * (e + 1) / (e + 2) + 1 / r, l = o(e * r), -l);
-  }
-  return v = $, v;
-}
-export {
-  F as __require
-};
-//# sourceMappingURL=inverse_students_t_hill.js.map

package/dist/node_modules/.pnpm/@stdlib_math-base-special-kernel-betaincinv@0.2.2/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/inverse_students_t_hill.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"inverse_students_t_hill.js","sources":["../../../../../../../../../../node_modules/.pnpm/@stdlib+math-base-special-kernel-betaincinv@0.2.2/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/inverse_students_t_hill.js"],"sourcesContent":["/**\n* @license Apache-2.0\n*\n* Copyright (c) 2018 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*\n*\n* ## Notice\n*\n* The original C++ code and copyright notice are from the [Boost library]{@link http://www.boost.org/doc/libs/1_61_0/boost/math/special_functions/detail/t_distribution_inv.hpp}. The implementation has been modified for JavaScript.\n*\n* ```text\n* Copyright John Maddock 2006.\n* Copyright Paul A. Bristow 2007.\n*\n* Use, modification and distribution are subject to the\n* Boost Software License, Version 1.0. (See accompanying file\n* LICENSE or copy at http://www.boost.org/LICENSE_1_0.txt)\n* ```\n*/\n\n/* eslint-disable no-mixed-operators */\n\n'use strict';\n\n// MODULES //\n\nvar erfcinv = require( '@stdlib/math-base-special-erfcinv' );\nvar expm1 = require( '@stdlib/math-base-special-expm1' );\nvar sqrt = require( '@stdlib/math-base-special-sqrt' );\nvar pow = require( '@stdlib/math-base-special-pow' );\nvar HALF_PI = require( '@stdlib/constants-float64-half-pi' );\nvar SQRT2 = require( '@stdlib/constants-float64-sqrt-two' );\n\n\n// MAIN //\n\n/**\n* Evaluates Student's t quantiles via a method due to Hill.\n*\n* ## References\n*\n* -   Hill, G. W. 1970. \"Algorithm 396: Student's T-Quantiles.\" _Communications of the ACM_ 13 (10). New York, NY, USA: ACM: 619–20. doi:[10.1145/355598.355600](https://doi.org/10.1145/355598.355600).\n*\n* @private\n* @param {PositiveNumber} ndf - degrees of freedom\n* @param {Probability} u - input probability\n* @returns {number} function value\n*/\nfunction inverseStudentsTHill( ndf, u ) {\n\tvar a;\n\tvar b;\n\tvar c;\n\tvar d;\n\tvar q;\n\tvar x;\n\tvar y;\n\n\tif ( ndf > 1e20 ) {\n\t\treturn -erfcinv( 2 * u ) * SQRT2;\n\t}\n\ta = 1.0 / ( ndf - 0.5 );\n\tb = 48.0 / (a * a);\n\tc = ( ( ( ( (20700.0*a/b) - 98.0 ) * a ) - 16.0 ) * a ) + 96.36;\n\td = ( ( ( (94.5/(b+c)) - 3.0 ) / b ) + 1.0 ) * sqrt( a * HALF_PI ) * ndf;\n\ty = pow( d * 2.0 * u, 2.0 / ndf );\n\n\tif ( y > ( 0.05 + a ) ) {\n\t\t// Asymptotic inverse expansion about normal:\n\t\tx = -erfcinv( 2.0 * u ) * SQRT2;\n\t\ty = x * x;\n\n\t\tif ( ndf < 5.0 ) {\n\t\t\tc += 0.3 * ( ndf-4.5 ) * ( x + 0.6 );\n\t\t}\n\t\tc += ( ( ( ( ( ( (0.05*d*x)-5.0 ) * x ) - 7.0 ) * x )- 2.0 ) * x ) + b;\n\t\ty = ((((((0.4*y+6.3)*y)+36.0) * y + 94.5) / c - y - 3.0) / b + 1.0) * x;\n\t\ty = expm1( a * y * y );\n\t} else {\n\t\ty = ((1.0 / ( ( (ndf+6.0) / (ndf*y) - 0.089 * d - 0.822 ) *\n\t\t(ndf+2.0) * 3.0 ) + 0.5 / (ndf+4.0)) * y - 1.0) *\n\t\t(ndf+1.0) / (ndf+2.0) + 1.0 / y;\n\t}\n\tq = sqrt( ndf * y );\n\treturn -q;\n}\n\n\n// EXPORTS //\n\nmodule.exports = inverseStudentsTHill;\n"],"names":["erfcinv","require$$0","expm1","require$$1","sqrt","require$$2","pow","require$$3","HALF_PI","require$$4","SQRT2","require$$5","inverseStudentsTHill","ndf","u","a","b","c","d","q","x","y","inverse_students_t_hill"],"mappings":";;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;AAsCA,MAAIA,IAAUC,EAAA,GACVC,IAAQC,EAAA,GACRC,IAAOC,EAAA,GACPC,IAAMC,EAAA,GACNC,IAAUC,EAAA,GACVC,IAAQC,EAAA;AAiBZ,WAASC,EAAsBC,GAAKC,GAAI;AACvC,QAAIC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC;AAEJ,WAAKR,IAAM,OACH,CAACb,EAAS,IAAIc,CAAC,IAAKJ,KAE5BK,IAAI,KAAQF,IAAM,MAClBG,IAAI,MAAQD,IAAIA,IAChBE,MAAa,QAAQF,IAAEC,IAAK,MAASD,IAAM,MAASA,IAAM,OAC1DG,MAAW,QAAMF,IAAEC,KAAM,KAAQD,IAAM,KAAQZ,EAAMW,IAAIP,CAAO,IAAKK,GACrEQ,IAAIf,EAAKY,IAAI,IAAMJ,GAAG,IAAMD,CAAG,GAE1BQ,IAAM,OAAON,KAEjBK,IAAI,CAACpB,EAAS,IAAMc,CAAC,IAAKJ,GAC1BW,IAAID,IAAIA,GAEHP,IAAM,MACVI,KAAK,OAAQJ,IAAI,QAAUO,IAAI,OAEhCH,QAAkB,OAAKC,IAAEE,IAAG,KAAQA,IAAM,KAAQA,IAAK,KAAQA,IAAMJ,GACrEK,SAAU,MAAIA,IAAE,OAAKA,IAAG,MAAQA,IAAI,QAAQJ,IAAII,IAAI,KAAOL,IAAI,KAAOI,GACtEC,IAAInB,EAAOa,IAAIM,IAAIA,CAAC,KAEpBA,MAAM,OAAWR,IAAI,MAAQA,IAAIQ,KAAK,QAAQH,IAAI,UACjDL,IAAI,KAAO,KAAQ,OAAOA,IAAI,MAAQQ,IAAI,MAC1CR,IAAI,MAAQA,IAAI,KAAO,IAAMQ,GAE/BF,IAAIf,EAAMS,IAAMQ,CAAC,GACV,CAACF;AAAA,EACT;AAKA,SAAAG,IAAiBV;;","x_google_ignoreList":[0]}

package/dist/node_modules/.pnpm/@stdlib_math-base-special-kernel-betaincinv@0.2.2/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/inverse_students_t_tail_series.js

@@ -1,51 +0,0 @@
-import { __require as I } from "../../../../../@stdlib_math-base-special-gamma-delta-ratio@0.2.2/node_modules/@stdlib/math-base-special-gamma-delta-ratio/lib/index.js";
-import { __require as c } from "../../../../../@stdlib_math-base-tools-evalpoly@0.2.2/node_modules/@stdlib/math-base-tools-evalpoly/lib/index.js";
-import { __require as R } from "../../../../../@stdlib_math-base-special-sqrt@0.2.2/node_modules/@stdlib/math-base-special-sqrt/lib/index.js";
-import { __require as S } from "../../../../../@stdlib_math-base-special-pow@0.3.0/node_modules/@stdlib/math-base-special-pow/lib/index.js";
-import { __require as T } from "../../../../../@stdlib_constants-float64-pi@0.2.2/node_modules/@stdlib/constants-float64-pi/lib/index.js";
-/**
-* @license Apache-2.0
-*
-* Copyright (c) 2018 The Stdlib Authors.
-*
-* Licensed under the Apache License, Version 2.0 (the "License");
-* you may not use this file except in compliance with the License.
-* You may obtain a copy of the License at
-*
-*    http://www.apache.org/licenses/LICENSE-2.0
-*
-* Unless required by applicable law or agreed to in writing, software
-* distributed under the License is distributed on an "AS IS" BASIS,
-* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
-* See the License for the specific language governing permissions and
-* limitations under the License.
-*
-*
-* ## Notice
-*
-* The original C++ code and copyright notice are from the [Boost library]{@link http://www.boost.org/doc/libs/1_61_0/boost/math/special_functions/detail/t_distribution_inv.hpp}. The implementation has been modified for JavaScript.
-*
-* ```text
-* Copyright John Maddock 2006.
-* Copyright Paul A. Bristow 2007.
-*
-* Use, modification and distribution are subject to the
-* Boost Software License, Version 1.0. (See accompanying file
-* LICENSE or copy at http://www.boost.org/LICENSE_1_0.txt)
-* ```
-*/
-var v, m;
-function P() {
-  if (m) return v;
-  m = 1;
-  var p = I(), l = c(), n = R(), b = S(), L = T(), i = [0, 0, 0, 0, 0, 0, 0];
-  function $(r, w) {
-    var s, o, _, e, a, t, u, q;
-    return q = p(r / 2, 0.5) * n(r * L) * w, e = r + 2, a = r + 4, t = r + 6, i[0] = 1, i[1] = -(r + 1) / (2 * e), e *= r + 2, i[2] = -r * (r + 1) * (r + 3) / (8 * e * a), e *= r + 2, i[3] = -r * (r + 1) * (r + 5) * ((3 * r + 7) * r - 2) / (48 * e * a * t), e *= r + 2, a *= r + 4, i[4] = -r * (r + 1) * (r + 7) * (((((15 * r + 154) * r + 465) * r + 286) * r - 336) * r + 64) / (384 * e * a * t * (r + 8)), e *= r + 2, i[5] = -r * (r + 1) * (r + 3) * (r + 9) * ((((((35 * r + 452) * r + 1573) * r + 600) * r - 2020) * r + 928) * r - 128) / (1280 * e * a * t * (r + 8) * (r + 10)), e *= r + 2, a *= r + 4, t *= r + 6, i[6] = -r * (r + 1) * (r + 11) * (((((((((((945 * r + 31506) * r + 425858) * r + 2980236) * r + 11266745) * r + 20675018) * r + 7747124) * r - 22574632) * r - 8565600) * r + 18108416) * r - 7099392) * r + 884736) / (46080 * e * a * t * (r + 8) * (r + 10) * (r + 12)), u = n(r), _ = b(u * q, 1 / r), o = _ * _, s = l(i, o), s *= u, s /= _, -s;
-  }
-  return v = $, v;
-}
-export {
-  P as __require
-};
-//# sourceMappingURL=inverse_students_t_tail_series.js.map

package/dist/node_modules/.pnpm/@stdlib_math-base-special-kernel-betaincinv@0.2.2/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/inverse_students_t_tail_series.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"inverse_students_t_tail_series.js","sources":["../../../../../../../../../../node_modules/.pnpm/@stdlib+math-base-special-kernel-betaincinv@0.2.2/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/inverse_students_t_tail_series.js"],"sourcesContent":["/**\n* @license Apache-2.0\n*\n* Copyright (c) 2018 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*\n*\n* ## Notice\n*\n* The original C++ code and copyright notice are from the [Boost library]{@link http://www.boost.org/doc/libs/1_61_0/boost/math/special_functions/detail/t_distribution_inv.hpp}. The implementation has been modified for JavaScript.\n*\n* ```text\n* Copyright John Maddock 2006.\n* Copyright Paul A. Bristow 2007.\n*\n* Use, modification and distribution are subject to the\n* Boost Software License, Version 1.0. (See accompanying file\n* LICENSE or copy at http://www.boost.org/LICENSE_1_0.txt)\n* ```\n*/\n\n'use strict';\n\n// MODULES //\n\nvar gammaDeltaRatio = require( '@stdlib/math-base-special-gamma-delta-ratio' );\nvar evalpoly = require( '@stdlib/math-base-tools-evalpoly' );\nvar sqrt = require( '@stdlib/math-base-special-sqrt' );\nvar pow = require( '@stdlib/math-base-special-pow' );\nvar PI = require( '@stdlib/constants-float64-pi' );\n\n\n// VARIABLES //\n\n// Array for the coefficients d(k), these depend only on the number of degrees of freedom df, so at least in theory we could tabulate these for fixed df, see p15 of Shaw:\nvar d = [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ]; // WARNING: not thread safe\n\n\n// MAIN //\n\n/**\n* Evaluates Student's t quantiles via a tail series expansion. Tail and body series are due to Shaw.\n*\n* ## References\n*\n* -   Shaw, William T. 2006. \"Sampling Student's T distribution – use of the inverse cumulative distribution function.\" _Journal of Computational Finance_ 9 (4): 37–73. [www.mth.kcl.ac.uk/~shaww/web\\_page/papers/Tdistribution06.pdf](www.mth.kcl.ac.uk/~shaww/web_page/papers/Tdistribution06.pdf).\n*\n* @private\n* @param {number} df - degrees of freedom\n* @param {number} v - function value\n* @returns {number} tail value\n*/\nfunction inverseStudentsTTailSeries( df, v ) { // eslint-disable-line id-length\n\tvar result;\n\tvar power;\n\tvar div;\n\tvar np2;\n\tvar np4;\n\tvar np6;\n\tvar rn;\n\tvar w;\n\n\t// Tail series expansion, see section 6 of Shaw's paper. `w` is calculated using Eq 60:\n\tw = gammaDeltaRatio( df/2.0, 0.5 ) * sqrt( df*PI ) * v;\n\n\t// Define some variables:\n\tnp2 = df + 2.0;\n\tnp4 = df + 4.0;\n\tnp6 = df + 6.0;\n\n\td[ 0 ] = 1.0;\n\td[ 1 ] = -(df+1.0) / (2.0*np2);\n\tnp2 *= (df + 2.0);\n\td[ 2 ] = -df * (df+1.0) * (df+3.0) / (8.0*np2*np4);\n\tnp2 *= df + 2.0;\n\td[ 3 ] = -df * (df+1.0) * (df+5.0) * (((3.0*df) + 7.0) * df - 2.0) / (48.0*np2*np4*np6); // eslint-disable-line max-len, no-mixed-operators\n\tnp2 *= (df + 2.0);\n\tnp4 *= (df + 4.0);\n\td[ 4 ] = -df * (df+1.0) * (df+7.0) * ( (((((15.0*df) + 154.0) * df + 465.0) * df + 286.0) * df - 336.0) * df + 64.0) / (384.0*np2*np4*np6*(df+8.0)); // eslint-disable-line max-len, no-mixed-operators\n\tnp2 *= (df + 2.0);\n\td[ 5 ] = -df * (df+1.0) * (df+3.0) * (df+9.0) * (((((((35.0 * df + 452.0) * df+1573.0) * df + 600.0) * df - 2020.0) * df) + 928.0) * df - 128.0) / (1280.0*np2*np4*np6*(df+8.0) * (df+10.0)); // eslint-disable-line max-len, no-mixed-operators\n\tnp2 *= (df + 2.0);\n\tnp4 *= (df + 4.0);\n\tnp6 *= (df + 6.0);\n\td[ 6 ] = -df * (df+1.0) * (df+11.0) * ((((((((((((945.0*df) + 31506.0) * df + 425858.0) * df + 2980236.0) * df + 11266745.0) * df + 20675018.0) * df + 7747124.0) * df - 22574632.0) * df - 8565600.0) * df + 18108416.0) * df - 7099392.0) * df + 884736.0) / (46080.0*np2*np4*np6*(df+8.0) * (df+10.0) * (df+12.0)); // eslint-disable-line max-len, no-mixed-operators\n\n\t// Now bring everything together to provide the result this is Eq 62 of Shaw:\n\trn = sqrt( df );\n\tdiv = pow( rn*w, 1.0/df );\n\tpower = div * div;\n\tresult = evalpoly( d, power );\n\tresult *= rn;\n\tresult /= div;\n\treturn -result;\n}\n\n\n// EXPORTS //\n\nmodule.exports = inverseStudentsTTailSeries;\n"],"names":["gammaDeltaRatio","require$$0","evalpoly","require$$1","sqrt","require$$2","pow","require$$3","PI","require$$4","d","inverseStudentsTTailSeries","df","v","result","power","div","np2","np4","np6","rn","w","inverse_students_t_tail_series"],"mappings":";;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;AAoCA,MAAIA,IAAkBC,EAAA,GAClBC,IAAWC,EAAA,GACXC,IAAOC,EAAA,GACPC,IAAMC,EAAA,GACNC,IAAKC,EAAA,GAMLC,IAAI,CAAE,GAAK,GAAK,GAAK,GAAK,GAAK,GAAK;AAiBxC,WAASC,EAA4BC,GAAIC,GAAI;AAC5C,QAAIC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC;AAGJ,WAAAA,IAAIrB,EAAiBY,IAAG,GAAK,GAAG,IAAKR,EAAMQ,IAAGJ,CAAE,IAAKK,GAGrDI,IAAML,IAAK,GACXM,IAAMN,IAAK,GACXO,IAAMP,IAAK,GAEXF,EAAG,CAAC,IAAK,GACTA,EAAG,CAAC,IAAK,EAAEE,IAAG,MAAQ,IAAIK,IAC1BA,KAAQL,IAAK,GACbF,EAAG,CAAC,IAAK,CAACE,KAAMA,IAAG,MAAQA,IAAG,MAAQ,IAAIK,IAAIC,IAC9CD,KAAOL,IAAK,GACZF,EAAG,CAAC,IAAK,CAACE,KAAMA,IAAG,MAAQA,IAAG,OAAU,IAAIA,IAAM,KAAOA,IAAK,MAAQ,KAAKK,IAAIC,IAAIC,IACnFF,KAAQL,IAAK,GACbM,KAAQN,IAAK,GACbF,EAAG,CAAC,IAAK,CAACE,KAAMA,IAAG,MAAQA,IAAG,UAAc,KAAKA,IAAM,OAASA,IAAK,OAASA,IAAK,OAASA,IAAK,OAASA,IAAK,OAAS,MAAMK,IAAIC,IAAIC,KAAKP,IAAG,KAC9IK,KAAQL,IAAK,GACbF,EAAG,CAAC,IAAK,CAACE,KAAMA,IAAG,MAAQA,IAAG,MAAQA,IAAG,WAAc,KAAOA,IAAK,OAASA,IAAG,QAAUA,IAAK,OAASA,IAAK,QAAUA,IAAM,OAASA,IAAK,QAAU,OAAOK,IAAIC,IAAIC,KAAKP,IAAG,MAAQA,IAAG,MACtLK,KAAQL,IAAK,GACbM,KAAQN,IAAK,GACbO,KAAQP,IAAK,GACbF,EAAG,CAAC,IAAK,CAACE,KAAMA,IAAG,MAAQA,IAAG,iBAAoB,MAAMA,IAAM,SAAWA,IAAK,UAAYA,IAAK,WAAaA,IAAK,YAAcA,IAAK,YAAcA,IAAK,WAAaA,IAAK,YAAcA,IAAK,WAAaA,IAAK,YAAcA,IAAK,WAAaA,IAAK,WAAa,QAAQK,IAAIC,IAAIC,KAAKP,IAAG,MAAQA,IAAG,OAASA,IAAG,MAG/SQ,IAAKhB,EAAMQ,CAAE,GACbI,IAAMV,EAAKc,IAAGC,GAAG,IAAIT,CAAE,GACvBG,IAAQC,IAAMA,GACdF,IAASZ,EAAUQ,GAAGK,CAAK,GAC3BD,KAAUM,GACVN,KAAUE,GACH,CAACF;AAAA,EACT;AAKA,SAAAQ,IAAiBX;;","x_google_ignoreList":[0]}

package/dist/node_modules/.pnpm/@stdlib_math-base-special-kernel-betaincinv@0.2.2/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/main.js

@@ -1,65 +0,0 @@
-import { __require as or } from "../../../../../@stdlib_math-base-tools-evalpoly@0.2.2/node_modules/@stdlib/math-base-tools-evalpoly/lib/index.js";
-import { __require as _r } from "../../../../../@stdlib_math-base-special-betainc@0.2.2/node_modules/@stdlib/math-base-special-betainc/lib/index.js";
-import { __require as sr } from "../../../../../@stdlib_math-base-special-expm1@0.2.3/node_modules/@stdlib/math-base-special-expm1/lib/index.js";
-import { __require as ar } from "../../../../../@stdlib_math-base-special-log1p@0.2.3/node_modules/@stdlib/math-base-special-log1p/lib/index.js";
-import { __require as nr } from "../../../../../@stdlib_math-base-special-asin@0.2.3/node_modules/@stdlib/math-base-special-asin/lib/index.js";
-import { __require as lr } from "../../../../../@stdlib_math-base-special-beta@0.3.0/node_modules/@stdlib/math-base-special-beta/lib/index.js";
-import { __require as qr } from "../../../../../@stdlib_math-base-special-sqrt@0.2.2/node_modules/@stdlib/math-base-special-sqrt/lib/index.js";
-import { __require as Lr } from "../../../../../@stdlib_math-base-special-abs@0.2.2/node_modules/@stdlib/math-base-special-abs/lib/index.js";
-import { __require as $r } from "../../../../../@stdlib_math-base-special-exp@0.2.4/node_modules/@stdlib/math-base-special-exp/lib/index.js";
-import { __require as Ir } from "../../../../../@stdlib_math-base-special-pow@0.3.0/node_modules/@stdlib/math-base-special-pow/lib/index.js";
-import { __require as xr } from "../../../../../@stdlib_math-base-special-sin@0.3.0/node_modules/@stdlib/math-base-special-sin/lib/index.js";
-import { __require as Tr } from "../../../../../@stdlib_math-base-special-max@0.3.0/node_modules/@stdlib/math-base-special-max/lib/index.js";
-import { __require as dr } from "../../../../../@stdlib_math-base-special-min@0.2.3/node_modules/@stdlib/math-base-special-min/lib/index.js";
-import { __require as yr } from "../../../../../@stdlib_math-base-special-ln@0.2.4/node_modules/@stdlib/math-base-special-ln/lib/index.js";
-import { __require as Mr } from "../../../../../@stdlib_constants-float64-smallest-normal@0.2.2/node_modules/@stdlib/constants-float64-smallest-normal/lib/index.js";
-import { __require as cr } from "../../../../../@stdlib_constants-float64-half-pi@0.2.2/node_modules/@stdlib/constants-float64-half-pi/lib/index.js";
-import { __require as gr } from "../../../../../@stdlib_constants-float64-eps@0.2.2/node_modules/@stdlib/constants-float64-eps/lib/index.js";
-import { __require as Ar } from "./find_ibeta_inv_from_t_dist.js";
-import { __require as Fr } from "./temme1.js";
-import { __require as Nr } from "./temme2.js";
-import { __require as Or } from "./temme3.js";
-import { __require as Rr } from "./halley_iterate.js";
-import { __require as hr } from "./ibeta_roots.js";
-var w, U;
-function re() {
-  if (U) return w;
-  U = 1;
-  var V = or(), D = _r(), c = sr(), E = ar(), W = nr(), q = lr(), H = qr(), Y = Lr(), P = $r(), s = Ir(), B = xr(), Z = Tr(), p = dr(), G = yr(), g = Mr(), X = cr(), j = gr(), b = Ar(), rr = Fr(), er = Nr(), ir = Or(), fr = Rr(), mr = hr(), tr = 32, ur = 1e3, n = [0, 0, 0, 0, 0];
-  function vr(r, e, f, u) {
-    var x, v, A, _, k, L, z, C, $, T, F, m, N, l, d, O, R, J, y, o, h, a, S, M, I, K, Q, i, t;
-    if (v = !1, u === 0)
-      return [1, 0];
-    if (f === 0)
-      return [0, 1];
-    if (r === 1) {
-      if (e === 1)
-        return [f, 1 - f];
-      m = e, e = r, r = m, m = u, u = f, f = m, v = !0;
-    }
-    if (i = 0, _ = 0, L = 1, r === 0.5) {
-      if (e === 0.5)
-        return i = B(f * X), i *= i, t = B(u * X), t *= t, [i, t];
-      e > 0.5 && (m = e, e = r, r = m, m = u, u = f, f = m, v = !v);
-    }
-    if (e === 0.5 && r >= 0.5 && f !== 1)
-      h = {}, i = b(r, f, h), t = h.value;
-    else {
-      if (e === 1)
-        return f < u ? r > 1 ? (i = s(f, 1 / r), t = -c(G(f) / r)) : (i = s(f, 1 / r), t = 1 - i) : (i = P(E(-u) / r), t = -c(E(-u) / r)), v && (m = t, t = i, i = m), [i, t];
-      if (r + e > 5)
-        f > 0.5 && (m = e, e = r, r = m, m = u, u = f, f = m, v = !v), $ = p(r, e), C = Z(r, e), H($) > C - $ && $ > 5 ? (i = rr(r, e, f), t = 1 - i) : (I = r + e, k = W(H(r / I)), A = $ / I, A >= 0.2 && A <= 0.8 && I >= 10 ? (F = s(f, 1 / r), F < 25e-4 && r + e < 200 ? i = F * s(r * q(r, e), 1 / r) : i = er(f, I, k), t = 1 - i) : (r < e && (m = e, e = r, r = m, m = u, u = f, f = m, v = !v), T = 0, e < 2 && (T = q(r, e)), T === 0 ? t = 1 : (t = s(e * u * T, 1 / e), i = 1 - t)), t > 1e-5 && (i = ir(r, e, f, u), t = 1 - i));
-      else if (r < 1 && e < 1) {
-        if (o = (1 - r) / (2 - r - e), O = D(o, r, e) - f, Y(O) / f < j * 3)
-          return v ? [1 - o, o] : [o, 1 - o];
-        O < 0 && (m = e, e = r, r = m, m = u, u = f, f = m, v = !v, o = 1 - o), y = s(r * f * q(r, e), 1 / r), i = y / (1 + y), t = 1 / (1 + y), i > o && (i = o), L = o;
-      } else r > 1 && e > 1 ? (o = (r - 1) / (r + e - 2), N = (e - 1) / (r + e - 2), J = D(o, r, e) - f, J < 0 && (m = e, e = r, r = m, m = u, u = f, f = m, m = N, N = o, o = m, v = !v), R = G(f * r * q(r, e)) / r, i = P(R), t = i < 0.9 ? 1 - i : -c(R), e < r && i < 0.2 && (l = r - 1, d = e - 1, a = r * r, S = r * a, M = e * e, n[0] = 0, n[1] = 1, n[2] = d / l, l *= l, n[3] = d * (3 * r * e + 5 * e + a - r - 4) / (2 * (r + 2) * l), l *= r + 1, n[4] = d * (33 * r * M + 31 * M + 8 * a * M - 30 * r * e - 47 * e + 11 * a * e + 6 * S * e + 18 + 4 * r - S + a * a - 10 * a), n[4] /= 3 * (r + 3) * (r + 2) * l, i = V(n, i)), i > o && (i = o), L = o) : (e < r && (m = e, e = r, r = m, m = u, u = f, f = m, v = !v), s(f, 1 / r) < 0.5 ? (i = s(f * r * q(r, e), 1 / r), i === 0 && (i = g), t = 1 - i) : (t = s(1 - s(f, e * q(r, e)), 1 / e), t === 0 && (t = g), i = 1 - t));
-    }
-    return i > 0.5 && (m = e, e = r, r = m, m = u, u = f, f = m, m = t, t = i, i = m, v = !v, K = 1 - L, Q = 1 - _, _ = K, L = Q), _ === 0 && (v ? (_ = j, i < _ && (i = _)) : _ = g, i < _ && (i = _)), x = tr, i < 1e-50 && (r < 1 || e < 1) && (x *= 3, x /= 2), z = mr(r, e, f < u ? f : u, f >= u), i = fr(z, i, _, L, x, ur), i === _ && (i = 0), v ? [1 - i, i] : [i, 1 - i];
-  }
-  return w = vr, w;
-}
-export {
-  re as __require
-};
-//# sourceMappingURL=main.js.map

package/dist/node_modules/.pnpm/@stdlib_math-base-special-kernel-betaincinv@0.2.2/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/main.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"main.js","sources":["../../../../../../../../../../node_modules/.pnpm/@stdlib+math-base-special-kernel-betaincinv@0.2.2/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/main.js"],"sourcesContent":["/* eslint-disable max-statements, max-lines */\n\n/**\n* @license Apache-2.0\n*\n* Copyright (c) 2018 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*\n*\n* ## Notice\n*\n* The original C++ code and copyright notice are from the [Boost library]{@link http://www.boost.org/doc/libs/1_62_0/boost/math/special_functions/detail/ibeta_inverse.hpp}. The implementation has been modified for JavaScript.\n*\n* ```text\n* Copyright John Maddock 2006.\n* Copyright Paul A. Bristow 2007.\n*\n* Use, modification and distribution are subject to the\n* Boost Software License, Version 1.0. (See accompanying file\n* LICENSE or copy at http://www.boost.org/LICENSE_1_0.txt)\n* ```\n*/\n\n'use strict';\n\n// MODULES //\n\nvar evalpoly = require( '@stdlib/math-base-tools-evalpoly' );\nvar betainc = require( '@stdlib/math-base-special-betainc' );\nvar expm1 = require( '@stdlib/math-base-special-expm1' );\nvar log1p = require( '@stdlib/math-base-special-log1p' );\nvar asin = require( '@stdlib/math-base-special-asin' );\nvar beta = require( '@stdlib/math-base-special-beta' );\nvar sqrt = require( '@stdlib/math-base-special-sqrt' );\nvar abs = require( '@stdlib/math-base-special-abs' );\nvar exp = require( '@stdlib/math-base-special-exp' );\nvar pow = require( '@stdlib/math-base-special-pow' );\nvar sin = require( '@stdlib/math-base-special-sin' );\nvar max = require( '@stdlib/math-base-special-max' );\nvar min = require( '@stdlib/math-base-special-min' );\nvar ln = require( '@stdlib/math-base-special-ln' );\nvar FLOAT64_MIN_NORM = require( '@stdlib/constants-float64-smallest-normal' );\nvar HALF_PI = require( '@stdlib/constants-float64-half-pi' );\nvar EPSILON = require( '@stdlib/constants-float64-eps' );\nvar findIBetaInvFromTDist = require( './find_ibeta_inv_from_t_dist.js' );\nvar temme1 = require( './temme1.js' );\nvar temme2 = require( './temme2.js' );\nvar temme3 = require( './temme3.js' );\nvar halleyIterate = require( './halley_iterate.js' );\nvar ibetaRoots = require( './ibeta_roots.js' );\n\n\n// VARIABLES //\n\nvar DIGITS = 32;\nvar MAX_ITERATIONS = 1000;\n\n// Workspace for the polynomial coefficients:\nvar terms = [ 0.0, 0.0, 0.0, 0.0, 0.0 ]; // WARNING: not thread safe\n\n\n// MAIN //\n\n/**\n* Calculates the inverse of the incomplete beta function.\n*\n* @private\n* @param {PositiveNumber} a - function parameter\n* @param {PositiveNumber} b - function parameter\n* @param {Probability} p - function parameter\n* @param {Probability} q - probability equal to `1 - p`\n* @returns {Array} two-element array holding function value `y` and `1-y`\n*/\nfunction ibetaInvImp( a, b, p, q ) {\n\tvar digits;\n\tvar invert;\n\tvar lambda;\n\tvar lower;\n\tvar theta;\n\tvar upper;\n\tvar roots;\n\tvar maxv;\n\tvar minv;\n\tvar bet;\n\tvar ppa;\n\tvar tmp;\n\tvar xs2;\n\tvar ap1;\n\tvar bm1;\n\tvar fs;\n\tvar lx;\n\tvar ps;\n\tvar xg;\n\tvar xs;\n\tvar yp;\n\tvar a2;\n\tvar a3;\n\tvar b2;\n\tvar r;\n\tvar l;\n\tvar u;\n\tvar x;\n\tvar y;\n\n\t// The flag invert is set to true if we swap a for b and p for q, in which case the result has to be subtracted from 1:\n\tinvert = false;\n\n\t// Handle trivial cases first...\n\tif ( q === 0.0 ) {\n\t\treturn [ 1.0, 0.0 ];\n\t}\n\tif ( p === 0.0 ) {\n\t\treturn [ 0.0, 1.0 ];\n\t}\n\tif ( a === 1.0 ) {\n\t\tif ( b === 1.0 ) {\n\t\t\treturn [ p, 1.0-p ];\n\t\t}\n\t\t// Change things around so we can handle as b == 1 special case below:\n\t\ttmp = b;\n\t\tb = a;\n\t\ta = tmp;\n\n\t\ttmp = q;\n\t\tq = p;\n\t\tp = tmp;\n\n\t\tinvert = true;\n\t}\n\t// Depending upon which approximation method we use, we may end up calculating either x or y initially (where y = 1-x):\n\tx = 0.0; // Set to a safe zero to avoid a\n\n\t// For some of the methods we can put tighter bounds on the result than simply [0,1]:\n\tlower = 0.0;\n\tupper = 1.0;\n\n\t// Student's T with b = 0.5 gets handled as a special case, swap around if the arguments are in the \"wrong\" order:\n\tif ( a === 0.5 ) {\n\t\tif ( b === 0.5 ) {\n\t\t\tx = sin( p*HALF_PI );\n\t\t\tx *= x;\n\t\t\ty = sin( q*HALF_PI );\n\t\t\ty *= y;\n\t\t\treturn [ x, y ];\n\t\t}\n\t\tif ( b > 0.5 ) {\n\t\t\ttmp = b;\n\t\t\tb = a;\n\t\t\ta = tmp;\n\n\t\t\ttmp = q;\n\t\t\tq = p;\n\t\t\tp = tmp;\n\n\t\t\tinvert = !invert;\n\t\t}\n\t}\n\t// Select calculation method for the initial estimate:\n\tif ( b === 0.5 && a >= 0.5 && p !== 1.0 ) {\n\t\t// We have a Student's T distribution:\n\t\typ = {};\n\t\tx = findIBetaInvFromTDist( a, p, yp );\n\t\ty = yp.value;\n\t}\n\telse if ( b === 1.0 ) {\n\t\tif ( p < q ) {\n\t\t\tif ( a > 1.0 ) {\n\t\t\t\tx = pow( p, 1.0/a );\n\t\t\t\ty = -expm1( ln(p) / a );\n\t\t\t} else {\n\t\t\t\tx = pow( p, 1.0/a );\n\t\t\t\ty = 1.0 - x;\n\t\t\t}\n\t\t} else {\n\t\t\tx = exp( log1p(-q) / a );\n\t\t\ty = -expm1( log1p(-q) / a );\n\t\t}\n\t\tif ( invert ) {\n\t\t\ttmp = y;\n\t\t\ty = x;\n\t\t\tx = tmp;\n\t\t}\n\t\treturn [ x, y ];\n\t}\n\telse if ( a+b > 5.0 ) {\n\t\t// When a+b is large then we can use one of Prof Temme's asymptotic expansions, begin by swapping things around so that p < 0.5, we do this to avoid cancellations errors when p is large.\n\t\tif ( p > 0.5 ) {\n\t\t\ttmp = b;\n\t\t\tb = a;\n\t\t\ta = tmp;\n\n\t\t\ttmp = q;\n\t\t\tq = p;\n\t\t\tp = tmp;\n\n\t\t\tinvert = !invert;\n\t\t}\n\t\tminv = min( a, b );\n\t\tmaxv = max( a, b );\n\t\tif ( ( sqrt(minv) > (maxv-minv) ) && minv > 5.0 ) {\n\t\t\t// When a and b differ by a small amount the curve is quite symmetrical and we can use an error function to approximate the inverse. This is the cheapest of the three Temme expansions, and the calculated value for x will never be much larger than p, so we don't have to worry about cancellation as long as p is small.\n\t\t\tx = temme1( a, b, p );\n\t\t\ty = 1.0 - x;\n\t\t} else {\n\t\t\tr = a + b;\n\t\t\ttheta = asin( sqrt( a/r ) );\n\t\t\tlambda = minv / r;\n\t\t\tif (\n\t\t\t\tlambda >= 0.2 &&\n\t\t\t\tlambda <= 0.8 &&\n\t\t\t\tr >= 10\n\t\t\t) {\n\t\t\t\t// The second error function case is the next cheapest to use, it breaks down when the result is likely to be very small, if `a+b` is also small, but we can use a cheaper expansion there in any case. As before `x` won't be much larger than `p`, so as long as `p` is small we should be free of cancellation error.\n\t\t\t\tppa = pow( p, 1.0/a );\n\t\t\t\tif ( ppa < 0.0025 && ( a+b ) < 200.0 ) {\n\t\t\t\t\tx = ppa * pow( a*beta( a, b ), 1.0/a );\n\t\t\t\t} else {\n\t\t\t\t\tx = temme2( p, r, theta );\n\t\t\t\t}\n\t\t\t\ty = 1.0 - x;\n\t\t\t} else {\n\t\t\t\t// If we get here then a and b are very different in magnitude and we need to use the third of Temme's methods which involves inverting the incomplete gamma.  This is much more expensive than the other methods.  We also can only use this method when a > b, which can lead to cancellation errors if we really want y (as we will when x is close to 1), so a different expansion is used in that case.\n\t\t\t\tif ( a < b ) {\n\t\t\t\t\ttmp = b;\n\t\t\t\t\tb = a;\n\t\t\t\t\ta = tmp;\n\n\t\t\t\t\ttmp = q;\n\t\t\t\t\tq = p;\n\t\t\t\t\tp = tmp;\n\t\t\t\t\tinvert = !invert;\n\t\t\t\t}\n\t\t\t\t// Try and compute the easy way first:\n\t\t\t\tbet = 0.0;\n\t\t\t\tif ( b < 2.0 ) {\n\t\t\t\t\tbet = beta( a, b );\n\t\t\t\t}\n\t\t\t\tif ( bet === 0.0 ) {\n\t\t\t\t\ty = 1.0;\n\t\t\t\t} else {\n\t\t\t\t\ty = pow( b*q*bet, 1.0/b );\n\t\t\t\t\tx = 1.0 - y;\n\t\t\t\t}\n\t\t\t}\n\t\t\tif ( y > 1.0e-5 ) {\n\t\t\t\tx = temme3( a, b, p, q );\n\t\t\t\ty = 1.0 - x;\n\t\t\t}\n\t\t}\n\t}\n\telse if ( a < 1.0 && b < 1.0 ) {\n\t\t// Both a and b less than 1, there is a point of inflection at xs:\n\t\txs = ( 1.0-a ) / ( 2.0-a-b );\n\n\t\t// Now we need to ensure that we start our iteration from the right side of the inflection point:\n\t\tfs = betainc( xs, a, b ) - p;\n\t\tif ( abs(fs)/p < EPSILON*3.0 ) {\n\t\t\t// The result is at the point of inflection, best just return it:\n\t\t\tif ( invert ) {\n\t\t\t\treturn [ 1.0-xs, xs ];\n\t\t\t}\n\t\t\treturn [ xs, 1.0-xs ];\n\t\t}\n\t\tif ( fs < 0.0 ) {\n\t\t\ttmp = b;\n\t\t\tb = a;\n\t\t\ta = tmp;\n\n\t\t\ttmp = q;\n\t\t\tq = p;\n\t\t\tp = tmp;\n\n\t\t\tinvert = !invert;\n\t\t\txs = 1.0 - xs;\n\t\t}\n\t\txg = pow( a*p*beta( a, b ), 1.0/a );\n\t\tx = xg / ( 1.0+xg );\n\t\ty = 1.0 / ( 1.0+xg );\n\n\t\t// And finally we know that our result is below the inflection point, so set an upper limit on our search:\n\t\tif ( x > xs ) {\n\t\t\tx = xs;\n\t\t}\n\t\tupper = xs;\n\t}\n\telse if ( a > 1.0 && b > 1.0 ) {\n\t\t// Small a and b, both greater than 1, there is a point of inflection at xs, and it's complement is xs2, we must always start our iteration from the right side of the point of inflection.\n\t\txs = ( a-1.0 ) / ( a+b-2.0 );\n\t\txs2 = ( b-1.0 ) / ( a+b-2.0 );\n\t\tps = betainc( xs, a, b ) - p;\n\n\t\tif ( ps < 0.0 ) {\n\t\t\ttmp = b;\n\t\t\tb = a;\n\t\t\ta = tmp;\n\n\t\t\ttmp = q;\n\t\t\tq = p;\n\t\t\tp = tmp;\n\n\t\t\ttmp = xs2;\n\t\t\txs2 = xs;\n\t\t\txs = tmp;\n\n\t\t\tinvert = !invert;\n\t\t}\n\t\t// Estimate x and y, using expm1 to get a good estimate for y when it's very small:\n\t\tlx = ln( p*a*beta( a, b ) ) / a;\n\t\tx = exp( lx );\n\t\ty = ( x < 0.9 ) ? 1.0-x : -expm1(lx);\n\n\t\tif ( b < a && x < 0.2 ) {\n\t\t\t// Under a limited range of circumstances we can improve our estimate for x...\n\t\t\tap1 = a - 1.0;\n\t\t\tbm1 = b - 1.0;\n\t\t\ta2 = a * a;\n\t\t\ta3 = a * a2;\n\t\t\tb2 = b * b;\n\t\t\tterms[ 0 ] = 0.0;\n\t\t\tterms[ 1 ] = 1.0;\n\t\t\tterms[ 2 ] = bm1 / ap1;\n\t\t\tap1 *= ap1;\n\t\t\tterms[ 3 ] = bm1 * (3.0*a*b + 5.0*b + a2 - a - 4.0) / (2.0 * (a+2.0) * ap1); // eslint-disable-line max-len, no-mixed-operators\n\t\t\tap1 *= (a + 1.0);\n\t\t\tterms[ 4 ] = bm1 * (33.0*a*b2 + 31.0*b2 + 8.0*a2*b2 - 30.0*a*b - 47.0*b + 11.0*a2*b + 6.0*a3*b + 18.0 + 4.0*a - a3 + a2*a2 - 10.0*a2); // eslint-disable-line max-len, no-mixed-operators\n\t\t\tterms[ 4 ] /= (3.0 * (a+3.0) * (a+2.0) * ap1);\n\t\t\tx = evalpoly( terms, x );\n\t\t}\n\t\t// Know that result is below the inflection point, so set an upper limit on search...\n\t\tif ( x > xs ) {\n\t\t\tx = xs;\n\t\t}\n\t\tupper = xs;\n\t} else {\n\t\t// Case: ( a <= 1 ) != ( b <= 1 ). If all else fails we get here, only one of a and b is above 1, and a+b is small.  Start by swapping things around so that we have a concave curve with b > a and no points of inflection in [0,1].  As long as we expect x to be small then we can use the simple (and cheap) power term to estimate x, but when we expect x to be large then this greatly underestimates x and leaves us trying to iterate \"round the corner\" which may take almost forever. We could use Temme's inverse gamma function case in that case, this works really rather well (albeit expensively) even though strictly speaking we're outside it's defined range. However it's expensive to compute, and an alternative approach which models the curve as a distorted quarter circle is much cheaper to compute, and still keeps the number of iterations required down to a reasonable level. With thanks to Prof. Temme for this suggestion.\n\t\tif ( b < a ) {\n\t\t\ttmp = b;\n\t\t\tb = a;\n\t\t\ta = tmp;\n\n\t\t\ttmp = q;\n\t\t\tq = p;\n\t\t\tp = tmp;\n\t\t\tinvert = !invert;\n\t\t}\n\t\tif ( pow( p, 1.0/a ) < 0.5 ) {\n\t\t\tx = pow( p*a*beta( a, b ), 1.0/a );\n\t\t\tif ( x === 0.0 ) {\n\t\t\t\tx = FLOAT64_MIN_NORM;\n\t\t\t}\n\t\t\ty = 1.0 - x;\n\t\t}\n\t\t// Case: pow(q, 1/b) < 0.1\n\t\telse {\n\t\t\t// Model a distorted quarter circle:\n\t\t\ty = pow( 1.0-pow( p, b*beta( a, b ) ), 1.0/b );\n\t\t\tif ( y === 0 ) {\n\t\t\t\ty = FLOAT64_MIN_NORM;\n\t\t\t}\n\t\t\tx = 1.0 - y;\n\t\t}\n\t}\n\t// Now we have a guess for x (and for y) we can set things up for iteration.  If x > 0.5 it pays to swap things round:\n\tif ( x > 0.5 ) {\n\t\ttmp = b;\n\t\tb = a;\n\t\ta = tmp;\n\n\t\ttmp = q;\n\t\tq = p;\n\t\tp = tmp;\n\n\t\ttmp = y;\n\t\ty = x;\n\t\tx = tmp;\n\n\t\tinvert = !invert;\n\t\tl = 1.0 - upper;\n\t\tu = 1.0 - lower;\n\t\tlower = l;\n\t\tupper = u;\n\t}\n\t// Lower bound for our search:  We're not interested in denormalized answers as these tend to take up lots of iterations, given that we can't get accurate derivatives in this area (they tend to be infinite).\n\tif ( lower === 0 ) {\n\t\tif ( invert ) {\n\t\t\t// We're not interested in answers smaller than machine epsilon:\n\t\t\tlower = EPSILON;\n\t\t\tif ( x < lower ) {\n\t\t\t\tx = lower;\n\t\t\t}\n\t\t} else {\n\t\t\tlower = FLOAT64_MIN_NORM;\n\t\t}\n\t\tif ( x < lower ) {\n\t\t\tx = lower;\n\t\t}\n\t}\n\t// Figure out how many digits to iterate towards:\n\tdigits = DIGITS;\n\tif ( x < 1.0e-50 && ( a < 1.0 || b < 1.0 ) ) {\n\t\t// If we're in a region where the first derivative is very large, then we have to take care that the root-finder doesn't terminate prematurely.  We'll bump the precision up to avoid this, but we have to take care not to set the precision too high or the last few iterations will just thrash around and convergence may be slow in this case. Try 3/4 of machine epsilon:\n\t\tdigits *= 3;\n\t\tdigits /= 2;\n\t}\n\t// Now iterate, we can use either p or q as the target here depending on which is smaller:\n\troots = ibetaRoots( a, b, ( (p < q) ? p : q ), p >= q );\n\tx = halleyIterate( roots, x, lower, upper, digits, MAX_ITERATIONS );\n\n\t// Tidy up, if we \"lower\" was too high then zero is the best answer we have:\n\tif ( x === lower ) {\n\t\tx = 0.0;\n\t}\n\tif ( invert ) {\n\t\treturn [ 1.0-x, x ];\n\t}\n\treturn [ x, 1.0-x ];\n}\n\n\n// EXPORTS //\n\nmodule.exports = ibetaInvImp;\n"],"names":["evalpoly","require$$0","betainc","require$$1","expm1","require$$2","log1p","require$$3","asin","require$$4","beta","require$$5","sqrt","require$$6","abs","require$$7","exp","require$$8","pow","require$$9","sin","require$$10","max","require$$11","min","require$$12","ln","require$$13","FLOAT64_MIN_NORM","require$$14","HALF_PI","require$$15","EPSILON","require$$16","findIBetaInvFromTDist","require$$17","temme1","require$$18","temme2","require$$19","temme3","require$$20","halleyIterate","require$$21","ibetaRoots","require$$22","DIGITS","MAX_ITERATIONS","terms","ibetaInvImp","a","b","p","q","digits","invert","lambda","lower","theta","upper","roots","maxv","minv","bet","ppa","tmp","xs2","ap1","bm1","fs","lx","ps","xg","xs","yp","a2","a3","b2","r","l","u","x","y","main"],"mappings":";;;;;;;;;;;;;;;;;;;;;;;;;;;AAsCA,MAAIA,IAAWC,GAAA,GACXC,IAAUC,GAAA,GACVC,IAAQC,GAAA,GACRC,IAAQC,GAAA,GACRC,IAAOC,GAAA,GACPC,IAAOC,GAAA,GACPC,IAAOC,GAAA,GACPC,IAAMC,GAAA,GACNC,IAAMC,GAAA,GACNC,IAAMC,GAAA,GACNC,IAAMC,GAAA,GACNC,IAAMC,GAAA,GACNC,IAAMC,GAAA,GACNC,IAAKC,GAAA,GACLC,IAAmBC,GAAA,GACnBC,IAAUC,GAAA,GACVC,IAAUC,GAAA,GACVC,IAAwBC,GAAA,GACxBC,KAASC,GAAA,GACTC,KAASC,GAAA,GACTC,KAASC,GAAA,GACTC,KAAgBC,GAAA,GAChBC,KAAaC,GAAA,GAKbC,KAAS,IACTC,KAAiB,KAGjBC,IAAQ,CAAE,GAAK,GAAK,GAAK,GAAK;AAelC,WAASC,GAAaC,GAAGC,GAAGC,GAAGC,GAAI;AAClC,QAAIC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC;AAMJ,QAHA3B,IAAS,IAGJF,MAAM;AACV,aAAO,CAAE,GAAK,CAAG;AAElB,QAAKD,MAAM;AACV,aAAO,CAAE,GAAK,CAAG;AAElB,QAAKF,MAAM,GAAM;AAChB,UAAKC,MAAM;AACV,eAAO,CAAEC,GAAG,IAAIA,CAAC;AAGlB,MAAAa,IAAMd,GACNA,IAAID,GACJA,IAAIe,GAEJA,IAAMZ,GACNA,IAAID,GACJA,IAAIa,GAEJV,IAAS;AAAA,IACX;AASC,QAPA0B,IAAI,GAGJxB,IAAQ,GACRE,IAAQ,GAGHT,MAAM,KAAM;AAChB,UAAKC,MAAM;AACV,eAAA8B,IAAI7D,EAAKgC,IAAEtB,CAAO,GAClBmD,KAAKA,GACLC,IAAI9D,EAAKiC,IAAEvB,CAAO,GAClBoD,KAAKA,GACE,CAAED,GAAGC,CAAC;AAEd,MAAK/B,IAAI,QACRc,IAAMd,GACNA,IAAID,GACJA,IAAIe,GAEJA,IAAMZ,GACNA,IAAID,GACJA,IAAIa,GAEJV,IAAS,CAACA;AAAA,IAEb;AAEC,QAAKJ,MAAM,OAAOD,KAAK,OAAOE,MAAM;AAEnC,MAAAsB,IAAK,CAAA,GACLO,IAAI/C,EAAuBgB,GAAGE,GAAGsB,CAAE,GACnCQ,IAAIR,EAAG;AAAA,SAEH;AAAA,UAAKvB,MAAM;AACf,eAAKC,IAAIC,IACHH,IAAI,KACR+B,IAAI/D,EAAKkC,GAAG,IAAIF,CAAC,GACjBgC,IAAI,CAAC9E,EAAOsB,EAAG0B,CAAC,IAAIF,CAAC,MAErB+B,IAAI/D,EAAKkC,GAAG,IAAIF,CAAC,GACjBgC,IAAI,IAAMD,MAGXA,IAAIjE,EAAKV,EAAM,CAAC+C,CAAC,IAAIH,CAAC,GACtBgC,IAAI,CAAC9E,EAAOE,EAAM,CAAC+C,CAAC,IAAIH,CAAC,IAErBK,MACJU,IAAMiB,GACNA,IAAID,GACJA,IAAIhB,IAEE,CAAEgB,GAAGC,CAAC;AAET,UAAKhC,IAAEC,IAAI;AAEf,QAAKC,IAAI,QACRa,IAAMd,GACNA,IAAID,GACJA,IAAIe,GAEJA,IAAMZ,GACNA,IAAID,GACJA,IAAIa,GAEJV,IAAS,CAACA,IAEXO,IAAOtC,EAAK0B,GAAGC,CAAC,GAChBU,IAAOvC,EAAK4B,GAAGC,CAAC,GACTvC,EAAKkD,CAAI,IAAKD,IAAKC,KAAWA,IAAO,KAE3CmB,IAAI7C,GAAQc,GAAGC,GAAGC,CAAC,GACnB8B,IAAI,IAAMD,MAEVH,IAAI5B,IAAIC,GACRO,IAAQlD,EAAMI,EAAMsC,IAAE4B,CAAC,CAAE,GACzBtB,IAASM,IAAOgB,GAEftB,KAAU,OACVA,KAAU,OACVsB,KAAK,MAGLd,IAAM9C,EAAKkC,GAAG,IAAIF,CAAC,GACdc,IAAM,SAAYd,IAAEC,IAAM,MAC9B8B,IAAIjB,IAAM9C,EAAKgC,IAAExC,EAAMwC,GAAGC,CAAC,GAAI,IAAID,CAAC,IAEpC+B,IAAI3C,GAAQc,GAAG0B,GAAGpB,CAAK,GAExBwB,IAAI,IAAMD,MAGL/B,IAAIC,MACRc,IAAMd,GACNA,IAAID,GACJA,IAAIe,GAEJA,IAAMZ,GACNA,IAAID,GACJA,IAAIa,GACJV,IAAS,CAACA,IAGXQ,IAAM,GACDZ,IAAI,MACRY,IAAMrD,EAAMwC,GAAGC,CAAC,IAEZY,MAAQ,IACZmB,IAAI,KAEJA,IAAIhE,EAAKiC,IAAEE,IAAEU,GAAK,IAAIZ,CAAC,GACvB8B,IAAI,IAAMC,KAGPA,IAAI,SACRD,IAAIzC,GAAQU,GAAGC,GAAGC,GAAGC,CAAC,GACtB6B,IAAI,IAAMD;AAAA,eAIH/B,IAAI,KAAOC,IAAI,GAAM;AAM9B,YAJAsB,KAAO,IAAIvB,MAAQ,IAAIA,IAAEC,IAGzBkB,IAAKnE,EAASuE,GAAIvB,GAAGC,CAAC,IAAKC,GACtBtC,EAAIuD,CAAE,IAAEjB,IAAIpB,IAAQ;AAExB,iBAAKuB,IACG,CAAE,IAAIkB,GAAIA,CAAE,IAEb,CAAEA,GAAI,IAAIA,CAAE;AAEpB,QAAKJ,IAAK,MACTJ,IAAMd,GACNA,IAAID,GACJA,IAAIe,GAEJA,IAAMZ,GACNA,IAAID,GACJA,IAAIa,GAEJV,IAAS,CAACA,GACVkB,IAAK,IAAMA,IAEZD,IAAKtD,EAAKgC,IAAEE,IAAE1C,EAAMwC,GAAGC,CAAC,GAAI,IAAID,CAAC,GACjC+B,IAAIT,KAAO,IAAIA,IACfU,IAAI,KAAQ,IAAIV,IAGXS,IAAIR,MACRQ,IAAIR,IAELd,IAAQc;AAAA,MACV,MACM,CAAKvB,IAAI,KAAOC,IAAI,KAExBsB,KAAOvB,IAAE,MAAUA,IAAEC,IAAE,IACvBe,KAAQf,IAAE,MAAUD,IAAEC,IAAE,IACxBoB,IAAKrE,EAASuE,GAAIvB,GAAGC,CAAC,IAAKC,GAEtBmB,IAAK,MACTN,IAAMd,GACNA,IAAID,GACJA,IAAIe,GAEJA,IAAMZ,GACNA,IAAID,GACJA,IAAIa,GAEJA,IAAMC,GACNA,IAAMO,GACNA,IAAKR,GAELV,IAAS,CAACA,IAGXe,IAAK5C,EAAI0B,IAAEF,IAAExC,EAAMwC,GAAGC,CAAC,CAAE,IAAKD,GAC9B+B,IAAIjE,EAAKsD,CAAE,GACXY,IAAMD,IAAI,MAAQ,IAAIA,IAAI,CAAC7E,EAAMkE,CAAE,GAE9BnB,IAAID,KAAK+B,IAAI,QAEjBd,IAAMjB,IAAI,GACVkB,IAAMjB,IAAI,GACVwB,IAAKzB,IAAIA,GACT0B,IAAK1B,IAAIyB,GACTE,IAAK1B,IAAIA,GACTH,EAAO,CAAC,IAAK,GACbA,EAAO,CAAC,IAAK,GACbA,EAAO,KAAMoB,IAAMD,GACnBA,KAAOA,GACPnB,EAAO,CAAC,IAAKoB,KAAO,IAAIlB,IAAEC,IAAI,IAAIA,IAAIwB,IAAKzB,IAAI,MAAQ,KAAOA,IAAE,KAAOiB,IACvEA,KAAQjB,IAAI,GACZF,EAAO,CAAC,IAAKoB,KAAO,KAAKlB,IAAE2B,IAAK,KAAKA,IAAK,IAAIF,IAAGE,IAAK,KAAK3B,IAAEC,IAAI,KAAKA,IAAI,KAAKwB,IAAGxB,IAAI,IAAIyB,IAAGzB,IAAI,KAAO,IAAID,IAAI0B,IAAKD,IAAGA,IAAK,KAAKA,IAClI3B,EAAO,MAAQ,KAAOE,IAAE,MAAQA,IAAE,KAAOiB,GACzCc,IAAIjF,EAAUgD,GAAOiC,CAAC,IAGlBA,IAAIR,MACRQ,IAAIR,IAELd,IAAQc,MAGHtB,IAAID,MACRe,IAAMd,GACNA,IAAID,GACJA,IAAIe,GAEJA,IAAMZ,GACNA,IAAID,GACJA,IAAIa,GACJV,IAAS,CAACA,IAENrC,EAAKkC,GAAG,IAAIF,CAAC,IAAK,OACtB+B,IAAI/D,EAAKkC,IAAEF,IAAExC,EAAMwC,GAAGC,CAAC,GAAI,IAAID,CAAC,GAC3B+B,MAAM,MACVA,IAAIrD,IAELsD,IAAI,IAAMD,MAKVC,IAAIhE,EAAK,IAAIA,EAAKkC,GAAGD,IAAEzC,EAAMwC,GAAGC,CAAC,IAAM,IAAIA,CAAC,GACvC+B,MAAM,MACVA,IAAItD,IAELqD,IAAI,IAAMC;AAAA;AAqDZ,WAjDKD,IAAI,QACRhB,IAAMd,GACNA,IAAID,GACJA,IAAIe,GAEJA,IAAMZ,GACNA,IAAID,GACJA,IAAIa,GAEJA,IAAMiB,GACNA,IAAID,GACJA,IAAIhB,GAEJV,IAAS,CAACA,GACVwB,IAAI,IAAMpB,GACVqB,IAAI,IAAMvB,GACVA,IAAQsB,GACRpB,IAAQqB,IAGJvB,MAAU,MACTF,KAEJE,IAAQzB,GACHiD,IAAIxB,MACRwB,IAAIxB,MAGLA,IAAQ7B,GAEJqD,IAAIxB,MACRwB,IAAIxB,KAINH,IAASR,IACJmC,IAAI,UAAa/B,IAAI,KAAOC,IAAI,OAEpCG,KAAU,GACVA,KAAU,IAGXM,IAAQhB,GAAYM,GAAGC,GAAMC,IAAIC,IAAKD,IAAIC,GAAKD,KAAKC,CAAC,GACrD4B,IAAIvC,GAAekB,GAAOqB,GAAGxB,GAAOE,GAAOL,GAAQP,EAAc,GAG5DkC,MAAMxB,MACVwB,IAAI,IAEA1B,IACG,CAAE,IAAI0B,GAAGA,CAAC,IAEX,CAAEA,GAAG,IAAIA,CAAC;AAAA,EAClB;AAKA,SAAAE,IAAiBlC;;","x_google_ignoreList":[0]}

package/dist/node_modules/.pnpm/@stdlib_math-base-special-kernel-betaincinv@0.2.2/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/newton_raphson.js

@@ -1,60 +0,0 @@
-import { __require as N } from "../../../../../@stdlib_math-base-special-signum@0.2.2/node_modules/@stdlib/math-base-special-signum/lib/index.js";
-import { __require as R } from "../../../../../@stdlib_math-base-special-abs@0.2.2/node_modules/@stdlib/math-base-special-abs/lib/index.js";
-import { __require as E } from "../../../../../@stdlib_math-base-special-ldexp@0.2.3/node_modules/@stdlib/math-base-special-ldexp/lib/index.js";
-import { __require as I } from "../../../../../@stdlib_constants-float64-max@0.2.2/node_modules/@stdlib/constants-float64-max/lib/index.js";
-/**
-* @license Apache-2.0
-*
-* Copyright (c) 2018 The Stdlib Authors.
-*
-* Licensed under the Apache License, Version 2.0 (the "License");
-* you may not use this file except in compliance with the License.
-* You may obtain a copy of the License at
-*
-*    http://www.apache.org/licenses/LICENSE-2.0
-*
-* Unless required by applicable law or agreed to in writing, software
-* distributed under the License is distributed on an "AS IS" BASIS,
-* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
-* See the License for the specific language governing permissions and
-* limitations under the License.
-*
-*
-* ## Notice
-*
-* The original C++ code and copyright notice are from the [Boost library]{@link http://www.boost.org/doc/libs/1_62_0/boost/math/tools/roots.hpp}. The implementation has been modified for JavaScript.
-*
-* ```text
-* Copyright John Maddock 2006.
-*
-* Use, modification and distribution are subject to the
-* Boost Software License, Version 1.0. (See accompanying file
-* LICENSE or copy at http://www.boost.org/LICENSE_1_0.txt)
-* ```
-*/
-var d, L;
-function j() {
-  if (L) return d;
-  L = 1;
-  var h = N(), v = R(), c = E(), l = I();
-  function k(u, a, i, f, $, A) {
-    var o, _, q, w, r, p, e, n, t, b;
-    t = 0, o = 0, r = a, w = c(1, 1 - $), e = l, _ = l, q = l, p = A;
-    do {
-      if (o = t, q = _, _ = e, n = u(r), t = n[0], b = n[1], p -= 1, t === 0)
-        break;
-      if (b === 0 ? (o === 0 && (r === i ? a = f : a = i, o = u(a), e = a - r), h(o) * h(t) < 0 ? e < 0 ? e = (r - i) / 2 : e = (r - f) / 2 : e < 0 ? e = (r - f) / 2 : e = (r - i) / 2) : e = t / b, v(e * 2) > v(q) && (e = e > 0 ? (r - i) / 2 : (r - f) / 2), a = r, r -= e, r <= i) {
-        if (e = 0.5 * (a - i), r = a - e, r === i || r === f)
-          break;
-      } else if (r >= f && (e = 0.5 * (a - f), r = a - e, r === i || r === f))
-        break;
-      e > 0 ? f = a : i = a;
-    } while (p && v(r * w) < v(e));
-    return r;
-  }
-  return d = k, d;
-}
-export {
-  j as __require
-};
-//# sourceMappingURL=newton_raphson.js.map