da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin/***********************************************************************
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin* *
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin* This software is part of the ast package *
3e14f97f673e8a630f076077de35afdd43dc1587Roger A. Faulkner* Copyright (c) 1985-2010 AT&T Intellectual Property *
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin* and is licensed under the *
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin* Common Public License, Version 1.0 *
7c2fbfb345896881c631598ee3852ce9ce33fb07April Chin* by AT&T Intellectual Property *
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin* *
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin* A copy of the License is available at *
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin* http://www.opensource.org/licenses/cpl1.0.txt *
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin* (with md5 checksum 059e8cd6165cb4c31e351f2b69388fd9) *
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin* *
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin* Information and Software Systems Research *
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin* AT&T Research *
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin* Florham Park NJ *
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin* *
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin* Glenn Fowler <gsf@research.att.com> *
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin* David Korn <dgk@research.att.com> *
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin* Phong Vo <kpv@research.att.com> *
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin* *
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin***********************************************************************/
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin#include "sfhdr.h"
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin/* Convert a Sfdouble_t value represented in an ASCII format into
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin** the internal Sfdouble_t representation.
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin**
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin** Written by Kiem-Phong Vo.
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin*/
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin#define BATCH (2*sizeof(int)) /* accumulate this many digits at a time */
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin#define IPART 0 /* doing integer part */
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin#define FPART 1 /* doing fractional part */
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin#define EPART 2 /* doing exponent part */
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin#if __STD_C
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chinstatic Sfdouble_t sfpow10(reg int n)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin#else
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chinstatic Sfdouble_t sfpow10(n)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chinreg int n;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin#endif
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin{
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin Sfdouble_t dval;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin switch(n)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin { case -3: return .001;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin case -2: return .01;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin case -1: return .1;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin case 0: return 1.;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin case 1: return 10.;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin case 2: return 100.;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin case 3: return 1000.;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin }
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin if(n < 0)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin { dval = .0001;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin for(n += 4; n < 0; n += 1)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin dval /= 10.;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin }
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin else
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin { dval = 10000.;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin for(n -= 4; n > 0; n -= 1)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin dval *= 10.;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin }
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin return dval;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin}
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin#if __STD_C
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chinSfdouble_t _sfstrtod(reg const char* s, char** retp)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin#else
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chinSfdouble_t _sfstrtod(s,retp)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chinreg char* s; /* string to convert */
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chinchar** retp; /* to return the remainder of string */
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin#endif
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin{
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin reg int n, c, m;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin reg int mode, fexp, sign, expsign;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin Sfdouble_t dval;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin#if _lib_locale
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin int decpoint = 0;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin int thousand = 0;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin SFSETLOCALE(&decpoint,&thousand);
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin#else
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin#define decpoint '.'
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin#endif
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin /* skip initial blanks */
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin while(isspace(*s))
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin ++s;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin /* get the sign */
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin if((sign = (*s == '-')) || *s == '+')
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin s += 1;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin mode = IPART;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin fexp = expsign = 0;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin dval = 0.;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin while(*s)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin { /* accumulate a handful of the digits */
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin for(m = BATCH, n = 0; m > 0; --m, ++s)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin { /* get and process a char */
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin c = *s;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin if(isdigit(c))
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin n = 10*n + (c - '0');
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin else break;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin }
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin /* number of digits accumulated */
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin m = BATCH-m;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin if(mode == IPART)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin { /* doing the integer part */
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin if(dval == 0.)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin dval = (Sfdouble_t)n;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin else dval = dval*sfpow10(m) + (Sfdouble_t)n;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin }
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin else if(mode == FPART)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin { /* doing the fractional part */
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin fexp -= m;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin if(n > 0)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin dval += n*sfpow10(fexp);
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin }
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin else if(n)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin { /* doing the exponent part */
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin if(expsign)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin n = -n;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin dval *= sfpow10(n);
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin }
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin if(!c)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin break;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin if(m < BATCH)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin { /* detected a non-digit */
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin if(c == decpoint)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin { /* start the fractional part or no match */
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin if(mode != IPART)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin break;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin mode = FPART;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin s += 1;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin }
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin else if(c == 'e' || c == 'E')
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin { if(mode == EPART)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin break;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin mode = EPART;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin c = *++s;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin if((expsign = (c == '-')) || c == '+')
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin s += 1;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin }
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin else break;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin }
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin }
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin if(retp)
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin *retp = (char*)s;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin return sign ? -dval : dval;
da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968chin}