]>
git.proxmox.com Git - mirror_edk2.git/blob - StdLib/LibC/gdtoa/ldtoa.c
1 /* $NetBSD: ldtoa.c,v 1.4.2.1 2007/05/07 19:49:06 pavel Exp $ */
4 * Copyright (c) 2003 David Schultz <das@FreeBSD.ORG>
7 * Redistribution and use in source and binary forms, with or without
8 * modification, are permitted provided that the following conditions
10 * 1. Redistributions of source code must retain the above copyright
11 * notice, this list of conditions and the following disclaimer.
12 * 2. Redistributions in binary form must reproduce the above copyright
13 * notice, this list of conditions and the following disclaimer in the
14 * documentation and/or other materials provided with the distribution.
16 * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND
17 * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
18 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
19 * ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
20 * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
21 * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
22 * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
23 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
24 * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
25 * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
28 #include <LibConfig.h>
29 #include <sys/EfiCdefs.h>
36 #include <machine/ieee.h>
40 /* Disable warnings about conversions to narrower data types,
41 primarily for the fpclassify() macro.
43 #pragma warning ( disable : 4244 )
44 // Squelch bogus warnings about uninitialized variable use.
45 #pragma warning ( disable : 4700 )
49 * ldtoa() is a wrapper for gdtoa() that makes it smell like dtoa(),
50 * except that the floating point argument is passed by reference.
51 * When dtoa() is passed a NaN or infinity, it sets expt to 9999.
52 * However, a long double could have a valid exponent of 9999, so we
53 * use INT_MAX in ldtoa() instead.
56 ldtoa(long double *ld
, int mode
, int ndigits
, int *decpt
, int *sign
, char **rve
)
60 LDBL_MANT_DIG
, /* nbits */
61 LDBL_MIN_EXP
- LDBL_MANT_DIG
, /* emin */
62 LDBL_MAX_EXP
- LDBL_MANT_DIG
, /* emax */
63 FPI_Round_near
, /* rounding */
64 #ifdef Sudden_Underflow /* unused, but correct anyway */
73 uint32_t bits
[(LDBL_MANT_DIG
+ 31) / 32];
76 *sign
= (int)(u
.extu_ext
.ext_sign
);
77 be
= (int)(u
.extu_ext
.ext_exp
- (LDBL_MAX_EXP
- 1) - (LDBL_MANT_DIG
- 1));
78 EXT_TO_ARRAY32(u
, bits
);
80 switch (fpclassify(u
.extu_ld
)) {
83 #ifdef LDBL_IMPLICIT_NBIT
84 bits
[LDBL_MANT_DIG
/ 32] |= 1 << ((LDBL_MANT_DIG
- 1) % 32);
85 #endif /* LDBL_IMPLICIT_NBIT */
91 kind
= STRTOG_Denormal
;
92 #ifdef LDBL_IMPLICIT_NBIT
97 kind
= STRTOG_Infinite
;
106 ret
= gdtoa(&fpi
, be
, (ULong
*)bits
, &kind
, mode
, ndigits
, decpt
, rve
);
107 if (*decpt
== -32768)
111 return dtoa((double)*ld
, mode
, ndigits
, decpt
, sign
, rve
);