Changeset 09b082d

Timestamp:

Nov 30, 2004, 11:28:53 PM (20 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:

33bf500

Parents:

65f1edc

Message:

src/pitchyin.c

Files:

• ChangeLog (modified) (1 diff)
• src/pitchyin.c (modified) (1 diff)

ChangeLog

@@ -1 +1 @@
 2004-11-39  Paul Brossier <piem@altern.org>
         * configure.ac: added -lmx on macosx
+        * python/aubiocut: seeks for local minima before peak
+        * src/pitchyinc.c: adds draft for all-in-one faster function
 
 2004-10-28  Paul Brossier <piem@altern.org>

src/pitchyin.c

@@ -87 +87 @@
 }
 
+
+/* all the above in one */
+void aubio_pitchyin_getpitchfast(fvec_t * input, fvec_t * yin){
+        uint_t c,j,tau;
+        smpl_t tmp, tmp2;
+        for (c=0;c<input->channels;c++)
+        {
+                for (tau=0;tau<yin->length;tau++)
+                {
+                        yin->data[c][tau] = 0.;
+                }
+                for (tau=1;tau<yin->length;tau++)
+                {
+                        for (j=0;j<yin->length;j++)
+                        {
+                                tmp = input->data[c][j] - input->data[c][j+tau];
+                                yin->data[c][tau] += SQR(tmp);
+                        }
+                }
+                tmp2 = 0.;
+                yin->data[c][0] = 1.;
+                for (tau=1;tau<yin->length;tau++)
+                {
+                        tmp += yin->data[c][tau];
+                        yin->data[c][tau] *= tau/tmp;
+                }
+        }
+        /* should merge the following with above */
+        do 
+        {
+                if(yin->data[c][tau] < 0.1) { 
+                        while (yin->data[c][tau+1] < yin->data[c][tau]) {
+                                tau++;
+                        }
+                        return tau;
+                }
+                tau++;
+        } while (tau<yin->length);
+        return 0;
+}
+