comparison src/EM.cpp @ 38:5164bccf3064

Return pitch activation matrix

     for (int i = 0; i < m_instruments; ++i) {
         for (int n = 0; n < m_notes; ++n) {
             newSources[i][n] = epsilon;
             if (inRange(i, n)) {
                 float *w = silvet_templates[i].data[n];
+                double pitch = m_pitches[n];
+                double source = m_sources[i][n];
                 for (int j = 0; j < m_bins; ++j) {
-                    newSources[i][n] +=
-                        w[j] * m_q[j] * m_pitches[n] * m_sources[i][n];
+                    newSources[i][n] += w[j] * m_q[j] * pitch * source;
                 }
             }
         }
         normalise(newSources[i]);
     }