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