79 	long double lo, xn;  in sqrtl()  local
123 	xn = sqrt(u.e);			/* 53-bit estimate of sqrtl(x). */  in sqrtl()
125 	xn = (xn + (u.e / xn)) * 0.5;	/* 106-bit estimate. */  in sqrtl()
129 	lo = (lo - u.e) / xn;		/* Low bits divided by xn. */  in sqrtl()
130 	xn = xn + (u.e / xn);		/* High portion of estimate. */  in sqrtl()
131 	u.e = xn + lo;			/* Combine everything. */  in sqrtl()
137 	xn = x / u.e;			/* Chopped quotient (inexact?). */  in sqrtl()
140 		if (xn == u.e) {  in sqrtl()
145 		xn = dec(xn);		/* xn = xn - ulp. */  in sqrtl()
149 		xn = inc(xn);		/* xn = xn + ulp. */  in sqrtl()
152 		xn = inc(xn);		/* xn  = xn + ulp. */  in sqrtl()
154 	u.e = u.e + xn;				/* Chopped sum. */  in sqrtl()