diff options
| author | Damien George <damien.p.george@gmail.com> | 2014-12-18 14:44:02 +0000 |
|---|---|---|
| committer | Damien George <damien.p.george@gmail.com> | 2014-12-18 14:44:02 +0000 |
| commit | f04329e93b7bfa5462f232832d64aac3dc52d5db (patch) | |
| tree | 96b06142670e0bb1ddfddcd353872f4612954bf8 /lib/libm/acoshf.c | |
| parent | 6936f4626c8748c077ad2c63827d4deeb43b867d (diff) | |
| download | micropython-f04329e93b7bfa5462f232832d64aac3dc52d5db.tar.gz micropython-f04329e93b7bfa5462f232832d64aac3dc52d5db.zip | |
lib/libm: Add acosh, asinh, atanh, tan; get working with stmhal.
acoshf, asinhf, atanhf were added from musl. mathsincos.c was
split up into its original, separate files (from newlibe-nano-2).
tan was added.
All of the important missing float functions are now implemented,
and pyboard now passes tests/float/math_fun.py (finally!).
Diffstat (limited to 'lib/libm/acoshf.c')
| -rw-r--r-- | lib/libm/acoshf.c | 32 |
1 files changed, 32 insertions, 0 deletions
diff --git a/lib/libm/acoshf.c b/lib/libm/acoshf.c new file mode 100644 index 0000000000..8a8409f583 --- /dev/null +++ b/lib/libm/acoshf.c @@ -0,0 +1,32 @@ +/*****************************************************************************/ +/*****************************************************************************/ +// acoshf from musl-0.9.15 +/*****************************************************************************/ +/*****************************************************************************/ + +#include "libm.h" + +#if FLT_EVAL_METHOD==2 +#undef sqrtf +#define sqrtf sqrtl +#elif FLT_EVAL_METHOD==1 +#undef sqrtf +#define sqrtf sqrt +#endif + +/* acosh(x) = log(x + sqrt(x*x-1)) */ +float acoshf(float x) +{ + union {float f; uint32_t i;} u = {x}; + uint32_t a = u.i & 0x7fffffff; + + if (a < 0x3f800000+(1<<23)) + /* |x| < 2, invalid if x < 1 or nan */ + /* up to 2ulp error in [1,1.125] */ + return log1pf(x-1 + sqrtf((x-1)*(x-1)+2*(x-1))); + if (a < 0x3f800000+(12<<23)) + /* |x| < 0x1p12 */ + return logf(2*x - 1/(x+sqrtf(x*x-1))); + /* x >= 0x1p12 */ + return logf(x) + 0.693147180559945309417232121458176568f; +} |