--- /dev/null
+/*
+ (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 */