comparison yeti/cqtkernel.yeti @ 6:3caae61fa53f

Sparse etc

 bf = load may.vector.blockfuncs;
 complex = load may.complex;
 window = load may.signal.window;
 fft = load may.transform.fft;
 pl = load may.plot;
-complexmatrix = load may.matrix.complex;
+cm = load may.matrix.complex;
 
 { pow, round, floor, ceil, nextPowerOfTwo } = load may.mathmisc;
 
 fs = 48000;
 

 
         shift = atomOffset + ((i-1) * atomHop);
 
         println "shift = \(shift)";
 
+        shiftedReals = vec.concat [vec.zeros shift, reals];
+        shiftedImags = vec.concat [vec.zeros shift, imags];
+
+//        println "shiftedReals = \(vec.list shiftedReals)";
+//        println "shiftedImags = \(vec.list shiftedImags)";
+
         specKernel = fftFunc
            (complex.complexArray
                (vec.concat [vec.zeros shift, reals])
                (vec.concat [vec.zeros shift, imags]));
 
         map do c:
-            if complex.magnitude c < thresh then complex.zero
-            else complex.scale (1/fftLen) c fi
+            if complex.magnitude c < thresh then complex.zero else c fi
             done specKernel;
 
     done [1..winNr];
 
 done [1..bins];
 
-kmat = complexmatrix.newComplexMatrix (ColumnMajor()) (concat kernels); // or row major?
-println "density = \(complexmatrix.density kmat)";
-complexmatrix.toSparse kmat;
+kmat = cm.toSparse
+   (cm.scaled (1/fftLen)
+       (cm.newComplexMatrix (ColumnMajor()) (concat kernels))); // or row major?
+
+println "density = \(cm.density kmat)";
+
+kmat;