comparison src/EM.cpp @ 95:853b2d750688 bqvec

Fix silly bug

     double *norm = allocate<double>(m_binCount);
     v_copy(norm, column, m_binCount);
     normaliseColumn(norm, m_binCount);
     expectation(norm);
     maximisation(norm);
+    deallocate(norm);
 }
 
 const double *
 EM::templateFor(int instrument, int note, int shift)
 {

                 const double source = m_sources[i][n];
                 const double factor = pitch * source * shift;
                 const double *w = templateFor(i, n, f);
 
                 v_copy(contributions, w, m_binCount);
-                v_add(contributions, m_q, m_binCount);
+                v_multiply(contributions, m_q, m_binCount);
                 v_scale(contributions, factor, m_binCount);
 
                 double total = v_sum(contributions, m_binCount);
 
                 if (n >= m_lowestPitch && n <= m_highestPitch) {