dsp-collection 0.2.5 → 0.2.7

package/README.md

@@ -1,9 +1,11 @@
-# The source-code.biz JavaScript DSP Collection
+# The JavaScript DSP Collection
 
 A collection of JavaScript modules for digital signal processing (written in TypeScript).
 
 (This package is currently in development)
 
+NPM package: [dsp-collection](https://www.npmjs.com/package/dsp-collection)
+
 ## Online test applications
 
 * [Window functions test](https://www.source-code.biz/dsp/js/simpleTests/testWindowFunctions.html)
@@ -13,12 +15,14 @@ A collection of JavaScript modules for digital signal processing (written in Typ
 
 * [HarmSyn - Harmonic synthesizer](https://www.source-code.biz/harmSyn)
 * [KlattSyn - Klatt formant synthesizer](https://www.source-code.biz/klattSyn)
+* [NoiseSyn - Spectral noise synthesizer](https://www.source-code.biz/noiseSyn)
 * [SinSyn - Sinusoidal synthesizer](https://www.source-code.biz/sinSyn)
 * [SpecFilt - Spectral filter tool](https://www.source-code.biz/specFilt)
 
 ## Associated JavaScript packages
 
-* [commons-math-interpolation - Interpolation and regression algorithms](https://github.com/chdh/commons-math-interpolation)
+* [commons-math-interpolation](https://github.com/chdh/commons-math-interpolation) - Interpolation and regression algorithms
 * [Function curve editor widget](https://www.source-code.biz/snippets/typescript/functionCurveEditor)
 * [Function curve viewer widget](https://www.source-code.biz/snippets/typescript/functionCurveViewer)
 * [WAV file encoder](https://github.com/chdh/wav-file-encoder)
+* [WAV file decoder](https://github.com/chdh/wav-file-decoder)

package/filter/FirFilterWin.d.ts

@@ -0,0 +1,8 @@
+import MutableComplex from "../math/MutableComplex.ts";
+import { WindowFunction } from "../signal/WindowFunctions.ts";
+export declare function createWindowKernel(f: WindowFunction, width: number, symetric?: boolean): Float64Array;
+export declare function createLpFilterKernel(windowFunctionId: string, normFirstMinFreq: number): Float64Array;
+export declare function calcFreqRespAt(kernel: ArrayLike<number>, normFreq: number): MutableComplex;
+export declare function applyFirKernelAt(signal: ArrayLike<number>, pos: number, kernel: ArrayLike<number>): number;
+export declare function applyFirKernel(signal: ArrayLike<number>, kernel: Float64Array): Float64Array;
+//# sourceMappingURL=FirFilterWin.d.ts.map

package/filter/FirFilterWin.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"FirFilterWin.d.ts","sourceRoot":"","sources":["../../src/filter/FirFilterWin.ts"],"names":[],"mappings":"AAEA,OAAO,cAAc,MAAM,2BAA2B,CAAC;AAEvD,OAAO,EAAC,cAAc,EAAC,MAAM,8BAA8B,CAAC;AAM5D,wBAAgB,kBAAkB,CAAE,CAAC,EAAE,cAAc,EAAE,KAAK,EAAE,MAAM,EAAE,QAAQ,GAAE,OAAc,GAAI,YAAY,CAK/F;AAKf,wBAAgB,oBAAoB,CAAE,gBAAgB,EAAE,MAAM,EAAE,gBAAgB,EAAE,MAAM,GAAI,YAAY,CAKzD;AAI/C,wBAAgB,cAAc,CAAE,MAAM,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,QAAQ,EAAE,MAAM,GAAG,cAAc,CAEhD;AAY5C,wBAAgB,gBAAgB,CAAE,MAAM,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,GAAG,EAAE,MAAM,EAAE,MAAM,EAAE,SAAS,CAAC,MAAM,CAAC,GAAI,MAAM,CAU7F;AAGhB,wBAAgB,cAAc,CAAE,MAAM,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,MAAM,EAAE,YAAY,GAAI,YAAY,CAI/E"}

package/filter/FirFilterWin.js

@@ -0,0 +1,55 @@
+import * as WindowFunctions from "../signal/WindowFunctions.js";
+import { goertzelSingle } from "../signal/Goertzel.js";
+import * as ArrayUtils from "../utils/ArrayUtils.js";
+export function createWindowKernel(f, width, symetric = true) {
+    const nudge = symetric ? 1 : 0;
+    const a1 = Float64Array.from({ length: width }, (_x, i) => f(i / (width - nudge)));
+    const sum = ArrayUtils.sum(a1);
+    const a2 = a1.map(x => x / sum);
+    return a2;
+}
+export function createLpFilterKernel(windowFunctionId, normFirstMinFreq) {
+    const descr = WindowFunctions.getFunctionDescrById(windowFunctionId);
+    const width = Math.round(descr.firstMinPos / normFirstMinFreq);
+    if (width < 3) {
+        throw new Error("Filter parameters out of range.");
+    }
+    return createWindowKernel(descr.f, width);
+}
+export function calcFreqRespAt(kernel, normFreq) {
+    const relFreq = normFreq * kernel.length;
+    return goertzelSingle(kernel, relFreq);
+}
+function reflectIndex(i, n) {
+    if (i >= 0 && i < n) {
+        return i;
+    }
+    if (n <= 1) {
+        return 0;
+    }
+    const period = 2 * (n - 1);
+    const t = (i % period + period) % period;
+    return (t < n) ? t : period - t;
+}
+export function applyFirKernelAt(signal, pos, kernel) {
+    if (pos < 0 || pos >= signal.length) {
+        return NaN;
+    }
+    const width = kernel.length;
+    const m = Math.floor(kernel.length / 2);
+    const p0 = pos - m;
+    let acc = 0;
+    for (let i = 0; i < width; i++) {
+        const p = reflectIndex(p0 + i, signal.length);
+        acc += signal[p] * kernel[i];
+    }
+    return acc;
+}
+export function applyFirKernel(signal, kernel) {
+    const out = new Float64Array(signal.length);
+    for (let i = 0; i < signal.length; i++) {
+        out[i] = applyFirKernelAt(signal, i, kernel);
+    }
+    return out;
+}
+//# sourceMappingURL=FirFilterWin.js.map

package/filter/FirFilterWin.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"FirFilterWin.js","sourceRoot":"","sources":["../../src/filter/FirFilterWin.ts"],"names":[],"mappings":"AAGA,OAAO,KAAK,eAAe,MAAM,8BAA8B,CAAC;AAEhE,OAAO,EAAC,cAAc,EAAC,MAAM,uBAAuB,CAAC;AACrD,OAAO,KAAK,UAAU,MAAM,wBAAwB,CAAC;AAIrD,MAAM,UAAU,kBAAkB,CAAE,CAAiB,EAAE,KAAa,EAAE,WAAoB,IAAI;IAC3F,MAAM,KAAK,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC/B,MAAM,EAAE,GAAG,YAAY,CAAC,IAAI,CAAC,EAAC,MAAM,EAAE,KAAK,EAAC,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IACjF,MAAM,GAAG,GAAG,UAAU,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;IAC/B,MAAM,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;IAChC,OAAO,EAAE,CAAC;AAAC,CAAC;AAKf,MAAM,UAAU,oBAAoB,CAAE,gBAAwB,EAAE,gBAAwB;IACrF,MAAM,KAAK,GAAG,eAAe,CAAC,oBAAoB,CAAC,gBAAgB,CAAC,CAAC;IACrE,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,WAAW,GAAG,gBAAgB,CAAC,CAAC;IAC/D,IAAI,KAAK,GAAG,CAAC,EAAE,CAAC;QACb,MAAM,IAAI,KAAK,CAAC,iCAAiC,CAAC,CAAC;IAAC,CAAC;IACxD,OAAO,kBAAkB,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;AAAC,CAAC;AAI/C,MAAM,UAAU,cAAc,CAAE,MAAyB,EAAE,QAAgB;IACxE,MAAM,OAAO,GAAG,QAAQ,GAAG,MAAM,CAAC,MAAM,CAAC;IACzC,OAAO,cAAc,CAAC,MAAM,EAAE,OAAO,CAAC,CAAC;AAAC,CAAC;AAE5C,SAAS,YAAY,CAAE,CAAS,EAAE,CAAS;IACxC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QACnB,OAAO,CAAC,CAAC;IAAC,CAAC;IACd,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QACV,OAAO,CAAC,CAAC;IAAC,CAAC;IACd,MAAM,MAAM,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC3B,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,GAAG,MAAM,CAAC,GAAG,MAAM,CAAC;IACzC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;AAAC,CAAC;AAGrC,MAAM,UAAU,gBAAgB,CAAE,MAAyB,EAAE,GAAW,EAAE,MAAyB;IAChG,IAAI,GAAG,GAAG,CAAC,IAAI,GAAG,IAAI,MAAM,CAAC,MAAM,EAAE,CAAC;QACnC,OAAO,GAAG,CAAC;IAAC,CAAC;IAChB,MAAM,KAAK,GAAG,MAAM,CAAC,MAAM,CAAC;IAC5B,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IACxC,MAAM,EAAE,GAAG,GAAG,GAAG,CAAC,CAAC;IACnB,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC9B,MAAM,CAAC,GAAG,YAAY,CAAC,EAAE,GAAG,CAAC,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC;QAC9C,GAAG,IAAI,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;IAAC,CAAC;IAClC,OAAO,GAAG,CAAC;AAAC,CAAC;AAGhB,MAAM,UAAU,cAAc,CAAE,MAAyB,EAAE,MAAoB;IAC5E,MAAM,GAAG,GAAG,IAAI,YAAY,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC;IAC5C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACtC,GAAG,CAAC,CAAC,CAAC,GAAG,gBAAgB,CAAC,MAAM,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;IAAC,CAAC;IAClD,OAAO,GAAG,CAAC;AAAC,CAAC"}

package/filter/SpecFilt.d.ts

@@ -1,10 +1,11 @@
-export declare const enum FilterType {
-    lowPass = "LP",
-    highPass = "HP",
-    bandPass = "BP",
-    bandStop = "BS"
-}
-export declare type FilterCurveFunction = (frequency: number) => number;
-export declare function getFilterCurveFunction(filterType: FilterType, filterFreq1: number, filterFreq2: number, smoothingWidth: number): FilterCurveFunction;
-export declare function applyFilterCurveFunction(inAmplitudes: ArrayLike<number>, scalingFactor: number, filterCurveFunction: FilterCurveFunction): Float64Array;
-export declare function filterSignal(inSamples: ArrayLike<number>, sampleRate: number, filterType: FilterType, filterFreq1: number, filterFreq2: number, smoothingWidth: number): Float64Array;
+export declare const enum FilterType {
+    lowPass = "LP",
+    highPass = "HP",
+    bandPass = "BP",
+    bandStop = "BS"
+}
+export type FilterCurveFunction = (frequency: number) => number;
+export declare function getFilterCurveFunction(filterType: FilterType, filterFreq1: number, filterFreq2: number, smoothingWidth: number): FilterCurveFunction;
+export declare function applyFilterCurveFunction(inAmplitudes: ArrayLike<number>, scalingFactor: number, filterCurveFunction: FilterCurveFunction): Float64Array<ArrayBuffer>;
+export declare function filterSignal(inSamples: ArrayLike<number>, sampleRate: number, filterType: FilterType, filterFreq1: number, filterFreq2: number, smoothingWidth: number): Float64Array;
+//# sourceMappingURL=SpecFilt.d.ts.map

package/filter/SpecFilt.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"SpecFilt.d.ts","sourceRoot":"","sources":["../../src/filter/SpecFilt.ts"],"names":[],"mappings":"AASA,0BAAkB,UAAU;IACzB,OAAO,OAAQ;IACf,QAAQ,OAAO;IACf,QAAQ,OAAO;IACf,QAAQ,OAAO;CAAE;AAEpB,MAAM,MAAM,mBAAmB,GAAG,CAAC,SAAS,EAAE,MAAM,KAAK,MAAM,CAAC;AAoBhE,wBAAgB,sBAAsB,CAAE,UAAU,EAAE,UAAU,EAAE,WAAW,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,EAAE,cAAc,EAAE,MAAM,GAAI,mBAAmB,CA2C9F;AAKzD,wBAAgB,wBAAwB,CAAE,YAAY,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,aAAa,EAAE,MAAM,EAAE,mBAAmB,EAAE,mBAAmB,6BAOhH;AAwB1B,wBAAgB,YAAY,CAAE,SAAS,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,UAAU,EAAE,MAAM,EAAE,UAAU,EAAE,UAAU,EAAE,WAAW,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,EAAE,cAAc,EAAE,MAAM,GAAI,YAAY,CAUlK"}

package/filter/SpecFilt.js

@@ -1,99 +1,99 @@
-import ComplexArray from "../math/ComplexArray.js";
-import * as Fft from "../signal/Fft.js";
-function smoothingFunction(x) {
-    return (Math.sin(x * Math.PI / 2) + 1) / 2;
-}
-export function getFilterCurveFunction(filterType, filterFreq1, filterFreq2, smoothingWidth) {
-    switch (filterType) {
-        case "LP": {
-            return (freq) => {
-                if (freq < filterFreq1 - smoothingWidth) {
-                    return 1;
-                }
-                else if (freq < filterFreq1 + smoothingWidth) {
-                    return smoothingFunction((filterFreq1 - freq) / smoothingWidth);
-                }
-                else {
-                    return 0;
-                }
-            };
-        }
-        case "HP": {
-            return (freq) => {
-                if (freq < filterFreq1 - smoothingWidth) {
-                    return 0;
-                }
-                else if (freq < filterFreq1 + smoothingWidth) {
-                    return smoothingFunction((freq - filterFreq1) / smoothingWidth);
-                }
-                else {
-                    return 1;
-                }
-            };
-        }
-        case "BP": {
-            return (freq) => {
-                if (freq < filterFreq1 - smoothingWidth) {
-                    return 0;
-                }
-                else if (freq < filterFreq1 + smoothingWidth) {
-                    return smoothingFunction((freq - filterFreq1) / smoothingWidth);
-                }
-                else if (freq < filterFreq2 - smoothingWidth) {
-                    return 1;
-                }
-                else if (freq < filterFreq2 + smoothingWidth) {
-                    return smoothingFunction((filterFreq2 - freq) / smoothingWidth);
-                }
-                else {
-                    return 0;
-                }
-            };
-        }
-        case "BS": {
-            return (freq) => {
-                if (freq < filterFreq1 - smoothingWidth) {
-                    return 1;
-                }
-                else if (freq < filterFreq1 + smoothingWidth) {
-                    return smoothingFunction((filterFreq1 - freq) / smoothingWidth);
-                }
-                else if (freq < filterFreq2 - smoothingWidth) {
-                    return 0;
-                }
-                else if (freq < filterFreq2 + smoothingWidth) {
-                    return smoothingFunction((freq - filterFreq2) / smoothingWidth);
-                }
-                else {
-                    return 1;
-                }
-            };
-        }
-        default: {
-            throw new Error("Unsupported filter type.");
-        }
-    }
-}
-export function applyFilterCurveFunction(inAmplitudes, scalingFactor, filterCurveFunction) {
-    const n = inAmplitudes.length;
-    const outAmplitudes = new Float64Array(n);
-    for (let p = 0; p < n; p++) {
-        const frequency = p / scalingFactor;
-        const filterFactor = filterCurveFunction(frequency);
-        outAmplitudes[p] = filterFactor * inAmplitudes[p];
-    }
-    return outAmplitudes;
-}
-export function filterSignal(inSamples, sampleRate, filterType, filterFreq1, filterFreq2, smoothingWidth) {
-    const n = inSamples.length;
-    const inSpectrum = Fft.fftRealSpectrum(inSamples);
-    const inAmplitudes = inSpectrum.getAbsArray();
-    const inPhases = inSpectrum.getArgArray();
-    const filterCurveFunction = getFilterCurveFunction(filterType, filterFreq1, filterFreq2, smoothingWidth);
-    const outAmplitudes = applyFilterCurveFunction(inAmplitudes, n / sampleRate, filterCurveFunction);
-    const outPhases = inPhases;
-    const outSpectrum = ComplexArray.fromPolar(outAmplitudes, outPhases);
-    const outSignal = Fft.iFftRealHalf(outSpectrum, n);
-    return outSignal;
-}
+import ComplexArray from "../math/ComplexArray.js";
+import * as Fft from "../signal/Fft.js";
+function smoothingFunction(x) {
+    return (Math.sin(x * Math.PI / 2) + 1) / 2;
+}
+export function getFilterCurveFunction(filterType, filterFreq1, filterFreq2, smoothingWidth) {
+    switch (filterType) {
+        case "LP": {
+            return (freq) => {
+                if (freq < filterFreq1 - smoothingWidth) {
+                    return 1;
+                }
+                else if (freq < filterFreq1 + smoothingWidth) {
+                    return smoothingFunction((filterFreq1 - freq) / smoothingWidth);
+                }
+                else {
+                    return 0;
+                }
+            };
+        }
+        case "HP": {
+            return (freq) => {
+                if (freq < filterFreq1 - smoothingWidth) {
+                    return 0;
+                }
+                else if (freq < filterFreq1 + smoothingWidth) {
+                    return smoothingFunction((freq - filterFreq1) / smoothingWidth);
+                }
+                else {
+                    return 1;
+                }
+            };
+        }
+        case "BP": {
+            return (freq) => {
+                if (freq < filterFreq1 - smoothingWidth) {
+                    return 0;
+                }
+                else if (freq < filterFreq1 + smoothingWidth) {
+                    return smoothingFunction((freq - filterFreq1) / smoothingWidth);
+                }
+                else if (freq < filterFreq2 - smoothingWidth) {
+                    return 1;
+                }
+                else if (freq < filterFreq2 + smoothingWidth) {
+                    return smoothingFunction((filterFreq2 - freq) / smoothingWidth);
+                }
+                else {
+                    return 0;
+                }
+            };
+        }
+        case "BS": {
+            return (freq) => {
+                if (freq < filterFreq1 - smoothingWidth) {
+                    return 1;
+                }
+                else if (freq < filterFreq1 + smoothingWidth) {
+                    return smoothingFunction((filterFreq1 - freq) / smoothingWidth);
+                }
+                else if (freq < filterFreq2 - smoothingWidth) {
+                    return 0;
+                }
+                else if (freq < filterFreq2 + smoothingWidth) {
+                    return smoothingFunction((freq - filterFreq2) / smoothingWidth);
+                }
+                else {
+                    return 1;
+                }
+            };
+        }
+        default: {
+            throw new Error("Unsupported filter type.");
+        }
+    }
+}
+export function applyFilterCurveFunction(inAmplitudes, scalingFactor, filterCurveFunction) {
+    const n = inAmplitudes.length;
+    const outAmplitudes = new Float64Array(n);
+    for (let p = 0; p < n; p++) {
+        const frequency = p / scalingFactor;
+        const filterFactor = filterCurveFunction(frequency);
+        outAmplitudes[p] = filterFactor * inAmplitudes[p];
+    }
+    return outAmplitudes;
+}
+export function filterSignal(inSamples, sampleRate, filterType, filterFreq1, filterFreq2, smoothingWidth) {
+    const n = inSamples.length;
+    const inSpectrum = Fft.fftRealSpectrum(inSamples);
+    const inAmplitudes = inSpectrum.getAbsArray();
+    const inPhases = inSpectrum.getArgArray();
+    const filterCurveFunction = getFilterCurveFunction(filterType, filterFreq1, filterFreq2, smoothingWidth);
+    const outAmplitudes = applyFilterCurveFunction(inAmplitudes, n / sampleRate, filterCurveFunction);
+    const outPhases = inPhases;
+    const outSpectrum = ComplexArray.fromPolar(outAmplitudes, outPhases);
+    const outSignal = Fft.iFftRealHalf(outSpectrum, n);
+    return outSignal;
+}
 //# sourceMappingURL=SpecFilt.js.map

package/filter/SpecFilt.js.map

@@ -1 +1 @@
-{"version":3,"file":"SpecFilt.js","sourceRoot":"","sources":["../../src/filter/SpecFilt.ts"],"names":[],"mappings":"AAMA,OAAO,YAAY,MAAM,yBAAyB,CAAC;AACnD,OAAO,KAAK,GAAG,MAAM,kBAAkB,CAAC;AAWxC,SAAS,iBAAiB,CAAE,CAAS;IAClC,OAAO,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;AAAC,CAAC;AAgBhD,MAAM,UAAU,sBAAsB,CAAE,UAAsB,EAAE,WAAmB,EAAE,WAAmB,EAAE,cAAsB;IAC7H,QAAQ,UAAU,EAAE;QACjB,SAAuB,CAAC,CAAC;YACtB,OAAO,CAAC,IAAY,EAAE,EAAE;gBACrB,IAAI,IAAI,GAAG,WAAW,GAAG,cAAc,EAAE;oBACtC,OAAO,CAAC,CAAC;iBAAE;qBACR,IAAI,IAAI,GAAG,WAAW,GAAG,cAAc,EAAE;oBAC5C,OAAO,iBAAiB,CAAC,CAAC,WAAW,GAAG,IAAI,CAAC,GAAG,cAAc,CAAC,CAAC;iBAAE;qBAC/D;oBACH,OAAO,CAAC,CAAC;iBAAE;YAAA,CAAC,CAAC;SAAE;QACxB,SAAwB,CAAC,CAAC;YACvB,OAAO,CAAC,IAAY,EAAE,EAAE;gBACrB,IAAI,IAAI,GAAG,WAAW,GAAG,cAAc,EAAE;oBACtC,OAAO,CAAC,CAAC;iBAAE;qBACR,IAAI,IAAI,GAAG,WAAW,GAAG,cAAc,EAAE;oBAC5C,OAAO,iBAAiB,CAAC,CAAC,IAAI,GAAG,WAAW,CAAC,GAAG,cAAc,CAAC,CAAC;iBAAE;qBAC/D;oBACH,OAAO,CAAC,CAAC;iBAAE;YAAA,CAAC,CAAC;SAAE;QACxB,SAAwB,CAAC,CAAC;YACvB,OAAO,CAAC,IAAY,EAAE,EAAE;gBACrB,IAAI,IAAI,GAAG,WAAW,GAAG,cAAc,EAAE;oBACtC,OAAO,CAAC,CAAC;iBAAE;qBACR,IAAI,IAAI,GAAG,WAAW,GAAG,cAAc,EAAE;oBAC5C,OAAO,iBAAiB,CAAC,CAAC,IAAI,GAAG,WAAW,CAAC,GAAG,cAAc,CAAC,CAAC;iBAAE;qBAC/D,IAAI,IAAI,GAAG,WAAW,GAAG,cAAc,EAAE;oBAC5C,OAAO,CAAC,CAAC;iBAAE;qBACR,IAAI,IAAI,GAAG,WAAW,GAAG,cAAc,EAAE;oBAC5C,OAAO,iBAAiB,CAAC,CAAC,WAAW,GAAG,IAAI,CAAC,GAAG,cAAc,CAAC,CAAC;iBAAE;qBAC/D;oBACH,OAAO,CAAC,CAAC;iBAAE;YAAA,CAAC,CAAC;SAAE;QACxB,SAAwB,CAAC,CAAC;YACvB,OAAO,CAAC,IAAY,EAAE,EAAE;gBACrB,IAAI,IAAI,GAAG,WAAW,GAAG,cAAc,EAAE;oBACtC,OAAO,CAAC,CAAC;iBAAE;qBACR,IAAI,IAAI,GAAG,WAAW,GAAG,cAAc,EAAE;oBAC5C,OAAO,iBAAiB,CAAC,CAAC,WAAW,GAAG,IAAI,CAAC,GAAG,cAAc,CAAC,CAAC;iBAAE;qBAC/D,IAAI,IAAI,GAAG,WAAW,GAAG,cAAc,EAAE;oBAC5C,OAAO,CAAC,CAAC;iBAAE;qBACR,IAAI,IAAI,GAAG,WAAW,GAAG,cAAc,EAAE;oBAC5C,OAAO,iBAAiB,CAAC,CAAC,IAAI,GAAG,WAAW,CAAC,GAAG,cAAc,CAAC,CAAC;iBAAE;qBAC/D;oBACH,OAAO,CAAC,CAAC;iBAAE;YAAA,CAAC,CAAC;SAAE;QACxB,OAAO,CAAC,CAAC;YACN,MAAM,IAAI,KAAK,CAAC,0BAA0B,CAAC,CAAC;SAAE;KAAC;AAAA,CAAC;AAKzD,MAAM,UAAU,wBAAwB,CAAE,YAA+B,EAAE,aAAqB,EAAE,mBAAwC;IACvI,MAAM,CAAC,GAAG,YAAY,CAAC,MAAM,CAAC;IAC9B,MAAM,aAAa,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC1C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;QACzB,MAAM,SAAS,GAAG,CAAC,GAAG,aAAa,CAAC;QACpC,MAAM,YAAY,GAAG,mBAAmB,CAAC,SAAS,CAAC,CAAC;QACpD,aAAa,CAAC,CAAC,CAAC,GAAG,YAAY,GAAG,YAAY,CAAC,CAAC,CAAC,CAAC;KAAE;IACvD,OAAO,aAAa,CAAC;AAAC,CAAC;AAwB1B,MAAM,UAAU,YAAY,CAAE,SAA4B,EAAE,UAAkB,EAAE,UAAsB,EAAE,WAAmB,EAAE,WAAmB,EAAE,cAAsB;IACrK,MAAM,CAAC,GAAG,SAAS,CAAC,MAAM,CAAC;IAC3B,MAAM,UAAU,GAAG,GAAG,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC;IAClD,MAAM,YAAY,GAAG,UAAU,CAAC,WAAW,EAAE,CAAC;IAC9C,MAAM,QAAQ,GAAG,UAAU,CAAC,WAAW,EAAE,CAAC;IAC1C,MAAM,mBAAmB,GAAG,sBAAsB,CAAC,UAAU,EAAE,WAAW,EAAE,WAAW,EAAE,cAAc,CAAC,CAAC;IACzG,MAAM,aAAa,GAAG,wBAAwB,CAAC,YAAY,EAAE,CAAC,GAAG,UAAU,EAAE,mBAAmB,CAAC,CAAC;IAClG,MAAM,SAAS,GAAG,QAAQ,CAAC;IAC3B,MAAM,WAAW,GAAG,YAAY,CAAC,SAAS,CAAC,aAAa,EAAE,SAAS,CAAC,CAAC;IACrE,MAAM,SAAS,GAAG,GAAG,CAAC,YAAY,CAAC,WAAW,EAAE,CAAC,CAAC,CAAC;IACnD,OAAO,SAAS,CAAC;AAAC,CAAC"}
+{"version":3,"file":"SpecFilt.js","sourceRoot":"","sources":["../../src/filter/SpecFilt.ts"],"names":[],"mappings":"AAMA,OAAO,YAAY,MAAM,yBAAyB,CAAC;AACnD,OAAO,KAAK,GAAG,MAAM,kBAAkB,CAAC;AAWxC,SAAS,iBAAiB,CAAE,CAAS;IAClC,OAAO,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;AAAC,CAAC;AAgBhD,MAAM,UAAU,sBAAsB,CAAE,UAAsB,EAAE,WAAmB,EAAE,WAAmB,EAAE,cAAsB;IAC7H,QAAQ,UAAU,EAAE,CAAC;QAClB,SAAuB,CAAC,CAAC,CAAC;YACvB,OAAO,CAAC,IAAY,EAAE,EAAE;gBACrB,IAAI,IAAI,GAAG,WAAW,GAAG,cAAc,EAAE,CAAC;oBACvC,OAAO,CAAC,CAAC;gBAAC,CAAC;qBACR,IAAI,IAAI,GAAG,WAAW,GAAG,cAAc,EAAE,CAAC;oBAC7C,OAAO,iBAAiB,CAAC,CAAC,WAAW,GAAG,IAAI,CAAC,GAAG,cAAc,CAAC,CAAC;gBAAC,CAAC;qBAC/D,CAAC;oBACJ,OAAO,CAAC,CAAC;gBAAC,CAAC;YAAA,CAAC,CAAC;QAAC,CAAC;QACxB,SAAwB,CAAC,CAAC,CAAC;YACxB,OAAO,CAAC,IAAY,EAAE,EAAE;gBACrB,IAAI,IAAI,GAAG,WAAW,GAAG,cAAc,EAAE,CAAC;oBACvC,OAAO,CAAC,CAAC;gBAAC,CAAC;qBACR,IAAI,IAAI,GAAG,WAAW,GAAG,cAAc,EAAE,CAAC;oBAC7C,OAAO,iBAAiB,CAAC,CAAC,IAAI,GAAG,WAAW,CAAC,GAAG,cAAc,CAAC,CAAC;gBAAC,CAAC;qBAC/D,CAAC;oBACJ,OAAO,CAAC,CAAC;gBAAC,CAAC;YAAA,CAAC,CAAC;QAAC,CAAC;QACxB,SAAwB,CAAC,CAAC,CAAC;YACxB,OAAO,CAAC,IAAY,EAAE,EAAE;gBACrB,IAAI,IAAI,GAAG,WAAW,GAAG,cAAc,EAAE,CAAC;oBACvC,OAAO,CAAC,CAAC;gBAAC,CAAC;qBACR,IAAI,IAAI,GAAG,WAAW,GAAG,cAAc,EAAE,CAAC;oBAC7C,OAAO,iBAAiB,CAAC,CAAC,IAAI,GAAG,WAAW,CAAC,GAAG,cAAc,CAAC,CAAC;gBAAC,CAAC;qBAC/D,IAAI,IAAI,GAAG,WAAW,GAAG,cAAc,EAAE,CAAC;oBAC7C,OAAO,CAAC,CAAC;gBAAC,CAAC;qBACR,IAAI,IAAI,GAAG,WAAW,GAAG,cAAc,EAAE,CAAC;oBAC7C,OAAO,iBAAiB,CAAC,CAAC,WAAW,GAAG,IAAI,CAAC,GAAG,cAAc,CAAC,CAAC;gBAAC,CAAC;qBAC/D,CAAC;oBACJ,OAAO,CAAC,CAAC;gBAAC,CAAC;YAAA,CAAC,CAAC;QAAC,CAAC;QACxB,SAAwB,CAAC,CAAC,CAAC;YACxB,OAAO,CAAC,IAAY,EAAE,EAAE;gBACrB,IAAI,IAAI,GAAG,WAAW,GAAG,cAAc,EAAE,CAAC;oBACvC,OAAO,CAAC,CAAC;gBAAC,CAAC;qBACR,IAAI,IAAI,GAAG,WAAW,GAAG,cAAc,EAAE,CAAC;oBAC7C,OAAO,iBAAiB,CAAC,CAAC,WAAW,GAAG,IAAI,CAAC,GAAG,cAAc,CAAC,CAAC;gBAAC,CAAC;qBAC/D,IAAI,IAAI,GAAG,WAAW,GAAG,cAAc,EAAE,CAAC;oBAC7C,OAAO,CAAC,CAAC;gBAAC,CAAC;qBACR,IAAI,IAAI,GAAG,WAAW,GAAG,cAAc,EAAE,CAAC;oBAC7C,OAAO,iBAAiB,CAAC,CAAC,IAAI,GAAG,WAAW,CAAC,GAAG,cAAc,CAAC,CAAC;gBAAC,CAAC;qBAC/D,CAAC;oBACJ,OAAO,CAAC,CAAC;gBAAC,CAAC;YAAA,CAAC,CAAC;QAAC,CAAC;QACxB,OAAO,CAAC,CAAC,CAAC;YACP,MAAM,IAAI,KAAK,CAAC,0BAA0B,CAAC,CAAC;QAAC,CAAC;IAAA,CAAC;AAAA,CAAC;AAKzD,MAAM,UAAU,wBAAwB,CAAE,YAA+B,EAAE,aAAqB,EAAE,mBAAwC;IACvI,MAAM,CAAC,GAAG,YAAY,CAAC,MAAM,CAAC;IAC9B,MAAM,aAAa,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC1C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC1B,MAAM,SAAS,GAAG,CAAC,GAAG,aAAa,CAAC;QACpC,MAAM,YAAY,GAAG,mBAAmB,CAAC,SAAS,CAAC,CAAC;QACpD,aAAa,CAAC,CAAC,CAAC,GAAG,YAAY,GAAG,YAAY,CAAC,CAAC,CAAC,CAAC;IAAC,CAAC;IACvD,OAAO,aAAa,CAAC;AAAC,CAAC;AAwB1B,MAAM,UAAU,YAAY,CAAE,SAA4B,EAAE,UAAkB,EAAE,UAAsB,EAAE,WAAmB,EAAE,WAAmB,EAAE,cAAsB;IACrK,MAAM,CAAC,GAAG,SAAS,CAAC,MAAM,CAAC;IAC3B,MAAM,UAAU,GAAG,GAAG,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC;IAClD,MAAM,YAAY,GAAG,UAAU,CAAC,WAAW,EAAE,CAAC;IAC9C,MAAM,QAAQ,GAAG,UAAU,CAAC,WAAW,EAAE,CAAC;IAC1C,MAAM,mBAAmB,GAAG,sBAAsB,CAAC,UAAU,EAAE,WAAW,EAAE,WAAW,EAAE,cAAc,CAAC,CAAC;IACzG,MAAM,aAAa,GAAG,wBAAwB,CAAC,YAAY,EAAE,CAAC,GAAG,UAAU,EAAE,mBAAmB,CAAC,CAAC;IAClG,MAAM,SAAS,GAAG,QAAQ,CAAC;IAC3B,MAAM,WAAW,GAAG,YAAY,CAAC,SAAS,CAAC,aAAa,EAAE,SAAS,CAAC,CAAC;IACrE,MAAM,SAAS,GAAG,GAAG,CAAC,YAAY,CAAC,WAAW,EAAE,CAAC,CAAC,CAAC;IACnD,OAAO,SAAS,CAAC;AAAC,CAAC"}

package/math/Complex.d.ts

@@ -1,42 +1,43 @@
-export default class Complex {
-    static readonly I: Complex;
-    static readonly ZERO: Complex;
-    static readonly ONE: Complex;
-    static readonly TWO: Complex;
-    static readonly NaN: Complex;
-    static readonly INFINITY: Complex;
-    readonly re: number;
-    readonly im: number;
-    constructor(re: number, im?: number);
-    toString(): string;
-    toNumber(eps: number): number;
-    isNaN(): boolean;
-    isInfinite(): boolean;
-    isFinite(): boolean;
-    equals(x: Complex): boolean;
-    fuzzyEquals(x: Complex, eps: number): boolean;
-    static expj(arg: number): Complex;
-    static fromPolar(abs: number, arg: number): Complex;
-    abs(): number;
-    arg(): number;
-    conj(): Complex;
-    neg(): Complex;
-    reciprocal(): Complex;
-    exp(): Complex;
-    log(): Complex;
-    sqr(): Complex;
-    sqrt(): Complex;
-    addReal(x: number): Complex;
-    add(x: Complex): Complex;
-    subReal(x: number): Complex;
-    static subFromReal(x: number, y: Complex): Complex;
-    sub(x: Complex): Complex;
-    mulReal(x: number): Complex;
-    mul(x: Complex): Complex;
-    divReal(x: number): Complex;
-    div(x: Complex): Complex;
-    static divFromReal(x: number, y: Complex): Complex;
-    powInt(x: number): Complex;
-    powReal(x: number): Complex;
-    pow(x: Complex): Complex;
-}
+export default class Complex {
+    static readonly I: Complex;
+    static readonly ZERO: Complex;
+    static readonly ONE: Complex;
+    static readonly TWO: Complex;
+    static readonly NaN: Complex;
+    static readonly INFINITY: Complex;
+    readonly re: number;
+    readonly im: number;
+    constructor(re: number, im?: number);
+    toString(): string;
+    toNumber(eps: number): number;
+    isNaN(): boolean;
+    isInfinite(): boolean;
+    isFinite(): boolean;
+    equals(x: Complex): boolean;
+    fuzzyEquals(x: Complex, eps: number): boolean;
+    static expj(arg: number): Complex;
+    static fromPolar(abs: number, arg: number): Complex;
+    abs(): number;
+    arg(): number;
+    conj(): Complex;
+    neg(): Complex;
+    reciprocal(): Complex;
+    exp(): Complex;
+    log(): Complex;
+    sqr(): Complex;
+    sqrt(): Complex;
+    addReal(x: number): Complex;
+    add(x: Complex): Complex;
+    subReal(x: number): Complex;
+    static subFromReal(x: number, y: Complex): Complex;
+    sub(x: Complex): Complex;
+    mulReal(x: number): Complex;
+    mul(x: Complex): Complex;
+    divReal(x: number): Complex;
+    div(x: Complex): Complex;
+    static divFromReal(x: number, y: Complex): Complex;
+    powInt(x: number): Complex;
+    powReal(x: number): Complex;
+    pow(x: Complex): Complex;
+}
+//# sourceMappingURL=Complex.d.ts.map

package/math/Complex.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"Complex.d.ts","sourceRoot":"","sources":["../../src/math/Complex.ts"],"names":[],"mappings":"AAKA,MAAM,CAAC,OAAO,OAAO,OAAO;IAOzB,gBAAuB,CAAC,UAAqB;IAG7C,gBAAuB,IAAI,UAAkB;IAG7C,gBAAuB,GAAG,UAAkB;IAG5C,gBAAuB,GAAG,UAAkB;IAG5C,gBAAuB,GAAG,UAAyB;IAGnD,gBAAuB,QAAQ,UAAmC;IAOlE,SAAgB,EAAE,EAAE,MAAM,CAAC;IAG3B,SAAgB,EAAE,EAAE,MAAM,CAAC;gBAYP,EAAE,EAAE,MAAM,EAAE,EAAE,GAAE,MAAU;IASvC,QAAQ,IAAK,MAAM;IAOnB,QAAQ,CAAE,GAAG,EAAE,MAAM,GAAI,MAAM;IAS/B,KAAK,IAAK,OAAO;IAMjB,UAAU;IAMV,QAAQ;IAMR,MAAM,CAAE,CAAC,EAAE,OAAO,GAAI,OAAO;IAO7B,WAAW,CAAE,CAAC,EAAE,OAAO,EAAE,GAAG,EAAE,MAAM,GAAI,OAAO;WAQxC,IAAI,CAAE,GAAG,EAAE,MAAM,GAAI,OAAO;WAM5B,SAAS,CAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAI,OAAO;IAQrD,GAAG,IAAK,MAAM;IAMd,GAAG,IAAK,MAAM;IAMd,IAAI,IAAK,OAAO;IAMhB,GAAG,IAAK,OAAO;IAMf,UAAU,IAAK,OAAO;IActB,GAAG,IAAK,OAAO;IAMf,GAAG,IAAK,OAAO;IAMf,GAAG,IAAK,OAAO;IAMf,IAAI,IAAK,OAAO;IAWhB,OAAO,CAAE,CAAC,EAAE,MAAM,GAAI,OAAO;IAM7B,GAAG,CAAE,CAAC,EAAE,OAAO,GAAI,OAAO;IAM1B,OAAO,CAAE,CAAC,EAAE,MAAM,GAAI,OAAO;WAMtB,WAAW,CAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,OAAO,GAAI,OAAO;IAMpD,GAAG,CAAE,CAAC,EAAE,OAAO,GAAI,OAAO;IAM1B,OAAO,CAAE,CAAC,EAAE,MAAM,GAAI,OAAO;IAM7B,GAAG,CAAE,CAAC,EAAE,OAAO,GAAI,OAAO;IAM1B,OAAO,CAAE,CAAC,EAAE,MAAM,GAAI,OAAO;IAM7B,GAAG,CAAE,CAAC,EAAE,OAAO,GAAI,OAAO;WAOnB,WAAW,CAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,OAAO,GAAI,OAAO;IAOpD,MAAM,CAAE,CAAC,EAAE,MAAM,GAAI,OAAO;IAQ5B,OAAO,CAAE,CAAC,EAAE,MAAM,GAAI,OAAO;IAM7B,GAAG,CAAE,CAAC,EAAE,OAAO,GAAI,OAAO;CAGhC"}