strtopQ.c revision 1.1.1.2 1 1.1 kleink /****************************************************************
2 1.1 kleink
3 1.1 kleink The author of this software is David M. Gay.
4 1.1 kleink
5 1.1 kleink Copyright (C) 1998, 2000 by Lucent Technologies
6 1.1 kleink All Rights Reserved
7 1.1 kleink
8 1.1 kleink Permission to use, copy, modify, and distribute this software and
9 1.1 kleink its documentation for any purpose and without fee is hereby
10 1.1 kleink granted, provided that the above copyright notice appear in all
11 1.1 kleink copies and that both that the copyright notice and this
12 1.1 kleink permission notice and warranty disclaimer appear in supporting
13 1.1 kleink documentation, and that the name of Lucent or any of its entities
14 1.1 kleink not be used in advertising or publicity pertaining to
15 1.1 kleink distribution of the software without specific, written prior
16 1.1 kleink permission.
17 1.1 kleink
18 1.1 kleink LUCENT DISCLAIMS ALL WARRANTIES WITH REGARD TO THIS SOFTWARE,
19 1.1 kleink INCLUDING ALL IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS.
20 1.1 kleink IN NO EVENT SHALL LUCENT OR ANY OF ITS ENTITIES BE LIABLE FOR ANY
21 1.1 kleink SPECIAL, INDIRECT OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES
22 1.1 kleink WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER
23 1.1 kleink IN AN ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION,
24 1.1 kleink ARISING OUT OF OR IN CONNECTION WITH THE USE OR PERFORMANCE OF
25 1.1 kleink THIS SOFTWARE.
26 1.1 kleink
27 1.1 kleink ****************************************************************/
28 1.1 kleink
29 1.1 kleink /* Please send bug reports to David M. Gay (dmg at acm dot org,
30 1.1 kleink * with " at " changed at "@" and " dot " changed to "."). */
31 1.1 kleink
32 1.1 kleink #include "gdtoaimp.h"
33 1.1 kleink
34 1.1 kleink #undef _0
35 1.1 kleink #undef _1
36 1.1 kleink
37 1.1 kleink /* one or the other of IEEE_MC68k or IEEE_8087 should be #defined */
38 1.1 kleink
39 1.1 kleink #ifdef IEEE_MC68k
40 1.1 kleink #define _0 0
41 1.1 kleink #define _1 1
42 1.1 kleink #define _2 2
43 1.1 kleink #define _3 3
44 1.1 kleink #endif
45 1.1 kleink #ifdef IEEE_8087
46 1.1 kleink #define _0 3
47 1.1 kleink #define _1 2
48 1.1 kleink #define _2 1
49 1.1 kleink #define _3 0
50 1.1 kleink #endif
51 1.1 kleink
52 1.1 kleink int
53 1.1 kleink #ifdef KR_headers
54 1.1 kleink strtopQ(s, sp, V) CONST char *s; char **sp; void *V;
55 1.1 kleink #else
56 1.1 kleink strtopQ(CONST char *s, char **sp, void *V)
57 1.1 kleink #endif
58 1.1 kleink {
59 1.1.1.2 christos static FPI fpi0 = { 113, 1-16383-113+1, 32766 - 16383 - 113 + 1, 1, SI };
60 1.1 kleink ULong bits[4];
61 1.1 kleink Long exp;
62 1.1 kleink int k;
63 1.1 kleink ULong *L = (ULong*)V;
64 1.1.1.2 christos #ifdef Honor_FLT_ROUNDS
65 1.1.1.2 christos #include "gdtoa_fltrnds.h"
66 1.1.1.2 christos #else
67 1.1.1.2 christos #define fpi &fpi0
68 1.1.1.2 christos #endif
69 1.1 kleink
70 1.1.1.2 christos k = strtodg(s, sp, fpi, &exp, bits);
71 1.1 kleink switch(k & STRTOG_Retmask) {
72 1.1 kleink case STRTOG_NoNumber:
73 1.1 kleink case STRTOG_Zero:
74 1.1 kleink L[0] = L[1] = L[2] = L[3] = 0;
75 1.1 kleink break;
76 1.1 kleink
77 1.1 kleink case STRTOG_Normal:
78 1.1 kleink case STRTOG_NaNbits:
79 1.1 kleink L[_3] = bits[0];
80 1.1 kleink L[_2] = bits[1];
81 1.1 kleink L[_1] = bits[2];
82 1.1 kleink L[_0] = (bits[3] & ~0x10000) | ((exp + 0x3fff + 112) << 16);
83 1.1 kleink break;
84 1.1 kleink
85 1.1 kleink case STRTOG_Denormal:
86 1.1 kleink L[_3] = bits[0];
87 1.1 kleink L[_2] = bits[1];
88 1.1 kleink L[_1] = bits[2];
89 1.1 kleink L[_0] = bits[3];
90 1.1 kleink break;
91 1.1 kleink
92 1.1 kleink case STRTOG_Infinite:
93 1.1 kleink L[_0] = 0x7fff0000;
94 1.1 kleink L[_1] = L[_2] = L[_3] = 0;
95 1.1 kleink break;
96 1.1 kleink
97 1.1 kleink case STRTOG_NaN:
98 1.1 kleink L[0] = ld_QNAN0;
99 1.1 kleink L[1] = ld_QNAN1;
100 1.1 kleink L[2] = ld_QNAN2;
101 1.1 kleink L[3] = ld_QNAN3;
102 1.1 kleink }
103 1.1 kleink if (k & STRTOG_Neg)
104 1.1 kleink L[_0] |= 0x80000000L;
105 1.1 kleink return k;
106 1.1 kleink }
107