Mercurial > hg > constant-q-cpp

comparison yeti/cqtkernel.yeti @ 72:642df7b3346f

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Support returning a magnitude spectrum (dense) etc
author Chris Cannam <c.cannam@qmul.ac.uk>
date Thu, 20 Mar 2014 16:15:43 +0000
parents 27007f8302f4
children c2e9c91ed3f7
comparison
equal deleted inserted replaced
71:4767faa6726f 72:642df7b3346f
122 // Normalisation 122 // Normalisation
123 123
124 wx1 = vec.maxindex (complex.magnitudes (cm.getRow 0 kmat)); 124 wx1 = vec.maxindex (complex.magnitudes (cm.getRow 0 kmat));
125 wx2 = vec.maxindex (complex.magnitudes (cm.getRow (cm.height kmat - 1) kmat)); 125 wx2 = vec.maxindex (complex.magnitudes (cm.getRow (cm.height kmat - 1) kmat));
126 126
127 subset = cm.columnSlice kmat wx1 (wx2+1); 127 subset = cm.flipped (cm.columnSlice kmat wx1 (wx2+1));
128 square = cm.product (cm.conjugateTransposed subset) subset; 128 square = cm.product (cm.conjugateTransposed subset) subset;
129 129
130 diag = complex.magnitudes (cm.getDiagonal 0 square); 130 diag = complex.magnitudes (cm.getDiagonal 0 square);
131 wK = vec.slice diag (round(1/q)) (vec.length diag - round(1/q) - 2); 131 wK = vec.slice diag (round(1/q)) (vec.length diag - round(1/q) - 2);
132 132