comparison src/Silvet.cpp @ 305:04a3c152e590

Now that we are returning the pitch activation matrix, we don't want to skip silent columns in the postprocess loop

     if (wantShifts) {
         localBestShifts = 
             vector<vector<int> >(width, vector<int>(pack.templateNoteCount, 0));
     }
 
-    vector<bool> present(width, false);
-
+    double columnThreshold = 1e-5;
+    
 #pragma omp parallel for
     for (int i = 0; i < width; ++i) {
 
         double sum = 0.0;
         for (int j = 0; j < pack.templateHeight; ++j) {
             sum += filtered.at(i).at(j);
         }
-        if (sum < 1e-5) continue;
-
-        present[i] = true;
+        if (sum < columnThreshold) continue;
 
         EM em(&pack, m_hqMode);
 
         em.setPitchSparsity(pack.pitchSparsity);
         em.setSourceSparsity(pack.sourceSparsity);

             }                
         }
     }
         
     for (int i = 0; i < width; ++i) {
-
-        if (!present[i]) {
-            // silent column
-            for (int j = 0; j < pack.templateNoteCount; ++j) {
-                m_postFilter[j]->push(0.0);
-            }
-            m_pianoRoll.push_back(map<int, double>());
-            if (wantShifts) {
-                m_pianoRollShifts.push_back(map<int, int>());
-            }
-            continue;
-        }
 
         vector<double> filtered = postProcess
             (localPitches[i], localBestShifts[i], wantShifts);
 
         Feature f;