| Index: fusl/src/math/tanl.c
|
| diff --git a/fusl/src/math/tanl.c b/fusl/src/math/tanl.c
|
| new file mode 100644
|
| index 0000000000000000000000000000000000000000..6af067127ade990b98d26696594e5737fbf50ade
|
| --- /dev/null
|
| +++ b/fusl/src/math/tanl.c
|
| @@ -0,0 +1,29 @@
|
| +#include "libm.h"
|
| +
|
| +#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
|
| +long double tanl(long double x)
|
| +{
|
| + return tan(x);
|
| +}
|
| +#elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384
|
| +long double tanl(long double x)
|
| +{
|
| + union ldshape u = {x};
|
| + long double y[2];
|
| + unsigned n;
|
| +
|
| + u.i.se &= 0x7fff;
|
| + if (u.i.se == 0x7fff)
|
| + return x - x;
|
| + if (u.f < M_PI_4) {
|
| + if (u.i.se < 0x3fff - LDBL_MANT_DIG/2) {
|
| + /* raise inexact if x!=0 and underflow if subnormal */
|
| + FORCE_EVAL(u.i.se == 0 ? x*0x1p-120f : x+0x1p120f);
|
| + return x;
|
| + }
|
| + return __tanl(x, 0, 0);
|
| + }
|
| + n = __rem_pio2l(x, y);
|
| + return __tanl(y[0], y[1], n&1);
|
| +}
|
| +#endif
|
|
|