Changeset d8604ac for src

Timestamp:

Mar 16, 2006, 4:35:51 PM (18 years ago)

Author:

Paul Brossier <piem@altern.org>

Branches:

feature/autosink, feature/cnn, feature/cnn_org, feature/constantq, feature/crepe, feature/crepe_org, feature/pitchshift, feature/pydocstrings, feature/timestretch, fix/ffmpeg5, master, pitchshift, sampler, timestretch, yinfft+

Children:

8595fe3

Parents:

5d1c070

Message:

activate adapt thresh, reduce numb partials for high frequencies, add support for samplerate
activate adapt thresh, reduce numb partials for high frequencies, add support for samplerate

Location:

src

Files:

• pitchmcomb.c (modified) (7 diffs)
• pitchmcomb.h (modified) (1 diff)

src/pitchmcomb.c

@@ -94 +94 @@
     newmag->data[i][j]=fftgrain->norm[i][j];
   /* detect only if local energy > 10. */ 
-  if (vec_local_energy(newmag)>10.) {
+  //if (vec_local_energy(newmag)>10.) {
     //hfc = vec_local_hfc(newmag); //not used
     aubio_pitchmcomb_spectral_pp(p, newmag);
     aubio_pitchmcomb_combdet(p,newmag);
+    //aubio_pitchmcomb_sort_cand_freq(p->candidates,p->ncand);
     // print picked candidate
     //AUBIO_DBG("%f\n",p->candidates[p->goodcandidate]->ebin);
     //AUBIO_DBG("%f\n",p->candidates[p->goodcandidate]->ebin);
     return p->candidates[p->goodcandidate]->ebin;
-  } else {
+  /*} else {
     return -1.;
-  }
+  }*/
 }
 
@@ -139 +140 @@
 void aubio_pitchmcomb_spectral_pp(aubio_pitchmcomb_t * p, fvec_t * newmag) {
   fvec_t * mag = (fvec_t *)p->scratch;
-  //fvec_t * tmp = (fvec_t *)p->scratch2;
+  fvec_t * tmp = (fvec_t *)p->scratch2;
   uint_t i=0,j;
   uint_t length = mag->length;
@@ -150 +151 @@
   /* skipped */                      /* low pass filtering   */
   /** \bug: vec_movind_thres writes out of bounds */
-  //vec_adapt_thres(mag,tmp,p->win_post,p->win_pre); /* adaptative threshold */
+  vec_adapt_thres(mag,tmp,p->win_post,p->win_pre); /* adaptative threshold */
   vec_add(mag,-p->threshold);        /* fixed threshold      */
   {
@@ -192 +193 @@
   /* get the biggest peak in the spectrum */
   root_peak = aubio_pitchmcomb_get_root_peak(peaks,count);
+  if (peaks[root_peak].ebin > aubio_miditofreq(80.)/p->tau) N = 3;
+  if (peaks[root_peak].ebin > aubio_miditofreq(85.)/p->tau) N = 2;
+  if (peaks[root_peak].ebin > aubio_miditofreq(90.)/p->tau) N = 1;
   /* now calculate the energy of each of the 5 combs */
   for (l=0;l<M;l++) {
@@ -208 +212 @@
     /* for each in candidate[l]->ecomb[k] */
     for (k=0;k<curlen;k++) {
-      xx = 10000.;
+      xx = 100000.;
       /** get the candidate->ecomb the closer to peaks.ebin 
        * (to cope with the inharmonicity)*/
@@ -255 +259 @@
         count += ispeak;
         spectral_peaks[count-1].bin = j;
-        spectral_peaks[count-1].ebin = vec_quadint(X,j);
+        spectral_peaks[count-1].ebin = vec_quadint(X,j) - 1.;
       }
     }
@@ -307 +311 @@
 }
 
-aubio_pitchmcomb_t * new_aubio_pitchmcomb(uint_t size, uint_t channels) {
+aubio_pitchmcomb_t * new_aubio_pitchmcomb(uint_t size, uint_t channels, uint_t samplerate) {
   aubio_pitchmcomb_t * p = AUBIO_NEW(aubio_pitchmcomb_t);
   /** \bug should check if size / 8 > post+pre+1 */
   uint_t i;
   uint_t spec_size;
-  p->spec_partition   = 2;
+  p->spec_partition   = 4;
   p->ncand            = 5;
-  p->npartials        = 10;
+  p->npartials        = 5;
   p->cutoff           = 1.;
   p->threshold        = 0.01;
   p->win_post         = 8;
   p->win_pre          = 7;
-  p->tau              = 44100./size;
+  p->tau              = samplerate/size;
   p->alpha            = 9.;
   p->goodcandidate    = 0;

src/pitchmcomb.h

@@ -35 +35 @@
 smpl_t aubio_pitchmcomb_detect(aubio_pitchmcomb_t * p, cvec_t * fftgrain);
 uint_t aubio_pitch_cands(aubio_pitchmcomb_t * p, cvec_t * fftgrain, smpl_t * cands);
-aubio_pitchmcomb_t * new_aubio_pitchmcomb(uint_t size, uint_t channels);
+aubio_pitchmcomb_t * new_aubio_pitchmcomb(uint_t size, uint_t channels, uint_t samplerate);
 void del_aubio_pitchmcomb(aubio_pitchmcomb_t *p);