nerdamer 1.1.13 → 2.0.0-rc.1

package/output/core/functions/decimal.js

@@ -0,0 +1,664 @@
+"use strict";
+var __importDefault = (this && this.__importDefault) || function (mod) {
+    return (mod && mod.__esModule) ? mod : { "default": mod };
+};
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.EULER_GAMMA = void 0;
+exports.erf = erf;
+exports.gamma = gamma;
+exports.sinc = sinc;
+exports.Si = Si;
+exports.Shi = Shi;
+exports.Ci = Ci;
+exports.Chi = Chi;
+exports.Ei = Ei;
+exports.Li = Li;
+const decimal_js_1 = __importDefault(require("decimal.js"));
+const bigint_1 = require("./bigint/bigint");
+/*
+ * TODO: Currently only up to a 100 digits. This will need to be calculated in the future.
+ */
+exports.EULER_GAMMA = new decimal_js_1.default('0.5772156649015328606065120900824024310421593359399235988057672348848677267776646709369470632917467495');
+const DEFAULT_SERIES_GUARD_DIGITS = 8;
+const MIN_PRECISION_SERIES_ITERATIONS = 1000;
+const SERIES_ITERATIONS_PER_PRECISION_DIGIT = 8;
+const DECIMAL_SINC_SERIES_CUTOFF = '1e-8';
+const OSCILLATORY_ASYMPTOTIC_CUTOFF = 1e9;
+const RATIONAL_APPROXIMATION_CUTOFF = 4;
+const ASYMPTOTIC_APPROXIMATION_CUTOFF = 8;
+const EXPONENTIAL_ASYMPTOTIC_CUTOFF = 40;
+const SERIES_MAX_ITERATIONS = 200;
+const ASYMPTOTIC_MAX_TERMS = 100;
+const POWER_SERIES_MAX_TERMS = 500;
+function getSeriesError(guardDigits = DEFAULT_SERIES_GUARD_DIGITS) {
+    return new decimal_js_1.default(10).pow(-(decimal_js_1.default.precision + guardDigits));
+}
+function gammaSpouge(x) {
+    const precision = decimal_js_1.default.precision;
+    decimal_js_1.default.set({ precision: Math.floor(precision * 1.45) });
+    try {
+        const n = precision;
+        const coeffs = calcCoeffs(n);
+        let accm = coeffs[0];
+        for (let k = 1; k < n; k++) {
+            accm = accm.plus(coeffs[k].div(x.plus(k)));
+        }
+        accm = accm.times(decimal_js_1.default.exp(x.plus(n).neg())).times(decimal_js_1.default.pow(x.plus(n), x.plus(0.5)));
+        return accm.div(x);
+    }
+    finally {
+        decimal_js_1.default.set({ precision });
+    }
+}
+function erf(x) {
+    x = new decimal_js_1.default(x);
+    const pi = decimal_js_1.default.acos(-1);
+    const twoOverSqrtPi = new decimal_js_1.default(2).div(pi.sqrt());
+    const error = getSeriesError();
+    const sgn = decimal_js_1.default.sign(x);
+    x = x.abs();
+    if (x.isZero()) {
+        return new decimal_js_1.default(0);
+    }
+    if (x.lt(2)) {
+        let sum = new decimal_js_1.default(0);
+        let k = 0;
+        let dx = x;
+        const max = Math.max(MIN_PRECISION_SERIES_ITERATIONS, decimal_js_1.default.precision * SERIES_ITERATIONS_PER_PRECISION_DIGIT);
+        while (dx.gte(error)) {
+            const num = new decimal_js_1.default(-1).pow(k).times(x.pow(2 * k + 1));
+            const den = new decimal_js_1.default(2 * k + 1).times(String((0, bigint_1.factorial)(k)));
+            const result = sum.plus(num.div(den));
+            dx = result.minus(sum).abs();
+            sum = result;
+            k++;
+            if (k > max) {
+                throw new Error('erf not converging. Exiting!');
+            }
+        }
+        let retval = decimal_js_1.default.min(twoOverSqrtPi.times(sum), 1);
+        if (sgn === -1) {
+            retval = retval.neg();
+        }
+        return retval;
+    }
+    // Asymptotic expansion of erfc(x) for larger |x|:
+    // erfc(x) ~ exp(-x^2)/(sqrt(pi)*x) * (1 - 1/(2x^2) + 3/(4x^4) - 15/(8x^6) + ...)
+    const xx = x.mul(x);
+    const twoX2 = xx.mul(2);
+    let term = new decimal_js_1.default(1);
+    let series = term;
+    let best = term.abs();
+    const max = Math.max(MIN_PRECISION_SERIES_ITERATIONS, decimal_js_1.default.precision * SERIES_ITERATIONS_PER_PRECISION_DIGIT);
+    for (let n = 1; n < max; n++) {
+        term = term.mul(-(2 * n - 1)).div(twoX2);
+        const at = term.abs();
+        if (at.gt(best)) {
+            break;
+        }
+        best = at;
+        series = series.plus(term);
+        if (at.lt(error)) {
+            break;
+        }
+    }
+    const erfc = x.neg().mul(x).exp().div(pi.sqrt().mul(x)).times(series);
+    let retval = new decimal_js_1.default(1).minus(erfc);
+    if (sgn === -1) {
+        retval = retval.neg();
+    }
+    return retval;
+}
+function calcCoeffs(n) {
+    const c = [];
+    let kFact = new decimal_js_1.default(1);
+    // Get pi
+    const pi = decimal_js_1.default.acos(-1);
+    // sqrt(2*pi);
+    c[0] = pi.times(2).sqrt();
+    for (let k = 1; k < n; k++) {
+        c[k] = decimal_js_1.default.exp(n - k)
+            .times(decimal_js_1.default.pow(n - k, k - 0.5))
+            .div(kFact);
+        kFact = kFact.times(-k);
+    }
+    return c;
+}
+// https://math.stackexchange.com/questions/2204020/definition-of-the-gamma-function-for-non-integer-negative-values
+/**
+ * The gamma function using Spouge approximation.
+ * This function is slow converging but appropriate for arbitrary precision. It relies on Decimal.js.
+ *
+ * Ported from: https://rosettacode.org/wiki/Gamma_function#C++ (the implementation)
+ * https://math.stackexchange.com/questions/2204020/definition-of-the-gamma-function-for-non-integer-negative-values (negative gamma)
+ * https://en.wikipedia.org/wiki/Spouge%27s_approximation (Spouge's Approximation)
+ *
+ * @param n The number of elements to be calculated in the coefficient array
+ * @param x The number being evaluated
+ * @returns
+ */
+function gamma(x) {
+    x = new decimal_js_1.default(x);
+    if (x.isInteger()) {
+        if (x.lte(0)) {
+            return new decimal_js_1.default(Infinity);
+        }
+        return new decimal_js_1.default(String((0, bigint_1.factorial)(Number(x.minus(1)))));
+    }
+    if (x.lt(0)) {
+        const pi = decimal_js_1.default.acos(-1);
+        return pi.div(pi
+            .times(x)
+            .sin()
+            .times(gammaSpouge(new decimal_js_1.default(1).minus(x))));
+    }
+    return gammaSpouge(x);
+}
+function sinc(x) {
+    const d = new decimal_js_1.default(x);
+    if (d.isZero()) {
+        return new decimal_js_1.default(1);
+    }
+    const ax = d.abs();
+    if (ax.lt(DECIMAL_SINC_SERIES_CUTOFF)) {
+        const x2 = d.mul(d);
+        return new decimal_js_1.default(1)
+            .minus(x2.div(6))
+            .plus(x2.mul(x2).div(120))
+            .minus(x2.mul(x2).mul(x2).div(5040));
+    }
+    return d.sin().div(d);
+}
+function Si(x) {
+    let d = new decimal_js_1.default(x);
+    if (d.isZero()) {
+        return new decimal_js_1.default(0);
+    }
+    let sign = new decimal_js_1.default(1);
+    if (d.isNeg()) {
+        sign = new decimal_js_1.default(-1);
+        d = d.neg();
+    }
+    function polevl(y, coef) {
+        let ans = coef[0];
+        for (let i = 1; i < coef.length; i++) {
+            ans = ans.mul(y).plus(coef[i]);
+        }
+        return ans;
+    }
+    function p1evl(y, coef) {
+        let ans = y.plus(coef[0]);
+        for (let i = 1; i < coef.length; i++) {
+            ans = ans.mul(y).plus(coef[i]);
+        }
+        return ans;
+    }
+    if (d.gt(OSCILLATORY_ASYMPTOTIC_CUTOFF)) {
+        const halfPi = decimal_js_1.default.acos(-1).div(2);
+        return sign.mul(halfPi.minus(d.cos().div(d)));
+    }
+    if (d.lte(RATIONAL_APPROXIMATION_CUTOFF)) {
+        const SN = [
+            new decimal_js_1.default('-8.39167827910303881427e-11'),
+            new decimal_js_1.default('4.62591714427012837309e-8'),
+            new decimal_js_1.default('-9.75759303843632795789e-6'),
+            new decimal_js_1.default('9.76945438170435310816e-4'),
+            new decimal_js_1.default('-4.13470316229406538752e-2'),
+            new decimal_js_1.default('1.00000000000000000302'),
+        ];
+        const SD = [
+            new decimal_js_1.default('2.03269266195951942049e-12'),
+            new decimal_js_1.default('1.27997891179943299903e-9'),
+            new decimal_js_1.default('4.41827842801218905784e-7'),
+            new decimal_js_1.default('9.96412122043875552487e-5'),
+            new decimal_js_1.default('1.42085239326149893930e-2'),
+            new decimal_js_1.default('9.99999999999999996984e-1'),
+        ];
+        const z = d.mul(d);
+        const s = d.mul(polevl(z, SN)).div(polevl(z, SD));
+        return sign.mul(s);
+    }
+    const s = d.sin();
+    const c = d.cos();
+    const z = new decimal_js_1.default(1).div(d.mul(d));
+    let f;
+    let g;
+    if (d.lt(ASYMPTOTIC_APPROXIMATION_CUTOFF)) {
+        const FN4 = [
+            new decimal_js_1.default('4.23612862892216586994'),
+            new decimal_js_1.default('5.45937717161812843388'),
+            new decimal_js_1.default('1.62083287701538329132'),
+            new decimal_js_1.default('1.67006611831323023771e-1'),
+            new decimal_js_1.default('6.81020132472518137426e-3'),
+            new decimal_js_1.default('1.08936580650328664411e-4'),
+            new decimal_js_1.default('5.48900223421373614008e-7'),
+        ];
+        const FD4 = [
+            new decimal_js_1.default('8.16496634205391016773'),
+            new decimal_js_1.default('7.30828822505564552187'),
+            new decimal_js_1.default('1.86792257950184183883'),
+            new decimal_js_1.default('1.78792052963149907262e-1'),
+            new decimal_js_1.default('7.01710668322789753610e-3'),
+            new decimal_js_1.default('1.10034357153915731354e-4'),
+            new decimal_js_1.default('5.48900252756255700982e-7'),
+        ];
+        const GN4 = [
+            new decimal_js_1.default('8.71001698973114191777e-2'),
+            new decimal_js_1.default('6.11379109952219284151e-1'),
+            new decimal_js_1.default('3.97180296392337498885e-1'),
+            new decimal_js_1.default('7.48527737628469092119e-2'),
+            new decimal_js_1.default('5.38868681462177273157e-3'),
+            new decimal_js_1.default('1.61999794598934024525e-4'),
+            new decimal_js_1.default('1.97963874140963632189e-6'),
+            new decimal_js_1.default('7.82579040744090311069e-9'),
+        ];
+        const GD4 = [
+            new decimal_js_1.default('1.64402202413355338886'),
+            new decimal_js_1.default('6.66296701268987968381e-1'),
+            new decimal_js_1.default('9.88771761277688796203e-2'),
+            new decimal_js_1.default('6.22396345441768420760e-3'),
+            new decimal_js_1.default('1.73221081474177119497e-4'),
+            new decimal_js_1.default('2.02659182086343991969e-6'),
+            new decimal_js_1.default('7.82579218933534490868e-9'),
+        ];
+        f = polevl(z, FN4).div(d.mul(p1evl(z, FD4)));
+        g = z.mul(polevl(z, GN4)).div(p1evl(z, GD4));
+    }
+    else {
+        const FN8 = [
+            new decimal_js_1.default('4.55880873470465315206e-1'),
+            new decimal_js_1.default('7.13715274100146711374e-1'),
+            new decimal_js_1.default('1.60300158222319456320e-1'),
+            new decimal_js_1.default('1.16064229408124407915e-2'),
+            new decimal_js_1.default('3.49556442447859055605e-4'),
+            new decimal_js_1.default('4.86215430826454749482e-6'),
+            new decimal_js_1.default('3.20092790091004902806e-8'),
+            new decimal_js_1.default('9.41779576128512936592e-11'),
+            new decimal_js_1.default('9.70507110881952024631e-14'),
+        ];
+        const FD8 = [
+            new decimal_js_1.default('9.17463611873684053703e-1'),
+            new decimal_js_1.default('1.78685545332074536321e-1'),
+            new decimal_js_1.default('1.22253594771971293032e-2'),
+            new decimal_js_1.default('3.58696481881851580297e-4'),
+            new decimal_js_1.default('4.92435064317881464393e-6'),
+            new decimal_js_1.default('3.21956939101046018377e-8'),
+            new decimal_js_1.default('9.43720590350276732376e-11'),
+            new decimal_js_1.default('9.70507110881952025725e-14'),
+        ];
+        const GN8 = [
+            new decimal_js_1.default('6.97359953443276214934e-1'),
+            new decimal_js_1.default('3.30410979305632063225e-1'),
+            new decimal_js_1.default('3.84878767649974295920e-2'),
+            new decimal_js_1.default('1.71718239052347903558e-3'),
+            new decimal_js_1.default('3.48941165502279436777e-5'),
+            new decimal_js_1.default('3.47131167084116673800e-7'),
+            new decimal_js_1.default('1.70404452782044526189e-9'),
+            new decimal_js_1.default('3.85945925430276600453e-12'),
+            new decimal_js_1.default('3.14040098946363334640e-15'),
+        ];
+        const GD8 = [
+            new decimal_js_1.default('1.68548898811011640017'),
+            new decimal_js_1.default('4.87852258695304967486e-1'),
+            new decimal_js_1.default('4.67913194259625806320e-2'),
+            new decimal_js_1.default('1.90284426674399523638e-3'),
+            new decimal_js_1.default('3.68475504442561108162e-5'),
+            new decimal_js_1.default('3.57043223443740838771e-7'),
+            new decimal_js_1.default('1.72693748966316146736e-9'),
+            new decimal_js_1.default('3.87830166023954706752e-12'),
+            new decimal_js_1.default('3.14040098946363335242e-15'),
+        ];
+        f = polevl(z, FN8).div(d.mul(p1evl(z, FD8)));
+        g = z.mul(polevl(z, GN8)).div(p1evl(z, GD8));
+    }
+    const halfPi = decimal_js_1.default.acos(-1).div(2);
+    return sign.mul(halfPi.minus(f.mul(c)).minus(g.mul(s)));
+}
+function Shi(x) {
+    const d = new decimal_js_1.default(x);
+    if (d.isZero()) {
+        return new decimal_js_1.default(0);
+    }
+    const sign = d.isNeg() ? new decimal_js_1.default(-1) : new decimal_js_1.default(1);
+    const ax = d.abs();
+    if (ax.lt(ASYMPTOTIC_APPROXIMATION_CUTOFF)) {
+        let term = ax;
+        let sum = term;
+        let n = 0;
+        for (let iter = 0; iter < SERIES_MAX_ITERATIONS; iter++) {
+            const a = 2 * n + 1;
+            term = term
+                .mul(ax)
+                .mul(ax)
+                .mul(a)
+                .div(new decimal_js_1.default(a + 2).mul(a + 2).mul(a + 1));
+            sum = sum.plus(term);
+            n++;
+            if (term.abs().lt(getSeriesError())) {
+                break;
+            }
+        }
+        return sign.mul(sum);
+    }
+    function EiLocal(y) {
+        if (y.isZero()) {
+            return new decimal_js_1.default(-Infinity);
+        }
+        if (y.gt(EXPONENTIAL_ASYMPTOTIC_CUTOFF) || y.lt(-EXPONENTIAL_ASYMPTOTIC_CUTOFF)) {
+            let term = new decimal_js_1.default(1);
+            let sum = new decimal_js_1.default(1);
+            let best = term.abs();
+            for (let n = 1; n < ASYMPTOTIC_MAX_TERMS; n++) {
+                term = term.mul(n).div(y);
+                const at = term.abs();
+                if (at.gt(best)) {
+                    break;
+                }
+                best = at;
+                sum = sum.plus(term);
+                if (at.lt(getSeriesError())) {
+                    break;
+                }
+            }
+            return y.exp().mul(sum).div(y);
+        }
+        let term = y;
+        let sum = term;
+        for (let k = 1; k < POWER_SERIES_MAX_TERMS; k++) {
+            term = term
+                .mul(y)
+                .mul(k)
+                .div(new decimal_js_1.default(k + 1).mul(k + 1));
+            sum = sum.plus(term);
+            if (term.abs().lt(getSeriesError())) {
+                break;
+            }
+        }
+        return exports.EULER_GAMMA.plus(y.abs().ln()).plus(sum);
+    }
+    return EiLocal(d).minus(EiLocal(d.neg())).div(2);
+}
+function Ci(x) {
+    let d = new decimal_js_1.default(x);
+    if (d.isZero()) {
+        return new decimal_js_1.default(-Infinity);
+    }
+    if (d.isNeg()) {
+        d = d.neg();
+    }
+    function polevl(y, coef) {
+        let ans = coef[0];
+        for (let i = 1; i < coef.length; i++) {
+            ans = ans.mul(y).plus(coef[i]);
+        }
+        return ans;
+    }
+    function p1evl(y, coef) {
+        let ans = y.plus(coef[0]);
+        for (let i = 1; i < coef.length; i++) {
+            ans = ans.mul(y).plus(coef[i]);
+        }
+        return ans;
+    }
+    if (d.gt(OSCILLATORY_ASYMPTOTIC_CUTOFF)) {
+        return d.sin().div(d);
+    }
+    if (d.lte(RATIONAL_APPROXIMATION_CUTOFF)) {
+        const CN = [
+            new decimal_js_1.default('2.02524002389102268789e-11'),
+            new decimal_js_1.default('-1.35249504915790756375e-8'),
+            new decimal_js_1.default('3.59325051419993077021e-6'),
+            new decimal_js_1.default('-4.74007206873407909465e-4'),
+            new decimal_js_1.default('2.89159652607555242092e-2'),
+            new decimal_js_1.default('-1.00000000000000000080'),
+        ];
+        const CD = [
+            new decimal_js_1.default('4.07746040061880559506e-12'),
+            new decimal_js_1.default('3.06780997581887812692e-9'),
+            new decimal_js_1.default('1.23210355685883423679e-6'),
+            new decimal_js_1.default('3.17442024775032769882e-4'),
+            new decimal_js_1.default('5.10028056236446052392e-2'),
+            new decimal_js_1.default('4.00000000000000000080'),
+        ];
+        const z = d.mul(d);
+        const c0 = z.mul(polevl(z, CN)).div(polevl(z, CD));
+        return exports.EULER_GAMMA.plus(d.ln()).plus(c0);
+    }
+    const s = d.sin();
+    const c = d.cos();
+    const z = new decimal_js_1.default(1).div(d.mul(d));
+    let f;
+    let g;
+    if (d.lt(ASYMPTOTIC_APPROXIMATION_CUTOFF)) {
+        const FN4 = [
+            new decimal_js_1.default('4.23612862892216586994'),
+            new decimal_js_1.default('5.45937717161812843388'),
+            new decimal_js_1.default('1.62083287701538329132'),
+            new decimal_js_1.default('1.67006611831323023771e-1'),
+            new decimal_js_1.default('6.81020132472518137426e-3'),
+            new decimal_js_1.default('1.08936580650328664411e-4'),
+            new decimal_js_1.default('5.48900223421373614008e-7'),
+        ];
+        const FD4 = [
+            new decimal_js_1.default('8.16496634205391016773'),
+            new decimal_js_1.default('7.30828822505564552187'),
+            new decimal_js_1.default('1.86792257950184183883'),
+            new decimal_js_1.default('1.78792052963149907262e-1'),
+            new decimal_js_1.default('7.01710668322789753610e-3'),
+            new decimal_js_1.default('1.10034357153915731354e-4'),
+            new decimal_js_1.default('5.48900252756255700982e-7'),
+        ];
+        const GN4 = [
+            new decimal_js_1.default('8.71001698973114191777e-2'),
+            new decimal_js_1.default('6.11379109952219284151e-1'),
+            new decimal_js_1.default('3.97180296392337498885e-1'),
+            new decimal_js_1.default('7.48527737628469092119e-2'),
+            new decimal_js_1.default('5.38868681462177273157e-3'),
+            new decimal_js_1.default('1.61999794598934024525e-4'),
+            new decimal_js_1.default('1.97963874140963632189e-6'),
+            new decimal_js_1.default('7.82579040744090311069e-9'),
+        ];
+        const GD4 = [
+            new decimal_js_1.default('1.64402202413355338886'),
+            new decimal_js_1.default('6.66296701268987968381e-1'),
+            new decimal_js_1.default('9.88771761277688796203e-2'),
+            new decimal_js_1.default('6.22396345441768420760e-3'),
+            new decimal_js_1.default('1.73221081474177119497e-4'),
+            new decimal_js_1.default('2.02659182086343991969e-6'),
+            new decimal_js_1.default('7.82579218933534490868e-9'),
+        ];
+        f = polevl(z, FN4).div(d.mul(p1evl(z, FD4)));
+        g = z.mul(polevl(z, GN4)).div(p1evl(z, GD4));
+    }
+    else {
+        const FN8 = [
+            new decimal_js_1.default('4.55880873470465315206e-1'),
+            new decimal_js_1.default('7.13715274100146711374e-1'),
+            new decimal_js_1.default('1.60300158222319456320e-1'),
+            new decimal_js_1.default('1.16064229408124407915e-2'),
+            new decimal_js_1.default('3.49556442447859055605e-4'),
+            new decimal_js_1.default('4.86215430826454749482e-6'),
+            new decimal_js_1.default('3.20092790091004902806e-8'),
+            new decimal_js_1.default('9.41779576128512936592e-11'),
+            new decimal_js_1.default('9.70507110881952024631e-14'),
+        ];
+        const FD8 = [
+            new decimal_js_1.default('9.17463611873684053703e-1'),
+            new decimal_js_1.default('1.78685545332074536321e-1'),
+            new decimal_js_1.default('1.22253594771971293032e-2'),
+            new decimal_js_1.default('3.58696481881851580297e-4'),
+            new decimal_js_1.default('4.92435064317881464393e-6'),
+            new decimal_js_1.default('3.21956939101046018377e-8'),
+            new decimal_js_1.default('9.43720590350276732376e-11'),
+            new decimal_js_1.default('9.70507110881952025725e-14'),
+        ];
+        const GN8 = [
+            new decimal_js_1.default('6.97359953443276214934e-1'),
+            new decimal_js_1.default('3.30410979305632063225e-1'),
+            new decimal_js_1.default('3.84878767649974295920e-2'),
+            new decimal_js_1.default('1.71718239052347903558e-3'),
+            new decimal_js_1.default('3.48941165502279436777e-5'),
+            new decimal_js_1.default('3.47131167084116673800e-7'),
+            new decimal_js_1.default('1.70404452782044526189e-9'),
+            new decimal_js_1.default('3.85945925430276600453e-12'),
+            new decimal_js_1.default('3.14040098946363334640e-15'),
+        ];
+        const GD8 = [
+            new decimal_js_1.default('1.68548898811011640017'),
+            new decimal_js_1.default('4.87852258695304967486e-1'),
+            new decimal_js_1.default('4.67913194259625806320e-2'),
+            new decimal_js_1.default('1.90284426674399523638e-3'),
+            new decimal_js_1.default('3.68475504442561108162e-5'),
+            new decimal_js_1.default('3.57043223443740838771e-7'),
+            new decimal_js_1.default('1.72693748966316146736e-9'),
+            new decimal_js_1.default('3.87830166023954706752e-12'),
+            new decimal_js_1.default('3.14040098946363335242e-15'),
+        ];
+        f = polevl(z, FN8).div(d.mul(p1evl(z, FD8)));
+        g = z.mul(polevl(z, GN8)).div(p1evl(z, GD8));
+    }
+    return f.mul(s).minus(g.mul(c));
+}
+function Chi(x) {
+    const d = new decimal_js_1.default(x);
+    const ax = d.abs();
+    if (ax.isZero()) {
+        return new decimal_js_1.default(-Infinity);
+    }
+    if (ax.lt(ASYMPTOTIC_APPROXIMATION_CUTOFF)) {
+        let n = 1;
+        let term = ax.mul(ax).div(4);
+        let sum = term;
+        for (let iter = 0; iter < SERIES_MAX_ITERATIONS; iter++) {
+            term = term
+                .mul(ax)
+                .mul(ax)
+                .mul(n)
+                .div(new decimal_js_1.default(n + 1).mul(2 * n + 1).mul(2 * n + 2));
+            sum = sum.plus(term);
+            n++;
+            if (term.abs().lt(getSeriesError())) {
+                break;
+            }
+        }
+        return exports.EULER_GAMMA.plus(ax.ln()).plus(sum);
+    }
+    function EiLocal(y) {
+        if (y.isZero()) {
+            return new decimal_js_1.default(-Infinity);
+        }
+        if (y.gt(EXPONENTIAL_ASYMPTOTIC_CUTOFF) || y.lt(-EXPONENTIAL_ASYMPTOTIC_CUTOFF)) {
+            let term = new decimal_js_1.default(1);
+            let sum = new decimal_js_1.default(1);
+            let best = term.abs();
+            for (let n = 1; n < ASYMPTOTIC_MAX_TERMS; n++) {
+                term = term.mul(n).div(y);
+                const at = term.abs();
+                if (at.gt(best)) {
+                    break;
+                }
+                best = at;
+                sum = sum.plus(term);
+                if (at.lt(getSeriesError())) {
+                    break;
+                }
+            }
+            return y.exp().mul(sum).div(y);
+        }
+        let term = y;
+        let sum = term;
+        for (let k = 1; k < POWER_SERIES_MAX_TERMS; k++) {
+            term = term
+                .mul(y)
+                .mul(k)
+                .div(new decimal_js_1.default(k + 1).mul(k + 1));
+            sum = sum.plus(term);
+            if (term.abs().lt(getSeriesError())) {
+                break;
+            }
+        }
+        return exports.EULER_GAMMA.plus(y.abs().ln()).plus(sum);
+    }
+    return EiLocal(d).plus(EiLocal(d.neg())).div(2);
+}
+function Ei(x) {
+    const d = new decimal_js_1.default(x);
+    if (d.isZero()) {
+        return new decimal_js_1.default(-Infinity);
+    }
+    if (d.gt(EXPONENTIAL_ASYMPTOTIC_CUTOFF) || d.lt(-EXPONENTIAL_ASYMPTOTIC_CUTOFF)) {
+        let term = new decimal_js_1.default(1);
+        let sum = new decimal_js_1.default(1);
+        let best = term.abs();
+        for (let n = 1; n < ASYMPTOTIC_MAX_TERMS; n++) {
+            term = term.mul(n).div(d);
+            const at = term.abs();
+            if (at.gt(best)) {
+                break;
+            }
+            best = at;
+            sum = sum.plus(term);
+            if (at.lt(getSeriesError())) {
+                break;
+            }
+        }
+        return d.exp().mul(sum).div(d);
+    }
+    let term = d;
+    let sum = term;
+    for (let k = 1; k < POWER_SERIES_MAX_TERMS; k++) {
+        term = term
+            .mul(d)
+            .mul(k)
+            .div(new decimal_js_1.default(k + 1).mul(k + 1));
+        sum = sum.plus(term);
+        if (term.abs().lt(getSeriesError())) {
+            break;
+        }
+    }
+    return exports.EULER_GAMMA.plus(d.abs().ln()).plus(sum);
+}
+function Li(x) {
+    const d = new decimal_js_1.default(x);
+    if (d.lte(0)) {
+        return new decimal_js_1.default(NaN);
+    }
+    if (d.eq(1)) {
+        return new decimal_js_1.default(-Infinity);
+    }
+    const y = d.ln();
+    if (y.isZero()) {
+        return new decimal_js_1.default(-Infinity);
+    }
+    if (y.gt(EXPONENTIAL_ASYMPTOTIC_CUTOFF) || y.lt(-EXPONENTIAL_ASYMPTOTIC_CUTOFF)) {
+        let term = new decimal_js_1.default(1);
+        let sum = new decimal_js_1.default(1);
+        let best = term.abs();
+        for (let n = 1; n < ASYMPTOTIC_MAX_TERMS; n++) {
+            term = term.mul(n).div(y);
+            const at = term.abs();
+            if (at.gt(best)) {
+                break;
+            }
+            best = at;
+            sum = sum.plus(term);
+            if (at.lt(getSeriesError())) {
+                break;
+            }
+        }
+        return y.exp().mul(sum).div(y);
+    }
+    let term = y;
+    let sum = term;
+    for (let k = 1; k < POWER_SERIES_MAX_TERMS; k++) {
+        term = term
+            .mul(y)
+            .mul(k)
+            .div(new decimal_js_1.default(k + 1).mul(k + 1));
+        sum = sum.plus(term);
+        if (term.abs().lt(getSeriesError())) {
+            break;
+        }
+    }
+    return exports.EULER_GAMMA.plus(y.abs().ln()).plus(sum);
+}