dsp-collection 0.2.5 → 0.2.7

package/math/PolyReal.js

@@ -1,227 +1,227 @@
-import MutableComplex from "./MutableComplex.js";
-export function evaluateReal(a, x) {
-    if (a.length == 0) {
-        throw new Error("Zero length array.");
-    }
-    const n = a.length - 1;
-    let r = a[n];
-    for (let i = n - 1; i >= 0; i--) {
-        r *= x;
-        r += a[i];
-    }
-    return r;
-}
-export function evaluateComplex(a, x) {
-    if (a.length == 0) {
-        throw new Error("Zero length array.");
-    }
-    const n = a.length - 1;
-    const r = new MutableComplex(a[n]);
-    for (let i = n - 1; i >= 0; i--) {
-        r.mulBy(x);
-        r.addRealTo(a[i]);
-    }
-    return r;
-}
-export function expand(zeros) {
-    const n = zeros.length;
-    if (n == 0) {
-        return Float64Array.of(1);
-    }
-    let a = Float64Array.of(-zeros[0], 1);
-    for (let i = 1; i < n; i++) {
-        const a2 = Float64Array.of(-zeros[i], 1);
-        a = multiply(a, a2);
-    }
-    return a;
-}
-export function compareEqual(a1, a2, eps = 0) {
-    const n1 = a1.length - 1;
-    const n2 = a2.length - 1;
-    const n = Math.max(n1, n2);
-    for (let i = 0; i <= n; i++) {
-        const v1 = (i <= n1) ? a1[i] : 0;
-        const v2 = (i <= n2) ? a2[i] : 0;
-        if (Math.abs(v1 - v2) > eps) {
-            return false;
-        }
-    }
-    return true;
-}
-export function add(a1, a2, eps = 0) {
-    const n1 = a1.length - 1;
-    const n2 = a2.length - 1;
-    const n3 = Math.max(n1, n2);
-    const a3 = new Float64Array(n3 + 1);
-    for (let i = 0; i <= n3; i++) {
-        const v1 = (i <= n1) ? a1[i] : 0;
-        const v2 = (i <= n2) ? a2[i] : 0;
-        a3[i] = v1 + v2;
-    }
-    return trim(a3, eps);
-}
-export function multiply(a1, a2, eps = 0) {
-    if (a1.length == 0 || a2.length == 0) {
-        throw new Error("Zero length arrays.");
-    }
-    if (a1.length == 1 && a1[0] == 0 || a2.length == 1 && a2[0] == 0) {
-        return Float64Array.of(0);
-    }
-    const n1 = a1.length - 1;
-    const n2 = a2.length - 1;
-    const n3 = n1 + n2;
-    const a3 = new Float64Array(n3 + 1);
-    for (let i = 0; i <= n3; i++) {
-        let t = 0;
-        const p1 = Math.max(0, i - n2);
-        const p2 = Math.min(n1, i);
-        for (let j = p1; j <= p2; j++) {
-            t += a1[j] * a2[i - j];
-        }
-        a3[i] = t;
-    }
-    return trim(a3, eps);
-}
-export function divide(a1r, a2r, eps = 0) {
-    if (a1r.length == 0 || a2r.length == 0) {
-        throw new Error("Zero length arrays.");
-    }
-    const a1 = trim(a1r, eps);
-    const a2 = trim(a2r, eps);
-    if (a2.length == 1) {
-        if (a2[0] == 0) {
-            throw new Error("Polynomial division by zero.");
-        }
-        if (a2[0] == 1) {
-            return [Float64Array.from(a1), Float64Array.of(0)];
-        }
-        return [divByReal(a1, a2[0]), Float64Array.of(0)];
-    }
-    const n1 = a1.length - 1;
-    const n2 = a2.length - 1;
-    if (n1 < n2) {
-        return [Float64Array.of(0), Float64Array.from(a1)];
-    }
-    const a = Float64Array.from(a1);
-    const lc2 = a2[n2];
-    for (let i = n1 - n2; i >= 0; i--) {
-        const r = a[n2 + i] / lc2;
-        a[n2 + i] = r;
-        for (let j = 0; j < n2; ++j) {
-            a[i + j] -= r * a2[j];
-        }
-    }
-    const quotient = trim(a.subarray(n2), eps);
-    const remainder = trim(a.subarray(0, n2), eps);
-    return [quotient, remainder];
-}
-export function gcd(a1, a2, eps = 0) {
-    let r1 = trim(a1, eps);
-    let r2 = trim(a2, eps);
-    makeMonic(r1);
-    makeMonic(r2);
-    if (r1.length < r2.length) {
-        [r1, r2] = [r2, r1];
-    }
-    while (true) {
-        if (r2.length < 2) {
-            return Float64Array.of(1);
-        }
-        const r = divide(r1, r2, eps)[1];
-        if (r.length == 1 && r[0] == 0) {
-            return r2;
-        }
-        makeMonic(r);
-        r1 = r2;
-        r2 = r;
-    }
-}
-function trim(a, eps = 0) {
-    if (a.length == 0) {
-        throw new Error("Zero length array.");
-    }
-    if (Math.abs(a[a.length - 1]) > eps) {
-        return Float64Array.from(a);
-    }
-    let len = a.length - 1;
-    while (len > 0 && Math.abs(a[len - 1]) <= eps) {
-        len--;
-    }
-    if (len == 0) {
-        return Float64Array.of(0);
-    }
-    const a2 = new Float64Array(len);
-    for (let i = 0; i < len; i++) {
-        a2[i] = a[i];
-    }
-    return a2;
-}
-function makeMonic(a) {
-    const len = a.length;
-    if (len == 0) {
-        throw new Error("Zero length array.");
-    }
-    const lc = a[len - 1];
-    if (lc == 1) {
-        return;
-    }
-    if (lc == 0) {
-        throw new Error("Leading coefficient is zero.");
-    }
-    a[len - 1] = 1;
-    for (let i = 0; i < len - 1; i++) {
-        a[i] /= lc;
-    }
-}
-function divByReal(a, b) {
-    const a2 = new Float64Array(a.length);
-    for (let i = 0; i < a.length; i++) {
-        a2[i] = a[i] / b;
-    }
-    return a2;
-}
-function divByRealInPlace(a, b) {
-    for (let i = 0; i < a.length; i++) {
-        a[i] /= b;
-    }
-}
-export function evaluateFractionComplex(f, x) {
-    const v1 = evaluateComplex(f[0], x);
-    const v2 = evaluateComplex(f[1], x);
-    return v1.div(v2);
-}
-export function addFractions(f1, f2, eps = 0) {
-    if (compareEqual(f1[1], f2[1], eps)) {
-        return [add(f1[0], f2[0], eps), Float64Array.from(f1[1])];
-    }
-    const g = gcd(f1[1], f2[1], eps);
-    if (g.length == 1 && g[0] == 1) {
-        const top = add(multiply(f1[0], f2[1], eps), multiply(f2[0], f1[1], eps));
-        const bottom = multiply(f1[1], f2[1], eps);
-        return [top, bottom];
-    }
-    const q1 = divide(f1[1], g, eps);
-    const q2 = divide(f2[1], g, eps);
-    const m1 = q1[0];
-    const m2 = q2[0];
-    const top = add(multiply(f1[0], m2, eps), multiply(f2[0], m1, eps));
-    const bottom = multiply(f1[1], m2, eps);
-    return [top, bottom];
-}
-export function multiplyFractions(f1, f2, eps = 0) {
-    const top = multiply(f1[0], f2[0], eps);
-    const bottom = multiply(f1[1], f2[1], eps);
-    return [top, bottom];
-}
-export function normalizeFraction(f, eps = 0) {
-    const top = trim(f[0], eps);
-    const bottom = trim(f[1], eps);
-    const lc = bottom[bottom.length - 1];
-    if (lc == 0) {
-        throw new Error("Fraction denominator is zero.");
-    }
-    divByRealInPlace(top, lc);
-    divByRealInPlace(bottom, lc);
-    return [top, bottom];
-}
+import MutableComplex from "./MutableComplex.js";
+export function evaluateReal(a, x) {
+    if (a.length == 0) {
+        throw new Error("Zero length array.");
+    }
+    const n = a.length - 1;
+    let r = a[n];
+    for (let i = n - 1; i >= 0; i--) {
+        r *= x;
+        r += a[i];
+    }
+    return r;
+}
+export function evaluateComplex(a, x) {
+    if (a.length == 0) {
+        throw new Error("Zero length array.");
+    }
+    const n = a.length - 1;
+    const r = new MutableComplex(a[n]);
+    for (let i = n - 1; i >= 0; i--) {
+        r.mulBy(x);
+        r.addRealTo(a[i]);
+    }
+    return r;
+}
+export function expand(zeros) {
+    const n = zeros.length;
+    if (n == 0) {
+        return Float64Array.of(1);
+    }
+    let a = Float64Array.of(-zeros[0], 1);
+    for (let i = 1; i < n; i++) {
+        const a2 = Float64Array.of(-zeros[i], 1);
+        a = multiply(a, a2);
+    }
+    return a;
+}
+export function compareEqual(a1, a2, eps = 0) {
+    const n1 = a1.length - 1;
+    const n2 = a2.length - 1;
+    const n = Math.max(n1, n2);
+    for (let i = 0; i <= n; i++) {
+        const v1 = (i <= n1) ? a1[i] : 0;
+        const v2 = (i <= n2) ? a2[i] : 0;
+        if (Math.abs(v1 - v2) > eps) {
+            return false;
+        }
+    }
+    return true;
+}
+export function add(a1, a2, eps = 0) {
+    const n1 = a1.length - 1;
+    const n2 = a2.length - 1;
+    const n3 = Math.max(n1, n2);
+    const a3 = new Float64Array(n3 + 1);
+    for (let i = 0; i <= n3; i++) {
+        const v1 = (i <= n1) ? a1[i] : 0;
+        const v2 = (i <= n2) ? a2[i] : 0;
+        a3[i] = v1 + v2;
+    }
+    return trim(a3, eps);
+}
+export function multiply(a1, a2, eps = 0) {
+    if (a1.length == 0 || a2.length == 0) {
+        throw new Error("Zero length arrays.");
+    }
+    if (a1.length == 1 && a1[0] == 0 || a2.length == 1 && a2[0] == 0) {
+        return Float64Array.of(0);
+    }
+    const n1 = a1.length - 1;
+    const n2 = a2.length - 1;
+    const n3 = n1 + n2;
+    const a3 = new Float64Array(n3 + 1);
+    for (let i = 0; i <= n3; i++) {
+        let t = 0;
+        const p1 = Math.max(0, i - n2);
+        const p2 = Math.min(n1, i);
+        for (let j = p1; j <= p2; j++) {
+            t += a1[j] * a2[i - j];
+        }
+        a3[i] = t;
+    }
+    return trim(a3, eps);
+}
+export function divide(a1r, a2r, eps = 0) {
+    if (a1r.length == 0 || a2r.length == 0) {
+        throw new Error("Zero length arrays.");
+    }
+    const a1 = trim(a1r, eps);
+    const a2 = trim(a2r, eps);
+    if (a2.length == 1) {
+        if (a2[0] == 0) {
+            throw new Error("Polynomial division by zero.");
+        }
+        if (a2[0] == 1) {
+            return [Float64Array.from(a1), Float64Array.of(0)];
+        }
+        return [divByReal(a1, a2[0]), Float64Array.of(0)];
+    }
+    const n1 = a1.length - 1;
+    const n2 = a2.length - 1;
+    if (n1 < n2) {
+        return [Float64Array.of(0), Float64Array.from(a1)];
+    }
+    const a = Float64Array.from(a1);
+    const lc2 = a2[n2];
+    for (let i = n1 - n2; i >= 0; i--) {
+        const r = a[n2 + i] / lc2;
+        a[n2 + i] = r;
+        for (let j = 0; j < n2; ++j) {
+            a[i + j] -= r * a2[j];
+        }
+    }
+    const quotient = trim(a.subarray(n2), eps);
+    const remainder = trim(a.subarray(0, n2), eps);
+    return [quotient, remainder];
+}
+export function gcd(a1, a2, eps = 0) {
+    let r1 = trim(a1, eps);
+    let r2 = trim(a2, eps);
+    makeMonic(r1);
+    makeMonic(r2);
+    if (r1.length < r2.length) {
+        [r1, r2] = [r2, r1];
+    }
+    while (true) {
+        if (r2.length < 2) {
+            return Float64Array.of(1);
+        }
+        const r = divide(r1, r2, eps)[1];
+        if (r.length == 1 && r[0] == 0) {
+            return r2;
+        }
+        makeMonic(r);
+        r1 = r2;
+        r2 = r;
+    }
+}
+function trim(a, eps = 0) {
+    if (a.length == 0) {
+        throw new Error("Zero length array.");
+    }
+    if (Math.abs(a[a.length - 1]) > eps) {
+        return Float64Array.from(a);
+    }
+    let len = a.length - 1;
+    while (len > 0 && Math.abs(a[len - 1]) <= eps) {
+        len--;
+    }
+    if (len == 0) {
+        return Float64Array.of(0);
+    }
+    const a2 = new Float64Array(len);
+    for (let i = 0; i < len; i++) {
+        a2[i] = a[i];
+    }
+    return a2;
+}
+function makeMonic(a) {
+    const len = a.length;
+    if (len == 0) {
+        throw new Error("Zero length array.");
+    }
+    const lc = a[len - 1];
+    if (lc == 1) {
+        return;
+    }
+    if (lc == 0) {
+        throw new Error("Leading coefficient is zero.");
+    }
+    a[len - 1] = 1;
+    for (let i = 0; i < len - 1; i++) {
+        a[i] /= lc;
+    }
+}
+function divByReal(a, b) {
+    const a2 = new Float64Array(a.length);
+    for (let i = 0; i < a.length; i++) {
+        a2[i] = a[i] / b;
+    }
+    return a2;
+}
+function divByRealInPlace(a, b) {
+    for (let i = 0; i < a.length; i++) {
+        a[i] /= b;
+    }
+}
+export function evaluateFractionComplex(f, x) {
+    const v1 = evaluateComplex(f[0], x);
+    const v2 = evaluateComplex(f[1], x);
+    return v1.div(v2);
+}
+export function addFractions(f1, f2, eps = 0) {
+    if (compareEqual(f1[1], f2[1], eps)) {
+        return [add(f1[0], f2[0], eps), Float64Array.from(f1[1])];
+    }
+    const g = gcd(f1[1], f2[1], eps);
+    if (g.length == 1 && g[0] == 1) {
+        const top = add(multiply(f1[0], f2[1], eps), multiply(f2[0], f1[1], eps));
+        const bottom = multiply(f1[1], f2[1], eps);
+        return [top, bottom];
+    }
+    const q1 = divide(f1[1], g, eps);
+    const q2 = divide(f2[1], g, eps);
+    const m1 = q1[0];
+    const m2 = q2[0];
+    const top = add(multiply(f1[0], m2, eps), multiply(f2[0], m1, eps));
+    const bottom = multiply(f1[1], m2, eps);
+    return [top, bottom];
+}
+export function multiplyFractions(f1, f2, eps = 0) {
+    const top = multiply(f1[0], f2[0], eps);
+    const bottom = multiply(f1[1], f2[1], eps);
+    return [top, bottom];
+}
+export function normalizeFraction(f, eps = 0) {
+    const top = trim(f[0], eps);
+    const bottom = trim(f[1], eps);
+    const lc = bottom[bottom.length - 1];
+    if (lc == 0) {
+        throw new Error("Fraction denominator is zero.");
+    }
+    divByRealInPlace(top, lc);
+    divByRealInPlace(bottom, lc);
+    return [top, bottom];
+}
 //# sourceMappingURL=PolyReal.js.map

package/math/PolyReal.js.map

@@ -1 +1 @@
-{"version":3,"file":"PolyReal.js","sourceRoot":"","sources":["../../src/math/PolyReal.ts"],"names":[],"mappings":"AAiBA,OAAO,cAAc,MAAM,qBAAqB,CAAC;AAYjD,MAAM,UAAU,YAAY,CAAE,CAAoB,EAAE,CAAS;IAC1D,IAAI,CAAC,CAAC,MAAM,IAAI,CAAC,EAAE;QAChB,MAAM,IAAI,KAAK,CAAC,oBAAoB,CAAC,CAAC;KAAE;IAC3C,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;IACvB,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACb,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE;QAC9B,CAAC,IAAI,CAAC,CAAC;QACP,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;KAAE;IACf,OAAO,CAAC,CAAC;AAAC,CAAC;AAYd,MAAM,UAAU,eAAe,CAAE,CAAoB,EAAE,CAAU;IAC9D,IAAI,CAAC,CAAC,MAAM,IAAI,CAAC,EAAE;QAChB,MAAM,IAAI,KAAK,CAAC,oBAAoB,CAAC,CAAC;KAAE;IAC3C,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;IACvB,MAAM,CAAC,GAAG,IAAI,cAAc,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACnC,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE;QAC9B,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QACX,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;KAAE;IACvB,OAAO,CAAC,CAAC;AAAC,CAAC;AAKd,MAAM,UAAU,MAAM,CAAE,KAAwB;IAC7C,MAAM,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC;IACvB,IAAI,CAAC,IAAI,CAAC,EAAE;QACT,OAAO,YAAY,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;KAAE;IAC/B,IAAI,CAAC,GAAG,YAAY,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACtC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;QACzB,MAAM,EAAE,GAAG,YAAY,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACzC,CAAC,GAAG,QAAQ,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;KAAE;IACzB,OAAO,CAAC,CAAC;AAAC,CAAC;AAKd,MAAM,UAAU,YAAY,CAAE,EAAqB,EAAE,EAAqB,EAAE,GAAG,GAAG,CAAC;IAChF,MAAM,EAAE,GAAG,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC;IACzB,MAAM,EAAE,GAAG,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC;IACzB,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE;QAC1B,MAAM,EAAE,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACjC,MAAM,EAAE,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACjC,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,GAAG,EAAE;YAC1B,OAAO,KAAK,CAAC;SAAE;KAAC;IACtB,OAAO,IAAI,CAAC;AAAC,CAAC;AAKjB,MAAM,UAAU,GAAG,CAAE,EAAqB,EAAE,EAAqB,EAAE,GAAG,GAAG,CAAC;IACvE,MAAM,EAAE,GAAG,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC;IACzB,MAAM,EAAE,GAAG,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC;IACzB,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IAC5B,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;IACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE,EAAE;QAC3B,MAAM,EAAE,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACjC,MAAM,EAAE,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACjC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC;KAAE;IACrB,OAAO,IAAI,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;AAAC,CAAC;AAK1B,MAAM,UAAU,QAAQ,CAAE,EAAqB,EAAE,EAAqB,EAAE,GAAG,GAAG,CAAC;IAC5E,IAAI,EAAE,CAAC,MAAM,IAAI,CAAC,IAAI,EAAE,CAAC,MAAM,IAAI,CAAC,EAAE;QACnC,MAAM,IAAI,KAAK,CAAC,qBAAqB,CAAC,CAAC;KAAE;IAC5C,IAAI,EAAE,CAAC,MAAM,IAAI,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,EAAE,CAAC,MAAM,IAAI,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE;QAC/D,OAAO,YAAY,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;KAAE;IAC/B,MAAM,EAAE,GAAG,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC;IACzB,MAAM,EAAE,GAAG,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC;IACzB,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;IACnB,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;IACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE,EAAE;QAC3B,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC;QAC/B,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;QAC3B,KAAK,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE,EAAE;YAC5B,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;SAAE;QAC5B,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;KAAE;IACf,OAAO,IAAI,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;AAAC,CAAC;AAM1B,MAAM,UAAU,MAAM,CAAE,GAAsB,EAAE,GAAsB,EAAE,GAAG,GAAG,CAAC;IAC5E,IAAI,GAAG,CAAC,MAAM,IAAI,CAAC,IAAI,GAAG,CAAC,MAAM,IAAI,CAAC,EAAE;QACrC,MAAM,IAAI,KAAK,CAAC,qBAAqB,CAAC,CAAC;KAAE;IAC5C,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;IAC1B,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;IAC1B,IAAI,EAAE,CAAC,MAAM,IAAI,CAAC,EAAE;QACjB,IAAI,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE;YACb,MAAM,IAAI,KAAK,CAAC,8BAA8B,CAAC,CAAC;SAAE;QACrD,IAAI,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE;YACb,OAAO,CAAC,YAAY,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,YAAY,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;SAAE;QACxD,OAAO,CAAC,SAAS,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,YAAY,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;KAAE;IACvD,MAAM,EAAE,GAAG,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC;IACzB,MAAM,EAAE,GAAG,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC;IACzB,IAAI,EAAE,GAAG,EAAE,EAAE;QACV,OAAO,CAAC,YAAY,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,YAAY,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;KAAE;IACxD,MAAM,CAAC,GAAG,YAAY,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IAChC,MAAM,GAAG,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC;IACnB,KAAK,IAAI,CAAC,GAAG,EAAE,GAAG,EAAE,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE;QAChC,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;QAC1B,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,EAAE,CAAC,EAAE;YAC1B,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;SAAE;KAAC;IAC/B,MAAM,QAAQ,GAAG,IAAI,CAAC,CAAC,CAAC,QAAQ,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;IAC3C,MAAM,SAAS,GAAG,IAAI,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;IAC/C,OAAO,CAAC,QAAQ,EAAE,SAAS,CAAC,CAAC;AAAC,CAAC;AAKlC,MAAM,UAAU,GAAG,CAAE,EAAqB,EAAE,EAAqB,EAAE,GAAG,GAAG,CAAC;IACvE,IAAI,EAAE,GAAG,IAAI,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;IACvB,IAAI,EAAE,GAAG,IAAI,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;IACvB,SAAS,CAAC,EAAE,CAAC,CAAC;IACd,SAAS,CAAC,EAAE,CAAC,CAAC;IACd,IAAI,EAAE,CAAC,MAAM,GAAG,EAAE,CAAC,MAAM,EAAE;QACxB,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;KAAE;IACzB,OAAO,IAAI,EAAE;QACV,IAAI,EAAE,CAAC,MAAM,GAAG,CAAC,EAAE;YAChB,OAAO,YAAY,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;SAAE;QAC/B,MAAM,CAAC,GAAG,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACjC,IAAI,CAAC,CAAC,MAAM,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE;YAC7B,OAAO,EAAE,CAAC;SAAE;QACf,SAAS,CAAC,CAAC,CAAC,CAAC;QACb,EAAE,GAAG,EAAE,CAAC;QACR,EAAE,GAAG,CAAC,CAAC;KAAE;AAAA,CAAC;AAGhB,SAAS,IAAI,CAAE,CAAoB,EAAE,GAAG,GAAG,CAAC;IACzC,IAAI,CAAC,CAAC,MAAM,IAAI,CAAC,EAAE;QAChB,MAAM,IAAI,KAAK,CAAC,oBAAoB,CAAC,CAAC;KAAE;IAC3C,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,EAAE;QAClC,OAAO,YAAY,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;KAAE;IACjC,IAAI,GAAG,GAAG,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;IACvB,OAAO,GAAG,GAAG,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,IAAI,GAAG,EAAE;QAC5C,GAAG,EAAE,CAAC;KAAE;IACX,IAAI,GAAG,IAAI,CAAC,EAAE;QACX,OAAO,YAAY,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;KAAE;IAC/B,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC;IACjC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE;QAC3B,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;KAAE;IAClB,OAAO,EAAE,CAAC;AAAC,CAAC;AAGf,SAAS,SAAS,CAAE,CAAe;IAChC,MAAM,GAAG,GAAG,CAAC,CAAC,MAAM,CAAC;IACrB,IAAI,GAAG,IAAI,CAAC,EAAE;QACX,MAAM,IAAI,KAAK,CAAC,oBAAoB,CAAC,CAAC;KAAE;IAC3C,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;IACtB,IAAI,EAAE,IAAI,CAAC,EAAE;QACV,OAAO;KAAE;IACZ,IAAI,EAAE,IAAI,CAAC,EAAE;QACV,MAAM,IAAI,KAAK,CAAC,8BAA8B,CAAC,CAAC;KAAE;IACrD,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;IACf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;QAC/B,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC;KAAE;AAAA,CAAC;AAEpB,SAAS,SAAS,CAAE,CAAoB,EAAE,CAAS;IAChD,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;IACtC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QAChC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;KAAE;IACtB,OAAO,EAAE,CAAC;AAAC,CAAC;AAEf,SAAS,gBAAgB,CAAE,CAAe,EAAE,CAAS;IAClD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QAChC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;KAAE;AAAA,CAAC;AAOnB,MAAM,UAAU,uBAAuB,CAAE,CAAsB,EAAE,CAAU;IACxE,MAAM,EAAE,GAAG,eAAe,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACpC,MAAM,EAAE,GAAG,eAAe,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACpC,OAAO,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;AAAC,CAAC;AAKvB,MAAM,UAAU,YAAY,CAAE,EAAuB,EAAE,EAAuB,EAAE,GAAG,GAAG,CAAC;IACpF,IAAI,YAAY,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE;QAClC,OAAO,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE,YAAY,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;KAAE;IAC/D,MAAM,CAAC,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;IACjC,IAAI,CAAC,CAAC,MAAM,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE;QAC7B,MAAM,GAAG,GAAG,GAAG,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QAC1E,MAAM,MAAM,GAAG,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;QAC3C,OAAO,CAAC,GAAG,EAAE,MAAM,CAAC,CAAC;KAAE;IAC1B,MAAM,EAAE,GAAG,MAAM,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;IACjC,MAAM,EAAE,GAAG,MAAM,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;IAKjC,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IACjB,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IACjB,MAAM,GAAG,GAAG,GAAG,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,GAAG,CAAC,EAAE,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC;IACpE,MAAM,MAAM,GAAG,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC;IACxC,OAAO,CAAC,GAAG,EAAE,MAAM,CAAC,CAAC;AAAC,CAAC;AAK1B,MAAM,UAAU,iBAAiB,CAAE,EAAuB,EAAE,EAAuB,EAAE,GAAG,GAAG,CAAC;IACzF,MAAM,GAAG,GAAM,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;IAC3C,MAAM,MAAM,GAAG,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;IAC3C,OAAO,CAAC,GAAG,EAAE,MAAM,CAAC,CAAC;AAAC,CAAC;AAK1B,MAAM,UAAU,iBAAiB,CAAE,CAAsB,EAAE,GAAG,GAAG,CAAC;IAC/D,MAAM,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;IAC5B,MAAM,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;IAC/B,MAAM,EAAE,GAAG,MAAM,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IACrC,IAAI,EAAE,IAAI,CAAC,EAAE;QACV,MAAM,IAAI,KAAK,CAAC,+BAA+B,CAAC,CAAC;KAAE;IACtD,gBAAgB,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;IAC1B,gBAAgB,CAAC,MAAM,EAAE,EAAE,CAAC,CAAC;IAC7B,OAAO,CAAC,GAAG,EAAE,MAAM,CAAC,CAAC;AAAC,CAAC"}
+{"version":3,"file":"PolyReal.js","sourceRoot":"","sources":["../../src/math/PolyReal.ts"],"names":[],"mappings":"AAiBA,OAAO,cAAc,MAAM,qBAAqB,CAAC;AAYjD,MAAM,UAAU,YAAY,CAAE,CAAoB,EAAE,CAAS;IAC1D,IAAI,CAAC,CAAC,MAAM,IAAI,CAAC,EAAE,CAAC;QACjB,MAAM,IAAI,KAAK,CAAC,oBAAoB,CAAC,CAAC;IAAC,CAAC;IAC3C,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;IACvB,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACb,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,CAAC,IAAI,CAAC,CAAC;QACP,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IAAC,CAAC;IACf,OAAO,CAAC,CAAC;AAAC,CAAC;AAYd,MAAM,UAAU,eAAe,CAAE,CAAoB,EAAE,CAAU;IAC9D,IAAI,CAAC,CAAC,MAAM,IAAI,CAAC,EAAE,CAAC;QACjB,MAAM,IAAI,KAAK,CAAC,oBAAoB,CAAC,CAAC;IAAC,CAAC;IAC3C,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;IACvB,MAAM,CAAC,GAAG,IAAI,cAAc,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACnC,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QACX,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAAC,CAAC;IACvB,OAAO,CAAC,CAAC;AAAC,CAAC;AAKd,MAAM,UAAU,MAAM,CAAE,KAAwB;IAC7C,MAAM,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC;IACvB,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QACV,OAAO,YAAY,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAAC,CAAC;IAC/B,IAAI,CAAC,GAAiB,YAAY,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACpD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC1B,MAAM,EAAE,GAAG,YAAY,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACzC,CAAC,GAAG,QAAQ,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;IAAC,CAAC;IACzB,OAAO,CAAC,CAAC;AAAC,CAAC;AAKd,MAAM,UAAU,YAAY,CAAE,EAAqB,EAAE,EAAqB,EAAE,GAAG,GAAG,CAAC;IAChF,MAAM,EAAE,GAAG,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC;IACzB,MAAM,EAAE,GAAG,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC;IACzB,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,EAAE,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACjC,MAAM,EAAE,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACjC,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC;YAC3B,OAAO,KAAK,CAAC;QAAC,CAAC;IAAA,CAAC;IACtB,OAAO,IAAI,CAAC;AAAC,CAAC;AAKjB,MAAM,UAAU,GAAG,CAAE,EAAqB,EAAE,EAAqB,EAAE,GAAG,GAAG,CAAC;IACvE,MAAM,EAAE,GAAG,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC;IACzB,MAAM,EAAE,GAAG,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC;IACzB,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IAC5B,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;IACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;QAC5B,MAAM,EAAE,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACjC,MAAM,EAAE,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACjC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC;IAAC,CAAC;IACrB,OAAO,IAAI,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;AAAC,CAAC;AAK1B,MAAM,UAAU,QAAQ,CAAE,EAAqB,EAAE,EAAqB,EAAE,GAAG,GAAG,CAAC;IAC5E,IAAI,EAAE,CAAC,MAAM,IAAI,CAAC,IAAI,EAAE,CAAC,MAAM,IAAI,CAAC,EAAE,CAAC;QACpC,MAAM,IAAI,KAAK,CAAC,qBAAqB,CAAC,CAAC;IAAC,CAAC;IAC5C,IAAI,EAAE,CAAC,MAAM,IAAI,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,EAAE,CAAC,MAAM,IAAI,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,CAAC;QAChE,OAAO,YAAY,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAAC,CAAC;IAC/B,MAAM,EAAE,GAAG,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC;IACzB,MAAM,EAAE,GAAG,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC;IACzB,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;IACnB,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;IACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;QAC5B,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC;QAC/B,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;QAC3B,KAAK,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;YAC7B,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAAC,CAAC;QAC5B,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAAC,CAAC;IACf,OAAO,IAAI,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;AAAC,CAAC;AAM1B,MAAM,UAAU,MAAM,CAAE,GAAsB,EAAE,GAAsB,EAAE,GAAG,GAAG,CAAC;IAC5E,IAAI,GAAG,CAAC,MAAM,IAAI,CAAC,IAAI,GAAG,CAAC,MAAM,IAAI,CAAC,EAAE,CAAC;QACtC,MAAM,IAAI,KAAK,CAAC,qBAAqB,CAAC,CAAC;IAAC,CAAC;IAC5C,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;IAC1B,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;IAC1B,IAAI,EAAE,CAAC,MAAM,IAAI,CAAC,EAAE,CAAC;QAClB,IAAI,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,CAAC;YACd,MAAM,IAAI,KAAK,CAAC,8BAA8B,CAAC,CAAC;QAAC,CAAC;QACrD,IAAI,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,CAAC;YACd,OAAO,CAAC,YAAY,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,YAAY,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAAC,CAAC;QACxD,OAAO,CAAC,SAAS,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,YAAY,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAAC,CAAC;IACvD,MAAM,EAAE,GAAG,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC;IACzB,MAAM,EAAE,GAAG,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC;IACzB,IAAI,EAAE,GAAG,EAAE,EAAE,CAAC;QACX,OAAO,CAAC,YAAY,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,YAAY,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;IAAC,CAAC;IACxD,MAAM,CAAC,GAAG,YAAY,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IAChC,MAAM,GAAG,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC;IACnB,KAAK,IAAI,CAAC,GAAG,EAAE,GAAG,EAAE,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QACjC,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;QAC1B,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC;YAC3B,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;QAAC,CAAC;IAAA,CAAC;IAC/B,MAAM,QAAQ,GAAG,IAAI,CAAC,CAAC,CAAC,QAAQ,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;IAC3C,MAAM,SAAS,GAAG,IAAI,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;IAC/C,OAAO,CAAC,QAAQ,EAAE,SAAS,CAAC,CAAC;AAAC,CAAC;AAKlC,MAAM,UAAU,GAAG,CAAE,EAAqB,EAAE,EAAqB,EAAE,GAAG,GAAG,CAAC;IACvE,IAAI,EAAE,GAAG,IAAI,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;IACvB,IAAI,EAAE,GAAG,IAAI,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;IACvB,SAAS,CAAC,EAAE,CAAC,CAAC;IACd,SAAS,CAAC,EAAE,CAAC,CAAC;IACd,IAAI,EAAE,CAAC,MAAM,GAAG,EAAE,CAAC,MAAM,EAAE,CAAC;QACzB,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IAAC,CAAC;IACzB,OAAO,IAAI,EAAE,CAAC;QACX,IAAI,EAAE,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;YACjB,OAAO,YAAY,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAAC,CAAC;QAC/B,MAAM,CAAC,GAAG,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACjC,IAAI,CAAC,CAAC,MAAM,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,CAAC;YAC9B,OAAO,EAAE,CAAC;QAAC,CAAC;QACf,SAAS,CAAC,CAAC,CAAC,CAAC;QACb,EAAE,GAAG,EAAE,CAAC;QACR,EAAE,GAAG,CAAC,CAAC;IAAC,CAAC;AAAA,CAAC;AAGhB,SAAS,IAAI,CAAE,CAAoB,EAAE,GAAG,GAAG,CAAC;IACzC,IAAI,CAAC,CAAC,MAAM,IAAI,CAAC,EAAE,CAAC;QACjB,MAAM,IAAI,KAAK,CAAC,oBAAoB,CAAC,CAAC;IAAC,CAAC;IAC3C,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC;QACnC,OAAO,YAAY,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAAC,CAAC;IACjC,IAAI,GAAG,GAAG,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;IACvB,OAAO,GAAG,GAAG,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,IAAI,GAAG,EAAE,CAAC;QAC7C,GAAG,EAAE,CAAC;IAAC,CAAC;IACX,IAAI,GAAG,IAAI,CAAC,EAAE,CAAC;QACZ,OAAO,YAAY,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAAC,CAAC;IAC/B,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,CAAC,CAAC,CAAC,CAAC,CAAC;IAAC,CAAC;IAClB,OAAO,EAAE,CAAC;AAAC,CAAC;AAGf,SAAS,SAAS,CAAE,CAAe;IAChC,MAAM,GAAG,GAAG,CAAC,CAAC,MAAM,CAAC;IACrB,IAAI,GAAG,IAAI,CAAC,EAAE,CAAC;QACZ,MAAM,IAAI,KAAK,CAAC,oBAAoB,CAAC,CAAC;IAAC,CAAC;IAC3C,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;IACtB,IAAI,EAAE,IAAI,CAAC,EAAE,CAAC;QACX,OAAO;IAAC,CAAC;IACZ,IAAI,EAAE,IAAI,CAAC,EAAE,CAAC;QACX,MAAM,IAAI,KAAK,CAAC,8BAA8B,CAAC,CAAC;IAAC,CAAC;IACrD,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;IACf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAChC,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC;IAAC,CAAC;AAAA,CAAC;AAEpB,SAAS,SAAS,CAAE,CAAoB,EAAE,CAAS;IAChD,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;IACtC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACjC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAAC,CAAC;IACtB,OAAO,EAAE,CAAC;AAAC,CAAC;AAEf,SAAS,gBAAgB,CAAE,CAAe,EAAE,CAAS;IAClD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACjC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;IAAC,CAAC;AAAA,CAAC;AAOnB,MAAM,UAAU,uBAAuB,CAAE,CAAsB,EAAE,CAAU;IACxE,MAAM,EAAE,GAAG,eAAe,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACpC,MAAM,EAAE,GAAG,eAAe,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACpC,OAAO,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;AAAC,CAAC;AAKvB,MAAM,UAAU,YAAY,CAAE,EAAuB,EAAE,EAAuB,EAAE,GAAG,GAAG,CAAC;IACpF,IAAI,YAAY,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE,CAAC;QACnC,OAAO,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE,YAAY,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAAC,CAAC;IAC/D,MAAM,CAAC,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;IACjC,IAAI,CAAC,CAAC,MAAM,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,CAAC;QAC9B,MAAM,GAAG,GAAG,GAAG,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QAC1E,MAAM,MAAM,GAAG,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;QAC3C,OAAO,CAAC,GAAG,EAAE,MAAM,CAAC,CAAC;IAAC,CAAC;IAC1B,MAAM,EAAE,GAAG,MAAM,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;IACjC,MAAM,EAAE,GAAG,MAAM,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;IAKjC,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IACjB,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IACjB,MAAM,GAAG,GAAG,GAAG,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,GAAG,CAAC,EAAE,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC;IACpE,MAAM,MAAM,GAAG,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC;IACxC,OAAO,CAAC,GAAG,EAAE,MAAM,CAAC,CAAC;AAAC,CAAC;AAK1B,MAAM,UAAU,iBAAiB,CAAE,EAAuB,EAAE,EAAuB,EAAE,GAAG,GAAG,CAAC;IACzF,MAAM,GAAG,GAAM,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;IAC3C,MAAM,MAAM,GAAG,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;IAC3C,OAAO,CAAC,GAAG,EAAE,MAAM,CAAC,CAAC;AAAC,CAAC;AAK1B,MAAM,UAAU,iBAAiB,CAAE,CAAsB,EAAE,GAAG,GAAG,CAAC;IAC/D,MAAM,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;IAC5B,MAAM,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;IAC/B,MAAM,EAAE,GAAG,MAAM,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IACrC,IAAI,EAAE,IAAI,CAAC,EAAE,CAAC;QACX,MAAM,IAAI,KAAK,CAAC,+BAA+B,CAAC,CAAC;IAAC,CAAC;IACtD,gBAAgB,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;IAC1B,gBAAgB,CAAC,MAAM,EAAE,EAAE,CAAC,CAAC;IAC7B,OAAO,CAAC,GAAG,EAAE,MAAM,CAAC,CAAC;AAAC,CAAC"}

package/package.json

@@ -1,6 +1,6 @@
 {
    "name": "dsp-collection",
-   "version": "0.2.5",
+   "version": "0.2.7",
    "description": "A collection of JavaScript modules for digital signal processing (written in TypeScript)",
    "keywords": [
       "DSP",
@@ -13,11 +13,16 @@
       "Spectral filtering",
       "Window functions",
       "Instantaneous frequency",
-      "Pitch estimation"
+      "Pitch estimation",
+      "Simple moving average (SMA)",
+      "triangular moving average (TMA)"
    ],
    "license": "MIT",
    "homepage": "http://www.source-code.biz/dsp/js",
-   "repository": "github:chdh/dsp-collection-js",
+   "repository": {
+      "type": "git",
+      "url": "git+https://github.com/chdh/dsp-collection-js.git"
+   },
    "contributors": [
       {
          "name": "Christian d'Heureuse",
@@ -32,5 +37,8 @@
       "doPack": "cd dist && npm pack",
       "doPublish": "cd dist && npm publish"
    },
+   "devDependencies": {
+      "@types/node": "^22"
+   },
    "type": "module"
 }

package/signal/AdaptiveStft.d.ts

@@ -1,12 +1,13 @@
-import Complex from "../math/Complex.js";
-import * as WindowFunctions from "./WindowFunctions.js";
-export interface ComponentResult {
-    component: Complex;
-    frequency: number;
-    windowStartPosition: number;
-    windowWidth: number;
-}
-export declare function getSingle_relWindow(samples: Float64Array | Float32Array, roughFrequency: number, roughWindowCenterPosition: number, relWindowWidth: number, windowFunction: WindowFunctions.WindowFunction | undefined): ComponentResult | undefined;
-export declare function getSingle_maxWindow(samples: Float64Array | Float32Array, roughFrequency: number, roughWindowCenterPosition: number, maxWindowWidth: number, windowFunction: WindowFunctions.WindowFunction | undefined): ComponentResult | undefined;
-export declare function getSingle(samples: Float64Array | Float32Array, roughFrequency: number, windowFunction: WindowFunctions.WindowFunction | undefined): Complex;
-export declare function getHarmonicAmplitudes(samples: Float64Array | Float32Array, windowCenterPosition: number, f0: number, harmonics: number, relWindowWidth?: number, windowFunction?: typeof WindowFunctions.flatTopWindowNorm): Float64Array | undefined;
+import Complex from "../math/Complex.ts";
+import * as WindowFunctions from "./WindowFunctions.ts";
+export interface ComponentResult {
+    component: Complex;
+    frequency: number;
+    windowStartPosition: number;
+    windowWidth: number;
+}
+export declare function getSingle_relWindow(samples: Float64Array | Float32Array, roughFrequency: number, roughWindowCenterPosition: number, relWindowWidth: number, windowFunction: WindowFunctions.WindowFunction | undefined): ComponentResult | undefined;
+export declare function getSingle_maxWindow(samples: Float64Array | Float32Array, roughFrequency: number, roughWindowCenterPosition: number, maxWindowWidth: number, windowFunction: WindowFunctions.WindowFunction | undefined): ComponentResult | undefined;
+export declare function getSingle(samples: Float64Array | Float32Array, roughFrequency: number, windowFunction: WindowFunctions.WindowFunction | undefined): Complex;
+export declare function getHarmonicAmplitudes(samples: Float64Array | Float32Array, windowCenterPosition: number, f0: number, harmonics: number, relWindowWidth?: number, windowFunction?: typeof WindowFunctions.flatTopWindowNorm): Float64Array | undefined;
+//# sourceMappingURL=AdaptiveStft.d.ts.map

package/signal/AdaptiveStft.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"AdaptiveStft.d.ts","sourceRoot":"","sources":["../../src/signal/AdaptiveStft.ts"],"names":[],"mappings":"AASA,OAAO,OAAO,MAAM,oBAAoB,CAAC;AAEzC,OAAO,KAAK,eAAe,MAAM,sBAAsB,CAAC;AAExD,MAAM,WAAW,eAAe;IAC7B,SAAS,EAAiB,OAAO,CAAC;IAClC,SAAS,EAAiB,MAAM,CAAC;IACjC,mBAAmB,EAAO,MAAM,CAAC;IACjC,WAAW,EAAe,MAAM,CAAC;CAAE;AAyBtC,wBAAgB,mBAAmB,CAAE,OAAO,EAAE,YAAY,GAAG,YAAY,EAAE,cAAc,EAAE,MAAM,EAAE,yBAAyB,EAAE,MAAM,EAC9H,cAAc,EAAE,MAAM,EAAE,cAAc,EAAE,eAAe,CAAC,cAAc,GAAG,SAAS,GAAI,eAAe,GAAG,SAAS,CAkB1E;AAwB7C,wBAAgB,mBAAmB,CAAE,OAAO,EAAE,YAAY,GAAG,YAAY,EAAE,cAAc,EAAE,MAAM,EAAE,yBAAyB,EAAE,MAAM,EAC9H,cAAc,EAAE,MAAM,EAAE,cAAc,EAAE,eAAe,CAAC,cAAc,GAAG,SAAS,GAAI,eAAe,GAAG,SAAS,CAGH;AAiBpH,wBAAgB,SAAS,CAAE,OAAO,EAAE,YAAY,GAAG,YAAY,EAAE,cAAc,EAAE,MAAM,EAAE,cAAc,EAAE,eAAe,CAAC,cAAc,GAAG,SAAS,GAAI,OAAO,CAEpH;AAuB1C,wBAAgB,qBAAqB,CAAE,OAAO,EAAE,YAAY,GAAG,YAAY,EAAE,oBAAoB,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,cAAc,SAAI,EAAE,cAAc,2CAAoC,GAAI,YAAY,GAAG,SAAS,CAgBrN"}