# HG changeset patch
# User Chris Cannam <c.cannam@qmul.ac.uk>
# Date 1398789172 -3600
# Node ID 31a1271596c495a86192a65a872fb142446ef183
# Parent  24ff66e793fd2955a8d0ace6779f7aa712bd93e8
Another step in icqt

diff -r 24ff66e793fd -r 31a1271596c4 yeti/cqt.yeti
--- a/yeti/cqt.yeti	Tue Apr 29 17:27:04 2014 +0100
+++ b/yeti/cqt.yeti	Tue Apr 29 17:32:52 2014 +0100
@@ -62,7 +62,8 @@
 
     streams = manipulate.duplicated octaves str;
 
-    //!!! can't be right!
+    // forward transform uses the conjugate-transposed kernel, inverse
+    // uses the original
     kernel = cm.transposed (cm.conjugateTransposed kdata.kernel);
 
 //    eprintln "have kernel";
diff -r 24ff66e793fd -r 31a1271596c4 yeti/icqt.yeti
--- a/yeti/icqt.yeti	Tue Apr 29 17:27:04 2014 +0100
+++ b/yeti/icqt.yeti	Tue Apr 29 17:32:52 2014 +0100
@@ -105,6 +105,17 @@
         map pile (list (flattenOctaves (map decomposeOctaves cq.cqComplex)));
 
     for octaveColumnLists do l: println "octave column list length: \(length l)" done;
+
+    kernel = cm.transposed kdata.kernel; // right way around for the multiply
+    
+    spectra = 
+        map do l:
+            map do col:
+                cm.product kernel (cm.newComplexColumnVector col);
+            done l;
+        done octaveColumnLists;
+
+    eprintln "calculated spectra, now to ifft, overlap-add..."
     
 );