numparser2.h

/*
 * libtu/numparser2.h
 *
 * Copyright (c) Tuomo Valkonen 1999-2000.
 * 
 * This file is distributed under the terms of the "Artistic License".
 * See the included file LICENSE for details.
 */

#define MAX_MANTISSA 10  /* should be enough for our needs */
#define ULONG_SIZE (sizeof(ulong)*8)

enum{
	NPNUM_INT,
	NPNUM_FLOAT
};

#ifdef NP_SIMPLE_IMPL

typedef struct _NPNum{
	int type;
	int base;
 	bool negative;
	long double fval;
	ulong ival;
} NPNum;

#define NUM_INIT {0, 0, 0, 0.0, 0}

static int npnum_mulbase_add(NPNum *num, long base, long v)
{
	long double iold=num->ival;
	
	num->fval=num->fval*base+(long double)v;
	
	num->ival*=base;
	
	if(num->ival<iold)
		num->type=NPNUM_FLOAT;
	
	num->ival+=v;
	
	return 0;
}

#else /* NP_SIMPLE_IMPL */

typedef struct _NPNum{
	unsigned char nmantissa;
	int type;
	int base;
 	bool negative;
	ulong mantissa[MAX_MANTISSA];
	long exponent;
} NPNum;

#define NUM_INIT {0, 0, 0, 0, {0,}, 0}

#define ADD_EXP(NUM, X) (NUM)->exponent+=(X);

#if defined(__GNUG__) && defined(i386) && !defined(NP_NO_I386_ASM)
 #define NP_I386_ASM
#endif

static int npnum_mulbase_add(NPNum *num, long base, long v)
{
	long i, j;
	ulong overflow;
#ifndef NP_I386_ASM
	ulong val;
#endif
	
	for(i=num->nmantissa;i>=0;i--){
#ifdef NP_I386_ASM
		__asm__("mul %4\n"
				: "=a" (num->mantissa[i]), "=d" (overflow)
				: "0" (num->mantissa[i]), "1" (0), "q" (base)
				: "eax", "edx");
#else
		overflow=0;
		val=num->mantissa[i];
		
		if(val<ULONG_MAX/base){
			val*=base;
		}else if(val){
			ulong tmp=val;
			ulong old=val;
			for(j=0; j<base; j++){
				val+=tmp;
				if(val<=old)
					overflow++;
				old=val;
			}
		}
		num->mantissa[i]=val;
#endif
		if(overflow){
			if(i==num->nmantissa){
				if(num->nmantissa==MAX_MANTISSA)
					return E_TOKZ_TOOBIG;
				num->nmantissa++;
			}
			num->mantissa[i+1]+=overflow;
		}
	}
	num->mantissa[0]+=v;
	
	return 0;
}

#undef NP_I386_ASM

#endif /* NP_SIMPLE_IMPL */


/* */


static bool is_end(int c)
{
	/* oops... EOF was missing */
	return (c==EOF || (c!='.' && ispunct(c)) || isspace(c) || iscntrl(c));
}


/* */


static int parse_exponent(NPNum *num, Tokenizer *tokz, int c)
{
	long exp=0;
	bool neg=FALSE;
	int err=0;
	
	c=GETCH();
	
	if(c=='-' || c=='+'){
		if(c=='-')
			neg=TRUE;
		c=GETCH();
	}
	
	for(; 1; c=GETCH()){
		if(isdigit(c)){
			exp*=10;
			exp+=c-'0';
		}else if(is_end(c)){
			UNGETCH(c);
			break;
		}else{
			err=E_TOKZ_NUMFMT;
		}
	}

	if(neg)
		exp*=-1;

#ifndef NP_SIMPLE_IMPL
	ADD_EXP(num, exp);
#else
	num->fval*=pow(num->base, exp);
#endif	
	return err;
}


static int parse_number(NPNum *num, Tokenizer *tokz, int c)
{
	int base=10;
	int dm=1;
	int err=0;
	int tmp;
#ifdef NP_SIMPLE_IMPL
	long double divisor=base;
#endif	
	
	if(c=='-' || c=='+'){
		if(c=='-')
			num->negative=TRUE;
		c=GETCH();
		if(!isdigit(c))
			err=E_TOKZ_NUMFMT;
	}
	
	if(c=='0'){
		dm=0;
		c=GETCH();
		if(c=='x' || c=='X'){
			base=16;
			c=GETCH();
		}else if(c=='b' || c=='B'){
			base=2;
			c=GETCH();
		}else if('0'<=c && c<='7'){
			base=8;
		}else{
			dm=2;
		}
	}
	
	num->base=base;
	
	for(; 1; c=GETCH()){
		if((c=='e' || c=='E') && dm!=0){
			if(dm<2){
				err=E_TOKZ_NUMFMT;
				continue;
			}
			tmp=parse_exponent(num, tokz, c);
			if(err==0)
				err=tmp;
			break;
		}
		
		if(isxdigit(c)){
			if('0'<=c && c<='9')
				c-='0';
			else if(isupper(c))
				c-='A'-10;
			else
				c-='a'-10;
			
			if(c>=base)
				err=E_TOKZ_NUMFMT;

#ifdef NP_SIMPLE_IMPL
			if(dm==3){
				num->fval+=(long double)c/divisor;
				divisor*=base;
			}else
#endif
			{
				tmp=npnum_mulbase_add(num, base, c);
				if(err==0)
					err=tmp;
			}
			
			if(dm==1)
				dm=2;
#ifndef NP_SIMPLE_IMPL			
			else if(dm==3)
				ADD_EXP(num, -1);
#endif			
			continue;
		}
		
		if(c=='.'){
			if(dm!=2){
				err=E_TOKZ_NUMFMT;
			}
			dm=3;
#ifdef NP_SIMPLE_IMPL
			num->type=NPNUM_FLOAT;
			divisor=base;
#endif
			continue;
		}
		
		if(is_end(c)){
			UNGETCH(c);
			break;
		}
		
		err=E_TOKZ_NUMFMT;
	}
	
#ifndef NP_SIMPLE_IMPL			
	num->type=(num->exponent==0 ? NPNUM_INT : NPNUM_FLOAT);
#endif

	return err;
}