typos fixed
[libfirm] / ir / tv / tv.c
index 38834a2..b665375 100644 (file)
@@ -1,5 +1,5 @@
 /*
- * Copyright (C) 1995-2007 University of Karlsruhe.  All right reserved.
+ * Copyright (C) 1995-2008 University of Karlsruhe.  All right reserved.
  *
  * This file is part of libFirm.
  *
@@ -158,6 +158,29 @@ static int hash_val(const void *value, unsigned int length) {
        return hash;
 }
 
+static int cmp_tv(const void *p1, const void *p2, size_t n) {
+       const tarval *tv1 = p1;
+       const tarval *tv2 = p2;
+       (void) n;
+
+       assert(tv1->kind == k_tarval);
+       assert(tv2->kind == k_tarval);
+       if(tv1->mode < tv2->mode)
+               return -1;
+       if(tv1->mode > tv2->mode)
+               return 1;
+       if(tv1->length < tv2->length)
+               return -1;
+       if(tv1->length > tv2->length)
+               return 1;
+       if(tv1->value < tv2->value)
+               return -1;
+       if(tv1->value > tv2->value)
+               return 1;
+
+       return 0;
+}
+
 /** finds tarval with value/mode or creates new tarval */
 static tarval *get_tarval(const void *value, int length, ir_mode *mode) {
        tarval tv;
@@ -168,7 +191,12 @@ static tarval *get_tarval(const void *value, int length, ir_mode *mode) {
        if (length > 0) {
                /* if there already is such a value, it is returned, else value
                 * is copied into the set */
-               tv.value = INSERT_VALUE(value, length);
+               char *temp = alloca(length);
+               memcpy(temp, value, length);
+               if (get_mode_arithmetic(mode) == irma_twos_complement) {
+                       sign_extend(temp, mode);
+               }
+               tv.value = INSERT_VALUE(temp, length);
        } else {
                tv.value = value;
        }
@@ -188,8 +216,8 @@ static tarval *get_tarval_overflow(const void *value, int length, ir_mode *mode)
        case irms_reference:
                /* addresses always wrap around */
                temp = alloca(sc_get_buffer_length());
-               sc_val_from_ulong(-1, temp);
-               sc_and(temp, value, temp);
+               memcpy(temp, value, sc_get_buffer_length());
+               sc_truncate(get_mode_size_bits(mode), temp);
                /* the sc_ module expects that all bits are set ... */
                sign_extend(temp, mode);
                return get_tarval(temp, length, mode);
@@ -201,8 +229,8 @@ static tarval *get_tarval_overflow(const void *value, int length, ir_mode *mode)
                                return get_mode_max(mode);
                        case TV_OVERFLOW_WRAP:
                                temp = alloca(sc_get_buffer_length());
-                               sc_val_from_ulong(-1, temp);
-                               sc_and(temp, value, temp);
+                               memcpy(temp, value, sc_get_buffer_length());
+                               sc_truncate(get_mode_size_bits(mode), temp);
                                /* the sc_ module expects that all bits are set ... */
                                sign_extend(temp, mode);
                                return get_tarval(temp, length, mode);
@@ -218,8 +246,8 @@ static tarval *get_tarval_overflow(const void *value, int length, ir_mode *mode)
                                return get_mode_min(mode);
                        case TV_OVERFLOW_WRAP: {
                                char *temp = alloca(sc_get_buffer_length());
-                               sc_val_from_ulong(-1, temp);
-                               sc_and(temp, value, temp);
+                               memcpy(temp, value, sc_get_buffer_length());
+                               sc_truncate(get_mode_size_bits(mode), temp);
                                return get_tarval(temp, length, mode);
                        }
                        case TV_OVERFLOW_BAD:
@@ -232,10 +260,12 @@ static tarval *get_tarval_overflow(const void *value, int length, ir_mode *mode)
 
        case irms_float_number:
                if (SWITCH_NOINFINITY && fc_is_inf(value)) {
-                       return fc_is_negative(value)?get_mode_min(mode):get_mode_max(mode);
+                       /* clip infinity to maximum value */
+                       return fc_is_negative(value) ? get_mode_min(mode) : get_mode_max(mode);
                }
 
                if (SWITCH_NODENORMALS && fc_is_subnormal(value)) {
+                       /* clip denormals to zero */
                        return get_mode_null(mode);
                }
                break;
@@ -293,8 +323,11 @@ tarval *new_tarval_from_str(const char *str, size_t len, ir_mode *mode)
                        fc_val_from_str(str, len, 11, 52, NULL);
                        break;
                case 80:
+               case 96:
                        fc_val_from_str(str, len, 15, 64, NULL);
                        break;
+               default:
+                       panic("Unsupported mode in new_tarval_from_str()");
                }
                return get_tarval(fc_get_buffer(), fc_get_buffer_length(), mode);
 
@@ -324,7 +357,6 @@ tarval *new_tarval_from_long(long l, ir_mode *mode) {
                /* same as integer modes */
        case irms_int_number:
                sc_val_from_long(l, NULL);
-               sign_extend(sc_get_buffer(), mode);
                return get_tarval(sc_get_buffer(), sc_get_buffer_length(), mode);
 
        case irms_float_number:
@@ -370,8 +402,11 @@ tarval *new_tarval_from_double(long double d, ir_mode *mode) {
                fc_val_from_ieee754(d, 11, 52, NULL);
                break;
        case 80:
+       case 96:
                fc_val_from_ieee754(d, 15, 64, NULL);
                break;
+       default:
+               panic("Unsupported mode in new_tarval_from_double()");
        }
        return get_tarval(fc_get_buffer(), fc_get_buffer_length(), mode);
 }
@@ -452,8 +487,11 @@ tarval *get_tarval_max(ir_mode *mode) {
                        fc_get_max(11, 52, NULL);
                        break;
                case 80:
+               case 96:
                        fc_get_max(15, 64, NULL);
                        break;
+               default:
+                       panic("Unsupported mode in get_tarval_max()");
                }
                return get_tarval(fc_get_buffer(), fc_get_buffer_length(), mode);
 
@@ -492,8 +530,11 @@ tarval *get_tarval_min(ir_mode *mode) {
                        fc_get_min(11, 52, NULL);
                        break;
                case 80:
+               case 96:
                        fc_get_min(15, 64, NULL);
                        break;
+               default:
+                       panic("Unsupported mode in get_tarval_min()");
                }
                return get_tarval(fc_get_buffer(), fc_get_buffer_length(), mode);
 
@@ -520,13 +561,13 @@ tarval *get_tarval_null(ir_mode *mode) {
        case irms_control_flow:
        case irms_memory:
        case irms_auxiliary:
-       case irms_internal_boolean:
                assert(0);
                break;
 
        case irms_float_number:
                return new_tarval_from_double(0.0, mode);
 
+       case irms_internal_boolean:
        case irms_int_number:
                return new_tarval_from_long(0l,  mode);
 
@@ -541,6 +582,7 @@ tarval *get_tarval_one(ir_mode *mode) {
 
        if (get_mode_n_vector_elems(mode) > 1) {
                /* vector arithmetic not implemented yet */
+               assert(0);
                return tarval_bad;
        }
 
@@ -548,26 +590,28 @@ tarval *get_tarval_one(ir_mode *mode) {
        case irms_control_flow:
        case irms_memory:
        case irms_auxiliary:
-       case irms_internal_boolean:
-       case irms_reference:
                assert(0);
                break;
 
+       case irms_internal_boolean:
+               return tarval_b_true;
+
        case irms_float_number:
                return new_tarval_from_double(1.0, mode);
 
+       case irms_reference:
        case irms_int_number:
                return new_tarval_from_long(1l, mode);
-               break;
        }
        return tarval_bad;
 }
 
-tarval *get_tarval_minus_one(ir_mode *mode) {
+tarval *get_tarval_all_one(ir_mode *mode) {
        assert(mode);
 
        if (get_mode_n_vector_elems(mode) > 1) {
                /* vector arithmetic not implemented yet */
+               assert(0);
                return tarval_bad;
        }
 
@@ -575,11 +619,40 @@ tarval *get_tarval_minus_one(ir_mode *mode) {
        case irms_control_flow:
        case irms_memory:
        case irms_auxiliary:
+               assert(0);
+               return tarval_bad;
+
+       case irms_int_number:
        case irms_internal_boolean:
        case irms_reference:
+               return tarval_not(get_mode_null(mode));
+
+
+       case irms_float_number:
+               return new_tarval_from_double(1.0, mode);
+       }
+       return tarval_bad;
+}
+
+tarval *get_tarval_minus_one(ir_mode *mode) {
+       assert(mode);
+
+       if (get_mode_n_vector_elems(mode) > 1) {
+               /* vector arithmetic not implemented yet */
+               return tarval_bad;
+       }
+
+       switch(get_mode_sort(mode)) {
+       case irms_control_flow:
+       case irms_memory:
+       case irms_auxiliary:
+       case irms_internal_boolean:
                assert(0);
                break;
 
+       case irms_reference:
+               return tarval_bad;
+
        case irms_float_number:
                return mode_is_signed(mode) ? new_tarval_from_double(-1.0, mode) : tarval_bad;
 
@@ -606,8 +679,11 @@ tarval *get_tarval_nan(ir_mode *mode) {
                        fc_get_qnan(11, 52, NULL);
                        break;
                case 80:
+               case 96:
                        fc_get_qnan(15, 64, NULL);
                        break;
+               default:
+                       panic("Unsupported mode in get_tarval_nan()");
                }
                return get_tarval(fc_get_buffer(), fc_get_buffer_length(), mode);
        } else {
@@ -633,8 +709,11 @@ tarval *get_tarval_plus_inf(ir_mode *mode) {
                        fc_get_plusinf(11, 52, NULL);
                        break;
                case 80:
+               case 96:
                        fc_get_plusinf(15, 64, NULL);
                        break;
+               default:
+                       panic("Unsupported mode in get_tarval_plus_inf()");
                }
                return get_tarval(fc_get_buffer(), fc_get_buffer_length(), mode);
        } else {
@@ -660,8 +739,11 @@ tarval *get_tarval_minus_inf(ir_mode *mode) {
                        fc_get_minusinf(11, 52, NULL);
                        break;
                case 80:
+               case 96:
                        fc_get_minusinf(15, 64, NULL);
                        break;
+               default:
+                       panic("Unsupported mode in get_tarval_minus_inf()");
                }
                return get_tarval(fc_get_buffer(), fc_get_buffer_length(), mode);
        } else {
@@ -693,7 +775,7 @@ int tarval_is_negative(tarval *a) {
                        return sc_comp(a->value, get_mode_null(a->mode)->value) == -1 ? 1 : 0;
 
        case irms_float_number:
-               return fc_comp(a->value, get_mode_null(a->mode)->value) == -1 ? 1 : 0;
+               return fc_is_negative(a->value);
 
        default:
                assert(0 && "not implemented");
@@ -705,18 +787,33 @@ int tarval_is_negative(tarval *a) {
  * test if null, 1 means 'yes'
  */
 int tarval_is_null(tarval *a) {
-       ir_mode *m = get_tarval_mode(a);
-
-       return a == get_tarval_null(m);
+       return
+               a != tarval_bad &&
+               a == get_mode_null(get_tarval_mode(a));
 }
 
 /*
  * test if one, 1 means 'yes'
  */
 int tarval_is_one(tarval *a) {
-       ir_mode *m = get_tarval_mode(a);
+       return
+               a != tarval_bad &&
+               a == get_mode_one(get_tarval_mode(a));
+}
+
+int tarval_is_all_one(tarval *tv) {
+       return
+               tv != tarval_bad &&
+               tv == get_mode_all_one(get_tarval_mode(tv));
+}
 
-       return a == get_tarval_one(m);
+/*
+ * test if one, 1 means 'yes'
+ */
+int tarval_is_minus_one(tarval *a) {
+       return
+               a != tarval_bad &&
+               a == get_mode_minus_one(get_tarval_mode(a));
 }
 
 /*
@@ -752,8 +849,9 @@ pn_Cmp tarval_cmp(tarval *a, tarval *b) {
                return pn_Cmp_False;
 
        case irms_float_number:
-               if (no_float)
-                       return pn_Cmp_False;
+               /* it should be safe to enable this even if other arithmetic is disabled */
+               /*if (no_float)
+                       return pn_Cmp_False;*/
                /*
                 * BEWARE: we cannot compare a == b here, because
                 * a NaN is always Unordered to any other value, even to itself!
@@ -782,13 +880,14 @@ pn_Cmp tarval_cmp(tarval *a, tarval *b) {
 /*
  * convert to other mode
  */
-tarval *tarval_convert_to(tarval *src, ir_mode *m) {
+tarval *tarval_convert_to(tarval *src, ir_mode *dst_mode) {
        char *buffer;
+       fp_value *res;
 
        assert(src);
-       assert(m);
+       assert(dst_mode);
 
-       if (src->mode == m) return src;
+       if (src->mode == dst_mode) return src;
 
        if (get_mode_n_vector_elems(src->mode) > 1) {
                /* vector arithmetic not implemented yet */
@@ -803,9 +902,9 @@ tarval *tarval_convert_to(tarval *src, ir_mode *m) {
 
                /* cast float to something */
        case irms_float_number:
-               switch (get_mode_sort(m)) {
+               switch (get_mode_sort(dst_mode)) {
                case irms_float_number:
-                       switch (get_mode_size_bits(m)) {
+                       switch (get_mode_size_bits(dst_mode)) {
                        case 32:
                                fc_cast(src->value, 8, 23, NULL);
                                break;
@@ -813,30 +912,30 @@ tarval *tarval_convert_to(tarval *src, ir_mode *m) {
                                fc_cast(src->value, 11, 52, NULL);
                                break;
                        case 80:
+                       case 96:
                                fc_cast(src->value, 15, 64, NULL);
                                break;
                        default:
-                               break;
+                               panic("Unsupported mode in tarval_convert_to()");
                        }
-                       return get_tarval(fc_get_buffer(), fc_get_buffer_length(), m);
+                       return get_tarval(fc_get_buffer(), fc_get_buffer_length(), dst_mode);
 
                case irms_int_number:
                        switch (GET_FLOAT_TO_INT_MODE()) {
                        case TRUNCATE:
-                               fc_int(src->value, NULL);
+                               res = fc_int(src->value, NULL);
                                break;
                        case ROUND:
-                               fc_rnd(src->value, NULL);
+                               res = fc_rnd(src->value, NULL);
                                break;
                        default:
-                               assert(0);
+                               panic("Unsupported float to int conversion mode in tarval_convert_to()");
                                break;
                        }
-                       /* FIXME: floating point unit can't produce a value in integer
-                        * representation
-                        * an intermediate representation is needed here first. */
-                       /*  return get_tarval(); */
-                       return tarval_bad;
+                       buffer = alloca(sc_get_buffer_length());
+                       if (! fc_flt2int(res, buffer, dst_mode))
+                               return tarval_bad;
+                       return get_tarval(buffer, sc_get_buffer_length(), dst_mode);
 
                default:
                        /* the rest can't be converted */
@@ -846,14 +945,14 @@ tarval *tarval_convert_to(tarval *src, ir_mode *m) {
 
        /* cast int/characters to something */
        case irms_int_number:
-               switch (get_mode_sort(m)) {
+               switch (get_mode_sort(dst_mode)) {
 
                case irms_reference:
                case irms_int_number:
                        buffer = alloca(sc_get_buffer_length());
                        memcpy(buffer, src->value, sc_get_buffer_length());
-                       sign_extend(buffer, m);
-                       return get_tarval_overflow(buffer, src->length, m);
+                       sign_extend(buffer, dst_mode);
+                       return get_tarval_overflow(buffer, src->length, dst_mode);
 
                case irms_internal_boolean:
                        /* XXX C semantics */
@@ -870,7 +969,7 @@ tarval *tarval_convert_to(tarval *src, ir_mode *m) {
                        snprintf(buffer, 100, "%s",
                                sc_print(src->value, get_mode_size_bits(src->mode), SC_DEC, mode_is_signed(src->mode)));
                        buffer[100 - 1] = '\0';
-                       switch (get_mode_size_bits(m)) {
+                       switch (get_mode_size_bits(dst_mode)) {
                        case 32:
                                fc_val_from_str(buffer, 0, 8, 23, NULL);
                                break;
@@ -878,10 +977,13 @@ tarval *tarval_convert_to(tarval *src, ir_mode *m) {
                                fc_val_from_str(buffer, 0, 11, 52, NULL);
                                break;
                        case 80:
+                       case 96:
                                fc_val_from_str(buffer, 0, 15, 64, NULL);
                                break;
+                       default:
+                               panic("Unsupported mode in tarval_convert_to()");
                        }
-                       return get_tarval(fc_get_buffer(), fc_get_buffer_length(), m);
+                       return get_tarval(fc_get_buffer(), fc_get_buffer_length(), dst_mode);
 
                default:
                        break;
@@ -890,16 +992,16 @@ tarval *tarval_convert_to(tarval *src, ir_mode *m) {
 
        case irms_internal_boolean:
                /* beware: this is C semantic for the INTERNAL boolean mode */
-               if (get_mode_sort(m) == irms_int_number)
-                       return src == tarval_b_true ? get_mode_one(m) : get_mode_null(m);
+               if (get_mode_sort(dst_mode) == irms_int_number)
+                       return src == tarval_b_true ? get_mode_one(dst_mode) : get_mode_null(dst_mode);
                break;
 
        case irms_reference:
-               if (get_mode_sort(m) == irms_int_number) {
+               if (get_mode_sort(dst_mode) == irms_int_number) {
                        buffer = alloca(sc_get_buffer_length());
                        memcpy(buffer, src->value, sc_get_buffer_length());
                        sign_extend(buffer, src->mode);
-                       return get_tarval_overflow(buffer, src->length, m);
+                       return get_tarval_overflow(buffer, src->length, dst_mode);
                }
                break;
        }
@@ -918,6 +1020,7 @@ tarval *tarval_not(tarval *a) {
        /* works for vector mode without changes */
 
        switch (get_mode_sort(a->mode)) {
+       case irms_reference:
        case irms_int_number:
                buffer = alloca(sc_get_buffer_length());
                sc_not(a->value, buffer);
@@ -959,8 +1062,9 @@ tarval *tarval_neg(tarval *a) {
                return get_tarval_overflow(buffer, a->length, a->mode);
 
        case irms_float_number:
-               if(no_float)
-                       return tarval_bad;
+               /* it should be safe to enable this even if other arithmetic is disabled */
+               /*if (no_float)
+                       return tarval_bad;*/
 
                fc_neg(a->value, NULL);
                return get_tarval_overflow(fc_get_buffer(), fc_get_buffer_length(), a->mode);
@@ -993,7 +1097,7 @@ tarval *tarval_add(tarval *a, tarval *b) {
                return get_tarval_overflow(buffer, a->length, a->mode);
 
        case irms_float_number:
-               if(no_float)
+               if (no_float)
                        return tarval_bad;
 
                fc_add(a->value, b->value, NULL);
@@ -1026,7 +1130,7 @@ tarval *tarval_sub(tarval *a, tarval *b) {
                return get_tarval_overflow(buffer, a->length, a->mode);
 
        case irms_float_number:
-               if(no_float)
+               if (no_float)
                        return tarval_bad;
 
                fc_sub(a->value, b->value, NULL);
@@ -1060,7 +1164,7 @@ tarval *tarval_mul(tarval *a, tarval *b) {
                return get_tarval_overflow(buffer, a->length, a->mode);
 
        case irms_float_number:
-               if(no_float)
+               if (no_float)
                        return tarval_bad;
 
                fc_mul(a->value, b->value, NULL);
@@ -1079,7 +1183,7 @@ tarval *tarval_quo(tarval *a, tarval *b) {
        assert(b);
        assert((a->mode == b->mode) && mode_is_float(a->mode));
 
-       if(no_float)
+       if (no_float)
                return tarval_bad;
 
        if (get_mode_n_vector_elems(a->mode) > 1) {
@@ -1184,8 +1288,9 @@ tarval *tarval_abs(tarval *a) {
                return a;
 
        case irms_float_number:
-               if(no_float)
-                       return tarval_bad;
+               /* it should be safe to enable this even if other arithmetic is disabled */
+               /*if (no_float)
+                       return tarval_bad;*/
 
                if (fc_comp(a->value, get_mode_null(a->mode)->value) == -1) {
                        fc_neg(a->value, NULL);
@@ -1383,6 +1488,7 @@ tarval *tarval_rot(tarval *a, tarval *b) {
  * carry flag of the last operation
  */
 int tarval_carry(void) {
+       panic("tarval_carry() requetsed: not implemented on all operations");
        return sc_had_carry();
 }
 
@@ -1535,29 +1641,6 @@ const tarval_mode_info *get_tarval_mode_output_option(ir_mode *mode) {
        return mode->tv_priv;
 }
 
-/*
- * Identifying tarvals values for algebraic simplifications.
- *
- * Returns:
- *   - TV_CLASSIFY_NULL    for additive neutral or the NULL tarval for reference modes,
- *   - TV_CLASSIFY_ONE     for multiplicative neutral,
- *   - TV_CLASSIFY_ALL_ONE for bitwise-and neutral
- *   - TV_CLASSIFY_OTHER   else
- */
-tarval_classification_t classify_tarval(tarval *tv) {
-       if (!tv || tv == tarval_bad) return TV_CLASSIFY_OTHER;
-
-       if (tv == get_mode_null(tv->mode))
-               return TV_CLASSIFY_NULL;
-       else if (tv == get_mode_one(tv->mode))
-               return TV_CLASSIFY_ONE;
-       else if ((get_mode_sort(tv->mode) == irms_int_number)
-               && (tv == get_mode_minus_one(tv->mode)))
-               return TV_CLASSIFY_ALL_ONE;
-
-       return TV_CLASSIFY_OTHER;
-}
-
 /*
  * Returns non-zero if a given (integer) tarval has only one single bit
  * set.
@@ -1599,6 +1682,30 @@ int tarval_ieee754_get_exponent(tarval *tv) {
        return fc_get_exponent(tv->value);
 }
 
+/*
+ * Check if the tarval can be converted to the given mode without
+ * precision loss.
+ */
+int tarval_ieee754_can_conv_lossless(tarval *tv, ir_mode *mode) {
+       char exp_size, mant_size;
+       switch (get_mode_size_bits(mode)) {
+       case 32:
+               exp_size = 8; mant_size = 23;
+               break;
+       case 64:
+               exp_size = 11; mant_size = 52;
+               break;
+       case 80:
+       case 96:
+               exp_size = 15; mant_size = 64;
+               break;
+       default:
+               panic("Unsupported mode in tarval_ieee754_can_conv_lossless()");
+               return 0;
+       }
+       return fc_can_lossless_conv_to(tv->value, exp_size, mant_size);
+}
+
 /* Set the immediate precision for IEEE-754 results. */
 unsigned tarval_ieee754_set_immediate_precision(unsigned bits) {
        return fc_set_immediate_precision(bits);
@@ -1609,6 +1716,34 @@ unsigned tarval_ieee754_get_exact(void) {
        return fc_is_exact();
 }
 
+/* check if its the a floating point NaN */
+int tarval_is_NaN(tarval *tv) {
+       if (! mode_is_float(tv->mode))
+               return 0;
+       return fc_is_nan(tv->value);
+}
+
+/* check if its the a floating point +inf */
+int tarval_is_plus_inf(tarval *tv) {
+       if (! mode_is_float(tv->mode))
+               return 0;
+       return fc_is_inf(tv->value) && !fc_is_negative(tv->value);
+}
+
+/* check if its the a floating point -inf */
+int tarval_is_minus_inf(tarval *tv) {
+       if (! mode_is_float(tv->mode))
+               return 0;
+       return fc_is_inf(tv->value) && fc_is_negative(tv->value);
+}
+
+/* check if the tarval represents a finite value */
+int tarval_is_finite(tarval *tv) {
+       if (mode_is_float(tv->mode))
+               return !fc_is_nan(tv->value) && !fc_is_inf(tv->value);
+       return 1;
+}
+
 /*
  * Sets the overflow mode for integer operations.
  */
@@ -1646,7 +1781,7 @@ void init_tarval_1(long null_value) {
 
        /* initialize the sets holding the tarvals with a comparison function and
         * an initial size, which is the expected number of constants */
-       tarvals = new_set(memcmp, N_CONSTANTS);
+       tarvals = new_set(cmp_tv, N_CONSTANTS);
        values  = new_set(memcmp, N_CONSTANTS);
        /* init strcalc with precision of 68 to support floating point values with 64
         * bit mantissa (needs extra bits for rounding and overflow) */