strtorQ.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 void
53 1.1 kleink #ifdef KR_headers
54 1.1 kleink ULtoQ(L, bits, exp, k) ULong *L; ULong *bits; Long exp; int k;
55 1.1 kleink #else
56 1.1 kleink ULtoQ(ULong *L, ULong *bits, Long exp, int k)
57 1.1 kleink #endif
58 1.1 kleink {
59 1.1 kleink switch(k & STRTOG_Retmask) {
60 1.1 kleink case STRTOG_NoNumber:
61 1.1 kleink case STRTOG_Zero:
62 1.1 kleink L[0] = L[1] = L[2] = L[3] = 0;
63 1.1 kleink break;
64 1.1 kleink
65 1.1 kleink case STRTOG_Normal:
66 1.1 kleink case STRTOG_NaNbits:
67 1.1 kleink L[_3] = bits[0];
68 1.1 kleink L[_2] = bits[1];
69 1.1 kleink L[_1] = bits[2];
70 1.1 kleink L[_0] = (bits[3] & ~0x10000) | ((exp + 0x3fff + 112) << 16);
71 1.1 kleink break;
72 1.1 kleink
73 1.1 kleink case STRTOG_Denormal:
74 1.1 kleink L[_3] = bits[0];
75 1.1 kleink L[_2] = bits[1];
76 1.1 kleink L[_1] = bits[2];
77 1.1 kleink L[_0] = bits[3];
78 1.1 kleink break;
79 1.1 kleink
80 1.1 kleink case STRTOG_Infinite:
81 1.1 kleink L[_0] = 0x7fff0000;
82 1.1 kleink L[_1] = L[_2] = L[_3] = 0;
83 1.1 kleink break;
84 1.1 kleink
85 1.1 kleink case STRTOG_NaN:
86 1.1 kleink L[0] = ld_QNAN0;
87 1.1 kleink L[1] = ld_QNAN1;
88 1.1 kleink L[2] = ld_QNAN2;
89 1.1 kleink L[3] = ld_QNAN3;
90 1.1 kleink }
91 1.1 kleink if (k & STRTOG_Neg)
92 1.1 kleink L[_0] |= 0x80000000L;
93 1.1 kleink }
94 1.1 kleink
95 1.1 kleink int
96 1.1 kleink #ifdef KR_headers
97 1.1 kleink strtorQ(s, sp, rounding, L) CONST char *s; char **sp; int rounding; void *L;
98 1.1 kleink #else
99 1.1 kleink strtorQ(CONST char *s, char **sp, int rounding, void *L)
100 1.1 kleink #endif
101 1.1 kleink {
102 1.1 kleink static FPI fpi0 = { 113, 1-16383-113+1, 32766-16383-113+1, 1, SI };
103 1.1 kleink FPI *fpi, fpi1;
104 1.1 kleink ULong bits[4];
105 1.1 kleink Long exp;
106 1.1 kleink int k;
107 1.1 kleink
108 1.1 kleink fpi = &fpi0;
109 1.1 kleink if (rounding != FPI_Round_near) {
110 1.1 kleink fpi1 = fpi0;
111 1.1 kleink fpi1.rounding = rounding;
112 1.1 kleink fpi = &fpi1;
113 1.1 kleink }
114 1.1 kleink k = strtodg(s, sp, fpi, &exp, bits);
115 1.1 kleink ULtoQ((ULong*)L, bits, exp, k);
116 1.1 kleink return k;
117 1.1 kleink }
118