ack/h/con_float

/*
  (c) copyright 1988 by the Vrije Universiteit, Amsterdam, The Netherlands.
  See the copyright notice in the ACK home directory, in the file "Copyright".
*/

/* $Id$ */

/* 
   #define CODE_GENERATOR for code generator
   #define CODE_EXPANDER for code expander

   #define IEEEFLOAT for machines using IEEE floating point format
   #define PDPFLOAT for machines using the PDP-11 floating point format
   If none of these are defined, the format of the machine on which the
   code generator runs is used.
   Returns 1 if sz has an illegal value, 2 in case of overflow,
   and 0 if all went well.
   If neither IEEEFLOAT nor PDPFLOAT are defined, the return value is not
   trustworthy.

   Unfortunately, the IEEE standard does not define the byte-order.
   depends on the #defines
	FL_MSL_AT_LOW_ADDRESS	1 if most significant long is at low address
	FL_MSW_AT_LOW_ADDRESS	1 if most significant word is at low address
	FL_MSB_AT_LOW_ADDRESS	1 if most significant byte is at low address
*/
#ifdef IEEEFLOAT
#define USE_FLT
#endif
#ifdef PDPFLOAT
#define USE_FLT
#undef FL_MSL_AT_LOW_ADDRESS
#define FL_MSL_AT_LOW_ADDRESS 1
#undef FL_MSW_AT_LOW_ADDRESS
#define FL_MSW_AT_LOW_ADDRESS 1
#undef FL_MSB_AT_LOW_ADDRESS
#define FL_MSB_AT_LOW_ADDRESS 0
#endif

#define I0 ((FL_MSL_AT_LOW_ADDRESS ? 0 : 4) + (FL_MSW_AT_LOW_ADDRESS ? 0 : 2) \
		+ (FL_MSB_AT_LOW_ADDRESS ? 0 : 1))
#define I1 ((FL_MSL_AT_LOW_ADDRESS ? 0 : 4) + (FL_MSW_AT_LOW_ADDRESS ? 0 : 2) \
		+ (FL_MSB_AT_LOW_ADDRESS ? 1 : 0))
#define I2 ((FL_MSL_AT_LOW_ADDRESS ? 0 : 4) + (FL_MSW_AT_LOW_ADDRESS ? 2 : 0) \
		+ (FL_MSB_AT_LOW_ADDRESS ? 0 : 1))
#define I3 ((FL_MSL_AT_LOW_ADDRESS ? 0 : 4) + (FL_MSW_AT_LOW_ADDRESS ? 2 : 0) \
		+ (FL_MSB_AT_LOW_ADDRESS ? 1 : 0))
#define I4 ((FL_MSL_AT_LOW_ADDRESS ? 4 : 0) + (FL_MSW_AT_LOW_ADDRESS ? 0 : 2) \
		+ (FL_MSB_AT_LOW_ADDRESS ? 0 : 1))
#define I5 ((FL_MSL_AT_LOW_ADDRESS ? 4 : 0) + (FL_MSW_AT_LOW_ADDRESS ? 0 : 2) \
		+ (FL_MSB_AT_LOW_ADDRESS ? 1 : 0))
#define I6 ((FL_MSL_AT_LOW_ADDRESS ? 4 : 0) + (FL_MSW_AT_LOW_ADDRESS ? 2 : 0) \
		+ (FL_MSB_AT_LOW_ADDRESS ? 0 : 1))
#define I7 ((FL_MSL_AT_LOW_ADDRESS ? 4 : 0) + (FL_MSW_AT_LOW_ADDRESS ? 2 : 0) \
		+ (FL_MSB_AT_LOW_ADDRESS ? 1 : 0))

#ifndef USE_FLT
static int
float_cst(str, sz, buf)
	char *str, *buf;
	int sz;
{
	int i;
	char *p;
	float fl;
	double f;
	double atof();

	if (sz!= 4 && sz!= 8)	{
		return 1;
	}
	f = atof(str);
	if (sz == 4) {
		fl = f;
		p = (char *) &fl;
	}
	else {
		p = (char *) &f;
	}
	for (i = sz; i; i--) {
		*buf++ = *p++;
	}
	return 0;
}
#else /* USE_FLT */

#include <ctype.h>
#include <flt_arith.h>

int
float_cst(str, sz, buf)
	char *str, *buf;
	int sz;
{
	int overflow = 0;
	flt_arith e;

	if (sz!= 4 && sz!= 8)	{
		return 1;
	}
	flt_str2flt(str, &e);
#ifdef IEEEFLOAT
	if (sz == 4) {
#endif
#ifdef PDPFLOAT
		e.flt_exp += 129;
#else
		e.flt_exp += 127;
#endif
		if (e.flt_mantissa.flt_h_32 == 0) e.flt_exp = 0;
#ifdef IEEEFLOAT
		if (e.flt_mantissa.flt_h_32 & 0x80) {
			/* rounding */
			if ((e.flt_mantissa.flt_h_32 & 0xffffff00) == 0xffffff00) {
				e.flt_exp++;
				e.flt_mantissa.flt_h_32 = 0x80000000;
			}
			else {
				e.flt_mantissa.flt_h_32 += 0x80;
			}
		}
		if (e.flt_exp >= 255) {
			overflow = 1;
			e.flt_exp = 255;
			e.flt_mantissa.flt_h_32 = e.flt_mantissa.flt_l_32 = 0;
		}
		if (e.flt_exp <= 0) {
			flt_b64_sft(&(e.flt_mantissa), 1);
			if (e.flt_exp < 0) {
				flt_b64_sft(&(e.flt_mantissa), -e.flt_exp);
				e.flt_exp = 0;
			}
		}
#endif
#ifndef IEEEFLOAT
		if (sz == 4 && (e.flt_mantissa.flt_h_32 & 0x80)) {
			/* rounding */
			if ((e.flt_mantissa.flt_h_32 & 0xffffff00) == 0xffffff00) {
				e.flt_exp++;
				e.flt_mantissa.flt_h_32 = 0x80000000;
			}
			else {
				e.flt_mantissa.flt_h_32 += 0x80;
			}
		}
		if (sz == 8 && (e.flt_mantissa.flt_l_32 & 0x80)) {
			/* rounding */
			if ((e.flt_mantissa.flt_l_32 & 0xffffff00) == 0xffffff00) {
				e.flt_mantissa.flt_l_32 = 0;
				if (e.flt_mantissa.flt_h_32 == 0xffffffff) {
					e.flt_exp++;
					e.flt_mantissa.flt_h_32 = 0x80000000;
				}
				else e.flt_mantissa.flt_h_32++;
			}
			else {
				e.flt_mantissa.flt_l_32 += 0x80;
			}
		}
		if (e.flt_exp > 255) {
			overflow = 1;
			e.flt_exp = 255;
			e.flt_mantissa.flt_h_32 = e.flt_mantissa.flt_l_32 = 0xffffffff;
		}
#endif
		buf[I0] = (e.flt_sign << 7) | (e.flt_exp >> 1);
		buf[I1] = ((e.flt_exp&1) << 7) |
				((e.flt_mantissa.flt_h_32 & 0x7fffffff) >> 24);
		buf[I2] = e.flt_mantissa.flt_h_32 >> 16;
		buf[I3] = e.flt_mantissa.flt_h_32 >> 8;
#ifndef IEEEFLOAT
		if (sz == 8) {
			buf[I4] = e.flt_mantissa.flt_h_32;
			buf[I5] = e.flt_mantissa.flt_l_32 >> 24;
			buf[I6] = e.flt_mantissa.flt_l_32 >> 16;
			buf[I7] = e.flt_mantissa.flt_l_32 >> 8;
			flt_b64_sft(&(e.flt_mantissa), -56);
		}
		else
#endif
			flt_b64_sft(&(e.flt_mantissa), -24);
#ifdef IEEEFLOAT
	}
	else {
		e.flt_exp += 1023;
		if (e.flt_mantissa.flt_h_32 == 0) e.flt_exp = 0;
		if (e.flt_mantissa.flt_l_32 & 0x400) {
			/* rounding */
			if ((e.flt_mantissa.flt_l_32 & 0xfffff800) == 0xfffff800) {
				e.flt_mantissa.flt_l_32 = 0;
				if (e.flt_mantissa.flt_h_32 == 0xffffffff) {
					e.flt_exp++;
					e.flt_mantissa.flt_h_32 = 0x80000000;
				}
				else e.flt_mantissa.flt_h_32++;
			}
			else {
				e.flt_mantissa.flt_l_32 += 0x400;
			}
		}
		if (e.flt_exp >= 2047) {
			overflow = 1;
			e.flt_exp = 2047;
			e.flt_mantissa.flt_h_32 = e.flt_mantissa.flt_l_32 = 0;
		}
		if (e.flt_exp <= 0) {
			flt_b64_sft(&(e.flt_mantissa), 1);
			if (e.flt_exp < 0) {
				flt_b64_sft(&(e.flt_mantissa), -e.flt_exp);
				e.flt_exp = 0;
			}
		}
		buf[I0] = (e.flt_sign << 7) | (e.flt_exp >> 4);
		buf[I1] = ((e.flt_exp & 017)<< 4) | ((e.flt_mantissa.flt_h_32 >> 27) & 017);
		buf[I2] = e.flt_mantissa.flt_h_32 >> 19;
		buf[I3] = e.flt_mantissa.flt_h_32 >> 11;
		buf[I4] = e.flt_mantissa.flt_h_32 >> 3;
		buf[I5] = (e.flt_mantissa.flt_h_32 << 5) | ((e.flt_mantissa.flt_l_32 >> 27) & 037);
		buf[I6] = e.flt_mantissa.flt_l_32 >> 19;
		buf[I7] = e.flt_mantissa.flt_l_32 >> 11;
		flt_b64_sft(&(e.flt_mantissa), -53);
	}
#endif
#if ! FL_MSL_AT_LOW_ADDRESS
	if (sz == 4) {
		buf[I4] = buf[I0];
		buf[I5] = buf[I1];
		buf[I6] = buf[I2];
		buf[I7] = buf[I3];
	}
#endif
	if (overflow) {
		return 2;
	}
	return 0;
}
#endif /* USE_FLT */

#ifdef CODE_GENERATOR
con_float()
{
	char buf[8];
	int rval = float_cst(str, (int)argval, buf);
	int i;

	if (rval == 1) {
		fprintf(stderr,"float constant size = %d\n",(int)argval);
		fatal("bad fcon size");
	}
	fprintf(codefile,"!float %s sz %d\n", str, (int)argval);
	if (rval == 2) {
		fprintf(stderr, "Warning: overflow in floating point constant %s\n", str);
	}
	fprintf(codefile, ".data1 0%o", buf[0] & 0377);
	for (i = 1; i < (int)argval; i++) {
		fprintf(codefile, ",0%o", buf[i] & 0377);
	}
	putc('\n', codefile);
}
#endif /* CODE_GENERATOR */

#ifdef CODE_EXPANDER
con_float(str, argval)
	char *str;
	arith argval;
{
	char buf[8];
	int rval = float_cst(str, (int)argval, buf);
	int i;

	if (rval == 1) {
		argval = 8;
		rval = float_cst(str, 8, buf);
	}
	for (i = 0; i < (int)argval; i++) {
		gen1(buf[i]);
	}
}
#endif /* CODE_EXPANDER */