gcc/fixed-value.c

   1 /* Fixed-point arithmetic support.
   2    Copyright (C) 2006, 2007, 2008, 2009 Free Software Foundation, Inc.
   3
   4 This file is part of GCC.
   5
   6 GCC is free software; you can redistribute it and/or modify it under
   7 the terms of the GNU General Public License as published by the Free
   8 Software Foundation; either version 3, or (at your option) any later
   9 version.
  10
  11 GCC is distributed in the hope that it will be useful, but WITHOUT ANY
  12 WARRANTY; without even the implied warranty of MERCHANTABILITY or
  13 FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
  14 for more details.
  15
  16 You should have received a copy of the GNU General Public License
  17 along with GCC; see the file COPYING3.  If not see
  18 <http://www.gnu.org/licenses/>.  */
  19
  20 #include "config.h"
  21 #include "system.h"
  22 #include "coretypes.h"
  23 #include "tm.h"
  24 #include "tree.h"
  25 #include "toplev.h"
  26 #include "fixed-value.h"
  27
  28 /* Compare two fixed objects for bitwise identity.  */
  29
  30 bool
  31 fixed_identical (const FIXED_VALUE_TYPE *a, const FIXED_VALUE_TYPE *b)
  32 {
  33   return (a->mode == b->mode
  34           && a->data.high == b->data.high
  35           && a->data.low == b->data.low);
  36 }
  37
  38 /* Calculate a hash value.  */
  39
  40 unsigned int
  41 fixed_hash (const FIXED_VALUE_TYPE *f)
  42 {
  43   return (unsigned int) (f->data.low ^ f->data.high);
  44 }
  45
  46 /* Define the enum code for the range of the fixed-point value.  */
  47 enum fixed_value_range_code {
  48   FIXED_OK,             /* The value is within the range.  */
  49   FIXED_UNDERFLOW,      /* The value is less than the minimum.  */
  50   FIXED_GT_MAX_EPS,     /* The value is greater than the maximum, but not equal
  51                            to the maximum plus the epsilon.  */
  52   FIXED_MAX_EPS         /* The value equals the maximum plus the epsilon.  */
  53 };
  54
  55 /* Check REAL_VALUE against the range of the fixed-point mode.
  56    Return FIXED_OK, if it is within the range.
  57           FIXED_UNDERFLOW, if it is less than the minimum.
  58           FIXED_GT_MAX_EPS, if it is greater than the maximum, but not equal to
  59             the maximum plus the epsilon.
  60           FIXED_MAX_EPS, if it is equal to the maximum plus the epsilon.  */
  61
  62 static enum fixed_value_range_code
  63 check_real_for_fixed_mode (REAL_VALUE_TYPE *real_value, enum machine_mode mode)
  64 {
  65   REAL_VALUE_TYPE max_value, min_value, epsilon_value;
  66
  67   real_2expN (&max_value, GET_MODE_IBIT (mode), mode);
  68   real_2expN (&epsilon_value, -GET_MODE_FBIT (mode), mode);
  69
  70   if (SIGNED_FIXED_POINT_MODE_P (mode))
  71     min_value = REAL_VALUE_NEGATE (max_value);
  72   else
  73     real_from_string (&min_value, "0.0");
  74
  75   if (real_compare (LT_EXPR, real_value, &min_value))
  76     return FIXED_UNDERFLOW;
  77   if (real_compare (EQ_EXPR, real_value, &max_value))
  78     return FIXED_MAX_EPS;
  79   real_arithmetic (&max_value, MINUS_EXPR, &max_value, &epsilon_value);
  80   if (real_compare (GT_EXPR, real_value, &max_value))
  81     return FIXED_GT_MAX_EPS;
  82   return FIXED_OK;
  83 }
  84
  85 /* Initialize from a decimal or hexadecimal string.  */
  86
  87 void
  88 fixed_from_string (FIXED_VALUE_TYPE *f, const char *str, enum machine_mode mode)
  89 {
  90   REAL_VALUE_TYPE real_value, fixed_value, base_value;
  91   unsigned int fbit;
  92   enum fixed_value_range_code temp;
  93
  94   f->mode = mode;
  95   fbit = GET_MODE_FBIT (mode);
  96
  97   real_from_string (&real_value, str);
  98   temp = check_real_for_fixed_mode (&real_value, f->mode);
  99   /* We don't want to warn the case when the _Fract value is 1.0.  */
 100   if (temp == FIXED_UNDERFLOW
 101       || temp == FIXED_GT_MAX_EPS
 102       || (temp == FIXED_MAX_EPS && ALL_ACCUM_MODE_P (f->mode)))
 103     warning (OPT_Woverflow,
 104              "large fixed-point constant implicitly truncated to fixed-point type");
 105   real_2expN (&base_value, fbit, mode);
 106   real_arithmetic (&fixed_value, MULT_EXPR, &real_value, &base_value);
 107   real_to_integer2 ((HOST_WIDE_INT *)&f->data.low, &f->data.high,
 108                     &fixed_value);
 109
 110   if (temp == FIXED_MAX_EPS && ALL_FRACT_MODE_P (f->mode))
 111     {
 112       /* From the spec, we need to evaluate 1 to the maximal value.  */
 113       f->data.low = -1;
 114       f->data.high = -1;
 115       f->data = double_int_ext (f->data,
 116                                 GET_MODE_FBIT (f->mode)
 117                                 + GET_MODE_IBIT (f->mode), 1);
 118     }
 119   else
 120     f->data = double_int_ext (f->data,
 121                               SIGNED_FIXED_POINT_MODE_P (f->mode)
 122                               + GET_MODE_FBIT (f->mode)
 123                               + GET_MODE_IBIT (f->mode),
 124                               UNSIGNED_FIXED_POINT_MODE_P (f->mode));
 125 }
 126
 127 /* Render F as a decimal floating point constant.  */
 128
 129 void
 130 fixed_to_decimal (char *str, const FIXED_VALUE_TYPE *f_orig,
 131                   size_t buf_size)
 132 {
 133   REAL_VALUE_TYPE real_value, base_value, fixed_value;
 134
 135   real_2expN (&base_value, GET_MODE_FBIT (f_orig->mode), f_orig->mode);
 136   real_from_integer (&real_value, VOIDmode, f_orig->data.low, f_orig->data.high,
 137                      UNSIGNED_FIXED_POINT_MODE_P (f_orig->mode));
 138   real_arithmetic (&fixed_value, RDIV_EXPR, &real_value, &base_value);
 139   real_to_decimal (str, &fixed_value, buf_size, 0, 1);
 140 }
 141
 142 /* If SAT_P, saturate A to the maximum or the minimum, and save to *F based on
 143    the machine mode MODE.
 144    Do not modify *F otherwise.
 145    This function assumes the width of double_int is greater than the width
 146    of the fixed-point value (the sum of a possible sign bit, possible ibits,
 147    and fbits).
 148    Return true, if !SAT_P and overflow.  */
 149
 150 static bool
 151 fixed_saturate1 (enum machine_mode mode, double_int a, double_int *f,
 152                  bool sat_p)
 153 {
 154   bool overflow_p = false;
 155   bool unsigned_p = UNSIGNED_FIXED_POINT_MODE_P (mode);
 156   int i_f_bits = GET_MODE_IBIT (mode) + GET_MODE_FBIT (mode);
 157
 158   if (unsigned_p) /* Unsigned type.  */
 159     {
 160       double_int max;
 161       max.low = -1;
 162       max.high = -1;
 163       max = double_int_ext (max, i_f_bits, 1);
 164       if (double_int_cmp (a, max, 1) == 1)
 165         {
 166           if (sat_p)
 167             *f = max;
 168           else
 169             overflow_p = true;
 170         }
 171     }
 172   else /* Signed type.  */
 173     {
 174       double_int max, min;
 175       max.high = -1;
 176       max.low = -1;
 177       max = double_int_ext (max, i_f_bits, 1);
 178       min.high = 0;
 179       min.low = 1;
 180       lshift_double (min.low, min.high, i_f_bits,
 181                      2 * HOST_BITS_PER_WIDE_INT,
 182                      &min.low, &min.high, 1);
 183       min = double_int_ext (min, 1 + i_f_bits, 0);
 184       if (double_int_cmp (a, max, 0) == 1)
 185         {
 186           if (sat_p)
 187             *f = max;
 188           else
 189             overflow_p = true;
 190         }
 191       else if (double_int_cmp (a, min, 0) == -1)
 192         {
 193           if (sat_p)
 194             *f = min;
 195           else
 196             overflow_p = true;
 197         }
 198     }
 199   return overflow_p;
 200 }
 201
 202 /* If SAT_P, saturate {A_HIGH, A_LOW} to the maximum or the minimum, and
 203    save to *F based on the machine mode MODE.
 204    Do not modify *F otherwise.
 205    This function assumes the width of two double_int is greater than the width
 206    of the fixed-point value (the sum of a possible sign bit, possible ibits,
 207    and fbits).
 208    Return true, if !SAT_P and overflow.  */
 209
 210 static bool
 211 fixed_saturate2 (enum machine_mode mode, double_int a_high, double_int a_low,
 212                  double_int *f, bool sat_p)
 213 {
 214   bool overflow_p = false;
 215   bool unsigned_p = UNSIGNED_FIXED_POINT_MODE_P (mode);
 216   int i_f_bits = GET_MODE_IBIT (mode) + GET_MODE_FBIT (mode);
 217
 218   if (unsigned_p) /* Unsigned type.  */
 219     {
 220       double_int max_r, max_s;
 221       max_r.high = 0;
 222       max_r.low = 0;
 223       max_s.high = -1;
 224       max_s.low = -1;
 225       max_s = double_int_ext (max_s, i_f_bits, 1);
 226       if (double_int_cmp (a_high, max_r, 1) == 1
 227           || (double_int_equal_p (a_high, max_r) &&
 228               double_int_cmp (a_low, max_s, 1) == 1))
 229         {
 230           if (sat_p)
 231             *f = max_s;
 232           else
 233             overflow_p = true;
 234         }
 235     }
 236   else /* Signed type.  */
 237     {
 238       double_int max_r, max_s, min_r, min_s;
 239       max_r.high = 0;
 240       max_r.low = 0;
 241       max_s.high = -1;
 242       max_s.low = -1;
 243       max_s = double_int_ext (max_s, i_f_bits, 1);
 244       min_r.high = -1;
 245       min_r.low = -1;
 246       min_s.high = 0;
 247       min_s.low = 1;
 248       lshift_double (min_s.low, min_s.high, i_f_bits,
 249                      2 * HOST_BITS_PER_WIDE_INT,
 250                      &min_s.low, &min_s.high, 1);
 251       min_s = double_int_ext (min_s, 1 + i_f_bits, 0);
 252       if (double_int_cmp (a_high, max_r, 0) == 1
 253           || (double_int_equal_p (a_high, max_r) &&
 254               double_int_cmp (a_low, max_s, 1) == 1))
 255         {
 256           if (sat_p)
 257             *f = max_s;
 258           else
 259             overflow_p = true;
 260         }
 261       else if (double_int_cmp (a_high, min_r, 0) == -1
 262                || (double_int_equal_p (a_high, min_r) &&
 263                    double_int_cmp (a_low, min_s, 1) == -1))
 264         {
 265           if (sat_p)
 266             *f = min_s;
 267           else
 268             overflow_p = true;
 269         }
 270     }
 271   return overflow_p;
 272 }
 273
 274 /* Return the sign bit based on I_F_BITS.  */
 275
 276 static inline int
 277 get_fixed_sign_bit (double_int a, int i_f_bits)
 278 {
 279   if (i_f_bits < HOST_BITS_PER_WIDE_INT)
 280     return (a.low >> i_f_bits) & 1;
 281   else
 282     return (a.high >> (i_f_bits - HOST_BITS_PER_WIDE_INT)) & 1;
 283 }
 284
 285 /* Calculate F = A + (SUBTRACT_P ? -B : B).
 286    If SAT_P, saturate the result to the max or the min.
 287    Return true, if !SAT_P and overflow.  */
 288
 289 static bool
 290 do_fixed_add (FIXED_VALUE_TYPE *f, const FIXED_VALUE_TYPE *a,
 291               const FIXED_VALUE_TYPE *b, bool subtract_p, bool sat_p)
 292 {
 293   bool overflow_p = false;
 294   bool unsigned_p;
 295   double_int temp;
 296   int i_f_bits;
 297
 298   /* This was a conditional expression but it triggered a bug in
 299      Sun C 5.5.  */
 300   if (subtract_p)
 301     temp = double_int_neg (b->data);
 302   else
 303     temp = b->data;
 304
 305   unsigned_p = UNSIGNED_FIXED_POINT_MODE_P (a->mode);
 306   i_f_bits = GET_MODE_IBIT (a->mode) + GET_MODE_FBIT (a->mode);
 307   f->mode = a->mode;
 308   f->data = double_int_add (a->data, temp);
 309   if (unsigned_p) /* Unsigned type.  */
 310     {
 311       if (subtract_p) /* Unsigned subtraction.  */
 312         {
 313           if (double_int_cmp (a->data, b->data, 1) == -1)
 314             {
 315               if (sat_p)
 316                 {
 317                   f->data.high = 0;
 318                   f->data.low = 0;
 319                  }
 320               else
 321                 overflow_p = true;
 322             }
 323         }
 324       else /* Unsigned addition.  */
 325         {
 326           f->data = double_int_ext (f->data, i_f_bits, 1);
 327           if (double_int_cmp (f->data, a->data, 1) == -1
 328               || double_int_cmp (f->data, b->data, 1) == -1)
 329             {
 330               if (sat_p)
 331                 {
 332                   f->data.high = -1;
 333                   f->data.low = -1;
 334                 }
 335               else
 336                 overflow_p = true;
 337             }
 338         }
 339     }
 340   else /* Signed type.  */
 341     {
 342       if ((!subtract_p
 343            && (get_fixed_sign_bit (a->data, i_f_bits)
 344                == get_fixed_sign_bit (b->data, i_f_bits))
 345            && (get_fixed_sign_bit (a->data, i_f_bits)
 346                != get_fixed_sign_bit (f->data, i_f_bits)))
 347           || (subtract_p
 348               && (get_fixed_sign_bit (a->data, i_f_bits)
 349                   != get_fixed_sign_bit (b->data, i_f_bits))
 350               && (get_fixed_sign_bit (a->data, i_f_bits)
 351                   != get_fixed_sign_bit (f->data, i_f_bits))))
 352         {
 353           if (sat_p)
 354             {
 355               f->data.low = 1;
 356               f->data.high = 0;
 357               lshift_double (f->data.low, f->data.high, i_f_bits,
 358                              2 * HOST_BITS_PER_WIDE_INT,
 359                              &f->data.low, &f->data.high, 1);
 360               if (get_fixed_sign_bit (a->data, i_f_bits) == 0)
 361                 {
 362                   double_int one;
 363                   one.low = 1;
 364                   one.high = 0;
 365                   f->data = double_int_add (f->data, double_int_neg (one));
 366                 }
 367             }
 368           else
 369             overflow_p = true;
 370         }
 371     }
 372   f->data = double_int_ext (f->data, (!unsigned_p) + i_f_bits, unsigned_p);
 373   return overflow_p;
 374 }
 375
 376 /* Calculate F = A * B.
 377    If SAT_P, saturate the result to the max or the min.
 378    Return true, if !SAT_P and overflow.  */
 379
 380 static bool
 381 do_fixed_multiply (FIXED_VALUE_TYPE *f, const FIXED_VALUE_TYPE *a,
 382                    const FIXED_VALUE_TYPE *b, bool sat_p)
 383 {
 384   bool overflow_p = false;
 385   bool unsigned_p = UNSIGNED_FIXED_POINT_MODE_P (a->mode);
 386   int i_f_bits = GET_MODE_IBIT (a->mode) + GET_MODE_FBIT (a->mode);
 387   f->mode = a->mode;
 388   if (GET_MODE_PRECISION (f->mode) <= HOST_BITS_PER_WIDE_INT)
 389     {
 390       f->data = double_int_mul (a->data, b->data);
 391       lshift_double (f->data.low, f->data.high,
 392                      (-GET_MODE_FBIT (f->mode)),
 393                      2 * HOST_BITS_PER_WIDE_INT,
 394                      &f->data.low, &f->data.high, !unsigned_p);
 395       overflow_p = fixed_saturate1 (f->mode, f->data, &f->data, sat_p);
 396     }
 397   else
 398     {
 399       /* The result of multiplication expands to two double_int.  */
 400       double_int a_high, a_low, b_high, b_low;
 401       double_int high_high, high_low, low_high, low_low;
 402       double_int r, s, temp1, temp2;
 403       int carry = 0;
 404
 405       /* Decompose a and b to four double_int.  */
 406       a_high.low = a->data.high;
 407       a_high.high = 0;
 408       a_low.low = a->data.low;
 409       a_low.high = 0;
 410       b_high.low = b->data.high;
 411       b_high.high = 0;
 412       b_low.low = b->data.low;
 413       b_low.high = 0;
 414
 415       /* Perform four multiplications.  */
 416       low_low = double_int_mul (a_low, b_low);
 417       low_high = double_int_mul (a_low, b_high);
 418       high_low = double_int_mul (a_high, b_low);
 419       high_high = double_int_mul (a_high, b_high);
 420
 421       /* Accumulate four results to {r, s}.  */
 422       temp1.high = high_low.low;
 423       temp1.low = 0;
 424       s = double_int_add (low_low, temp1);
 425       if (double_int_cmp (s, low_low, 1) == -1
 426           || double_int_cmp (s, temp1, 1) == -1)
 427         carry ++; /* Carry */
 428       temp1.high = s.high;
 429       temp1.low = s.low;
 430       temp2.high = low_high.low;
 431       temp2.low = 0;
 432       s = double_int_add (temp1, temp2);
 433       if (double_int_cmp (s, temp1, 1) == -1
 434           || double_int_cmp (s, temp2, 1) == -1)
 435         carry ++; /* Carry */
 436
 437       temp1.low = high_low.high;
 438       temp1.high = 0;
 439       r = double_int_add (high_high, temp1);
 440       temp1.low = low_high.high;
 441       temp1.high = 0;
 442       r = double_int_add (r, temp1);
 443       temp1.low = carry;
 444       temp1.high = 0;
 445       r = double_int_add (r, temp1);
 446
 447       /* We need to add neg(b) to r, if a < 0.  */
 448       if (!unsigned_p && a->data.high < 0)
 449         r = double_int_add (r, double_int_neg (b->data));
 450       /* We need to add neg(a) to r, if b < 0.  */
 451       if (!unsigned_p && b->data.high < 0)
 452         r = double_int_add (r, double_int_neg (a->data));
 453
 454       /* Shift right the result by FBIT.  */
 455       if (GET_MODE_FBIT (f->mode) == 2 * HOST_BITS_PER_WIDE_INT)
 456         {
 457           s.low = r.low;
 458           s.high = r.high;
 459           if (unsigned_p)
 460             {
 461               r.low = 0;
 462               r.high = 0;
 463             }
 464           else
 465             {
 466               r.low = -1;
 467               r.high = -1;
 468             }
 469           f->data.low = s.low;
 470           f->data.high = s.high;
 471         }
 472       else
 473         {
 474           lshift_double (s.low, s.high,
 475                          (-GET_MODE_FBIT (f->mode)),
 476                          2 * HOST_BITS_PER_WIDE_INT,
 477                          &s.low, &s.high, 0);
 478           lshift_double (r.low, r.high,
 479                          (2 * HOST_BITS_PER_WIDE_INT
 480                           - GET_MODE_FBIT (f->mode)),
 481                          2 * HOST_BITS_PER_WIDE_INT,
 482                          &f->data.low, &f->data.high, 0);
 483           f->data.low = f->data.low | s.low;
 484           f->data.high = f->data.high | s.high;
 485           s.low = f->data.low;
 486           s.high = f->data.high;
 487           lshift_double (r.low, r.high,
 488                          (-GET_MODE_FBIT (f->mode)),
 489                          2 * HOST_BITS_PER_WIDE_INT,
 490                          &r.low, &r.high, !unsigned_p);
 491         }
 492
 493       overflow_p = fixed_saturate2 (f->mode, r, s, &f->data, sat_p);
 494     }
 495
 496   f->data = double_int_ext (f->data, (!unsigned_p) + i_f_bits, unsigned_p);
 497   return overflow_p;
 498 }
 499
 500 /* Calculate F = A / B.
 501    If SAT_P, saturate the result to the max or the min.
 502    Return true, if !SAT_P and overflow.  */
 503
 504 static bool
 505 do_fixed_divide (FIXED_VALUE_TYPE *f, const FIXED_VALUE_TYPE *a,
 506                  const FIXED_VALUE_TYPE *b, bool sat_p)
 507 {
 508   bool overflow_p = false;
 509   bool unsigned_p = UNSIGNED_FIXED_POINT_MODE_P (a->mode);
 510   int i_f_bits = GET_MODE_IBIT (a->mode) + GET_MODE_FBIT (a->mode);
 511   f->mode = a->mode;
 512   if (GET_MODE_PRECISION (f->mode) <= HOST_BITS_PER_WIDE_INT)
 513     {
 514       lshift_double (a->data.low, a->data.high,
 515                      GET_MODE_FBIT (f->mode),
 516                      2 * HOST_BITS_PER_WIDE_INT,
 517                      &f->data.low, &f->data.high, !unsigned_p);
 518       f->data = double_int_div (f->data, b->data, unsigned_p, TRUNC_DIV_EXPR);
 519       overflow_p = fixed_saturate1 (f->mode, f->data, &f->data, sat_p);
 520     }
 521   else
 522     {
 523       double_int pos_a, pos_b, r, s;
 524       double_int quo_r, quo_s, mod, temp;
 525       int num_of_neg = 0;
 526       int i;
 527
 528       /* If a < 0, negate a.  */
 529       if (!unsigned_p && a->data.high < 0)
 530         {
 531           pos_a = double_int_neg (a->data);
 532           num_of_neg ++;
 533         }
 534       else
 535         pos_a = a->data;
 536
 537       /* If b < 0, negate b.  */
 538       if (!unsigned_p && b->data.high < 0)
 539         {
 540           pos_b = double_int_neg (b->data);
 541           num_of_neg ++;
 542         }
 543       else
 544         pos_b = b->data;
 545
 546       /* Left shift pos_a to {r, s} by FBIT.  */
 547       if (GET_MODE_FBIT (f->mode) == 2 * HOST_BITS_PER_WIDE_INT)
 548         {
 549           r = pos_a;
 550           s.high = 0;
 551           s.low = 0;
 552         }
 553       else
 554         {
 555           lshift_double (pos_a.low, pos_a.high,
 556                          GET_MODE_FBIT (f->mode),
 557                          2 * HOST_BITS_PER_WIDE_INT,
 558                          &s.low, &s.high, 0);
 559           lshift_double (pos_a.low, pos_a.high,
 560                          - (2 * HOST_BITS_PER_WIDE_INT
 561                             - GET_MODE_FBIT (f->mode)),
 562                          2 * HOST_BITS_PER_WIDE_INT,
 563                          &r.low, &r.high, 0);
 564         }
 565
 566       /* Divide r by pos_b to quo_r.  The remainder is in mod.  */
 567       div_and_round_double (TRUNC_DIV_EXPR, 1, r.low, r.high, pos_b.low,
 568                             pos_b.high, &quo_r.low, &quo_r.high, &mod.low,
 569                             &mod.high);
 570
 571       quo_s.high = 0;
 572       quo_s.low = 0;
 573
 574       for (i = 0; i < 2 * HOST_BITS_PER_WIDE_INT; i++)
 575         {
 576           /* Record the leftmost bit of mod.  */
 577           int leftmost_mod = (mod.high < 0);
 578
 579           /* Shift left mod by 1 bit.  */
 580           lshift_double (mod.low, mod.high, 1, 2 * HOST_BITS_PER_WIDE_INT,
 581                          &mod.low, &mod.high, 0);
 582
 583           /* Test the leftmost bit of s to add to mod.  */
 584           if (s.high < 0)
 585             mod.low += 1;
 586
 587           /* Shift left quo_s by 1 bit.  */
 588           lshift_double (quo_s.low, quo_s.high, 1, 2 * HOST_BITS_PER_WIDE_INT,
 589                          &quo_s.low, &quo_s.high, 0);
 590
 591           /* Try to calculate (mod - pos_b).  */
 592           temp = double_int_add (mod, double_int_neg (pos_b));
 593
 594           if (leftmost_mod == 1 || double_int_cmp (mod, pos_b, 1) != -1)
 595             {
 596               quo_s.low += 1;
 597               mod = temp;
 598             }
 599
 600           /* Shift left s by 1 bit.  */
 601           lshift_double (s.low, s.high, 1, 2 * HOST_BITS_PER_WIDE_INT,
 602                          &s.low, &s.high, 0);
 603
 604         }
 605
 606       if (num_of_neg == 1)
 607         {
 608           quo_s = double_int_neg (quo_s);
 609           if (quo_s.high == 0 && quo_s.low == 0)
 610             quo_r = double_int_neg (quo_r);
 611           else
 612             {
 613               quo_r.low = ~quo_r.low;
 614               quo_r.high = ~quo_r.high;
 615             }
 616         }
 617
 618       f->data = quo_s;
 619       overflow_p = fixed_saturate2 (f->mode, quo_r, quo_s, &f->data, sat_p);
 620     }
 621
 622   f->data = double_int_ext (f->data, (!unsigned_p) + i_f_bits, unsigned_p);
 623   return overflow_p;
 624 }
 625
 626 /* Calculate F = A << B if LEFT_P.  Otherwise, F = A >> B.
 627    If SAT_P, saturate the result to the max or the min.
 628    Return true, if !SAT_P and overflow.  */
 629
 630 static bool
 631 do_fixed_shift (FIXED_VALUE_TYPE *f, const FIXED_VALUE_TYPE *a,
 632               const FIXED_VALUE_TYPE *b, bool left_p, bool sat_p)
 633 {
 634   bool overflow_p = false;
 635   bool unsigned_p = UNSIGNED_FIXED_POINT_MODE_P (a->mode);
 636   int i_f_bits = GET_MODE_IBIT (a->mode) + GET_MODE_FBIT (a->mode);
 637   f->mode = a->mode;
 638
 639   if (b->data.low == 0)
 640     {
 641       f->data = a->data;
 642       return overflow_p;
 643     }
 644
 645   if (GET_MODE_PRECISION (f->mode) <= HOST_BITS_PER_WIDE_INT || (!left_p))
 646     {
 647       lshift_double (a->data.low, a->data.high,
 648                      left_p ? b->data.low : (-b->data.low),
 649                      2 * HOST_BITS_PER_WIDE_INT,
 650                      &f->data.low, &f->data.high, !unsigned_p);
 651       if (left_p) /* Only left shift saturates.  */
 652         overflow_p = fixed_saturate1 (f->mode, f->data, &f->data, sat_p);
 653     }
 654   else /* We need two double_int to store the left-shift result.  */
 655     {
 656       double_int temp_high, temp_low;
 657       if (b->data.low == 2 * HOST_BITS_PER_WIDE_INT)
 658         {
 659           temp_high = a->data;
 660           temp_low.high = 0;
 661           temp_low.low = 0;
 662         }
 663       else
 664         {
 665           lshift_double (a->data.low, a->data.high,
 666                          b->data.low,
 667                          2 * HOST_BITS_PER_WIDE_INT,
 668                          &temp_low.low, &temp_low.high, !unsigned_p);
 669           /* Logical shift right to temp_high.  */
 670           lshift_double (a->data.low, a->data.high,
 671                          b->data.low - 2 * HOST_BITS_PER_WIDE_INT,
 672                          2 * HOST_BITS_PER_WIDE_INT,
 673                          &temp_high.low, &temp_high.high, 0);
 674         }
 675       if (!unsigned_p && a->data.high < 0) /* Signed-extend temp_high.  */
 676         temp_high = double_int_ext (temp_high, b->data.low, unsigned_p);
 677       f->data = temp_low;
 678       overflow_p = fixed_saturate2 (f->mode, temp_high, temp_low, &f->data,
 679                                     sat_p);
 680     }
 681   f->data = double_int_ext (f->data, (!unsigned_p) + i_f_bits, unsigned_p);
 682   return overflow_p;
 683 }
 684
 685 /* Calculate F = -A.
 686    If SAT_P, saturate the result to the max or the min.
 687    Return true, if !SAT_P and overflow.  */
 688
 689 static bool
 690 do_fixed_neg (FIXED_VALUE_TYPE *f, const FIXED_VALUE_TYPE *a, bool sat_p)
 691 {
 692   bool overflow_p = false;
 693   bool unsigned_p = UNSIGNED_FIXED_POINT_MODE_P (a->mode);
 694   int i_f_bits = GET_MODE_IBIT (a->mode) + GET_MODE_FBIT (a->mode);
 695   f->mode = a->mode;
 696   f->data = double_int_neg (a->data);
 697   f->data = double_int_ext (f->data, (!unsigned_p) + i_f_bits, unsigned_p);
 698
 699   if (unsigned_p) /* Unsigned type.  */
 700     {
 701       if (f->data.low != 0 || f->data.high != 0)
 702         {
 703           if (sat_p)
 704             {
 705               f->data.low = 0;
 706               f->data.high = 0;
 707             }
 708           else
 709             overflow_p = true;
 710         }
 711     }
 712   else /* Signed type.  */
 713     {
 714       if (!(f->data.high == 0 && f->data.low == 0)
 715           && f->data.high == a->data.high && f->data.low == a->data.low )
 716         {
 717           if (sat_p)
 718             {
 719               /* Saturate to the maximum by subtracting f->data by one.  */
 720               f->data.low = -1;
 721               f->data.high = -1;
 722               f->data = double_int_ext (f->data, i_f_bits, 1);
 723             }
 724           else
 725             overflow_p = true;
 726         }
 727     }
 728   return overflow_p;
 729 }
 730
 731 /* Perform the binary or unary operation described by CODE.
 732    Note that OP0 and OP1 must have the same mode for binary operators.
 733    For a unary operation, leave OP1 NULL.
 734    Return true, if !SAT_P and overflow.  */
 735
 736 bool
 737 fixed_arithmetic (FIXED_VALUE_TYPE *f, int icode, const FIXED_VALUE_TYPE *op0,
 738                   const FIXED_VALUE_TYPE *op1, bool sat_p)
 739 {
 740   switch (icode)
 741     {
 742     case NEGATE_EXPR:
 743       return do_fixed_neg (f, op0, sat_p);
 744       break;
 745
 746     case PLUS_EXPR:
 747       gcc_assert (op0->mode == op1->mode);
 748       return do_fixed_add (f, op0, op1, false, sat_p);
 749       break;
 750
 751     case MINUS_EXPR:
 752       gcc_assert (op0->mode == op1->mode);
 753       return do_fixed_add (f, op0, op1, true, sat_p);
 754       break;
 755
 756     case MULT_EXPR:
 757       gcc_assert (op0->mode == op1->mode);
 758       return do_fixed_multiply (f, op0, op1, sat_p);
 759       break;
 760
 761     case TRUNC_DIV_EXPR:
 762       gcc_assert (op0->mode == op1->mode);
 763       return do_fixed_divide (f, op0, op1, sat_p);
 764       break;
 765
 766     case LSHIFT_EXPR:
 767       return do_fixed_shift (f, op0, op1, true, sat_p);
 768       break;
 769
 770     case RSHIFT_EXPR:
 771       return do_fixed_shift (f, op0, op1, false, sat_p);
 772       break;
 773
 774     default:
 775       gcc_unreachable ();
 776     }
 777   return false;
 778 }
 779
 780 /* Compare fixed-point values by tree_code.
 781    Note that OP0 and OP1 must have the same mode.  */
 782
 783 bool
 784 fixed_compare (int icode, const FIXED_VALUE_TYPE *op0,
 785                const FIXED_VALUE_TYPE *op1)
 786 {
 787   enum tree_code code = (enum tree_code) icode;
 788   gcc_assert (op0->mode == op1->mode);
 789
 790   switch (code)
 791     {
 792     case NE_EXPR:
 793       return !double_int_equal_p (op0->data, op1->data);
 794
 795     case EQ_EXPR:
 796       return double_int_equal_p (op0->data, op1->data);
 797
 798     case LT_EXPR:
 799       return double_int_cmp (op0->data, op1->data,
 800                              UNSIGNED_FIXED_POINT_MODE_P (op0->mode)) == -1;
 801
 802     case LE_EXPR:
 803       return double_int_cmp (op0->data, op1->data,
 804                              UNSIGNED_FIXED_POINT_MODE_P (op0->mode)) != 1;
 805
 806     case GT_EXPR:
 807       return double_int_cmp (op0->data, op1->data,
 808                              UNSIGNED_FIXED_POINT_MODE_P (op0->mode)) == 1;
 809
 810     case GE_EXPR:
 811       return double_int_cmp (op0->data, op1->data,
 812                              UNSIGNED_FIXED_POINT_MODE_P (op0->mode)) != -1;
 813
 814     default:
 815       gcc_unreachable ();
 816     }
 817 }
 818
 819 /* Extend or truncate to a new mode.
 820    If SAT_P, saturate the result to the max or the min.
 821    Return true, if !SAT_P and overflow.  */
 822
 823 bool
 824 fixed_convert (FIXED_VALUE_TYPE *f, enum machine_mode mode,
 825                const FIXED_VALUE_TYPE *a, bool sat_p)
 826 {
 827   bool overflow_p = false;
 828   if (mode == a->mode)
 829     {
 830       *f = *a;
 831       return overflow_p;
 832     }
 833
 834   if (GET_MODE_FBIT (mode) > GET_MODE_FBIT (a->mode))
 835     {
 836       /* Left shift a to temp_high, temp_low based on a->mode.  */
 837       double_int temp_high, temp_low;
 838       int amount = GET_MODE_FBIT (mode) - GET_MODE_FBIT (a->mode);
 839       lshift_double (a->data.low, a->data.high,
 840                      amount,
 841                      2 * HOST_BITS_PER_WIDE_INT,
 842                      &temp_low.low, &temp_low.high,
 843                      SIGNED_FIXED_POINT_MODE_P (a->mode));
 844       /* Logical shift right to temp_high.  */
 845       lshift_double (a->data.low, a->data.high,
 846                      amount - 2 * HOST_BITS_PER_WIDE_INT,
 847                      2 * HOST_BITS_PER_WIDE_INT,
 848                      &temp_high.low, &temp_high.high, 0);
 849       if (SIGNED_FIXED_POINT_MODE_P (a->mode)
 850           && a->data.high < 0) /* Signed-extend temp_high.  */
 851         temp_high = double_int_ext (temp_high, amount, 0);
 852       f->mode = mode;
 853       f->data = temp_low;
 854       if (SIGNED_FIXED_POINT_MODE_P (a->mode) ==
 855           SIGNED_FIXED_POINT_MODE_P (f->mode))
 856         overflow_p = fixed_saturate2 (f->mode, temp_high, temp_low, &f->data,
 857                                       sat_p);
 858       else
 859         {
 860           /* Take care of the cases when converting between signed and
 861              unsigned.  */
 862           if (SIGNED_FIXED_POINT_MODE_P (a->mode))
 863             {
 864               /* Signed -> Unsigned.  */
 865               if (a->data.high < 0)
 866                 {
 867                   if (sat_p)
 868                     {
 869                       f->data.low = 0;  /* Set to zero.  */
 870                       f->data.high = 0;  /* Set to zero.  */
 871                     }
 872                   else
 873                     overflow_p = true;
 874                 }
 875               else
 876                 overflow_p = fixed_saturate2 (f->mode, temp_high, temp_low,
 877                                               &f->data, sat_p);
 878             }
 879           else
 880             {
 881               /* Unsigned -> Signed.  */
 882               if (temp_high.high < 0)
 883                 {
 884                   if (sat_p)
 885                     {
 886                       /* Set to maximum.  */
 887                       f->data.low = -1;  /* Set to all ones.  */
 888                       f->data.high = -1;  /* Set to all ones.  */
 889                       f->data = double_int_ext (f->data,
 890                                                 GET_MODE_FBIT (f->mode)
 891                                                 + GET_MODE_IBIT (f->mode),
 892                                                 1); /* Clear the sign.  */
 893                     }
 894                   else
 895                     overflow_p = true;
 896                 }
 897               else
 898                 overflow_p = fixed_saturate2 (f->mode, temp_high, temp_low,
 899                                               &f->data, sat_p);
 900             }
 901         }
 902     }
 903   else
 904     {
 905       /* Right shift a to temp based on a->mode.  */
 906       double_int temp;
 907       lshift_double (a->data.low, a->data.high,
 908                      GET_MODE_FBIT (mode) - GET_MODE_FBIT (a->mode),
 909                      2 * HOST_BITS_PER_WIDE_INT,
 910                      &temp.low, &temp.high,
 911                      SIGNED_FIXED_POINT_MODE_P (a->mode));
 912       f->mode = mode;
 913       f->data = temp;
 914       if (SIGNED_FIXED_POINT_MODE_P (a->mode) ==
 915           SIGNED_FIXED_POINT_MODE_P (f->mode))
 916         overflow_p = fixed_saturate1 (f->mode, f->data, &f->data, sat_p);
 917       else
 918         {
 919           /* Take care of the cases when converting between signed and
 920              unsigned.  */
 921           if (SIGNED_FIXED_POINT_MODE_P (a->mode))
 922             {
 923               /* Signed -> Unsigned.  */
 924               if (a->data.high < 0)
 925                 {
 926                   if (sat_p)
 927                     {
 928                       f->data.low = 0;  /* Set to zero.  */
 929                       f->data.high = 0;  /* Set to zero.  */
 930                     }
 931                   else
 932                     overflow_p = true;
 933                 }
 934               else
 935                 overflow_p = fixed_saturate1 (f->mode, f->data, &f->data,
 936                                               sat_p);
 937             }
 938           else
 939             {
 940               /* Unsigned -> Signed.  */
 941               if (temp.high < 0)
 942                 {
 943                   if (sat_p)
 944                     {
 945                       /* Set to maximum.  */
 946                       f->data.low = -1;  /* Set to all ones.  */
 947                       f->data.high = -1;  /* Set to all ones.  */
 948                       f->data = double_int_ext (f->data,
 949                                                 GET_MODE_FBIT (f->mode)
 950                                                 + GET_MODE_IBIT (f->mode),
 951                                                 1); /* Clear the sign.  */
 952                     }
 953                   else
 954                     overflow_p = true;
 955                 }
 956               else
 957                 overflow_p = fixed_saturate1 (f->mode, f->data, &f->data,
 958                                               sat_p);
 959             }
 960         }
 961     }
 962
 963   f->data = double_int_ext (f->data,
 964                             SIGNED_FIXED_POINT_MODE_P (f->mode)
 965                             + GET_MODE_FBIT (f->mode)
 966                             + GET_MODE_IBIT (f->mode),
 967                             UNSIGNED_FIXED_POINT_MODE_P (f->mode));
 968   return overflow_p;
 969 }
 970
 971 /* Convert to a new fixed-point mode from an integer.
 972    If UNSIGNED_P, this integer is unsigned.
 973    If SAT_P, saturate the result to the max or the min.
 974    Return true, if !SAT_P and overflow.  */
 975
 976 bool
 977 fixed_convert_from_int (FIXED_VALUE_TYPE *f, enum machine_mode mode,
 978                         double_int a, bool unsigned_p, bool sat_p)
 979 {
 980   bool overflow_p = false;
 981   /* Left shift a to temp_high, temp_low.  */
 982   double_int temp_high, temp_low;
 983   int amount = GET_MODE_FBIT (mode);
 984   if (amount == 2 * HOST_BITS_PER_WIDE_INT)
 985     {
 986        temp_high = a;
 987        temp_low.low = 0;
 988        temp_low.high = 0;
 989     }
 990   else
 991     {
 992       lshift_double (a.low, a.high,
 993                      amount,
 994                      2 * HOST_BITS_PER_WIDE_INT,
 995                      &temp_low.low, &temp_low.high, 0);
 996
 997       /* Logical shift right to temp_high.  */
 998       lshift_double (a.low, a.high,
 999                      amount - 2 * HOST_BITS_PER_WIDE_INT,
1000                      2 * HOST_BITS_PER_WIDE_INT,
1001                      &temp_high.low, &temp_high.high, 0);
1002     }
1003   if (!unsigned_p && a.high < 0) /* Signed-extend temp_high.  */
1004     temp_high = double_int_ext (temp_high, amount, 0);
1005
1006   f->mode = mode;
1007   f->data = temp_low;
1008
1009   if (unsigned_p == UNSIGNED_FIXED_POINT_MODE_P (f->mode))
1010     overflow_p = fixed_saturate2 (f->mode, temp_high, temp_low, &f->data,
1011                                   sat_p);
1012   else
1013     {
1014       /* Take care of the cases when converting between signed and unsigned.  */
1015       if (!unsigned_p)
1016         {
1017           /* Signed -> Unsigned.  */
1018           if (a.high < 0)
1019             {
1020               if (sat_p)
1021                 {
1022                   f->data.low = 0;  /* Set to zero.  */
1023                   f->data.high = 0;  /* Set to zero.  */
1024                 }
1025               else
1026                 overflow_p = true;
1027             }
1028           else
1029             overflow_p = fixed_saturate2 (f->mode, temp_high, temp_low,
1030                                           &f->data, sat_p);
1031         }
1032       else
1033         {
1034           /* Unsigned -> Signed.  */
1035           if (temp_high.high < 0)
1036             {
1037               if (sat_p)
1038                 {
1039                   /* Set to maximum.  */
1040                   f->data.low = -1;  /* Set to all ones.  */
1041                   f->data.high = -1;  /* Set to all ones.  */
1042                   f->data = double_int_ext (f->data,
1043                                             GET_MODE_FBIT (f->mode)
1044                                             + GET_MODE_IBIT (f->mode),
1045                                             1); /* Clear the sign.  */
1046                 }
1047               else
1048                 overflow_p = true;
1049             }
1050           else
1051             overflow_p = fixed_saturate2 (f->mode, temp_high, temp_low,
1052                                           &f->data, sat_p);
1053         }
1054     }
1055   f->data = double_int_ext (f->data,
1056                             SIGNED_FIXED_POINT_MODE_P (f->mode)
1057                             + GET_MODE_FBIT (f->mode)
1058                             + GET_MODE_IBIT (f->mode),
1059                             UNSIGNED_FIXED_POINT_MODE_P (f->mode));
1060   return overflow_p;
1061 }
1062
1063 /* Convert to a new fixed-point mode from a real.
1064    If SAT_P, saturate the result to the max or the min.
1065    Return true, if !SAT_P and overflow.  */
1066
1067 bool
1068 fixed_convert_from_real (FIXED_VALUE_TYPE *f, enum machine_mode mode,
1069                          const REAL_VALUE_TYPE *a, bool sat_p)
1070 {
1071   bool overflow_p = false;
1072   REAL_VALUE_TYPE real_value, fixed_value, base_value;
1073   bool unsigned_p = UNSIGNED_FIXED_POINT_MODE_P (mode);
1074   int i_f_bits = GET_MODE_IBIT (mode) + GET_MODE_FBIT (mode);
1075   unsigned int fbit = GET_MODE_FBIT (mode);
1076   enum fixed_value_range_code temp;
1077
1078   real_value = *a;
1079   f->mode = mode;
1080   real_2expN (&base_value, fbit, mode);
1081   real_arithmetic (&fixed_value, MULT_EXPR, &real_value, &base_value);
1082   real_to_integer2 ((HOST_WIDE_INT *)&f->data.low, &f->data.high, &fixed_value);
1083   temp = check_real_for_fixed_mode (&real_value, mode);
1084   if (temp == FIXED_UNDERFLOW) /* Minimum.  */
1085     {
1086       if (sat_p)
1087         {
1088           if (unsigned_p)
1089             {
1090               f->data.low = 0;
1091               f->data.high = 0;
1092             }
1093           else
1094             {
1095               f->data.low = 1;
1096               f->data.high = 0;
1097               lshift_double (f->data.low, f->data.high, i_f_bits,
1098                              2 * HOST_BITS_PER_WIDE_INT,
1099                              &f->data.low, &f->data.high, 1);
1100               f->data = double_int_ext (f->data, 1 + i_f_bits, 0);
1101             }
1102         }
1103       else
1104         overflow_p = true;
1105     }
1106   else if (temp == FIXED_GT_MAX_EPS || temp == FIXED_MAX_EPS) /* Maximum.  */
1107     {
1108       if (sat_p)
1109         {
1110           f->data.low = -1;
1111           f->data.high = -1;
1112           f->data = double_int_ext (f->data, i_f_bits, 1);
1113         }
1114       else
1115         overflow_p = true;
1116     }
1117   f->data = double_int_ext (f->data, (!unsigned_p) + i_f_bits, unsigned_p);
1118   return overflow_p;
1119 }
1120
1121 /* Convert to a new real mode from a fixed-point.  */
1122
1123 void
1124 real_convert_from_fixed (REAL_VALUE_TYPE *r, enum machine_mode mode,
1125                          const FIXED_VALUE_TYPE *f)
1126 {
1127   REAL_VALUE_TYPE base_value, fixed_value, real_value;
1128
1129   real_2expN (&base_value, GET_MODE_FBIT (f->mode), f->mode);
1130   real_from_integer (&fixed_value, VOIDmode, f->data.low, f->data.high,
1131                      UNSIGNED_FIXED_POINT_MODE_P (f->mode));
1132   real_arithmetic (&real_value, RDIV_EXPR, &fixed_value, &base_value);
1133   real_convert (r, mode, &real_value);
1134 }
1135
1136 /* Determine whether a fixed-point value F is negative.  */
1137
1138 bool
1139 fixed_isneg (const FIXED_VALUE_TYPE *f)
1140 {
1141   if (SIGNED_FIXED_POINT_MODE_P (f->mode))
1142     {
1143       int i_f_bits = GET_MODE_IBIT (f->mode) + GET_MODE_FBIT (f->mode);
1144       int sign_bit = get_fixed_sign_bit (f->data, i_f_bits);
1145       if (sign_bit == 1)
1146         return true;
1147     }
1148
1149   return false;
1150 }