gonzui

    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 David M. Gay (dmg at acm dot org,
   30:  * with " at " changed at "@" and " dot " changed to ".").      */
   31: 
   32: #include "gdtoaimp.h"
   33: 
   34:  void
   35: #ifdef KR_headers
   36: ULtod(L, bits, exp, k) ULong *L; ULong *bits; Long exp; int k;
   37: #else
   38: ULtod(ULong *L, ULong *bits, Long exp, int k)
   39: #endif
   40: {
   41:         switch(k & STRTOG_Retmask) {
   42:           case STRTOG_NoNumber:
   43:           case STRTOG_Zero:
   44:                 L[0] = L[1] = 0;
   45:                 break;
   46: 
   47:           case STRTOG_Denormal:
   48:                 L[_1] = bits[0];
   49:                 L[_0] = bits[1];
   50:                 break;
   51: 
   52:           case STRTOG_Normal:
   53:           case STRTOG_NaNbits:
   54:                 L[_1] = bits[0];
   55:                 L[_0] = (bits[1] & ~0x100000) | ((exp + 0x3ff + 52) << 20);
   56:                 break;
   57: 
   58:           case STRTOG_NoMemory:
   59:                 /*errno = ERANGE*/;
   60:                 /* FALLTHROUGH */
   61:           case STRTOG_Infinite:
   62:                 L[_0] = 0x7ff00000;
   63:                 L[_1] = 0;
   64:                 break;
   65: 
   66:           case STRTOG_NaN:
   67:                 L[0] = d_QNAN0;
   68:                 L[1] = d_QNAN1;
   69:           }
   70:         if (k & STRTOG_Neg)
   71:                 L[_0] |= 0x80000000L;
   72:         }
   73: 
   74:  int
   75: #ifdef KR_headers
   76: strtord(s, sp, rounding, d) CONST char *s; char **sp; int rounding; double *d;
   77: #else
   78: strtord(CONST char *s, char **sp, int rounding, double *d)
   79: #endif
   80: {
   81:         static FPI fpi0 = { 53, 1-1023-53+1, 2046-1023-53+1, 1, SI };
   82:         FPI *fpi, fpi1;
   83:         ULong bits[2];
   84:         Long exp;
   85:         int k;
   86: 
   87:         fpi = &fpi0;
   88:         if (rounding != FPI_Round_near) {
   89:                 fpi1 = fpi0;
   90:                 fpi1.rounding = rounding;
   91:                 fpi = &fpi1;
   92:                 }
   93:         k = strtodg(s, sp, fpi, &exp, bits);
   94:         ULtod((ULong*)d, bits, exp, k);
   95:         return k;
   96:         }