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