* Copyright (c) 1992 The Regents of the University of California.
* This software was developed by the Computer Systems Engineering group
* at Lawrence Berkeley Laboratory under DARPA contract BG 91-66 and
* contributed to Berkeley.
* %sccs.include.redist.c%
#if defined(LIBC_SCCS) && !defined(lint)
static char sccsid
[] = "@(#)fixunsdfdi.c 5.4 (Berkeley) %G%";
#endif /* LIBC_SCCS and not lint */
#define ONE_FOURTH (1 << (LONG_BITS - 2))
#define ONE_HALF (ONE_FOURTH * 2.0)
#define ONE (ONE_FOURTH * 4.0)
* Convert double to (unsigned) quad.
* Not sure what to do with negative numbers---for now, anything out
* of range becomes UQUAD_MAX.
return (UQUAD_MAX
); /* ??? should be 0? ERANGE??? */
* Get the upper part of the result. Note that the divide
* may round up; we want to avoid this if possible, so we
toppart
= (x
- ONE_HALF
) / ONE
;
* Now build a u_quad_t out of the top part. The difference
* between x and this is the bottom part (this may introduce
* a few fuzzy bits, but what the heck). With any luck this
* difference will be nonnegative: x should wind up in the
* range [0..ULONG_MAX]. For paranoia, we assume [LONG_MIN..
t
.ul
[H
] = (unsigned long)toppart
;