@jtarrio/signals 0.9.0

package/dist/dsp/converters.js

@@ -0,0 +1,43 @@
+// Copyright 2024 Jacobo Tarrio Barreiro. All rights reserved.
+// Copyright 2013 Google Inc. All rights reserved.
+//
+// Licensed under the Apache License, Version 2.0 (the "License");
+// you may not use this file except in compliance with the License.
+// You may obtain a copy of the License at
+//
+//     http://www.apache.org/licenses/LICENSE-2.0
+//
+// Unless required by applicable law or agreed to in writing, software
+// distributed under the License is distributed on an "AS IS" BASIS,
+// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+// See the License for the specific language governing permissions and
+// limitations under the License.
+import { IqPool } from "./buffers.js";
+/**
+ * Converts the given buffer of unsigned 8-bit samples into a pair of
+ * 32-bit floating-point sample streams.
+ */
+export class U8ToFloat32 {
+    /** @param length The expected length of each sample block. */
+    constructor(length) {
+        this.pool = new IqPool(4, length);
+    }
+    pool;
+    /**
+     * @param input A buffer containing the unsigned 8-bit samples.
+     * @returns An array that contains first the I stream and then the Q stream.
+     */
+    convert(input) {
+        let u8 = new Uint8Array(input);
+        const len = u8.length / 2;
+        let out = this.pool.get(len);
+        const outI = out[0];
+        const outQ = out[1];
+        for (let i = 0; i < len; ++i) {
+            outI[i] = u8[2 * i] / 128 - 0.995;
+            outQ[i] = u8[2 * i + 1] / 128 - 0.995;
+        }
+        return out;
+    }
+}
+//# sourceMappingURL=converters.js.map

package/dist/dsp/converters.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"converters.js","sourceRoot":"","sources":["../../src/dsp/converters.ts"],"names":[],"mappings":"AAAA,8DAA8D;AAC9D,kDAAkD;AAClD,EAAE;AACF,kEAAkE;AAClE,mEAAmE;AACnE,0CAA0C;AAC1C,EAAE;AACF,iDAAiD;AACjD,EAAE;AACF,sEAAsE;AACtE,oEAAoE;AACpE,2EAA2E;AAC3E,sEAAsE;AACtE,iCAAiC;AAEjC,OAAO,EAAE,MAAM,EAAE,MAAM,cAAc,CAAC;AAEtC;;;GAGG;AACH,MAAM,OAAO,WAAW;IACtB,8DAA8D;IAC9D,YAAY,MAAe;QACzB,IAAI,CAAC,IAAI,GAAG,IAAI,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;IACpC,CAAC;IAEO,IAAI,CAAS;IAErB;;;OAGG;IACH,OAAO,CAAC,KAAkB;QACxB,IAAI,EAAE,GAAG,IAAI,UAAU,CAAC,KAAK,CAAC,CAAC;QAC/B,MAAM,GAAG,GAAG,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC;QAC1B,IAAI,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAC7B,MAAM,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QACpB,MAAM,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QACpB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC,EAAE,CAAC;YAC7B,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,KAAK,CAAC;YAClC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,KAAK,CAAC;QACxC,CAAC;QACD,OAAO,GAAG,CAAC;IACb,CAAC;CACF"}

package/dist/dsp/demodulators.d.ts

@@ -0,0 +1,68 @@
+/** The sideband to demodulate. */
+export declare enum Sideband {
+    Upper = 0,
+    Lower = 1
+}
+/** A class to demodulate a USB or LSB signal. */
+export declare class SSBDemodulator {
+    /**
+     * @param sideband The sideband to demodulate.
+     * @param kernelLen The length of the Hilbert filter kernel to use.
+     */
+    constructor(sideband: Sideband, kernelLen: number);
+    private filterDelay;
+    private filterHilbert;
+    private hilbertMul;
+    /** Switches the demodulator's sideband on the fly. */
+    setSideband(sideband: Sideband): void;
+    /** Demodulates the given I/Q samples into the real output. */
+    demodulate(I: Float32Array, Q: Float32Array, out: Float32Array): void;
+}
+/** A class to demodulate an AM signal. */
+export declare class AMDemodulator {
+    /**
+     * @param sampleRate The signal's sample rate.
+     */
+    constructor(sampleRate: number);
+    private alpha;
+    private carrierAmplitude;
+    /** Demodulates the given I/Q samples into the real output. */
+    demodulate(I: Float32Array, Q: Float32Array, out: Float32Array): void;
+}
+/** A class to demodulate an FM signal. */
+export declare class FMDemodulator {
+    /**
+     * @param maxDeviation The maximum deviation for the signal, as a fraction of the sample rate.
+     */
+    constructor(maxDeviation: number);
+    private mul;
+    private lI;
+    private lQ;
+    /** Changes the maximum deviation. */
+    setMaxDeviation(maxDeviation: number): void;
+    /** Demodulates the given I/Q samples into the real output. */
+    demodulate(I: Float32Array, Q: Float32Array, out: Float32Array): void;
+}
+/** A class to demodulate the stereo signal in a demodulated FM signal. */
+export declare class StereoSeparator {
+    /**
+     * @param sampleRate The sample rate for the input signal.
+     * @param pilotFreq The frequency of the pilot tone.
+     */
+    constructor(sampleRate: number, pilotFreq: number);
+    private pool;
+    private pll;
+    /**
+     * Locks on to the pilot tone and uses it to demodulate the stereo audio.
+     * @param samples The original audio stream.
+     * @returns An object with a key 'found' that tells whether a
+     *     consistent stereo pilot tone was detected and a key 'diff'
+     *     that contains the original stream demodulated with the
+     *     reconstructed stereo carrier.
+     */
+    separate(samples: Float32Array): {
+        found: boolean;
+        diff: Float32Array;
+    };
+}
+//# sourceMappingURL=demodulators.d.ts.map

package/dist/dsp/demodulators.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"demodulators.d.ts","sourceRoot":"","sources":["../../src/dsp/demodulators.ts"],"names":[],"mappings":"AAoBA,kCAAkC;AAClC,oBAAY,QAAQ;IAClB,KAAK,IAAA;IACL,KAAK,IAAA;CACN;AAED,iDAAiD;AACjD,qBAAa,cAAc;IACzB;;;OAGG;gBACS,QAAQ,EAAE,QAAQ,EAAE,SAAS,EAAE,MAAM;IAOjD,OAAO,CAAC,WAAW,CAAY;IAC/B,OAAO,CAAC,aAAa,CAAY;IACjC,OAAO,CAAC,UAAU,CAAS;IAE3B,sDAAsD;IACtD,WAAW,CAAC,QAAQ,EAAE,QAAQ;IAI9B,8DAA8D;IAC9D,UAAU,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,EAAE,GAAG,EAAE,YAAY;CAU/D;AAED,0CAA0C;AAC1C,qBAAa,aAAa;IACxB;;OAEG;gBACS,UAAU,EAAE,MAAM;IAK9B,OAAO,CAAC,KAAK,CAAS;IACtB,OAAO,CAAC,gBAAgB,CAAS;IAEjC,8DAA8D;IAC9D,UAAU,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,EAAE,GAAG,EAAE,YAAY;CAa/D;AAED,0CAA0C;AAC1C,qBAAa,aAAa;IACxB;;OAEG;gBACS,YAAY,EAAE,MAAM;IAMhC,OAAO,CAAC,GAAG,CAAS;IACpB,OAAO,CAAC,EAAE,CAAS;IACnB,OAAO,CAAC,EAAE,CAAS;IAEnB,qCAAqC;IACrC,eAAe,CAAC,YAAY,EAAE,MAAM;IAIpC,8DAA8D;IAC9D,UAAU,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,EAAE,GAAG,EAAE,YAAY;CAe/D;AAED,0EAA0E;AAC1E,qBAAa,eAAe;IAC1B;;;OAGG;gBACS,UAAU,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM;IAKjD,OAAO,CAAC,IAAI,CAAc;IAC1B,OAAO,CAAC,GAAG,CAAM;IAEjB;;;;;;;OAOG;IACH,QAAQ,CAAC,OAAO,EAAE,YAAY,GAAG;QAAE,KAAK,EAAE,OAAO,CAAC;QAAC,IAAI,EAAE,YAAY,CAAA;KAAE;CAaxE"}

package/dist/dsp/demodulators.js

@@ -0,0 +1,149 @@
+// Copyright 2024 Jacobo Tarrio Barreiro. All rights reserved.
+// Copyright 2013 Google Inc. All rights reserved.
+//
+// Licensed under the Apache License, Version 2.0 (the "License");
+// you may not use this file except in compliance with the License.
+// You may obtain a copy of the License at
+//
+//     http://www.apache.org/licenses/LICENSE-2.0
+//
+// Unless required by applicable law or agreed to in writing, software
+// distributed under the License is distributed on an "AS IS" BASIS,
+// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+// See the License for the specific language governing permissions and
+// limitations under the License.
+import { makeHilbertKernel } from "./coefficients.js";
+import { decay, FIRFilter, PLL } from "./filters.js";
+import { Float32Pool } from "./buffers.js";
+import { atan2 } from "./math.js";
+/** The sideband to demodulate. */
+export var Sideband;
+(function (Sideband) {
+    Sideband[Sideband["Upper"] = 0] = "Upper";
+    Sideband[Sideband["Lower"] = 1] = "Lower";
+})(Sideband || (Sideband = {}));
+/** A class to demodulate a USB or LSB signal. */
+export class SSBDemodulator {
+    /**
+     * @param sideband The sideband to demodulate.
+     * @param kernelLen The length of the Hilbert filter kernel to use.
+     */
+    constructor(sideband, kernelLen) {
+        let hilbert = makeHilbertKernel(kernelLen);
+        this.filterDelay = new FIRFilter(hilbert);
+        this.filterHilbert = new FIRFilter(hilbert);
+        this.hilbertMul = sideband == Sideband.Upper ? -1 : 1;
+    }
+    filterDelay;
+    filterHilbert;
+    hilbertMul;
+    /** Switches the demodulator's sideband on the fly. */
+    setSideband(sideband) {
+        this.hilbertMul = sideband == Sideband.Upper ? -1 : 1;
+    }
+    /** Demodulates the given I/Q samples into the real output. */
+    demodulate(I, Q, out) {
+        this.filterDelay.loadSamples(I);
+        this.filterHilbert.loadSamples(Q);
+        for (let i = 0; i < out.length; ++i) {
+            out[i] =
+                (this.filterDelay.getDelayed(i) +
+                    this.filterHilbert.get(i) * this.hilbertMul) /
+                    2;
+        }
+    }
+}
+/** A class to demodulate an AM signal. */
+export class AMDemodulator {
+    /**
+     * @param sampleRate The signal's sample rate.
+     */
+    constructor(sampleRate) {
+        this.alpha = decay(sampleRate, 0.5);
+        this.carrierAmplitude = 0;
+    }
+    alpha;
+    carrierAmplitude;
+    /** Demodulates the given I/Q samples into the real output. */
+    demodulate(I, Q, out) {
+        const alpha = this.alpha;
+        let carrierAmplitude = this.carrierAmplitude;
+        for (let i = 0; i < out.length; ++i) {
+            const vI = I[i];
+            const vQ = Q[i];
+            const power = vI * vI + vQ * vQ;
+            const amplitude = Math.sqrt(power);
+            carrierAmplitude += alpha * (amplitude - carrierAmplitude);
+            out[i] = carrierAmplitude == 0 ? 0 : amplitude / carrierAmplitude - 1;
+        }
+        this.carrierAmplitude = carrierAmplitude;
+    }
+}
+/** A class to demodulate an FM signal. */
+export class FMDemodulator {
+    /**
+     * @param maxDeviation The maximum deviation for the signal, as a fraction of the sample rate.
+     */
+    constructor(maxDeviation) {
+        this.mul = 1 / (2 * Math.PI * maxDeviation);
+        this.lI = 0;
+        this.lQ = 0;
+    }
+    mul;
+    lI;
+    lQ;
+    /** Changes the maximum deviation. */
+    setMaxDeviation(maxDeviation) {
+        this.mul = 1 / (2 * Math.PI * maxDeviation);
+    }
+    /** Demodulates the given I/Q samples into the real output. */
+    demodulate(I, Q, out) {
+        const mul = this.mul;
+        let lI = this.lI;
+        let lQ = this.lQ;
+        for (let i = 0; i < I.length; ++i) {
+            // s = s * conj(l)
+            let real = lI * I[i] + lQ * Q[i];
+            let imag = lI * Q[i] - I[i] * lQ;
+            lI = I[i];
+            lQ = Q[i];
+            out[i] = atan2(imag, real) * mul;
+        }
+        this.lI = lI;
+        this.lQ = lQ;
+    }
+}
+/** A class to demodulate the stereo signal in a demodulated FM signal. */
+export class StereoSeparator {
+    /**
+     * @param sampleRate The sample rate for the input signal.
+     * @param pilotFreq The frequency of the pilot tone.
+     */
+    constructor(sampleRate, pilotFreq) {
+        this.pool = new Float32Pool(4);
+        this.pll = new PLL(sampleRate, pilotFreq, 10);
+    }
+    pool;
+    pll;
+    /**
+     * Locks on to the pilot tone and uses it to demodulate the stereo audio.
+     * @param samples The original audio stream.
+     * @returns An object with a key 'found' that tells whether a
+     *     consistent stereo pilot tone was detected and a key 'diff'
+     *     that contains the original stream demodulated with the
+     *     reconstructed stereo carrier.
+     */
+    separate(samples) {
+        let out = this.pool.get(samples.length);
+        for (let i = 0; i < samples.length; ++i) {
+            this.pll.add(samples[i]);
+            // Multiply by 4 instead of 2 so 'out' has the same level as the input samples.
+            out[i] = samples[i] * this.pll.sin * this.pll.cos * 4;
+        }
+        return {
+            found: this.pll.locked,
+            diff: out,
+        };
+    }
+}
+//# sourceMappingURL=demodulators.js.map

package/dist/dsp/demodulators.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"demodulators.js","sourceRoot":"","sources":["../../src/dsp/demodulators.ts"],"names":[],"mappings":"AAAA,8DAA8D;AAC9D,kDAAkD;AAClD,EAAE;AACF,kEAAkE;AAClE,mEAAmE;AACnE,0CAA0C;AAC1C,EAAE;AACF,iDAAiD;AACjD,EAAE;AACF,sEAAsE;AACtE,oEAAoE;AACpE,2EAA2E;AAC3E,sEAAsE;AACtE,iCAAiC;AAEjC,OAAO,EAAE,iBAAiB,EAAE,MAAM,mBAAmB,CAAC;AACtD,OAAO,EAAE,KAAK,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,cAAc,CAAC;AACrD,OAAO,EAAE,WAAW,EAAE,MAAM,cAAc,CAAC;AAC3C,OAAO,EAAE,KAAK,EAAE,MAAM,WAAW,CAAC;AAElC,kCAAkC;AAClC,MAAM,CAAN,IAAY,QAGX;AAHD,WAAY,QAAQ;IAClB,yCAAK,CAAA;IACL,yCAAK,CAAA;AACP,CAAC,EAHW,QAAQ,KAAR,QAAQ,QAGnB;AAED,iDAAiD;AACjD,MAAM,OAAO,cAAc;IACzB;;;OAGG;IACH,YAAY,QAAkB,EAAE,SAAiB;QAC/C,IAAI,OAAO,GAAG,iBAAiB,CAAC,SAAS,CAAC,CAAC;QAC3C,IAAI,CAAC,WAAW,GAAG,IAAI,SAAS,CAAC,OAAO,CAAC,CAAC;QAC1C,IAAI,CAAC,aAAa,GAAG,IAAI,SAAS,CAAC,OAAO,CAAC,CAAC;QAC5C,IAAI,CAAC,UAAU,GAAG,QAAQ,IAAI,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACxD,CAAC;IAEO,WAAW,CAAY;IACvB,aAAa,CAAY;IACzB,UAAU,CAAS;IAE3B,sDAAsD;IACtD,WAAW,CAAC,QAAkB;QAC5B,IAAI,CAAC,UAAU,GAAG,QAAQ,IAAI,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACxD,CAAC;IAED,8DAA8D;IAC9D,UAAU,CAAC,CAAe,EAAE,CAAe,EAAE,GAAiB;QAC5D,IAAI,CAAC,WAAW,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC;QAChC,IAAI,CAAC,aAAa,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC;QAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,EAAE,CAAC,EAAE,CAAC;YACpC,GAAG,CAAC,CAAC,CAAC;gBACJ,CAAC,IAAI,CAAC,WAAW,CAAC,UAAU,CAAC,CAAC,CAAC;oBAC7B,IAAI,CAAC,aAAa,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC;oBAC9C,CAAC,CAAC;QACN,CAAC;IACH,CAAC;CACF;AAED,0CAA0C;AAC1C,MAAM,OAAO,aAAa;IACxB;;OAEG;IACH,YAAY,UAAkB;QAC5B,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC,UAAU,EAAE,GAAG,CAAC,CAAC;QACpC,IAAI,CAAC,gBAAgB,GAAG,CAAC,CAAC;IAC5B,CAAC;IAEO,KAAK,CAAS;IACd,gBAAgB,CAAS;IAEjC,8DAA8D;IAC9D,UAAU,CAAC,CAAe,EAAE,CAAe,EAAE,GAAiB;QAC5D,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC;QACzB,IAAI,gBAAgB,GAAG,IAAI,CAAC,gBAAgB,CAAC;QAC7C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,EAAE,CAAC,EAAE,CAAC;YACpC,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YAChB,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YAChB,MAAM,KAAK,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;YAChC,MAAM,SAAS,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;YACnC,gBAAgB,IAAI,KAAK,GAAG,CAAC,SAAS,GAAG,gBAAgB,CAAC,CAAC;YAC3D,GAAG,CAAC,CAAC,CAAC,GAAG,gBAAgB,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,SAAS,GAAG,gBAAgB,GAAG,CAAC,CAAC;QACxE,CAAC;QACD,IAAI,CAAC,gBAAgB,GAAG,gBAAgB,CAAC;IAC3C,CAAC;CACF;AAED,0CAA0C;AAC1C,MAAM,OAAO,aAAa;IACxB;;OAEG;IACH,YAAY,YAAoB;QAC9B,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,YAAY,CAAC,CAAC;QAC5C,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;QACZ,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;IACd,CAAC;IAEO,GAAG,CAAS;IACZ,EAAE,CAAS;IACX,EAAE,CAAS;IAEnB,qCAAqC;IACrC,eAAe,CAAC,YAAoB;QAClC,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,YAAY,CAAC,CAAC;IAC9C,CAAC;IAED,8DAA8D;IAC9D,UAAU,CAAC,CAAe,EAAE,CAAe,EAAE,GAAiB;QAC5D,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC;QACrB,IAAI,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;QACjB,IAAI,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;QACjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,EAAE,CAAC,EAAE,CAAC;YAClC,kBAAkB;YAClB,IAAI,IAAI,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YACjC,IAAI,IAAI,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;YACjC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YACV,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YACV,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,GAAG,CAAC;QACnC,CAAC;QACD,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;QACb,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;IACf,CAAC;CACF;AAED,0EAA0E;AAC1E,MAAM,OAAO,eAAe;IAC1B;;;OAGG;IACH,YAAY,UAAkB,EAAE,SAAiB;QAC/C,IAAI,CAAC,IAAI,GAAG,IAAI,WAAW,CAAC,CAAC,CAAC,CAAC;QAC/B,IAAI,CAAC,GAAG,GAAG,IAAI,GAAG,CAAC,UAAU,EAAE,SAAS,EAAE,EAAE,CAAC,CAAC;IAChD,CAAC;IAEO,IAAI,CAAc;IAClB,GAAG,CAAM;IAEjB;;;;;;;OAOG;IACH,QAAQ,CAAC,OAAqB;QAC5B,IAAI,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;QACxC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,EAAE,CAAC,EAAE,CAAC;YACxC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;YACzB,+EAA+E;YAC/E,GAAG,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC;QACxD,CAAC;QAED,OAAO;YACL,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,MAAM;YACtB,IAAI,EAAE,GAAG;SACV,CAAC;IACJ,CAAC;CACF"}

package/dist/dsp/fft.d.ts

@@ -0,0 +1,54 @@
+import { Float32RingBuffer } from "./buffers.js";
+/** Fast Fourier Transform implementation. */
+/**
+ * Returns the length of the FFT for a given array length.
+ *
+ * This FFT implementation only works in power-of-2 lengths,
+ * so this function returns the next available length.
+ */
+export declare function actualLength(minimumLength: number): number;
+/**
+ * The output of the transform functions.
+ *
+ * The first array contains the real parts, and the second array contains the imaginary parts.
+ */
+export type FFTOutput = [Float32Array, Float32Array];
+/** Fast Fourier Transform and reverse transform with a given length. */
+export declare class FFT {
+    length: number;
+    /**
+     * Returns an FFT instance that fits the given length.
+     *
+     * The actual length may be greater than the given length if it
+     * is not a power of 2.
+     */
+    static ofLength(minimumLength: number): FFT;
+    private constructor();
+    private revIndex;
+    private fwd;
+    private bwd;
+    private copy;
+    private out;
+    private window;
+    /** Sets the window function for this FFT. */
+    setWindow(window: Float32Array): void;
+    /**
+     * Transforms the given time-domain input.
+     * @param real An array of real parts.
+     * @param imag An array of imaginary parts.
+     * @return The output of the transform.
+     */
+    transform(real: Float32Array, imag: Float32Array): FFTOutput;
+    transform(real: number[], imag: number[]): FFTOutput;
+    transformCircularBuffers(real: Float32RingBuffer, imag: Float32RingBuffer): FFTOutput;
+    /**
+     * Does a reverse transform of the given frequency-domain input.
+     * The input and output arrays must be the same length as the FFT.
+     * @param real An array of real parts.
+     * @param imag An array of imaginary parts.
+     * @return The output of the reverse transform.
+     */
+    reverse(real: Float32Array, imag: Float32Array): FFTOutput;
+    reverse(real: number[], imag: number[]): FFTOutput;
+}
+//# sourceMappingURL=fft.d.ts.map

package/dist/dsp/fft.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"fft.d.ts","sourceRoot":"","sources":["../../src/dsp/fft.ts"],"names":[],"mappings":"AAcA,OAAO,EAAE,iBAAiB,EAAU,MAAM,cAAc,CAAC;AAEzD,6CAA6C;AAE7C;;;;;GAKG;AACH,wBAAgB,YAAY,CAAC,aAAa,EAAE,MAAM,GAAG,MAAM,CAM1D;AAED;;;;GAIG;AACH,MAAM,MAAM,SAAS,GAAG,CAAC,YAAY,EAAE,YAAY,CAAC,CAAC;AAErD,wEAAwE;AACxE,qBAAa,GAAG;IAWa,MAAM,EAAE,MAAM;IAVzC;;;;;OAKG;IACH,MAAM,CAAC,QAAQ,CAAC,aAAa,EAAE,MAAM,GAAG,GAAG;IAI3C,OAAO;IAWP,OAAO,CAAC,QAAQ,CAAa;IAC7B,OAAO,CAAC,GAAG,CAAiB;IAC5B,OAAO,CAAC,GAAG,CAAiB;IAC5B,OAAO,CAAC,IAAI,CAAS;IACrB,OAAO,CAAC,GAAG,CAAS;IACpB,OAAO,CAAC,MAAM,CAAe;IAE7B,6CAA6C;IAC7C,SAAS,CAAC,MAAM,EAAE,YAAY;IAI9B;;;;;OAKG;IACH,SAAS,CAAC,IAAI,EAAE,YAAY,EAAE,IAAI,EAAE,YAAY,GAAG,SAAS;IAC5D,SAAS,CAAC,IAAI,EAAE,MAAM,EAAE,EAAE,IAAI,EAAE,MAAM,EAAE,GAAG,SAAS;IAepD,wBAAwB,CACtB,IAAI,EAAE,iBAAiB,EACvB,IAAI,EAAE,iBAAiB,GACtB,SAAS;IAQZ;;;;;;OAMG;IACH,OAAO,CAAC,IAAI,EAAE,YAAY,EAAE,IAAI,EAAE,YAAY,GAAG,SAAS;IAC1D,OAAO,CAAC,IAAI,EAAE,MAAM,EAAE,EAAE,IAAI,EAAE,MAAM,EAAE,GAAG,SAAS;CAcnD"}

package/dist/dsp/fft.js

@@ -0,0 +1,175 @@
+// Copyright 2024 Jacobo Tarrio Barreiro. All rights reserved.
+//
+// Licensed under the Apache License, Version 2.0 (the "License");
+// you may not use this file except in compliance with the License.
+// You may obtain a copy of the License at
+//
+//     http://www.apache.org/licenses/LICENSE-2.0
+//
+// Unless required by applicable law or agreed to in writing, software
+// distributed under the License is distributed on an "AS IS" BASIS,
+// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+// See the License for the specific language governing permissions and
+// limitations under the License.
+import { IqPool } from "./buffers.js";
+/** Fast Fourier Transform implementation. */
+/**
+ * Returns the length of the FFT for a given array length.
+ *
+ * This FFT implementation only works in power-of-2 lengths,
+ * so this function returns the next available length.
+ */
+export function actualLength(minimumLength) {
+    if (minimumLength < 2)
+        return 0;
+    if (((minimumLength - 1) & minimumLength) == 0)
+        return minimumLength;
+    let realLength = 1;
+    while (realLength < minimumLength)
+        realLength <<= 1;
+    return realLength;
+}
+/** Fast Fourier Transform and reverse transform with a given length. */
+export class FFT {
+    length;
+    /**
+     * Returns an FFT instance that fits the given length.
+     *
+     * The actual length may be greater than the given length if it
+     * is not a power of 2.
+     */
+    static ofLength(minimumLength) {
+        return new FFT(actualLength(minimumLength));
+    }
+    constructor(length) {
+        this.length = length;
+        this.revIndex = reversedBitIndices(length);
+        let [fwd, bwd] = makeFftCoefficients(length);
+        this.fwd = fwd;
+        this.bwd = bwd;
+        this.copy = new IqPool(2, length);
+        this.out = new IqPool(2, length);
+        this.window = new Float32Array(length);
+        this.window.fill(1);
+    }
+    revIndex;
+    fwd;
+    bwd;
+    copy;
+    out;
+    window;
+    /** Sets the window function for this FFT. */
+    setWindow(window) {
+        this.window.set(window);
+    }
+    transform(real, imag) {
+        const length = this.length;
+        let [outReal, outImag] = this.out.get(length);
+        outReal.fill(0);
+        outImag.fill(0);
+        for (let i = 0; i < length && i < real.length && i < imag.length; ++i) {
+            const ri = this.revIndex[i];
+            outReal[ri] = (this.window[i] * real[i]) / length;
+            outImag[ri] = (this.window[i] * imag[i]) / length;
+        }
+        doFastTransform(this.length, this.fwd, outReal, outImag);
+        return [outReal, outImag];
+    }
+    transformCircularBuffers(real, imag) {
+        const length = this.length;
+        let [copyReal, copyImag] = this.copy.get(length);
+        real.copyTo(copyReal);
+        imag.copyTo(copyImag);
+        return this.transform(copyReal, copyImag);
+    }
+    reverse(real, imag) {
+        const length = this.length;
+        let [outReal, outImag] = this.out.get(length);
+        outReal.fill(0);
+        outImag.fill(0);
+        for (let i = 0; i < length && i < real.length && i < imag.length; ++i) {
+            const ri = this.revIndex[i];
+            outReal[ri] = this.window[i] * real[i];
+            outImag[ri] = this.window[i] * imag[i];
+        }
+        doFastTransform(this.length, this.bwd, outReal, outImag);
+        return [outReal, outImag];
+    }
+}
+/** Performs a fast direct or reverse transform in place. */
+function doFastTransform(length, coefs, real, imag) {
+    for (let dftSize = 2, coeffBin = 0; dftSize <= length; dftSize *= 2, ++coeffBin) {
+        const binCoefficients = coefs[coeffBin];
+        const halfDftSize = dftSize / 2;
+        for (let dftStart = 0; dftStart < length; dftStart += dftSize) {
+            for (let i = 0; i < halfDftSize; ++i) {
+                const near = dftStart + i;
+                const far = near + halfDftSize;
+                const evenReal = real[near];
+                const evenImag = imag[near];
+                const cr = binCoefficients.real[i];
+                const ci = binCoefficients.imag[i];
+                const or = real[far];
+                const oi = imag[far];
+                const oddReal = cr * or - ci * oi;
+                const oddImag = cr * oi + ci * or;
+                real[near] = evenReal + oddReal;
+                imag[near] = evenImag + oddImag;
+                real[far] = evenReal - oddReal;
+                imag[far] = evenImag - oddImag;
+            }
+        }
+    }
+}
+/** Builds a triangle of direct and reverse FFT coefficients for the given length. */
+function makeFftCoefficients(length) {
+    let numBits = getNumBits(length);
+    let fwd = [];
+    let bwd = [];
+    for (let bin = 0, halfSize = 1; bin < numBits; ++bin, halfSize *= 2) {
+        fwd.push({
+            real: new Float32Array(halfSize),
+            imag: new Float32Array(halfSize),
+        });
+        bwd.push({
+            real: new Float32Array(halfSize),
+            imag: new Float32Array(halfSize),
+        });
+        for (let i = 0; i < halfSize; ++i) {
+            const fwdAngle = (-1 * Math.PI * i) / halfSize;
+            fwd[bin].real[i] = Math.cos(fwdAngle);
+            fwd[bin].imag[i] = Math.sin(fwdAngle);
+            const bwdAngle = (Math.PI * i) / halfSize;
+            bwd[bin].real[i] = Math.cos(bwdAngle);
+            bwd[bin].imag[i] = Math.sin(bwdAngle);
+        }
+    }
+    return [fwd, bwd];
+}
+/** Builds an array of numbers with their bits reversed. */
+function reversedBitIndices(length) {
+    const numBits = getNumBits(length);
+    let output = new Int32Array(length);
+    for (let i = 0; i < length; ++i) {
+        output[i] = reverseBits(i, numBits);
+    }
+    return output;
+}
+/** Returns how many bits we need to fit 'length' distinct values. */
+function getNumBits(length) {
+    let numBits = 0;
+    for (let shifted = length - 1; shifted > 0; shifted >>= 1)
+        ++numBits;
+    return numBits;
+}
+/** Reverses the bits in a number. */
+function reverseBits(num, bits) {
+    let output = 0;
+    for (let b = 0; b < bits; ++b) {
+        output <<= 1;
+        output |= num & 1;
+        num >>= 1;
+    }
+    return output;
+}
+//# sourceMappingURL=fft.js.map

package/dist/dsp/fft.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"fft.js","sourceRoot":"","sources":["../../src/dsp/fft.ts"],"names":[],"mappings":"AAAA,8DAA8D;AAC9D,EAAE;AACF,kEAAkE;AAClE,mEAAmE;AACnE,0CAA0C;AAC1C,EAAE;AACF,iDAAiD;AACjD,EAAE;AACF,sEAAsE;AACtE,oEAAoE;AACpE,2EAA2E;AAC3E,sEAAsE;AACtE,iCAAiC;AAEjC,OAAO,EAAqB,MAAM,EAAE,MAAM,cAAc,CAAC;AAEzD,6CAA6C;AAE7C;;;;;GAKG;AACH,MAAM,UAAU,YAAY,CAAC,aAAqB;IAChD,IAAI,aAAa,GAAG,CAAC;QAAE,OAAO,CAAC,CAAC;IAChC,IAAI,CAAC,CAAC,aAAa,GAAG,CAAC,CAAC,GAAG,aAAa,CAAC,IAAI,CAAC;QAAE,OAAO,aAAa,CAAC;IACrE,IAAI,UAAU,GAAG,CAAC,CAAC;IACnB,OAAO,UAAU,GAAG,aAAa;QAAE,UAAU,KAAK,CAAC,CAAC;IACpD,OAAO,UAAU,CAAC;AACpB,CAAC;AASD,wEAAwE;AACxE,MAAM,OAAO,GAAG;IAWa;IAV3B;;;;;OAKG;IACH,MAAM,CAAC,QAAQ,CAAC,aAAqB;QACnC,OAAO,IAAI,GAAG,CAAC,YAAY,CAAC,aAAa,CAAC,CAAC,CAAC;IAC9C,CAAC;IAED,YAA2B,MAAc;QAAd,WAAM,GAAN,MAAM,CAAQ;QACvC,IAAI,CAAC,QAAQ,GAAG,kBAAkB,CAAC,MAAM,CAAC,CAAC;QAC3C,IAAI,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,mBAAmB,CAAC,MAAM,CAAC,CAAC;QAC7C,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC;QACf,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,IAAI,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;QAClC,IAAI,CAAC,GAAG,GAAG,IAAI,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;QACjC,IAAI,CAAC,MAAM,GAAG,IAAI,YAAY,CAAC,MAAM,CAAC,CAAC;QACvC,IAAI,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACtB,CAAC;IAEO,QAAQ,CAAa;IACrB,GAAG,CAAiB;IACpB,GAAG,CAAiB;IACpB,IAAI,CAAS;IACb,GAAG,CAAS;IACZ,MAAM,CAAe;IAE7B,6CAA6C;IAC7C,SAAS,CAAC,MAAoB;QAC5B,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC;IAC1B,CAAC;IAUD,SAAS,CAA0B,IAAO,EAAE,IAAO;QACjD,MAAM,MAAM,GAAG,IAAI,CAAC,MAAM,CAAC;QAC3B,IAAI,CAAC,OAAO,EAAE,OAAO,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC;QAC9C,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAChB,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAChB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,IAAI,CAAC,GAAG,IAAI,CAAC,MAAM,IAAI,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,EAAE,CAAC,EAAE,CAAC;YACtE,MAAM,EAAE,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;YAC5B,OAAO,CAAC,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;YAClD,OAAO,CAAC,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;QACpD,CAAC;QACD,eAAe,CAAC,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC,GAAG,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QACzD,OAAO,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;IAC5B,CAAC;IAED,wBAAwB,CACtB,IAAuB,EACvB,IAAuB;QAEvB,MAAM,MAAM,GAAG,IAAI,CAAC,MAAM,CAAC;QAC3B,IAAI,CAAC,QAAQ,EAAE,QAAQ,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC;QACjD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC;QACtB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC;QACtB,OAAO,IAAI,CAAC,SAAS,CAAC,QAAQ,EAAE,QAAQ,CAAC,CAAC;IAC5C,CAAC;IAWD,OAAO,CAA0B,IAAO,EAAE,IAAO;QAC/C,MAAM,MAAM,GAAG,IAAI,CAAC,MAAM,CAAC;QAC3B,IAAI,CAAC,OAAO,EAAE,OAAO,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC;QAC9C,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAChB,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAChB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,IAAI,CAAC,GAAG,IAAI,CAAC,MAAM,IAAI,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,EAAE,CAAC,EAAE,CAAC;YACtE,MAAM,EAAE,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;YAC5B,OAAO,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;YACvC,OAAO,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;QACzC,CAAC;QACD,eAAe,CAAC,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC,GAAG,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QACzD,OAAO,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;IAC5B,CAAC;CACF;AAED,4DAA4D;AAC5D,SAAS,eAAe,CACtB,MAAc,EACd,KAAqB,EACrB,IAAkB,EAClB,IAAkB;IAElB,KACE,IAAI,OAAO,GAAG,CAAC,EAAE,QAAQ,GAAG,CAAC,EAC7B,OAAO,IAAI,MAAM,EACjB,OAAO,IAAI,CAAC,EAAE,EAAE,QAAQ,EACxB,CAAC;QACD,MAAM,eAAe,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC;QACxC,MAAM,WAAW,GAAG,OAAO,GAAG,CAAC,CAAC;QAChC,KAAK,IAAI,QAAQ,GAAG,CAAC,EAAE,QAAQ,GAAG,MAAM,EAAE,QAAQ,IAAI,OAAO,EAAE,CAAC;YAC9D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,WAAW,EAAE,EAAE,CAAC,EAAE,CAAC;gBACrC,MAAM,IAAI,GAAG,QAAQ,GAAG,CAAC,CAAC;gBAC1B,MAAM,GAAG,GAAG,IAAI,GAAG,WAAW,CAAC;gBAC/B,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC;gBAC5B,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC;gBAC5B,MAAM,EAAE,GAAG,eAAe,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;gBACnC,MAAM,EAAE,GAAG,eAAe,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;gBACnC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;gBACrB,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;gBACrB,MAAM,OAAO,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;gBAClC,MAAM,OAAO,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;gBAClC,IAAI,CAAC,IAAI,CAAC,GAAG,QAAQ,GAAG,OAAO,CAAC;gBAChC,IAAI,CAAC,IAAI,CAAC,GAAG,QAAQ,GAAG,OAAO,CAAC;gBAChC,IAAI,CAAC,GAAG,CAAC,GAAG,QAAQ,GAAG,OAAO,CAAC;gBAC/B,IAAI,CAAC,GAAG,CAAC,GAAG,QAAQ,GAAG,OAAO,CAAC;YACjC,CAAC;QACH,CAAC;IACH,CAAC;AACH,CAAC;AAKD,qFAAqF;AACrF,SAAS,mBAAmB,CAAC,MAAc;IACzC,IAAI,OAAO,GAAG,UAAU,CAAC,MAAM,CAAC,CAAC;IACjC,IAAI,GAAG,GAAmB,EAAE,CAAC;IAC7B,IAAI,GAAG,GAAmB,EAAE,CAAC;IAE7B,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,QAAQ,GAAG,CAAC,EAAE,GAAG,GAAG,OAAO,EAAE,EAAE,GAAG,EAAE,QAAQ,IAAI,CAAC,EAAE,CAAC;QACpE,GAAG,CAAC,IAAI,CAAC;YACP,IAAI,EAAE,IAAI,YAAY,CAAC,QAAQ,CAAC;YAChC,IAAI,EAAE,IAAI,YAAY,CAAC,QAAQ,CAAC;SACjC,CAAC,CAAC;QACH,GAAG,CAAC,IAAI,CAAC;YACP,IAAI,EAAE,IAAI,YAAY,CAAC,QAAQ,CAAC;YAChC,IAAI,EAAE,IAAI,YAAY,CAAC,QAAQ,CAAC;SACjC,CAAC,CAAC;QACH,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,EAAE,CAAC,EAAE,CAAC;YAClC,MAAM,QAAQ,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,QAAQ,CAAC;YAC/C,GAAG,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC;YACtC,GAAG,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC;YACtC,MAAM,QAAQ,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,QAAQ,CAAC;YAC1C,GAAG,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC;YACtC,GAAG,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC;QACxC,CAAC;IACH,CAAC;IAED,OAAO,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;AACpB,CAAC;AAED,2DAA2D;AAC3D,SAAS,kBAAkB,CAAC,MAAc;IACxC,MAAM,OAAO,GAAG,UAAU,CAAC,MAAM,CAAC,CAAC;IACnC,IAAI,MAAM,GAAG,IAAI,UAAU,CAAC,MAAM,CAAC,CAAC;IACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,EAAE,CAAC,EAAE,CAAC;QAChC,MAAM,CAAC,CAAC,CAAC,GAAG,WAAW,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;IACtC,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAED,qEAAqE;AACrE,SAAS,UAAU,CAAC,MAAc;IAChC,IAAI,OAAO,GAAG,CAAC,CAAC;IAChB,KAAK,IAAI,OAAO,GAAG,MAAM,GAAG,CAAC,EAAE,OAAO,GAAG,CAAC,EAAE,OAAO,KAAK,CAAC;QAAE,EAAE,OAAO,CAAC;IACrE,OAAO,OAAO,CAAC;AACjB,CAAC;AAED,qCAAqC;AACrC,SAAS,WAAW,CAAC,GAAW,EAAE,IAAY;IAC5C,IAAI,MAAM,GAAG,CAAC,CAAC;IACf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,EAAE,CAAC,EAAE,CAAC;QAC9B,MAAM,KAAK,CAAC,CAAC;QACb,MAAM,IAAI,GAAG,GAAG,CAAC,CAAC;QAClB,GAAG,KAAK,CAAC,CAAC;IACZ,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC"}