]>
Commit | Line | Data |
---|---|---|
2aa62f2b | 1 | /* $NetBSD: hexnan.c,v 1.3 2006/03/11 18:38:14 kleink Exp $ */\r |
2 | \r | |
3 | /****************************************************************\r | |
4 | \r | |
5 | The author of this software is David M. Gay.\r | |
6 | \r | |
7 | Copyright (C) 2000 by Lucent Technologies\r | |
8 | All Rights Reserved\r | |
9 | \r | |
10 | Permission to use, copy, modify, and distribute this software and\r | |
11 | its documentation for any purpose and without fee is hereby\r | |
12 | granted, provided that the above copyright notice appear in all\r | |
13 | copies and that both that the copyright notice and this\r | |
14 | permission notice and warranty disclaimer appear in supporting\r | |
15 | documentation, and that the name of Lucent or any of its entities\r | |
16 | not be used in advertising or publicity pertaining to\r | |
17 | distribution of the software without specific, written prior\r | |
18 | permission.\r | |
19 | \r | |
20 | LUCENT DISCLAIMS ALL WARRANTIES WITH REGARD TO THIS SOFTWARE,\r | |
21 | INCLUDING ALL IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS.\r | |
22 | IN NO EVENT SHALL LUCENT OR ANY OF ITS ENTITIES BE LIABLE FOR ANY\r | |
23 | SPECIAL, INDIRECT OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES\r | |
24 | WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER\r | |
25 | IN AN ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION,\r | |
26 | ARISING OUT OF OR IN CONNECTION WITH THE USE OR PERFORMANCE OF\r | |
27 | THIS SOFTWARE.\r | |
28 | \r | |
29 | ****************************************************************/\r | |
30 | \r | |
31 | /* Please send bug reports to David M. Gay (dmg at acm dot org,\r | |
32 | * with " at " changed at "@" and " dot " changed to "."). */\r | |
33 | #include <LibConfig.h>\r | |
34 | \r | |
35 | #include "gdtoaimp.h"\r | |
36 | \r | |
37 | static void\r | |
38 | #ifdef KR_headers\r | |
39 | L_shift(x, x1, i) ULong *x; ULong *x1; int i;\r | |
40 | #else\r | |
41 | L_shift(ULong *x, ULong *x1, int i)\r | |
42 | #endif\r | |
43 | {\r | |
44 | int j;\r | |
45 | \r | |
46 | i = 8 - i;\r | |
47 | i <<= 2;\r | |
48 | j = ULbits - i;\r | |
49 | do {\r | |
50 | *x |= x[1] << j;\r | |
51 | x[1] >>= i;\r | |
52 | } while(++x < x1);\r | |
53 | }\r | |
54 | \r | |
55 | int\r | |
56 | #ifdef KR_headers\r | |
57 | hexnan(sp, fpi, x0)\r | |
58 | CONST char **sp; CONST FPI *fpi; ULong *x0;\r | |
59 | #else\r | |
60 | hexnan( CONST char **sp, CONST FPI *fpi, ULong *x0)\r | |
61 | #endif\r | |
62 | {\r | |
63 | ULong c, h, *x, *x1, *xe;\r | |
64 | CONST char *s;\r | |
65 | int havedig, hd0, i, nbits;\r | |
66 | \r | |
67 | if (!hexdig['0'])\r | |
68 | hexdig_init_D2A();\r | |
69 | nbits = fpi->nbits;\r | |
70 | x = x0 + ((unsigned int)nbits >> kshift);\r | |
71 | if (nbits & kmask)\r | |
72 | x++;\r | |
73 | *--x = 0;\r | |
74 | x1 = xe = x;\r | |
75 | havedig = hd0 = i = 0;\r | |
76 | s = *sp;\r | |
77 | while((c = *(CONST unsigned char*)++s) != 0) {\r | |
78 | if ((h = hexdig[c]) == 0) {\r | |
79 | if (c <= ' ') {\r | |
80 | if (hd0 < havedig) {\r | |
81 | if (x < x1 && i < 8)\r | |
82 | L_shift(x, x1, i);\r | |
83 | if (x <= x0) {\r | |
84 | i = 8;\r | |
85 | continue;\r | |
86 | }\r | |
87 | hd0 = havedig;\r | |
88 | *--x = 0;\r | |
89 | x1 = x;\r | |
90 | i = 0;\r | |
91 | }\r | |
92 | continue;\r | |
93 | }\r | |
94 | if (/*(*/ c == ')' && havedig) {\r | |
95 | *sp = s + 1;\r | |
96 | break;\r | |
97 | }\r | |
98 | return STRTOG_NaN;\r | |
99 | }\r | |
100 | havedig++;\r | |
101 | if (++i > 8) {\r | |
102 | if (x <= x0)\r | |
103 | continue;\r | |
104 | i = 1;\r | |
105 | *--x = 0;\r | |
106 | }\r | |
107 | *x = (*x << 4) | (h & 0xf);\r | |
108 | }\r | |
109 | if (!havedig)\r | |
110 | return STRTOG_NaN;\r | |
111 | if (x < x1 && i < 8)\r | |
112 | L_shift(x, x1, i);\r | |
113 | if (x > x0) {\r | |
114 | x1 = x0;\r | |
115 | do *x1++ = *x++;\r | |
116 | while(x <= xe);\r | |
117 | do *x1++ = 0;\r | |
118 | while(x1 <= xe);\r | |
119 | }\r | |
120 | else {\r | |
121 | /* truncate high-order word if necessary */\r | |
122 | if ( (i = nbits & (ULbits-1)) !=0)\r | |
123 | *xe &= ((ULong)0xffffffff) >> (ULbits - i);\r | |
124 | }\r | |
125 | for(x1 = xe;; --x1) {\r | |
126 | if (*x1 != 0)\r | |
127 | break;\r | |
128 | if (x1 == x0) {\r | |
129 | *x1 = 1;\r | |
130 | break;\r | |
131 | }\r | |
132 | }\r | |
133 | return STRTOG_NaNbits;\r | |
134 | }\r |