@milaboratories/graph-maker 1.1.131 → 1.1.133

package/dist/node_modules/@milaboratories/miplots4/dist/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/inverse_students_t.js

@@ -0,0 +1,80 @@
+import { l as T } from "../../math-base-special-erfcinv/lib/index.js";
+import { l as U } from "../../math-base-special-floor/lib/index.js";
+import { l as V } from "../../math-base-special-ldexp/lib/index.js";
+import { l as W } from "../../math-base-special-round/lib/index.js";
+import { l as Y } from "../../math-base-special-acos/lib/index.js";
+import { l as Z } from "../../math-base-special-sqrt/lib/index.js";
+import { l as _ } from "../../math-base-special-abs/lib/index.js";
+import { l as g } from "../node_modules/@stdlib/math-base-special-cos/lib/index.js";
+import { l as j } from "../../math-base-special-pow/lib/index.js";
+import { l as q } from "../node_modules/@stdlib/math-base-special-sin/lib/index.js";
+import { l as z } from "../../constants-float64-sqrt-two/lib/index.js";
+import { l as A } from "../../constants-float64-pi/lib/index.js";
+import { i as B } from "./inverse_students_t_body_series.js";
+import { i as C } from "./inverse_students_t_tail_series.js";
+import { i as F } 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 R = T, y = U, k = V, D = W, G = Y, $ = Z, J = _, E = g, K = j, L = q, H = z, I = A, M = B, l = C, n = F, S = 268435456, N = 1 / 3, Q = 2 * 53 / 3, X = 0.8549879733383485;
+function e(o, i, t) {
+  var s, v, c, f, m, d, h, w, p, u, x, a, O, b, r, P;
+  if (m = 0, i > t ? (h = t, t = i, i = h, f = !0) : f = !1, y(o) === o && o < 20)
+    switch (v = k(1, Q), y(o)) {
+      case 1:
+        i === 0.5 ? m = 0 : m = -E(I * i) / L(I * i);
+        break;
+      case 2:
+        m = (2 * i - 1) / $(2 * i * t);
+        break;
+      case 4:
+        d = 4 * i * t, c = $(d), O = 4 * E(G(c) / 3) / c, b = $(O - 4), m = i - 0.5 < 0 ? -b : b;
+        break;
+      case 6:
+        if (i < 1e-150)
+          return (f ? -1 : 1) * n(o, i);
+        r = 4 * (i - i * i), P = K(r, N), a = 6 * (1 + X * (1 / P - 1));
+        do
+          p = a * a, u = p * p, x = a * u, w = a, a = 2 * (8 * r * x - 270 * p + 2187) / (5 * (4 * r * u - 216 * a - 243));
+        while (J((a - w) / a) > v);
+        a = $(a - o), m = i - 0.5 < 0 ? -a : a;
+        break;
+      default:
+        o > S ? m = R(2 * i) * H : o < 3 ? (s = 0.2742 - o * 0.0242143, i > s ? m = M(o, i) : m = l(o, i)) : (s = k(1, D(o / -0.654)), i > s ? m = n(o, i) : m = l(o, i));
+    }
+  else o > S ? m = -R(2 * i) * H : o < 3 ? (s = 0.2742 - o * 0.0242143, i > s ? m = M(o, i) : m = l(o, i)) : (s = k(1, D(o / -0.654)), i > s ? m = n(o, i) : m = l(o, i));
+  return f ? -m : m;
+}
+var ni = e;
+export {
+  ni as i
+};
+//# sourceMappingURL=inverse_students_t.js.map

package/dist/node_modules/@milaboratories/miplots4/dist/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/inverse_students_t.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"inverse_students_t.js","sources":["../../../../../../../../../node_modules/@milaboratories/miplots4/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,IAAIA,IAAUC,GACVC,IAAQC,GACRC,IAAQC,GACRC,IAAQC,GACRC,IAAOC,GACPC,IAAOC,GACPC,IAAMC,GACNC,IAAMC,GACNC,IAAMC,GACNC,IAAMC,GACNC,IAAQC,GACRC,IAAKC,GACLC,IAA6BC,GAC7BC,IAA6BC,GAC7BC,IAAuBC,GAKvBC,IAAe,WACfC,IAAY,IAAM,GAClBC,IAAQ,IAAM,KAAS,GACvBC,IAAI;AAcR,SAASC,EAAkBC,GAAIC,GAAGC,GAAI;AACrC,MAAIC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC;AAYJ,MAVAX,IAAS,GACJN,IAAIC,KAERO,IAAMP,GACNA,IAAID,GACJA,IAAIQ,GACJH,IAAS,MAETA,IAAS,IAELvC,EAAMiC,CAAE,MAAMA,KAAMA,IAAK;AAI7B,YAFAI,IAAYnC,EAAO,GAAK4B,CAAG,GAElB9B,EAAOiC,CAAE,GAAA;AAAA,MAClB,KAAK;AAECC,QAAAA,MAAM,MACVM,IAAS,IAETA,IAAS,CAAC5B,EAAKQ,IAAKc,CAAC,IAAKlB,EAAKI,IAAKc,CAAC;AAEtC;AAAA,MACD,KAAK;AAEJM,QAAAA,KAAY,IAAIN,IAAK,KAAQ1B,EAAM,IAAM0B,IAAIC,CAAC;AAC9C;AAAA,MACD,KAAK;AAEJM,QAAAA,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,MACD,KAAK;AAEJ,YAAKf,IAAI;AACR,kBAAWK,IAAW,KAAK,KAAMb,EAAsBO,GAAIC,CAAC;AAG7DgB,QAAAA,IAAI,KAAQhB,IAAKA,IAAEA,IACnBiB,IAAIrC,EAAKoC,GAAGrB,CAAS,GACrBkB,IAAI,KAAQ,IAAQhB,KAAO,IAAIoB,IAAK;AACpC;AACCP,UAAAA,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,eAC3BrC,GAAMqC,IAAIJ,KAAMI,CAAC,IAAKV;AAGhCU,QAAAA,IAAIvC,EAAMuC,IAAId,CAAE,GAChBO,IAAWN,IAAI,MAAM,IAAQ,CAACa,IAAIA;AAClC;AAAA,MACD;AACMd,QAAAA,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,IAAA;AAAA,MAIlCD,CAAAA,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,SAASK,IAAW,CAACC,IAASA;AAC/B;AAKA,IAAAY,KAAiBpB;","x_google_ignoreList":[0]}

package/dist/node_modules/@milaboratories/miplots4/dist/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/inverse_students_t_body_series.js

@@ -0,0 +1,54 @@
+import { l as i } from "../../math-base-special-gamma-delta-ratio/lib/index.js";
+import { l as t } from "../../math-base-tools-evalpoly/lib/index.js";
+import { l as f } from "../../math-base-special-sqrt/lib/index.js";
+import { l as s } from "../../constants-float64-pi/lib/index.js";
+import { p as $ } from "./polyval_co14.js";
+import { p as n } from "./polyval_co15.js";
+import { p as v } from "./polyval_co16.js";
+import { p as e } from "./polyval_co17.js";
+import { p as l } from "./polyval_co18.js";
+import { p as u } from "./polyval_co19.js";
+import { p as x } from "./polyval_co20.js";
+import { p as c } from "./polyval_co21.js";
+import { p as d } 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 g = i, h = t, j = f, k = s, q = $, B = n, D = v, I = e, O = l, P = u, R = x, S = c, T = d, b = 0, r = [1, 0, 0, 0, 0, 0, 0, 0, 0, 0];
+function w(p, a) {
+  var o, m;
+  return m = g(p / 2, 0.5) * j(p * k) * (a - 0.5), o = 1 / p, r[1] = q(o), r[2] = B(o), r[3] = D(o), r[4] = I(o), r[5] = O(o), r[6] = P(o), r[7] = R(o), r[8] = S(o), r[9] = T(o), b + m * h(r, m * m);
+}
+var Q = w;
+export {
+  Q as i
+};
+//# sourceMappingURL=inverse_students_t_body_series.js.map

package/dist/node_modules/@milaboratories/miplots4/dist/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/inverse_students_t_body_series.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"inverse_students_t_body_series.js","sources":["../../../../../../../../../node_modules/@milaboratories/miplots4/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,IAAIA,IAAkBC,GAClBC,IAAWC,GACXC,IAAOC,GACPC,IAAKC,GACLC,IAAWC,GACXC,IAAWC,GACXC,IAAWC,GACXC,IAAWC,GACXC,IAAWC,GACXC,IAAWC,GACXC,IAAWC,GACXC,IAAWC,GACXC,IAAWC,GAKXC,IAAK,GAGLC,IAAI,CAAE,GAAK,GAAK,GAAK,GAAK,GAAK,GAAK,GAAK,GAAK,GAAK,CAAG;AAiB1D,SAASC,EAA4BC,GAAIC,GAAI;AAC5C,MAAIC,GACAC;AAGJ,SAAAA,IAAIhC,EAAiB6B,IAAG,GAAG,GAAG,IAAKzB,EAAMyB,IAAGvB,CAAAA,KAASwB,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,CAAAA;AAChC;AAKA,IAAAC,IAAiBL;","x_google_ignoreList":[0]}

package/dist/node_modules/@milaboratories/miplots4/dist/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/inverse_students_t_hill.js

@@ -0,0 +1,47 @@
+import { l as u } from "../../math-base-special-erfcinv/lib/index.js";
+import { l as b } from "../../math-base-special-expm1/lib/index.js";
+import { l as c } from "../../math-base-special-sqrt/lib/index.js";
+import { l as e } from "../../math-base-special-pow/lib/index.js";
+import { l as h } from "../../constants-float64-half-pi/lib/index.js";
+import { l as x } from "../../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 = u, H = b, $ = c, P = e, S = h, n = x;
+function T(o, f) {
+  var i, t, p, s, l, m, r;
+  return o > 1e20 ? -v(2 * f) * n : (i = 1 / (o - 0.5), t = 48 / (i * i), p = ((20700 * i / t - 98) * i - 16) * i + 96.36, s = ((94.5 / (t + p) - 3) / t + 1) * $(i * S) * o, r = P(s * 2 * f, 2 / o), r > 0.05 + i ? (m = -v(2 * f) * n, r = m * m, o < 5 && (p += 0.3 * (o - 4.5) * (m + 0.6)), p += (((0.05 * s * m - 5) * m - 7) * m - 2) * m + t, r = (((((0.4 * r + 6.3) * r + 36) * r + 94.5) / p - r - 3) / t + 1) * m, r = H(i * r * r)) : r = ((1 / (((o + 6) / (o * r) - 0.089 * s - 0.822) * (o + 2) * 3) + 0.5 / (o + 4)) * r - 1) * (o + 1) / (o + 2) + 1 / r, l = $(o * r), -l);
+}
+var y = T;
+export {
+  y as i
+};
+//# sourceMappingURL=inverse_students_t_hill.js.map

package/dist/node_modules/@milaboratories/miplots4/dist/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/inverse_students_t_hill.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"inverse_students_t_hill.js","sources":["../../../../../../../../../node_modules/@milaboratories/miplots4/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,IAAIA,IAAUC,GACVC,IAAQC,GACRC,IAAOC,GACPC,IAAMC,GACNC,IAAUC,GACVC,IAAQC;AAiBZ,SAASC,EAAsBC,GAAKC,GAAI;AACvC,MAAIC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC;AAEJ,SAAKR,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;AACT;AAKA,IAAAG,IAAiBV;","x_google_ignoreList":[0]}

package/dist/node_modules/@milaboratories/miplots4/dist/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/inverse_students_t_tail_series.js

@@ -0,0 +1,46 @@
+import { l as $ } from "../../math-base-special-gamma-delta-ratio/lib/index.js";
+import { l as c } from "../../math-base-tools-evalpoly/lib/index.js";
+import { l as u } from "../../math-base-special-sqrt/lib/index.js";
+import { l as g } from "../../math-base-special-pow/lib/index.js";
+import { l as q } from "../../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 x = $, R = c, e = u, S = g, T = q, a = [0, 0, 0, 0, 0, 0, 0];
+function b(o, n) {
+  var p, l, s, i, m, t, f, v;
+  return v = x(o / 2, 0.5) * e(o * T) * n, i = o + 2, m = o + 4, t = o + 6, a[0] = 1, a[1] = -(o + 1) / (2 * i), i *= o + 2, a[2] = -o * (o + 1) * (o + 3) / (8 * i * m), i *= o + 2, a[3] = -o * (o + 1) * (o + 5) * ((3 * o + 7) * o - 2) / (48 * i * m * t), i *= o + 2, m *= o + 4, a[4] = -o * (o + 1) * (o + 7) * (((((15 * o + 154) * o + 465) * o + 286) * o - 336) * o + 64) / (384 * i * m * t * (o + 8)), i *= o + 2, a[5] = -o * (o + 1) * (o + 3) * (o + 9) * ((((((35 * o + 452) * o + 1573) * o + 600) * o - 2020) * o + 928) * o - 128) / (1280 * i * m * t * (o + 8) * (o + 10)), i *= o + 2, m *= o + 4, t *= o + 6, a[6] = -o * (o + 1) * (o + 11) * (((((((((((945 * o + 31506) * o + 425858) * o + 2980236) * o + 11266745) * o + 20675018) * o + 7747124) * o - 22574632) * o - 8565600) * o + 18108416) * o - 7099392) * o + 884736) / (46080 * i * m * t * (o + 8) * (o + 10) * (o + 12)), f = e(o), s = S(f * v, 1 / o), l = s * s, p = R(a, l), p *= f, p /= s, -p;
+}
+var y = b;
+export {
+  y as i
+};
+//# sourceMappingURL=inverse_students_t_tail_series.js.map

package/dist/node_modules/@milaboratories/miplots4/dist/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/inverse_students_t_tail_series.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"inverse_students_t_tail_series.js","sources":["../../../../../../../../../node_modules/@milaboratories/miplots4/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,IAAIA,IAAkBC,GAClBC,IAAWC,GACXC,IAAOC,GACPC,IAAMC,GACNC,IAAKC,GAMLC,IAAI,CAAE,GAAK,GAAK,GAAK,GAAK,GAAK,GAAK,CAAA;AAiBxC,SAASC,EAA4BC,GAAIC,GAAI;AAC5C,MAAIC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC;AAGJ,SAAAA,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;AACT;AAKA,IAAAQ,IAAiBX;","x_google_ignoreList":[0]}

package/dist/node_modules/@milaboratories/miplots4/dist/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/main.js

@@ -0,0 +1,60 @@
+import { l as U } from "../../math-base-tools-evalpoly/lib/index.js";
+import { l as V } from "../../math-base-special-betainc/lib/index.js";
+import { l as W } from "../../math-base-special-expm1/lib/index.js";
+import { l as Y } from "../../math-base-special-log1p/lib/index.js";
+import { l as f } from "../../math-base-special-asin/lib/index.js";
+import { l as r } from "../../math-base-special-beta/lib/index.js";
+import { l as i } from "../../math-base-special-sqrt/lib/index.js";
+import { l as v } from "../../math-base-special-abs/lib/index.js";
+import { l as o1 } from "../../math-base-special-exp/lib/index.js";
+import { l as m1 } from "../../math-base-special-pow/lib/index.js";
+import { l as t1 } from "../node_modules/@stdlib/math-base-special-sin/lib/index.js";
+import { l as a1 } from "../../math-base-special-max/lib/index.js";
+import { l as s1 } from "../../math-base-special-min/lib/index.js";
+import { l as p1 } from "../../math-base-special-ln/lib/index.js";
+import { l as l1 } from "../../constants-float64-smallest-normal/lib/index.js";
+import { l as e1 } from "../../constants-float64-half-pi/lib/index.js";
+import { l as $1 } from "../../constants-float64-eps/lib/index.js";
+import { f as u1 } from "./find_ibeta_inv_from_t_dist.js";
+import { t as n1 } from "./temme1.js";
+import { t as c1 } from "./temme2.js";
+import { t as d1 } from "./temme3.js";
+import { h as h1 } from "./halley_iterate.js";
+import { i as x1 } from "./ibeta_roots.js";
+var M1 = U, H = V, q = W, J = Y, g1 = f, M = r, P = i, y1 = v, Q = o1, n = m1, S = t1, A1 = a1, F1 = s1, X = p1, w = l1, Z = e1, _ = $1, I1 = u1, K1 = n1, L1 = c1, N1 = d1, O1 = h1, R1 = x1, T1 = 32, b1 = 1e3, h = [0, 0, 0, 0, 0];
+function j1(o, m, a, l) {
+  var A, e, N, u, z, x, B, C, g, F, O, s, R, d, I, T, b, D, K, $, j, c, k, L, y, E, G, t, p;
+  if (e = !1, l === 0)
+    return [1, 0];
+  if (a === 0)
+    return [0, 1];
+  if (o === 1) {
+    if (m === 1)
+      return [a, 1 - a];
+    s = m, m = o, o = s, s = l, l = a, a = s, e = !0;
+  }
+  if (t = 0, u = 0, x = 1, o === 0.5) {
+    if (m === 0.5)
+      return t = S(a * Z), t *= t, p = S(l * Z), p *= p, [t, p];
+    m > 0.5 && (s = m, m = o, o = s, s = l, l = a, a = s, e = !e);
+  }
+  if (m === 0.5 && o >= 0.5 && a !== 1)
+    j = {}, t = I1(o, a, j), p = j.value;
+  else {
+    if (m === 1)
+      return a < l ? o > 1 ? (t = n(a, 1 / o), p = -q(X(a) / o)) : (t = n(a, 1 / o), p = 1 - t) : (t = Q(J(-l) / o), p = -q(J(-l) / o)), e && (s = p, p = t, t = s), [t, p];
+    if (o + m > 5)
+      a > 0.5 && (s = m, m = o, o = s, s = l, l = a, a = s, e = !e), g = F1(o, m), C = A1(o, m), P(g) > C - g && g > 5 ? (t = K1(o, m, a), p = 1 - t) : (y = o + m, z = g1(P(o / y)), N = g / y, N >= 0.2 && N <= 0.8 && y >= 10 ? (O = n(a, 1 / o), O < 25e-4 && o + m < 200 ? t = O * n(o * M(o, m), 1 / o) : t = L1(a, y, z), p = 1 - t) : (o < m && (s = m, m = o, o = s, s = l, l = a, a = s, e = !e), F = 0, m < 2 && (F = M(o, m)), F === 0 ? p = 1 : (p = n(m * l * F, 1 / m), t = 1 - p)), p > 1e-5 && (t = N1(o, m, a, l), p = 1 - t));
+    else if (o < 1 && m < 1) {
+      if ($ = (1 - o) / (2 - o - m), T = H($, o, m) - a, y1(T) / a < _ * 3)
+        return e ? [1 - $, $] : [$, 1 - $];
+      T < 0 && (s = m, m = o, o = s, s = l, l = a, a = s, e = !e, $ = 1 - $), K = n(o * a * M(o, m), 1 / o), t = K / (1 + K), p = 1 / (1 + K), t > $ && (t = $), x = $;
+    } else o > 1 && m > 1 ? ($ = (o - 1) / (o + m - 2), R = (m - 1) / (o + m - 2), D = H($, o, m) - a, D < 0 && (s = m, m = o, o = s, s = l, l = a, a = s, s = R, R = $, $ = s, e = !e), b = X(a * o * M(o, m)) / o, t = Q(b), p = t < 0.9 ? 1 - t : -q(b), m < o && t < 0.2 && (d = o - 1, I = m - 1, c = o * o, k = o * c, L = m * m, h[0] = 0, h[1] = 1, h[2] = I / d, d *= d, h[3] = I * (3 * o * m + 5 * m + c - o - 4) / (2 * (o + 2) * d), d *= o + 1, h[4] = I * (33 * o * L + 31 * L + 8 * c * L - 30 * o * m - 47 * m + 11 * c * m + 6 * k * m + 18 + 4 * o - k + c * c - 10 * c), h[4] /= 3 * (o + 3) * (o + 2) * d, t = M1(h, t)), t > $ && (t = $), x = $) : (m < o && (s = m, m = o, o = s, s = l, l = a, a = s, e = !e), n(a, 1 / o) < 0.5 ? (t = n(a * o * M(o, m), 1 / o), t === 0 && (t = w), p = 1 - t) : (p = n(1 - n(a, m * M(o, m)), 1 / m), p === 0 && (p = w), t = 1 - p));
+  }
+  return t > 0.5 && (s = m, m = o, o = s, s = l, l = a, a = s, s = p, p = t, t = s, e = !e, E = 1 - x, G = 1 - u, u = E, x = G), u === 0 && (e ? (u = _, t < u && (t = u)) : u = w, t < u && (t = u)), A = T1, t < 1e-50 && (o < 1 || m < 1) && (A *= 3, A /= 2), B = R1(o, m, a < l ? a : l, a >= l), t = O1(B, t, u, x, A, b1), t === u && (t = 0), e ? [1 - t, t] : [t, 1 - t];
+}
+var i1 = j1;
+export {
+  i1 as m
+};
+//# sourceMappingURL=main.js.map

package/dist/node_modules/@milaboratories/miplots4/dist/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/main.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"main.js","sources":["../../../../../../../../../node_modules/@milaboratories/miplots4/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,IAAIA,KAAWC,GACXC,IAAUC,GACVC,IAAQC,GACRC,IAAQC,GACRC,KAAOC,GACPC,IAAOC,GACPC,IAAOC,GACPC,KAAMC,GACNC,IAAMC,IACNC,IAAMC,IACNC,IAAMC,IACNC,KAAMC,IACNC,KAAMC,IACNC,IAAKC,IACLC,IAAmBC,IACnBC,IAAUC,IACVC,IAAUC,IACVC,KAAwBC,IACxBC,KAASC,IACTC,KAASC,IACTC,KAASC,IACTC,KAAgBC,IAChBC,KAAaC,IAKbC,KAAS,IACTC,KAAiB,KAGjBC,IAAQ,CAAE,GAAK,GAAK,GAAK,GAAK,CAAA;AAelC,SAASC,GAAaC,GAAGC,GAAGC,GAAGC,GAAI;AAClC,MAAIC,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,MAHA3B,IAAS,IAGJF,MAAM;AACV,WAAO,CAAE,GAAK,CAAG;AAElB,MAAKD,MAAM;AACV,WAAO,CAAE,GAAK,CAAG;AAElB,MAAKF,MAAM,GAAM;AAChB,QAAKC,MAAM;AACV,aAAO,CAAEC,GAAG,IAAIA,CAAC;AAGlBa,IAAAA,IAAMd,GACNA,IAAID,GACJA,IAAIe,GAEJA,IAAMZ,GACNA,IAAID,GACJA,IAAIa,GAEJV,IAAS;AAAA,EACX;AASC,MAPA0B,IAAI,GAGJxB,IAAQ,GACRE,IAAQ,GAGHT,MAAM,KAAM;AAChB,QAAKC,MAAM;AACV,aAAA8B,IAAI7D,EAAKgC,IAAEtB,CAAO,GAClBmD,KAAKA,GACLC,IAAI9D,EAAKiC,IAAEvB,CAAO,GAClBoD,KAAKA,GACE,CAAED,GAAGC,CAAC;AAET/B,IAAAA,IAAI,QACRc,IAAMd,GACNA,IAAID,GACJA,IAAIe,GAEJA,IAAMZ,GACNA,IAAID,GACJA,IAAIa,GAEJV,IAAS,CAACA;AAAAA,EAEb;AAEC,MAAKJ,MAAM,OAAOD,KAAK,OAAOE,MAAM;AAEnCsB,IAAAA,IAAK,CAAA,GACLO,IAAI/C,GAAuBgB,GAAGE,GAAGsB,CAAE,GACnCQ,IAAIR,EAAG;AAAA,OAEH;AAAA,QAAKvB,MAAM;AACf,aAAKC,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,QAAKhC,IAAEC,IAAI;AAEVC,MAAAA,IAAI,QACRa,IAAMd,GACNA,IAAID,GACJA,IAAIe,GAEJA,IAAMZ,GACNA,IAAID,GACJA,IAAIa,GAEJV,IAAS,CAACA,IAEXO,IAAOtC,GAAK0B,GAAGC,CAAC,GAChBU,IAAOvC,GAAK4B,GAAGC,CAAC,GACTvC,EAAKkD,CAAI,IAAKD,IAAKC,KAAWA,IAAO,KAE3CmB,IAAI7C,GAAQc,GAAGC,GAAGC,CAAC,GACnB8B,IAAI,IAAMD,MAEVH,IAAI5B,IAAIC,GACRO,IAAQlD,GAAMI,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;AAAAA,aAIH/B,IAAI,KAAOC,IAAI,GAAM;AAM9B,UAJAsB,KAAO,IAAIvB,MAAQ,IAAIA,IAAEC,IAGzBkB,IAAKnE,EAASuE,GAAIvB,GAAGC,CAAC,IAAKC,GACtBtC,GAAIuD,CAAE,IAAEjB,IAAIpB,IAAQ;AAExB,eAAKuB,IACG,CAAE,IAAIkB,GAAIA,CAAE,IAEb,CAAEA,GAAI,IAAIA,CAAE;AAEfJ,MAAAA,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;AAAAA,IACV,MACWvB,CAAAA,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,CAAA,IAAMoB,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,CAAA,KAAQ,KAAOE,IAAE,MAAQA,IAAE,KAAOiB,GACzCc,IAAIjF,GAAUgD,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,CAAA,GAAM,IAAIA,CAAC,GACvC+B,MAAM,MACVA,IAAItD,IAELqD,IAAI,IAAMC;AAAAA,EAAA;AAqDZ,SAjDKD,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;AAClB;AAKA,IAAAE,KAAiBlC;","x_google_ignoreList":[0]}

package/dist/node_modules/@milaboratories/miplots4/dist/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/newton_raphson.js

@@ -0,0 +1,55 @@
+import { l as x } from "../../math-base-special-signum/lib/index.js";
+import { l as M } from "../../math-base-special-abs/lib/index.js";
+import { l as R } from "../../math-base-special-ldexp/lib/index.js";
+import { l as e } from "../../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 u = x, s = M, f = R, d = e;
+function g(h, a, o, m, w, q) {
+  var p, b, n, k, r, v, i, $, t, c;
+  t = 0, p = 0, r = a, k = f(1, 1 - w), i = d, b = d, n = d, v = q;
+  do {
+    if (p = t, n = b, b = i, $ = h(r), t = $[0], c = $[1], v -= 1, t === 0)
+      break;
+    if (c === 0 ? (p === 0 && (r === o ? a = m : a = o, p = h(a), i = a - r), u(p) * u(t) < 0 ? i < 0 ? i = (r - o) / 2 : i = (r - m) / 2 : i < 0 ? i = (r - m) / 2 : i = (r - o) / 2) : i = t / c, s(i * 2) > s(n) && (i = i > 0 ? (r - o) / 2 : (r - m) / 2), a = r, r -= i, r <= o) {
+      if (i = 0.5 * (a - o), r = a - i, r === o || r === m)
+        break;
+    } else if (r >= m && (i = 0.5 * (a - m), r = a - i, r === o || r === m))
+      break;
+    i > 0 ? m = a : o = a;
+  } while (v && s(r * k) < s(i));
+  return r;
+}
+var B = g;
+export {
+  B as n
+};
+//# sourceMappingURL=newton_raphson.js.map

package/dist/node_modules/@milaboratories/miplots4/dist/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/newton_raphson.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"newton_raphson.js","sources":["../../../../../../../../../node_modules/@milaboratories/miplots4/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/newton_raphson.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/tools/roots.hpp}. The implementation has been modified for JavaScript.\n*\n* ```text\n* Copyright John Maddock 2006.\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 sign = require( '@stdlib/math-base-special-signum' );\nvar abs = require( '@stdlib/math-base-special-abs' );\nvar ldexp = require( '@stdlib/math-base-special-ldexp' );\nvar MAX_VALUE = require( '@stdlib/constants-float64-max' );\n\n\n// MAIN //\n\n/**\n* Performs root finding via second order Newton-Raphson iteration.\n*\n* @private\n* @param {Array} fun - two-element array of the function and its first derivative\n* @param {number} guess - initial starting value.\n* @param {number} min - minimum possible value for the result,used as initial lower bracket.\n* @param {number} max - maximum possible value for the result, used as initial upper bracket.\n* @param {PositiveInteger} digits - desired number of binary digits\n* @param {PositiveInteger} maxIter - maximum number of iterations\n* @returns {number} function value\n*/\nfunction newtonRaphsonIterate( fun, guess, min, max, digits, maxIter ) {\n\tvar f0last;\n\tvar delta1;\n\tvar delta2;\n\tvar factor;\n\tvar result;\n\tvar count;\n\tvar delta;\n\tvar res;\n\tvar f0;\n\tvar f1;\n\n\tf0 = 0.0;\n\tf0last = 0.0;\n\tresult = guess;\n\n\tfactor = ldexp( 1.0, 1.0 - digits );\n\tdelta = MAX_VALUE;\n\tdelta1 = MAX_VALUE;\n\tdelta2 = MAX_VALUE;\n\n\tcount = maxIter;\n\tdo {\n\t\tf0last = f0;\n\t\tdelta2 = delta1;\n\t\tdelta1 = delta;\n\t\tres = fun(result);\n\t\tf0 = res[ 0 ];\n\t\tf1 = res[ 1 ];\n\t\tcount -= 1;\n\t\tif ( f0 === 0.0 ) {\n\t\t\tbreak;\n\t\t}\n\t\tif ( f1 === 0.0 ) {\n\t\t\t// Oops zero derivative!!!\n\t\t\tif ( f0last === 0.0 ) {\n\t\t\t\t// Must be the first iteration, pretend that we had a previous one at either min or max:\n\t\t\t\tif ( result === min ) {\n\t\t\t\t\tguess = max;\n\t\t\t\t} else {\n\t\t\t\t\tguess = min;\n\t\t\t\t}\n\t\t\t\tf0last = fun( guess );\n\t\t\t\tdelta = guess - result;\n\t\t\t}\n\t\t\tif ( sign(f0last) * sign(f0) < 0 ) {\n\t\t\t\t// We've crossed over so move in opposite direction to last step:\n\t\t\t\tif ( delta < 0 ) {\n\t\t\t\t\tdelta = (result - min) / 2.0;\n\t\t\t\t} else {\n\t\t\t\t\tdelta = (result - max) / 2.0;\n\t\t\t\t}\n\t\t\t} else if ( delta < 0 ) {\n\t\t\t\tdelta = (result - max) / 2.0;\n\t\t\t} else {\n\t\t\t\tdelta = (result - min) / 2.0;\n\t\t\t}\n\t\t} else {\n\t\t\tdelta = f0 / f1;\n\t\t}\n\t\tif ( abs(delta * 2.0) > abs(delta2) ) {\n\t\t\t// Last two steps haven't converged, try bisection:\n\t\t\tdelta = ( delta > 0.0 ) ? (result-min) / 2.0 : (result-max) / 2.0;\n\t\t}\n\t\tguess = result;\n\t\tresult -= delta;\n\t\tif ( result <= min ) {\n\t\t\tdelta = 0.5 * (guess - min);\n\t\t\tresult = guess - delta;\n\t\t\tif ( result === min || result === max ) {\n\t\t\t\tbreak;\n\t\t\t}\n\t\t} else if ( result >= max ) {\n\t\t\tdelta = 0.5 * (guess - max);\n\t\t\tresult = guess - delta;\n\t\t\tif ( result === min || result === max ) {\n\t\t\t\tbreak;\n\t\t\t}\n\t\t}\n\t\t// Update brackets:\n\t\tif ( delta > 0.0 ) {\n\t\t\tmax = guess;\n\t\t} else {\n\t\t\tmin = guess;\n\t\t}\n\t}\n\twhile ( count && ( abs(result * factor) < abs(delta) ) );\n\n\treturn result;\n}\n\n\n// EXPORTS //\n\nmodule.exports = newtonRaphsonIterate;\n"],"names":["sign","require$$0","abs","require$$1","ldexp","require$$2","MAX_VALUE","require$$3","newtonRaphsonIterate","fun","guess","min","max","digits","maxIter","f0last","delta1","delta2","factor","result","count","delta","res","f0","f1","newton_raphson"],"mappings":";;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;AAmCA,IAAIA,IAAOC,GACPC,IAAMC,GACNC,IAAQC,GACRC,IAAYC;AAiBhB,SAASC,EAAsBC,GAAKC,GAAOC,GAAKC,GAAKC,GAAQC,GAAU;AACtE,MAAIC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC,GACAC;AAEJD,MAAK,GACLR,IAAS,GACTI,IAAST,GAETQ,IAASd,EAAO,GAAK,IAAMS,CAAM,GACjCQ,IAAQf,GACRU,IAASV,GACTW,IAASX,GAETc,IAAQN;AACR,KAAG;AAQF,QAPAC,IAASQ,GACTN,IAASD,GACTA,IAASK,GACTC,IAAMb,EAAIU,CAAM,GAChBI,IAAKD,EAAK,CAAC,GACXE,IAAKF,EAAK,CAAC,GACXF,KAAS,GACJG,MAAO;AACX;AAmCD,QAjCKC,MAAO,KAENT,MAAW,MAEVI,MAAWR,IACfD,IAAQE,IAERF,IAAQC,GAETI,IAASN,EAAKC,CAAK,GACnBW,IAAQX,IAAQS,IAEZnB,EAAKe,CAAM,IAAIf,EAAKuB,CAAE,IAAI,IAEzBF,IAAQ,IACZA,KAASF,IAASR,KAAO,IAEzBU,KAASF,IAASP,KAAO,IAEfS,IAAQ,IACnBA,KAASF,IAASP,KAAO,IAEzBS,KAASF,IAASR,KAAO,KAG1BU,IAAQE,IAAKC,GAETtB,EAAImB,IAAQ,CAAG,IAAInB,EAAIe,CAAM,MAEjCI,IAAUA,IAAQ,KAASF,IAAOR,KAAO,KAAOQ,IAAOP,KAAO,IAE/DF,IAAQS,GACRA,KAAUE,GACLF,KAAUR;AAGd,UAFAU,IAAQ,OAAOX,IAAQC,IACvBQ,IAAST,IAAQW,GACZF,MAAWR,KAAOQ,MAAWP;AACjC;AAAA,eAEUO,KAAUP,MACrBS,IAAQ,OAAOX,IAAQE,IACvBO,IAAST,IAAQW,GACZF,MAAWR,KAAOQ,MAAWP;AACjC;AAIGS,IAAAA,IAAQ,IACZT,IAAMF,IAENC,IAAMD;AAAAA,EAET,SACSU,KAAWlB,EAAIiB,IAASD,CAAM,IAAIhB,EAAImB,CAAK;AAEnD,SAAOF;AACR;AAKA,IAAAM,IAAiBjB;","x_google_ignoreList":[0]}

package/dist/node_modules/@milaboratories/miplots4/dist/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/polyval_co1.js

@@ -0,0 +1,25 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2022 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.
+*/
+function n(a) {
+  return a === 0 ? -1 : -1 + a * (-5 + a * 5);
+}
+var o = n;
+export {
+  o as p
+};
+//# sourceMappingURL=polyval_co1.js.map

package/dist/node_modules/@milaboratories/miplots4/dist/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/polyval_co1.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"polyval_co1.js","sources":["../../../../../../../../../node_modules/@milaboratories/miplots4/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/polyval_co1.js"],"sourcesContent":["/**\n* @license Apache-2.0\n*\n* Copyright (c) 2022 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/* This is a generated file. Do not edit directly. */\n'use strict';\n\n// MAIN //\n\n/**\n* Evaluates a polynomial.\n*\n* ## Notes\n*\n* -   The implementation uses [Horner's rule][horners-method] for efficient computation.\n*\n* [horners-method]: https://en.wikipedia.org/wiki/Horner%27s_method\n*\n* @private\n* @param {number} x - value at which to evaluate the polynomial\n* @returns {number} evaluated polynomial\n*/\nfunction evalpoly( x ) {\n\tif ( x === 0.0 ) {\n\t\treturn -1.0;\n\t}\n\treturn -1.0 + (x * (-5.0 + (x * 5.0)));\n}\n\n\n// EXPORTS //\n\nmodule.exports = evalpoly;\n"],"names":["evalpoly","x","polyval_co1"],"mappings":";;;;;;;;;;;;;;;;;AAoCA,SAASA,EAAUC,GAAI;AACtB,SAAKA,MAAM,IACH,KAED,KAAQA,KAAK,KAAQA,IAAI;AACjC;AAKA,IAAAC,IAAiBF;","x_google_ignoreList":[0]}

package/dist/node_modules/@milaboratories/miplots4/dist/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/polyval_co10.js

@@ -0,0 +1,25 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2022 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.
+*/
+function n(a) {
+  return a === 0 ? 63149 : 63149 + a * (-151557 + a * (140052 + a * (-727469 + a * (2239932 + a * (-2251437 + a * 750479)))));
+}
+var o = n;
+export {
+  o as p
+};
+//# sourceMappingURL=polyval_co10.js.map

package/dist/node_modules/@milaboratories/miplots4/dist/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/polyval_co10.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"polyval_co10.js","sources":["../../../../../../../../../node_modules/@milaboratories/miplots4/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/polyval_co10.js"],"sourcesContent":["/**\n* @license Apache-2.0\n*\n* Copyright (c) 2022 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/* This is a generated file. Do not edit directly. */\n'use strict';\n\n// MAIN //\n\n/**\n* Evaluates a polynomial.\n*\n* ## Notes\n*\n* -   The implementation uses [Horner's rule][horners-method] for efficient computation.\n*\n* [horners-method]: https://en.wikipedia.org/wiki/Horner%27s_method\n*\n* @private\n* @param {number} x - value at which to evaluate the polynomial\n* @returns {number} evaluated polynomial\n*/\nfunction evalpoly( x ) {\n\tif ( x === 0.0 ) {\n\t\treturn 63149.0;\n\t}\n\treturn 63149.0 + (x * (-151557.0 + (x * (140052.0 + (x * (-727469.0 + (x * (2239932.0 + (x * (-2251437.0 + (x * 750479.0))))))))))); // eslint-disable-line max-len\n}\n\n\n// EXPORTS //\n\nmodule.exports = evalpoly;\n"],"names":["evalpoly","x","polyval_co10"],"mappings":";;;;;;;;;;;;;;;;;AAoCA,SAASA,EAAUC,GAAI;AACtB,SAAKA,MAAM,IACH,QAED,QAAWA,KAAK,UAAaA,KAAK,SAAYA,KAAK,UAAaA,KAAK,UAAaA,KAAK,WAAcA,IAAI;AACjH;AAKA,IAAAC,IAAiBF;","x_google_ignoreList":[0]}

package/dist/node_modules/@milaboratories/miplots4/dist/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/polyval_co11.js

@@ -0,0 +1,25 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2022 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.
+*/
+function n(a) {
+  return a === 0 ? 29233 : 29233 + a * (-78755 + a * (105222 + a * (146879 + a * (-1602610 + a * (3195183 + a * (-2554139 + a * 729754))))));
+}
+var o = n;
+export {
+  o as p
+};
+//# sourceMappingURL=polyval_co11.js.map

package/dist/node_modules/@milaboratories/miplots4/dist/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/polyval_co11.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"polyval_co11.js","sources":["../../../../../../../../../node_modules/@milaboratories/miplots4/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/polyval_co11.js"],"sourcesContent":["/**\n* @license Apache-2.0\n*\n* Copyright (c) 2022 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/* This is a generated file. Do not edit directly. */\n'use strict';\n\n// MAIN //\n\n/**\n* Evaluates a polynomial.\n*\n* ## Notes\n*\n* -   The implementation uses [Horner's rule][horners-method] for efficient computation.\n*\n* [horners-method]: https://en.wikipedia.org/wiki/Horner%27s_method\n*\n* @private\n* @param {number} x - value at which to evaluate the polynomial\n* @returns {number} evaluated polynomial\n*/\nfunction evalpoly( x ) {\n\tif ( x === 0.0 ) {\n\t\treturn 29233.0;\n\t}\n\treturn 29233.0 + (x * (-78755.0 + (x * (105222.0 + (x * (146879.0 + (x * (-1602610.0 + (x * (3195183.0 + (x * (-2554139.0 + (x * 729754.0))))))))))))); // eslint-disable-line max-len\n}\n\n\n// EXPORTS //\n\nmodule.exports = evalpoly;\n"],"names":["evalpoly","x","polyval_co11"],"mappings":";;;;;;;;;;;;;;;;;AAoCA,SAASA,EAAUC,GAAI;AACtB,SAAKA,MAAM,IACH,QAED,QAAWA,KAAK,SAAYA,KAAK,SAAYA,KAAK,SAAYA,KAAK,WAAcA,KAAK,UAAaA,KAAK,WAAcA,IAAI;AAClI;AAKA,IAAAC,IAAiBF;","x_google_ignoreList":[0]}

package/dist/node_modules/@milaboratories/miplots4/dist/node_modules/@stdlib/math-base-special-kernel-betaincinv/lib/polyval_co12.js

@@ -0,0 +1,25 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2022 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.
+*/
+function n(a) {
+  return a === 0 ? 1 : 1 + a * (-13 + a * 13);
+}
+var o = n;
+export {
+  o as p
+};
+//# sourceMappingURL=polyval_co12.js.map