shithub: aubio

Info  •  Files  •  Log  •  Branches

src/spectral/specdesc.c: add more noise to make sure log doesn't explode when fftnorm is very small and old fftnorm is null

--- a/src/spectral/specdesc.c

+++ b/src/spectral/specdesc.c

@@ -193,7 +193,7 @@

     onset->data[0] = 0.;

     for (j=0;j<fftgrain->length;j++) {

       onset->data[0] += fftgrain->norm[j]

-        *LOG(1.+fftgrain->norm[j]/(o->oldmag->data[j]+1.e-10));

+        *LOG(1.+fftgrain->norm[j]/(o->oldmag->data[j]+1.e-1));

       o->oldmag->data[j] = fftgrain->norm[j];

     }

     if (isnan(onset->data[0])) onset->data[0] = 0.;

@@ -206,7 +206,7 @@

   uint_t j;

     onset->data[0] = 0.;

     for (j=0;j<fftgrain->length;j++) {

-      onset->data[0] += LOG(1.+fftgrain->norm[j]/(o->oldmag->data[j]+1.e-10));

+      onset->data[0] += LOG(1.+fftgrain->norm[j]/(o->oldmag->data[j]+1.e-1));

       o->oldmag->data[j] = fftgrain->norm[j];

     }

     if (isnan(onset->data[0])) onset->data[0] = 0.;