comparison src/EM.cpp @ 164:629c9525b815

Some measurements suggest we get slightly better results with more iterations, but not drastically worse with fewer -- so go lower for draft and higher for hq mode

                 v_add_with_gain(m_estimate, w, factor, m_binCount);
             }
         }
     }
 
-    //!!! try normalising so as to sum to unity and then taking L2
-    //!!! norm of the two (for measure of how close they are)
-
     for (int i = 0; i < m_binCount; ++i) {
         m_q[i] = column[i] / m_estimate[i];
     }
+
+/*
+    double l2norm = 0.0;
+    
+    for (int i = 0; i < m_binCount; ++i) {
+        l2norm += (column[i] - m_estimate[i]) * (column[i] - m_estimate[i]);
+    }
+
+    l2norm = sqrt(l2norm);
+    cerr << "l2norm = " << l2norm << endl;
+*/
 }
 
 void
 EM::maximisation(const float *column)
 {