36 	u.e = x;  in inc()
43 	return (u.e);  in inc()
52 	u.e = x;  in dec()
59 	return (u.e);  in dec()
83 	u.e = x;  in sqrtl()
106 		u.e *= 0x1.0p514L;  in sqrtl()
133 	xn = (long double)sqrt((double)u.e);			/* 53-bit estimate of sqrtl(x). */  in sqrtl()
134 	xn = (xn + (u.e / xn)) * 0.5L;	/* 106-bit estimate. */  in sqrtl()
136 	lo = u.e;  in sqrtl()
142 	lo = (lo - u.e) / xn;		/* Low bits divided by xn. */  in sqrtl()
143 	xn = xn + (u.e / xn);		/* High portion of estimate. */  in sqrtl()
144 	u.e = xn + lo;			/* Combine everything. */  in sqrtl()
157                 xn = x / u.e;			/* Chopped quotient (inexact?). */  in sqrtl()
160                         if (xn == u.e) {  in sqrtl()
162                                 return (u.e);  in sqrtl()
171                         u.e = inc(u.e);		/* u.e = u.e + ulp. */  in sqrtl()
174                 u.e = u.e + xn;				/* Chopped sum. */  in sqrtl()
179 	return (u.e);  in sqrtl()