wavescope-mcp 1.0.2

package/dist/wavelet.d.ts

@@ -0,0 +1,46 @@
+/**
+ * Ricker (Mexican hat) wavelet: ψ(t) = (1 - t²) · exp(-t²/2)
+ */
+export declare function rickerWavelet(t: number): number;
+export interface WaveletCoefficients {
+    scales: number[];
+    coefficients: number[][];
+}
+export interface Peak {
+    position: number;
+    coefficient: number;
+    scale: number;
+}
+export type Boundary = "reflect" | "zero";
+export interface CWTOptions {
+    boundary?: Boundary;
+}
+/**
+ * Compute the Continuous Wavelet Transform (Ricker) over the signal.
+ *
+ * For each scale a, the wavelet ψ_a(t) = (1/√a)·ψ(t/a) is convolved
+ * with the signal. The result for scale a at position b is:
+ *   W(a, b) = Σ_t ψ_a(t-b) · signal[t]
+ *
+ * Boundary handling defaults to symmetric reflection; pass
+ * `{ boundary: "zero" }` for the older zero-pad behavior.
+ */
+export declare function computeCWT(signal: number[], scales?: number[], options?: CWTOptions): WaveletCoefficients;
+/**
+ * Detect local maxima in wavelet coefficient magnitudes across all scales.
+ *
+ * Returns peaks sorted by |coefficient| descending. Each peak is a
+ * local maximum in its scale band — meaning it's larger than its
+ * immediate neighbors at the same scale.
+ *
+ * Plateau handling: >= left, > right — selects the rightmost element
+ * of a flat plateau region.
+ *
+ * Cross-scale ridge collapse: a single structural feature produces local
+ * maxima at the same position across multiple scales. After magnitude
+ * sorting, peaks whose position is within `ridgeWindow` of an already-kept
+ * stronger peak are dropped, so a single spike yields one peak (the
+ * dominant scale) rather than one per scale.
+ */
+export declare function detectPeaks(cwt: WaveletCoefficients, threshold: number, maxPeaks?: number, ridgeWindow?: number): Peak[];
+//# sourceMappingURL=wavelet.d.ts.map

package/dist/wavelet.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"wavelet.d.ts","sourceRoot":"","sources":["../src/wavelet.ts"],"names":[],"mappings":"AAAA;;GAEG;AACH,wBAAgB,aAAa,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAG/C;AAqBD,MAAM,WAAW,mBAAmB;IAClC,MAAM,EAAE,MAAM,EAAE,CAAC;IACjB,YAAY,EAAE,MAAM,EAAE,EAAE,CAAC;CAC1B;AAED,MAAM,WAAW,IAAI;IACnB,QAAQ,EAAE,MAAM,CAAC;IACjB,WAAW,EAAE,MAAM,CAAC;IACpB,KAAK,EAAE,MAAM,CAAC;CACf;AAED,MAAM,MAAM,QAAQ,GAAG,SAAS,GAAG,MAAM,CAAC;AAE1C,MAAM,WAAW,UAAU;IACzB,QAAQ,CAAC,EAAE,QAAQ,CAAC;CACrB;AAiBD;;;;;;;;;GASG;AACH,wBAAgB,UAAU,CACxB,MAAM,EAAE,MAAM,EAAE,EAChB,MAAM,GAAE,MAAM,EAAmB,EACjC,OAAO,GAAE,UAAe,GACvB,mBAAmB,CAmCrB;AAED;;;;;;;;;;;;;;;GAeG;AACH,wBAAgB,WAAW,CACzB,GAAG,EAAE,mBAAmB,EACxB,SAAS,EAAE,MAAM,EACjB,QAAQ,GAAE,MAAY,EACtB,WAAW,GAAE,MAAU,GACtB,IAAI,EAAE,CA2CR"}

package/dist/wavelet.js

@@ -0,0 +1,141 @@
+/**
+ * Ricker (Mexican hat) wavelet: ψ(t) = (1 - t²) · exp(-t²/2)
+ */
+export function rickerWavelet(t) {
+    const t2 = t * t;
+    return (1 - t2) * Math.exp(-t2 / 2);
+}
+/**
+ * Generate wavelet kernel values for scale a, centered at 0.
+ * Truncated to ±5a (covers ~99.99% of the Ricker energy) but bounded
+ * by half the signal length to keep the kernel finite on short inputs.
+ * Includes 1/√a normalization to keep coefficient magnitudes comparable
+ * across scales.
+ */
+function makeKernel(a, numPoints) {
+    if (!Number.isFinite(a) || a <= 0)
+        throw new Error(`Invalid scale: ${a}`);
+    const halfWidth = Math.ceil(5 * a);
+    const half = Math.min(halfWidth, Math.ceil(numPoints / 2));
+    const invSqrtA = 1 / Math.sqrt(a);
+    const kernel = [];
+    for (let t = -half; t <= half; t++) {
+        kernel.push(invSqrtA * rickerWavelet(t / a));
+    }
+    return kernel;
+}
+const DEFAULT_SCALES = [1, 2, 4, 8, 16, 32, 64, 128];
+/**
+ * Reflect-index: mirror out-of-range indices back into [0, N-1].
+ * Used to suppress boundary artifacts where the wavelet's negative
+ * lobes would otherwise be clipped by zero-padding.
+ */
+function reflectIndex(idx, N) {
+    if (N === 1)
+        return 0;
+    const period = 2 * (N - 1);
+    let i = idx % period;
+    if (i < 0)
+        i += period;
+    return i >= N ? period - i : i;
+}
+/**
+ * Compute the Continuous Wavelet Transform (Ricker) over the signal.
+ *
+ * For each scale a, the wavelet ψ_a(t) = (1/√a)·ψ(t/a) is convolved
+ * with the signal. The result for scale a at position b is:
+ *   W(a, b) = Σ_t ψ_a(t-b) · signal[t]
+ *
+ * Boundary handling defaults to symmetric reflection; pass
+ * `{ boundary: "zero" }` for the older zero-pad behavior.
+ */
+export function computeCWT(signal, scales = DEFAULT_SCALES, options = {}) {
+    const boundary = options.boundary ?? "reflect";
+    const N = signal.length;
+    const usedScales = [];
+    for (const a of scales) {
+        if (!usedScales.includes(a))
+            usedScales.push(a);
+    }
+    const coefficients = [];
+    if (N === 0) {
+        return { scales: usedScales, coefficients: usedScales.map(() => []) };
+    }
+    for (const a of usedScales) {
+        const kernel = makeKernel(a, N);
+        const halfKernel = Math.floor(kernel.length / 2);
+        const coeffs = new Array(N);
+        for (let pos = 0; pos < N; pos++) {
+            let sum = 0;
+            for (let k = 0; k < kernel.length; k++) {
+                const signalIdx = pos + k - halfKernel;
+                if (signalIdx >= 0 && signalIdx < N) {
+                    sum += kernel[k] * signal[signalIdx];
+                }
+                else if (boundary === "reflect") {
+                    sum += kernel[k] * signal[reflectIndex(signalIdx, N)];
+                }
+            }
+            coeffs[pos] = sum;
+        }
+        coefficients.push(coeffs);
+    }
+    return { scales: usedScales, coefficients };
+}
+/**
+ * Detect local maxima in wavelet coefficient magnitudes across all scales.
+ *
+ * Returns peaks sorted by |coefficient| descending. Each peak is a
+ * local maximum in its scale band — meaning it's larger than its
+ * immediate neighbors at the same scale.
+ *
+ * Plateau handling: >= left, > right — selects the rightmost element
+ * of a flat plateau region.
+ *
+ * Cross-scale ridge collapse: a single structural feature produces local
+ * maxima at the same position across multiple scales. After magnitude
+ * sorting, peaks whose position is within `ridgeWindow` of an already-kept
+ * stronger peak are dropped, so a single spike yields one peak (the
+ * dominant scale) rather than one per scale.
+ */
+export function detectPeaks(cwt, threshold, maxPeaks = 250, ridgeWindow = 2) {
+    if (cwt.coefficients.length === 0)
+        return [];
+    const peaks = [];
+    for (let si = 0; si < cwt.scales.length; si++) {
+        const scale = cwt.scales[si];
+        const coeffs = cwt.coefficients[si];
+        const N = coeffs.length;
+        for (let pos = 0; pos < N; pos++) {
+            const mag = Math.abs(coeffs[pos]);
+            if (mag < threshold)
+                continue;
+            const leftOk = pos === 0 || mag >= Math.abs(coeffs[pos - 1]);
+            const rightOk = pos === N - 1 || mag > Math.abs(coeffs[pos + 1]);
+            if (leftOk && rightOk) {
+                peaks.push({
+                    position: pos,
+                    coefficient: coeffs[pos],
+                    scale,
+                });
+            }
+        }
+    }
+    peaks.sort((a, b) => Math.abs(b.coefficient) - Math.abs(a.coefficient));
+    const kept = [];
+    for (const peak of peaks) {
+        let overlap = false;
+        for (const k of kept) {
+            if (Math.abs(k.position - peak.position) <= ridgeWindow) {
+                overlap = true;
+                break;
+            }
+        }
+        if (!overlap)
+            kept.push(peak);
+        if (kept.length >= maxPeaks)
+            break;
+    }
+    return kept;
+}
+//# sourceMappingURL=wavelet.js.map

package/dist/wavelet.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"wavelet.js","sourceRoot":"","sources":["../src/wavelet.ts"],"names":[],"mappings":"AAAA;;GAEG;AACH,MAAM,UAAU,aAAa,CAAC,CAAS;IACrC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACjB,OAAO,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;AACtC,CAAC;AAED;;;;;;GAMG;AACH,SAAS,UAAU,CAAC,CAAS,EAAE,SAAiB;IAC9C,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,kBAAkB,CAAC,EAAE,CAAC,CAAC;IAC1E,MAAM,SAAS,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACnC,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,EAAE,IAAI,CAAC,IAAI,CAAC,SAAS,GAAG,CAAC,CAAC,CAAC,CAAC;IAC3D,MAAM,QAAQ,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAClC,MAAM,MAAM,GAAa,EAAE,CAAC;IAC5B,KAAK,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,CAAC,IAAI,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;QACnC,MAAM,CAAC,IAAI,CAAC,QAAQ,GAAG,aAAa,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAC/C,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAmBD,MAAM,cAAc,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC;AAErD;;;;GAIG;AACH,SAAS,YAAY,CAAC,GAAW,EAAE,CAAS;IAC1C,IAAI,CAAC,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC;IACtB,MAAM,MAAM,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC3B,IAAI,CAAC,GAAG,GAAG,GAAG,MAAM,CAAC;IACrB,IAAI,CAAC,GAAG,CAAC;QAAE,CAAC,IAAI,MAAM,CAAC;IACvB,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AACjC,CAAC;AAED;;;;;;;;;GASG;AACH,MAAM,UAAU,UAAU,CACxB,MAAgB,EAChB,SAAmB,cAAc,EACjC,UAAsB,EAAE;IAExB,MAAM,QAAQ,GAAG,OAAO,CAAC,QAAQ,IAAI,SAAS,CAAC;IAC/C,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;IACxB,MAAM,UAAU,GAAa,EAAE,CAAC;IAChC,KAAK,MAAM,CAAC,IAAI,MAAM,EAAE,CAAC;QACvB,IAAI,CAAC,UAAU,CAAC,QAAQ,CAAC,CAAC,CAAC;YAAE,UAAU,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAClD,CAAC;IACD,MAAM,YAAY,GAAe,EAAE,CAAC;IAEpC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QACZ,OAAO,EAAE,MAAM,EAAE,UAAU,EAAE,YAAY,EAAE,UAAU,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;IACxE,CAAC;IAED,KAAK,MAAM,CAAC,IAAI,UAAU,EAAE,CAAC;QAC3B,MAAM,MAAM,GAAG,UAAU,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAChC,MAAM,UAAU,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;QACjD,MAAM,MAAM,GAAG,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC;QAEpC,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;YACjC,IAAI,GAAG,GAAG,CAAC,CAAC;YACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;gBACvC,MAAM,SAAS,GAAG,GAAG,GAAG,CAAC,GAAG,UAAU,CAAC;gBACvC,IAAI,SAAS,IAAI,CAAC,IAAI,SAAS,GAAG,CAAC,EAAE,CAAC;oBACpC,GAAG,IAAI,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,SAAS,CAAC,CAAC;gBACvC,CAAC;qBAAM,IAAI,QAAQ,KAAK,SAAS,EAAE,CAAC;oBAClC,GAAG,IAAI,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,YAAY,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC,CAAC;gBACxD,CAAC;YACH,CAAC;YACD,MAAM,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC;QACpB,CAAC;QAED,YAAY,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;IAC5B,CAAC;IAED,OAAO,EAAE,MAAM,EAAE,UAAU,EAAE,YAAY,EAAE,CAAC;AAC9C,CAAC;AAED;;;;;;;;;;;;;;;GAeG;AACH,MAAM,UAAU,WAAW,CACzB,GAAwB,EACxB,SAAiB,EACjB,WAAmB,GAAG,EACtB,cAAsB,CAAC;IAEvB,IAAI,GAAG,CAAC,YAAY,CAAC,MAAM,KAAK,CAAC;QAAE,OAAO,EAAE,CAAC;IAE7C,MAAM,KAAK,GAAW,EAAE,CAAC;IAEzB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,CAAC,MAAM,CAAC,MAAM,EAAE,EAAE,EAAE,EAAE,CAAC;QAC9C,MAAM,KAAK,GAAG,GAAG,CAAC,MAAM,CAAC,EAAE,CAAC,CAAC;QAC7B,MAAM,MAAM,GAAG,GAAG,CAAC,YAAY,CAAC,EAAE,CAAC,CAAC;QACpC,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;QAExB,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;YACjC,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC;YAClC,IAAI,GAAG,GAAG,SAAS;gBAAE,SAAS;YAE9B,MAAM,MAAM,GAAG,GAAG,KAAK,CAAC,IAAI,GAAG,IAAI,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;YAC7D,MAAM,OAAO,GAAG,GAAG,KAAK,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;YAEjE,IAAI,MAAM,IAAI,OAAO,EAAE,CAAC;gBACtB,KAAK,CAAC,IAAI,CAAC;oBACT,QAAQ,EAAE,GAAG;oBACb,WAAW,EAAE,MAAM,CAAC,GAAG,CAAC;oBACxB,KAAK;iBACN,CAAC,CAAC;YACL,CAAC;QACH,CAAC;IACH,CAAC;IAED,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,WAAW,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,WAAW,CAAC,CAAC,CAAC;IAExE,MAAM,IAAI,GAAW,EAAE,CAAC;IACxB,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACzB,IAAI,OAAO,GAAG,KAAK,CAAC;QACpB,KAAK,MAAM,CAAC,IAAI,IAAI,EAAE,CAAC;YACrB,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,QAAQ,GAAG,IAAI,CAAC,QAAQ,CAAC,IAAI,WAAW,EAAE,CAAC;gBACxD,OAAO,GAAG,IAAI,CAAC;gBACf,MAAM;YACR,CAAC;QACH,CAAC;QACD,IAAI,CAAC,OAAO;YAAE,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QAC9B,IAAI,IAAI,CAAC,MAAM,IAAI,QAAQ;YAAE,MAAM;IACrC,CAAC;IAED,OAAO,IAAI,CAAC;AACd,CAAC"}

package/package.json

@@ -0,0 +1,43 @@
+{
+  "name": "wavescope-mcp",
+  "version": "1.0.2",
+  "description": "Wavelet-based multi-resolution context management for MCP",
+  "type": "module",
+  "main": "./dist/index.js",
+  "bin": {
+    "wavescope-mcp": "./dist/index.js"
+  },
+  "files": [
+    "dist",
+    "README.md",
+    "LICENSE"
+  ],
+  "keywords": [
+    "mcp",
+    "wavelet",
+    "context",
+    "llm"
+  ],
+  "author": "",
+  "license": "ISC",
+  "devDependencies": {
+    "@types/node": "^25.9.1",
+    "@vitest/coverage-v8": "^4.1.7",
+    "tsx": "^4.22.3",
+    "typescript": "^6.0.3",
+    "vitest": "^4.1.7"
+  },
+  "dependencies": {
+    "@modelcontextprotocol/sdk": "^1.29.0",
+    "zod": "^4.4.3"
+  },
+  "scripts": {
+    "build": "tsc",
+    "dev": "tsx src/index.ts",
+    "start": "node dist/index.js",
+    "test": "vitest",
+    "test:run": "vitest run",
+    "test:coverage": "vitest run --coverage",
+    "typecheck": "tsc --noEmit"
+  }
+}