36 	double dr, dt, dx;  in cbrtl()  local
95 	dt = ft;  in cbrtl()
96 	dr = dt * dt * dt;  in cbrtl()
97 	dt = dt * (dx + dx + dr) / (dx + dr + dr);  in cbrtl()
100 	dr = dt * dt * dt;  in cbrtl()
101 	dt = dt * (dx + dx + dr) / (dx + dr + dr);  in cbrtl()
105 	 * dt is cbrtl(x) to ~47 bits (after x has been reduced to 1 <= x < 8).  in cbrtl()
113 	t = dt + vd - 0x1.0p32;  in cbrtl()
116 	 * Round dt away from zero to 47 bits.  Since we don't trust the 47,  in cbrtl()
118 	 * might be avoidable in this case, since on most machines dt will  in cbrtl()
121 	 * dt is much more accurate than needed.  in cbrtl()
123 	t = dt + 0x2.0p-46 + 0x1.0p60L - 0x1.0p60;  in cbrtl()