goscript 0.0.34 → 0.0.36

package/dist/gs/math/const.gs.js

@@ -0,0 +1,61 @@
+// https://oeis.org/A001113
+export let E = Math.E;
+// https://oeis.org/A000796
+export let Pi = Math.PI;
+// https://oeis.org/A001622
+export let Phi = 1.61803398874989484820458683436563811772030917980576286213544862;
+// https://oeis.org/A002193
+export let Sqrt2 = Math.SQRT2;
+// https://oeis.org/A019774
+export let SqrtE = Math.sqrt(Math.E);
+// https://oeis.org/A002161
+export let SqrtPi = Math.sqrt(Math.PI);
+// https://oeis.org/A139339
+export let SqrtPhi = Math.sqrt(1.61803398874989484820458683436563811772030917980576286213544862);
+// https://oeis.org/A002162
+export let Ln2 = Math.LN2;
+export let Log2E = Math.LOG2E;
+// https://oeis.org/A002392
+export let Ln10 = Math.LN10;
+export let Log10E = Math.LOG10E;
+// 3.40282346638528859811704183484516925440e+38
+export let MaxFloat32 = 3.4028234663852886e+38;
+// 1.401298464324817070923729583289916131280e-45
+export let SmallestNonzeroFloat32 = 1.401298464324817e-45;
+// 1.79769313486231570814527423731704356798070e+308
+export let MaxFloat64 = Number.MAX_VALUE;
+// 4.9406564584124654417656879286822137236505980e-324
+export let SmallestNonzeroFloat64 = Number.MIN_VALUE;
+// 32 or 64
+let intSize = (32 << ((~0 >> 63)));
+// MaxInt32 or MaxInt64 depending on intSize.
+export let MaxInt = 9223372036854775807n;
+// MinInt32 or MinInt64 depending on intSize.
+export let MinInt = -9223372036854775808n;
+// 127
+export let MaxInt8 = 127;
+// -128
+export let MinInt8 = -128;
+// 32767
+export let MaxInt16 = 32767;
+// -32768
+export let MinInt16 = -32768;
+// 2147483647
+export let MaxInt32 = 2147483647;
+// -2147483648
+export let MinInt32 = -2147483648;
+// 9223372036854775807
+export let MaxInt64 = 9223372036854775807n;
+// -9223372036854775808
+export let MinInt64 = -9223372036854775808n;
+// MaxUint32 or MaxUint64 depending on intSize.
+export let MaxUint = 0xffffffffffffffffn;
+// 255
+export let MaxUint8 = 255;
+// 65535
+export let MaxUint16 = 65535;
+// 4294967295
+export let MaxUint32 = 4294967295;
+// 18446744073709551615
+export let MaxUint64 = 0xffffffffffffffffn;
+//# sourceMappingURL=const.gs.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"const.gs.js","sourceRoot":"","sources":["../../../gs/math/const.gs.ts"],"names":[],"mappings":"AAEA,2BAA2B;AAC3B,MAAM,CAAC,IAAI,CAAC,GAAW,IAAI,CAAC,CAAC,CAAA;AAE7B,2BAA2B;AAC3B,MAAM,CAAC,IAAI,EAAE,GAAW,IAAI,CAAC,EAAE,CAAA;AAE/B,2BAA2B;AAC3B,MAAM,CAAC,IAAI,GAAG,GAAW,gEAAgE,CAAA;AAEzF,2BAA2B;AAC3B,MAAM,CAAC,IAAI,KAAK,GAAW,IAAI,CAAC,KAAK,CAAA;AAErC,2BAA2B;AAC3B,MAAM,CAAC,IAAI,KAAK,GAAW,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAA;AAE5C,2BAA2B;AAC3B,MAAM,CAAC,IAAI,MAAM,GAAW,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAA;AAE9C,2BAA2B;AAC3B,MAAM,CAAC,IAAI,OAAO,GAAW,IAAI,CAAC,IAAI,CAAC,gEAAgE,CAAC,CAAA;AAExG,2BAA2B;AAC3B,MAAM,CAAC,IAAI,GAAG,GAAW,IAAI,CAAC,GAAG,CAAA;AAEjC,MAAM,CAAC,IAAI,KAAK,GAAW,IAAI,CAAC,KAAK,CAAA;AAErC,2BAA2B;AAC3B,MAAM,CAAC,IAAI,IAAI,GAAW,IAAI,CAAC,IAAI,CAAA;AAEnC,MAAM,CAAC,IAAI,MAAM,GAAW,IAAI,CAAC,MAAM,CAAA;AAEvC,+CAA+C;AAC/C,MAAM,CAAC,IAAI,UAAU,GAAW,sBAAsB,CAAA;AAEtD,gDAAgD;AAChD,MAAM,CAAC,IAAI,sBAAsB,GAAW,qBAAqB,CAAA;AAEjE,mDAAmD;AACnD,MAAM,CAAC,IAAI,UAAU,GAAW,MAAM,CAAC,SAAS,CAAA;AAEhD,qDAAqD;AACrD,MAAM,CAAC,IAAI,sBAAsB,GAAW,MAAM,CAAC,SAAS,CAAA;AAE5D,WAAW;AACX,IAAI,OAAO,GAAW,CAAC,EAAE,IAAI,CAAC,CAAC,CAAE,CAAY,IAAI,EAAE,CAAC,CAAC,CAAC,CAAA;AAEtD,6CAA6C;AAC7C,MAAM,CAAC,IAAI,MAAM,GAAW,oBAAoB,CAAA;AAEhD,6CAA6C;AAC7C,MAAM,CAAC,IAAI,MAAM,GAAW,CAAC,oBAAoB,CAAA;AAEjD,MAAM;AACN,MAAM,CAAC,IAAI,OAAO,GAAW,GAAG,CAAA;AAEhC,OAAO;AACP,MAAM,CAAC,IAAI,OAAO,GAAW,CAAC,GAAG,CAAA;AAEjC,QAAQ;AACR,MAAM,CAAC,IAAI,QAAQ,GAAW,KAAK,CAAA;AAEnC,SAAS;AACT,MAAM,CAAC,IAAI,QAAQ,GAAW,CAAC,KAAK,CAAA;AAEpC,aAAa;AACb,MAAM,CAAC,IAAI,QAAQ,GAAW,UAAU,CAAA;AAExC,cAAc;AACd,MAAM,CAAC,IAAI,QAAQ,GAAW,CAAC,UAAU,CAAA;AAEzC,sBAAsB;AACtB,MAAM,CAAC,IAAI,QAAQ,GAAW,oBAAoB,CAAA;AAElD,uBAAuB;AACvB,MAAM,CAAC,IAAI,QAAQ,GAAW,CAAC,oBAAoB,CAAA;AAEnD,+CAA+C;AAC/C,MAAM,CAAC,IAAI,OAAO,GAAW,mBAAmB,CAAA;AAEhD,MAAM;AACN,MAAM,CAAC,IAAI,QAAQ,GAAW,GAAG,CAAA;AAEjC,QAAQ;AACR,MAAM,CAAC,IAAI,SAAS,GAAW,KAAK,CAAA;AAEpC,aAAa;AACb,MAAM,CAAC,IAAI,SAAS,GAAW,UAAU,CAAA;AAEzC,uBAAuB;AACvB,MAAM,CAAC,IAAI,SAAS,GAAW,mBAAmB,CAAA"}

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

@@ -0,0 +1 @@
+export declare function Copysign(f: number, sign: number): number;

package/dist/gs/math/copysign.gs.js

@@ -0,0 +1,20 @@
+// Copysign returns a value with the magnitude of f
+// and the sign of sign.
+export function Copysign(f, sign) {
+    // Handle special cases for zero
+    if (f === 0) {
+        return sign < 0 || Object.is(sign, -0) ? -0 : 0;
+    }
+    const magnitude = Math.abs(f);
+    const signValue = Math.sign(sign);
+    // Handle NaN case
+    if (Number.isNaN(sign)) {
+        return signValue < 0 ? -magnitude : magnitude;
+    }
+    // Handle negative zero case
+    if (Object.is(sign, -0)) {
+        return -magnitude;
+    }
+    return signValue < 0 ? -magnitude : magnitude;
+}
+//# sourceMappingURL=copysign.gs.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"copysign.gs.js","sourceRoot":"","sources":["../../../gs/math/copysign.gs.ts"],"names":[],"mappings":"AAGA,mDAAmD;AACnD,wBAAwB;AACxB,MAAM,UAAU,QAAQ,CAAC,CAAS,EAAE,IAAY;IAC/C,gCAAgC;IAChC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QACb,OAAO,IAAI,GAAG,CAAC,IAAI,MAAM,CAAC,EAAE,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;IAChD,CAAC;IAED,MAAM,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAA;IAC7B,MAAM,SAAS,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAA;IAEjC,kBAAkB;IAClB,IAAI,MAAM,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC;QACxB,OAAO,SAAS,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,SAAS,CAAA;IAC9C,CAAC;IAED,4BAA4B;IAC5B,IAAI,MAAM,CAAC,EAAE,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;QACzB,OAAO,CAAC,SAAS,CAAA;IAClB,CAAC;IAED,OAAO,SAAS,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,SAAS,CAAA;AAC9C,CAAC"}

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

@@ -0,0 +1,5 @@
+export declare function Dim(x: number, y: number): number;
+export declare function Max(x: number, y: number): number;
+export declare function max(x: number, y: number): number;
+export declare function Min(x: number, y: number): number;
+export declare function min(x: number, y: number): number;

package/dist/gs/math/dim.gs.js

@@ -0,0 +1,69 @@
+// Dim returns the maximum of x-y or 0.
+//
+// Special cases are:
+//
+//	Dim(+Inf, +Inf) = NaN
+//	Dim(-Inf, -Inf) = NaN
+//	Dim(x, NaN) = Dim(NaN, x) = NaN
+export function Dim(x, y) {
+    const v = x - y;
+    return v <= 0 ? 0 : v;
+}
+// Max returns the larger of x or y.
+//
+// Special cases are:
+//
+//	Max(x, +Inf) = Max(+Inf, x) = +Inf
+//	Max(x, NaN) = Max(NaN, x) = NaN
+//	Max(+0, ±0) = Max(±0, +0) = +0
+//	Max(-0, -0) = -0
+//
+// Note that this differs from the built-in function max when called
+// with NaN and +Inf.
+export function Max(x, y) {
+    return max(x, y);
+}
+export function max(x, y) {
+    // Handle NaN cases
+    if (Number.isNaN(x) || Number.isNaN(y)) {
+        return Number.NaN;
+    }
+    // Handle infinity cases
+    if (x === Number.POSITIVE_INFINITY || y === Number.POSITIVE_INFINITY) {
+        return Number.POSITIVE_INFINITY;
+    }
+    // Handle zero cases - prefer +0 over -0
+    if (x === 0 && y === 0) {
+        return Object.is(x, -0) ? y : x;
+    }
+    return Math.max(x, y);
+}
+// Min returns the smaller of x or y.
+//
+// Special cases are:
+//
+//	Min(x, -Inf) = Min(-Inf, x) = -Inf
+//	Min(x, NaN) = Min(NaN, x) = NaN
+//	Min(-0, ±0) = Min(±0, -0) = -0
+//
+// Note that this differs from the built-in function min when called
+// with NaN and -Inf.
+export function Min(x, y) {
+    return min(x, y);
+}
+export function min(x, y) {
+    // Handle NaN cases
+    if (Number.isNaN(x) || Number.isNaN(y)) {
+        return Number.NaN;
+    }
+    // Handle infinity cases
+    if (x === Number.NEGATIVE_INFINITY || y === Number.NEGATIVE_INFINITY) {
+        return Number.NEGATIVE_INFINITY;
+    }
+    // Handle zero cases - prefer -0 over +0
+    if (x === 0 && y === 0) {
+        return Object.is(x, -0) ? x : y;
+    }
+    return Math.min(x, y);
+}
+//# sourceMappingURL=dim.gs.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"dim.gs.js","sourceRoot":"","sources":["../../../gs/math/dim.gs.ts"],"names":[],"mappings":"AAKA,uCAAuC;AACvC,EAAE;AACF,qBAAqB;AACrB,EAAE;AACF,wBAAwB;AACxB,wBAAwB;AACxB,kCAAkC;AAClC,MAAM,UAAU,GAAG,CAAC,CAAS,EAAE,CAAS;IACvC,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAA;IACf,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;AACtB,CAAC;AAED,oCAAoC;AACpC,EAAE;AACF,qBAAqB;AACrB,EAAE;AACF,qCAAqC;AACrC,kCAAkC;AAClC,iCAAiC;AACjC,mBAAmB;AACnB,EAAE;AACF,oEAAoE;AACpE,qBAAqB;AACrB,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,mBAAmB;IACnB,IAAI,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC;QACxC,OAAO,MAAM,CAAC,GAAG,CAAA;IAClB,CAAC;IAED,wBAAwB;IACxB,IAAI,CAAC,KAAK,MAAM,CAAC,iBAAiB,IAAI,CAAC,KAAK,MAAM,CAAC,iBAAiB,EAAE,CAAC;QACtE,OAAO,MAAM,CAAC,iBAAiB,CAAA;IAChC,CAAC;IAED,wCAAwC;IACxC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QACxB,OAAO,MAAM,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;IAChC,CAAC;IAED,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAA;AACtB,CAAC;AAED,qCAAqC;AACrC,EAAE;AACF,qBAAqB;AACrB,EAAE;AACF,qCAAqC;AACrC,kCAAkC;AAClC,iCAAiC;AACjC,EAAE;AACF,oEAAoE;AACpE,qBAAqB;AACrB,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,mBAAmB;IACnB,IAAI,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC;QACxC,OAAO,MAAM,CAAC,GAAG,CAAA;IAClB,CAAC;IAED,wBAAwB;IACxB,IAAI,CAAC,KAAK,MAAM,CAAC,iBAAiB,IAAI,CAAC,KAAK,MAAM,CAAC,iBAAiB,EAAE,CAAC;QACtE,OAAO,MAAM,CAAC,iBAAiB,CAAA;IAChC,CAAC;IAED,wCAAwC;IACxC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QACxB,OAAO,MAAM,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;IAChC,CAAC;IAED,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAA;AACtB,CAAC"}

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

@@ -0,0 +1,4 @@
+export declare function Erf(x: number): number;
+export declare function erf(x: number): number;
+export declare function Erfc(x: number): number;
+export declare function erfc(x: number): number;

package/dist/gs/math/erf.gs.js

@@ -0,0 +1,336 @@
+import { IsInf, IsNaN, NaN } from "./bits.gs.js";
+import { Exp } from "./exp.gs.js";
+import { Float64bits, Float64frombits } from "./unsafe.gs.js";
+// 0x3FEB0AC160000000
+let erx = 8.45062911510467529297e-01;
+// Coefficients for approximation to  erf in [0, 0.84375]
+// 0x3FC06EBA8214DB69
+let efx = 1.28379167095512586316e-01;
+// 0x3FF06EBA8214DB69
+let efx8 = 1.02703333676410069053e+00;
+// 0x3FC06EBA8214DB68
+let pp0 = 1.28379167095512558561e-01;
+// 0xBFD4CD7D691CB913
+let pp1 = -3.25042107247001499370e-01;
+// 0xBF9D2A51DBD7194F
+let pp2 = -2.84817495755985104766e-02;
+// 0xBF77A291236668E4
+let pp3 = -5.77027029648944159157e-03;
+// 0xBEF8EAD6120016AC
+let pp4 = -2.37630166566501626084e-05;
+// 0x3FD97779CDDADC09
+let qq1 = 3.97917223959155352819e-01;
+// 0x3FB0A54C5536CEBA
+let qq2 = 6.50222499887672944485e-02;
+// 0x3F74D022C4D36B0F
+let qq3 = 5.08130628187576562776e-03;
+// 0x3F215DC9221C1A10
+let qq4 = 1.32494738004321644526e-04;
+// 0xBED09C4342A26120
+let qq5 = -3.96022827877536812320e-06;
+// Coefficients for approximation to  erf  in [0.84375, 1.25]
+// 0xBF6359B8BEF77538
+let pa0 = -2.36211856075265944077e-03;
+// 0x3FDA8D00AD92B34D
+let pa1 = 4.14856118683748331666e-01;
+// 0xBFD7D240FBB8C3F1
+let pa2 = -3.72207876035701323847e-01;
+// 0x3FD45FCA805120E4
+let pa3 = 3.18346619901161753674e-01;
+// 0xBFBC63983D3E28EC
+let pa4 = -1.10894694282396677476e-01;
+// 0x3FA22A36599795EB
+let pa5 = 3.54783043256182359371e-02;
+// 0xBF61BF380A96073F
+let pa6 = -2.16637559486879084300e-03;
+// 0x3FBB3E6618EEE323
+let qa1 = 1.06420880400844228286e-01;
+// 0x3FE14AF092EB6F33
+let qa2 = 5.40397917702171048937e-01;
+// 0x3FB2635CD99FE9A7
+let qa3 = 7.18286544141962662868e-02;
+// 0x3FC02660E763351F
+let qa4 = 1.26171219808761642112e-01;
+// 0x3F8BEDC26B51DD1C
+let qa5 = 1.36370839120290507362e-02;
+// 0x3F888B545735151D
+let qa6 = 1.19844998467991074170e-02;
+// Coefficients for approximation to  erfc in [1.25, 1/0.35]
+// 0xBF843412600D6435
+let ra0 = -9.86494403484714822705e-03;
+// 0xBFE63416E4BA7360
+let ra1 = -6.93858572707181764372e-01;
+// 0xC0251E0441B0E726
+let ra2 = -1.05586262253232909814e+01;
+// 0xC04F300AE4CBA38D
+let ra3 = -6.23753324503260060396e+01;
+// 0xC0644CB184282266
+let ra4 = -1.62396669462573470355e+02;
+// 0xC067135CEBCCABB2
+let ra5 = -1.84605092906711035994e+02;
+// 0xC054526557E4D2F2
+let ra6 = -8.12874355063065934246e+01;
+// 0xC023A0EFC69AC25C
+let ra7 = -9.81432934416914548592e+00;
+// 0x4033A6B9BD707687
+let sa1 = 1.96512716674392571292e+01;
+// 0x4061350C526AE721
+let sa2 = 1.37657754143519042600e+02;
+// 0x407B290DD58A1A71
+let sa3 = 4.34565877475229228821e+02;
+// 0x40842B1921EC2868
+let sa4 = 6.45387271733267880336e+02;
+// 0x407AD02157700314
+let sa5 = 4.29008140027567833386e+02;
+// 0x405B28A3EE48AE2C
+let sa6 = 1.08635005541779435134e+02;
+// 0x401A47EF8E484A93
+let sa7 = 6.57024977031928170135e+00;
+// 0xBFAEEFF2EE749A62
+let sa8 = -6.04244152148580987438e-02;
+// Coefficients for approximation to  erfc in [1/.35, 28]
+// 0xBF84341239E86F4A
+let rb0 = -9.86494292470009928597e-03;
+// 0xBFE993BA70C285DE
+let rb1 = -7.99283237680523006574e-01;
+// 0xC031C209555F995A
+let rb2 = -1.77579549177547519889e+01;
+// 0xC064145D43C5ED98
+let rb3 = -1.60636384855821916062e+02;
+// 0xC083EC881375F228
+let rb4 = -6.37566443368389627722e+02;
+// 0xC09004616A2E5992
+let rb5 = -1.02509513161107724954e+03;
+// 0xC07E384E9BDC383F
+let rb6 = -4.83519191608651397019e+02;
+// 0x403E568B261D5190
+let sb1 = 3.03380607434824582924e+01;
+// 0x40745CAE221B9F0A
+let sb2 = 3.25792512996573918826e+02;
+// 0x409802EB189D5118
+let sb3 = 1.53672958608443695994e+03;
+// 0x40A8FFB7688C246A
+let sb4 = 3.19985821950859553908e+03;
+// 0x40A3F219CEDF3BE6
+let sb5 = 2.55305040643316442583e+03;
+// 0x407DA874E79FE763
+let sb6 = 4.74528541206955367215e+02;
+// 0xC03670E242712D62
+let sb7 = -2.24409524465858183362e+01;
+// Erf returns the error function of x.
+//
+// Special cases are:
+//
+//	Erf(+Inf) = 1
+//	Erf(-Inf) = -1
+//	Erf(NaN) = NaN
+export function Erf(x) {
+    return erf(x);
+}
+export function erf(x) {
+    // 0x0080000000000000
+    // 2**-28
+    // 0x0080000000000000
+    let VeryTiny = 2.848094538889218e-306;
+    // 2**-28
+    let Small = 1.0 / ((1 << 28));
+    // special cases
+    switch (true) {
+        case IsNaN(x):
+            return NaN();
+            break;
+        case IsInf(x, 1):
+            return 1;
+            break;
+        case IsInf(x, -1):
+            return -1;
+            break;
+    }
+    let sign = false;
+    if (x < 0) {
+        x = -x;
+        sign = true;
+    }
+    // |x| < 0.84375
+    // |x| < 2**-28
+    // avoid underflow
+    if (x < 0.84375) {
+        // |x| < 0.84375
+        let temp = 0;
+        // |x| < 2**-28
+        // avoid underflow
+        if (x < 3.72529e-09) {
+            // |x| < 2**-28
+            // avoid underflow
+            if (x < 2.84809e-306) {
+                temp = 0.125 * (8.0 * x + 1.02703 * x); // avoid underflow
+            }
+            else {
+                temp = x + 0.128379 * x;
+            }
+        }
+        else {
+            let z = x * x;
+            let r = 0.128379 + z * (-0.325042 + z * (-0.0284817 + z * (-0.00577027 + z * -2.3763e-05)));
+            let s = 1 + z * (0.397917 + z * (0.0650222 + z * (0.00508131 + z * (0.000132495 + z * -3.96023e-06))));
+            let y = r / s;
+            temp = x + x * y;
+        }
+        if (sign) {
+            return -temp;
+        }
+        return temp;
+    }
+    // 0.84375 <= |x| < 1.25
+    if (x < 1.25) {
+        // 0.84375 <= |x| < 1.25
+        let s = x - 1;
+        let P = -0.00236212 + s * (0.414856 + s * (-0.372208 + s * (0.318347 + s * (-0.110895 + s * (0.0354783 + s * -0.00216638)))));
+        let Q = 1 + s * (0.106421 + s * (0.540398 + s * (0.0718287 + s * (0.126171 + s * (0.0136371 + s * 0.0119845)))));
+        if (sign) {
+            return -0.845063 - P / Q;
+        }
+        return 0.845063 + P / Q;
+    }
+    // inf > |x| >= 6
+    if (x >= 6) {
+        // inf > |x| >= 6
+        if (sign) {
+            return -1;
+        }
+        return 1;
+    }
+    let s = 1 / (x * x);
+    let R;
+    let S;
+    // |x| < 1 / 0.35  ~ 2.857143
+    // |x| >= 1 / 0.35  ~ 2.857143
+    if (x < 1 / 0.35) {
+        // |x| < 1 / 0.35  ~ 2.857143
+        R = -0.00986494 + s * (-0.693859 + s * (-10.5586 + s * (-62.3753 + s * (-162.397 + s * (-184.605 + s * (-81.2874 + s * -9.81433))))));
+        S = 1 + s * (19.6513 + s * (137.658 + s * (434.566 + s * (645.387 + s * (429.008 + s * (108.635 + s * (6.57025 + s * -0.0604244)))))));
+    }
+    else {
+        // |x| >= 1 / 0.35  ~ 2.857143
+        R = -0.00986494 + s * (-0.799283 + s * (-17.758 + s * (-160.636 + s * (-637.566 + s * (-1025.1 + s * -483.519)))));
+        S = 1 + s * (30.3381 + s * (325.793 + s * (1536.73 + s * (3199.86 + s * (2553.05 + s * (474.529 + s * -22.441))))));
+    }
+    let z = Float64frombits((Float64bits(x) & 0xffffffff00000000)); // pseudo-single (20-bit) precision x
+    let r = Exp(-z * z - 0.5625) * Exp((z - x) * (z + x) + R / S);
+    if (sign) {
+        return r / x - 1;
+    }
+    return 1 - r / x;
+}
+// Erfc returns the complementary error function of x.
+//
+// Special cases are:
+//
+//	Erfc(+Inf) = 0
+//	Erfc(-Inf) = 2
+//	Erfc(NaN) = NaN
+export function Erfc(x) {
+    return erfc(x);
+}
+export function erfc(x) {
+    // 2**-56
+    let Tiny = 1.0 / ((1 << 56));
+    // special cases
+    switch (true) {
+        case IsNaN(x):
+            return NaN();
+            break;
+        case IsInf(x, 1):
+            return 0;
+            break;
+        case IsInf(x, -1):
+            return 2;
+            break;
+    }
+    let sign = false;
+    if (x < 0) {
+        x = -x;
+        sign = true;
+    }
+    // |x| < 0.84375
+    // |x| < 2**-56
+    // |x| < 1/4
+    if (x < 0.84375) {
+        // |x| < 0.84375
+        let temp = 0;
+        // |x| < 2**-56
+        // |x| < 1/4
+        if (x < 1.38778e-17) {
+            // |x| < 2**-56
+            temp = x;
+        }
+        else {
+            let z = x * x;
+            let r = 0.128379 + z * (-0.325042 + z * (-0.0284817 + z * (-0.00577027 + z * -2.3763e-05)));
+            let s = 1 + z * (0.397917 + z * (0.0650222 + z * (0.00508131 + z * (0.000132495 + z * -3.96023e-06))));
+            let y = r / s;
+            // |x| < 1/4
+            if (x < 0.25) {
+                // |x| < 1/4
+                temp = x + x * y;
+            }
+            else {
+                temp = 0.5 + (x * y + (x - 0.5));
+            }
+        }
+        if (sign) {
+            return 1 + temp;
+        }
+        return 1 - temp;
+    }
+    // 0.84375 <= |x| < 1.25
+    if (x < 1.25) {
+        // 0.84375 <= |x| < 1.25
+        let s = x - 1;
+        let P = -0.00236212 + s * (0.414856 + s * (-0.372208 + s * (0.318347 + s * (-0.110895 + s * (0.0354783 + s * -0.00216638)))));
+        let Q = 1 + s * (0.106421 + s * (0.540398 + s * (0.0718287 + s * (0.126171 + s * (0.0136371 + s * 0.0119845)))));
+        if (sign) {
+            return 1 + 0.845063 + P / Q;
+        }
+        return 1 - 0.845063 - P / Q;
+    }
+    // |x| < 28
+    // |x| < 1 / 0.35 ~ 2.857143
+    // |x| >= 1 / 0.35 ~ 2.857143
+    // x < -6
+    // pseudo-single (20-bit) precision x
+    if (x < 28) {
+        // |x| < 28
+        let s = 1 / (x * x);
+        let R;
+        let S;
+        // |x| < 1 / 0.35 ~ 2.857143
+        // |x| >= 1 / 0.35 ~ 2.857143
+        // x < -6
+        if (x < 1 / 0.35) {
+            // |x| < 1 / 0.35 ~ 2.857143
+            R = -0.00986494 + s * (-0.693859 + s * (-10.5586 + s * (-62.3753 + s * (-162.397 + s * (-184.605 + s * (-81.2874 + s * -9.81433))))));
+            S = 1 + s * (19.6513 + s * (137.658 + s * (434.566 + s * (645.387 + s * (429.008 + s * (108.635 + s * (6.57025 + s * -0.0604244)))))));
+        }
+        else {
+            // |x| >= 1 / 0.35 ~ 2.857143
+            // x < -6
+            if (sign && x > 6) {
+                return 2;
+            }
+            R = -0.00986494 + s * (-0.799283 + s * (-17.758 + s * (-160.636 + s * (-637.566 + s * (-1025.1 + s * -483.519)))));
+            S = 1 + s * (30.3381 + s * (325.793 + s * (1536.73 + s * (3199.86 + s * (2553.05 + s * (474.529 + s * -22.441))))));
+        }
+        let z = Float64frombits((Float64bits(x) & 0xffffffff00000000)); // pseudo-single (20-bit) precision x
+        let r = Exp(-z * z - 0.5625) * Exp((z - x) * (z + x) + R / S);
+        if (sign) {
+            return 2 - r / x;
+        }
+        return r / x;
+    }
+    if (sign) {
+        return 2;
+    }
+    return 0;
+}
+//# sourceMappingURL=erf.gs.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"erf.gs.js","sourceRoot":"","sources":["../../../gs/math/erf.gs.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,KAAK,EAAE,KAAK,EAAE,GAAG,EAAE,MAAM,cAAc,CAAC;AACjD,OAAO,EAAE,GAAG,EAAE,MAAM,aAAa,CAAC;AAElC,OAAO,EAAE,WAAW,EAAE,eAAe,EAAE,MAAM,gBAAgB,CAAC;AAE9D,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,yDAAyD;AACzD,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,IAAI,GAAW,0BAA0B,CAAA;AAE7C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,6DAA6D;AAC7D,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,4DAA4D;AAC5D,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,yDAAyD;AACzD,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,0BAA0B,CAAA;AAE5C,qBAAqB;AACrB,IAAI,GAAG,GAAW,CAAC,0BAA0B,CAAA;AAE7C,uCAAuC;AACvC,EAAE;AACF,qBAAqB;AACrB,EAAE;AACF,gBAAgB;AAChB,iBAAiB;AACjB,iBAAiB;AACjB,MAAM,UAAU,GAAG,CAAC,CAAS;IAC5B,OAAO,GAAG,CAAC,CAAC,CAAC,CAAA;AACd,CAAC;AAED,MAAM,UAAU,GAAG,CAAC,CAAS;IAE5B,qBAAqB;IACrB,SAAS;IACT,qBAAqB;IACrB,IAAI,QAAQ,GAAW,sBAAsB,CAAA;IAC7C,SAAS;IACT,IAAI,KAAK,GAAW,GAAG,GAAG,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAA;IACrC,gBAAgB;IAChB,QAAQ,IAAI,EAAE,CAAC;QACd,KAAK,KAAK,CAAC,CAAC,CAAC;YACZ,OAAO,GAAG,EAAE,CAAA;YACZ,MAAK;QACN,KAAK,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC;YACf,OAAO,CAAC,CAAA;YACR,MAAK;QACN,KAAK,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YAChB,OAAO,CAAC,CAAC,CAAA;YACT,MAAK;IACP,CAAC;IACD,IAAI,IAAI,GAAG,KAAK,CAAA;IAChB,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QACX,CAAC,GAAG,CAAC,CAAC,CAAA;QACN,IAAI,GAAG,IAAI,CAAA;IACZ,CAAC;IACD,gBAAgB;IAEhB,eAAe;IAEf,kBAAkB;IAClB,IAAI,CAAC,GAAG,OAAO,EAAE,CAAC;QACjB,gBAAgB;QAChB,IAAI,IAAI,GAAW,CAAC,CAAA;QACpB,eAAe;QAEf,kBAAkB;QAClB,IAAI,CAAC,GAAG,WAAW,EAAE,CAAC;YACrB,eAAe;YAEf,kBAAkB;YAClB,IAAI,CAAC,GAAG,YAAY,EAAE,CAAC;gBACtB,IAAI,GAAG,KAAK,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,OAAO,GAAG,CAAC,CAAC,CAAA,CAAC,kBAAkB;YAC1D,CAAC;iBAAM,CAAC;gBACP,IAAI,GAAG,CAAC,GAAG,QAAQ,GAAG,CAAC,CAAA;YACxB,CAAC;QACF,CAAC;aAAM,CAAC;YACP,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAA;YACb,IAAI,CAAC,GAAG,QAAQ,GAAG,CAAC,GAAG,CAAC,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,CAAC,SAAS,GAAG,CAAC,GAAG,CAAC,CAAC,UAAU,GAAG,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,CAAA;YAC3F,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,SAAS,GAAG,CAAC,GAAG,CAAC,UAAU,GAAG,CAAC,GAAG,CAAC,WAAW,GAAG,CAAC,GAAG,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC,CAAA;YACtG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAA;YACb,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAA;QACjB,CAAC;QACD,IAAI,IAAI,EAAE,CAAC;YACV,OAAO,CAAC,IAAI,CAAA;QACb,CAAC;QACD,OAAO,IAAI,CAAA;IACZ,CAAC;IACD,wBAAwB;IACxB,IAAI,CAAC,GAAG,IAAI,EAAE,CAAC;QACd,wBAAwB;QACxB,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAA;QACb,IAAI,CAAC,GAAG,CAAC,UAAU,GAAG,CAAC,GAAG,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,SAAS,GAAG,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;QAC7H,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,SAAS,GAAG,CAAC,GAAG,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,SAAS,GAAG,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;QAChH,IAAI,IAAI,EAAE,CAAC;YACV,OAAO,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,CAAA;QACzB,CAAC;QACD,OAAO,QAAQ,GAAG,CAAC,GAAG,CAAC,CAAA;IACxB,CAAC;IACD,iBAAiB;IACjB,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QACZ,iBAAiB;QACjB,IAAI,IAAI,EAAE,CAAC;YACV,OAAO,CAAC,CAAC,CAAA;QACV,CAAC;QACD,OAAO,CAAC,CAAA;IACT,CAAC;IACD,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;IACnB,IAAI,CAAS,CAAA;IACb,IAAI,CAAS,CAAA;IACb,6BAA6B;IAE7B,8BAA8B;IAC9B,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,EAAE,CAAC;QAClB,6BAA6B;QAC7B,CAAC,GAAG,CAAC,UAAU,GAAG,CAAC,GAAG,CAAC,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;QACrI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;IACvI,CAAC;SAAM,CAAC;QACP,8BAA8B;QAC9B,CAAC,GAAG,CAAC,UAAU,GAAG,CAAC,GAAG,CAAC,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,CAAC,MAAM,GAAG,CAAC,GAAG,CAAC,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,CAAC,MAAM,GAAG,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;QAClH,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;IACpH,CAAC;IACD,IAAI,CAAC,GAAG,eAAe,CAAC,CAAC,WAAW,CAAC,CAAC,CAAC,GAAG,kBAAkB,CAAC,CAAC,CAAA,CAAC,qCAAqC;IACpG,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAA;IAC7D,IAAI,IAAI,EAAE,CAAC;QACV,OAAO,CAAC,GAAG,CAAC,GAAG,CAAC,CAAA;IACjB,CAAC;IACD,OAAO,CAAC,GAAG,CAAC,GAAG,CAAC,CAAA;AACjB,CAAC;AAED,sDAAsD;AACtD,EAAE;AACF,qBAAqB;AACrB,EAAE;AACF,iBAAiB;AACjB,iBAAiB;AACjB,kBAAkB;AAClB,MAAM,UAAU,IAAI,CAAC,CAAS;IAC7B,OAAO,IAAI,CAAC,CAAC,CAAC,CAAA;AACf,CAAC;AAED,MAAM,UAAU,IAAI,CAAC,CAAS;IAC7B,SAAS;IACT,IAAI,IAAI,GAAW,GAAG,GAAG,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAA;IACpC,gBAAgB;IAChB,QAAQ,IAAI,EAAE,CAAC;QACd,KAAK,KAAK,CAAC,CAAC,CAAC;YACZ,OAAO,GAAG,EAAE,CAAA;YACZ,MAAK;QACN,KAAK,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC;YACf,OAAO,CAAC,CAAA;YACR,MAAK;QACN,KAAK,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YAChB,OAAO,CAAC,CAAA;YACR,MAAK;IACP,CAAC;IACD,IAAI,IAAI,GAAG,KAAK,CAAA;IAChB,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QACX,CAAC,GAAG,CAAC,CAAC,CAAA;QACN,IAAI,GAAG,IAAI,CAAA;IACZ,CAAC;IACD,gBAAgB;IAEhB,eAAe;IAEf,YAAY;IACZ,IAAI,CAAC,GAAG,OAAO,EAAE,CAAC;QACjB,gBAAgB;QAChB,IAAI,IAAI,GAAW,CAAC,CAAA;QACpB,eAAe;QAEf,YAAY;QACZ,IAAI,CAAC,GAAG,WAAW,EAAE,CAAC;YACrB,eAAe;YACf,IAAI,GAAG,CAAC,CAAA;QACT,CAAC;aAAM,CAAC;YACP,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAA;YACb,IAAI,CAAC,GAAG,QAAQ,GAAG,CAAC,GAAG,CAAC,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,CAAC,SAAS,GAAG,CAAC,GAAG,CAAC,CAAC,UAAU,GAAG,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,CAAA;YAC3F,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,SAAS,GAAG,CAAC,GAAG,CAAC,UAAU,GAAG,CAAC,GAAG,CAAC,WAAW,GAAG,CAAC,GAAG,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC,CAAA;YACtG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAA;YACb,YAAY;YACZ,IAAI,CAAC,GAAG,IAAI,EAAE,CAAC;gBACd,YAAY;gBACZ,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAA;YACjB,CAAC;iBAAM,CAAC;gBACP,IAAI,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAA;YACjC,CAAC;QACF,CAAC;QACD,IAAI,IAAI,EAAE,CAAC;YACV,OAAO,CAAC,GAAG,IAAI,CAAA;QAChB,CAAC;QACD,OAAO,CAAC,GAAG,IAAI,CAAA;IAChB,CAAC;IACD,wBAAwB;IACxB,IAAI,CAAC,GAAG,IAAI,EAAE,CAAC;QACd,wBAAwB;QACxB,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAA;QACb,IAAI,CAAC,GAAG,CAAC,UAAU,GAAG,CAAC,GAAG,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,SAAS,GAAG,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;QAC7H,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,SAAS,GAAG,CAAC,GAAG,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,SAAS,GAAG,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;QAChH,IAAI,IAAI,EAAE,CAAC;YACV,OAAO,CAAC,GAAG,QAAQ,GAAG,CAAC,GAAG,CAAC,CAAA;QAC5B,CAAC;QACD,OAAO,CAAC,GAAG,QAAQ,GAAG,CAAC,GAAG,CAAC,CAAA;IAE5B,CAAC;IACD,WAAW;IAEX,4BAA4B;IAE5B,6BAA6B;IAE7B,SAAS;IAET,qCAAqC;IACrC,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC;QACZ,WAAW;QACX,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;QACnB,IAAI,CAAS,CAAA;QACb,IAAI,CAAS,CAAA;QACb,4BAA4B;QAE5B,6BAA6B;QAE7B,SAAS;QACT,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,EAAE,CAAC;YAClB,4BAA4B;YAC5B,CAAC,GAAG,CAAC,UAAU,GAAG,CAAC,GAAG,CAAC,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;YACrI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;QACvI,CAAC;aAAM,CAAC;YACP,6BAA6B;YAE7B,SAAS;YACT,IAAI,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;gBACnB,OAAO,CAAC,CAAA;YACT,CAAC;YACD,CAAC,GAAG,CAAC,UAAU,GAAG,CAAC,GAAG,CAAC,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,CAAC,MAAM,GAAG,CAAC,GAAG,CAAC,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,CAAC,MAAM,GAAG,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;YAClH,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;QACpH,CAAC;QACD,IAAI,CAAC,GAAG,eAAe,CAAC,CAAC,WAAW,CAAC,CAAC,CAAC,GAAG,kBAAkB,CAAC,CAAC,CAAA,CAAC,qCAAqC;QACpG,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAA;QAC7D,IAAI,IAAI,EAAE,CAAC;YACV,OAAO,CAAC,GAAG,CAAC,GAAG,CAAC,CAAA;QACjB,CAAC;QACD,OAAO,CAAC,GAAG,CAAC,CAAA;IACb,CAAC;IACD,IAAI,IAAI,EAAE,CAAC;QACV,OAAO,CAAC,CAAA;IACT,CAAC;IACD,OAAO,CAAC,CAAA;AACT,CAAC"}

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

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

package/dist/gs/math/erfinv.gs.js

@@ -0,0 +1,118 @@
+import * as $ from "@goscript/builtin/builtin.js";
+import { Inf, IsNaN, NaN } from "./bits.gs.js";
+import { Log } from "./log.gs.js";
+import { Sqrt } from "./sqrt.gs.js";
+// Coefficients for approximation to erf in |x| <= 0.85
+let a0 = 1.1975323115670912564578e0;
+let a1 = 4.7072688112383978012285e1;
+let a2 = 6.9706266534389598238465e2;
+let a3 = 4.8548868893843886794648e3;
+let a4 = 1.6235862515167575384252e4;
+let a5 = 2.3782041382114385731252e4;
+let a6 = 1.1819493347062294404278e4;
+let a7 = 8.8709406962545514830200e2;
+let b0 = 1.0000000000000000000e0;
+let b1 = 4.2313330701600911252e1;
+let b2 = 6.8718700749205790830e2;
+let b3 = 5.3941960214247511077e3;
+let b4 = 2.1213794301586595867e4;
+let b5 = 3.9307895800092710610e4;
+let b6 = 2.8729085735721942674e4;
+let b7 = 5.2264952788528545610e3;
+// Coefficients for approximation to erf in 0.85 < |x| <= 1-2*exp(-25)
+let c0 = 1.42343711074968357734e0;
+let c1 = 4.63033784615654529590e0;
+let c2 = 5.76949722146069140550e0;
+let c3 = 3.64784832476320460504e0;
+let c4 = 1.27045825245236838258e0;
+let c5 = 2.41780725177450611770e-1;
+let c6 = 2.27238449892691845833e-2;
+let c7 = 7.74545014278341407640e-4;
+let d0 = 1.4142135623730950488016887e0;
+let d1 = 2.9036514445419946173133295e0;
+let d2 = 2.3707661626024532365971225e0;
+let d3 = 9.7547832001787427186894837e-1;
+let d4 = 2.0945065210512749128288442e-1;
+let d5 = 2.1494160384252876777097297e-2;
+let d6 = 7.7441459065157709165577218e-4;
+let d7 = 1.4859850019840355905497876e-9;
+// Coefficients for approximation to erf in 1-2*exp(-25) < |x| < 1
+let e0 = 6.65790464350110377720e0;
+let e1 = 5.46378491116411436990e0;
+let e2 = 1.78482653991729133580e0;
+let e3 = 2.96560571828504891230e-1;
+let e4 = 2.65321895265761230930e-2;
+let e5 = 1.24266094738807843860e-3;
+let e6 = 2.71155556874348757815e-5;
+let e7 = 2.01033439929228813265e-7;
+let f0 = 1.414213562373095048801689e0;
+let f1 = 8.482908416595164588112026e-1;
+let f2 = 1.936480946950659106176712e-1;
+let f3 = 2.103693768272068968719679e-2;
+let f4 = 1.112800997078859844711555e-3;
+let f5 = 2.611088405080593625138020e-5;
+let f6 = 2.010321207683943062279931e-7;
+let f7 = 2.891024605872965461538222e-15;
+// Erfinv returns the inverse error function of x.
+//
+// Special cases are:
+//
+//	Erfinv(1) = +Inf
+//	Erfinv(-1) = -Inf
+//	Erfinv(x) = NaN if x < -1 or x > 1
+//	Erfinv(NaN) = NaN
+export function Erfinv(x) {
+    // special cases
+    if (IsNaN(x) || x <= -1 || x >= 1) {
+        if (x == -1 || x == 1) {
+            return Inf($.int(x));
+        }
+        return NaN();
+    }
+    let sign = false;
+    if (x < 0) {
+        x = -x;
+        sign = true;
+    }
+    let ans = 0;
+    // |x| <= 0.85
+    if (x <= 0.85) {
+        // |x| <= 0.85
+        let r = 0.180625 - 0.25 * x * x;
+        let z1 = ((((((887.094 * r + 11819.5) * r + 23782.041382114385) * r + 16235.9) * r + 4854.89) * r + 697.063) * r + 47.0727) * r + 1.19753;
+        let z2 = ((((((5226.5 * r + 28729.1) * r + 39307.9) * r + 21213.8) * r + 5394.2) * r + 687.187) * r + 42.3133) * r + 1;
+        ans = (x * z1) / z2;
+    }
+    else {
+        let z1;
+        let z2;
+        let r = Sqrt(0.693147 - Log(1.0 - x));
+        if (r <= 5.0) {
+            r -= 1.6;
+            z1 = ((((((0.000774545 * r + 0.0227238) * r + 0.241781) * r + 1.27046) * r + 3.64785) * r + 5.7695) * r + 4.63034) * r + 1.42344;
+            z2 = ((((((1.48599e-09 * r + 0.000774415) * r + 0.0214942) * r + 0.209451) * r + 0.975478) * r + 2.37077) * r + 2.90365) * r + 1.41421;
+        }
+        else {
+            r -= 5.0;
+            z1 = ((((((2.01033e-07 * r + 2.71156e-05) * r + 0.00124266) * r + 0.0265322) * r + 0.296561) * r + 1.78483) * r + 5.46378) * r + 6.6579;
+            z2 = ((((((2.89102e-15 * r + 2.01032e-07) * r + 2.61109e-05) * r + 0.0011128) * r + 0.0210369) * r + 0.193648) * r + 0.848291) * r + 1.41421;
+        }
+        ans = z1 / z2;
+    }
+    if (sign) {
+        return -ans;
+    }
+    return ans;
+}
+// Erfcinv returns the inverse of [Erfc](x).
+//
+// Special cases are:
+//
+//	Erfcinv(0) = +Inf
+//	Erfcinv(2) = -Inf
+//	Erfcinv(x) = NaN if x < 0 or x > 2
+//	Erfcinv(NaN) = NaN
+export function Erfcinv(x) {
+    return Erfinv(1 - x);
+}
+//# sourceMappingURL=erfinv.gs.js.map