comparison src/EM.cpp @ 86:6075e92d63ab timing

Reorder internal loops

                 newSources[i][n] = epsilon;
 
                 const float *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) {
+                        newPitches[n] += w[j] * m_q[j] * factor;
+                    }
+
+                    if (inRange(i, n)) {
+                        for (int j = 0; j < m_binCount; ++j) {
+                            newSources[i][n] += w[j] * m_q[j] * factor;
+                        }
+                    }
+                }
+
                 for (int j = 0; j < m_binCount; ++j) {
-
-                    const double contribution = w[j] * m_q[j] * factor;
-
-                    if (n >= m_lowestPitch && n <= m_highestPitch) {
-                        newPitches[n] += contribution;
-                    }
-
-                    newShifts[f][n] += contribution;
-        
-                    if (inRange(i, n)) {
-                        newSources[i][n] += contribution;
-                    }
+                    newShifts[f][n] += w[j] * m_q[j] * factor;
                 }
             }
         }
     }