diff options
| -rwxr-xr-x | configure | 2 | ||||
| -rw-r--r-- | libavutil/libm.h | 19 |
2 files changed, 21 insertions, 0 deletions
@@ -1815,6 +1815,8 @@ MATH_FUNCS=" copysign cosf erf + exp10 + exp10f exp2 exp2f expf diff --git a/libavutil/libm.h b/libavutil/libm.h index 146768a..6f9ac1b 100644 --- a/libavutil/libm.h +++ b/libavutil/libm.h @@ -29,6 +29,7 @@ #include "config.h" #include "attributes.h" #include "intfloat.h" +#include "mathematics.h" #if HAVE_MIPSFPU && HAVE_INLINE_ASM #include "libavutil/mips/libm_mips.h" @@ -292,6 +293,24 @@ static inline double erf(double z) #define exp2f(x) ((float)exp2(x)) #endif /* HAVE_EXP2F */ +/* Somewhat inaccurate fallbacks, relative error ~ 1e-13 concentrated on very +small and very large values. For perfection accuracy-wise, should use pow. +Speed benefits (>2x average, with no super slow paths) deemed to be worth the +accuracy tradeoff */ +#if !HAVE_EXP10 +static av_always_inline double exp10(double x) +{ + return exp2(M_LOG2_10 * x); +} +#endif /* HAVE_EXP10 */ + +#if !HAVE_EXP10F +static av_always_inline float exp10f(float x) +{ + return exp2f(M_LOG2_10 * x); +} +#endif /* HAVE_EXP10F */ + #if !HAVE_ISINF #undef isinf /* Note: these do not follow the BSD/Apple/GNU convention of returning -1 for |
