comparison src/EM.cpp @ 88:97b77e7cb94c timing

Store templates as doubles instead of floats

     normaliseColumn(column);
     expectation(column);
     maximisation(column);
 }
 
-const float *
+const double *
 EM::templateFor(int instrument, int note, int shift)
 {
     return silvet_templates[instrument].data[note] + shift;
 }
 

     for (int i = 0; i < m_instrumentCount; ++i) {
         for (int n = 0; n < m_noteCount; ++n) {
             const double pitch = m_pitches[n];
             const double source = m_sources[i][n];
             for (int f = 0; f < m_shiftCount; ++f) {
-                const float *w = templateFor(i, n, f);
+                const double *w = templateFor(i, n, f);
                 const double shift = m_shifts[f][n];
                 const double factor = pitch * source * shift;
                 for (int j = 0; j < m_binCount; ++j) {
                     m_estimate[j] += w[j] * factor;
                 }

             for (int i = 0; i < m_instrumentCount; ++i) {
 
                 const double source = m_sources[i][n];
                 newSources[i][n] = epsilon;
 
-                const float *w = templateFor(i, n, f);
+                const double *w = templateFor(i, n, f);
                 const double factor = pitch * source * shift;
 
                 if (n >= m_lowestPitch && n <= m_highestPitch) {
 
                     for (int j = 0; j < m_binCount; ++j) {