Changeset 422a2a9

Timestamp:

Feb 2, 2017, 2:29:43 AM (8 years ago)

Author:

Paul Brossier <piem@piem.org>

Branches:

yinfft+

Children:

5bdbb83

Parents:

09558b8

Message:

src/pitch/pitchyinfft.c: lookahead only, stop at valley low, improve comments

File:

• src/pitch/pitchyinfft.c (modified) (3 diffs)

src/pitch/pitchyinfft.c

@@ -115 +115 @@
   fvec_t *yin = p->yinfft;
   smpl_t tmp = 0., sum = 0.;
-  uint_t startbin = 0, endbin = 0, before = 0, after = 0;
+  // empirically derived peak width to look for
+  uint_t lookafter = yin->length / 128, endbin;
 
   // window the input
@@ -159 +160 @@
   }
 
-  // choose lowest confident candidate first, to avoid choosing harmonics
+  // choose first confident candidate, to avoid choosing lower harmonics
   tau = 0;
   for (l = 1; l < yin->length; l++) {
     // is this candidate "roughly" as good as the lowest one?
+    // the constant 0.1 is empirically derived
     if (ABS (yin->data[l] - tmp) < 0.1) {
       tau = l;
@@ -168 +170 @@
     }
   }
-  // find local min around current pick to sharpen the results
-  before = 2, after = 5;
-  // for low frequency, the peak may be large, we use a larger after
-  if (tau > yin->length/2) { before = 1; after = 20; }
-  startbin = tau > before ? tau - before : 0;
-  endbin = tau < yin->length - 1 - after ? tau + after : yin->length - 1;
+  // find local min around current peak to sharpen the results
+  endbin = tau + lookafter < yin->length - 1 ? tau + lookafter : yin->length - 1;
   tmp = yin->data[tau];
-  for (l = startbin; l < endbin; l++) {
+  for (l = tau; l < endbin; l++) {
     if (yin->data[l] < tmp ) {
       tmp = yin->data[l];
       tau = l;
+    }
+    // stop as soon as we start going up again
+    if (yin->data[l] > tmp && l > tau) {
+      break;
     }
   }