949f079c9c027cc6b57d3cb3e65bc5040f870aa0
[Faustine.git] / dsp_files / fft_pattern_matching.dsp
1 import ( "complex.lib" ) ;
2 import ( "fft.lib.dsp" ) ;
3
4 concats = case {
5 (1) => _, 1 : vectorize ;
6 (m) => concats(m - 1), ( _, 1 : vectorize) : #;
7 };
8
9 //subvect_even(n) = par(i, (n/2), (_, (i, 2 : *) : [])) : concats(n/2);
10 //subvect_odd(n) = par(i, (n/2), (_, (i, 2 : * : _, 1 : +) : [])) : concats(n/2);
11
12 evens = case {
13 (2) => [0] : vectorize(1);
14 (n) => evens(n-2) # ((_, (n-2) : []) : vectorize(1));
15 };
16
17 odds = case {
18 (2) => [1] : vectorize(1);
19 (n) => odds(n - 2) # ((_, (n-1) : []) : vectorize(1));
20 };
21
22 subvect_even(n) = _ <: evens(n);
23 subvect_odd(n) = _ <: odds(n);
24
25 eo(n) = _ <: subvect_even(n), subvect_odd(n);
26
27 output = serialize, serialize : vectorize(1), vectorize(1) : #;
28
29 //process = (1,2,3,4,5,6,7,8) : evens(8);
30
31 process = _ : vectorize(8) : eo(8) : output;