1#include "libm.h"
2
3#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
4long double frexpl(long double x, int *e)
5{
6	return frexp(x, e);
7}
8#elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384
9long double frexpl(long double x, int *e)
10{
11	union ldshape u = {x};
12	int ee = u.i.se & 0x7fff;
13
14	if (!ee) {
15		if (x) {
16			x = frexpl(x*0x1p120, e);
17			*e -= 120;
18		} else *e = 0;
19		return x;
20	} else if (ee == 0x7fff) {
21		return x;
22	}
23
24	*e = ee - 0x3ffe;
25	u.i.se &= 0x8000;
26	u.i.se |= 0x3ffe;
27	return u.f;
28}
29#endif
30