* [PATCH] The 2nd patch from the fixed-point branch
@ 2007-06-26 1:17 Fu, Chao-Ying
2007-06-26 15:13 ` Richard Sandiford
0 siblings, 1 reply; 4+ messages in thread
From: Fu, Chao-Ying @ 2007-06-26 1:17 UTC (permalink / raw)
To: gcc-patches; +Cc: Mark Mitchell, Thekkath, Radhika, Stephens, Nigel
[-- Attachment #1: Type: text/plain, Size: 751 bytes --]
Hi,
> 2. Merge in fixed-value.h and fixed-value.c to handle fixed-point values.
This is the 2nd patch from the fixed-point branch. The structure to
hold fixed-point values and its associated functions are added.
"double_int_scmp" is fixed to correctly compare two double_int values.
Please review this patch. Thanks a lot!
Regards,
Chao-ying
2007-06-25 Chao-ying Fu <fu@mips.com>
* gengtype.c (main): Handle FIXED_VALUE_TYPE type as scalar typedef.
* double-int.c (double_int_scmp): Use casts of unsigned HOST_WIDE_INT
to compare a.low and b.low.
* fixed-value.c: New file.
* fixed-value.h: New file.
* Makefile.in (OBJS-common): Add fixed-value.o.
(fixed-value.o): New rule.
(GTFILES): Add fixed-value.h.
[-- Attachment #2: gcc.diff --]
[-- Type: application/octet-stream, Size: 39252 bytes --]
Index: gengtype.c
===================================================================
--- gengtype.c (revision 126007)
+++ gengtype.c (working copy)
@@ -3535,6 +3535,7 @@
pos.line = __LINE__ + 1;
do_scalar_typedef ("CUMULATIVE_ARGS", &pos); pos.line++;
do_scalar_typedef ("REAL_VALUE_TYPE", &pos); pos.line++;
+ do_scalar_typedef ("FIXED_VALUE_TYPE", &pos); pos.line++;
do_scalar_typedef ("double_int", &pos); pos.line++;
do_scalar_typedef ("uint8", &pos); pos.line++;
do_scalar_typedef ("jword", &pos); pos.line++;
Index: double-int.c
===================================================================
--- double-int.c (revision 126007)
+++ double-int.c (working copy)
@@ -363,9 +363,9 @@
return -1;
if (a.high > b.high)
return 1;
- if ((HOST_WIDE_INT) a.low < (HOST_WIDE_INT) b.low)
+ if ((unsigned HOST_WIDE_INT) a.low < (unsigned HOST_WIDE_INT) b.low)
return -1;
- if ((HOST_WIDE_INT) a.low > (HOST_WIDE_INT) b.low)
+ if ((unsigned HOST_WIDE_INT) a.low > (unsigned HOST_WIDE_INT) b.low)
return 1;
return 0;
Index: fixed-value.c
===================================================================
--- fixed-value.c (revision 0)
+++ fixed-value.c (revision 0)
@@ -0,0 +1,1144 @@
+/* Fixed-point arithmetic support.
+ Copyright (C) 2006, 2007 Free Software Foundation, Inc.
+
+This file is part of GCC.
+
+GCC is free software; you can redistribute it and/or modify it
+under the terms of the GNU General Public License as published by the
+Free Software Foundation; either version 2, or (at your option) any
+later version.
+
+GCC is distributed in the hope that it will be useful, but WITHOUT
+ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
+for more details.
+
+You should have received a copy of the GNU General Public License
+along with GCC; see the file COPYING. If not, write to the Free
+Software Foundation, 51 Franklin Street, Fifth Floor, Boston, MA
+02110-1301, USA. */
+
+#include "config.h"
+#include "system.h"
+#include "coretypes.h"
+#include "tm.h"
+#include "tree.h"
+#include "toplev.h"
+#include "fixed-value.h"
+
+static int check_real_for_fixed_mode (REAL_VALUE_TYPE *, unsigned int);
+static int get_fixed_sign_bit (double_int, int);
+static bool do_fixed_add (FIXED_VALUE_TYPE *, const FIXED_VALUE_TYPE *,
+ const FIXED_VALUE_TYPE *, int, int);
+static bool do_fixed_multiply (FIXED_VALUE_TYPE *, const FIXED_VALUE_TYPE *,
+ const FIXED_VALUE_TYPE *, int);
+static bool do_fixed_divide (FIXED_VALUE_TYPE *, const FIXED_VALUE_TYPE *,
+ const FIXED_VALUE_TYPE *, int);
+static bool do_fixed_shift (FIXED_VALUE_TYPE *, const FIXED_VALUE_TYPE *,
+ const FIXED_VALUE_TYPE *, int, int);
+static bool do_fixed_neg (FIXED_VALUE_TYPE *, const FIXED_VALUE_TYPE *, int);
+static bool fixed_saturate1 (unsigned int, double_int, double_int *, int);
+static bool fixed_saturate2 (unsigned int, double_int, double_int,
+ double_int *, int);
+
+/* Compare two fixed objects for bitwise identity. */
+
+bool
+fixed_identical (const FIXED_VALUE_TYPE *a, const FIXED_VALUE_TYPE *b)
+{
+ return a->mode == b->mode && a->data.high == b->data.high
+ && a->data.low == b->data.low;
+}
+
+/* Calculate a hash value. */
+
+unsigned int
+fixed_hash (const FIXED_VALUE_TYPE *f)
+{
+ return (unsigned int) (f->data.low ^ f->data.high);
+}
+
+/* Check REAL_VALUE against the range of the fixed-point mode.
+ Return 0, if it is within the range.
+ 1, if it is less than the minimum.
+ 2, if it is greater than the maximum, but not equal to
+ the maximum + the epsilon.
+ 3, if it is equal to the maximum + the epsilon. */
+
+int
+check_real_for_fixed_mode (REAL_VALUE_TYPE *real_value, unsigned int mode)
+{
+ REAL_VALUE_TYPE max_value, min_value, epsilon_value;
+ char max_string[20], min_string[20], epsilon_string[20];
+
+ sprintf (max_string, "0x1.0p%d", GET_MODE_IBIT (mode));
+ sprintf (epsilon_string, "0x1.0p-%d", GET_MODE_FBIT (mode));
+ real_from_string (&max_value, max_string);
+ real_from_string (&epsilon_value, epsilon_string);
+
+ if (SIGNED_FIXED_POINT_MODE_P (mode))
+ {
+ sprintf (min_string, "-0x1.0p%d", GET_MODE_IBIT (mode));
+ real_from_string (&min_value, min_string);
+ }
+ else
+ real_from_string (&min_value, "0.0");
+
+ if (real_compare (LT_EXPR, real_value, &min_value))
+ return 1;
+ if (real_compare (EQ_EXPR, real_value, &max_value))
+ return 3;
+ real_arithmetic (&max_value, MINUS_EXPR, &max_value, &epsilon_value);
+ if (real_compare (GT_EXPR, real_value, &max_value))
+ return 2;
+ return 0;
+}
+
+/* Initialize from a decimal or hexadecimal string. */
+
+void
+fixed_from_string (FIXED_VALUE_TYPE *f, const char *str, tree type)
+{
+ REAL_VALUE_TYPE real_value, fixed_value, base_value;
+ char base_string[20];
+ unsigned int fbit;
+ int temp;
+
+ f->mode = TYPE_MODE (type);
+ fbit = TYPE_FBIT (type);
+
+ sprintf (base_string, "0x1.0p%d", fbit);
+ real_from_string (&real_value, str);
+ temp = check_real_for_fixed_mode (&real_value, f->mode);
+ /* We don't want to warn the case when the _Fract value is 1.0. */
+ if (temp == 1 || temp == 2 || (temp == 3 && ALL_ACCUM_MODE_P (f->mode)))
+ warning (OPT_Woverflow,
+ "large fixed-point constant implicitly truncated to fixed-point type");
+ real_from_string (&base_value, base_string);
+ real_arithmetic (&fixed_value, MULT_EXPR, &real_value, &base_value);
+ real_to_integer2 ((HOST_WIDE_INT *)&f->data.low, &f->data.high,
+ &fixed_value);
+
+ if (temp == 3 && ALL_FRACT_MODE_P (f->mode))
+ {
+ /* From the spec, we need to evaluate 1 to the maximal value. */
+ f->data.low = -1;
+ f->data.high = -1;
+ f->data = double_int_ext (f->data,
+ GET_MODE_FBIT (f->mode)
+ + GET_MODE_IBIT (f->mode), 1);
+ }
+ else
+ f->data = double_int_ext (f->data,
+ SIGNED_FIXED_POINT_MODE_P (f->mode)
+ + GET_MODE_FBIT (f->mode)
+ + GET_MODE_IBIT (f->mode),
+ UNSIGNED_FIXED_POINT_MODE_P (f->mode));
+}
+
+/* Render F as a decimal floating point constant. */
+
+void fixed_to_decimal (char *str, const FIXED_VALUE_TYPE *f_orig,
+ size_t buf_size)
+{
+ char base_string[20];
+ REAL_VALUE_TYPE real_value, base_value, fixed_value;
+
+ sprintf (base_string, "0x1.0p%d", GET_MODE_FBIT (f_orig->mode));
+ real_from_string (&base_value, base_string);
+ real_from_integer (&real_value, VOIDmode, f_orig->data.low, f_orig->data.high,
+ UNSIGNED_FIXED_POINT_MODE_P (f_orig->mode));
+ real_arithmetic (&fixed_value, RDIV_EXPR, &real_value, &base_value);
+ real_to_decimal (str, &fixed_value, buf_size, 0, 1);
+}
+
+/* If SATP, saturate A to the maximum or the minimum, and save to *F based on
+ the machine mode MODE.
+ This function assumes the width of double_int is greater than the width
+ of the fixed-point value at the fixed-point mode.
+ Return true, if !SATP and overflow. */
+
+static bool
+fixed_saturate1 (unsigned int mode, double_int a, double_int *f, int satp)
+{
+ bool overflow = false;
+ int unsignedp = UNSIGNED_FIXED_POINT_MODE_P (mode);
+ int i_f_bits = GET_MODE_IBIT (mode) + GET_MODE_FBIT (mode);
+
+ if (unsignedp) /* Unsigned type. */
+ {
+ double_int max;
+ max.low = -1;
+ max.high = -1;
+ max = double_int_ext (max, i_f_bits, 1);
+ if (double_int_cmp (a, max, 1) == 1)
+ {
+ if (satp)
+ *f = max;
+ else
+ overflow = true;
+ }
+ }
+ else /* Signed type. */
+ {
+ double_int max, min;
+ max.high = -1;
+ max.low = -1;
+ max = double_int_ext (max, i_f_bits, 1);
+ min.high = 0;
+ min.low = 1;
+ lshift_double (min.low, min.high, i_f_bits,
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &min.low, &min.high, 1);
+ min = double_int_ext (min, 1 + i_f_bits, 0);
+ if (double_int_cmp (a, max, 0) == 1)
+ {
+ if (satp)
+ *f = max;
+ else
+ overflow = true;
+ }
+ else if (double_int_cmp (a, min, 0) == -1)
+ {
+ if (satp)
+ *f = min;
+ else
+ overflow = true;
+ }
+ }
+ return overflow;
+}
+
+/* If SATP, saturate {A_HIGH, A_LOW} to the maximum or the minimum, and
+ save to *F based on the machine mode MODE.
+ This function assumes the width of two double_int is greater than the width
+ of the fixed-point value at the fixed-point mode.
+ Return true, if !SATP and overflow. */
+
+static bool
+fixed_saturate2 (unsigned int mode, double_int a_high, double_int a_low,
+ double_int *f, int satp)
+{
+ bool overflow = false;
+ int unsignedp = UNSIGNED_FIXED_POINT_MODE_P (mode);
+ int i_f_bits = GET_MODE_IBIT (mode) + GET_MODE_FBIT (mode);
+
+ if (unsignedp) /* Unsigned type. */
+ {
+ double_int max_r, max_s;
+ max_r.high = 0;
+ max_r.low = 0;
+ max_s.high = -1;
+ max_s.low = -1;
+ max_s = double_int_ext (max_s, i_f_bits, 1);
+ if (double_int_cmp (a_high, max_r, 1) == 1
+ || (double_int_equal_p (a_high, max_r) &&
+ double_int_cmp (a_low, max_s, 1) == 1))
+ {
+ if (satp)
+ *f = max_s;
+ else
+ overflow = true;
+ }
+ }
+ else /* Signed type. */
+ {
+ double_int max_r, max_s, min_r, min_s;
+ max_r.high = 0;
+ max_r.low = 0;
+ max_s.high = -1;
+ max_s.low = -1;
+ max_s = double_int_ext (max_s, i_f_bits, 1);
+ min_r.high = -1;
+ min_r.low = -1;
+ min_s.high = 0;
+ min_s.low = 1;
+ lshift_double (min_s.low, min_s.high, i_f_bits,
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &min_s.low, &min_s.high, 1);
+ min_s = double_int_ext (min_s, 1 + i_f_bits, 0);
+ if (double_int_cmp (a_high, max_r, 0) == 1
+ || (double_int_equal_p (a_high, max_r) &&
+ double_int_cmp (a_low, max_s, 1) == 1))
+ {
+ if (satp)
+ *f = max_s;
+ else
+ overflow = true;
+ }
+ else if (double_int_cmp (a_high, min_r, 0) == -1
+ || (double_int_equal_p (a_high, min_r) &&
+ double_int_cmp (a_low, min_s, 1) == -1))
+ {
+ if (satp)
+ *f = min_s;
+ else
+ overflow = true;
+ }
+ }
+ return overflow;
+}
+
+/* Return the sign bit based on I_F_BITS. */
+
+inline int get_fixed_sign_bit (double_int a, int i_f_bits)
+{
+ if (i_f_bits < HOST_BITS_PER_WIDE_INT)
+ return (a.low >> i_f_bits) & 1;
+ else
+ return (a.high >> (i_f_bits - HOST_BITS_PER_WIDE_INT)) & 1;
+}
+
+/* Calculate F = A + (SUBTRACT_P ? -B : B).
+ If SATP, saturate the result to the max or the min.
+ Return true, if !SATP and overflow. */
+
+static bool
+do_fixed_add (FIXED_VALUE_TYPE *f, const FIXED_VALUE_TYPE *a,
+ const FIXED_VALUE_TYPE *b, int subtract_p, int satp)
+{
+ bool overflow = false;
+ double_int temp = subtract_p ? double_int_neg (b->data) : b->data;
+ int unsignedp = UNSIGNED_FIXED_POINT_MODE_P (a->mode);
+ int i_f_bits = GET_MODE_IBIT (a->mode) + GET_MODE_FBIT (a->mode);
+ f->mode = a->mode;
+ f->data = double_int_add (a->data, temp);
+ if (unsignedp) /* Unsigned type. */
+ {
+ if (subtract_p) /* Unsigned subtraction. */
+ {
+ if (double_int_cmp (a->data, b->data, 1) == -1)
+ {
+ if (satp)
+ {
+ f->data.high = 0;
+ f->data.low = 0;
+ }
+ else
+ overflow = true;
+ }
+ }
+ else /* Unsigned addition. */
+ {
+ f->data = double_int_ext (f->data, i_f_bits, 1);
+ if (double_int_cmp (f->data, a->data, 1) == -1
+ || double_int_cmp (f->data, b->data, 1) == -1)
+ {
+ if (satp)
+ {
+ f->data.high = -1;
+ f->data.low = -1;
+ }
+ else
+ overflow = true;
+ }
+ }
+ }
+ else /* Signed type. */
+ {
+ if (((!subtract_p)
+ && get_fixed_sign_bit (a->data, i_f_bits)
+ == get_fixed_sign_bit (b->data, i_f_bits)
+ && get_fixed_sign_bit (a->data, i_f_bits)
+ != get_fixed_sign_bit (f->data, i_f_bits))
+ || (subtract_p
+ && get_fixed_sign_bit (a->data, i_f_bits)
+ != get_fixed_sign_bit (b->data, i_f_bits)
+ && get_fixed_sign_bit (a->data, i_f_bits)
+ != get_fixed_sign_bit (f->data, i_f_bits)))
+ {
+ if (satp)
+ {
+ f->data.low = 1;
+ f->data.high = 0;
+ lshift_double (f->data.low, f->data.high, i_f_bits,
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &f->data.low, &f->data.high, 1);
+ if (get_fixed_sign_bit (a->data, i_f_bits) == 0)
+ {
+ double_int one;
+ one.low = 1;
+ one.high = 0;
+ f->data = double_int_add (f->data, double_int_neg (one));
+ }
+ }
+ else
+ overflow = true;
+ }
+ }
+ f->data = double_int_ext (f->data, (!unsignedp) + i_f_bits, unsignedp);
+ return overflow;
+}
+
+/* Calculate F = A * B.
+ If SATP, saturate the result to the max or the min.
+ Return true, if !SATP and overflow. */
+
+static bool
+do_fixed_multiply (FIXED_VALUE_TYPE *f, const FIXED_VALUE_TYPE *a,
+ const FIXED_VALUE_TYPE *b, int satp)
+{
+ bool overflow = false;
+ int unsignedp = UNSIGNED_FIXED_POINT_MODE_P (a->mode);
+ int i_f_bits = GET_MODE_IBIT (a->mode) + GET_MODE_FBIT (a->mode);
+ f->mode = a->mode;
+ if (GET_MODE_PRECISION (f->mode) <= HOST_BITS_PER_WIDE_INT)
+ {
+ f->data = double_int_mul (a->data, b->data);
+ lshift_double (f->data.low, f->data.high,
+ (-GET_MODE_FBIT (f->mode)),
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &f->data.low, &f->data.high, !unsignedp);
+ overflow = fixed_saturate1 (f->mode, f->data, &f->data, satp);
+ }
+ else
+ {
+ /* The result of multiplication expands to two double_int. */
+ double_int a_high, a_low, b_high, b_low;
+ double_int high_high, high_low, low_high, low_low;
+ double_int r, s, temp1, temp2;
+ int carry = 0;
+
+ /* Decompose a and b to four double_int. */
+ a_high.low = a->data.high;
+ a_high.high = 0;
+ a_low.low = a->data.low;
+ a_low.high = 0;
+ b_high.low = b->data.high;
+ b_high.high = 0;
+ b_low.low = b->data.low;
+ b_low.high = 0;
+
+ /* Perform four multiplications. */
+ low_low = double_int_mul (a_low, b_low);
+ low_high = double_int_mul (a_low, b_high);
+ high_low = double_int_mul (a_high, b_low);
+ high_high = double_int_mul (a_high, b_high);
+
+ /* Accumulate four results to {r, s}. */
+ temp1.high = high_low.low;
+ temp1.low = 0;
+ s = double_int_add (low_low, temp1);
+ if (double_int_cmp (s, low_low, 1) == -1
+ || double_int_cmp (s, temp1, 1) == -1)
+ carry ++; /* Carry */
+ temp1.high = s.high;
+ temp1.low = s.low;
+ temp2.high = low_high.low;
+ temp2.low = 0;
+ s = double_int_add (temp1, temp2);
+ if (double_int_cmp (s, temp1, 1) == -1
+ || double_int_cmp (s, temp2, 1) == -1)
+ carry ++; /* Carry */
+
+ temp1.low = high_low.high;
+ temp1.high = 0;
+ r = double_int_add (high_high, temp1);
+ temp1.low = low_high.high;
+ temp1.high = 0;
+ r = double_int_add (r, temp1);
+ temp1.low = carry;
+ temp1.high = 0;
+ r = double_int_add (r, temp1);
+
+ /* We need to add neg(b) to r, if a < 0. */
+ if (!unsignedp && a->data.high < 0)
+ r = double_int_add (r, double_int_neg (b->data));
+ /* We need to add neg(a) to r, if b < 0. */
+ if (!unsignedp && b->data.high < 0)
+ r = double_int_add (r, double_int_neg (a->data));
+
+ /* Shift right the result by FBIT. */
+ if (GET_MODE_FBIT (f->mode) == 2 * HOST_BITS_PER_WIDE_INT)
+ {
+ s.low = r.low;
+ s.high = r.high;
+ if (unsignedp)
+ {
+ r.low = 0;
+ r.high = 0;
+ }
+ else
+ {
+ r.low = -1;
+ r.high = -1;
+ }
+ f->data.low = s.low;
+ f->data.high = s.high;
+ }
+ else
+ {
+ lshift_double (s.low, s.high,
+ (-GET_MODE_FBIT (f->mode)),
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &s.low, &s.high, 0);
+ lshift_double (r.low, r.high,
+ (2 * HOST_BITS_PER_WIDE_INT - GET_MODE_FBIT (f->mode)),
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &f->data.low, &f->data.high, 0);
+ f->data.low = f->data.low | s.low;
+ f->data.high = f->data.high | s.high;
+ s.low = f->data.low;
+ s.high = f->data.high;
+ lshift_double (r.low, r.high,
+ (-GET_MODE_FBIT (f->mode)),
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &r.low, &r.high, !unsignedp);
+ }
+
+ overflow = fixed_saturate2 (f->mode, r, s, &f->data, satp);
+ }
+
+ f->data = double_int_ext (f->data, (!unsignedp) + i_f_bits, unsignedp);
+ return false;
+}
+
+/* Calculate F = A / B.
+ If SATP, saturate the result to the max or the min.
+ Return true, if !SATP and overflow. */
+
+static bool
+do_fixed_divide (FIXED_VALUE_TYPE *f, const FIXED_VALUE_TYPE *a,
+ const FIXED_VALUE_TYPE *b, int satp)
+{
+ bool overflow = false;
+ int unsignedp = UNSIGNED_FIXED_POINT_MODE_P (a->mode);
+ int i_f_bits = GET_MODE_IBIT (a->mode) + GET_MODE_FBIT (a->mode);
+ f->mode = a->mode;
+ if (GET_MODE_PRECISION (f->mode) <= HOST_BITS_PER_WIDE_INT)
+ {
+ lshift_double (a->data.low, a->data.high,
+ GET_MODE_FBIT (f->mode),
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &f->data.low, &f->data.high, !unsignedp);
+ f->data = double_int_div (f->data, b->data, unsignedp, TRUNC_DIV_EXPR);
+ overflow = fixed_saturate1 (f->mode, f->data, &f->data, satp);
+ }
+ else
+ {
+ double_int pos_a, pos_b, r, s;
+ double_int quo_r, quo_s, mod, temp;
+ int num_of_neg = 0;
+ int i;
+
+ /* If a < 0, negate a. */
+ if (!unsignedp && a->data.high < 0)
+ {
+ pos_a = double_int_neg (a->data);
+ num_of_neg ++;
+ }
+ else
+ pos_a = a->data;
+
+ /* If b < 0, negate b. */
+ if (!unsignedp && b->data.high < 0)
+ {
+ pos_b = double_int_neg (b->data);
+ num_of_neg ++;
+ }
+ else
+ pos_b = b->data;
+
+ /* Left shift pos_a to {r, s} by FBIT. */
+ if (GET_MODE_FBIT (f->mode) == 2 * HOST_BITS_PER_WIDE_INT)
+ {
+ r = pos_a;
+ s.high = 0;
+ s.low = 0;
+ }
+ else
+ {
+ lshift_double (pos_a.low, pos_a.high,
+ GET_MODE_FBIT (f->mode),
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &s.low, &s.high, 0);
+ lshift_double (pos_a.low, pos_a.high,
+ - (2 * HOST_BITS_PER_WIDE_INT - GET_MODE_FBIT (f->mode)),
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &r.low, &r.high, 0);
+ }
+
+ /* Divide r by pos_b to quo_r. The remanider is in mod. */
+ div_and_round_double (TRUNC_DIV_EXPR, 1, r.low, r.high, pos_b.low,
+ pos_b.high, &quo_r.low, &quo_r.high, &mod.low,
+ &mod.high);
+
+ quo_s.high = 0;
+ quo_s.low = 0;
+
+ for (i = 0; i < 2 * HOST_BITS_PER_WIDE_INT; i++)
+ {
+ /* Record the leftmost bit of mod. */
+ int leftmost_mod = (mod.high < 0);
+
+ /* Shift left mod by 1 bit. */
+ lshift_double (mod.low, mod.high, 1, 2 * HOST_BITS_PER_WIDE_INT,
+ &mod.low, &mod.high, 0);
+
+ /* Test the leftmost bit of s to add to mod. */
+ if (s.high < 0)
+ mod.low += 1;
+
+ /* Shift left quo_s by 1 bit. */
+ lshift_double (quo_s.low, quo_s.high, 1, 2 * HOST_BITS_PER_WIDE_INT,
+ &quo_s.low, &quo_s.high, 0);
+
+ /* Try to calculate (mod - pos_b). */
+ temp = double_int_add (mod, double_int_neg (pos_b));
+
+ if (leftmost_mod == 1 || double_int_cmp (mod, pos_b, 1) != -1)
+ {
+ quo_s.low += 1;
+ mod = temp;
+ }
+
+ /* Shift left s by 1 bit. */
+ lshift_double (s.low, s.high, 1, 2 * HOST_BITS_PER_WIDE_INT,
+ &s.low, &s.high, 0);
+
+ }
+
+ if (num_of_neg == 1)
+ {
+ quo_s = double_int_neg (quo_s);
+ if (quo_s.high == 0 && quo_s.low == 0)
+ quo_r = double_int_neg (quo_r);
+ else
+ {
+ quo_r.low = ~quo_r.low;
+ quo_r.high = ~quo_r.high;
+ }
+ }
+
+ f->data = quo_s;
+ overflow = fixed_saturate2 (f->mode, quo_r, quo_s, &f->data, satp);
+ }
+
+ f->data = double_int_ext (f->data, (!unsignedp) + i_f_bits, unsignedp);
+ return overflow;
+}
+
+/* Calculate F = A << B if LEFT_P. Otherwies, F = A >> B.
+ If SATP, saturate the result to the max or the min.
+ Return true, if !SATP and overflow. */
+
+static bool
+do_fixed_shift (FIXED_VALUE_TYPE *f, const FIXED_VALUE_TYPE *a,
+ const FIXED_VALUE_TYPE *b, int left_p, int satp)
+{
+ bool overflow = false;
+ int unsignedp = UNSIGNED_FIXED_POINT_MODE_P (a->mode);
+ int i_f_bits = GET_MODE_IBIT (a->mode) + GET_MODE_FBIT (a->mode);
+ f->mode = a->mode;
+
+ if (b->data.low == 0)
+ {
+ f->data = a->data;
+ return overflow;
+ }
+
+ if (GET_MODE_PRECISION (f->mode) <= HOST_BITS_PER_WIDE_INT || (!left_p))
+ {
+ lshift_double (a->data.low, a->data.high,
+ left_p ? b->data.low : (-b->data.low),
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &f->data.low, &f->data.high, !unsignedp);
+ if (left_p) /* Only left shift saturates. */
+ overflow = fixed_saturate1 (f->mode, f->data, &f->data, satp);
+ }
+ else /* We need two double_int to store the left-shift result. */
+ {
+ double_int temp_high, temp_low;
+ if (b->data.low == 2 * HOST_BITS_PER_WIDE_INT)
+ {
+ temp_high = a->data;
+ temp_low.high = 0;
+ temp_low.low = 0;
+ }
+ else
+ {
+ lshift_double (a->data.low, a->data.high,
+ b->data.low,
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &temp_low.low, &temp_low.high, !unsignedp);
+ /* Logical shift right to temp_high. */
+ lshift_double (a->data.low, a->data.high,
+ b->data.low - 2 * HOST_BITS_PER_WIDE_INT,
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &temp_high.low, &temp_high.high, 0);
+ }
+ if (!unsignedp && a->data.high < 0) /* Signed-extend temp_high. */
+ temp_high = double_int_ext (temp_high, b->data.low, unsignedp);
+ f->data = temp_low;
+ overflow = fixed_saturate2 (f->mode, temp_high, temp_low, &f->data, satp);
+ }
+ f->data = double_int_ext (f->data, (!unsignedp) + i_f_bits, unsignedp);
+ return overflow;
+}
+
+/* Calculate F = -A.
+ If SATP, saturate the result to the max or the min.
+ Return true, if !SATP and overflow. */
+
+static bool
+do_fixed_neg (FIXED_VALUE_TYPE *f, const FIXED_VALUE_TYPE *a, int satp)
+{
+ bool overflow = false;
+ int unsignedp = UNSIGNED_FIXED_POINT_MODE_P (a->mode);
+ int i_f_bits = GET_MODE_IBIT (a->mode) + GET_MODE_FBIT (a->mode);
+ f->mode = a->mode;
+ f->data = double_int_neg (a->data);
+ f->data = double_int_ext (f->data, (!unsignedp) + i_f_bits, unsignedp);
+
+ if (unsignedp) /* Unsigned type. */
+ {
+ if (f->data.low != 0 || f->data.high != 0)
+ {
+ if (satp)
+ {
+ f->data.low = 0;
+ f->data.high = 0;
+ }
+ else
+ overflow = true;
+ }
+ }
+ else /* Signed type. */
+ {
+ if (!(f->data.high == 0 && f->data.low == 0)
+ && f->data.high == a->data.high && f->data.low == a->data.low )
+ {
+ if (satp)
+ {
+ /* Saturate to the maximum by subtracting f->data by one. */
+ f->data.low = -1;
+ f->data.high = -1;
+ f->data = double_int_ext (f->data, i_f_bits, 1);
+ }
+ else
+ overflow = true;
+ }
+ }
+ return overflow;
+}
+
+/* Perform the binary or unary operation described by CODE.
+ For a unary operation, leave OP1 NULL.
+ Return true, if !SATP and overflow. */
+
+bool
+fixed_arithmetic (FIXED_VALUE_TYPE *f, int icode, const FIXED_VALUE_TYPE *op0,
+ const FIXED_VALUE_TYPE *op1, int satp)
+{
+ switch (icode)
+ {
+ case NEGATE_EXPR:
+ return do_fixed_neg (f, op0, satp);
+ break;
+
+ case PLUS_EXPR:
+ return do_fixed_add (f, op0, op1, 0, satp);
+ break;
+
+ case MINUS_EXPR:
+ return do_fixed_add (f, op0, op1, 1, satp);
+ break;
+
+ case MULT_EXPR:
+ return do_fixed_multiply (f, op0, op1, satp);
+ break;
+
+ case TRUNC_DIV_EXPR:
+ return do_fixed_divide (f, op0, op1, satp);
+ break;
+
+ case LSHIFT_EXPR:
+ return do_fixed_shift (f, op0, op1, 1, satp);
+ break;
+
+ case RSHIFT_EXPR:
+ return do_fixed_shift (f, op0, op1, 0, satp);
+ break;
+
+ default:
+ gcc_unreachable ();
+ }
+ return false;
+}
+
+/* Compare fixed-point values by tree_code. */
+
+bool
+fixed_compare (int icode, const FIXED_VALUE_TYPE *op0,
+ const FIXED_VALUE_TYPE *op1)
+{
+ enum tree_code code = icode;
+
+ switch (code)
+ {
+ case NE_EXPR:
+ return op0->mode != op1->mode
+ || !double_int_equal_p (op0->data, op1->data);
+
+ case EQ_EXPR:
+ return op0->mode == op1->mode
+ && double_int_equal_p (op0->data, op1->data);
+
+ case LT_EXPR:
+ return op0->mode == op1->mode
+ && double_int_cmp (op0->data, op1->data,
+ UNSIGNED_FIXED_POINT_MODE_P (op0->mode)) == -1;
+
+ case LE_EXPR:
+ return op0->mode == op1->mode
+ && double_int_cmp (op0->data, op1->data,
+ UNSIGNED_FIXED_POINT_MODE_P (op0->mode)) != 1;
+
+ case GT_EXPR:
+ return op0->mode == op1->mode
+ && double_int_cmp (op0->data, op1->data,
+ UNSIGNED_FIXED_POINT_MODE_P (op0->mode)) == 1;
+
+ case GE_EXPR:
+ return op0->mode == op1->mode
+ && double_int_cmp (op0->data, op1->data,
+ UNSIGNED_FIXED_POINT_MODE_P (op0->mode)) != -1;
+
+ default:
+ gcc_unreachable ();
+ }
+}
+
+/* Extend or truncate to a new mode.
+ If SATP, saturate the result to the max or the min.
+ Return true, if !SATP and overflow. */
+
+bool
+fixed_convert (FIXED_VALUE_TYPE *f, enum machine_mode mode,
+ const FIXED_VALUE_TYPE *a, int satp)
+{
+ bool overflow = false;
+ if (mode == a->mode)
+ {
+ *f = *a;
+ return overflow;
+ }
+
+ if (GET_MODE_FBIT (mode) > GET_MODE_FBIT (a->mode))
+ {
+ /* Left shift a to temp_high, temp_low based on a->mode. */
+ double_int temp_high, temp_low;
+ int amount = GET_MODE_FBIT (mode) - GET_MODE_FBIT (a->mode);
+ lshift_double (a->data.low, a->data.high,
+ amount,
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &temp_low.low, &temp_low.high,
+ SIGNED_FIXED_POINT_MODE_P (a->mode));
+ /* Logical shift right to temp_high. */
+ lshift_double (a->data.low, a->data.high,
+ amount - 2 * HOST_BITS_PER_WIDE_INT,
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &temp_high.low, &temp_high.high, 0);
+ if (SIGNED_FIXED_POINT_MODE_P (a->mode)
+ && a->data.high < 0) /* Signed-extend temp_high. */
+ temp_high = double_int_ext (temp_high, amount, 0);
+ f->mode = mode;
+ f->data = temp_low;
+ if (SIGNED_FIXED_POINT_MODE_P (a->mode) ==
+ SIGNED_FIXED_POINT_MODE_P (f->mode))
+ overflow = fixed_saturate2 (f->mode, temp_high, temp_low, &f->data,
+ satp);
+ else
+ {
+ /* Take care of the cases when converting between signed and
+ unsigned. */
+ if (SIGNED_FIXED_POINT_MODE_P (a->mode))
+ {
+ /* Signed -> Unsigned. */
+ if (a->data.high < 0)
+ {
+ if (satp)
+ {
+ f->data.low = 0; /* Set to zero. */
+ f->data.high = 0; /* Set to zero. */
+ }
+ else
+ overflow = true;
+ }
+ else
+ overflow = fixed_saturate2 (f->mode, temp_high, temp_low,
+ &f->data, satp);
+ }
+ else
+ {
+ /* Unsigned -> Signed. */
+ if (temp_high.high < 0)
+ {
+ if (satp)
+ {
+ /* Set to maximum. */
+ f->data.low = -1; /* Set to all ones. */
+ f->data.high = -1; /* Set to all ones. */
+ f->data = double_int_ext (f->data,
+ GET_MODE_FBIT (f->mode) +
+ GET_MODE_IBIT (f->mode),
+ 1); /* Clear the sign. */
+ }
+ else
+ overflow = true;
+ }
+ else
+ overflow = fixed_saturate2 (f->mode, temp_high, temp_low,
+ &f->data, satp);
+ }
+ }
+ }
+ else
+ {
+ /* Right shift a to temp based on a->mode. */
+ double_int temp;
+ lshift_double (a->data.low, a->data.high,
+ GET_MODE_FBIT (mode) - GET_MODE_FBIT (a->mode),
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &temp.low, &temp.high,
+ SIGNED_FIXED_POINT_MODE_P (a->mode));
+ f->mode = mode;
+ f->data = temp;
+ if (SIGNED_FIXED_POINT_MODE_P (a->mode) ==
+ SIGNED_FIXED_POINT_MODE_P (f->mode))
+ overflow = fixed_saturate1 (f->mode, f->data, &f->data, satp);
+ else
+ {
+ /* Take care of the cases when converting between signed and
+ unsigned. */
+ if (SIGNED_FIXED_POINT_MODE_P (a->mode))
+ {
+ /* Signed -> Unsigned. */
+ if (a->data.high < 0)
+ {
+ if (satp)
+ {
+ f->data.low = 0; /* Set to zero. */
+ f->data.high = 0; /* Set to zero. */
+ }
+ else
+ overflow = true;
+ }
+ else
+ overflow = fixed_saturate1 (f->mode, f->data, &f->data, satp);
+ }
+ else
+ {
+ /* Unsigned -> Signed. */
+ if (temp.high < 0)
+ {
+ if (satp)
+ {
+ /* Set to maximum. */
+ f->data.low = -1; /* Set to all ones. */
+ f->data.high = -1; /* Set to all ones. */
+ f->data = double_int_ext (f->data,
+ GET_MODE_FBIT (f->mode) +
+ GET_MODE_IBIT (f->mode),
+ 1); /* Clear the sign. */
+ }
+ else
+ overflow = true;
+ }
+ else
+ overflow = fixed_saturate1 (f->mode, f->data, &f->data, satp);
+ }
+ }
+ }
+
+ f->data = double_int_ext (f->data,
+ SIGNED_FIXED_POINT_MODE_P (f->mode)
+ + GET_MODE_FBIT (f->mode)
+ + GET_MODE_IBIT (f->mode),
+ UNSIGNED_FIXED_POINT_MODE_P (f->mode));
+ return overflow;
+}
+
+/* Convert to a new fixed-point mode from an integer.
+ If UNSIGNEDP, this integer is unsigned.
+ If SATP, saturate the result to the max or the min.
+ Return true, if !SATP and overflow. */
+
+bool
+fixed_convert_from_int (FIXED_VALUE_TYPE *f, enum machine_mode mode,
+ double_int a, int unsignedp, int satp)
+{
+ bool overflow = false;
+ /* Left shift a to temp_high, temp_low. */
+ double_int temp_high, temp_low;
+ int amount = GET_MODE_FBIT (mode);
+ if (amount == 2 * HOST_BITS_PER_WIDE_INT)
+ {
+ temp_high = a;
+ temp_low.low = 0;
+ temp_low.high = 0;
+ }
+ else
+ {
+ lshift_double (a.low, a.high,
+ amount,
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &temp_low.low, &temp_low.high, 0);
+
+ /* Logical shift right to temp_high. */
+ lshift_double (a.low, a.high,
+ amount - 2 * HOST_BITS_PER_WIDE_INT,
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &temp_high.low, &temp_high.high, 0);
+ }
+ if (!unsignedp && a.high < 0) /* Signed-extend temp_high. */
+ temp_high = double_int_ext (temp_high, amount, 0);
+
+ f->mode = mode;
+ f->data = temp_low;
+
+ if (unsignedp == UNSIGNED_FIXED_POINT_MODE_P (f->mode))
+ overflow = fixed_saturate2 (f->mode, temp_high, temp_low, &f->data, satp);
+ else
+ {
+ /* Take care of the cases when converting between signed and unsigned. */
+ if (!unsignedp)
+ {
+ /* Signed -> Unsigned. */
+ if (a.high < 0)
+ {
+ if (satp)
+ {
+ f->data.low = 0; /* Set to zero. */
+ f->data.high = 0; /* Set to zero. */
+ }
+ else
+ overflow = true;
+ }
+ else
+ overflow = fixed_saturate2 (f->mode, temp_high, temp_low, &f->data,
+ satp);
+ }
+ else
+ {
+ /* Unsigned -> Signed. */
+ if (temp_high.high < 0)
+ {
+ if (satp)
+ {
+ /* Set to maximum. */
+ f->data.low = -1; /* Set to all ones. */
+ f->data.high = -1; /* Set to all ones. */
+ f->data = double_int_ext (f->data,
+ GET_MODE_FBIT (f->mode) +
+ GET_MODE_IBIT (f->mode),
+ 1); /* Clear the sign. */
+ }
+ else
+ overflow = true;
+ }
+ else
+ overflow = fixed_saturate2 (f->mode, temp_high, temp_low,
+ &f->data, satp);
+ }
+ }
+ f->data = double_int_ext (f->data,
+ SIGNED_FIXED_POINT_MODE_P (f->mode)
+ + GET_MODE_FBIT (f->mode)
+ + GET_MODE_IBIT (f->mode),
+ UNSIGNED_FIXED_POINT_MODE_P (f->mode));
+ return overflow;
+}
+
+/* Convert to a new fixed-point mode from a real.
+ If SATP, saturate the result to the max or the min.
+ Return true, if !SATP and overflow. */
+
+bool
+fixed_convert_from_real (FIXED_VALUE_TYPE *f, enum machine_mode mode,
+ const REAL_VALUE_TYPE *a, int satp)
+{
+ bool overflow = false;
+ REAL_VALUE_TYPE real_value, fixed_value, base_value;
+ char base_string[20];
+ int unsignedp = UNSIGNED_FIXED_POINT_MODE_P (mode);
+ int i_f_bits = GET_MODE_IBIT (mode) + GET_MODE_FBIT (mode);
+ unsigned int fbit = GET_MODE_FBIT (mode);
+ int temp;
+
+ real_value = *a;
+ f->mode = mode;
+ sprintf (base_string, "0x1.0p%d", fbit);
+ real_from_string (&base_value, base_string);
+ real_arithmetic (&fixed_value, MULT_EXPR, &real_value, &base_value);
+ real_to_integer2 ((HOST_WIDE_INT *)&f->data.low, &f->data.high, &fixed_value);
+ temp = check_real_for_fixed_mode (&real_value, mode);
+ if (temp == 1) /* Minimum. */
+ {
+ if (satp)
+ {
+ if (unsignedp)
+ {
+ f->data.low = 0;
+ f->data.high = 0;
+ }
+ else
+ {
+ f->data.low = 1;
+ f->data.high = 0;
+ lshift_double (f->data.low, f->data.high, i_f_bits,
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &f->data.low, &f->data.high, 1);
+ f->data = double_int_ext (f->data, 1 + i_f_bits, 0);
+ }
+ }
+ else
+ overflow = true;
+ }
+ else if (temp == 2 || temp == 3) /* Maximum. */
+ {
+ if (satp)
+ {
+ f->data.low = -1;
+ f->data.high = -1;
+ f->data = double_int_ext (f->data, i_f_bits, 1);
+ }
+ else
+ overflow = true;
+ }
+ f->data = double_int_ext (f->data, (!unsignedp) + i_f_bits, unsignedp);
+ return overflow;
+}
+
+/* Convert to a new real mode from a fixed-point. */
+
+void
+real_convert_from_fixed (REAL_VALUE_TYPE *r, enum machine_mode mode,
+ const FIXED_VALUE_TYPE *f)
+{
+ char base_string[20];
+ REAL_VALUE_TYPE base_value, fixed_value, real_value;
+
+ sprintf (base_string, "0x1.0p%d", GET_MODE_FBIT (f->mode));
+ real_from_string (&base_value, base_string);
+ real_from_integer (&fixed_value, VOIDmode, f->data.low, f->data.high,
+ UNSIGNED_FIXED_POINT_MODE_P (f->mode));
+ real_arithmetic (&real_value, RDIV_EXPR, &fixed_value, &base_value);
+ real_convert (r, mode, &real_value);
+}
+
+/* Determine whether a fixed-point value F is negative. */
+
+bool
+fixed_isneg (const FIXED_VALUE_TYPE *f)
+{
+ if (SIGNED_FIXED_POINT_MODE_P (f->mode))
+ {
+ int i_f_bits = GET_MODE_IBIT (f->mode) + GET_MODE_FBIT (f->mode);
+ int sign_bit = get_fixed_sign_bit (f->data, i_f_bits);
+ if (sign_bit == 1)
+ return true;
+ }
+
+ return false;
+}
Index: fixed-value.h
===================================================================
--- fixed-value.h (revision 0)
+++ fixed-value.h (revision 0)
@@ -0,0 +1,92 @@
+/* Fixed-point arithmetic support.
+ Copyright (C) 2006, 2007 Free Software Foundation, Inc.
+
+This file is part of GCC.
+
+GCC is free software; you can redistribute it and/or modify it
+under the terms of the GNU General Public License as published by the
+Free Software Foundation; either version 2, or (at your option) any
+later version.
+
+GCC is distributed in the hope that it will be useful, but WITHOUT
+ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
+for more details.
+
+You should have received a copy of the GNU General Public License
+along with GCC; see the file COPYING. If not, write to the Free
+Software Foundation, 51 Franklin Street, Fifth Floor, Boston, MA
+02110-1301, USA. */
+
+#ifndef GCC_FIXED_VALUE_H
+#define GCC_FIXED_VALUE_H
+
+#include "machmode.h"
+#include "real.h"
+#include "double-int.h"
+
+struct fixed_value GTY(())
+{
+ double_int data; /* Store data up to 2 wide integers. */
+ unsigned int mode; /* Use machine mode to know IBIT and FBIT. */
+};
+
+#define FIXED_VALUE_TYPE struct fixed_value
+
+/* Constant fixed-point values 0 and 1. */
+extern FIXED_VALUE_TYPE fconst0[];
+extern FIXED_VALUE_TYPE fconst1[];
+
+/* Return a CONST_FIXED with value R and mode M. */
+#define CONST_FIXED_FROM_FIXED_VALUE(r, m) \
+ const_fixed_from_fixed_value (r, m)
+extern rtx const_fixed_from_fixed_value (FIXED_VALUE_TYPE, enum machine_mode);
+
+/* Initialize from a decimal or hexadecimal string. */
+extern void fixed_from_string (FIXED_VALUE_TYPE *, const char *, tree);
+
+/* In tree.c: wrap up a FIXED_VALUE_TYPE in a tree node. */
+extern tree build_fixed (tree, FIXED_VALUE_TYPE);
+
+/* Extend or truncate to a new mode. */
+extern bool fixed_convert (FIXED_VALUE_TYPE *, enum machine_mode,
+ const FIXED_VALUE_TYPE *, int);
+
+/* Convert to a fixed-point mode from an integer. */
+extern bool fixed_convert_from_int (FIXED_VALUE_TYPE *, enum machine_mode,
+ double_int, int, int);
+
+/* Convert to a fixed-point mode from a real. */
+extern bool fixed_convert_from_real (FIXED_VALUE_TYPE *, enum machine_mode,
+ const REAL_VALUE_TYPE *, int);
+
+/* Convert to a real mode from a fixed-point. */
+extern void real_convert_from_fixed (REAL_VALUE_TYPE *, enum machine_mode,
+ const FIXED_VALUE_TYPE *);
+
+/* Compare two fixed-point objects for bitwise identity. */
+extern bool fixed_identical (const FIXED_VALUE_TYPE *, const FIXED_VALUE_TYPE *);
+
+/* Calculate a hash value. */
+extern unsigned int fixed_hash (const FIXED_VALUE_TYPE *);
+
+#define FIXED_VALUES_IDENTICAL(x, y) fixed_identical (&(x), &(y))
+
+/* Determine whether a fixed-point value X is negative. */
+#define FIXED_VALUE_NEGATIVE(x) fixed_isneg (&(x))
+
+/* Render F as a decimal floating point constant. */
+extern void fixed_to_decimal (char *str, const FIXED_VALUE_TYPE *, size_t);
+
+/* Binary or unary arithmetic on tree_code. */
+extern bool fixed_arithmetic (FIXED_VALUE_TYPE *, int, const FIXED_VALUE_TYPE *,
+ const FIXED_VALUE_TYPE *, int);
+
+/* Compare fixed-point values by tree_code. */
+extern bool fixed_compare (int, const FIXED_VALUE_TYPE *,
+ const FIXED_VALUE_TYPE *);
+
+/* Determine whether a fixed-point value X is negative. */
+extern bool fixed_isneg (const FIXED_VALUE_TYPE *);
+
+#endif /* GCC_FIXED_VALUE_H */
Index: Makefile.in
===================================================================
--- Makefile.in (revision 126007)
+++ Makefile.in (working copy)
@@ -999,6 +999,7 @@
expmed.o \
expr.o \
final.o \
+ fixed-value.o \
fold-const.o \
function.o \
fwprop.o \
@@ -2377,6 +2378,8 @@
toplev.h $(TM_P_H) $(REAL_H)
dfp.o : dfp.c dfp.h $(CONFIG_H) $(SYSTEM_H) coretypes.h $(TM_H) $(TREE_H) \
toplev.h $(TM_P_H) $(REAL_H) $(DECNUM_H)
+fixed-value.o: fixed-value.c $(CONFIG_H) $(SYSTEM_H) coretypes.h $(TM_H) \
+ $(TREE_H) fixed-value.h $(REAL_H) toplev.h
integrate.o : integrate.c $(CONFIG_H) $(SYSTEM_H) coretypes.h $(TM_H) \
$(RTL_H) $(TREE_H) $(FLAGS_H) debug.h $(INTEGRATE_H) insn-config.h \
$(EXPR_H) $(REAL_H) $(REGS_H) intl.h $(FUNCTION_H) output.h $(RECOG_H) \
@@ -3000,6 +3003,7 @@
$(srcdir)/coverage.c $(srcdir)/rtl.h \
$(srcdir)/optabs.h $(srcdir)/tree.h $(srcdir)/function.h $(srcdir)/libfuncs.h $(SYMTAB_H) \
$(srcdir)/real.h $(srcdir)/varray.h $(srcdir)/insn-addr.h $(srcdir)/hwint.h \
+ $(srcdir)/fixed-value.h \
$(srcdir)/ipa-reference.h $(srcdir)/output.h $(srcdir)/cfgloop.h \
$(srcdir)/cselib.h $(srcdir)/basic-block.h $(srcdir)/cgraph.h \
$(srcdir)/reload.h \
^ permalink raw reply [flat|nested] 4+ messages in thread
* Re: [PATCH] The 2nd patch from the fixed-point branch
2007-06-26 1:17 [PATCH] The 2nd patch from the fixed-point branch Fu, Chao-Ying
@ 2007-06-26 15:13 ` Richard Sandiford
2007-06-27 1:27 ` Fu, Chao-Ying
0 siblings, 1 reply; 4+ messages in thread
From: Richard Sandiford @ 2007-06-26 15:13 UTC (permalink / raw)
To: Fu, Chao-Ying
Cc: gcc-patches, Mark Mitchell, Thekkath, Radhika, Stephens, Nigel
"Fu, Chao-Ying" <fu@mips.com> writes:
> +static int check_real_for_fixed_mode (REAL_VALUE_TYPE *, unsigned int);
> +static int get_fixed_sign_bit (double_int, int);
> +static bool do_fixed_add (FIXED_VALUE_TYPE *, const FIXED_VALUE_TYPE *,
> + const FIXED_VALUE_TYPE *, int, int);
> +static bool do_fixed_multiply (FIXED_VALUE_TYPE *, const FIXED_VALUE_TYPE *,
> + const FIXED_VALUE_TYPE *, int);
> +static bool do_fixed_divide (FIXED_VALUE_TYPE *, const FIXED_VALUE_TYPE *,
> + const FIXED_VALUE_TYPE *, int);
> +static bool do_fixed_shift (FIXED_VALUE_TYPE *, const FIXED_VALUE_TYPE *,
> + const FIXED_VALUE_TYPE *, int, int);
> +static bool do_fixed_neg (FIXED_VALUE_TYPE *, const FIXED_VALUE_TYPE *, int);
> +static bool fixed_saturate1 (unsigned int, double_int, double_int *, int);
> +static bool fixed_saturate2 (unsigned int, double_int, double_int,
> + double_int *, int);
New style is not to have these forward declarations unless they are
really needed (due to mutual recursion, etc.)
> +/* Compare two fixed objects for bitwise identity. */
> +
> +bool
> +fixed_identical (const FIXED_VALUE_TYPE *a, const FIXED_VALUE_TYPE *b)
> +{
> + return a->mode == b->mode && a->data.high == b->data.high
> + && a->data.low == b->data.low;
Formatting nit, should be:
return (a->mode == b->mode
&& a->data.high == b->data.high
&& a->data.low == b->data.low);
> +
> +/* Check REAL_VALUE against the range of the fixed-point mode.
> + Return 0, if it is within the range.
> + 1, if it is less than the minimum.
> + 2, if it is greater than the maximum, but not equal to
> + the maximum + the epsilon.
> + 3, if it is equal to the maximum + the epsilon. */
Would the code be clearer if this were an enum? Or is this mirroring
some well-known library function?
> +int
> +check_real_for_fixed_mode (REAL_VALUE_TYPE *real_value, unsigned int mode)
Why is MODE not an enum machine_mode? (I'm sure there's a good reason,
I just wasn't sure what it was.)
> +{
> + REAL_VALUE_TYPE max_value, min_value, epsilon_value;
> + char max_string[20], min_string[20], epsilon_string[20];
> +
> + sprintf (max_string, "0x1.0p%d", GET_MODE_IBIT (mode));
> + sprintf (epsilon_string, "0x1.0p-%d", GET_MODE_FBIT (mode));
> + real_from_string (&max_value, max_string);
> + real_from_string (&epsilon_value, epsilon_string);
Couldn't you use real_2expN here instead? Other cases later in file.
> +void fixed_to_decimal (char *str, const FIXED_VALUE_TYPE *f_orig,
> + size_t buf_size)
Formatting; function name should start a new line.
> +/* If SATP, saturate A to the maximum or the minimum, and save to *F based on
> + the machine mode MODE.
> + This function assumes the width of double_int is greater than the width
> + of the fixed-point value at the fixed-point mode.
> + Return true, if !SATP and overflow. */
> +
> +static bool
> +fixed_saturate1 (unsigned int mode, double_int a, double_int *f, int satp)
Maybe satp should be a bool too? I don't quite understand what you
mean by "the fixed-point value at the fixed-point mode".
It might be clearer if you add "Do not modify *F otherwise." after the
first sentence. Same comments apply to fixed_saturate2.
> +/* Return the sign bit based on I_F_BITS. */
> +
> +inline int get_fixed_sign_bit (double_int a, int i_f_bits)
Must be "static inline". Formatting: function on new line.
> + if (((!subtract_p)
Redundant brackets.
> + && get_fixed_sign_bit (a->data, i_f_bits)
> + == get_fixed_sign_bit (b->data, i_f_bits)
Multiline conditions should be bracketed (so that emacs will indent
the "==" correctly).
> +/* Calculate F = A << B if LEFT_P. Otherwies, F = A >> B.
> + If SATP, saturate the result to the max or the min.
> + Return true, if !SATP and overflow. */
Mixture of naming styles: LEFT_P vs. SATP. FWIW, I think LEFT_P
is the usual convention in gcc. (Same thing with SUBTRACT_P earlier.)
> +/* Perform the binary or unary operation described by CODE.
> + For a unary operation, leave OP1 NULL.
> + Return true, if !SATP and overflow. */
> +
> +bool
> +fixed_arithmetic (FIXED_VALUE_TYPE *f, int icode, const FIXED_VALUE_TYPE *op0,
> + const FIXED_VALUE_TYPE *op1, int satp)
Am I right in thinking that OP0 and OP1 must have the same mode?
If so, it might be worth a comment and a gcc_assert().
The reason I ask is that...
> +/* Compare fixed-point values by tree_code. */
> +
> +bool
> +fixed_compare (int icode, const FIXED_VALUE_TYPE *op0,
> + const FIXED_VALUE_TYPE *op1)
> +{
> + enum tree_code code = icode;
> +
> + switch (code)
> + {
> + case NE_EXPR:
> + return op0->mode != op1->mode
> + || !double_int_equal_p (op0->data, op1->data);
...I was somewhat surprised that fixed_arithmentic appeared to need
operands of the same mode, but fixed_compare didn't. It also wasn't
immediately obvious that values of different modes were treated as
unequal, even if they have the same conceptual value. Might be worth
adding a comment.
> +/* Constant fixed-point values 0 and 1. */
> +extern FIXED_VALUE_TYPE fconst0[];
> +extern FIXED_VALUE_TYPE fconst1[];
What are the indexes?
> +/* Determine whether a fixed-point value X is negative. */
> +#define FIXED_VALUE_NEGATIVE(x) fixed_isneg (&(x))
Odd spacing (in context).
Richard
^ permalink raw reply [flat|nested] 4+ messages in thread
* RE: [PATCH] The 2nd patch from the fixed-point branch
2007-06-26 15:13 ` Richard Sandiford
@ 2007-06-27 1:27 ` Fu, Chao-Ying
2007-07-02 21:10 ` Fu, Chao-Ying
0 siblings, 1 reply; 4+ messages in thread
From: Fu, Chao-Ying @ 2007-06-27 1:27 UTC (permalink / raw)
To: Richard Sandiford
Cc: gcc-patches, Mark Mitchell, Thekkath, Radhika, Stephens, Nigel
[-- Attachment #1: Type: text/plain, Size: 7553 bytes --]
Richard Sandiford wrote:
> > +static int check_real_for_fixed_mode (REAL_VALUE_TYPE *,
> unsigned int);
> > +static int get_fixed_sign_bit (double_int, int);
> > +static bool do_fixed_add (FIXED_VALUE_TYPE *, const
> FIXED_VALUE_TYPE *,
> > + const FIXED_VALUE_TYPE *, int, int);
> > +static bool do_fixed_multiply (FIXED_VALUE_TYPE *, const
> FIXED_VALUE_TYPE *,
> > + const FIXED_VALUE_TYPE *, int);
> > +static bool do_fixed_divide (FIXED_VALUE_TYPE *, const
> FIXED_VALUE_TYPE *,
> > + const FIXED_VALUE_TYPE *, int);
> > +static bool do_fixed_shift (FIXED_VALUE_TYPE *, const
> FIXED_VALUE_TYPE *,
> > + const FIXED_VALUE_TYPE *, int, int);
> > +static bool do_fixed_neg (FIXED_VALUE_TYPE *, const
> FIXED_VALUE_TYPE *, int);
> > +static bool fixed_saturate1 (unsigned int, double_int,
> double_int *, int);
> > +static bool fixed_saturate2 (unsigned int, double_int, double_int,
> > + double_int *, int);
>
> New style is not to have these forward declarations unless they are
> really needed (due to mutual recursion, etc.)
Yes. These declarations are deleted.
>
> > +/* Compare two fixed objects for bitwise identity. */
> > +
> > +bool
> > +fixed_identical (const FIXED_VALUE_TYPE *a, const
> FIXED_VALUE_TYPE *b)
> > +{
> > + return a->mode == b->mode && a->data.high == b->data.high
> > + && a->data.low == b->data.low;
>
> Formatting nit, should be:
>
> return (a->mode == b->mode
> && a->data.high == b->data.high
> && a->data.low == b->data.low);
Yes.
>
> > +
> > +/* Check REAL_VALUE against the range of the fixed-point mode.
> > + Return 0, if it is within the range.
> > + 1, if it is less than the minimum.
> > + 2, if it is greater than the maximum, but not equal to
> > + the maximum + the epsilon.
> > + 3, if it is equal to the maximum + the epsilon. */
>
> Would the code be clearer if this were an enum? Or is this mirroring
> some well-known library function?
Yes. I created an enum for the returned values.
>
> > +int
> > +check_real_for_fixed_mode (REAL_VALUE_TYPE *real_value,
> unsigned int mode)
>
> Why is MODE not an enum machine_mode? (I'm sure there's a
> good reason,
> I just wasn't sure what it was.)
I just changed to enum machine_mode.
>
> > +{
> > + REAL_VALUE_TYPE max_value, min_value, epsilon_value;
> > + char max_string[20], min_string[20], epsilon_string[20];
> > +
> > + sprintf (max_string, "0x1.0p%d", GET_MODE_IBIT (mode));
> > + sprintf (epsilon_string, "0x1.0p-%d", GET_MODE_FBIT (mode));
> > + real_from_string (&max_value, max_string);
> > + real_from_string (&epsilon_value, epsilon_string);
>
> Couldn't you use real_2expN here instead? Other cases later in file.
Yes. I replaced them with real_2expN.
>
> > +void fixed_to_decimal (char *str, const FIXED_VALUE_TYPE *f_orig,
> > + size_t buf_size)
>
> Formatting; function name should start a new line.
Yes.
>
> > +/* If SATP, saturate A to the maximum or the minimum, and
> save to *F based on
> > + the machine mode MODE.
> > + This function assumes the width of double_int is
> greater than the width
> > + of the fixed-point value at the fixed-point mode.
> > + Return true, if !SATP and overflow. */
> > +
> > +static bool
> > +fixed_saturate1 (unsigned int mode, double_int a,
> double_int *f, int satp)
>
> Maybe satp should be a bool too? I don't quite understand what you
> mean by "the fixed-point value at the fixed-point mode".
SATP is a bool now.
I try to say that the width is the sum of a possible sign bit, possible ibits
and fbits => without any padding bits.
>
> It might be clearer if you add "Do not modify *F otherwise." after the
> first sentence. Same comments apply to fixed_saturate2.
Yes.
>
> > +/* Return the sign bit based on I_F_BITS. */
> > +
> > +inline int get_fixed_sign_bit (double_int a, int i_f_bits)
>
> Must be "static inline". Formatting: function on new line.
>
Yes.
> > + if (((!subtract_p)
>
> Redundant brackets.
Yes.
>
> > + && get_fixed_sign_bit (a->data, i_f_bits)
> > + == get_fixed_sign_bit (b->data, i_f_bits)
>
> Multiline conditions should be bracketed (so that emacs will indent
> the "==" correctly).
Yes.
>
> > +/* Calculate F = A << B if LEFT_P. Otherwies, F = A >> B.
> > + If SATP, saturate the result to the max or the min.
> > + Return true, if !SATP and overflow. */
>
> Mixture of naming styles: LEFT_P vs. SATP. FWIW, I think LEFT_P
> is the usual convention in gcc. (Same thing with SUBTRACT_P earlier.)
I changed SATP to SAT_P.
>
> > +/* Perform the binary or unary operation described by CODE.
> > + For a unary operation, leave OP1 NULL.
> > + Return true, if !SATP and overflow. */
> > +
> > +bool
> > +fixed_arithmetic (FIXED_VALUE_TYPE *f, int icode, const
> FIXED_VALUE_TYPE *op0,
> > + const FIXED_VALUE_TYPE *op1, int satp)
>
>
> Am I right in thinking that OP0 and OP1 must have the same mode?
> If so, it might be worth a comment and a gcc_assert().
Yes.
> The reason I ask is that...
>
> > +/* Compare fixed-point values by tree_code. */
> > +
> > +bool
> > +fixed_compare (int icode, const FIXED_VALUE_TYPE *op0,
> > + const FIXED_VALUE_TYPE *op1)
> > +{
> > + enum tree_code code = icode;
> > +
> > + switch (code)
> > + {
> > + case NE_EXPR:
> > + return op0->mode != op1->mode
> > + || !double_int_equal_p (op0->data, op1->data);
>
> ...I was somewhat surprised that fixed_arithmentic appeared to need
> operands of the same mode, but fixed_compare didn't. It also wasn't
> immediately obvious that values of different modes were treated as
> unequal, even if they have the same conceptual value. Might be worth
> adding a comment.
"fixed_compare" must have the same mode as well. I changed the comment
and the code to have gcc_assert.
>
> > +/* Constant fixed-point values 0 and 1. */
> > +extern FIXED_VALUE_TYPE fconst0[];
> > +extern FIXED_VALUE_TYPE fconst1[];
>
> What are the indexes?
Ok. I created two defines and two new macros to access them.
#define MAX_FCONST0 18 /* For storing 18 fixed-point zeros per
fract, ufract, accum, and uaccum modes . */
#define MAX_FCONST1 8 /* For storing 8 fixed-point ones per accum
and uaccum modes. */
/* Constant fixed-point values 0 and 1. */
extern FIXED_VALUE_TYPE fconst0[MAX_FCONST0];
extern FIXED_VALUE_TYPE fconst1[MAX_FCONST1];
/* Macros to access fconst0 and fconst1 via machine modes. */
#define FCONST0(mode) fconst0[mode - QQmode]
#define FCONST1(mode) fconst1[mode - HAmode]
>
> > +/* Determine whether a fixed-point value X is negative. */
> > +#define FIXED_VALUE_NEGATIVE(x) fixed_isneg (&(x))
>
> Odd spacing (in context).
Yes.
The new patch is attached. Thanks a lot!
Regards,
Chao-ying
2007-06-25 Chao-ying Fu <fu@mips.com>
* gengtype.c (main): Handle FIXED_VALUE_TYPE type as scalar typedef.
* double-int.c (double_int_scmp): Use casts of unsigned HOST_WIDE_INT
to compare a.low and b.low.
* fixed-value.c: New file.
* fixed-value.h: New file.
* Makefile.in (OBJS-common): Add fixed-value.o.
(fixed-value.o): New rule.
(GTFILES): Add fixed-value.h.
[-- Attachment #2: fixed-value.diff --]
[-- Type: application/octet-stream, Size: 39506 bytes --]
Index: gengtype.c
===================================================================
--- gengtype.c (revision 126045)
+++ gengtype.c (working copy)
@@ -3535,6 +3535,7 @@
pos.line = __LINE__ + 1;
do_scalar_typedef ("CUMULATIVE_ARGS", &pos); pos.line++;
do_scalar_typedef ("REAL_VALUE_TYPE", &pos); pos.line++;
+ do_scalar_typedef ("FIXED_VALUE_TYPE", &pos); pos.line++;
do_scalar_typedef ("double_int", &pos); pos.line++;
do_scalar_typedef ("uint8", &pos); pos.line++;
do_scalar_typedef ("jword", &pos); pos.line++;
Index: double-int.c
===================================================================
--- double-int.c (revision 126045)
+++ double-int.c (working copy)
@@ -363,9 +363,9 @@
return -1;
if (a.high > b.high)
return 1;
- if ((HOST_WIDE_INT) a.low < (HOST_WIDE_INT) b.low)
+ if ((unsigned HOST_WIDE_INT) a.low < (unsigned HOST_WIDE_INT) b.low)
return -1;
- if ((HOST_WIDE_INT) a.low > (HOST_WIDE_INT) b.low)
+ if ((unsigned HOST_WIDE_INT) a.low > (unsigned HOST_WIDE_INT) b.low)
return 1;
return 0;
Index: fixed-value.c
===================================================================
--- fixed-value.c (revision 0)
+++ fixed-value.c (revision 0)
@@ -0,0 +1,1138 @@
+/* Fixed-point arithmetic support.
+ Copyright (C) 2006, 2007 Free Software Foundation, Inc.
+
+This file is part of GCC.
+
+GCC is free software; you can redistribute it and/or modify it
+under the terms of the GNU General Public License as published by the
+Free Software Foundation; either version 2, or (at your option) any
+later version.
+
+GCC is distributed in the hope that it will be useful, but WITHOUT
+ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
+for more details.
+
+You should have received a copy of the GNU General Public License
+along with GCC; see the file COPYING. If not, write to the Free
+Software Foundation, 51 Franklin Street, Fifth Floor, Boston, MA
+02110-1301, USA. */
+
+#include "config.h"
+#include "system.h"
+#include "coretypes.h"
+#include "tm.h"
+#include "tree.h"
+#include "toplev.h"
+#include "fixed-value.h"
+
+/* Compare two fixed objects for bitwise identity. */
+
+bool
+fixed_identical (const FIXED_VALUE_TYPE *a, const FIXED_VALUE_TYPE *b)
+{
+ return (a->mode == b->mode
+ && a->data.high == b->data.high
+ && a->data.low == b->data.low);
+}
+
+/* Calculate a hash value. */
+
+unsigned int
+fixed_hash (const FIXED_VALUE_TYPE *f)
+{
+ return (unsigned int) (f->data.low ^ f->data.high);
+}
+
+/* Define the enum code for the range of the fixed-point value. */
+enum fixed_value_range_code {
+ FIXED_OK, /* The value is within the range. */
+ FIXED_UNDERFLOW, /* The value is less than the minimum. */
+ FIXED_GT_MAX_EPS, /* The value is greater than the maximum, but not equal
+ to the maximum plus the epsilon. */
+ FIXED_MAX_EPS /* The value equals the maximum plus the epsilon. */
+};
+
+/* Check REAL_VALUE against the range of the fixed-point mode.
+ Return FIXED_OK, if it is within the range.
+ FIXED_UNDERFLOW, if it is less than the minimum.
+ FIXED_GT_MAX_EPS, if it is greater than the maximum, but not equal to
+ the maximum plus the epsilon.
+ FIXED_MAX_EPS, if it is equal to the maximum plus the epsilon. */
+
+static enum fixed_value_range_code
+check_real_for_fixed_mode (REAL_VALUE_TYPE *real_value, enum machine_mode mode)
+{
+ REAL_VALUE_TYPE max_value, min_value, epsilon_value;
+
+ real_2expN (&max_value, GET_MODE_IBITS (mode));
+ real_2expN (&epsilon_value, -GET_MODE_FBITS (mode));
+
+ if (SIGNED_FIXED_POINT_MODE_P (mode))
+ min_value = REAL_VALUE_NEGATE (max_value);
+ else
+ real_from_string (&min_value, "0.0");
+
+ if (real_compare (LT_EXPR, real_value, &min_value))
+ return FIXED_UNDERFLOW;
+ if (real_compare (EQ_EXPR, real_value, &max_value))
+ return FIXED_MAX_EPS;
+ real_arithmetic (&max_value, MINUS_EXPR, &max_value, &epsilon_value);
+ if (real_compare (GT_EXPR, real_value, &max_value))
+ return FIXED_GT_MAX_EPS;
+ return FIXED_OK;
+}
+
+/* Initialize from a decimal or hexadecimal string. */
+
+void
+fixed_from_string (FIXED_VALUE_TYPE *f, const char *str, enum machine_mode mode)
+{
+ REAL_VALUE_TYPE real_value, fixed_value, base_value;
+ unsigned int fbit;
+ enum fixed_value_range_code temp;
+
+ f->mode = mode;
+ fbit = GET_MODE_FBITS (mode);
+
+ real_from_string (&real_value, str);
+ temp = check_real_for_fixed_mode (&real_value, f->mode);
+ /* We don't want to warn the case when the _Fract value is 1.0. */
+ if (temp == FIXED_UNDERFLOW
+ || temp == FIXED_GT_MAX_EPS
+ || (temp == FIXED_MAX_EPS && ALL_ACCUM_MODE_P (f->mode)))
+ warning (OPT_Woverflow,
+ "large fixed-point constant implicitly truncated to fixed-point type");
+ real_2expN (&base_value, fbit);
+ real_arithmetic (&fixed_value, MULT_EXPR, &real_value, &base_value);
+ real_to_integer2 ((HOST_WIDE_INT *)&f->data.low, &f->data.high,
+ &fixed_value);
+
+ if (temp == FIXED_MAX_EPS && ALL_FRACT_MODE_P (f->mode))
+ {
+ /* From the spec, we need to evaluate 1 to the maximal value. */
+ f->data.low = -1;
+ f->data.high = -1;
+ f->data = double_int_ext (f->data,
+ GET_MODE_FBITS (f->mode)
+ + GET_MODE_IBITS (f->mode), 1);
+ }
+ else
+ f->data = double_int_ext (f->data,
+ SIGNED_FIXED_POINT_MODE_P (f->mode)
+ + GET_MODE_FBITS (f->mode)
+ + GET_MODE_IBITS (f->mode),
+ UNSIGNED_FIXED_POINT_MODE_P (f->mode));
+}
+
+/* Render F as a decimal floating point constant. */
+
+void
+fixed_to_decimal (char *str, const FIXED_VALUE_TYPE *f_orig,
+ size_t buf_size)
+{
+ REAL_VALUE_TYPE real_value, base_value, fixed_value;
+
+ real_2expN (&base_value, GET_MODE_FBITS (f_orig->mode));
+ real_from_integer (&real_value, VOIDmode, f_orig->data.low, f_orig->data.high,
+ UNSIGNED_FIXED_POINT_MODE_P (f_orig->mode));
+ real_arithmetic (&fixed_value, RDIV_EXPR, &real_value, &base_value);
+ real_to_decimal (str, &fixed_value, buf_size, 0, 1);
+}
+
+/* If SAT_P, saturate A to the maximum or the minimum, and save to *F based on
+ the machine mode MODE.
+ Do not modify *F otherwise.
+ This function assumes the width of double_int is greater than the width
+ of the fixed-point value (the sum of a possible sign bit, possible ibits,
+ and fbits).
+ Return true, if !SAT_P and overflow. */
+
+static bool
+fixed_saturate1 (enum machine_mode mode, double_int a, double_int *f,
+ bool sat_p)
+{
+ bool overflow = false;
+ bool unsigned_p = UNSIGNED_FIXED_POINT_MODE_P (mode);
+ int i_f_bits = GET_MODE_IBITS (mode) + GET_MODE_FBITS (mode);
+
+ if (unsigned_p) /* Unsigned type. */
+ {
+ double_int max;
+ max.low = -1;
+ max.high = -1;
+ max = double_int_ext (max, i_f_bits, 1);
+ if (double_int_cmp (a, max, 1) == 1)
+ {
+ if (sat_p)
+ *f = max;
+ else
+ overflow = true;
+ }
+ }
+ else /* Signed type. */
+ {
+ double_int max, min;
+ max.high = -1;
+ max.low = -1;
+ max = double_int_ext (max, i_f_bits, 1);
+ min.high = 0;
+ min.low = 1;
+ lshift_double (min.low, min.high, i_f_bits,
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &min.low, &min.high, 1);
+ min = double_int_ext (min, 1 + i_f_bits, 0);
+ if (double_int_cmp (a, max, 0) == 1)
+ {
+ if (sat_p)
+ *f = max;
+ else
+ overflow = true;
+ }
+ else if (double_int_cmp (a, min, 0) == -1)
+ {
+ if (sat_p)
+ *f = min;
+ else
+ overflow = true;
+ }
+ }
+ return overflow;
+}
+
+/* If SAT_P, saturate {A_HIGH, A_LOW} to the maximum or the minimum, and
+ save to *F based on the machine mode MODE.
+ Do not modify *F otherwise.
+ This function assumes the width of two double_int is greater than the width
+ of the fixed-point value (the sum of a possible sign bit, possible ibits,
+ and fbits).
+ Return true, if !SAT_P and overflow. */
+
+static bool
+fixed_saturate2 (enum machine_mode mode, double_int a_high, double_int a_low,
+ double_int *f, bool sat_p)
+{
+ bool overflow = false;
+ bool unsigned_p = UNSIGNED_FIXED_POINT_MODE_P (mode);
+ int i_f_bits = GET_MODE_IBITS (mode) + GET_MODE_FBITS (mode);
+
+ if (unsigned_p) /* Unsigned type. */
+ {
+ double_int max_r, max_s;
+ max_r.high = 0;
+ max_r.low = 0;
+ max_s.high = -1;
+ max_s.low = -1;
+ max_s = double_int_ext (max_s, i_f_bits, 1);
+ if (double_int_cmp (a_high, max_r, 1) == 1
+ || (double_int_equal_p (a_high, max_r) &&
+ double_int_cmp (a_low, max_s, 1) == 1))
+ {
+ if (sat_p)
+ *f = max_s;
+ else
+ overflow = true;
+ }
+ }
+ else /* Signed type. */
+ {
+ double_int max_r, max_s, min_r, min_s;
+ max_r.high = 0;
+ max_r.low = 0;
+ max_s.high = -1;
+ max_s.low = -1;
+ max_s = double_int_ext (max_s, i_f_bits, 1);
+ min_r.high = -1;
+ min_r.low = -1;
+ min_s.high = 0;
+ min_s.low = 1;
+ lshift_double (min_s.low, min_s.high, i_f_bits,
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &min_s.low, &min_s.high, 1);
+ min_s = double_int_ext (min_s, 1 + i_f_bits, 0);
+ if (double_int_cmp (a_high, max_r, 0) == 1
+ || (double_int_equal_p (a_high, max_r) &&
+ double_int_cmp (a_low, max_s, 1) == 1))
+ {
+ if (sat_p)
+ *f = max_s;
+ else
+ overflow = true;
+ }
+ else if (double_int_cmp (a_high, min_r, 0) == -1
+ || (double_int_equal_p (a_high, min_r) &&
+ double_int_cmp (a_low, min_s, 1) == -1))
+ {
+ if (sat_p)
+ *f = min_s;
+ else
+ overflow = true;
+ }
+ }
+ return overflow;
+}
+
+/* Return the sign bit based on I_F_BITS. */
+
+static inline int
+get_fixed_sign_bit (double_int a, int i_f_bits)
+{
+ if (i_f_bits < HOST_BITS_PER_WIDE_INT)
+ return (a.low >> i_f_bits) & 1;
+ else
+ return (a.high >> (i_f_bits - HOST_BITS_PER_WIDE_INT)) & 1;
+}
+
+/* Calculate F = A + (SUBTRACT_P ? -B : B).
+ If SAT_P, saturate the result to the max or the min.
+ Return true, if !SAT_P and overflow. */
+
+static bool
+do_fixed_add (FIXED_VALUE_TYPE *f, const FIXED_VALUE_TYPE *a,
+ const FIXED_VALUE_TYPE *b, bool subtract_p, bool sat_p)
+{
+ bool overflow = false;
+ double_int temp = subtract_p ? double_int_neg (b->data) : b->data;
+ bool unsigned_p = UNSIGNED_FIXED_POINT_MODE_P (a->mode);
+ int i_f_bits = GET_MODE_IBITS (a->mode) + GET_MODE_FBITS (a->mode);
+ f->mode = a->mode;
+ f->data = double_int_add (a->data, temp);
+ if (unsigned_p) /* Unsigned type. */
+ {
+ if (subtract_p) /* Unsigned subtraction. */
+ {
+ if (double_int_cmp (a->data, b->data, 1) == -1)
+ {
+ if (sat_p)
+ {
+ f->data.high = 0;
+ f->data.low = 0;
+ }
+ else
+ overflow = true;
+ }
+ }
+ else /* Unsigned addition. */
+ {
+ f->data = double_int_ext (f->data, i_f_bits, 1);
+ if (double_int_cmp (f->data, a->data, 1) == -1
+ || double_int_cmp (f->data, b->data, 1) == -1)
+ {
+ if (sat_p)
+ {
+ f->data.high = -1;
+ f->data.low = -1;
+ }
+ else
+ overflow = true;
+ }
+ }
+ }
+ else /* Signed type. */
+ {
+ if ((!subtract_p
+ && (get_fixed_sign_bit (a->data, i_f_bits)
+ == get_fixed_sign_bit (b->data, i_f_bits))
+ && (get_fixed_sign_bit (a->data, i_f_bits)
+ != get_fixed_sign_bit (f->data, i_f_bits)))
+ || (subtract_p
+ && (get_fixed_sign_bit (a->data, i_f_bits)
+ != get_fixed_sign_bit (b->data, i_f_bits))
+ && (get_fixed_sign_bit (a->data, i_f_bits)
+ != get_fixed_sign_bit (f->data, i_f_bits))))
+ {
+ if (sat_p)
+ {
+ f->data.low = 1;
+ f->data.high = 0;
+ lshift_double (f->data.low, f->data.high, i_f_bits,
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &f->data.low, &f->data.high, 1);
+ if (get_fixed_sign_bit (a->data, i_f_bits) == 0)
+ {
+ double_int one;
+ one.low = 1;
+ one.high = 0;
+ f->data = double_int_add (f->data, double_int_neg (one));
+ }
+ }
+ else
+ overflow = true;
+ }
+ }
+ f->data = double_int_ext (f->data, (!unsigned_p) + i_f_bits, unsigned_p);
+ return overflow;
+}
+
+/* Calculate F = A * B.
+ If SAT_P, saturate the result to the max or the min.
+ Return true, if !SAT_P and overflow. */
+
+static bool
+do_fixed_multiply (FIXED_VALUE_TYPE *f, const FIXED_VALUE_TYPE *a,
+ const FIXED_VALUE_TYPE *b, bool sat_p)
+{
+ bool overflow = false;
+ bool unsigned_p = UNSIGNED_FIXED_POINT_MODE_P (a->mode);
+ int i_f_bits = GET_MODE_IBITS (a->mode) + GET_MODE_FBITS (a->mode);
+ f->mode = a->mode;
+ if (GET_MODE_PRECISION (f->mode) <= HOST_BITS_PER_WIDE_INT)
+ {
+ f->data = double_int_mul (a->data, b->data);
+ lshift_double (f->data.low, f->data.high,
+ (-GET_MODE_FBITS (f->mode)),
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &f->data.low, &f->data.high, !unsigned_p);
+ overflow = fixed_saturate1 (f->mode, f->data, &f->data, sat_p);
+ }
+ else
+ {
+ /* The result of multiplication expands to two double_int. */
+ double_int a_high, a_low, b_high, b_low;
+ double_int high_high, high_low, low_high, low_low;
+ double_int r, s, temp1, temp2;
+ int carry = 0;
+
+ /* Decompose a and b to four double_int. */
+ a_high.low = a->data.high;
+ a_high.high = 0;
+ a_low.low = a->data.low;
+ a_low.high = 0;
+ b_high.low = b->data.high;
+ b_high.high = 0;
+ b_low.low = b->data.low;
+ b_low.high = 0;
+
+ /* Perform four multiplications. */
+ low_low = double_int_mul (a_low, b_low);
+ low_high = double_int_mul (a_low, b_high);
+ high_low = double_int_mul (a_high, b_low);
+ high_high = double_int_mul (a_high, b_high);
+
+ /* Accumulate four results to {r, s}. */
+ temp1.high = high_low.low;
+ temp1.low = 0;
+ s = double_int_add (low_low, temp1);
+ if (double_int_cmp (s, low_low, 1) == -1
+ || double_int_cmp (s, temp1, 1) == -1)
+ carry ++; /* Carry */
+ temp1.high = s.high;
+ temp1.low = s.low;
+ temp2.high = low_high.low;
+ temp2.low = 0;
+ s = double_int_add (temp1, temp2);
+ if (double_int_cmp (s, temp1, 1) == -1
+ || double_int_cmp (s, temp2, 1) == -1)
+ carry ++; /* Carry */
+
+ temp1.low = high_low.high;
+ temp1.high = 0;
+ r = double_int_add (high_high, temp1);
+ temp1.low = low_high.high;
+ temp1.high = 0;
+ r = double_int_add (r, temp1);
+ temp1.low = carry;
+ temp1.high = 0;
+ r = double_int_add (r, temp1);
+
+ /* We need to add neg(b) to r, if a < 0. */
+ if (!unsigned_p && a->data.high < 0)
+ r = double_int_add (r, double_int_neg (b->data));
+ /* We need to add neg(a) to r, if b < 0. */
+ if (!unsigned_p && b->data.high < 0)
+ r = double_int_add (r, double_int_neg (a->data));
+
+ /* Shift right the result by FBIT. */
+ if (GET_MODE_FBITS (f->mode) == 2 * HOST_BITS_PER_WIDE_INT)
+ {
+ s.low = r.low;
+ s.high = r.high;
+ if (unsigned_p)
+ {
+ r.low = 0;
+ r.high = 0;
+ }
+ else
+ {
+ r.low = -1;
+ r.high = -1;
+ }
+ f->data.low = s.low;
+ f->data.high = s.high;
+ }
+ else
+ {
+ lshift_double (s.low, s.high,
+ (-GET_MODE_FBITS (f->mode)),
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &s.low, &s.high, 0);
+ lshift_double (r.low, r.high,
+ (2 * HOST_BITS_PER_WIDE_INT
+ - GET_MODE_FBITS (f->mode)),
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &f->data.low, &f->data.high, 0);
+ f->data.low = f->data.low | s.low;
+ f->data.high = f->data.high | s.high;
+ s.low = f->data.low;
+ s.high = f->data.high;
+ lshift_double (r.low, r.high,
+ (-GET_MODE_FBITS (f->mode)),
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &r.low, &r.high, !unsigned_p);
+ }
+
+ overflow = fixed_saturate2 (f->mode, r, s, &f->data, sat_p);
+ }
+
+ f->data = double_int_ext (f->data, (!unsigned_p) + i_f_bits, unsigned_p);
+ return false;
+}
+
+/* Calculate F = A / B.
+ If SAT_P, saturate the result to the max or the min.
+ Return true, if !SAT_P and overflow. */
+
+static bool
+do_fixed_divide (FIXED_VALUE_TYPE *f, const FIXED_VALUE_TYPE *a,
+ const FIXED_VALUE_TYPE *b, bool sat_p)
+{
+ bool overflow = false;
+ bool unsigned_p = UNSIGNED_FIXED_POINT_MODE_P (a->mode);
+ int i_f_bits = GET_MODE_IBITS (a->mode) + GET_MODE_FBITS (a->mode);
+ f->mode = a->mode;
+ if (GET_MODE_PRECISION (f->mode) <= HOST_BITS_PER_WIDE_INT)
+ {
+ lshift_double (a->data.low, a->data.high,
+ GET_MODE_FBITS (f->mode),
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &f->data.low, &f->data.high, !unsigned_p);
+ f->data = double_int_div (f->data, b->data, unsigned_p, TRUNC_DIV_EXPR);
+ overflow = fixed_saturate1 (f->mode, f->data, &f->data, sat_p);
+ }
+ else
+ {
+ double_int pos_a, pos_b, r, s;
+ double_int quo_r, quo_s, mod, temp;
+ int num_of_neg = 0;
+ int i;
+
+ /* If a < 0, negate a. */
+ if (!unsigned_p && a->data.high < 0)
+ {
+ pos_a = double_int_neg (a->data);
+ num_of_neg ++;
+ }
+ else
+ pos_a = a->data;
+
+ /* If b < 0, negate b. */
+ if (!unsigned_p && b->data.high < 0)
+ {
+ pos_b = double_int_neg (b->data);
+ num_of_neg ++;
+ }
+ else
+ pos_b = b->data;
+
+ /* Left shift pos_a to {r, s} by FBIT. */
+ if (GET_MODE_FBITS (f->mode) == 2 * HOST_BITS_PER_WIDE_INT)
+ {
+ r = pos_a;
+ s.high = 0;
+ s.low = 0;
+ }
+ else
+ {
+ lshift_double (pos_a.low, pos_a.high,
+ GET_MODE_FBITS (f->mode),
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &s.low, &s.high, 0);
+ lshift_double (pos_a.low, pos_a.high,
+ - (2 * HOST_BITS_PER_WIDE_INT
+ - GET_MODE_FBITS (f->mode)),
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &r.low, &r.high, 0);
+ }
+
+ /* Divide r by pos_b to quo_r. The remanider is in mod. */
+ div_and_round_double (TRUNC_DIV_EXPR, 1, r.low, r.high, pos_b.low,
+ pos_b.high, &quo_r.low, &quo_r.high, &mod.low,
+ &mod.high);
+
+ quo_s.high = 0;
+ quo_s.low = 0;
+
+ for (i = 0; i < 2 * HOST_BITS_PER_WIDE_INT; i++)
+ {
+ /* Record the leftmost bit of mod. */
+ int leftmost_mod = (mod.high < 0);
+
+ /* Shift left mod by 1 bit. */
+ lshift_double (mod.low, mod.high, 1, 2 * HOST_BITS_PER_WIDE_INT,
+ &mod.low, &mod.high, 0);
+
+ /* Test the leftmost bit of s to add to mod. */
+ if (s.high < 0)
+ mod.low += 1;
+
+ /* Shift left quo_s by 1 bit. */
+ lshift_double (quo_s.low, quo_s.high, 1, 2 * HOST_BITS_PER_WIDE_INT,
+ &quo_s.low, &quo_s.high, 0);
+
+ /* Try to calculate (mod - pos_b). */
+ temp = double_int_add (mod, double_int_neg (pos_b));
+
+ if (leftmost_mod == 1 || double_int_cmp (mod, pos_b, 1) != -1)
+ {
+ quo_s.low += 1;
+ mod = temp;
+ }
+
+ /* Shift left s by 1 bit. */
+ lshift_double (s.low, s.high, 1, 2 * HOST_BITS_PER_WIDE_INT,
+ &s.low, &s.high, 0);
+
+ }
+
+ if (num_of_neg == 1)
+ {
+ quo_s = double_int_neg (quo_s);
+ if (quo_s.high == 0 && quo_s.low == 0)
+ quo_r = double_int_neg (quo_r);
+ else
+ {
+ quo_r.low = ~quo_r.low;
+ quo_r.high = ~quo_r.high;
+ }
+ }
+
+ f->data = quo_s;
+ overflow = fixed_saturate2 (f->mode, quo_r, quo_s, &f->data, sat_p);
+ }
+
+ f->data = double_int_ext (f->data, (!unsigned_p) + i_f_bits, unsigned_p);
+ return overflow;
+}
+
+/* Calculate F = A << B if LEFT_P. Otherwies, F = A >> B.
+ If SAT_P, saturate the result to the max or the min.
+ Return true, if !SAT_P and overflow. */
+
+static bool
+do_fixed_shift (FIXED_VALUE_TYPE *f, const FIXED_VALUE_TYPE *a,
+ const FIXED_VALUE_TYPE *b, bool left_p, bool sat_p)
+{
+ bool overflow = false;
+ bool unsigned_p = UNSIGNED_FIXED_POINT_MODE_P (a->mode);
+ int i_f_bits = GET_MODE_IBITS (a->mode) + GET_MODE_FBITS (a->mode);
+ f->mode = a->mode;
+
+ if (b->data.low == 0)
+ {
+ f->data = a->data;
+ return overflow;
+ }
+
+ if (GET_MODE_PRECISION (f->mode) <= HOST_BITS_PER_WIDE_INT || (!left_p))
+ {
+ lshift_double (a->data.low, a->data.high,
+ left_p ? b->data.low : (-b->data.low),
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &f->data.low, &f->data.high, !unsigned_p);
+ if (left_p) /* Only left shift saturates. */
+ overflow = fixed_saturate1 (f->mode, f->data, &f->data, sat_p);
+ }
+ else /* We need two double_int to store the left-shift result. */
+ {
+ double_int temp_high, temp_low;
+ if (b->data.low == 2 * HOST_BITS_PER_WIDE_INT)
+ {
+ temp_high = a->data;
+ temp_low.high = 0;
+ temp_low.low = 0;
+ }
+ else
+ {
+ lshift_double (a->data.low, a->data.high,
+ b->data.low,
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &temp_low.low, &temp_low.high, !unsigned_p);
+ /* Logical shift right to temp_high. */
+ lshift_double (a->data.low, a->data.high,
+ b->data.low - 2 * HOST_BITS_PER_WIDE_INT,
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &temp_high.low, &temp_high.high, 0);
+ }
+ if (!unsigned_p && a->data.high < 0) /* Signed-extend temp_high. */
+ temp_high = double_int_ext (temp_high, b->data.low, unsigned_p);
+ f->data = temp_low;
+ overflow = fixed_saturate2 (f->mode, temp_high, temp_low, &f->data,
+ sat_p);
+ }
+ f->data = double_int_ext (f->data, (!unsigned_p) + i_f_bits, unsigned_p);
+ return overflow;
+}
+
+/* Calculate F = -A.
+ If SAT_P, saturate the result to the max or the min.
+ Return true, if !SAT_P and overflow. */
+
+static bool
+do_fixed_neg (FIXED_VALUE_TYPE *f, const FIXED_VALUE_TYPE *a, bool sat_p)
+{
+ bool overflow = false;
+ bool unsigned_p = UNSIGNED_FIXED_POINT_MODE_P (a->mode);
+ int i_f_bits = GET_MODE_IBITS (a->mode) + GET_MODE_FBITS (a->mode);
+ f->mode = a->mode;
+ f->data = double_int_neg (a->data);
+ f->data = double_int_ext (f->data, (!unsigned_p) + i_f_bits, unsigned_p);
+
+ if (unsigned_p) /* Unsigned type. */
+ {
+ if (f->data.low != 0 || f->data.high != 0)
+ {
+ if (sat_p)
+ {
+ f->data.low = 0;
+ f->data.high = 0;
+ }
+ else
+ overflow = true;
+ }
+ }
+ else /* Signed type. */
+ {
+ if (!(f->data.high == 0 && f->data.low == 0)
+ && f->data.high == a->data.high && f->data.low == a->data.low )
+ {
+ if (sat_p)
+ {
+ /* Saturate to the maximum by subtracting f->data by one. */
+ f->data.low = -1;
+ f->data.high = -1;
+ f->data = double_int_ext (f->data, i_f_bits, 1);
+ }
+ else
+ overflow = true;
+ }
+ }
+ return overflow;
+}
+
+/* Perform the binary or unary operation described by CODE.
+ Note that OP0 and OP1 must have the same mode for binary operators.
+ For a unary operation, leave OP1 NULL.
+ Return true, if !SAT_P and overflow. */
+
+bool
+fixed_arithmetic (FIXED_VALUE_TYPE *f, int icode, const FIXED_VALUE_TYPE *op0,
+ const FIXED_VALUE_TYPE *op1, bool sat_p)
+{
+ switch (icode)
+ {
+ case NEGATE_EXPR:
+ return do_fixed_neg (f, op0, sat_p);
+ break;
+
+ case PLUS_EXPR:
+ gcc_assert (op0->mode == op1->mode);
+ return do_fixed_add (f, op0, op1, false, sat_p);
+ break;
+
+ case MINUS_EXPR:
+ gcc_assert (op0->mode == op1->mode);
+ return do_fixed_add (f, op0, op1, true, sat_p);
+ break;
+
+ case MULT_EXPR:
+ gcc_assert (op0->mode == op1->mode);
+ return do_fixed_multiply (f, op0, op1, sat_p);
+ break;
+
+ case TRUNC_DIV_EXPR:
+ gcc_assert (op0->mode == op1->mode);
+ return do_fixed_divide (f, op0, op1, sat_p);
+ break;
+
+ case LSHIFT_EXPR:
+ return do_fixed_shift (f, op0, op1, true, sat_p);
+ break;
+
+ case RSHIFT_EXPR:
+ return do_fixed_shift (f, op0, op1, false, sat_p);
+ break;
+
+ default:
+ gcc_unreachable ();
+ }
+ return false;
+}
+
+/* Compare fixed-point values by tree_code.
+ Note that OP0 and OP1 must have the same mode. */
+
+bool
+fixed_compare (int icode, const FIXED_VALUE_TYPE *op0,
+ const FIXED_VALUE_TYPE *op1)
+{
+ enum tree_code code = icode;
+ gcc_assert (op0->mode == op1->mode);
+
+ switch (code)
+ {
+ case NE_EXPR:
+ return !double_int_equal_p (op0->data, op1->data);
+
+ case EQ_EXPR:
+ return double_int_equal_p (op0->data, op1->data);
+
+ case LT_EXPR:
+ return double_int_cmp (op0->data, op1->data,
+ UNSIGNED_FIXED_POINT_MODE_P (op0->mode)) == -1;
+
+ case LE_EXPR:
+ return double_int_cmp (op0->data, op1->data,
+ UNSIGNED_FIXED_POINT_MODE_P (op0->mode)) != 1;
+
+ case GT_EXPR:
+ return double_int_cmp (op0->data, op1->data,
+ UNSIGNED_FIXED_POINT_MODE_P (op0->mode)) == 1;
+
+ case GE_EXPR:
+ return double_int_cmp (op0->data, op1->data,
+ UNSIGNED_FIXED_POINT_MODE_P (op0->mode)) != -1;
+
+ default:
+ gcc_unreachable ();
+ }
+}
+
+/* Extend or truncate to a new mode.
+ If SAT_P, saturate the result to the max or the min.
+ Return true, if !SAT_P and overflow. */
+
+bool
+fixed_convert (FIXED_VALUE_TYPE *f, enum machine_mode mode,
+ const FIXED_VALUE_TYPE *a, bool sat_p)
+{
+ bool overflow = false;
+ if (mode == a->mode)
+ {
+ *f = *a;
+ return overflow;
+ }
+
+ if (GET_MODE_FBITS (mode) > GET_MODE_FBITS (a->mode))
+ {
+ /* Left shift a to temp_high, temp_low based on a->mode. */
+ double_int temp_high, temp_low;
+ int amount = GET_MODE_FBITS (mode) - GET_MODE_FBITS (a->mode);
+ lshift_double (a->data.low, a->data.high,
+ amount,
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &temp_low.low, &temp_low.high,
+ SIGNED_FIXED_POINT_MODE_P (a->mode));
+ /* Logical shift right to temp_high. */
+ lshift_double (a->data.low, a->data.high,
+ amount - 2 * HOST_BITS_PER_WIDE_INT,
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &temp_high.low, &temp_high.high, 0);
+ if (SIGNED_FIXED_POINT_MODE_P (a->mode)
+ && a->data.high < 0) /* Signed-extend temp_high. */
+ temp_high = double_int_ext (temp_high, amount, 0);
+ f->mode = mode;
+ f->data = temp_low;
+ if (SIGNED_FIXED_POINT_MODE_P (a->mode) ==
+ SIGNED_FIXED_POINT_MODE_P (f->mode))
+ overflow = fixed_saturate2 (f->mode, temp_high, temp_low, &f->data,
+ sat_p);
+ else
+ {
+ /* Take care of the cases when converting between signed and
+ unsigned. */
+ if (SIGNED_FIXED_POINT_MODE_P (a->mode))
+ {
+ /* Signed -> Unsigned. */
+ if (a->data.high < 0)
+ {
+ if (sat_p)
+ {
+ f->data.low = 0; /* Set to zero. */
+ f->data.high = 0; /* Set to zero. */
+ }
+ else
+ overflow = true;
+ }
+ else
+ overflow = fixed_saturate2 (f->mode, temp_high, temp_low,
+ &f->data, sat_p);
+ }
+ else
+ {
+ /* Unsigned -> Signed. */
+ if (temp_high.high < 0)
+ {
+ if (sat_p)
+ {
+ /* Set to maximum. */
+ f->data.low = -1; /* Set to all ones. */
+ f->data.high = -1; /* Set to all ones. */
+ f->data = double_int_ext (f->data,
+ GET_MODE_FBITS (f->mode)
+ + GET_MODE_IBITS (f->mode),
+ 1); /* Clear the sign. */
+ }
+ else
+ overflow = true;
+ }
+ else
+ overflow = fixed_saturate2 (f->mode, temp_high, temp_low,
+ &f->data, sat_p);
+ }
+ }
+ }
+ else
+ {
+ /* Right shift a to temp based on a->mode. */
+ double_int temp;
+ lshift_double (a->data.low, a->data.high,
+ GET_MODE_FBITS (mode) - GET_MODE_FBITS (a->mode),
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &temp.low, &temp.high,
+ SIGNED_FIXED_POINT_MODE_P (a->mode));
+ f->mode = mode;
+ f->data = temp;
+ if (SIGNED_FIXED_POINT_MODE_P (a->mode) ==
+ SIGNED_FIXED_POINT_MODE_P (f->mode))
+ overflow = fixed_saturate1 (f->mode, f->data, &f->data, sat_p);
+ else
+ {
+ /* Take care of the cases when converting between signed and
+ unsigned. */
+ if (SIGNED_FIXED_POINT_MODE_P (a->mode))
+ {
+ /* Signed -> Unsigned. */
+ if (a->data.high < 0)
+ {
+ if (sat_p)
+ {
+ f->data.low = 0; /* Set to zero. */
+ f->data.high = 0; /* Set to zero. */
+ }
+ else
+ overflow = true;
+ }
+ else
+ overflow = fixed_saturate1 (f->mode, f->data, &f->data, sat_p);
+ }
+ else
+ {
+ /* Unsigned -> Signed. */
+ if (temp.high < 0)
+ {
+ if (sat_p)
+ {
+ /* Set to maximum. */
+ f->data.low = -1; /* Set to all ones. */
+ f->data.high = -1; /* Set to all ones. */
+ f->data = double_int_ext (f->data,
+ GET_MODE_FBITS (f->mode)
+ + GET_MODE_IBITS (f->mode),
+ 1); /* Clear the sign. */
+ }
+ else
+ overflow = true;
+ }
+ else
+ overflow = fixed_saturate1 (f->mode, f->data, &f->data, sat_p);
+ }
+ }
+ }
+
+ f->data = double_int_ext (f->data,
+ SIGNED_FIXED_POINT_MODE_P (f->mode)
+ + GET_MODE_FBITS (f->mode)
+ + GET_MODE_IBITS (f->mode),
+ UNSIGNED_FIXED_POINT_MODE_P (f->mode));
+ return overflow;
+}
+
+/* Convert to a new fixed-point mode from an integer.
+ If UNSIGNED_P, this integer is unsigned.
+ If SAT_P, saturate the result to the max or the min.
+ Return true, if !SAT_P and overflow. */
+
+bool
+fixed_convert_from_int (FIXED_VALUE_TYPE *f, enum machine_mode mode,
+ double_int a, bool unsigned_p, bool sat_p)
+{
+ bool overflow = false;
+ /* Left shift a to temp_high, temp_low. */
+ double_int temp_high, temp_low;
+ int amount = GET_MODE_FBITS (mode);
+ if (amount == 2 * HOST_BITS_PER_WIDE_INT)
+ {
+ temp_high = a;
+ temp_low.low = 0;
+ temp_low.high = 0;
+ }
+ else
+ {
+ lshift_double (a.low, a.high,
+ amount,
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &temp_low.low, &temp_low.high, 0);
+
+ /* Logical shift right to temp_high. */
+ lshift_double (a.low, a.high,
+ amount - 2 * HOST_BITS_PER_WIDE_INT,
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &temp_high.low, &temp_high.high, 0);
+ }
+ if (!unsigned_p && a.high < 0) /* Signed-extend temp_high. */
+ temp_high = double_int_ext (temp_high, amount, 0);
+
+ f->mode = mode;
+ f->data = temp_low;
+
+ if (unsigned_p == UNSIGNED_FIXED_POINT_MODE_P (f->mode))
+ overflow = fixed_saturate2 (f->mode, temp_high, temp_low, &f->data, sat_p);
+ else
+ {
+ /* Take care of the cases when converting between signed and unsigned. */
+ if (!unsigned_p)
+ {
+ /* Signed -> Unsigned. */
+ if (a.high < 0)
+ {
+ if (sat_p)
+ {
+ f->data.low = 0; /* Set to zero. */
+ f->data.high = 0; /* Set to zero. */
+ }
+ else
+ overflow = true;
+ }
+ else
+ overflow = fixed_saturate2 (f->mode, temp_high, temp_low, &f->data,
+ sat_p);
+ }
+ else
+ {
+ /* Unsigned -> Signed. */
+ if (temp_high.high < 0)
+ {
+ if (sat_p)
+ {
+ /* Set to maximum. */
+ f->data.low = -1; /* Set to all ones. */
+ f->data.high = -1; /* Set to all ones. */
+ f->data = double_int_ext (f->data,
+ GET_MODE_FBITS (f->mode)
+ + GET_MODE_IBITS (f->mode),
+ 1); /* Clear the sign. */
+ }
+ else
+ overflow = true;
+ }
+ else
+ overflow = fixed_saturate2 (f->mode, temp_high, temp_low,
+ &f->data, sat_p);
+ }
+ }
+ f->data = double_int_ext (f->data,
+ SIGNED_FIXED_POINT_MODE_P (f->mode)
+ + GET_MODE_FBITS (f->mode)
+ + GET_MODE_IBITS (f->mode),
+ UNSIGNED_FIXED_POINT_MODE_P (f->mode));
+ return overflow;
+}
+
+/* Convert to a new fixed-point mode from a real.
+ If SAT_P, saturate the result to the max or the min.
+ Return true, if !SAT_P and overflow. */
+
+bool
+fixed_convert_from_real (FIXED_VALUE_TYPE *f, enum machine_mode mode,
+ const REAL_VALUE_TYPE *a, bool sat_p)
+{
+ bool overflow = false;
+ REAL_VALUE_TYPE real_value, fixed_value, base_value;
+ bool unsigned_p = UNSIGNED_FIXED_POINT_MODE_P (mode);
+ int i_f_bits = GET_MODE_IBITS (mode) + GET_MODE_FBITS (mode);
+ unsigned int fbit = GET_MODE_FBITS (mode);
+ enum fixed_value_range_code temp;
+
+ real_value = *a;
+ f->mode = mode;
+ real_2expN (&base_value, fbit);
+ real_arithmetic (&fixed_value, MULT_EXPR, &real_value, &base_value);
+ real_to_integer2 ((HOST_WIDE_INT *)&f->data.low, &f->data.high, &fixed_value);
+ temp = check_real_for_fixed_mode (&real_value, mode);
+ if (temp == FIXED_UNDERFLOW) /* Minimum. */
+ {
+ if (sat_p)
+ {
+ if (unsigned_p)
+ {
+ f->data.low = 0;
+ f->data.high = 0;
+ }
+ else
+ {
+ f->data.low = 1;
+ f->data.high = 0;
+ lshift_double (f->data.low, f->data.high, i_f_bits,
+ 2 * HOST_BITS_PER_WIDE_INT,
+ &f->data.low, &f->data.high, 1);
+ f->data = double_int_ext (f->data, 1 + i_f_bits, 0);
+ }
+ }
+ else
+ overflow = true;
+ }
+ else if (temp == FIXED_GT_MAX_EPS || temp == FIXED_MAX_EPS) /* Maximum. */
+ {
+ if (sat_p)
+ {
+ f->data.low = -1;
+ f->data.high = -1;
+ f->data = double_int_ext (f->data, i_f_bits, 1);
+ }
+ else
+ overflow = true;
+ }
+ f->data = double_int_ext (f->data, (!unsigned_p) + i_f_bits, unsigned_p);
+ return overflow;
+}
+
+/* Convert to a new real mode from a fixed-point. */
+
+void
+real_convert_from_fixed (REAL_VALUE_TYPE *r, enum machine_mode mode,
+ const FIXED_VALUE_TYPE *f)
+{
+ REAL_VALUE_TYPE base_value, fixed_value, real_value;
+
+ real_2expN (&base_value, GET_MODE_FBITS (f->mode));
+ real_from_integer (&fixed_value, VOIDmode, f->data.low, f->data.high,
+ UNSIGNED_FIXED_POINT_MODE_P (f->mode));
+ real_arithmetic (&real_value, RDIV_EXPR, &fixed_value, &base_value);
+ real_convert (r, mode, &real_value);
+}
+
+/* Determine whether a fixed-point value F is negative. */
+
+bool
+fixed_isneg (const FIXED_VALUE_TYPE *f)
+{
+ if (SIGNED_FIXED_POINT_MODE_P (f->mode))
+ {
+ int i_f_bits = GET_MODE_IBITS (f->mode) + GET_MODE_FBITS (f->mode);
+ int sign_bit = get_fixed_sign_bit (f->data, i_f_bits);
+ if (sign_bit == 1)
+ return true;
+ }
+
+ return false;
+}
Index: fixed-value.h
===================================================================
--- fixed-value.h (revision 0)
+++ fixed-value.h (revision 0)
@@ -0,0 +1,101 @@
+/* Fixed-point arithmetic support.
+ Copyright (C) 2006, 2007 Free Software Foundation, Inc.
+
+This file is part of GCC.
+
+GCC is free software; you can redistribute it and/or modify it
+under the terms of the GNU General Public License as published by the
+Free Software Foundation; either version 2, or (at your option) any
+later version.
+
+GCC is distributed in the hope that it will be useful, but WITHOUT
+ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
+for more details.
+
+You should have received a copy of the GNU General Public License
+along with GCC; see the file COPYING. If not, write to the Free
+Software Foundation, 51 Franklin Street, Fifth Floor, Boston, MA
+02110-1301, USA. */
+
+#ifndef GCC_FIXED_VALUE_H
+#define GCC_FIXED_VALUE_H
+
+#include "machmode.h"
+#include "real.h"
+#include "double-int.h"
+
+struct fixed_value GTY(())
+{
+ double_int data; /* Store data up to 2 wide integers. */
+ unsigned int mode; /* Use machine mode to know IBIT and FBIT. */
+};
+
+#define FIXED_VALUE_TYPE struct fixed_value
+
+#define MAX_FCONST0 18 /* For storing 18 fixed-point zeros per
+ fract, ufract, accum, and uaccum modes . */
+#define MAX_FCONST1 8 /* For storing 8 fixed-point ones per accum
+ and uaccum modes. */
+/* Constant fixed-point values 0 and 1. */
+extern FIXED_VALUE_TYPE fconst0[MAX_FCONST0];
+extern FIXED_VALUE_TYPE fconst1[MAX_FCONST1];
+
+/* Macros to access fconst0 and fconst1 via machine modes. */
+#define FCONST0(mode) fconst0[mode - QQmode]
+#define FCONST1(mode) fconst1[mode - HAmode]
+
+/* Return a CONST_FIXED with value R and mode M. */
+#define CONST_FIXED_FROM_FIXED_VALUE(r, m) \
+ const_fixed_from_fixed_value (r, m)
+extern rtx const_fixed_from_fixed_value (FIXED_VALUE_TYPE, enum machine_mode);
+
+/* Initialize from a decimal or hexadecimal string. */
+extern void fixed_from_string (FIXED_VALUE_TYPE *, const char *,
+ enum machine_mode);
+
+/* In tree.c: wrap up a FIXED_VALUE_TYPE in a tree node. */
+extern tree build_fixed (tree, FIXED_VALUE_TYPE);
+
+/* Extend or truncate to a new mode. */
+extern bool fixed_convert (FIXED_VALUE_TYPE *, enum machine_mode,
+ const FIXED_VALUE_TYPE *, bool);
+
+/* Convert to a fixed-point mode from an integer. */
+extern bool fixed_convert_from_int (FIXED_VALUE_TYPE *, enum machine_mode,
+ double_int, bool, bool);
+
+/* Convert to a fixed-point mode from a real. */
+extern bool fixed_convert_from_real (FIXED_VALUE_TYPE *, enum machine_mode,
+ const REAL_VALUE_TYPE *, bool);
+
+/* Convert to a real mode from a fixed-point. */
+extern void real_convert_from_fixed (REAL_VALUE_TYPE *, enum machine_mode,
+ const FIXED_VALUE_TYPE *);
+
+/* Compare two fixed-point objects for bitwise identity. */
+extern bool fixed_identical (const FIXED_VALUE_TYPE *, const FIXED_VALUE_TYPE *);
+
+/* Calculate a hash value. */
+extern unsigned int fixed_hash (const FIXED_VALUE_TYPE *);
+
+#define FIXED_VALUES_IDENTICAL(x, y) fixed_identical (&(x), &(y))
+
+/* Determine whether a fixed-point value X is negative. */
+#define FIXED_VALUE_NEGATIVE(x) fixed_isneg (&(x))
+
+/* Render F as a decimal floating point constant. */
+extern void fixed_to_decimal (char *str, const FIXED_VALUE_TYPE *, size_t);
+
+/* Binary or unary arithmetic on tree_code. */
+extern bool fixed_arithmetic (FIXED_VALUE_TYPE *, int, const FIXED_VALUE_TYPE *,
+ const FIXED_VALUE_TYPE *, bool);
+
+/* Compare fixed-point values by tree_code. */
+extern bool fixed_compare (int, const FIXED_VALUE_TYPE *,
+ const FIXED_VALUE_TYPE *);
+
+/* Determine whether a fixed-point value X is negative. */
+extern bool fixed_isneg (const FIXED_VALUE_TYPE *);
+
+#endif /* GCC_FIXED_VALUE_H */
Index: Makefile.in
===================================================================
--- Makefile.in (revision 126045)
+++ Makefile.in (working copy)
@@ -999,6 +999,7 @@
expmed.o \
expr.o \
final.o \
+ fixed-value.o \
fold-const.o \
function.o \
fwprop.o \
@@ -2377,6 +2378,8 @@
toplev.h $(TM_P_H) $(REAL_H)
dfp.o : dfp.c dfp.h $(CONFIG_H) $(SYSTEM_H) coretypes.h $(TM_H) $(TREE_H) \
toplev.h $(TM_P_H) $(REAL_H) $(DECNUM_H)
+fixed-value.o: fixed-value.c $(CONFIG_H) $(SYSTEM_H) coretypes.h $(TM_H) \
+ $(TREE_H) fixed-value.h $(REAL_H) toplev.h
integrate.o : integrate.c $(CONFIG_H) $(SYSTEM_H) coretypes.h $(TM_H) \
$(RTL_H) $(TREE_H) $(FLAGS_H) debug.h $(INTEGRATE_H) insn-config.h \
$(EXPR_H) $(REAL_H) $(REGS_H) intl.h $(FUNCTION_H) output.h $(RECOG_H) \
@@ -3000,6 +3003,7 @@
$(srcdir)/coverage.c $(srcdir)/rtl.h \
$(srcdir)/optabs.h $(srcdir)/tree.h $(srcdir)/function.h $(srcdir)/libfuncs.h $(SYMTAB_H) \
$(srcdir)/real.h $(srcdir)/varray.h $(srcdir)/insn-addr.h $(srcdir)/hwint.h \
+ $(srcdir)/fixed-value.h \
$(srcdir)/ipa-reference.h $(srcdir)/output.h $(srcdir)/cfgloop.h \
$(srcdir)/cselib.h $(srcdir)/basic-block.h $(srcdir)/cgraph.h \
$(srcdir)/reload.h \
^ permalink raw reply [flat|nested] 4+ messages in thread
* RE: [PATCH] The 2nd patch from the fixed-point branch
2007-06-27 1:27 ` Fu, Chao-Ying
@ 2007-07-02 21:10 ` Fu, Chao-Ying
0 siblings, 0 replies; 4+ messages in thread
From: Fu, Chao-Ying @ 2007-07-02 21:10 UTC (permalink / raw)
To: Richard Sandiford, Mark Mitchell
Cc: gcc-patches, Thekkath, Radhika, Stephens, Nigel
Hi,
Any update for this patch? Thanks!
(I have the 3rd patch almost ready to send out.)
Regards,
Chao-ying
> -----Original Message-----
> From: gcc-patches-owner@gcc.gnu.org
> [mailto:gcc-patches-owner@gcc.gnu.org]On Behalf Of Fu, Chao-Ying
> Sent: Tuesday, June 26, 2007 6:24 PM
> To: Richard Sandiford
> Cc: gcc-patches@gcc.gnu.org; Mark Mitchell; Thekkath,
> Radhika; Stephens,
> Nigel
> Subject: RE: [PATCH] The 2nd patch from the fixed-point branch
>
>
> Richard Sandiford wrote:
>
> > > +static int check_real_for_fixed_mode (REAL_VALUE_TYPE *,
> > unsigned int);
> > > +static int get_fixed_sign_bit (double_int, int);
> > > +static bool do_fixed_add (FIXED_VALUE_TYPE *, const
> > FIXED_VALUE_TYPE *,
> > > + const FIXED_VALUE_TYPE *, int, int);
> > > +static bool do_fixed_multiply (FIXED_VALUE_TYPE *, const
> > FIXED_VALUE_TYPE *,
> > > + const FIXED_VALUE_TYPE *, int);
> > > +static bool do_fixed_divide (FIXED_VALUE_TYPE *, const
> > FIXED_VALUE_TYPE *,
> > > + const FIXED_VALUE_TYPE *, int);
> > > +static bool do_fixed_shift (FIXED_VALUE_TYPE *, const
> > FIXED_VALUE_TYPE *,
> > > + const FIXED_VALUE_TYPE *, int, int);
> > > +static bool do_fixed_neg (FIXED_VALUE_TYPE *, const
> > FIXED_VALUE_TYPE *, int);
> > > +static bool fixed_saturate1 (unsigned int, double_int,
> > double_int *, int);
> > > +static bool fixed_saturate2 (unsigned int, double_int,
> double_int,
> > > + double_int *, int);
> >
> > New style is not to have these forward declarations unless they are
> > really needed (due to mutual recursion, etc.)
>
> Yes. These declarations are deleted.
>
> >
> > > +/* Compare two fixed objects for bitwise identity. */
> > > +
> > > +bool
> > > +fixed_identical (const FIXED_VALUE_TYPE *a, const
> > FIXED_VALUE_TYPE *b)
> > > +{
> > > + return a->mode == b->mode && a->data.high == b->data.high
> > > + && a->data.low == b->data.low;
> >
> > Formatting nit, should be:
> >
> > return (a->mode == b->mode
> > && a->data.high == b->data.high
> > && a->data.low == b->data.low);
>
> Yes.
>
> >
> > > +
> > > +/* Check REAL_VALUE against the range of the fixed-point mode.
> > > + Return 0, if it is within the range.
> > > + 1, if it is less than the minimum.
> > > + 2, if it is greater than the maximum, but not equal to
> > > + the maximum + the epsilon.
> > > + 3, if it is equal to the maximum + the epsilon. */
> >
> > Would the code be clearer if this were an enum? Or is this
> mirroring
> > some well-known library function?
>
> Yes. I created an enum for the returned values.
>
> >
> > > +int
> > > +check_real_for_fixed_mode (REAL_VALUE_TYPE *real_value,
> > unsigned int mode)
> >
> > Why is MODE not an enum machine_mode? (I'm sure there's a
> > good reason,
> > I just wasn't sure what it was.)
>
> I just changed to enum machine_mode.
>
> >
> > > +{
> > > + REAL_VALUE_TYPE max_value, min_value, epsilon_value;
> > > + char max_string[20], min_string[20], epsilon_string[20];
> > > +
> > > + sprintf (max_string, "0x1.0p%d", GET_MODE_IBIT (mode));
> > > + sprintf (epsilon_string, "0x1.0p-%d", GET_MODE_FBIT (mode));
> > > + real_from_string (&max_value, max_string);
> > > + real_from_string (&epsilon_value, epsilon_string);
> >
> > Couldn't you use real_2expN here instead? Other cases
> later in file.
>
> Yes. I replaced them with real_2expN.
>
> >
> > > +void fixed_to_decimal (char *str, const FIXED_VALUE_TYPE *f_orig,
> > > + size_t buf_size)
> >
> > Formatting; function name should start a new line.
>
> Yes.
>
> >
> > > +/* If SATP, saturate A to the maximum or the minimum, and
> > save to *F based on
> > > + the machine mode MODE.
> > > + This function assumes the width of double_int is
> > greater than the width
> > > + of the fixed-point value at the fixed-point mode.
> > > + Return true, if !SATP and overflow. */
> > > +
> > > +static bool
> > > +fixed_saturate1 (unsigned int mode, double_int a,
> > double_int *f, int satp)
> >
> > Maybe satp should be a bool too? I don't quite understand what you
> > mean by "the fixed-point value at the fixed-point mode".
>
> SATP is a bool now.
> I try to say that the width is the sum of a possible sign
> bit, possible ibits
> and fbits => without any padding bits.
>
> >
> > It might be clearer if you add "Do not modify *F
> otherwise." after the
> > first sentence. Same comments apply to fixed_saturate2.
>
> Yes.
>
> >
> > > +/* Return the sign bit based on I_F_BITS. */
> > > +
> > > +inline int get_fixed_sign_bit (double_int a, int i_f_bits)
> >
> > Must be "static inline". Formatting: function on new line.
> >
> Yes.
>
> > > + if (((!subtract_p)
> >
> > Redundant brackets.
>
> Yes.
>
> >
> > > + && get_fixed_sign_bit (a->data, i_f_bits)
> > > + == get_fixed_sign_bit (b->data, i_f_bits)
> >
> > Multiline conditions should be bracketed (so that emacs will indent
> > the "==" correctly).
>
> Yes.
>
> >
> > > +/* Calculate F = A << B if LEFT_P. Otherwies, F = A >> B.
> > > + If SATP, saturate the result to the max or the min.
> > > + Return true, if !SATP and overflow. */
> >
> > Mixture of naming styles: LEFT_P vs. SATP. FWIW, I think LEFT_P
> > is the usual convention in gcc. (Same thing with
> SUBTRACT_P earlier.)
>
> I changed SATP to SAT_P.
>
> >
> > > +/* Perform the binary or unary operation described by CODE.
> > > + For a unary operation, leave OP1 NULL.
> > > + Return true, if !SATP and overflow. */
> > > +
> > > +bool
> > > +fixed_arithmetic (FIXED_VALUE_TYPE *f, int icode, const
> > FIXED_VALUE_TYPE *op0,
> > > + const FIXED_VALUE_TYPE *op1, int satp)
> >
> >
> > Am I right in thinking that OP0 and OP1 must have the same mode?
> > If so, it might be worth a comment and a gcc_assert().
>
> Yes.
>
> > The reason I ask is that...
> >
> > > +/* Compare fixed-point values by tree_code. */
> > > +
> > > +bool
> > > +fixed_compare (int icode, const FIXED_VALUE_TYPE *op0,
> > > + const FIXED_VALUE_TYPE *op1)
> > > +{
> > > + enum tree_code code = icode;
> > > +
> > > + switch (code)
> > > + {
> > > + case NE_EXPR:
> > > + return op0->mode != op1->mode
> > > + || !double_int_equal_p (op0->data, op1->data);
> >
> > ...I was somewhat surprised that fixed_arithmentic appeared to need
> > operands of the same mode, but fixed_compare didn't. It also wasn't
> > immediately obvious that values of different modes were treated as
> > unequal, even if they have the same conceptual value.
> Might be worth
> > adding a comment.
>
> "fixed_compare" must have the same mode as well. I changed
> the comment
> and the code to have gcc_assert.
>
> >
> > > +/* Constant fixed-point values 0 and 1. */
> > > +extern FIXED_VALUE_TYPE fconst0[];
> > > +extern FIXED_VALUE_TYPE fconst1[];
> >
> > What are the indexes?
>
> Ok. I created two defines and two new macros to access them.
>
> #define MAX_FCONST0 18 /* For storing 18 fixed-point
> zeros per
> fract, ufract, accum, and
> uaccum modes . */
> #define MAX_FCONST1 8 /* For storing 8 fixed-point
> ones per accum
> and uaccum modes. */
> /* Constant fixed-point values 0 and 1. */
> extern FIXED_VALUE_TYPE fconst0[MAX_FCONST0];
> extern FIXED_VALUE_TYPE fconst1[MAX_FCONST1];
>
> /* Macros to access fconst0 and fconst1 via machine modes. */
> #define FCONST0(mode) fconst0[mode - QQmode]
> #define FCONST1(mode) fconst1[mode - HAmode]
>
>
> >
> > > +/* Determine whether a fixed-point value X is negative. */
> > > +#define FIXED_VALUE_NEGATIVE(x) fixed_isneg (&(x))
> >
> > Odd spacing (in context).
>
> Yes.
>
> The new patch is attached. Thanks a lot!
>
> Regards,
> Chao-ying
>
> 2007-06-25 Chao-ying Fu <fu@mips.com>
>
> * gengtype.c (main): Handle FIXED_VALUE_TYPE type as
> scalar typedef.
> * double-int.c (double_int_scmp): Use casts of unsigned
> HOST_WIDE_INT
> to compare a.low and b.low.
> * fixed-value.c: New file.
> * fixed-value.h: New file.
> * Makefile.in (OBJS-common): Add fixed-value.o.
> (fixed-value.o): New rule.
> (GTFILES): Add fixed-value.h.
>
^ permalink raw reply [flat|nested] 4+ messages in thread
end of thread, other threads:[~2007-07-02 20:30 UTC | newest]
Thread overview: 4+ messages (download: mbox.gz / follow: Atom feed)
-- links below jump to the message on this page --
2007-06-26 1:17 [PATCH] The 2nd patch from the fixed-point branch Fu, Chao-Ying
2007-06-26 15:13 ` Richard Sandiford
2007-06-27 1:27 ` Fu, Chao-Ying
2007-07-02 21:10 ` Fu, Chao-Ying
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).