numz 0.1.0 → 0.3.0

package/test/index.js

@@ -0,0 +1,52 @@
+// import { fft, ifft } from "../src/signal/fft/index.js";
+// import { conv, conv2 } from '../src/signal/conv/index.js'
+// import { Complex } from "ziko/math/complex";
+// // import { cos } from 'ziko/math/functions'
+
+// // const x = Array.from({ length: 100 }, (_, j)=> j/10)
+// // const y = cos(x)
+
+// // const ff = fft(x)
+// // // const ff = fft([1,2,3])
+// // // console.log(ff)
+
+// // // console.log(ifft(ff))
+// const signal = [1, 2, 3, 4, 5];
+// let kernel = [1, 0.5];
+// // kernel = [1, new Complex(0, 0.5)];
+// const result = conv(signal, kernel);
+// console.log(
+//    result
+// )
+
+// // console.log(ff)
+// // console.log(ifft(ff).map(n=>+ n.a.toFixed(7)))
+// // console.log(x)
+
+
+// // Example usage:
+// const image = [
+//     [1, 2, 3],
+//     [4, 5, 6],
+//     [7, 8, 9]
+// ];
+
+// const filter = [
+//     [1, 0],
+//     [0, 1]
+// ];
+
+// console.log(conv2(image, filter, 'full', 'fill'));
+
+function hanning(N, ArrayType = Float64Array) {
+  const w = new ArrayType(N);
+  const factor = 2 * Math.PI / N;
+
+  for (let n = 0; n < N; n++) {
+    w[n] = 0.5 - 0.5 * Math.cos(factor * n);
+  }
+  return w;
+}
+
+console.log(hanning(10).map(a=>a*20))
+

package/test/index.py

@@ -0,0 +1,14 @@
+from scipy import signal
+image = [
+    [1, 2, 3],
+    [4, 5, 6],
+    [7, 8, 9]
+]
+
+kernel = [
+    [1, 0],
+    [0, 1]
+]
+print(
+    signal.convolve2d(image, kernel, 'full')
+)

package/src/signal/fft/fft1d/index.js

@@ -1,18 +0,0 @@
-import { twiddle } from "../utils/index.js";
-export function fft1d(arr) {
-    const N = arr.length;
-    if (N === 1) return [arr[0]];
-
-    // Split even/odd
-    const even = fft1d(arr.filter((_, i) => i % 2 === 0));
-    const odd  = fft1d(arr.filter((_, i) => i % 2 === 1));
-
-    const out = new Array(N);
-    for (let k = 0; k < N/2; k++) {
-        const w = twiddle(N, k);
-        const t = w.mul(odd[k]);
-        out[k]        = even[k].add(t);
-        out[k + N/2]  = even[k].sub(t);
-    }
-    return out;
-}

package/src/signal/fft/fftnd/index.js

@@ -1,27 +0,0 @@
-import { normalizeComplexND } from '../utils/index.js'
-import { fft1d } from '../fft1d/index.js';
-function fftND(arr) {
-    arr = normalizeComplexND(arr);
-
-    // Base case: 1-D
-    if (!Array.isArray(arr[0])) 
-        return fft1d(arr);
-
-    const transformed = arr.map(sub => fftND(sub));
-
-    const rows = transformed.length;
-    const cols = transformed[0].length;
-
-    const colsData = Array.from({ length: cols }, (_, c) =>
-        Array.from({ length: rows }, (_, r) => transformed[r][c])
-    );
-
-    const colsTransformed = colsData.map(col => fft1d(col));
-
-    // Transpose back
-    const out = Array.from({ length: rows }, (_, r) =>
-        Array.from({ length: cols }, (_, c) => colsTransformed[c][r])
-    );
-
-    return out;
-}