1 #include "libm.h"
2 
3 #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
frexpl(long double x,int * e)4 long 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
frexpl(long double x,int * e)9 long 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