1112158Sdas/**************************************************************** 2112158Sdas 3112158SdasThe author of this software is David M. Gay. 4112158Sdas 5112158SdasCopyright (C) 1998, 2000 by Lucent Technologies 6112158SdasAll Rights Reserved 7112158Sdas 8112158SdasPermission to use, copy, modify, and distribute this software and 9112158Sdasits documentation for any purpose and without fee is hereby 10112158Sdasgranted, provided that the above copyright notice appear in all 11112158Sdascopies and that both that the copyright notice and this 12112158Sdaspermission notice and warranty disclaimer appear in supporting 13112158Sdasdocumentation, and that the name of Lucent or any of its entities 14112158Sdasnot be used in advertising or publicity pertaining to 15112158Sdasdistribution of the software without specific, written prior 16112158Sdaspermission. 17112158Sdas 18112158SdasLUCENT DISCLAIMS ALL WARRANTIES WITH REGARD TO THIS SOFTWARE, 19112158SdasINCLUDING ALL IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS. 20112158SdasIN NO EVENT SHALL LUCENT OR ANY OF ITS ENTITIES BE LIABLE FOR ANY 21112158SdasSPECIAL, INDIRECT OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES 22112158SdasWHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER 23112158SdasIN AN ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, 24112158SdasARISING OUT OF OR IN CONNECTION WITH THE USE OR PERFORMANCE OF 25112158SdasTHIS SOFTWARE. 26112158Sdas 27112158Sdas****************************************************************/ 28112158Sdas 29165743Sdas/* Please send bug reports to David M. Gay (dmg at acm dot org, 30165743Sdas * with " at " changed at "@" and " dot " changed to "."). */ 31112158Sdas 32174679Sdas/* $FreeBSD$ */ 33174679Sdas 34112158Sdas#include "gdtoaimp.h" 35112158Sdas 36112158Sdas#undef _0 37112158Sdas#undef _1 38112158Sdas 39112158Sdas/* one or the other of IEEE_MC68k or IEEE_8087 should be #defined */ 40112158Sdas 41112158Sdas#ifdef IEEE_MC68k 42112158Sdas#define _0 0 43112158Sdas#define _1 1 44112158Sdas#define _2 2 45112158Sdas#define _3 3 46112158Sdas#endif 47112158Sdas#ifdef IEEE_8087 48112158Sdas#define _0 3 49112158Sdas#define _1 2 50112158Sdas#define _2 1 51112158Sdas#define _3 0 52112158Sdas#endif 53112158Sdas 54112158Sdas void 55112158Sdas#ifdef KR_headers 56112158SdasULtoQ(L, bits, exp, k) ULong *L; ULong *bits; Long exp; int k; 57112158Sdas#else 58112158SdasULtoQ(ULong *L, ULong *bits, Long exp, int k) 59112158Sdas#endif 60112158Sdas{ 61112158Sdas switch(k & STRTOG_Retmask) { 62112158Sdas case STRTOG_NoNumber: 63112158Sdas case STRTOG_Zero: 64112158Sdas L[0] = L[1] = L[2] = L[3] = 0; 65112158Sdas break; 66112158Sdas 67112158Sdas case STRTOG_Normal: 68174690Sdas L[_3] = bits[0]; 69174690Sdas L[_2] = bits[1]; 70174690Sdas L[_1] = bits[2]; 71174690Sdas L[_0] = (bits[3] & ~0x10000) | ((exp + 0x3fff + 112) << 16); 72174690Sdas break; 73174690Sdas 74112158Sdas case STRTOG_NaNbits: 75112158Sdas L[_3] = bits[0]; 76112158Sdas L[_2] = bits[1]; 77112158Sdas L[_1] = bits[2]; 78174679Sdas L[_0] = (bits[3] & ~0x10000) 79174679Sdas | (((exp + 0x3fff + 112) << 16) | (1 << 15)); 80112158Sdas break; 81112158Sdas 82112158Sdas case STRTOG_Denormal: 83112158Sdas L[_3] = bits[0]; 84112158Sdas L[_2] = bits[1]; 85112158Sdas L[_1] = bits[2]; 86112158Sdas L[_0] = bits[3]; 87112158Sdas break; 88112158Sdas 89112158Sdas case STRTOG_Infinite: 90112158Sdas L[_0] = 0x7fff0000; 91112158Sdas L[_1] = L[_2] = L[_3] = 0; 92112158Sdas break; 93112158Sdas 94112158Sdas case STRTOG_NaN: 95165743Sdas L[0] = ld_QNAN0; 96165743Sdas L[1] = ld_QNAN1; 97165743Sdas L[2] = ld_QNAN2; 98165743Sdas L[3] = ld_QNAN3; 99112158Sdas } 100112158Sdas if (k & STRTOG_Neg) 101112158Sdas L[_0] |= 0x80000000L; 102112158Sdas } 103112158Sdas 104112158Sdas int 105112158Sdas#ifdef KR_headers 106227753StheravenstrtorQ_l(s, sp, rounding, L, locale) CONST char *s; char **sp; int rounding; 107227753Stheravenvoid *L; locale_t locale; 108112158Sdas#else 109227753StheravenstrtorQ_l(CONST char *s, char **sp, int rounding, void *L, locale_t locale) 110112158Sdas#endif 111112158Sdas{ 112112158Sdas static FPI fpi0 = { 113, 1-16383-113+1, 32766-16383-113+1, 1, SI }; 113112158Sdas FPI *fpi, fpi1; 114112158Sdas ULong bits[4]; 115112158Sdas Long exp; 116112158Sdas int k; 117112158Sdas 118112158Sdas fpi = &fpi0; 119112158Sdas if (rounding != FPI_Round_near) { 120112158Sdas fpi1 = fpi0; 121112158Sdas fpi1.rounding = rounding; 122112158Sdas fpi = &fpi1; 123112158Sdas } 124227753Stheraven k = strtodg_l(s, sp, fpi, &exp, bits, locale); 125112158Sdas ULtoQ((ULong*)L, bits, exp, k); 126112158Sdas return k; 127112158Sdas } 128