goscript 0.0.34 → 0.0.36

package/dist/gs/math/lgamma.gs.js

@@ -0,0 +1,243 @@
+import * as $ from "@goscript/builtin/builtin.js";
+import { Abs } from "./abs.gs.js";
+import { Inf, IsInf, IsNaN } from "./bits.gs.js";
+import { Floor } from "./floor.gs.js";
+import { Log } from "./log.gs.js";
+import { Mod } from "./mod.gs.js";
+import { Cos, Sin } from "./sin.gs.js";
+import { Float64bits } from "./unsafe.gs.js";
+let _lgamA = $.arrayToSlice([7.72156649015328655494e-02, 3.22467033424113591611e-01, 6.73523010531292681824e-02, 2.05808084325167332806e-02, 7.38555086081402883957e-03, 2.89051383673415629091e-03, 1.19270763183362067845e-03, 5.10069792153511336608e-04, 2.20862790713908385557e-04, 1.08011567247583939954e-04, 2.52144565451257326939e-05, 4.48640949618915160150e-05]);
+let _lgamR = $.arrayToSlice([1.0, 1.39200533467621045958e+00, 7.21935547567138069525e-01, 1.71933865632803078993e-01, 1.86459191715652901344e-02, 7.77942496381893596434e-04, 7.32668430744625636189e-06]);
+let _lgamS = $.arrayToSlice([-7.72156649015328655494e-02, 2.14982415960608852501e-01, 3.25778796408930981787e-01, 1.46350472652464452805e-01, 2.66422703033638609560e-02, 1.84028451407337715652e-03, 3.19475326584100867617e-05]);
+let _lgamT = $.arrayToSlice([4.83836122723810047042e-01, -1.47587722994593911752e-01, 6.46249402391333854778e-02, -3.27885410759859649565e-02, 1.79706750811820387126e-02, -1.03142241298341437450e-02, 6.10053870246291332635e-03, -3.68452016781138256760e-03, 2.25964780900612472250e-03, -1.40346469989232843813e-03, 8.81081882437654011382e-04, -5.38595305356740546715e-04, 3.15632070903625950361e-04, -3.12754168375120860518e-04, 3.35529192635519073543e-04]);
+let _lgamU = $.arrayToSlice([-7.72156649015328655494e-02, 6.32827064025093366517e-01, 1.45492250137234768737e+00, 9.77717527963372745603e-01, 2.28963728064692451092e-01, 1.33810918536787660377e-02]);
+let _lgamV = $.arrayToSlice([1.0, 2.45597793713041134822e+00, 2.12848976379893395361e+00, 7.69285150456672783825e-01, 1.04222645593369134254e-01, 3.21709242282423911810e-03]);
+let _lgamW = $.arrayToSlice([4.18938533204672725052e-01, 8.33333333333329678849e-02, -2.77777777728775536470e-03, 7.93650558643019558500e-04, -5.95187557450339963135e-04, 8.36339918996282139126e-04, -1.63092934096575273989e-03]);
+// Lgamma returns the natural logarithm and sign (-1 or +1) of [Gamma](x).
+//
+// Special cases are:
+//
+//	Lgamma(+Inf) = +Inf
+//	Lgamma(0) = +Inf
+//	Lgamma(-integer) = +Inf
+//	Lgamma(-Inf) = -Inf
+//	Lgamma(NaN) = NaN
+export function Lgamma(x) {
+    let lgamma = 0;
+    let sign = 1;
+    // Constants
+    const Ymin = 1.461632144968362245;
+    const Two52 = (1 << 52);
+    const Two53 = (1 << 53);
+    const Two58 = (1 << 58);
+    const Tiny = 1.0 / ((1 << 70));
+    const Tc = 1.46163214496836224576e+00;
+    const Tf = -1.21486290535849611461e-01;
+    const Tt = -3.63867699703950536541e-18;
+    // special cases
+    if (IsNaN(x)) {
+        lgamma = x;
+        return [lgamma, sign];
+    }
+    if (IsInf(x, 0)) {
+        lgamma = x;
+        return [lgamma, sign];
+    }
+    if (x == 0) {
+        lgamma = Inf(1);
+        return [lgamma, sign];
+    }
+    let neg = false;
+    if (x < 0) {
+        x = -x;
+        neg = true;
+    }
+    // if |x| < 2**-70, return -log(|x|)
+    if (x < 8.47033e-22) {
+        if (neg) {
+            sign = -1;
+        }
+        lgamma = -Log(x);
+        return [lgamma, sign];
+    }
+    let nadj = 0;
+    // Handle negative values
+    if (neg) {
+        // |x| >= 2**52, must be -integer
+        if (x >= 4503599627370496) {
+            lgamma = Inf(1);
+            return [lgamma, sign];
+        }
+        let t = sinPi(x);
+        // -integer
+        if (t == 0) {
+            lgamma = Inf(1);
+            return [lgamma, sign];
+        }
+        nadj = Log(3.14159 / Abs(t * x));
+        if (t < 0) {
+            sign = -1;
+        }
+    }
+    // Main computation
+    if (x == 1 || x == 2) {
+        lgamma = 0;
+    }
+    else if (x < 2) {
+        let y = 0;
+        let i = 0;
+        if (x <= 0.9) {
+            lgamma = -Log(x);
+            if (x >= (1.46163 - 1 + 0.27)) {
+                y = 1 - x;
+                i = 0;
+            }
+            else if (x >= (1.46163 - 1 - 0.27)) {
+                y = x - (1.46163 - 1);
+                i = 1;
+            }
+            else {
+                y = x;
+                i = 2;
+            }
+        }
+        else {
+            lgamma = 0;
+            if (x >= (1.46163 + 0.27)) {
+                y = 2 - x;
+                i = 0;
+            }
+            else if (x >= (1.46163 - 0.27)) {
+                y = x - 1.46163;
+                i = 1;
+            }
+            else {
+                y = x - 1;
+                i = 2;
+            }
+        }
+        if (i === 0) {
+            let z = y * y;
+            let p1 = _lgamA[0] + z * (_lgamA[2] + z * (_lgamA[4] + z * (_lgamA[6] + z * (_lgamA[8] + z * _lgamA[10]))));
+            let p2 = z * (_lgamA[1] + z * (+_lgamA[3] + z * (_lgamA[5] + z * (_lgamA[7] + z * (_lgamA[9] + z * _lgamA[11])))));
+            let p = y * p1 + p2;
+            lgamma += (p - 0.5 * y);
+        }
+        else if (i === 1) {
+            let z = y * y;
+            let w = z * y;
+            let p1 = _lgamT[0] + w * (_lgamT[3] + w * (_lgamT[6] + w * (_lgamT[9] + w * _lgamT[12])));
+            let p2 = _lgamT[1] + w * (_lgamT[4] + w * (_lgamT[7] + w * (_lgamT[10] + w * _lgamT[13])));
+            let p3 = _lgamT[2] + w * (_lgamT[5] + w * (_lgamT[8] + w * (_lgamT[11] + w * _lgamT[14])));
+            let p = z * p1 - (-3.63868e-18 - w * (p2 + y * p3));
+            lgamma += (-0.121486 + p);
+        }
+        else { // i === 2
+            let p1 = y * (_lgamU[0] + y * (_lgamU[1] + y * (_lgamU[2] + y * (_lgamU[3] + y * (_lgamU[4] + y * _lgamU[5])))));
+            let p2 = 1 + y * (_lgamV[1] + y * (_lgamV[2] + y * (_lgamV[3] + y * (_lgamV[4] + y * _lgamV[5]))));
+            lgamma += (-0.5 * y + p1 / p2);
+        }
+    }
+    else if (x < 8) {
+        let i = $.int(x);
+        let y = x - i;
+        let p = y * (_lgamS[0] + y * (_lgamS[1] + y * (_lgamS[2] + y * (_lgamS[3] + y * (_lgamS[4] + y * (_lgamS[5] + y * _lgamS[6]))))));
+        let q = 1 + y * (_lgamR[1] + y * (_lgamR[2] + y * (_lgamR[3] + y * (_lgamR[4] + y * (_lgamR[5] + y * _lgamR[6])))));
+        lgamma = 0.5 * y + p / q;
+        let z = 1.0;
+        // Handle fallthrough cases properly
+        if (i === 7) {
+            z *= (y + 6);
+            z *= (y + 5);
+            z *= (y + 4);
+            z *= (y + 3);
+            z *= (y + 2);
+            lgamma += Log(z);
+        }
+        else if (i === 6) {
+            z *= (y + 5);
+            z *= (y + 4);
+            z *= (y + 3);
+            z *= (y + 2);
+            lgamma += Log(z);
+        }
+        else if (i === 5) {
+            z *= (y + 4);
+            z *= (y + 3);
+            z *= (y + 2);
+            lgamma += Log(z);
+        }
+        else if (i === 4) {
+            z *= (y + 3);
+            z *= (y + 2);
+            lgamma += Log(z);
+        }
+        else if (i === 3) {
+            z *= (y + 2);
+            lgamma += Log(z);
+        }
+    }
+    else if (x < 288230376151711744) {
+        let t = Log(x);
+        let z = 1 / x;
+        let y = z * z;
+        let w = _lgamW[0] + z * (_lgamW[1] + y * (_lgamW[2] + y * (_lgamW[3] + y * (_lgamW[4] + y * (_lgamW[5] + y * _lgamW[6])))));
+        lgamma = (x - 0.5) * (t - 1) + w;
+    }
+    else {
+        lgamma = x * (Log(x) - 1);
+    }
+    if (neg) {
+        lgamma = nadj - lgamma;
+    }
+    return [lgamma, sign];
+}
+// sinPi(x) is a helper function for negative x
+export function sinPi(x) {
+    const Two52 = (1 << 52);
+    const Two53 = (1 << 53);
+    if (x < 0.25) {
+        return -Sin(3.14159 * x);
+    }
+    // argument reduction
+    let z = Floor(x);
+    let n = 0;
+    if (z != x) {
+        // inexact
+        x = Mod(x, 2);
+        n = $.int(x * 4);
+    }
+    else {
+        // exact
+        if (x >= 9007199254740992) {
+            // x must be even
+            x = 0;
+            n = 0;
+        }
+        else {
+            if (x < 4503599627370496) {
+                z = x + 4503599627370496; // exact
+            }
+            n = $.int((1 & Float64bits(z)));
+            x = n;
+            n <<= 2;
+        }
+    }
+    if (n === 0) {
+        x = Sin(3.14159 * x);
+    }
+    else if (n === 1 || n === 2) {
+        x = Cos(3.14159 * (0.5 - x));
+    }
+    else if (n === 3 || n === 4) {
+        x = Sin(3.14159 * (1 - x));
+    }
+    else if (n === 5 || n === 6) {
+        x = -Cos(3.14159 * (x - 1.5));
+    }
+    else {
+        x = Sin(3.14159 * (x - 2));
+    }
+    return -x;
+}
+//# sourceMappingURL=lgamma.gs.js.map

package/dist/gs/math/lgamma.gs.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"lgamma.gs.js","sourceRoot":"","sources":["../../../gs/math/lgamma.gs.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,CAAC,MAAM,8BAA8B,CAAC;AAClD,OAAO,EAAE,GAAG,EAAE,MAAM,aAAa,CAAC;AAClC,OAAO,EAAE,GAAG,EAAE,KAAK,EAAE,KAAK,EAAE,MAAM,cAAc,CAAC;AACjD,OAAO,EAAE,KAAK,EAAE,MAAM,eAAe,CAAC;AACtC,OAAO,EAAE,GAAG,EAAE,MAAM,aAAa,CAAC;AAClC,OAAO,EAAE,GAAG,EAAE,MAAM,aAAa,CAAC;AAClC,OAAO,EAAE,GAAG,EAAE,GAAG,EAAE,MAAM,aAAa,CAAC;AACvC,OAAO,EAAE,WAAW,EAAE,MAAM,gBAAgB,CAAC;AAE7C,IAAI,MAAM,GAAG,CAAC,CAAC,YAAY,CAAS,CAAC,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,CAAC,CAAC,CAAA;AAErX,IAAI,MAAM,GAAG,CAAC,CAAC,YAAY,CAAS,CAAC,GAAG,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,CAAC,CAAC,CAAA;AAElN,IAAI,MAAM,GAAG,CAAC,CAAC,YAAY,CAAS,CAAC,CAAC,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,CAAC,CAAC,CAAA;AAE1O,IAAI,MAAM,GAAG,CAAC,CAAC,YAAY,CAAS,CAAC,0BAA0B,EAAE,CAAC,0BAA0B,EAAE,0BAA0B,EAAE,CAAC,0BAA0B,EAAE,0BAA0B,EAAE,CAAC,0BAA0B,EAAE,0BAA0B,EAAE,CAAC,0BAA0B,EAAE,0BAA0B,EAAE,CAAC,0BAA0B,EAAE,0BAA0B,EAAE,CAAC,0BAA0B,EAAE,0BAA0B,EAAE,CAAC,0BAA0B,EAAE,0BAA0B,CAAC,CAAC,CAAA;AAEhd,IAAI,MAAM,GAAG,CAAC,CAAC,YAAY,CAAS,CAAC,CAAC,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,CAAC,CAAC,CAAA;AAE9M,IAAI,MAAM,GAAG,CAAC,CAAC,YAAY,CAAS,CAAC,GAAG,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,EAAE,0BAA0B,CAAC,CAAC,CAAA;AAEtL,IAAI,MAAM,GAAG,CAAC,CAAC,YAAY,CAAS,CAAC,0BAA0B,EAAE,0BAA0B,EAAE,CAAC,0BAA0B,EAAE,0BAA0B,EAAE,CAAC,0BAA0B,EAAE,0BAA0B,EAAE,CAAC,0BAA0B,CAAC,CAAC,CAAA;AAE5O,0EAA0E;AAC1E,EAAE;AACF,qBAAqB;AACrB,EAAE;AACF,sBAAsB;AACtB,mBAAmB;AACnB,0BAA0B;AAC1B,sBAAsB;AACtB,oBAAoB;AACpB,MAAM,UAAU,MAAM,CAAC,CAAS;IAC/B,IAAI,MAAM,GAAW,CAAC,CAAC;IACvB,IAAI,IAAI,GAAW,CAAC,CAAC;IAErB,YAAY;IACZ,MAAM,IAAI,GAAW,oBAAoB,CAAC;IAC1C,MAAM,KAAK,GAAW,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC;IAChC,MAAM,KAAK,GAAW,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC;IAChC,MAAM,KAAK,GAAW,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC;IAChC,MAAM,IAAI,GAAW,GAAG,GAAG,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;IACvC,MAAM,EAAE,GAAW,0BAA0B,CAAC;IAC9C,MAAM,EAAE,GAAW,CAAC,0BAA0B,CAAC;IAC/C,MAAM,EAAE,GAAW,CAAC,0BAA0B,CAAC;IAE/C,gBAAgB;IAChB,IAAI,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC;QACd,MAAM,GAAG,CAAC,CAAC;QACX,OAAO,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC;IACvB,CAAC;IACD,IAAI,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC;QACjB,MAAM,GAAG,CAAC,CAAC;QACX,OAAO,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC;IACvB,CAAC;IACD,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QACZ,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QAChB,OAAO,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC;IACvB,CAAC;IAED,IAAI,GAAG,GAAG,KAAK,CAAC;IAChB,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QACX,CAAC,GAAG,CAAC,CAAC,CAAC;QACP,GAAG,GAAG,IAAI,CAAC;IACZ,CAAC;IAED,oCAAoC;IACpC,IAAI,CAAC,GAAG,WAAW,EAAE,CAAC;QACrB,IAAI,GAAG,EAAE,CAAC;YACT,IAAI,GAAG,CAAC,CAAC,CAAC;QACX,CAAC;QACD,MAAM,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACjB,OAAO,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC;IACvB,CAAC;IAED,IAAI,IAAI,GAAW,CAAC,CAAC;IAErB,yBAAyB;IACzB,IAAI,GAAG,EAAE,CAAC;QACT,iCAAiC;QACjC,IAAI,CAAC,IAAI,gBAAgB,EAAE,CAAC;YAC3B,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;YAChB,OAAO,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC;QACvB,CAAC;QACD,IAAI,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QAEjB,WAAW;QACX,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;YACZ,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;YAChB,OAAO,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC;QACvB,CAAC;QACD,IAAI,GAAG,GAAG,CAAC,OAAO,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACjC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;YACX,IAAI,GAAG,CAAC,CAAC,CAAC;QACX,CAAC;IACF,CAAC;IAED,mBAAmB;IACnB,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QACtB,MAAM,GAAG,CAAC,CAAC;IACZ,CAAC;SAAM,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QAClB,IAAI,CAAC,GAAW,CAAC,CAAC;QAClB,IAAI,CAAC,GAAW,CAAC,CAAC;QAElB,IAAI,CAAC,IAAI,GAAG,EAAE,CAAC;YACd,MAAM,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAEjB,IAAI,CAAC,IAAI,CAAC,OAAO,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;gBAC/B,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;gBACV,CAAC,GAAG,CAAC,CAAC;YACP,CAAC;iBAAM,IAAI,CAAC,IAAI,CAAC,OAAO,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;gBACtC,CAAC,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,CAAC,CAAC;gBACtB,CAAC,GAAG,CAAC,CAAC;YACP,CAAC;iBAAM,CAAC;gBACP,CAAC,GAAG,CAAC,CAAC;gBACN,CAAC,GAAG,CAAC,CAAC;YACP,CAAC;QACF,CAAC;aAAM,CAAC;YACP,MAAM,GAAG,CAAC,CAAC;YAEX,IAAI,CAAC,IAAI,CAAC,OAAO,GAAG,IAAI,CAAC,EAAE,CAAC;gBAC3B,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;gBACV,CAAC,GAAG,CAAC,CAAC;YACP,CAAC;iBAAM,IAAI,CAAC,IAAI,CAAC,OAAO,GAAG,IAAI,CAAC,EAAE,CAAC;gBAClC,CAAC,GAAG,CAAC,GAAG,OAAO,CAAC;gBAChB,CAAC,GAAG,CAAC,CAAC;YACP,CAAC;iBAAM,CAAC;gBACP,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;gBACV,CAAC,GAAG,CAAC,CAAC;YACP,CAAC;QACF,CAAC;QAED,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;YACb,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YACd,IAAI,EAAE,GAAG,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,MAAO,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAClH,IAAI,EAAE,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,MAAO,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACzH,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC;YACpB,MAAM,IAAI,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC;QACzB,CAAC;aAAM,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;YACpB,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YACd,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YACd,IAAI,EAAE,GAAG,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,MAAO,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;YAC/F,IAAI,EAAE,GAAG,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,MAAO,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;YAChG,IAAI,EAAE,GAAG,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,MAAO,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;YAChG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,WAAW,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;YACpD,MAAM,IAAI,CAAC,CAAC,QAAQ,GAAG,CAAC,CAAC,CAAC;QAC3B,CAAC;aAAM,CAAC,CAAC,UAAU;YAClB,IAAI,EAAE,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,MAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACvH,IAAI,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,MAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACxG,MAAM,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;QAChC,CAAC;IACF,CAAC;SAAM,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QAClB,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACjB,IAAI,CAAC,GAAG,CAAC,GAAI,CAAY,CAAC;QAC1B,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,MAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACzI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,MAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC1H,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACzB,IAAI,CAAC,GAAG,GAAG,CAAC;QAEZ,oCAAoC;QACpC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;YACb,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACb,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACb,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACb,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACb,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACb,MAAM,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC;QAClB,CAAC;aAAM,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;YACpB,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACb,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACb,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACb,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACb,MAAM,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC;QAClB,CAAC;aAAM,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;YACpB,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACb,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACb,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACb,MAAM,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC;QAClB,CAAC;aAAM,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;YACpB,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACb,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACb,MAAM,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC;QAClB,CAAC;aAAM,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;YACpB,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACb,MAAM,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC;QAClB,CAAC;IACF,CAAC;SAAM,IAAI,CAAC,GAAG,kBAAkB,EAAE,CAAC;QACnC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QACf,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACd,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACd,IAAI,CAAC,GAAG,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,MAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACnI,MAAM,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;IAClC,CAAC;SAAM,CAAC;QACP,MAAM,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC3B,CAAC;IAED,IAAI,GAAG,EAAE,CAAC;QACT,MAAM,GAAG,IAAI,GAAG,MAAM,CAAC;IACxB,CAAC;IAED,OAAO,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC;AACvB,CAAC;AAED,+CAA+C;AAC/C,MAAM,UAAU,KAAK,CAAC,CAAS;IAC9B,MAAM,KAAK,GAAW,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC;IAChC,MAAM,KAAK,GAAW,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC;IAEhC,IAAI,CAAC,GAAG,IAAI,EAAE,CAAC;QACd,OAAO,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,CAAC,CAAC;IAC1B,CAAC;IAED,qBAAqB;IACrB,IAAI,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IACjB,IAAI,CAAC,GAAW,CAAC,CAAC;IAElB,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QACZ,UAAU;QACV,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACd,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAClB,CAAC;SAAM,CAAC;QACP,QAAQ;QACR,IAAI,CAAC,IAAI,gBAAgB,EAAE,CAAC;YAC3B,iBAAiB;YACjB,CAAC,GAAG,CAAC,CAAC;YACN,CAAC,GAAG,CAAC,CAAC;QACP,CAAC;aAAM,CAAC;YACP,IAAI,CAAC,GAAG,gBAAgB,EAAE,CAAC;gBAC1B,CAAC,GAAG,CAAC,GAAG,gBAAgB,CAAC,CAAC,QAAQ;YACnC,CAAC;YACD,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAChC,CAAC,GAAI,CAAY,CAAC;YAClB,CAAC,KAAK,CAAC,CAAC;QACT,CAAC;IACF,CAAC;IAED,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QACb,CAAC,GAAG,GAAG,CAAC,OAAO,GAAG,CAAC,CAAC,CAAC;IACtB,CAAC;SAAM,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QAC/B,CAAC,GAAG,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IAC9B,CAAC;SAAM,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QAC/B,CAAC,GAAG,GAAG,CAAC,OAAO,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAC5B,CAAC;SAAM,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QAC/B,CAAC,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;IAC/B,CAAC;SAAM,CAAC;QACP,CAAC,GAAG,GAAG,CAAC,OAAO,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAC5B,CAAC;IAED,OAAO,CAAC,CAAC,CAAC;AACX,CAAC"}

package/dist/gs/math/log.gs.d.ts

@@ -0,0 +1,2 @@
+export declare function Log(x: number): number;
+export declare function log(x: number): number;

package/dist/gs/math/log.gs.js

@@ -0,0 +1,16 @@
+// archLog import removed - using optimized implementation
+// Log returns the natural logarithm of x.
+//
+// Special cases are:
+//
+//	Log(+Inf) = +Inf
+//	Log(0) = -Inf
+//	Log(x < 0) = NaN
+//	Log(NaN) = NaN
+export function Log(x) {
+    return Math.log(x);
+}
+export function log(x) {
+    return Math.log(x);
+}
+//# sourceMappingURL=log.gs.js.map

package/dist/gs/math/log.gs.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"log.gs.js","sourceRoot":"","sources":["../../../gs/math/log.gs.ts"],"names":[],"mappings":"AAGA,0DAA0D;AAE1D,0CAA0C;AAC1C,EAAE;AACF,qBAAqB;AACrB,EAAE;AACF,mBAAmB;AACnB,gBAAgB;AAChB,mBAAmB;AACnB,iBAAiB;AACjB,MAAM,UAAU,GAAG,CAAC,CAAS;IAC5B,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAA;AACnB,CAAC;AAED,MAAM,UAAU,GAAG,CAAC,CAAS;IAC5B,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAA;AACnB,CAAC"}

package/dist/gs/math/log10.gs.d.ts

@@ -0,0 +1,4 @@
+export declare function Log10(x: number): number;
+export declare function log10(x: number): number;
+export declare function Log2(x: number): number;
+export declare function log2(x: number): number;

package/dist/gs/math/log10.gs.js

@@ -0,0 +1,17 @@
+// Log10 returns the decimal logarithm of x.
+// The special cases are the same as for [Log].
+export function Log10(x) {
+    return Math.log10(x);
+}
+export function log10(x) {
+    return Math.log10(x);
+}
+// Log2 returns the binary logarithm of x.
+// The special cases are the same as for [Log].
+export function Log2(x) {
+    return Math.log2(x);
+}
+export function log2(x) {
+    return Math.log2(x);
+}
+//# sourceMappingURL=log10.gs.js.map

package/dist/gs/math/log10.gs.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"log10.gs.js","sourceRoot":"","sources":["../../../gs/math/log10.gs.ts"],"names":[],"mappings":"AAKA,4CAA4C;AAC5C,+CAA+C;AAC/C,MAAM,UAAU,KAAK,CAAC,CAAS;IAC9B,OAAO,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAA;AACrB,CAAC;AAED,MAAM,UAAU,KAAK,CAAC,CAAS;IAC9B,OAAO,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAA;AACrB,CAAC;AAED,0CAA0C;AAC1C,+CAA+C;AAC/C,MAAM,UAAU,IAAI,CAAC,CAAS;IAC7B,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAA;AACpB,CAAC;AAED,MAAM,UAAU,IAAI,CAAC,CAAS;IAC7B,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAA;AACpB,CAAC"}

package/dist/gs/math/log1p.gs.d.ts

@@ -0,0 +1,2 @@
+export declare function Log1p(x: number): number;
+export declare function log1p(x: number): number;

package/dist/gs/math/log1p.gs.js

@@ -0,0 +1,17 @@
+// Log1p returns the natural logarithm of 1 plus its argument x.
+// It is more accurate than [Log](1 + x) when x is near zero.
+//
+// Special cases are:
+//
+//	Log1p(+Inf) = +Inf
+//	Log1p(±0) = ±0
+//	Log1p(-1) = -Inf
+//	Log1p(x < -1) = NaN
+//	Log1p(NaN) = NaN
+export function Log1p(x) {
+    return Math.log1p(x);
+}
+export function log1p(x) {
+    return Math.log1p(x);
+}
+//# sourceMappingURL=log1p.gs.js.map

package/dist/gs/math/log1p.gs.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"log1p.gs.js","sourceRoot":"","sources":["../../../gs/math/log1p.gs.ts"],"names":[],"mappings":"AAMA,gEAAgE;AAChE,6DAA6D;AAC7D,EAAE;AACF,qBAAqB;AACrB,EAAE;AACF,qBAAqB;AACrB,iBAAiB;AACjB,mBAAmB;AACnB,sBAAsB;AACtB,mBAAmB;AACnB,MAAM,UAAU,KAAK,CAAC,CAAS;IAC9B,OAAO,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAA;AACrB,CAAC;AAED,MAAM,UAAU,KAAK,CAAC,CAAS;IAC9B,OAAO,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAA;AACrB,CAAC"}

package/dist/gs/math/logb.gs.d.ts

@@ -0,0 +1,3 @@
+export declare function Logb(x: number): number;
+export declare function Ilogb(x: number): number;
+export declare function ilogb(x: number): number;

package/dist/gs/math/logb.gs.js

@@ -0,0 +1,43 @@
+// Logb returns the binary exponent of x.
+//
+// Special cases are:
+//
+//	Logb(±Inf) = +Inf
+//	Logb(0) = -Inf
+//	Logb(NaN) = NaN
+export function Logb(x) {
+    // Handle special cases
+    if (x === 0) {
+        return Number.NEGATIVE_INFINITY;
+    }
+    if (!Number.isFinite(x)) {
+        return Number.POSITIVE_INFINITY;
+    }
+    if (Number.isNaN(x)) {
+        return x;
+    }
+    return Math.floor(Math.log2(Math.abs(x)));
+}
+// Ilogb returns the binary exponent of x as an integer.
+//
+// Special cases are:
+//
+//	Ilogb(±Inf) = MaxInt32
+//	Ilogb(0) = MinInt32
+//	Ilogb(NaN) = MaxInt32
+export function Ilogb(x) {
+    // Handle special cases
+    if (x === 0) {
+        return -2147483648; // MinInt32
+    }
+    if (Number.isNaN(x) || !Number.isFinite(x)) {
+        return 2147483647; // MaxInt32
+    }
+    return Math.floor(Math.log2(Math.abs(x)));
+}
+// ilogb returns the binary exponent of x. It assumes x is finite and
+// non-zero.
+export function ilogb(x) {
+    return Math.floor(Math.log2(Math.abs(x)));
+}
+//# sourceMappingURL=logb.gs.js.map

package/dist/gs/math/logb.gs.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"logb.gs.js","sourceRoot":"","sources":["../../../gs/math/logb.gs.ts"],"names":[],"mappings":"AAIA,yCAAyC;AACzC,EAAE;AACF,qBAAqB;AACrB,EAAE;AACF,oBAAoB;AACpB,iBAAiB;AACjB,kBAAkB;AAClB,MAAM,UAAU,IAAI,CAAC,CAAS;IAC7B,uBAAuB;IACvB,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QACb,OAAO,MAAM,CAAC,iBAAiB,CAAA;IAChC,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;QACzB,OAAO,MAAM,CAAC,iBAAiB,CAAA;IAChC,CAAC;IAED,IAAI,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC;QACrB,OAAO,CAAC,CAAA;IACT,CAAC;IAED,OAAO,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;AAC1C,CAAC;AAED,wDAAwD;AACxD,EAAE;AACF,qBAAqB;AACrB,EAAE;AACF,yBAAyB;AACzB,sBAAsB;AACtB,wBAAwB;AACxB,MAAM,UAAU,KAAK,CAAC,CAAS;IAC9B,uBAAuB;IACvB,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QACb,OAAO,CAAC,UAAU,CAAA,CAAC,WAAW;IAC/B,CAAC;IAED,IAAI,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;QAC5C,OAAO,UAAU,CAAA,CAAC,WAAW;IAC9B,CAAC;IAED,OAAO,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;AAC1C,CAAC;AAED,qEAAqE;AACrE,YAAY;AACZ,MAAM,UAAU,KAAK,CAAC,CAAS;IAC9B,OAAO,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;AAC1C,CAAC"}

package/dist/gs/math/mod.gs.d.ts

@@ -0,0 +1,2 @@
+export declare function Mod(x: number, y: number): number;
+export declare function mod(x: number, y: number): number;

package/dist/gs/math/mod.gs.js

@@ -0,0 +1,26 @@
+// Mod returns the floating-point remainder of x/y.
+// The magnitude of the result is less than y and its
+// sign agrees with that of x.
+//
+// Special cases are:
+//
+//	Mod(±Inf, y) = NaN
+//	Mod(NaN, y) = NaN
+//	Mod(x, 0) = NaN
+//	Mod(x, ±Inf) = x
+//	Mod(x, NaN) = NaN
+export function Mod(x, y) {
+    return mod(x, y);
+}
+export function mod(x, y) {
+    // Handle special cases
+    if (y === 0 || !Number.isFinite(x) || Number.isNaN(x) || Number.isNaN(y)) {
+        return Number.NaN;
+    }
+    if (!Number.isFinite(y)) {
+        return x;
+    }
+    // JavaScript's % operator already implements the correct behavior for Mod
+    return x % y;
+}
+//# sourceMappingURL=mod.gs.js.map

package/dist/gs/math/mod.gs.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"mod.gs.js","sourceRoot":"","sources":["../../../gs/math/mod.gs.ts"],"names":[],"mappings":"AAOA,mDAAmD;AACnD,qDAAqD;AACrD,8BAA8B;AAC9B,EAAE;AACF,qBAAqB;AACrB,EAAE;AACF,qBAAqB;AACrB,oBAAoB;AACpB,kBAAkB;AAClB,mBAAmB;AACnB,oBAAoB;AACpB,MAAM,UAAU,GAAG,CAAC,CAAS,EAAE,CAAS;IACvC,OAAO,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAA;AACjB,CAAC;AAED,MAAM,UAAU,GAAG,CAAC,CAAS,EAAE,CAAS;IACvC,uBAAuB;IACvB,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC;QAC1E,OAAO,MAAM,CAAC,GAAG,CAAA;IAClB,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;QACzB,OAAO,CAAC,CAAA;IACT,CAAC;IAED,0EAA0E;IAC1E,OAAO,CAAC,GAAG,CAAC,CAAA;AACb,CAAC"}

package/dist/gs/math/modf.gs.d.ts

@@ -0,0 +1,2 @@
+export declare function Modf(f: number): [number, number];
+export declare function modf(f: number): [number, number];

package/dist/gs/math/modf.gs.js

@@ -0,0 +1,24 @@
+// Modf returns integer and fractional floating-point numbers
+// that sum to f. Both values have the same sign as f.
+//
+// Special cases are:
+//
+//	Modf(±Inf) = ±Inf, NaN
+//	Modf(NaN) = NaN, NaN
+export function Modf(f) {
+    return modf(f);
+}
+export function modf(f) {
+    // Handle special cases
+    if (Number.isNaN(f)) {
+        return [Number.NaN, Number.NaN];
+    }
+    if (!Number.isFinite(f)) {
+        return [f, Number.NaN];
+    }
+    // Get integer part using Math.trunc (preserves sign)
+    const intPart = Math.trunc(f);
+    const fracPart = f - intPart;
+    return [intPart, fracPart];
+}
+//# sourceMappingURL=modf.gs.js.map

package/dist/gs/math/modf.gs.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"modf.gs.js","sourceRoot":"","sources":["../../../gs/math/modf.gs.ts"],"names":[],"mappings":"AAIA,6DAA6D;AAC7D,sDAAsD;AACtD,EAAE;AACF,qBAAqB;AACrB,EAAE;AACF,yBAAyB;AACzB,uBAAuB;AACvB,MAAM,UAAU,IAAI,CAAC,CAAS;IAC7B,OAAO,IAAI,CAAC,CAAC,CAAC,CAAA;AACf,CAAC;AAED,MAAM,UAAU,IAAI,CAAC,CAAS;IAC7B,uBAAuB;IACvB,IAAI,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC;QACrB,OAAO,CAAC,MAAM,CAAC,GAAG,EAAE,MAAM,CAAC,GAAG,CAAC,CAAA;IAChC,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;QACzB,OAAO,CAAC,CAAC,EAAE,MAAM,CAAC,GAAG,CAAC,CAAA;IACvB,CAAC;IAED,qDAAqD;IACrD,MAAM,OAAO,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAA;IAC7B,MAAM,QAAQ,GAAG,CAAC,GAAG,OAAO,CAAA;IAE5B,OAAO,CAAC,OAAO,EAAE,QAAQ,CAAC,CAAA;AAC3B,CAAC"}

package/dist/gs/math/nextafter.gs.d.ts

@@ -0,0 +1,2 @@
+export declare function Nextafter32(x: number, y: number): number;
+export declare function Nextafter(x: number, y: number): number;

package/dist/gs/math/nextafter.gs.js

@@ -0,0 +1,66 @@
+import { IsNaN, NaN } from "./bits.gs.js";
+import { Copysign } from "./copysign.gs.js";
+import { Float32bits, Float32frombits, Float64bits, Float64frombits } from "./unsafe.gs.js";
+// Nextafter32 returns the next representable float32 value after x towards y.
+//
+// Special cases are:
+//
+//	Nextafter32(x, x)   = x
+//	Nextafter32(NaN, y) = NaN
+//	Nextafter32(x, NaN) = NaN
+export function Nextafter32(x, y) {
+    let r = 0;
+    {
+        // special case
+        switch (true) {
+            case IsNaN(x) || IsNaN(y):
+                r = NaN();
+                break;
+            case x == y:
+                r = x;
+                break;
+            case x == 0:
+                r = Copysign(Float32frombits(1), y);
+                break;
+            case (y > x) == (x > 0):
+                r = Float32frombits(Float32bits(x) + 1);
+                break;
+            default:
+                r = Float32frombits(Float32bits(x) - 1);
+                break;
+        }
+        return r;
+    }
+}
+// Nextafter returns the next representable float64 value after x towards y.
+//
+// Special cases are:
+//
+//	Nextafter(x, x)   = x
+//	Nextafter(NaN, y) = NaN
+//	Nextafter(x, NaN) = NaN
+export function Nextafter(x, y) {
+    let r = 0;
+    {
+        // special case
+        switch (true) {
+            case IsNaN(x) || IsNaN(y):
+                r = NaN();
+                break;
+            case x == y:
+                r = x;
+                break;
+            case x == 0:
+                r = Copysign(Float64frombits(1), y);
+                break;
+            case (y > x) == (x > 0):
+                r = Float64frombits(Float64bits(x) + 1);
+                break;
+            default:
+                r = Float64frombits(Float64bits(x) - 1);
+                break;
+        }
+        return r;
+    }
+}
+//# sourceMappingURL=nextafter.gs.js.map

package/dist/gs/math/nextafter.gs.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"nextafter.gs.js","sourceRoot":"","sources":["../../../gs/math/nextafter.gs.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,KAAK,EAAE,GAAG,EAAE,MAAM,cAAc,CAAC;AAC1C,OAAO,EAAE,QAAQ,EAAE,MAAM,kBAAkB,CAAC;AAC5C,OAAO,EAAE,WAAW,EAAE,eAAe,EAAE,WAAW,EAAE,eAAe,EAAE,MAAM,gBAAgB,CAAC;AAE5F,8EAA8E;AAC9E,EAAE;AACF,qBAAqB;AACrB,EAAE;AACF,0BAA0B;AAC1B,4BAA4B;AAC5B,4BAA4B;AAC5B,MAAM,UAAU,WAAW,CAAC,CAAS,EAAE,CAAS;IAC/C,IAAI,CAAC,GAAW,CAAC,CAAA;IACjB,CAAC;QAEA,eAAe;QACf,QAAQ,IAAI,EAAE,CAAC;YACd,KAAK,KAAK,CAAE,CAAY,CAAC,IAAI,KAAK,CAAE,CAAY,CAAC;gBAChD,CAAC,GAAI,GAAG,EAAa,CAAA;gBACrB,MAAK;YACN,KAAK,CAAC,IAAI,CAAC;gBACV,CAAC,GAAG,CAAC,CAAA;gBACL,MAAK;YACN,KAAK,CAAC,IAAI,CAAC;gBACV,CAAC,GAAI,QAAQ,CAAE,eAAe,CAAC,CAAC,CAAY,EAAG,CAAY,CAAY,CAAA;gBACvE,MAAK;YACN,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC;gBACtB,CAAC,GAAG,eAAe,CAAC,WAAW,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;gBACvC,MAAK;YACN;gBACC,CAAC,GAAG,eAAe,CAAC,WAAW,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;gBACvC,MAAK;QACP,CAAC;QACD,OAAO,CAAC,CAAA;IACT,CAAC;AACF,CAAC;AAED,4EAA4E;AAC5E,EAAE;AACF,qBAAqB;AACrB,EAAE;AACF,wBAAwB;AACxB,0BAA0B;AAC1B,0BAA0B;AAC1B,MAAM,UAAU,SAAS,CAAC,CAAS,EAAE,CAAS;IAC7C,IAAI,CAAC,GAAW,CAAC,CAAA;IACjB,CAAC;QAEA,eAAe;QACf,QAAQ,IAAI,EAAE,CAAC;YACd,KAAK,KAAK,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC,CAAC,CAAC;gBACxB,CAAC,GAAG,GAAG,EAAE,CAAA;gBACT,MAAK;YACN,KAAK,CAAC,IAAI,CAAC;gBACV,CAAC,GAAG,CAAC,CAAA;gBACL,MAAK;YACN,KAAK,CAAC,IAAI,CAAC;gBACV,CAAC,GAAG,QAAQ,CAAC,eAAe,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAA;gBACnC,MAAK;YACN,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC;gBACtB,CAAC,GAAG,eAAe,CAAC,WAAW,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;gBACvC,MAAK;YACN;gBACC,CAAC,GAAG,eAAe,CAAC,WAAW,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;gBACvC,MAAK;QACP,CAAC;QACD,OAAO,CAAC,CAAA;IACT,CAAC;AACF,CAAC"}

package/dist/gs/math/pow.gs.d.ts

@@ -0,0 +1,3 @@
+export declare function isOddInt(x: number): boolean;
+export declare function Pow(x: number, y: number): number;
+export declare function pow(x: number, y: number): number;

package/dist/gs/math/pow.gs.js

@@ -0,0 +1,40 @@
+export function isOddInt(x) {
+    // 1 << 53 is the largest exact integer in the float64 format.
+    if (Math.abs(x) >= (1 << 53)) {
+        return false;
+    }
+    const truncated = Math.trunc(x);
+    const fractional = x - truncated;
+    return fractional === 0 && (truncated & 1) === 1;
+}
+// Pow returns x**y, the base-x exponential of y.
+//
+// Special cases are (in order):
+//
+//	Pow(x, ±0) = 1 for any x
+//	Pow(1, y) = 1 for any y
+//	Pow(x, 1) = x for any x
+//	Pow(NaN, y) = NaN
+//	Pow(x, NaN) = NaN
+//	Pow(±0, y) = ±Inf for y an odd integer < 0
+//	Pow(±0, -Inf) = +Inf
+//	Pow(±0, +Inf) = +0
+//	Pow(±0, y) = +Inf for finite y < 0 and not an odd integer
+//	Pow(±0, y) = ±0 for y an odd integer > 0
+//	Pow(±0, y) = +0 for finite y > 0 and not an odd integer
+//	Pow(-1, ±Inf) = 1
+//	Pow(x, +Inf) = +Inf for |x| > 1
+//	Pow(x, -Inf) = +0 for |x| > 1
+//	Pow(x, +Inf) = +0 for |x| < 1
+//	Pow(x, -Inf) = +Inf for |x| < 1
+//	Pow(+Inf, y) = +Inf for y > 0
+//	Pow(+Inf, y) = +0 for y < 0
+//	Pow(-Inf, y) = Pow(-0, -y)
+//	Pow(x, y) = NaN for finite x < 0 and finite non-integer y
+export function Pow(x, y) {
+    return Math.pow(x, y);
+}
+export function pow(x, y) {
+    return Math.pow(x, y);
+}
+//# sourceMappingURL=pow.gs.js.map

package/dist/gs/math/pow.gs.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"pow.gs.js","sourceRoot":"","sources":["../../../gs/math/pow.gs.ts"],"names":[],"mappings":"AAYA,MAAM,UAAU,QAAQ,CAAC,CAAS;IACjC,8DAA8D;IAC9D,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC;QAC9B,OAAO,KAAK,CAAA;IACb,CAAC;IAED,MAAM,SAAS,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAA;IAC/B,MAAM,UAAU,GAAG,CAAC,GAAG,SAAS,CAAA;IAChC,OAAO,UAAU,KAAK,CAAC,IAAI,CAAC,SAAS,GAAG,CAAC,CAAC,KAAK,CAAC,CAAA;AACjD,CAAC;AAED,iDAAiD;AACjD,EAAE;AACF,gCAAgC;AAChC,EAAE;AACF,2BAA2B;AAC3B,0BAA0B;AAC1B,0BAA0B;AAC1B,oBAAoB;AACpB,oBAAoB;AACpB,6CAA6C;AAC7C,uBAAuB;AACvB,qBAAqB;AACrB,4DAA4D;AAC5D,2CAA2C;AAC3C,0DAA0D;AAC1D,oBAAoB;AACpB,kCAAkC;AAClC,gCAAgC;AAChC,gCAAgC;AAChC,kCAAkC;AAClC,gCAAgC;AAChC,8BAA8B;AAC9B,6BAA6B;AAC7B,4DAA4D;AAC5D,MAAM,UAAU,GAAG,CAAC,CAAS,EAAE,CAAS;IACvC,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAA;AACtB,CAAC;AAED,MAAM,UAAU,GAAG,CAAC,CAAS,EAAE,CAAS;IACvC,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAA;AACtB,CAAC"}

package/dist/gs/math/pow10.gs.d.ts

@@ -0,0 +1 @@
+export declare function Pow10(n: number): number;

package/dist/gs/math/pow10.gs.js

@@ -0,0 +1,14 @@
+import * as $ from "@goscript/builtin/builtin.js";
+let pow10tab = $.arrayToSlice([1e00, 1e01, 1e02, 1e03, 1e04, 1e05, 1e06, 1e07, 1e08, 1e09, 1e10, 1e11, 1e12, 1e13, 1e14, 1e15, 1e16, 1e17, 1e18, 1e19, 1e20, 1e21, 1e22, 1e23, 1e24, 1e25, 1e26, 1e27, 1e28, 1e29, 1e30, 1e31]);
+let pow10postab32 = $.arrayToSlice([1e00, 1e32, 1e64, 1e96, 1e128, 1e160, 1e192, 1e224, 1e256, 1e288]);
+let pow10negtab32 = $.arrayToSlice([1e-00, 1e-32, 1e-64, 1e-96, 1e-128, 1e-160, 1e-192, 1e-224, 1e-256, 1e-288, 1e-320]);
+// Pow10 returns 10**n, the base-10 exponential of n.
+//
+// Special cases are:
+//
+//	Pow10(n) =    0 for n < -323
+//	Pow10(n) = +Inf for n > 308
+export function Pow10(n) {
+    return Math.pow(10, n);
+}
+//# sourceMappingURL=pow10.gs.js.map

package/dist/gs/math/pow10.gs.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"pow10.gs.js","sourceRoot":"","sources":["../../../gs/math/pow10.gs.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,CAAC,MAAM,8BAA8B,CAAC;AAGlD,IAAI,QAAQ,GAAG,CAAC,CAAC,YAAY,CAAS,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC,CAAA;AAEvO,IAAI,aAAa,GAAG,CAAC,CAAC,YAAY,CAAS,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,CAAC,CAAC,CAAA;AAE9G,IAAI,aAAa,GAAG,CAAC,CAAC,YAAY,CAAS,CAAC,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC,CAAA;AAEhI,qDAAqD;AACrD,EAAE;AACF,qBAAqB;AACrB,EAAE;AACF,+BAA+B;AAC/B,8BAA8B;AAC9B,MAAM,UAAU,KAAK,CAAC,CAAS;IAC9B,OAAO,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,CAAC,CAAA;AACvB,CAAC"}

package/dist/gs/math/remainder.gs.d.ts

@@ -0,0 +1,2 @@
+export declare function Remainder(x: number, y: number): number;
+export declare function remainder(x: number, y: number): number;

package/dist/gs/math/remainder.gs.js

@@ -0,0 +1,25 @@
+// Remainder returns the IEEE 754 floating-point remainder of x/y.
+//
+// Special cases are:
+//
+//	Remainder(±Inf, y) = NaN
+//	Remainder(NaN, y) = NaN
+//	Remainder(x, 0) = NaN
+//	Remainder(x, ±Inf) = x
+//	Remainder(x, NaN) = NaN
+export function Remainder(x, y) {
+    return remainder(x, y);
+}
+export function remainder(x, y) {
+    // Handle special cases
+    if (Number.isNaN(x) || Number.isNaN(y) || !Number.isFinite(x) || y === 0) {
+        return Number.NaN;
+    }
+    if (!Number.isFinite(y)) {
+        return x;
+    }
+    // IEEE 754 remainder: x - n*y where n is the integer nearest to x/y
+    const n = Math.round(x / y);
+    return x - n * y;
+}
+//# sourceMappingURL=remainder.gs.js.map

package/dist/gs/math/remainder.gs.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"remainder.gs.js","sourceRoot":"","sources":["../../../gs/math/remainder.gs.ts"],"names":[],"mappings":"AAKA,kEAAkE;AAClE,EAAE;AACF,qBAAqB;AACrB,EAAE;AACF,2BAA2B;AAC3B,0BAA0B;AAC1B,wBAAwB;AACxB,yBAAyB;AACzB,0BAA0B;AAC1B,MAAM,UAAU,SAAS,CAAC,CAAS,EAAE,CAAS;IAC7C,OAAO,SAAS,CAAC,CAAC,EAAE,CAAC,CAAC,CAAA;AACvB,CAAC;AAED,MAAM,UAAU,SAAS,CAAC,CAAS,EAAE,CAAS;IAC7C,uBAAuB;IACvB,IAAI,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QAC1E,OAAO,MAAM,CAAC,GAAG,CAAA;IAClB,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;QACzB,OAAO,CAAC,CAAA;IACT,CAAC;IAED,oEAAoE;IACpE,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;IAC3B,OAAO,CAAC,GAAG,CAAC,GAAG,CAAC,CAAA;AACjB,CAAC"}