dsp-collection 0.2.5 → 0.2.7

package/math/ComplexArray.js

@@ -1,171 +1,171 @@
-import Complex from "./Complex.js";
-import MutableComplex from "./MutableComplex.js";
-import { assert } from "../utils/MiscUtils.js";
-const emptyFloat64Array = new Float64Array(0);
-export default class ComplexArray {
-    constructor(x = 0) {
-        if (typeof x == "number") {
-            this.constructByLength(x);
-        }
-        else if (Array.isArray(x) && x[0] instanceof Complex) {
-            this.constructByArrayOfComplex(x);
-        }
-        else if (x instanceof Object && x.length !== undefined) {
-            this.constructByArrayOfNumber(x);
-        }
-        else {
-            throw new Error("Invalid constructor argument.");
-        }
-    }
-    constructByLength(length) {
-        this.length = length;
-        if (length) {
-            this.re = new Float64Array(length);
-            this.im = new Float64Array(length);
-        }
-        else {
-            this.re = emptyFloat64Array;
-            this.im = emptyFloat64Array;
-        }
-    }
-    constructByArrayOfComplex(a) {
-        this.length = a.length;
-        this.re = new Float64Array(a.length);
-        this.im = new Float64Array(a.length);
-        for (let i = 0; i < a.length; i++) {
-            this.re[i] = a[i].re;
-            this.im[i] = a[i].im;
-        }
-    }
-    constructByArrayOfNumber(a) {
-        this.length = a.length;
-        this.re = new Float64Array(a);
-        this.im = new Float64Array(a.length);
-    }
-    static fromPolar(absArray, argArray) {
-        const n = absArray.length;
-        assert(n == argArray.length);
-        const a = new ComplexArray(n);
-        for (let i = 0; i < n; i++) {
-            a.setPolar(i, absArray[i], argArray[i]);
-        }
-        return a;
-    }
-    slice(begin, end) {
-        const a2 = new ComplexArray();
-        a2.re = this.re.slice(begin, end);
-        a2.im = this.im.slice(begin, end);
-        a2.length = a2.re.length;
-        return a2;
-    }
-    subarray(begin, end) {
-        const a2 = new ComplexArray();
-        a2.re = this.re.subarray(begin, end);
-        a2.im = this.im.subarray(begin, end);
-        a2.length = end - begin;
-        return a2;
-    }
-    set(i, c) {
-        this.re[i] = c.re;
-        this.im[i] = c.im;
-    }
-    setReIm(i, re, im) {
-        this.re[i] = re;
-        this.im[i] = im;
-    }
-    setPolar(i, abs, arg) {
-        this.re[i] = abs * Math.cos(arg);
-        this.im[i] = abs * Math.sin(arg);
-    }
-    static copy1(a1, i1, a2, i2) {
-        a2.re[i2] = a1.re[i1];
-        a2.im[i2] = a1.im[i1];
-    }
-    get(i) {
-        return new MutableComplex(this.re[i], this.im[i]);
-    }
-    getAbs(i) {
-        return Math.hypot(this.re[i], this.im[i]);
-    }
-    getArg(i) {
-        return Math.atan2(this.im[i], this.re[i]);
-    }
-    toString() {
-        let s = "[";
-        for (let i = 0; i < this.length; i++) {
-            if (i > 0) {
-                s += ", ";
-            }
-            s += "(" + this.re[i] + ", " + this.im[i] + ")";
-        }
-        s += "]";
-        return s;
-    }
-    getAbsArray() {
-        const n = this.length;
-        const a = new Float64Array(n);
-        for (let i = 0; i < n; i++) {
-            a[i] = this.getAbs(i);
-        }
-        return a;
-    }
-    getArgArray() {
-        const n = this.length;
-        const a = new Float64Array(n);
-        for (let i = 0; i < n; i++) {
-            a[i] = this.getArg(i);
-        }
-        return a;
-    }
-    addRealTo(i, x) {
-        this.re[i] += x;
-    }
-    addTo(i, x) {
-        this.re[i] += x.re;
-        this.im[i] += x.im;
-    }
-    subRealFrom(i, x) {
-        this.re[i] -= x;
-    }
-    subFrom(i, x) {
-        this.re[i] -= x.re;
-        this.im[i] -= x.im;
-    }
-    mulByReal(i, x) {
-        this.re[i] *= x;
-        this.im[i] *= x;
-    }
-    mulBy(i, x) {
-        this.setMul(i, this.re[i], this.im[i], x.re, x.im);
-    }
-    divByReal(i, x) {
-        this.re[i] /= x;
-        this.im[i] /= x;
-    }
-    divBy(i, x) {
-        this.setDiv(i, this.re[i], this.im[i], x.re, x.im);
-    }
-    mulByArray(a2) {
-        const n = this.length;
-        assert(a2.length == n);
-        for (let i = 0; i < n; i++) {
-            this.setMul(i, this.re[i], this.im[i], a2.re[i], a2.im[i]);
-        }
-    }
-    mulAllByReal(x) {
-        const n = this.length;
-        for (let i = 0; i < n; i++) {
-            this.mulByReal(i, x);
-        }
-    }
-    setMul(i, re1, im1, re2, im2) {
-        this.re[i] = re1 * re2 - im1 * im2;
-        this.im[i] = re1 * im2 + im1 * re2;
-    }
-    setDiv(i, re1, im1, re2, im2) {
-        const m = re1 * re1 + im1 * im1;
-        this.re[i] = (re1 * re2 + im1 * im2) / m;
-        this.im[i] = (im1 * re2 - re1 * im2) / m;
-    }
-}
+import Complex from "./Complex.js";
+import MutableComplex from "./MutableComplex.js";
+import { assert } from "../utils/MiscUtils.js";
+const emptyFloat64Array = new Float64Array(0);
+export default class ComplexArray {
+    constructor(x = 0) {
+        if (typeof x == "number") {
+            this.constructByLength(x);
+        }
+        else if (Array.isArray(x) && x[0] instanceof Complex) {
+            this.constructByArrayOfComplex(x);
+        }
+        else if (x instanceof Object && x.length !== undefined) {
+            this.constructByArrayOfNumber(x);
+        }
+        else {
+            throw new Error("Invalid constructor argument.");
+        }
+    }
+    constructByLength(length) {
+        this.length = length;
+        if (length) {
+            this.re = new Float64Array(length);
+            this.im = new Float64Array(length);
+        }
+        else {
+            this.re = emptyFloat64Array;
+            this.im = emptyFloat64Array;
+        }
+    }
+    constructByArrayOfComplex(a) {
+        this.length = a.length;
+        this.re = new Float64Array(a.length);
+        this.im = new Float64Array(a.length);
+        for (let i = 0; i < a.length; i++) {
+            this.re[i] = a[i].re;
+            this.im[i] = a[i].im;
+        }
+    }
+    constructByArrayOfNumber(a) {
+        this.length = a.length;
+        this.re = new Float64Array(a);
+        this.im = new Float64Array(a.length);
+    }
+    static fromPolar(absArray, argArray) {
+        const n = absArray.length;
+        assert(n == argArray.length);
+        const a = new ComplexArray(n);
+        for (let i = 0; i < n; i++) {
+            a.setPolar(i, absArray[i], argArray[i]);
+        }
+        return a;
+    }
+    slice(begin, end) {
+        const a2 = new ComplexArray();
+        a2.re = this.re.slice(begin, end);
+        a2.im = this.im.slice(begin, end);
+        a2.length = a2.re.length;
+        return a2;
+    }
+    subarray(begin, end) {
+        const a2 = new ComplexArray();
+        a2.re = this.re.subarray(begin, end);
+        a2.im = this.im.subarray(begin, end);
+        a2.length = end - begin;
+        return a2;
+    }
+    set(i, c) {
+        this.re[i] = c.re;
+        this.im[i] = c.im;
+    }
+    setReIm(i, re, im) {
+        this.re[i] = re;
+        this.im[i] = im;
+    }
+    setPolar(i, abs, arg) {
+        this.re[i] = abs * Math.cos(arg);
+        this.im[i] = abs * Math.sin(arg);
+    }
+    static copy1(a1, i1, a2, i2) {
+        a2.re[i2] = a1.re[i1];
+        a2.im[i2] = a1.im[i1];
+    }
+    get(i) {
+        return new MutableComplex(this.re[i], this.im[i]);
+    }
+    getAbs(i) {
+        return Math.hypot(this.re[i], this.im[i]);
+    }
+    getArg(i) {
+        return Math.atan2(this.im[i], this.re[i]);
+    }
+    toString() {
+        let s = "[";
+        for (let i = 0; i < this.length; i++) {
+            if (i > 0) {
+                s += ", ";
+            }
+            s += "(" + this.re[i] + ", " + this.im[i] + ")";
+        }
+        s += "]";
+        return s;
+    }
+    getAbsArray() {
+        const n = this.length;
+        const a = new Float64Array(n);
+        for (let i = 0; i < n; i++) {
+            a[i] = this.getAbs(i);
+        }
+        return a;
+    }
+    getArgArray() {
+        const n = this.length;
+        const a = new Float64Array(n);
+        for (let i = 0; i < n; i++) {
+            a[i] = this.getArg(i);
+        }
+        return a;
+    }
+    addRealTo(i, x) {
+        this.re[i] += x;
+    }
+    addTo(i, x) {
+        this.re[i] += x.re;
+        this.im[i] += x.im;
+    }
+    subRealFrom(i, x) {
+        this.re[i] -= x;
+    }
+    subFrom(i, x) {
+        this.re[i] -= x.re;
+        this.im[i] -= x.im;
+    }
+    mulByReal(i, x) {
+        this.re[i] *= x;
+        this.im[i] *= x;
+    }
+    mulBy(i, x) {
+        this.setMul(i, this.re[i], this.im[i], x.re, x.im);
+    }
+    divByReal(i, x) {
+        this.re[i] /= x;
+        this.im[i] /= x;
+    }
+    divBy(i, x) {
+        this.setDiv(i, this.re[i], this.im[i], x.re, x.im);
+    }
+    mulByArray(a2) {
+        const n = this.length;
+        assert(a2.length == n);
+        for (let i = 0; i < n; i++) {
+            this.setMul(i, this.re[i], this.im[i], a2.re[i], a2.im[i]);
+        }
+    }
+    mulAllByReal(x) {
+        const n = this.length;
+        for (let i = 0; i < n; i++) {
+            this.mulByReal(i, x);
+        }
+    }
+    setMul(i, re1, im1, re2, im2) {
+        this.re[i] = re1 * re2 - im1 * im2;
+        this.im[i] = re1 * im2 + im1 * re2;
+    }
+    setDiv(i, re1, im1, re2, im2) {
+        const m = re1 * re1 + im1 * im1;
+        this.re[i] = (re1 * re2 + im1 * im2) / m;
+        this.im[i] = (im1 * re2 - re1 * im2) / m;
+    }
+}
 //# sourceMappingURL=ComplexArray.js.map

package/math/ComplexArray.js.map

@@ -1 +1 @@
-{"version":3,"file":"ComplexArray.js","sourceRoot":"","sources":["../../src/math/ComplexArray.ts"],"names":[],"mappings":"AAAA,OAAO,OAAO,MAAM,cAAc,CAAC;AACnC,OAAO,cAAc,MAAM,qBAAqB,CAAC;AACjD,OAAO,EAAC,MAAM,EAAC,MAAM,uBAAuB,CAAC;AAE7C,MAAM,iBAAiB,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;AAO9C,MAAM,CAAC,OAAO,OAAO,YAAY;IAM9B,YAAoB,IAA4C,CAAC;QAC9D,IAAI,OAAO,CAAC,IAAI,QAAQ,EAAE;YACvB,IAAI,CAAC,iBAAiB,CAAC,CAAC,CAAC,CAAC;SAAE;aACzB,IAAI,KAAK,CAAC,OAAO,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,YAAY,OAAO,EAAE;YACpD,IAAI,CAAC,yBAAyB,CAAC,CAAC,CAAC,CAAC;SAAE;aACjC,IAAI,CAAC,YAAY,MAAM,IAAI,CAAC,CAAC,MAAM,KAAK,SAAS,EAAE;YACtD,IAAI,CAAC,wBAAwB,CAAoB,CAAC,CAAC,CAAC;SAAE;aACnD;YACH,MAAM,IAAI,KAAK,CAAC,+BAA+B,CAAC,CAAC;SAAE;IAAA,CAAC;IAElD,iBAAiB,CAAE,MAAc;QACtC,IAAI,CAAC,MAAM,GAAG,MAAM,CAAC;QACrB,IAAI,MAAM,EAAE;YACT,IAAI,CAAC,EAAE,GAAG,IAAI,YAAY,CAAC,MAAM,CAAC,CAAC;YACnC,IAAI,CAAC,EAAE,GAAG,IAAI,YAAY,CAAC,MAAM,CAAC,CAAC;SAAE;aAClC;YACH,IAAI,CAAC,EAAE,GAAG,iBAAiB,CAAC;YAC5B,IAAI,CAAC,EAAE,GAAG,iBAAiB,CAAC;SAAE;IAAA,CAAC;IAE7B,yBAAyB,CAAE,CAAY;QAC5C,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,MAAM,CAAC;QACvB,IAAI,CAAC,EAAE,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;QACrC,IAAI,CAAC,EAAE,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;QACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;YAChC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;YACrB,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;SAAE;IAAA,CAAC;IAEtB,wBAAwB,CAAE,CAAoB;QACnD,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,MAAM,CAAC;QACvB,IAAI,CAAC,EAAE,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QAC9B,IAAI,CAAC,EAAE,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;IAAC,CAAC;IAEnC,MAAM,CAAC,SAAS,CAAE,QAA2B,EAAE,QAA2B;QAC9E,MAAM,CAAC,GAAG,QAAQ,CAAC,MAAM,CAAC;QAC1B,MAAM,CAAC,CAAC,IAAI,QAAQ,CAAC,MAAM,CAAC,CAAC;QAC7B,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YACzB,CAAC,CAAC,QAAQ,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;SAAE;QAC7C,OAAO,CAAC,CAAC;IAAC,CAAC;IAEP,KAAK,CAAE,KAAc,EAAE,GAAY;QACvC,MAAM,EAAE,GAAG,IAAI,YAAY,EAAE,CAAC;QAC9B,EAAE,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,KAAK,CAAC,KAAK,EAAE,GAAG,CAAC,CAAC;QAClC,EAAE,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,KAAK,CAAC,KAAK,EAAE,GAAG,CAAC,CAAC;QAClC,EAAE,CAAC,MAAM,GAAG,EAAE,CAAC,EAAE,CAAC,MAAM,CAAC;QACzB,OAAO,EAAE,CAAC;IAAC,CAAC;IAER,QAAQ,CAAE,KAAa,EAAE,GAAW;QACxC,MAAM,EAAE,GAAG,IAAI,YAAY,EAAE,CAAC;QAC9B,EAAE,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,QAAQ,CAAC,KAAK,EAAE,GAAG,CAAC,CAAC;QACrC,EAAE,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,QAAQ,CAAC,KAAK,EAAE,GAAG,CAAC,CAAC;QACrC,EAAE,CAAC,MAAM,GAAG,GAAG,GAAG,KAAK,CAAC;QACxB,OAAO,EAAE,CAAC;IAAC,CAAC;IAIR,GAAG,CAAE,CAAS,EAAE,CAAU;QAC9B,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;QAClB,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;IAAC,CAAC;IAEhB,OAAO,CAAE,CAAS,EAAE,EAAU,EAAE,EAAU;QAC9C,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QAChB,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;IAAC,CAAC;IAEd,QAAQ,CAAE,CAAS,EAAE,GAAW,EAAE,GAAW;QACjD,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACjC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAAC,CAAC;IAE/B,MAAM,CAAC,KAAK,CAAE,EAAgB,EAAE,EAAU,EAAE,EAAgB,EAAE,EAAU;QAC5E,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;QACtB,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IAAC,CAAC;IAIpB,GAAG,CAAE,CAAS;QAClB,OAAO,IAAI,cAAc,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAAC,CAAC;IAEhD,MAAM,CAAE,CAAS;QACrB,OAAO,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAAC,CAAC;IAExC,MAAM,CAAE,CAAS;QACrB,OAAO,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAAC,CAAC;IAIxC,QAAQ;QACZ,IAAI,CAAC,GAAG,GAAG,CAAC;QACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;YACnC,IAAI,CAAC,GAAG,CAAC,EAAE;gBACR,CAAC,IAAI,IAAI,CAAC;aAAE;YACf,CAAC,IAAI,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;SAAE;QACrD,CAAC,IAAI,GAAG,CAAC;QACT,OAAO,CAAC,CAAC;IAAC,CAAC;IAEP,WAAW;QACf,MAAM,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC;QACtB,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YACzB,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;SAAE;QAC3B,OAAO,CAAC,CAAC;IAAC,CAAC;IAEP,WAAW;QACf,MAAM,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC;QACtB,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YACzB,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;SAAE;QAC3B,OAAO,CAAC,CAAC;IAAC,CAAC;IAIP,SAAS,CAAE,CAAS,EAAE,CAAS;QACnC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;IAAC,CAAC;IAEd,KAAK,CAAE,CAAS,EAAE,CAAU;QAChC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC;QACnB,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC;IAAC,CAAC;IAEjB,WAAW,CAAE,CAAS,EAAE,CAAS;QACrC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;IAAC,CAAC;IAEd,OAAO,CAAE,CAAS,EAAE,CAAU;QAClC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC;QACnB,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC;IAAC,CAAC;IAEjB,SAAS,CAAE,CAAS,EAAE,CAAS;QACnC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;QAChB,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;IAAC,CAAC;IAEd,KAAK,CAAE,CAAS,EAAE,CAAU;QAChC,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;IAAC,CAAC;IAEjD,SAAS,CAAE,CAAS,EAAE,CAAS;QACnC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;QAChB,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;IAAC,CAAC;IAEd,KAAK,CAAE,CAAS,EAAE,CAAU;QAChC,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;IAAC,CAAC;IAMjD,UAAU,CAAE,EAAgB;QAChC,MAAM,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC;QACtB,MAAM,CAAC,EAAE,CAAC,MAAM,IAAI,CAAC,CAAC,CAAC;QACvB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YACzB,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;SAAE;IAAA,CAAC;IAE7D,YAAY,CAAE,CAAS;QAC3B,MAAM,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC;QACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YACzB,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;SAAE;IAAA,CAAC;IAOvB,MAAM,CAAE,CAAS,EAAE,GAAW,EAAE,GAAW,EAAE,GAAW,EAAE,GAAW;QACzE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;QACnC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;IAAC,CAAC;IAKjC,MAAM,CAAE,CAAS,EAAE,GAAW,EAAE,GAAW,EAAE,GAAW,EAAE,GAAW;QACzE,MAAM,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;QAChC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;QACzC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAAC,CAAC;CAE7C"}
+{"version":3,"file":"ComplexArray.js","sourceRoot":"","sources":["../../src/math/ComplexArray.ts"],"names":[],"mappings":"AAAA,OAAO,OAAO,MAAM,cAAc,CAAC;AACnC,OAAO,cAAc,MAAM,qBAAqB,CAAC;AACjD,OAAO,EAAC,MAAM,EAAC,MAAM,uBAAuB,CAAC;AAE7C,MAAM,iBAAiB,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;AAO9C,MAAM,CAAC,OAAO,OAAO,YAAY;IAM9B,YAAoB,IAA4C,CAAC;QAC9D,IAAI,OAAO,CAAC,IAAI,QAAQ,EAAE,CAAC;YACxB,IAAI,CAAC,iBAAiB,CAAC,CAAC,CAAC,CAAC;QAAC,CAAC;aACzB,IAAI,KAAK,CAAC,OAAO,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,YAAY,OAAO,EAAE,CAAC;YACrD,IAAI,CAAC,yBAAyB,CAAC,CAAC,CAAC,CAAC;QAAC,CAAC;aACjC,IAAI,CAAC,YAAY,MAAM,IAAI,CAAC,CAAC,MAAM,KAAK,SAAS,EAAE,CAAC;YACvD,IAAI,CAAC,wBAAwB,CAAoB,CAAC,CAAC,CAAC;QAAC,CAAC;aACnD,CAAC;YACJ,MAAM,IAAI,KAAK,CAAC,+BAA+B,CAAC,CAAC;QAAC,CAAC;IAAA,CAAC;IAElD,iBAAiB,CAAE,MAAc;QACtC,IAAI,CAAC,MAAM,GAAG,MAAM,CAAC;QACrB,IAAI,MAAM,EAAE,CAAC;YACV,IAAI,CAAC,EAAE,GAAG,IAAI,YAAY,CAAC,MAAM,CAAC,CAAC;YACnC,IAAI,CAAC,EAAE,GAAG,IAAI,YAAY,CAAC,MAAM,CAAC,CAAC;QAAC,CAAC;aAClC,CAAC;YACJ,IAAI,CAAC,EAAE,GAAG,iBAAiB,CAAC;YAC5B,IAAI,CAAC,EAAE,GAAG,iBAAiB,CAAC;QAAC,CAAC;IAAA,CAAC;IAE7B,yBAAyB,CAAE,CAAY;QAC5C,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,MAAM,CAAC;QACvB,IAAI,CAAC,EAAE,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;QACrC,IAAI,CAAC,EAAE,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;QACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YACjC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;YACrB,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;QAAC,CAAC;IAAA,CAAC;IAEtB,wBAAwB,CAAE,CAAoB;QACnD,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,MAAM,CAAC;QACvB,IAAI,CAAC,EAAE,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QAC9B,IAAI,CAAC,EAAE,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;IAAC,CAAC;IAEnC,MAAM,CAAC,SAAS,CAAE,QAA2B,EAAE,QAA2B;QAC9E,MAAM,CAAC,GAAG,QAAQ,CAAC,MAAM,CAAC;QAC1B,MAAM,CAAC,CAAC,IAAI,QAAQ,CAAC,MAAM,CAAC,CAAC;QAC7B,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC1B,CAAC,CAAC,QAAQ,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;QAAC,CAAC;QAC7C,OAAO,CAAC,CAAC;IAAC,CAAC;IAEP,KAAK,CAAE,KAAc,EAAE,GAAY;QACvC,MAAM,EAAE,GAAG,IAAI,YAAY,EAAE,CAAC;QAC9B,EAAE,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,KAAK,CAAC,KAAK,EAAE,GAAG,CAAC,CAAC;QAClC,EAAE,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,KAAK,CAAC,KAAK,EAAE,GAAG,CAAC,CAAC;QAClC,EAAE,CAAC,MAAM,GAAG,EAAE,CAAC,EAAE,CAAC,MAAM,CAAC;QACzB,OAAO,EAAE,CAAC;IAAC,CAAC;IAER,QAAQ,CAAE,KAAa,EAAE,GAAW;QACxC,MAAM,EAAE,GAAG,IAAI,YAAY,EAAE,CAAC;QAC9B,EAAE,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,QAAQ,CAAC,KAAK,EAAE,GAAG,CAAC,CAAC;QACrC,EAAE,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,QAAQ,CAAC,KAAK,EAAE,GAAG,CAAC,CAAC;QACrC,EAAE,CAAC,MAAM,GAAG,GAAG,GAAG,KAAK,CAAC;QACxB,OAAO,EAAE,CAAC;IAAC,CAAC;IAIR,GAAG,CAAE,CAAS,EAAE,CAAU;QAC9B,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;QAClB,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;IAAC,CAAC;IAEhB,OAAO,CAAE,CAAS,EAAE,EAAU,EAAE,EAAU;QAC9C,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QAChB,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;IAAC,CAAC;IAEd,QAAQ,CAAE,CAAS,EAAE,GAAW,EAAE,GAAW;QACjD,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACjC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAAC,CAAC;IAE/B,MAAM,CAAC,KAAK,CAAE,EAAgB,EAAE,EAAU,EAAE,EAAgB,EAAE,EAAU;QAC5E,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;QACtB,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IAAC,CAAC;IAIpB,GAAG,CAAE,CAAS;QAClB,OAAO,IAAI,cAAc,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAAC,CAAC;IAEhD,MAAM,CAAE,CAAS;QACrB,OAAO,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAAC,CAAC;IAExC,MAAM,CAAE,CAAS;QACrB,OAAO,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAAC,CAAC;IAIxC,QAAQ;QACZ,IAAI,CAAC,GAAG,GAAG,CAAC;QACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YACpC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;gBACT,CAAC,IAAI,IAAI,CAAC;YAAC,CAAC;YACf,CAAC,IAAI,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QAAC,CAAC;QACrD,CAAC,IAAI,GAAG,CAAC;QACT,OAAO,CAAC,CAAC;IAAC,CAAC;IAEP,WAAW;QACf,MAAM,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC;QACtB,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC1B,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QAAC,CAAC;QAC3B,OAAO,CAAC,CAAC;IAAC,CAAC;IAEP,WAAW;QACf,MAAM,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC;QACtB,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC1B,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QAAC,CAAC;QAC3B,OAAO,CAAC,CAAC;IAAC,CAAC;IAIP,SAAS,CAAE,CAAS,EAAE,CAAS;QACnC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;IAAC,CAAC;IAEd,KAAK,CAAE,CAAS,EAAE,CAAU;QAChC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC;QACnB,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC;IAAC,CAAC;IAEjB,WAAW,CAAE,CAAS,EAAE,CAAS;QACrC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;IAAC,CAAC;IAEd,OAAO,CAAE,CAAS,EAAE,CAAU;QAClC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC;QACnB,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC;IAAC,CAAC;IAEjB,SAAS,CAAE,CAAS,EAAE,CAAS;QACnC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;QAChB,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;IAAC,CAAC;IAEd,KAAK,CAAE,CAAS,EAAE,CAAU;QAChC,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;IAAC,CAAC;IAEjD,SAAS,CAAE,CAAS,EAAE,CAAS;QACnC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;QAChB,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;IAAC,CAAC;IAEd,KAAK,CAAE,CAAS,EAAE,CAAU;QAChC,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;IAAC,CAAC;IAMjD,UAAU,CAAE,EAAgB;QAChC,MAAM,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC;QACtB,MAAM,CAAC,EAAE,CAAC,MAAM,IAAI,CAAC,CAAC,CAAC;QACvB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC1B,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAAC,CAAC;IAAA,CAAC;IAE7D,YAAY,CAAE,CAAS;QAC3B,MAAM,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC;QACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC1B,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAAC,CAAC;IAAA,CAAC;IAOvB,MAAM,CAAE,CAAS,EAAE,GAAW,EAAE,GAAW,EAAE,GAAW,EAAE,GAAW;QACzE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;QACnC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;IAAC,CAAC;IAKjC,MAAM,CAAE,CAAS,EAAE,GAAW,EAAE,GAAW,EAAE,GAAW,EAAE,GAAW;QACzE,MAAM,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;QAChC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;QACzC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAAC,CAAC;CAE7C"}

package/math/MathUtils.d.ts

@@ -1,7 +1,10 @@
-export declare function fuzzyEquals(a: number, b: number, eps: number): boolean;
-export declare function isFuzzyInteger(i: number, eps: number): boolean;
-export declare function isPowerOf2(i: number): boolean;
-export declare function getNextPowerOf2(x: number): number;
-export declare function floorLog2(x: number): number;
-export declare function hyperbolicDecline(t: number, a: number, b: number): number;
-export declare function movingAverage(a: ArrayLike<number>, windowWidth: number): Float64Array;
+export declare function fuzzyEquals(a: number, b: number, eps: number): boolean;
+export declare function isFuzzyInteger(i: number, eps: number): boolean;
+export declare function isPowerOf2(i: number): boolean;
+export declare function getNextPowerOf2(x: number): number;
+export declare function floorLog2(x: number): number;
+export declare function hyperbolicDecline(t: number, a: number, b: number): number;
+export declare function simpleMovingAverage(a: ArrayLike<number>, windowWidth: number, shift?: boolean): Float64Array;
+export declare function triangularMovingAverage(a: ArrayLike<number>, windowWidth: number): Float64Array;
+export declare function triangularMovingAverageRef(a: ArrayLike<number>, windowWidth: number): Float64Array;
+//# sourceMappingURL=MathUtils.d.ts.map

package/math/MathUtils.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"MathUtils.d.ts","sourceRoot":"","sources":["../../src/math/MathUtils.ts"],"names":[],"mappings":"AAAA,wBAAgB,WAAW,CAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAI,OAAO,CAS3C;AAE9B,wBAAgB,cAAc,CAAE,CAAC,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAI,OAAO,CAClB;AAK/C,wBAAgB,UAAU,CAAE,CAAC,EAAE,MAAM,GAAI,OAAO,CAGjB;AAK/B,wBAAgB,eAAe,CAAE,CAAC,EAAE,MAAM,GAAI,MAAM,CAMtC;AAMd,wBAAgB,SAAS,CAAE,CAAC,EAAE,MAAM,GAAI,MAAM,CAGf;AAqB/B,wBAAgB,iBAAiB,CAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAI,MAAM,CASzB;AASnD,wBAAgB,mBAAmB,CAAE,CAAC,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,WAAW,EAAE,MAAM,EAAE,KAAK,UAAQ,GAAI,YAAY,CAoB9F;AAUf,wBAAgB,uBAAuB,CAAE,CAAC,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,WAAW,EAAE,MAAM,GAAI,YAAY,CAQnF;AAQf,wBAAgB,0BAA0B,CAAE,CAAC,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,WAAW,EAAE,MAAM,GAAI,YAAY,CAOtF"}

package/math/MathUtils.js

@@ -1,82 +1,117 @@
-export function fuzzyEquals(a, b, eps) {
-    if (!isFinite(a) || !isFinite(b)) {
-        return false;
-    }
-    if (a == b) {
-        return true;
-    }
-    const diff = Math.abs(a - b);
-    if (diff <= eps) {
-        return true;
-    }
-    const mag = Math.max(Math.abs(a), Math.abs(b));
-    return diff <= mag * eps;
-}
-export function isFuzzyInteger(i, eps) {
-    return Math.abs(i - Math.round(i)) <= eps;
-}
-export function isPowerOf2(i) {
-    if (!Number.isSafeInteger(i) || i < 1 || i > 0x40000000) {
-        return false;
-    }
-    return (i & (i - 1)) == 0;
-}
-export function getNextPowerOf2(x) {
-    if (!isFinite(x)) {
-        return NaN;
-    }
-    let n = 1;
-    while (n <= x) {
-        n *= 2;
-    }
-    return n;
-}
-export function floorLog2(x) {
-    if (x > 0x7FFFFFFF || x < 1) {
-        throw new Error("Argument is not a valid integer.");
-    }
-    return 31 - Math.clz32(x);
-}
-export function hyperbolicDecline(t, a, b) {
-    switch (b) {
-        case 1: {
-            return 1 / (1 + a * t);
-        }
-        case 0: {
-            return Math.exp(-a * t);
-        }
-        case -1: {
-            return Math.max(0, 1 - a * t);
-        }
-        default: {
-            return 1 / (1 + b * a * t) ** (1 / b);
-        }
-    }
-}
-export function movingAverage(a, windowWidth) {
-    if (windowWidth < 2 || !Number.isSafeInteger(windowWidth)) {
-        throw new Error("Specified window width is not a valid integer value.");
-    }
-    const len = a.length;
-    const a2 = new Float64Array(len);
-    const halfWindowWidth = Math.floor(windowWidth / 2);
-    const extendedLen = len + halfWindowWidth;
-    let movingSum = 0;
-    let n = 0;
-    for (let p = 0; p < extendedLen; p++) {
-        if (p >= windowWidth) {
-            movingSum -= a[p - windowWidth];
-            n--;
-        }
-        if (p < len) {
-            movingSum += a[p];
-            n++;
-        }
-        const p2 = p - halfWindowWidth;
-        if (p2 >= 0) {
-            a2[p2] = movingSum / n;
-        }
-    }
-    return a2;
-}
+export function fuzzyEquals(a, b, eps) {
+    if (!isFinite(a) || !isFinite(b)) {
+        return false;
+    }
+    if (a == b) {
+        return true;
+    }
+    const diff = Math.abs(a - b);
+    if (diff <= eps) {
+        return true;
+    }
+    const mag = Math.max(Math.abs(a), Math.abs(b));
+    return diff <= mag * eps;
+}
+export function isFuzzyInteger(i, eps) {
+    return Math.abs(i - Math.round(i)) <= eps;
+}
+export function isPowerOf2(i) {
+    if (!Number.isSafeInteger(i) || i < 1 || i > 0x40000000) {
+        return false;
+    }
+    return (i & (i - 1)) == 0;
+}
+export function getNextPowerOf2(x) {
+    if (!isFinite(x)) {
+        return NaN;
+    }
+    let n = 1;
+    while (n <= x) {
+        n *= 2;
+    }
+    return n;
+}
+export function floorLog2(x) {
+    if (x > 0x7FFFFFFF || x < 1) {
+        throw new Error("Argument is not a valid integer.");
+    }
+    return 31 - Math.clz32(x);
+}
+export function hyperbolicDecline(t, a, b) {
+    switch (b) {
+        case 1: {
+            return 1 / (1 + a * t);
+        }
+        case 0: {
+            return Math.exp(-a * t);
+        }
+        case -1: {
+            return Math.max(0, 1 - a * t);
+        }
+        default: {
+            return 1 / (1 + b * a * t) ** (1 / b);
+        }
+    }
+}
+export function simpleMovingAverage(a, windowWidth, shift = false) {
+    if (windowWidth < 2 || !Number.isSafeInteger(windowWidth)) {
+        throw new Error("Specified window width is not valid for SMA.");
+    }
+    const len = a.length;
+    const a2 = new Float64Array(len);
+    const halfWindowWidth = Math.floor(windowWidth / 2);
+    const posShift = (shift && windowWidth % 2 == 0) ? 1 : 0;
+    const extendedLen = len + halfWindowWidth;
+    let movingSum = 0;
+    let n = 0;
+    for (let p = 0; p < extendedLen; p++) {
+        if (p >= windowWidth) {
+            movingSum -= a[p - windowWidth];
+            n--;
+        }
+        if (p < len) {
+            movingSum += a[p];
+            n++;
+        }
+        const p2 = p - halfWindowWidth + posShift;
+        if (p2 >= 0 && p2 < len) {
+            a2[p2] = movingSum / n;
+        }
+    }
+    return a2;
+}
+export function triangularMovingAverage(a, windowWidth) {
+    if (windowWidth < 4 || !Number.isSafeInteger(windowWidth)) {
+        throw new Error("Specified window width is not valid for TMA.");
+    }
+    const w1 = Math.floor(windowWidth / 2);
+    const w2 = windowWidth - w1;
+    const a1 = simpleMovingAverage(a, w1);
+    const a2 = simpleMovingAverage(a1, w2, true);
+    return a2;
+}
+export function triangularMovingAverageRef(a, windowWidth) {
+    if (windowWidth < 4 || !Number.isSafeInteger(windowWidth)) {
+        throw new Error("Specified window width is not valid for TMA.");
+    }
+    const len = a.length;
+    const a2 = new Float64Array(len);
+    for (let p = 0; p < len; p++) {
+        a2[p] = computeTriangularAverageAt(a, p, windowWidth);
+    }
+    return a2;
+}
+function computeTriangularAverageAt(a, p, windowWidth) {
+    const len = a.length;
+    const p1 = Math.max(0, Math.ceil(p - windowWidth / 2 + 0.1));
+    const p2 = Math.min(len - 1, Math.floor(p + windowWidth / 2 - 0.1));
+    let sum = 0;
+    let weightSum = 0;
+    for (let i = p1; i <= p2; i++) {
+        const weight = 1 - Math.abs(i - p) / (windowWidth / 2);
+        sum += a[i] * weight;
+        weightSum += weight;
+    }
+    return sum / weightSum;
+}
 //# sourceMappingURL=MathUtils.js.map

package/math/MathUtils.js.map

@@ -1 +1 @@
-{"version":3,"file":"MathUtils.js","sourceRoot":"","sources":["../../src/math/MathUtils.ts"],"names":[],"mappings":"AAAA,MAAM,UAAU,WAAW,CAAE,CAAS,EAAE,CAAS,EAAE,GAAW;IAC3D,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE;QAC/B,OAAO,KAAK,CAAC;KAAE;IAClB,IAAI,CAAC,IAAI,CAAC,EAAE;QACT,OAAO,IAAI,CAAC;KAAE;IACjB,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC7B,IAAI,IAAI,IAAI,GAAG,EAAE;QACd,OAAO,IAAI,CAAC;KAAE;IACjB,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAC/C,OAAO,IAAI,IAAI,GAAG,GAAG,GAAG,CAAC;AAAC,CAAC;AAE9B,MAAM,UAAU,cAAc,CAAE,CAAS,EAAE,GAAW;IACnD,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC;AAAC,CAAC;AAK/C,MAAM,UAAU,UAAU,CAAE,CAAS;IAClC,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,UAAU,EAAE;QACtD,OAAO,KAAK,CAAC;KAAE;IAClB,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;AAAC,CAAC;AAK/B,MAAM,UAAU,eAAe,CAAE,CAAS;IACvC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE;QACf,OAAO,GAAG,CAAC;KAAE;IAChB,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,OAAO,CAAC,IAAI,CAAC,EAAE;QACZ,CAAC,IAAI,CAAC,CAAC;KAAE;IACZ,OAAO,CAAC,CAAC;AAAC,CAAC;AAMd,MAAM,UAAU,SAAS,CAAE,CAAS;IACjC,IAAI,CAAC,GAAG,UAAU,IAAI,CAAC,GAAG,CAAC,EAAE;QAC1B,MAAM,IAAI,KAAK,CAAC,kCAAkC,CAAC,CAAC;KAAE;IACzD,OAAO,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;AAAC,CAAC;AAqB/B,MAAM,UAAU,iBAAiB,CAAE,CAAS,EAAE,CAAS,EAAE,CAAS;IAC/D,QAAQ,CAAC,EAAE;QACR,KAAK,CAAC,CAAC,CAAC;YACL,OAAO,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;SAAE;QAC5B,KAAK,CAAC,CAAC,CAAC;YACL,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;SAAE;QAC7B,KAAK,CAAC,CAAC,CAAC,CAAC;YACN,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;SAAE;QACnC,OAAO,CAAC,CAAC;YACN,OAAO,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;SAAE;KAAC;AAAA,CAAC;AAMnD,MAAM,UAAU,aAAa,CAAE,CAAoB,EAAE,WAAmB;IACrE,IAAI,WAAW,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,WAAW,CAAC,EAAE;QACxD,MAAM,IAAI,KAAK,CAAC,sDAAsD,CAAC,CAAC;KAAE;IAC7E,MAAM,GAAG,GAAG,CAAC,CAAC,MAAM,CAAC;IACrB,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC;IACjC,MAAM,eAAe,GAAG,IAAI,CAAC,KAAK,CAAC,WAAW,GAAG,CAAC,CAAC,CAAC;IACpD,MAAM,WAAW,GAAG,GAAG,GAAG,eAAe,CAAC;IAC1C,IAAI,SAAS,GAAG,CAAC,CAAC;IAClB,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,WAAW,EAAE,CAAC,EAAE,EAAE;QACnC,IAAI,CAAC,IAAI,WAAW,EAAE;YACnB,SAAS,IAAI,CAAC,CAAC,CAAC,GAAG,WAAW,CAAC,CAAC;YAChC,CAAC,EAAE,CAAC;SAAE;QACT,IAAI,CAAC,GAAG,GAAG,EAAE;YACV,SAAS,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;YAClB,CAAC,EAAE,CAAC;SAAE;QACT,MAAM,EAAE,GAAG,CAAC,GAAG,eAAe,CAAC;QAC/B,IAAI,EAAE,IAAI,CAAC,EAAE;YACV,EAAE,CAAC,EAAE,CAAC,GAAG,SAAS,GAAG,CAAC,CAAC;SAAE;KAAC;IAChC,OAAO,EAAE,CAAC;AAAC,CAAC"}
+{"version":3,"file":"MathUtils.js","sourceRoot":"","sources":["../../src/math/MathUtils.ts"],"names":[],"mappings":"AAAA,MAAM,UAAU,WAAW,CAAE,CAAS,EAAE,CAAS,EAAE,GAAW;IAC3D,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;QAChC,OAAO,KAAK,CAAC;IAAC,CAAC;IAClB,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QACV,OAAO,IAAI,CAAC;IAAC,CAAC;IACjB,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC7B,IAAI,IAAI,IAAI,GAAG,EAAE,CAAC;QACf,OAAO,IAAI,CAAC;IAAC,CAAC;IACjB,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAC/C,OAAO,IAAI,IAAI,GAAG,GAAG,GAAG,CAAC;AAAC,CAAC;AAE9B,MAAM,UAAU,cAAc,CAAE,CAAS,EAAE,GAAW;IACnD,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC;AAAC,CAAC;AAK/C,MAAM,UAAU,UAAU,CAAE,CAAS;IAClC,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,UAAU,EAAE,CAAC;QACvD,OAAO,KAAK,CAAC;IAAC,CAAC;IAClB,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;AAAC,CAAC;AAK/B,MAAM,UAAU,eAAe,CAAE,CAAS;IACvC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;QAChB,OAAO,GAAG,CAAC;IAAC,CAAC;IAChB,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,OAAO,CAAC,IAAI,CAAC,EAAE,CAAC;QACb,CAAC,IAAI,CAAC,CAAC;IAAC,CAAC;IACZ,OAAO,CAAC,CAAC;AAAC,CAAC;AAMd,MAAM,UAAU,SAAS,CAAE,CAAS;IACjC,IAAI,CAAC,GAAG,UAAU,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QAC3B,MAAM,IAAI,KAAK,CAAC,kCAAkC,CAAC,CAAC;IAAC,CAAC;IACzD,OAAO,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;AAAC,CAAC;AAqB/B,MAAM,UAAU,iBAAiB,CAAE,CAAS,EAAE,CAAS,EAAE,CAAS;IAC/D,QAAQ,CAAC,EAAE,CAAC;QACT,KAAK,CAAC,CAAC,CAAC,CAAC;YACN,OAAO,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAAC,CAAC;QAC5B,KAAK,CAAC,CAAC,CAAC,CAAC;YACN,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAAC,CAAC;QAC7B,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;YACP,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAAC,CAAC;QACnC,OAAO,CAAC,CAAC,CAAC;YACP,OAAO,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAAC,CAAC;IAAA,CAAC;AAAA,CAAC;AASnD,MAAM,UAAU,mBAAmB,CAAE,CAAoB,EAAE,WAAmB,EAAE,KAAK,GAAG,KAAK;IAC1F,IAAI,WAAW,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,WAAW,CAAC,EAAE,CAAC;QACzD,MAAM,IAAI,KAAK,CAAC,8CAA8C,CAAC,CAAC;IAAC,CAAC;IACrE,MAAM,GAAG,GAAG,CAAC,CAAC,MAAM,CAAC;IACrB,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC;IACjC,MAAM,eAAe,GAAG,IAAI,CAAC,KAAK,CAAC,WAAW,GAAG,CAAC,CAAC,CAAC;IACpD,MAAM,QAAQ,GAAG,CAAC,KAAK,IAAI,WAAW,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACzD,MAAM,WAAW,GAAG,GAAG,GAAG,eAAe,CAAC;IAC1C,IAAI,SAAS,GAAG,CAAC,CAAC;IAClB,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,WAAW,EAAE,CAAC,EAAE,EAAE,CAAC;QACpC,IAAI,CAAC,IAAI,WAAW,EAAE,CAAC;YACpB,SAAS,IAAI,CAAC,CAAC,CAAC,GAAG,WAAW,CAAC,CAAC;YAChC,CAAC,EAAE,CAAC;QAAC,CAAC;QACT,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC;YACX,SAAS,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;YAClB,CAAC,EAAE,CAAC;QAAC,CAAC;QACT,MAAM,EAAE,GAAG,CAAC,GAAG,eAAe,GAAG,QAAQ,CAAC;QAC1C,IAAI,EAAE,IAAI,CAAC,IAAI,EAAE,GAAG,GAAG,EAAE,CAAC;YACvB,EAAE,CAAC,EAAE,CAAC,GAAG,SAAS,GAAG,CAAC,CAAC;QAAC,CAAC;IAAA,CAAC;IAChC,OAAO,EAAE,CAAC;AAAC,CAAC;AAUf,MAAM,UAAU,uBAAuB,CAAE,CAAoB,EAAE,WAAmB;IAC/E,IAAI,WAAW,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,WAAW,CAAC,EAAE,CAAC;QACzD,MAAM,IAAI,KAAK,CAAC,8CAA8C,CAAC,CAAC;IAAC,CAAC;IAErE,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,WAAW,GAAG,CAAC,CAAC,CAAC;IACvC,MAAM,EAAE,GAAG,WAAW,GAAG,EAAE,CAAC;IAC5B,MAAM,EAAE,GAAG,mBAAmB,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;IACtC,MAAM,EAAE,GAAG,mBAAmB,CAAC,EAAE,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;IAC7C,OAAO,EAAE,CAAC;AAAC,CAAC;AAQf,MAAM,UAAU,0BAA0B,CAAE,CAAoB,EAAE,WAAmB;IAClF,IAAI,WAAW,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,WAAW,CAAC,EAAE,CAAC;QACzD,MAAM,IAAI,KAAK,CAAC,8CAA8C,CAAC,CAAC;IAAC,CAAC;IACrE,MAAM,GAAG,GAAG,CAAC,CAAC,MAAM,CAAC;IACrB,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC;IACjC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;QAC5B,EAAE,CAAC,CAAC,CAAC,GAAG,0BAA0B,CAAC,CAAC,EAAE,CAAC,EAAE,WAAW,CAAC,CAAC;IAAC,CAAC;IAC3D,OAAO,EAAE,CAAC;AAAC,CAAC;AAGf,SAAS,0BAA0B,CAAE,CAAoB,EAAE,CAAS,EAAE,WAAmB;IACtF,MAAM,GAAG,GAAG,CAAC,CAAC,MAAM,CAAC;IACrB,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,WAAW,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;IAC7D,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,WAAW,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;IACpE,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,IAAI,SAAS,GAAG,CAAC,CAAC;IAClB,KAAK,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;QAC7B,MAAM,MAAM,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,WAAW,GAAG,CAAC,CAAC,CAAC;QACvD,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;QACrB,SAAS,IAAI,MAAM,CAAC;IAAC,CAAC;IACzB,OAAO,GAAG,GAAG,SAAS,CAAC;AAAC,CAAC"}