5 	union {double f; uint64_t i;} u = {x};  in modf()  local
7 	int e = (int)(u.i>>52 & 0x7ff) - 0x3ff;  in modf()
12 		if (e == 0x400 && u.i<<12 != 0) /* nan */  in modf()
14 		u.i &= 1ULL<<63;  in modf()
15 		return u.f;  in modf()
20 		u.i &= 1ULL<<63;  in modf()
21 		*iptr = u.f;  in modf()
26 	if ((u.i & mask) == 0) {  in modf()
28 		u.i &= 1ULL<<63;  in modf()
29 		return u.f;  in modf()
31 	u.i &= ~mask;  in modf()
32 	*iptr = u.f;  in modf()
33 	return x - u.f;  in modf()