Home | History | Annotate | Line # | Download | only in gdtoa
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