33 	union IEEEl2bits u;  in inc()  local
35 	u.e = x;  in inc()
36 	if (++u.bits.manl == 0) {  in inc()
37 		if (++u.bits.manh == 0) {  in inc()
38 			u.bits.exp++;  in inc()
39 			u.bits.manh |= LDBL_NBIT;  in inc()
42 	return (u.e);  in inc()
49 	union IEEEl2bits u;  in dec()  local
51 	u.e = x;  in dec()
52 	if (u.bits.manl-- == 0) {  in dec()
53 		if (u.bits.manh-- == LDBL_NBIT) {  in dec()
54 			u.bits.exp--;  in dec()
55 			u.bits.manh |= LDBL_NBIT;  in dec()
58 	return (u.e);  in dec()
77 	union IEEEl2bits u;  in sqrtl()  local
82 	u.e = x;  in sqrtl()
87 	if (u.bits.exp == LDBL_INF_NAN_EXP) {  in sqrtl()
88                 if (u.bits.sign && u.bits.manh == LDBL_NBIT_INF && u.bits.manl == 0)  in sqrtl()
94 	if ((u.bits.manh | u.bits.manl | u.bits.exp) == 0)  in sqrtl()
98 	if (u.bits.sign)  in sqrtl()
103 	if (u.bits.exp == 0) {  in sqrtl()
105 		u.e *= 0x1.0p514L;  in sqrtl()
114 	if ((u.bits.exp - (BIAS-1)) & 1) {	/* n is odd.     */  in sqrtl()
115 		k += u.bits.exp - BIAS;	/* 2k = n - 1.   */  in sqrtl()
117                 u.dbits.dl = scalbn(u.dbits.dl, BIAS - u.bits.exp);  in sqrtl()
119 		u.bits.exp = BIAS;		/* u.e in [1,2). */  in sqrtl()
121                 k += u.bits.exp - (BIAS + 1);	/* 2k = n - 2.   */  in sqrtl()
123                 u.dbits.dl = scalbn(u.dbits.dl, (BIAS + 1) - u.bits.exp);  in sqrtl()
125 		u.bits.exp = (BIAS + 1);	/* u.e in [2,4). */  in sqrtl()
132 	xn = (long double)sqrt((double)u.e);			/* 53-bit estimate of sqrtl(x). */  in sqrtl()
133 	xn = (xn + (u.e / xn)) * 0.5L;	/* 106-bit estimate. */  in sqrtl()
135 	lo = u.e;  in sqrtl()
137         u.dbits.dl = 0.0;               /* Zero out lower double */  in sqrtl()
139 	u.bits.manl = 0;		/* Zero out lower bits. */  in sqrtl()
141 	lo = (lo - u.e) / xn;		/* Low bits divided by xn. */  in sqrtl()
142 	xn = xn + (u.e / xn);		/* High portion of estimate. */  in sqrtl()
143 	u.e = xn + lo;			/* Combine everything. */  in sqrtl()
145 	u.bits.exp += (k >> 1) - 1;  in sqrtl()
147         u.dbits.dl = scalbn(u.dbits.dl, (k>>1) -1);  in sqrtl()
156                 xn = x / u.e;			/* Chopped quotient (inexact?). */  in sqrtl()
159                         if (xn == u.e) {  in sqrtl()
161                                 return (u.e);  in sqrtl()
170                         u.e = inc(u.e);		/* u.e = u.e + ulp. */  in sqrtl()
173                 u.e = u.e + xn;				/* Chopped sum. */  in sqrtl()
175                 u.bits.exp--;  in sqrtl()
178 	return (u.e);  in sqrtl()