Home | History | Annotate | Line # | Download | only in gdtoa
      1  1.4  christos /* $NetBSD: strtorx.c,v 1.4 2011/03/20 23:15:35 christos 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 #define _3 3
     46  1.1    kleink #define _4 4
     47  1.1    kleink #endif
     48  1.2    kleink #ifdef IEEE_LITTLE_ENDIAN
     49  1.1    kleink #define _0 4
     50  1.1    kleink #define _1 3
     51  1.1    kleink #define _2 2
     52  1.1    kleink #define _3 1
     53  1.1    kleink #define _4 0
     54  1.1    kleink #endif
     55  1.1    kleink 
     56  1.1    kleink  void
     57  1.1    kleink #ifdef KR_headers
     58  1.4  christos ULtox(L, bits, expt, k) UShort *L; ULong *bits; Long expt; int k;
     59  1.1    kleink #else
     60  1.4  christos ULtox(UShort *L, ULong *bits, Long expt, int k)
     61  1.1    kleink #endif
     62  1.1    kleink {
     63  1.1    kleink 	switch(k & STRTOG_Retmask) {
     64  1.1    kleink 	  case STRTOG_NoNumber:
     65  1.1    kleink 	  case STRTOG_Zero:
     66  1.1    kleink 		L[0] = L[1] = L[2] = L[3] = L[4] = 0;
     67  1.1    kleink 		break;
     68  1.1    kleink 
     69  1.1    kleink 	  case STRTOG_Denormal:
     70  1.1    kleink 		L[_0] = 0;
     71  1.1    kleink 		goto normal_bits;
     72  1.1    kleink 
     73  1.1    kleink 	  case STRTOG_Normal:
     74  1.1    kleink 	  case STRTOG_NaNbits:
     75  1.4  christos 		L[_0] = expt + 0x3fff + 63;
     76  1.1    kleink  normal_bits:
     77  1.1    kleink 		L[_4] = (UShort)bits[0];
     78  1.1    kleink 		L[_3] = (UShort)(bits[0] >> 16);
     79  1.1    kleink 		L[_2] = (UShort)bits[1];
     80  1.1    kleink 		L[_1] = (UShort)(bits[1] >> 16);
     81  1.1    kleink 		break;
     82  1.1    kleink 
     83  1.1    kleink 	  case STRTOG_Infinite:
     84  1.1    kleink 		L[_0] = 0x7fff;
     85  1.4  christos 		L[_1] = 0x8000;
     86  1.4  christos 		L[_2] = L[_3] = L[_4] = 0;
     87  1.1    kleink 		break;
     88  1.1    kleink 
     89  1.1    kleink 	  case STRTOG_NaN:
     90  1.1    kleink 		L[0] = ldus_QNAN0;
     91  1.1    kleink 		L[1] = ldus_QNAN1;
     92  1.1    kleink 		L[2] = ldus_QNAN2;
     93  1.1    kleink 		L[3] = ldus_QNAN3;
     94  1.1    kleink 		L[4] = ldus_QNAN4;
     95  1.1    kleink 	  }
     96  1.1    kleink 	if (k & STRTOG_Neg)
     97  1.1    kleink 		L[_0] |= 0x8000;
     98  1.1    kleink 	}
     99  1.1    kleink 
    100  1.1    kleink  int
    101  1.1    kleink #ifdef KR_headers
    102  1.1    kleink strtorx(s, sp, rounding, L) CONST char *s; char **sp; int rounding; void *L;
    103  1.1    kleink #else
    104  1.1    kleink strtorx(CONST char *s, char **sp, int rounding, void *L)
    105  1.1    kleink #endif
    106  1.1    kleink {
    107  1.4  christos 	static CONST FPI fpi0 = { 64, 1-16383-64+1, 32766 - 16383 - 64 + 1, 1, SI };
    108  1.4  christos 	CONST FPI *fpi;
    109  1.4  christos 	FPI fpi1;
    110  1.1    kleink 	ULong bits[2];
    111  1.4  christos 	Long expt;
    112  1.1    kleink 	int k;
    113  1.1    kleink 
    114  1.1    kleink 	fpi = &fpi0;
    115  1.1    kleink 	if (rounding != FPI_Round_near) {
    116  1.1    kleink 		fpi1 = fpi0;
    117  1.1    kleink 		fpi1.rounding = rounding;
    118  1.1    kleink 		fpi = &fpi1;
    119  1.1    kleink 		}
    120  1.4  christos 	k = strtodg(s, sp, fpi, &expt, bits);
    121  1.3  christos 	if (k == STRTOG_NoMemory)
    122  1.3  christos 		return k;
    123  1.4  christos 	    ULtox((UShort*)L, bits, expt, k);
    124  1.1    kleink 	return k;
    125  1.1    kleink 	}
    126