dsp-collection 0.2.2 → 0.2.5

package/math/PolyReal.js

@@ -1,39 +1,76 @@
-import MutableComplex from "./MutableComplex";
-export function evaluate(a, x) {
+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;
-    const sum = new MutableComplex(a[n]);
+    let r = a[n];
     for (let i = n - 1; i >= 0; i--) {
-        sum.mulBy(x);
-        sum.addRealTo(a[i]);
+        r *= x;
+        r += a[i];
     }
-    return sum;
+    return r;
 }
-export function add(a1, a2) {
+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 Array(n3 + 1);
+    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 a3;
+    return trim(a3, eps);
 }
-export function multiply(a1, a2) {
+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 [0];
+        return Float64Array.of(0);
     }
     const n1 = a1.length - 1;
     const n2 = a2.length - 1;
     const n3 = n1 + n2;
-    const a3 = new Array(n3 + 1);
+    const a3 = new Float64Array(n3 + 1);
     for (let i = 0; i <= n3; i++) {
         let t = 0;
         const p1 = Math.max(0, i - n2);
@@ -43,37 +80,148 @@ export function multiply(a1, a2) {
         }
         a3[i] = t;
     }
-    return a3;
+    return trim(a3, eps);
 }
-export function equals(a1, a2, eps = 0) {
+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;
-    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;
+    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];
         }
     }
-    return true;
+    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;
+    }
 }
-export function evaluateFraction(f, x) {
-    const v1 = evaluate(f[0], x);
-    const v2 = evaluate(f[1], x);
+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) {
-    if (equals(f1[1], f2[1])) {
-        return [add(f1[0], f2[0]), f1[1].slice()];
+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 top = add(multiply(f1[0], f2[1]), multiply(f2[0], f1[1]));
-    const bottom = multiply(f1[1], f2[1]);
+    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) {
-    const top = multiply(f1[0], f2[0]);
-    const bottom = multiply(f1[1], f2[1]);
+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":"AAeA,OAAO,cAAc,MAAM,kBAAkB,CAAC;AAY9C,MAAM,UAAU,QAAQ,CAAE,CAAW,EAAE,CAAU;IAC9C,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,GAAG,GAAG,IAAI,cAAc,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACrC,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE;QAC9B,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QACb,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;KAAE;IACzB,OAAO,GAAG,CAAC;AAAC,CAAC;AAKhB,MAAM,UAAU,GAAG,CAAE,EAAY,EAAE,EAAY;IAC5C,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,KAAK,CAAS,EAAE,GAAG,CAAC,CAAC,CAAC;IACrC,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,EAAE,CAAC;AAAC,CAAC;AAKf,MAAM,UAAU,QAAQ,CAAE,EAAY,EAAE,EAAY;IACjD,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,CAAC,CAAC,CAAC,CAAC;KAAE;IAChB,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,KAAK,CAAS,EAAE,GAAG,CAAC,CAAC,CAAC;IACrC,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,EAAE,CAAC;AAAC,CAAC;AAKf,MAAM,UAAU,MAAM,CAAE,EAAY,EAAE,EAAY,EAAE,GAAG,GAAG,CAAC;IACxD,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;AAOjB,MAAM,UAAU,gBAAgB,CAAE,CAAa,EAAE,CAAU;IACxD,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAC7B,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAC7B,OAAO,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;AAAC,CAAC;AAKvB,MAAM,UAAU,YAAY,CAAE,EAAc,EAAE,EAAc;IACzD,IAAI,MAAM,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE;QACvB,OAAO,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC;KAAE;IAC/C,MAAM,GAAG,GAAG,GAAG,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAChE,MAAM,MAAM,GAAG,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACtC,OAAO,CAAC,GAAG,EAAE,MAAM,CAAC,CAAC;AAAC,CAAC;AAK1B,MAAM,UAAU,iBAAiB,CAAE,EAAc,EAAE,EAAc;IAC9D,MAAM,GAAG,GAAM,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACtC,MAAM,MAAM,GAAG,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACtC,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;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"}

package/package.json

@@ -1,6 +1,6 @@
 {
    "name": "dsp-collection",
-   "version": "0.2.2",
+   "version": "0.2.5",
    "description": "A collection of JavaScript modules for digital signal processing (written in TypeScript)",
    "keywords": [
       "DSP",
@@ -9,6 +9,8 @@
       "Goertzel",
       "FFT",
       "Fourier Transform",
+      "Filter",
+      "Spectral filtering",
       "Window functions",
       "Instantaneous frequency",
       "Pitch estimation"
@@ -19,7 +21,8 @@
    "contributors": [
       {
          "name": "Christian d'Heureuse",
-         "email": "chdh@inventec.ch"
+         "email": "chdh@inventec.ch",
+         "url": "http://www.source-code.biz"
       }
    ],
    "scripts": {
@@ -29,10 +32,5 @@
       "doPack": "cd dist && npm pack",
       "doPublish": "cd dist && npm publish"
    },
-   "devDependencies": {
-      "@types/node": "^11.13.4",
-      "rimraf": "^2.6.3",
-      "source-map-support": "^0.5.12"
-   },
    "type": "module"
 }

package/signal/AdaptiveStft.d.ts

@@ -1,12 +1,12 @@
-import Complex from "../math/Complex";
-import * as WindowFunctions from "./WindowFunctions";
+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, roughFrequency: number, roughWindowCenterPosition: number, relWindowWidth: number, windowFunction: WindowFunctions.WindowFunction | undefined): ComponentResult | undefined;
-export declare function getSingle_maxWindow(samples: Float64Array, roughFrequency: number, roughWindowCenterPosition: number, maxWindowWidth: number, windowFunction: WindowFunctions.WindowFunction | undefined): ComponentResult | undefined;
-export declare function getSingle(samples: Float64Array, roughFrequency: number, windowFunction: WindowFunctions.WindowFunction | undefined): Complex;
-export declare function getHarmonicAmplitudes(samples: Float64Array, windowCenterPosition: number, f0: number, harmonics: number, relWindowWidth?: number, windowFunction?: typeof WindowFunctions.flatTopWindowNorm): Float64Array | undefined;
+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;

package/signal/AdaptiveStft.js

@@ -1,6 +1,6 @@
-import Complex from "../math/Complex";
-import * as Goertzel from "./Goertzel";
-import * as WindowFunctions from "./WindowFunctions";
+import Complex from "../math/Complex.js";
+import * as Goertzel from "./Goertzel.js";
+import * as WindowFunctions from "./WindowFunctions.js";
 export function getSingle_relWindow(samples, roughFrequency, roughWindowCenterPosition, relWindowWidth, windowFunction) {
     if (roughFrequency <= 0 || roughFrequency >= 0.5 || relWindowWidth < 1) {
         return;

package/signal/AdaptiveStft.js.map

@@ -1 +1 @@
-{"version":3,"file":"AdaptiveStft.js","sourceRoot":"","sources":["../../src/signal/AdaptiveStft.ts"],"names":[],"mappings":"AASA,OAAO,OAAO,MAAM,iBAAiB,CAAC;AACtC,OAAO,KAAK,QAAQ,MAAM,YAAY,CAAC;AACvC,OAAO,KAAK,eAAe,MAAM,mBAAmB,CAAC;AA+BrD,MAAM,UAAU,mBAAmB,CAAE,OAAqB,EAAE,cAAsB,EAAE,yBAAiC,EAC/G,cAAsB,EAAE,cAA0D;IACrF,IAAI,cAAc,IAAI,CAAC,IAAI,cAAc,IAAI,GAAG,IAAI,cAAc,GAAG,CAAC,EAAE;QACrE,OAAO;KAAE;IACZ,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,cAAc,CAAC,EAAE;QACpC,MAAM,IAAI,KAAK,CAAC,6CAA6C,CAAC,CAAC;KAAE;IACpE,MAAM,WAAW,GAAG,IAAI,CAAC,KAAK,CAAC,cAAc,GAAG,cAAc,CAAC,CAAC;IAChE,MAAM,mBAAmB,GAAG,IAAI,CAAC,KAAK,CAAC,yBAAyB,GAAG,WAAW,GAAG,CAAC,CAAC,CAAC;IACpF,IAAI,mBAAmB,GAAG,CAAC,IAAI,mBAAmB,GAAG,WAAW,GAAG,OAAO,CAAC,MAAM,EAAE;QAChF,OAAO;KAAE;IACZ,MAAM,aAAa,GAAG,OAAO,CAAC,QAAQ,CAAC,mBAAmB,EAAE,mBAAmB,GAAG,WAAW,CAAC,CAAC;IAC/F,MAAM,eAAe,GAAG,cAAc,CAAC,CAAC,CAAC,eAAe,CAAC,WAAW,CAAC,aAAa,EAAE,cAAc,CAAC,CAAC,CAAC,CAAC,aAAa,CAAC;IACpH,MAAM,UAAU,GAAG,QAAQ,CAAC,cAAc,CAAC,eAAe,EAAE,cAAc,CAAC,CAAC;IAC5E,MAAM,mBAAmB,GAAG,UAAU,CAAC,OAAO,CAAC,CAAC,GAAG,WAAW,CAAC,CAAC;IAChE,MAAM,sBAAsB,GAAG,CAAC,cAAc,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,mBAAmB,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,mBAAmB,CAAC;IAC3G,OAAO;QACJ,SAAS,EAAc,sBAAsB;QAC7C,SAAS,EAAc,cAAc,GAAG,WAAW;QACnD,mBAAmB,EAAI,mBAAmB;QAC1C,WAAW,EAAY,WAAW;KAAE,CAAC;AAAC,CAAC;AAwB7C,MAAM,UAAU,mBAAmB,CAAE,OAAqB,EAAE,cAAsB,EAAE,yBAAiC,EAC/G,cAAsB,EAAE,cAA0D;IACrF,MAAM,cAAc,GAAG,IAAI,CAAC,KAAK,CAAC,cAAc,GAAG,cAAc,CAAC,CAAC;IAEnE,OAAO,mBAAmB,CAAC,OAAO,EAAE,cAAc,EAAE,yBAAyB,EAAE,cAAc,EAAE,cAAc,CAAC,CAAC;AAAC,CAAC;AAiBpH,MAAM,UAAU,SAAS,CAAE,OAAqB,EAAE,cAAsB,EAAE,cAA0D;IACjI,MAAM,CAAC,GAAG,mBAAmB,CAAC,OAAO,EAAE,cAAc,EAAE,OAAO,CAAC,MAAM,GAAG,CAAC,EAAE,OAAO,CAAC,MAAM,EAAE,cAAc,CAAC,CAAC;IAC3G,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,OAAO,CAAC,GAAG,CAAC;AAAC,CAAC;AAuB1C,MAAM,UAAU,qBAAqB,CAAE,OAAqB,EAAE,oBAA4B,EAAE,EAAU,EAAE,SAAiB,EAAE,cAAc,GAAG,CAAC,EAAE,cAAc,GAAG,eAAe,CAAC,iBAAiB;IAC9L,IAAI,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,IAAI,CAAC,IAAI,SAAS,IAAI,CAAC,EAAE;QAC7C,OAAO;KAAE;IACZ,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,cAAc,CAAC,EAAE;QACpC,MAAM,IAAI,KAAK,CAAC,6CAA6C,CAAC,CAAC;KAAE;IACpE,MAAM,WAAW,GAAG,IAAI,CAAC,KAAK,CAAC,cAAc,GAAG,EAAE,CAAC,CAAC;IACpD,MAAM,mBAAmB,GAAG,IAAI,CAAC,KAAK,CAAC,oBAAoB,GAAG,WAAW,GAAG,CAAC,CAAC,CAAC;IAC/E,IAAI,mBAAmB,GAAG,CAAC,IAAI,mBAAmB,GAAG,WAAW,GAAG,OAAO,CAAC,MAAM,IAAI,WAAW,IAAI,CAAC,EAAE;QACpG,OAAO;KAAE;IACZ,MAAM,aAAa,GAAG,OAAO,CAAC,QAAQ,CAAC,mBAAmB,EAAE,mBAAmB,GAAG,WAAW,CAAC,CAAC;IAC/F,MAAM,eAAe,GAAG,eAAe,CAAC,WAAW,CAAC,aAAa,EAAE,cAAc,CAAC,CAAC;IACnF,MAAM,UAAU,GAAG,IAAI,YAAY,CAAC,SAAS,CAAC,CAAC;IAC/C,KAAK,IAAI,QAAQ,GAAG,CAAC,EAAE,QAAQ,IAAI,SAAS,EAAE,QAAQ,EAAE,EAAE;QACvD,MAAM,CAAC,GAAG,QAAQ,CAAC,cAAc,CAAC,eAAe,EAAE,QAAQ,GAAG,cAAc,CAAC,CAAC;QAC9E,MAAM,SAAS,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,WAAW,CAAC;QAC5C,UAAU,CAAC,QAAQ,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC;KAAE;IAC1C,OAAO,UAAU,CAAC;AAAC,CAAC"}
+{"version":3,"file":"AdaptiveStft.js","sourceRoot":"","sources":["../../src/signal/AdaptiveStft.ts"],"names":[],"mappings":"AASA,OAAO,OAAO,MAAM,oBAAoB,CAAC;AACzC,OAAO,KAAK,QAAQ,MAAM,eAAe,CAAC;AAC1C,OAAO,KAAK,eAAe,MAAM,sBAAsB,CAAC;AA+BxD,MAAM,UAAU,mBAAmB,CAAE,OAAoC,EAAE,cAAsB,EAAE,yBAAiC,EAC9H,cAAsB,EAAE,cAA0D;IACrF,IAAI,cAAc,IAAI,CAAC,IAAI,cAAc,IAAI,GAAG,IAAI,cAAc,GAAG,CAAC,EAAE;QACrE,OAAO;KAAE;IACZ,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,cAAc,CAAC,EAAE;QACpC,MAAM,IAAI,KAAK,CAAC,6CAA6C,CAAC,CAAC;KAAE;IACpE,MAAM,WAAW,GAAG,IAAI,CAAC,KAAK,CAAC,cAAc,GAAG,cAAc,CAAC,CAAC;IAChE,MAAM,mBAAmB,GAAG,IAAI,CAAC,KAAK,CAAC,yBAAyB,GAAG,WAAW,GAAG,CAAC,CAAC,CAAC;IACpF,IAAI,mBAAmB,GAAG,CAAC,IAAI,mBAAmB,GAAG,WAAW,GAAG,OAAO,CAAC,MAAM,EAAE;QAChF,OAAO;KAAE;IACZ,MAAM,aAAa,GAAG,OAAO,CAAC,QAAQ,CAAC,mBAAmB,EAAE,mBAAmB,GAAG,WAAW,CAAC,CAAC;IAC/F,MAAM,eAAe,GAAG,cAAc,CAAC,CAAC,CAAC,eAAe,CAAC,WAAW,CAAC,aAAa,EAAE,cAAc,CAAC,CAAC,CAAC,CAAC,aAAa,CAAC;IACpH,MAAM,UAAU,GAAG,QAAQ,CAAC,cAAc,CAAC,eAAe,EAAE,cAAc,CAAC,CAAC;IAC5E,MAAM,mBAAmB,GAAG,UAAU,CAAC,OAAO,CAAC,CAAC,GAAG,WAAW,CAAC,CAAC;IAChE,MAAM,sBAAsB,GAAG,CAAC,cAAc,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,mBAAmB,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,mBAAmB,CAAC;IAC3G,OAAO;QACJ,SAAS,EAAc,sBAAsB;QAC7C,SAAS,EAAc,cAAc,GAAG,WAAW;QACnD,mBAAmB,EAAI,mBAAmB;QAC1C,WAAW,EAAY,WAAW;KAAE,CAAC;AAAC,CAAC;AAwB7C,MAAM,UAAU,mBAAmB,CAAE,OAAoC,EAAE,cAAsB,EAAE,yBAAiC,EAC9H,cAAsB,EAAE,cAA0D;IACrF,MAAM,cAAc,GAAG,IAAI,CAAC,KAAK,CAAC,cAAc,GAAG,cAAc,CAAC,CAAC;IAEnE,OAAO,mBAAmB,CAAC,OAAO,EAAE,cAAc,EAAE,yBAAyB,EAAE,cAAc,EAAE,cAAc,CAAC,CAAC;AAAC,CAAC;AAiBpH,MAAM,UAAU,SAAS,CAAE,OAAoC,EAAE,cAAsB,EAAE,cAA0D;IAChJ,MAAM,CAAC,GAAG,mBAAmB,CAAC,OAAO,EAAE,cAAc,EAAE,OAAO,CAAC,MAAM,GAAG,CAAC,EAAE,OAAO,CAAC,MAAM,EAAE,cAAc,CAAC,CAAC;IAC3G,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,OAAO,CAAC,GAAG,CAAC;AAAC,CAAC;AAuB1C,MAAM,UAAU,qBAAqB,CAAE,OAAoC,EAAE,oBAA4B,EAAE,EAAU,EAAE,SAAiB,EAAE,cAAc,GAAG,CAAC,EAAE,cAAc,GAAG,eAAe,CAAC,iBAAiB;IAC7M,IAAI,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,IAAI,CAAC,IAAI,SAAS,IAAI,CAAC,EAAE;QAC7C,OAAO;KAAE;IACZ,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,cAAc,CAAC,EAAE;QACpC,MAAM,IAAI,KAAK,CAAC,6CAA6C,CAAC,CAAC;KAAE;IACpE,MAAM,WAAW,GAAG,IAAI,CAAC,KAAK,CAAC,cAAc,GAAG,EAAE,CAAC,CAAC;IACpD,MAAM,mBAAmB,GAAG,IAAI,CAAC,KAAK,CAAC,oBAAoB,GAAG,WAAW,GAAG,CAAC,CAAC,CAAC;IAC/E,IAAI,mBAAmB,GAAG,CAAC,IAAI,mBAAmB,GAAG,WAAW,GAAG,OAAO,CAAC,MAAM,IAAI,WAAW,IAAI,CAAC,EAAE;QACpG,OAAO;KAAE;IACZ,MAAM,aAAa,GAAG,OAAO,CAAC,QAAQ,CAAC,mBAAmB,EAAE,mBAAmB,GAAG,WAAW,CAAC,CAAC;IAC/F,MAAM,eAAe,GAAG,eAAe,CAAC,WAAW,CAAC,aAAa,EAAE,cAAc,CAAC,CAAC;IACnF,MAAM,UAAU,GAAG,IAAI,YAAY,CAAC,SAAS,CAAC,CAAC;IAC/C,KAAK,IAAI,QAAQ,GAAG,CAAC,EAAE,QAAQ,IAAI,SAAS,EAAE,QAAQ,EAAE,EAAE;QACvD,MAAM,CAAC,GAAG,QAAQ,CAAC,cAAc,CAAC,eAAe,EAAE,QAAQ,GAAG,cAAc,CAAC,CAAC;QAC9E,MAAM,SAAS,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,WAAW,CAAC;QAC5C,UAAU,CAAC,QAAQ,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC;KAAE;IAC1C,OAAO,UAAU,CAAC;AAAC,CAAC"}

package/signal/Autocorrelation.d.ts

@@ -1,5 +1,5 @@
-import * as WindowFunctions from "./WindowFunctions";
-export declare function nonPeriodicAutocorrelationSingle(x: Float64Array, distance: number, compensate: boolean): number;
-export declare function findNonPeriodicAutocorrelationMaximum(x: Float64Array, minDistance: number, maxDistance: number, fixedOverlapWidth: boolean): number;
-export declare function nonPeriodicAutocorrelation(x: Float64Array, normalize: boolean): Float64Array;
-export declare function windowedNonPeriodicAutocorrelation(x: Float64Array, windowFunction: WindowFunctions.WindowFunction, normalize: boolean): Float64Array;
+import * as WindowFunctions from "./WindowFunctions.js";
+export declare function nonPeriodicAutocorrelationSingle(x: ArrayLike<number>, distance: number, compensate: boolean): number;
+export declare function findNonPeriodicAutocorrelationMaximum(x: ArrayLike<number>, minDistance: number, maxDistance: number, fixedOverlapWidth: boolean): number;
+export declare function nonPeriodicAutocorrelation(x: ArrayLike<number>, normalize: boolean): Float64Array;
+export declare function windowedNonPeriodicAutocorrelation(x: ArrayLike<number>, windowFunction: WindowFunctions.WindowFunction, normalize: boolean): Float64Array;

package/signal/Autocorrelation.js

@@ -1,5 +1,5 @@
-import * as WindowFunctions from "./WindowFunctions";
-import * as ArrayUtils from "../utils/ArrayUtils";
+import * as WindowFunctions from "./WindowFunctions.js";
+import * as ArrayUtils from "../utils/ArrayUtils.js";
 function autocorrelationKernel(x, distance, n) {
     let sum = 0;
     for (let i = 0; i < n; i++) {

package/signal/Autocorrelation.js.map

@@ -1 +1 @@
-{"version":3,"file":"Autocorrelation.js","sourceRoot":"","sources":["../../src/signal/Autocorrelation.ts"],"names":[],"mappings":"AAIA,OAAO,KAAK,eAAe,MAAM,mBAAmB,CAAC;AACrD,OAAO,KAAK,UAAU,MAAM,qBAAqB,CAAC;AAElD,SAAS,qBAAqB,CAAE,CAAe,EAAE,QAAgB,EAAE,CAAS;IACzE,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;QACzB,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC;KAAE;IACnC,OAAO,GAAG,CAAC;AAAC,CAAC;AAchB,MAAM,UAAU,gCAAgC,CAAE,CAAe,EAAE,QAAgB,EAAE,UAAmB;IACrG,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,GAAG,QAAQ,CAAC;IAC9B,IAAI,GAAG,GAAG,qBAAqB,CAAC,CAAC,EAAE,QAAQ,EAAE,CAAC,CAAC,CAAC;IAChD,IAAI,UAAU,IAAI,CAAC,GAAG,CAAC,EAAE;QACtB,GAAG,IAAI,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;KAAE;IACzB,OAAO,GAAG,CAAC;AAAC,CAAC;AAgBhB,MAAM,UAAU,qCAAqC,CAAE,CAAe,EAAE,WAAmB,EAAE,WAAmB,EAAE,iBAA0B;IACzI,IAAI,MAAM,GAAG,CAAC,QAAQ,CAAC;IACvB,IAAI,MAAM,GAAG,WAAW,CAAC;IACzB,KAAK,IAAI,QAAQ,GAAG,WAAW,EAAE,QAAQ,GAAG,WAAW,EAAE,QAAQ,EAAE,EAAE;QAClE,MAAM,GAAG,GAAG,iBAAiB,CAAC,CAAC;YAC5B,qBAAqB,CAAC,CAAC,EAAE,QAAQ,EAAE,CAAC,CAAC,MAAM,GAAG,WAAW,CAAC,CAAC,CAAC;YAC5D,gCAAgC,CAAC,CAAC,EAAE,QAAQ,EAAE,IAAI,CAAC,CAAC;QACvD,IAAI,GAAG,GAAG,MAAM,EAAE;YACf,MAAM,GAAG,GAAG,CAAC;YACb,MAAM,GAAG,QAAQ,CAAC;SAAE;KAAC;IAC3B,OAAO,MAAM,CAAC;AAAC,CAAC;AAgBnB,MAAM,UAAU,0BAA0B,CAAE,CAAe,EAAE,SAAkB;IAC5E,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC9B,KAAK,IAAI,QAAQ,GAAG,CAAC,EAAE,QAAQ,GAAG,CAAC,EAAE,QAAQ,EAAE,EAAE;QAC9C,CAAC,CAAC,QAAQ,CAAC,GAAG,gCAAgC,CAAC,CAAC,EAAE,QAAQ,EAAE,IAAI,CAAC,CAAC;KAAE;IACvE,IAAI,SAAS,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE;QAClC,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACjB,KAAK,IAAI,QAAQ,GAAG,CAAC,EAAE,QAAQ,GAAG,CAAC,EAAE,QAAQ,EAAE,EAAE;YAC9C,CAAC,CAAC,QAAQ,CAAC,IAAI,GAAG,CAAC;SAAE;KAAC;IAC5B,OAAO,CAAC,CAAC;AAAC,CAAC;AAoBd,MAAM,UAAU,kCAAkC,CAAE,CAAe,EAAE,cAA8C,EAAE,SAAkB;IACpI,MAAM,QAAQ,GAAG,eAAe,CAAC,WAAW,CAAC,CAAC,EAAE,cAAc,CAAC,CAAC;IAChE,MAAM,QAAQ,GAAG,0BAA0B,CAAC,QAAQ,EAAE,SAAS,CAAC,CAAC;IACjE,MAAM,MAAM,GAAG,eAAe,CAAC,cAAc,CAAC,cAAc,EAAE,CAAC,CAAC,MAAM,CAAC,CAAC;IACxE,MAAM,cAAc,GAAG,0BAA0B,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC;IAChE,MAAM,WAAW,GAAG,UAAU,CAAC,MAAM,CAAC,QAAQ,EAAE,cAAc,CAAC,CAAC;IAChE,OAAO,WAAW,CAAC;AAAC,CAAC"}
+{"version":3,"file":"Autocorrelation.js","sourceRoot":"","sources":["../../src/signal/Autocorrelation.ts"],"names":[],"mappings":"AAIA,OAAO,KAAK,eAAe,MAAM,sBAAsB,CAAC;AACxD,OAAO,KAAK,UAAU,MAAM,wBAAwB,CAAC;AAErD,SAAS,qBAAqB,CAAE,CAAoB,EAAE,QAAgB,EAAE,CAAS;IAC9E,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;QACzB,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC;KAAE;IACnC,OAAO,GAAG,CAAC;AAAC,CAAC;AAchB,MAAM,UAAU,gCAAgC,CAAE,CAAoB,EAAE,QAAgB,EAAE,UAAmB;IAC1G,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,GAAG,QAAQ,CAAC;IAC9B,IAAI,GAAG,GAAG,qBAAqB,CAAC,CAAC,EAAE,QAAQ,EAAE,CAAC,CAAC,CAAC;IAChD,IAAI,UAAU,IAAI,CAAC,GAAG,CAAC,EAAE;QACtB,GAAG,IAAI,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;KAAE;IACzB,OAAO,GAAG,CAAC;AAAC,CAAC;AAgBhB,MAAM,UAAU,qCAAqC,CAAE,CAAoB,EAAE,WAAmB,EAAE,WAAmB,EAAE,iBAA0B;IAC9I,IAAI,MAAM,GAAG,CAAC,QAAQ,CAAC;IACvB,IAAI,MAAM,GAAG,WAAW,CAAC;IACzB,KAAK,IAAI,QAAQ,GAAG,WAAW,EAAE,QAAQ,GAAG,WAAW,EAAE,QAAQ,EAAE,EAAE;QAClE,MAAM,GAAG,GAAG,iBAAiB,CAAC,CAAC;YAC5B,qBAAqB,CAAC,CAAC,EAAE,QAAQ,EAAE,CAAC,CAAC,MAAM,GAAG,WAAW,CAAC,CAAC,CAAC;YAC5D,gCAAgC,CAAC,CAAC,EAAE,QAAQ,EAAE,IAAI,CAAC,CAAC;QACvD,IAAI,GAAG,GAAG,MAAM,EAAE;YACf,MAAM,GAAG,GAAG,CAAC;YACb,MAAM,GAAG,QAAQ,CAAC;SAAE;KAAC;IAC3B,OAAO,MAAM,CAAC;AAAC,CAAC;AAgBnB,MAAM,UAAU,0BAA0B,CAAE,CAAoB,EAAE,SAAkB;IACjF,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC9B,KAAK,IAAI,QAAQ,GAAG,CAAC,EAAE,QAAQ,GAAG,CAAC,EAAE,QAAQ,EAAE,EAAE;QAC9C,CAAC,CAAC,QAAQ,CAAC,GAAG,gCAAgC,CAAC,CAAC,EAAE,QAAQ,EAAE,IAAI,CAAC,CAAC;KAAE;IACvE,IAAI,SAAS,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE;QAClC,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACjB,KAAK,IAAI,QAAQ,GAAG,CAAC,EAAE,QAAQ,GAAG,CAAC,EAAE,QAAQ,EAAE,EAAE;YAC9C,CAAC,CAAC,QAAQ,CAAC,IAAI,GAAG,CAAC;SAAE;KAAC;IAC5B,OAAO,CAAC,CAAC;AAAC,CAAC;AAoBd,MAAM,UAAU,kCAAkC,CAAE,CAAoB,EAAE,cAA8C,EAAE,SAAkB;IACzI,MAAM,QAAQ,GAAG,eAAe,CAAC,WAAW,CAAC,CAAC,EAAE,cAAc,CAAC,CAAC;IAChE,MAAM,QAAQ,GAAG,0BAA0B,CAAC,QAAQ,EAAE,SAAS,CAAC,CAAC;IACjE,MAAM,MAAM,GAAG,eAAe,CAAC,cAAc,CAAC,cAAc,EAAE,CAAC,CAAC,MAAM,CAAC,CAAC;IACxE,MAAM,cAAc,GAAG,0BAA0B,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC;IAChE,MAAM,WAAW,GAAG,UAAU,CAAC,MAAM,CAAC,QAAQ,EAAE,cAAc,CAAC,CAAC;IAChE,OAAO,WAAW,CAAC;AAAC,CAAC"}

package/signal/Dft.d.ts

@@ -1,9 +1,9 @@
-import MutableComplex from "../math/MutableComplex";
-import ComplexArray from "../math/ComplexArray";
-export declare function dftRealSingle(x: Float64Array, relativeFrequency: number): MutableComplex;
+import MutableComplex from "../math/MutableComplex.js";
+import ComplexArray from "../math/ComplexArray.js";
+export declare function dftRealSingle(x: ArrayLike<number>, relativeFrequency: number): MutableComplex;
 export declare function dftSingle(x: ComplexArray, relativeFrequency: number, direction: boolean): MutableComplex;
-export declare function dftReal(x: Float64Array): ComplexArray;
-export declare function dftRealHalf(x: Float64Array): ComplexArray;
+export declare function dftReal(x: ArrayLike<number>): ComplexArray;
+export declare function dftRealHalf(x: ArrayLike<number>): ComplexArray;
 export declare function dft(x: ComplexArray, direction: boolean): ComplexArray;
-export declare function dftRealSpectrum(x: Float64Array, inclNyquist?: boolean): ComplexArray;
+export declare function dftRealSpectrum(x: ArrayLike<number>, inclNyquist?: boolean): ComplexArray;
 export declare function iDftRealSpectrum(x: ComplexArray, len: number): Float64Array;

package/signal/Dft.js

@@ -1,6 +1,6 @@
-import Complex from "../math/Complex";
-import MutableComplex from "../math/MutableComplex";
-import ComplexArray from "../math/ComplexArray";
+import Complex from "../math/Complex.js";
+import MutableComplex from "../math/MutableComplex.js";
+import ComplexArray from "../math/ComplexArray.js";
 export function dftRealSingle(x, relativeFrequency) {
     const n = x.length;
     if (n == 0) {

package/signal/Dft.js.map

@@ -1 +1 @@
-{"version":3,"file":"Dft.js","sourceRoot":"","sources":["../../src/signal/Dft.ts"],"names":[],"mappings":"AAUA,OAAO,OAAO,MAAM,iBAAiB,CAAC;AACtC,OAAO,cAAc,MAAM,wBAAwB,CAAC;AACpD,OAAO,YAAY,MAAM,sBAAsB,CAAC;AAoBhD,MAAM,UAAU,aAAa,CAAE,CAAe,EAAE,iBAAyB;IACtE,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,IAAI,CAAC,IAAI,CAAC,EAAE;QACT,MAAM,IAAI,KAAK,CAAC,gCAAgC,CAAC,CAAC;KAAE;IACvD,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,iBAAiB,CAAC;IAC/C,MAAM,GAAG,GAAG,IAAI,cAAc,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;QACzB,MAAM,CAAC,GAAG,OAAO,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QACzC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;KAAE;IAClB,OAAO,GAAG,CAAC;AAAC,CAAC;AAiBhB,MAAM,UAAU,SAAS,CAAE,CAAe,EAAE,iBAAyB,EAAE,SAAkB;IACtF,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,IAAI,CAAC,IAAI,CAAC,EAAE;QACT,MAAM,IAAI,KAAK,CAAC,gCAAgC,CAAC,CAAC;KAAE;IACvD,MAAM,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,iBAAiB,CAAC;IACrE,MAAM,GAAG,GAAG,IAAI,cAAc,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;QACzB,MAAM,CAAC,GAAG,cAAc,CAAC,SAAS,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QAC7C,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAClB,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;KAAE;IAClB,OAAO,GAAG,CAAC;AAAC,CAAC;AAWhB,MAAM,UAAU,OAAO,CAAE,CAAe;IACrC,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC9B,KAAK,IAAI,SAAS,GAAG,CAAC,EAAE,SAAS,GAAG,CAAC,EAAE,SAAS,EAAE,EAAE;QACjD,MAAM,CAAC,GAAG,aAAa,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC;QACtC,CAAC,CAAC,GAAG,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;KAAE;IACzB,OAAO,CAAC,CAAC;AAAC,CAAC;AAUd,MAAM,UAAU,WAAW,CAAE,CAAe;IACzC,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IAClC,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC9B,KAAK,IAAI,SAAS,GAAG,CAAC,EAAE,SAAS,GAAG,CAAC,EAAE,SAAS,EAAE,EAAE;QACjD,MAAM,CAAC,GAAG,aAAa,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC;QACtC,CAAC,CAAC,GAAG,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;KAAE;IACzB,OAAO,CAAC,CAAC;AAAC,CAAC;AAYd,MAAM,UAAU,GAAG,CAAE,CAAe,EAAE,SAAkB;IACrD,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC9B,KAAK,IAAI,SAAS,GAAG,CAAC,EAAE,SAAS,GAAG,CAAC,EAAE,SAAS,EAAE,EAAE;QACjD,MAAM,CAAC,GAAG,SAAS,CAAC,CAAC,EAAE,SAAS,EAAE,SAAS,CAAC,CAAC;QAC7C,CAAC,CAAC,GAAG,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;KAAE;IACzB,OAAO,CAAC,CAAC;AAAC,CAAC;AAkBd,MAAM,UAAU,eAAe,CAAE,CAAe,EAAE,WAAW,GAAG,KAAK;IAClE,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,IAAI,CAAC,IAAI,CAAC,EAAE;QACT,MAAM,IAAI,KAAK,CAAC,gCAAgC,CAAC,CAAC;KAAE;IACvD,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACrE,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC9B,KAAK,IAAI,SAAS,GAAG,CAAC,EAAE,SAAS,GAAG,CAAC,EAAE,SAAS,EAAE,EAAE;QACjD,MAAM,CAAC,GAAG,aAAa,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC;QACtC,MAAM,CAAC,GAAG,CAAC,SAAS,IAAI,CAAC,IAAI,SAAS,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACjE,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;QACf,CAAC,CAAC,GAAG,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;KAAE;IACzB,OAAO,CAAC,CAAC;AAAC,CAAC;AAkBd,MAAM,UAAU,gBAAgB,CAAE,CAAe,EAAE,GAAW;IAC3D,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC;IAChC,KAAK,IAAI,SAAS,GAAG,CAAC,EAAE,SAAS,GAAG,CAAC,CAAC,MAAM,EAAE,SAAS,EAAE,EAAE;QACxD,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;QAC3B,oBAAoB,CAAC,CAAC,EAAE,SAAS,EAAE,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;KAAE;IAC1D,OAAO,CAAC,CAAC;AAAC,CAAC;AAEd,SAAS,oBAAoB,CAAE,CAAe,EAAE,SAAiB,EAAE,SAAiB,EAAE,KAAa;IAChG,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,CAAC,GAAG,SAAS,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;IACtC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;QACzB,CAAC,CAAC,CAAC,CAAC,IAAI,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;KAAE;AAAA,CAAC"}
+{"version":3,"file":"Dft.js","sourceRoot":"","sources":["../../src/signal/Dft.ts"],"names":[],"mappings":"AAUA,OAAO,OAAO,MAAM,oBAAoB,CAAC;AACzC,OAAO,cAAc,MAAM,2BAA2B,CAAC;AACvD,OAAO,YAAY,MAAM,yBAAyB,CAAC;AAoBnD,MAAM,UAAU,aAAa,CAAE,CAAoB,EAAE,iBAAyB;IAC3E,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,IAAI,CAAC,IAAI,CAAC,EAAE;QACT,MAAM,IAAI,KAAK,CAAC,gCAAgC,CAAC,CAAC;KAAE;IACvD,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,iBAAiB,CAAC;IAC/C,MAAM,GAAG,GAAG,IAAI,cAAc,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;QACzB,MAAM,CAAC,GAAG,OAAO,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QACzC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;KAAE;IAClB,OAAO,GAAG,CAAC;AAAC,CAAC;AAiBhB,MAAM,UAAU,SAAS,CAAE,CAAe,EAAE,iBAAyB,EAAE,SAAkB;IACtF,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,IAAI,CAAC,IAAI,CAAC,EAAE;QACT,MAAM,IAAI,KAAK,CAAC,gCAAgC,CAAC,CAAC;KAAE;IACvD,MAAM,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,iBAAiB,CAAC;IACrE,MAAM,GAAG,GAAG,IAAI,cAAc,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;QACzB,MAAM,CAAC,GAAG,cAAc,CAAC,SAAS,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QAC7C,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAClB,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;KAAE;IAClB,OAAO,GAAG,CAAC;AAAC,CAAC;AAWhB,MAAM,UAAU,OAAO,CAAE,CAAoB;IAC1C,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC9B,KAAK,IAAI,SAAS,GAAG,CAAC,EAAE,SAAS,GAAG,CAAC,EAAE,SAAS,EAAE,EAAE;QACjD,MAAM,CAAC,GAAG,aAAa,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC;QACtC,CAAC,CAAC,GAAG,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;KAAE;IACzB,OAAO,CAAC,CAAC;AAAC,CAAC;AAUd,MAAM,UAAU,WAAW,CAAE,CAAoB;IAC9C,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IAClC,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC9B,KAAK,IAAI,SAAS,GAAG,CAAC,EAAE,SAAS,GAAG,CAAC,EAAE,SAAS,EAAE,EAAE;QACjD,MAAM,CAAC,GAAG,aAAa,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC;QACtC,CAAC,CAAC,GAAG,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;KAAE;IACzB,OAAO,CAAC,CAAC;AAAC,CAAC;AAYd,MAAM,UAAU,GAAG,CAAE,CAAe,EAAE,SAAkB;IACrD,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC9B,KAAK,IAAI,SAAS,GAAG,CAAC,EAAE,SAAS,GAAG,CAAC,EAAE,SAAS,EAAE,EAAE;QACjD,MAAM,CAAC,GAAG,SAAS,CAAC,CAAC,EAAE,SAAS,EAAE,SAAS,CAAC,CAAC;QAC7C,CAAC,CAAC,GAAG,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;KAAE;IACzB,OAAO,CAAC,CAAC;AAAC,CAAC;AAkBd,MAAM,UAAU,eAAe,CAAE,CAAoB,EAAE,WAAW,GAAG,KAAK;IACvE,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,IAAI,CAAC,IAAI,CAAC,EAAE;QACT,MAAM,IAAI,KAAK,CAAC,gCAAgC,CAAC,CAAC;KAAE;IACvD,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACrE,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC9B,KAAK,IAAI,SAAS,GAAG,CAAC,EAAE,SAAS,GAAG,CAAC,EAAE,SAAS,EAAE,EAAE;QACjD,MAAM,CAAC,GAAG,aAAa,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC;QACtC,MAAM,CAAC,GAAG,CAAC,SAAS,IAAI,CAAC,IAAI,SAAS,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACjE,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;QACf,CAAC,CAAC,GAAG,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;KAAE;IACzB,OAAO,CAAC,CAAC;AAAC,CAAC;AAiBd,MAAM,UAAU,gBAAgB,CAAE,CAAe,EAAE,GAAW;IAC3D,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC;IAChC,KAAK,IAAI,SAAS,GAAG,CAAC,EAAE,SAAS,GAAG,CAAC,CAAC,MAAM,EAAE,SAAS,EAAE,EAAE;QACxD,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;QAC3B,oBAAoB,CAAC,CAAC,EAAE,SAAS,EAAE,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;KAAE;IAC1D,OAAO,CAAC,CAAC;AAAC,CAAC;AAEd,SAAS,oBAAoB,CAAE,CAAe,EAAE,SAAiB,EAAE,SAAiB,EAAE,KAAa;IAChG,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,CAAC,GAAG,SAAS,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;IACtC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;QACzB,CAAC,CAAC,CAAC,CAAC,IAAI,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;KAAE;AAAA,CAAC"}

package/signal/EnvelopeDetection.d.ts

@@ -1 +1 @@
-export declare function generateSignalEnvelope(signal: Float64Array, windowWidthDc: number, windowWidthEnvelope: number): Float64Array;
+export declare function generateSignalEnvelope(signal: ArrayLike<number>, windowWidthDc: number, windowWidthEnvelope: number): Float64Array;

package/signal/EnvelopeDetection.js

@@ -1,4 +1,4 @@
-import * as MathUtils from "../math/MathUtils";
+import * as MathUtils from "../math/MathUtils.js";
 export function generateSignalEnvelope(signal, windowWidthDc, windowWidthEnvelope) {
     const a1 = MathUtils.movingAverage(signal, windowWidthDc);
     for (let i = 0; i < a1.length; i++) {

package/signal/EnvelopeDetection.js.map

@@ -1 +1 @@
-{"version":3,"file":"EnvelopeDetection.js","sourceRoot":"","sources":["../../src/signal/EnvelopeDetection.ts"],"names":[],"mappings":"AAIA,OAAO,KAAK,SAAS,MAAM,mBAAmB,CAAC;AAoB/C,MAAM,UAAU,sBAAsB,CAAE,MAAoB,EAAE,aAAqB,EAAE,mBAA2B;IAC7G,MAAM,EAAE,GAAG,SAAS,CAAC,aAAa,CAAC,MAAM,EAAE,aAAa,CAAC,CAAC;IAC1D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QACjC,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;KAAE;IACzC,MAAM,EAAE,GAAG,SAAS,CAAC,aAAa,CAAC,EAAE,EAAE,mBAAmB,CAAC,CAAC;IAC5D,OAAO,EAAE,CAAC;AAAC,CAAC"}
+{"version":3,"file":"EnvelopeDetection.js","sourceRoot":"","sources":["../../src/signal/EnvelopeDetection.ts"],"names":[],"mappings":"AAIA,OAAO,KAAK,SAAS,MAAM,sBAAsB,CAAC;AAoBlD,MAAM,UAAU,sBAAsB,CAAE,MAAyB,EAAE,aAAqB,EAAE,mBAA2B;IAClH,MAAM,EAAE,GAAG,SAAS,CAAC,aAAa,CAAC,MAAM,EAAE,aAAa,CAAC,CAAC;IAC1D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QACjC,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;KAAE;IACzC,MAAM,EAAE,GAAG,SAAS,CAAC,aAAa,CAAC,EAAE,EAAE,mBAAmB,CAAC,CAAC;IAC5D,OAAO,EAAE,CAAC;AAAC,CAAC"}

package/signal/Fft.d.ts

@@ -1,6 +1,9 @@
-import ComplexArray from "../math/ComplexArray";
+import ComplexArray from "../math/ComplexArray.js";
 export declare function fft(x: ComplexArray, direction?: boolean): ComplexArray;
-export declare function fftReal(x: Float64Array): ComplexArray;
-export declare function fftRealHalf(x: Float64Array): ComplexArray;
-export declare function fftRealSpectrum(x: Float64Array, inclNyquist?: boolean): ComplexArray;
+export declare function fftReal(x: ArrayLike<number>): ComplexArray;
+export declare function fftRealHalf(x: ArrayLike<number>, inclNyquist?: boolean): ComplexArray;
+export declare function fftRealSpectrum(x: ArrayLike<number>, inclNyquist?: boolean): ComplexArray;
 export declare function fftShift(x: ComplexArray): ComplexArray;
+export declare function iFftRealHalfSimple(x: ComplexArray, len: number, inclNyquist?: boolean): Float64Array;
+export declare function iFftRealHalfOpt(x: ComplexArray, len: number, inclNyquist?: boolean): Float64Array;
+export declare function iFftRealHalf(x: ComplexArray, len: number, inclNyquist?: boolean): Float64Array;

package/signal/Fft.js

@@ -1,6 +1,7 @@
-import ComplexArray from "../math/ComplexArray";
-import MutableComplex from "../math/MutableComplex";
-import * as MathUtils from "../math/MathUtils";
+import ComplexArray from "../math/ComplexArray.js";
+import MutableComplex from "../math/MutableComplex.js";
+import * as MathUtils from "../math/MathUtils.js";
+var cooleyTukeySineTableCache;
 export function fft(x, direction = true) {
     const n = x.length;
     if (n <= 1) {
@@ -13,7 +14,7 @@ export function fft(x, direction = true) {
 }
 function fftCooleyTukey(x) {
     const n = x.length;
-    const sineTable = createSineTable(n / 2, n, false);
+    const sineTable = getCachedCooleyTukeySineTable(n);
     const a = copyBitReversed(x);
     applyButterflies(a, sineTable);
     return a;
@@ -95,6 +96,19 @@ function convolve(a1, a2) {
     a5.mulAllByReal(1 / n);
     return a5;
 }
+function getCachedCooleyTukeySineTable(n) {
+    if (!cooleyTukeySineTableCache) {
+        cooleyTukeySineTableCache = new Array(16);
+    }
+    const log2N = MathUtils.floorLog2(n);
+    if (!cooleyTukeySineTableCache[log2N]) {
+        cooleyTukeySineTableCache[log2N] = createCooleyTukeySineTable(n);
+    }
+    return cooleyTukeySineTableCache[log2N];
+}
+function createCooleyTukeySineTable(n) {
+    return createSineTable(n / 2, n, false);
+}
 function createSineTable(tableLength, waveLength, rotationalDirection = true) {
     const w = 2 * Math.PI / waveLength;
     const a = new ComplexArray(tableLength);
@@ -125,7 +139,7 @@ function swapReIm(a) {
 export function fftReal(x) {
     return fft(new ComplexArray(x));
 }
-export function fftRealHalf(x) {
+export function fftRealHalf(x, inclNyquist = false) {
     if (x.length <= 1) {
         return new ComplexArray(x);
     }
@@ -140,12 +154,16 @@ export function fftRealHalf(x) {
         a1.im[i] = x[2 * i + 1];
     }
     const a2 = fft(a1);
-    const a3 = new ComplexArray(n);
-    const w = Math.PI / n;
+    const a3 = new ComplexArray(n + (inclNyquist ? 1 : 0));
     a3.re[0] = a2.re[0] + a2.im[0];
     a3.im[0] = 0;
+    if (inclNyquist) {
+        a3.re[n] = a2.re[0] - a2.im[0];
+        a3.im[n] = 0;
+    }
     const temp1 = new MutableComplex();
     const temp2 = new MutableComplex();
+    const w = Math.PI / n;
     for (let i = 1; i < n; i++) {
         const sRe = Math.sin(i * w);
         const sIm = Math.cos(i * w);
@@ -162,8 +180,8 @@ export function fftRealSpectrum(x, inclNyquist = false) {
         throw new Error("Input array must not be empty.");
     }
     let a;
-    if (n % 2 == 0 && !inclNyquist) {
-        a = fftRealHalf(x);
+    if (n % 2 == 0) {
+        a = fftRealHalf(x, inclNyquist);
     }
     else {
         const a0 = fftReal(x);
@@ -184,4 +202,75 @@ export function fftShift(x) {
     }
     return a;
 }
+export function iFftRealHalfSimple(x, len, inclNyquist = false) {
+    if (x.length == 0 || len <= 0) {
+        return new Float64Array(0);
+    }
+    const x2 = createFullSpectrumFromHalfSpectrum(x, len, inclNyquist);
+    const a = fft(x2, false);
+    return a.re;
+}
+function createFullSpectrumFromHalfSpectrum(x, len, inclNyquist) {
+    const x2 = new ComplexArray(len);
+    ComplexArray.copy1(x, 0, x2, 0);
+    if (inclNyquist && len % 2 == 0 && x.length > len / 2) {
+        ComplexArray.copy1(x, len / 2, x2, len / 2);
+    }
+    const n2 = Math.min(x.length - 1, Math.floor((len - 1) / 2));
+    for (let i = 0; i < n2; i++) {
+        const p1 = 1 + i;
+        const p2 = len - 1 - i;
+        x2.re[p2] = x.re[p1];
+        x2.im[p2] = -x.im[p1];
+    }
+    return x2;
+}
+export function iFftRealHalfOpt(x, len, inclNyquist = false) {
+    if (len <= 0) {
+        return new Float64Array(0);
+    }
+    if (len % 2 != 0) {
+        throw new Error("output length is not even.");
+    }
+    const n = len / 2;
+    const a1 = new ComplexArray(n);
+    a1.re[0] = xRe(0);
+    a1.im[0] = xRe(0);
+    if (inclNyquist) {
+        a1.re[0] += xRe(n);
+        a1.im[0] -= xRe(n);
+    }
+    const temp1 = new MutableComplex();
+    const temp2 = new MutableComplex();
+    const w = Math.PI / n;
+    for (let i = 1; i < n; i++) {
+        const sRe = Math.sin(i * w);
+        const sIm = Math.cos(i * w);
+        temp1.setMul(xRe(i), xIm(i), (1 - sRe) / 2, sIm / 2);
+        temp2.setMul(xRe(n - i), xIm(n - i), (1 + sRe) / 2, sIm / 2);
+        a1.re[i] = temp1.re + temp2.re;
+        a1.im[i] = temp1.im - temp2.im;
+    }
+    const a2 = fft(a1, false);
+    const a3 = new Float64Array(2 * n);
+    for (let i = 0; i < n; i++) {
+        a3[2 * i] = a2.re[i];
+        a3[2 * i + 1] = a2.im[i];
+    }
+    return a3;
+    function xRe(i) {
+        return (i < x.length) ? x.re[i] : 0;
+    }
+    function xIm(i) {
+        return (i < x.length) ? x.im[i] : 0;
+    }
+}
+export function iFftRealHalf(x, len, inclNyquist = false) {
+    if (len % 2 == 0) {
+        return iFftRealHalfOpt(x, len, inclNyquist);
+    }
+    else {
+        return iFftRealHalfSimple(x, len, inclNyquist);
+    }
+}
 //# sourceMappingURL=Fft.js.map