]>
Commit | Line | Data |
---|---|---|
2aa62f2b | 1 | /* $NetBSD: ldtoa.c,v 1.4.2.1 2007/05/07 19:49:06 pavel Exp $ */\r |
2 | \r | |
3 | /*-\r | |
4 | * Copyright (c) 2003 David Schultz <das@FreeBSD.ORG>\r | |
5 | * All rights reserved.\r | |
6 | *\r | |
7 | * Redistribution and use in source and binary forms, with or without\r | |
8 | * modification, are permitted provided that the following conditions\r | |
9 | * are met:\r | |
10 | * 1. Redistributions of source code must retain the above copyright\r | |
11 | * notice, this list of conditions and the following disclaimer.\r | |
12 | * 2. Redistributions in binary form must reproduce the above copyright\r | |
13 | * notice, this list of conditions and the following disclaimer in the\r | |
14 | * documentation and/or other materials provided with the distribution.\r | |
15 | *\r | |
16 | * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND\r | |
17 | * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE\r | |
18 | * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE\r | |
19 | * ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE\r | |
20 | * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL\r | |
21 | * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS\r | |
22 | * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)\r | |
23 | * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT\r | |
24 | * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY\r | |
25 | * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF\r | |
26 | * SUCH DAMAGE.\r | |
27 | */\r | |
28 | #include <LibConfig.h>\r | |
29 | #include <sys/EfiCdefs.h>\r | |
30 | \r | |
31 | #include <float.h>\r | |
32 | #include <inttypes.h>\r | |
33 | #include <limits.h>\r | |
34 | #include <math.h>\r | |
35 | #include <stdlib.h>\r | |
36 | #include <machine/ieee.h>\r | |
37 | #include "gdtoaimp.h"\r | |
38 | \r | |
39 | #if defined(_MSC_VER)\r | |
40 | /* Disable warnings about conversions to narrower data types,\r | |
41 | primarily for the fpclassify() macro.\r | |
42 | */\r | |
43 | #pragma warning ( disable : 4244 )\r | |
44 | // Squelch bogus warnings about uninitialized variable use.\r | |
45 | #pragma warning ( disable : 4700 )\r | |
46 | #endif\r | |
47 | \r | |
48 | /*\r | |
49 | * ldtoa() is a wrapper for gdtoa() that makes it smell like dtoa(),\r | |
50 | * except that the floating point argument is passed by reference.\r | |
51 | * When dtoa() is passed a NaN or infinity, it sets expt to 9999.\r | |
52 | * However, a long double could have a valid exponent of 9999, so we\r | |
53 | * use INT_MAX in ldtoa() instead.\r | |
54 | */\r | |
55 | char *\r | |
56 | ldtoa(long double *ld, int mode, int ndigits, int *decpt, int *sign, char **rve)\r | |
57 | {\r | |
58 | #ifdef EXT_EXPBITS\r | |
59 | static FPI fpi = {\r | |
60 | LDBL_MANT_DIG, /* nbits */\r | |
61 | LDBL_MIN_EXP - LDBL_MANT_DIG, /* emin */\r | |
62 | LDBL_MAX_EXP - LDBL_MANT_DIG, /* emax */\r | |
63 | FPI_Round_near, /* rounding */\r | |
64 | #ifdef Sudden_Underflow /* unused, but correct anyway */\r | |
65 | 1\r | |
66 | #else\r | |
67 | 0\r | |
68 | #endif\r | |
69 | };\r | |
70 | int be, kind;\r | |
71 | char *ret;\r | |
72 | union ieee_ext_u u;\r | |
73 | uint32_t bits[(LDBL_MANT_DIG + 31) / 32];\r | |
74 | \r | |
75 | u.extu_ld = *ld;\r | |
76 | *sign = (int)(u.extu_ext.ext_sign);\r | |
77 | be = (int)(u.extu_ext.ext_exp - (LDBL_MAX_EXP - 1) - (LDBL_MANT_DIG - 1));\r | |
78 | EXT_TO_ARRAY32(u, bits);\r | |
79 | \r | |
80 | switch (fpclassify(u.extu_ld)) {\r | |
81 | case FP_NORMAL:\r | |
82 | kind = STRTOG_Normal;\r | |
83 | #ifdef LDBL_IMPLICIT_NBIT\r | |
84 | bits[LDBL_MANT_DIG / 32] |= 1 << ((LDBL_MANT_DIG - 1) % 32);\r | |
85 | #endif /* LDBL_IMPLICIT_NBIT */\r | |
86 | break;\r | |
87 | case FP_ZERO:\r | |
88 | kind = STRTOG_Zero;\r | |
89 | break;\r | |
90 | case FP_SUBNORMAL:\r | |
91 | kind = STRTOG_Denormal;\r | |
92 | #ifdef LDBL_IMPLICIT_NBIT\r | |
93 | be++;\r | |
94 | #endif\r | |
95 | break;\r | |
96 | case FP_INFINITE:\r | |
97 | kind = STRTOG_Infinite;\r | |
98 | break;\r | |
99 | case FP_NAN:\r | |
100 | kind = STRTOG_NaN;\r | |
101 | break;\r | |
102 | default:\r | |
103 | abort();\r | |
104 | }\r | |
105 | \r | |
106 | ret = gdtoa(&fpi, be, (ULong *)bits, &kind, mode, ndigits, decpt, rve);\r | |
107 | if (*decpt == -32768)\r | |
108 | *decpt = INT_MAX;\r | |
109 | return ret;\r | |
110 | #else\r | |
111 | return dtoa((double)*ld, mode, ndigits, decpt, sign, rve);\r | |
112 | #endif\r | |
113 | }\r |