Mercurial > hg > silvet

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

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix silly bug
author Chris Cannam
date Tue, 06 May 2014 14:37:42 +0100
parents 891cbcf1e4d2
children 9ecad4c9c2a2
line wrap: on
line diff
--- a/src/EM.cpp	Tue May 06 14:29:07 2014 +0100
+++ b/src/EM.cpp	Tue May 06 14:37:42 2014 +0100
@@ -123,6 +123,7 @@
     normaliseColumn(norm, m_binCount);
     expectation(norm);
     maximisation(norm);
+    deallocate(norm);
 }
 
 const double *
@@ -193,7 +194,7 @@
                 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);