shithub: opus

Info  •  Files  •  Log  •  Branches

Avoiding problems with x87

logE[b] can end up slightly larger than tonal->highE[b] due to rounding

--- a/src/analysis.c

+++ b/src/analysis.c

@@ -654,7 +654,7 @@

           tonal->lowE[b] = logE[b];

           tonal->highE[b] = MIN32(tonal->lowE[b]+15, tonal->highE[b]);

        }

-       relativeE += (logE[b]-tonal->lowE[b])/(1e-15f + (tonal->highE[b]-tonal->lowE[b]));

+       relativeE += (logE[b]-tonal->lowE[b])/(1e-5f + (tonal->highE[b]-tonal->lowE[b]));

        L1=L2=0;

        for (i=0;i<NB_FRAMES;i++)