Changeset 258195f

Timestamp:

Sep 28, 2009, 8:58:03 PM (15 years ago)

Author:

Paul Brossier <piem@piem.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:

916f8983

Parents:

17b7d66

Message:

src/aubio_priv.h: define VERY_SMALL_NUMBER here, add KILL_DENORMAL and CEIL_DENORMAL macros, SAFE_LOG10 and SAFE_LOG

Location:

src

Files:

• aubio_priv.h (modified) (1 diff)
• temporal/filter.c (modified) (1 diff)

src/aubio_priv.h

@@ -190 +190 @@
 #define ELEM_SWAP(a,b) { register smpl_t t=(a);(a)=(b);(b)=t; }
 
-#define ISDENORMAL(f) f < 1.e-37
+#define VERY_SMALL_NUMBER 2.e-42 //1.e-37
+
+#define IS_DENORMAL(f) f < VERY_SMALL_NUMBER
+
+#define KILL_DENORMAL(f)  IS_DENORMAL(f) ? 0. : f
+#define CEIL_DENORMAL(f)  IS_DENORMAL(f) ? VERY_SMALL_NUMBER : f
+
+#define SAFE_LOG10(f) LOG10(CEIL_DENORMAL(f))
+#define SAFE_LOG(f)   LOG(CEIL_DENORMAL(f))
 
 #define UNUSED __attribute__((unused))

src/temporal/filter.c

@@ -45 +45 @@
     for (j = 0; j < in->length; j++) {
       /* new input */
-      if (ISDENORMAL(in->data[i][j])) {
+      if (IS_DENORMAL(in->data[i][j])) {
         x[0] = y[0] = 0.;
       } else {