};
#undef DEFTREECODE
+#define ENUM_BITFIELD(X) enum X
+#define class klass
+
#include "real.h"
+
+#undef class
}
/* We never produce signals from the library. Thus setjmp need do nothing. */
/* ====================================================================== */
/* The implementation of the abstract floating point class based on gcc's
- real.c. I.e. the object of this excersize. Templated so that we can
+ real.c. I.e. the object of this exercise. Templated so that we can
all fp sizes. */
-template<int SIZE, enum machine_mode MODE>
class real_c_float
{
+ public:
+ static const enum machine_mode MODE = SFmode;
+
private:
- long image[SIZE / 32];
+ static const int external_max = 128 / 32;
+ static const int internal_max
+ = (sizeof (REAL_VALUE_TYPE) + sizeof (long) + 1) / sizeof (long);
+ long image[external_max < internal_max ? internal_max : external_max];
void from_long(long);
void from_str(const char *);
void ldexp (int);
};
-template<int SIZE, enum machine_mode MODE>
void
-real_c_float<SIZE, MODE>::from_long (long l)
+real_c_float::from_long (long l)
{
REAL_VALUE_TYPE f;
real_to_target (image, &f, MODE);
}
-template<int SIZE, enum machine_mode MODE>
void
-real_c_float<SIZE, MODE>::from_str (const char *s)
+real_c_float::from_str (const char *s)
{
REAL_VALUE_TYPE f;
- char *p = s;
+ const char *p = s;
if (*p == '-' || *p == '+')
p++;
real_to_target (image, &f, MODE);
}
-template<int SIZE, enum machine_mode MODE>
void
-real_c_float<SIZE, MODE>::binop (int code, const real_c_float &b)
+real_c_float::binop (int code, const real_c_float &b)
{
REAL_VALUE_TYPE ai, bi, ri;
if (verbose)
{
char ab[64], bb[64], rb[64];
- const int digits = int(SIZE / 4);
+ const real_format *fmt = real_format_for_mode[MODE - QFmode];
+ const int digits = (fmt->p * fmt->log2_b + 3) / 4;
char symbol_for_code;
real_from_target (&ri, image, MODE);
- real_to_hexadecimal (ab, &ai, digits);
- real_to_hexadecimal (bb, &bi, digits);
- real_to_hexadecimal (rb, &ri, digits);
+ real_to_hexadecimal (ab, &ai, sizeof(ab), digits, 0);
+ real_to_hexadecimal (bb, &bi, sizeof(bb), digits, 0);
+ real_to_hexadecimal (rb, &ri, sizeof(rb), digits, 0);
switch (code)
{
}
}
-template<int SIZE, enum machine_mode MODE>
void
-real_c_float<SIZE, MODE>::unop (int code)
+real_c_float::unop (int code)
{
REAL_VALUE_TYPE ai, ri;
if (verbose)
{
char ab[64], rb[64];
- const int digits = int(SIZE / 4);
+ const real_format *fmt = real_format_for_mode[MODE - QFmode];
+ const int digits = (fmt->p * fmt->log2_b + 3) / 4;
const char *symbol_for_code;
real_from_target (&ri, image, MODE);
- real_to_hexadecimal (ab, &ai, digits);
- real_to_hexadecimal (rb, &ri, digits);
+ real_to_hexadecimal (ab, &ai, sizeof(ab), digits, 0);
+ real_to_hexadecimal (rb, &ri, sizeof(rb), digits, 0);
switch (code)
{
}
}
-template<int SIZE, enum machine_mode MODE>
bool
-real_c_float<SIZE, MODE>::cmp (int code, const real_c_float &b) const
+real_c_float::cmp (int code, const real_c_float &b) const
{
REAL_VALUE_TYPE ai, bi;
bool ret;
if (verbose)
{
char ab[64], bb[64];
- const int digits = int(SIZE / 4);
+ const real_format *fmt = real_format_for_mode[MODE - QFmode];
+ const int digits = (fmt->p * fmt->log2_b + 3) / 4;
const char *symbol_for_code;
- real_to_hexadecimal (ab, &ai, digits);
- real_to_hexadecimal (bb, &bi, digits);
+ real_to_hexadecimal (ab, &ai, sizeof(ab), digits, 0);
+ real_to_hexadecimal (bb, &bi, sizeof(bb), digits, 0);
switch (code)
{
return ret;
}
-template<int SIZE, enum machine_mode MODE>
const char *
-real_c_float<SIZE, MODE>::str() const
+real_c_float::str() const
{
REAL_VALUE_TYPE f;
- const int digits = int(SIZE * .30102999566398119521 + 1);
+ const real_format *fmt = real_format_for_mode[MODE - QFmode];
+ const int digits = int(fmt->p * fmt->log2_b * .30102999566398119521 + 1);
real_from_target (&f, image, MODE);
char *buf = new char[digits + 10];
- real_to_decimal (buf, &f, digits);
+ real_to_decimal (buf, &f, digits+10, digits, 0);
return buf;
}
-template<int SIZE, enum machine_mode MODE>
const char *
-real_c_float<SIZE, MODE>::hex() const
+real_c_float::hex() const
{
REAL_VALUE_TYPE f;
- const int digits = int(SIZE / 4);
+ const real_format *fmt = real_format_for_mode[MODE - QFmode];
+ const int digits = (fmt->p * fmt->log2_b + 3) / 4;
real_from_target (&f, image, MODE);
char *buf = new char[digits + 10];
- real_to_hexadecimal (buf, &f, digits);
+ real_to_hexadecimal (buf, &f, digits+10, digits, 0);
return buf;
}
-template<int SIZE, enum machine_mode MODE>
long
-real_c_float<SIZE, MODE>::integer() const
+real_c_float::integer() const
{
REAL_VALUE_TYPE f;
real_from_target (&f, image, MODE);
return real_to_integer (&f);
}
-template<int SIZE, enum machine_mode MODE>
int
-real_c_float<SIZE, MODE>::exp() const
+real_c_float::exp() const
{
REAL_VALUE_TYPE f;
real_from_target (&f, image, MODE);
return real_exponent (&f);
}
-template<int SIZE, enum machine_mode MODE>
void
-real_c_float<SIZE, MODE>::ldexp (int exp)
+real_c_float::ldexp (int exp)
{
REAL_VALUE_TYPE ai;
int main(int ac, char **av)
{
- init_real_once ();
+ setbuf(stdout, NULL);
+ setbuf(stderr, NULL);
while (1)
switch (getopt (ac, av, "pvg:fdl"))
break;
case 'g':
{
- int size = strtol (optarg, 0, 0);
-
- switch (size)
- {
- case 32:
- Paranoia< real_c_float<32, SFmode> >().main();
- break;
-
- case 64:
- Paranoia< real_c_float<64, DFmode> >().main();
- break;
-
- case 96:
- Paranoia< real_c_float<96, XFmode> >().main();
- break;
-
- case 128:
- Paranoia< real_c_float<128, TFmode> >().main();
+ static const struct {
+ const char *name;
+ const struct real_format *fmt;
+ } fmts[] = {
+#define F(x) { #x, &x##_format }
+ F(ieee_single),
+ F(ieee_double),
+ F(ieee_extended_motorola),
+ F(ieee_extended_intel_96),
+ F(ieee_extended_intel_128),
+ F(ibm_extended),
+ F(ieee_quad),
+ F(vax_f),
+ F(vax_d),
+ F(vax_g),
+ F(i370_single),
+ F(i370_double),
+ F(real_internal),
+#undef F
+ };
+
+ int i, n = sizeof (fmts)/sizeof(*fmts);
+
+ for (i = 0; i < n; ++i)
+ if (strcmp (fmts[i].name, optarg) == 0)
break;
- default:
- puts ("Invalid gcc implementation size.");
+ if (i == n)
+ {
+ printf ("Unknown implementation \"%s\"; "
+ "available implementations:\n", optarg);
+ for (i = 0; i < n; ++i)
+ printf ("\t%s\n", fmts[i].name);
return 1;
}
+
+ // We cheat and use the same mode all the time, but vary
+ // the format used for that mode.
+ real_format_for_mode[int(real_c_float::MODE) - int(QFmode)]
+ = fmts[i].fmt;
+
+ Paranoia<real_c_float>().main();
break;
}
case '?':
puts ("-p\tpause between pages");
- puts ("-g<N>\treal.c implementation size N");
+ puts ("-g<FMT>\treal.c implementation FMT");
puts ("-f\tnative float");
puts ("-d\tnative double");
puts ("-l\tnative long double");
int target_flags = 0;
-extern "C"
-enum machine_mode
-mode_for_size (unsigned int size, enum mode_class, int)
+extern "C" int
+floor_log2_wide (unsigned HOST_WIDE_INT x)
{
- switch (size)
- {
- case 32:
- return SFmode;
- case 64:
- return DFmode;
- case 96:
- return XFmode;
- case 128:
- return TFmode;
- }
- abort ();
+ int log = -1;
+ while (x != 0)
+ log++,
+ x >>= 1;
+ return log;
}