X-Git-Url: http://nsz.repo.hu/git/?a=blobdiff_plain;f=ir%2Ftv%2Ftv.c;h=a26ed24c7979d531bc6c6291cd78322f51c49a5a;hb=844d1d6841247022150e8391b5aa65391bf41b9b;hp=b8858d9cfedd72d6a0cca1234d5e2b4abd7f777f;hpb=69eb2dc7a32de8cefa7a18d2e257c8ff59e782f5;p=libfirm

diff --git a/ir/tv/tv.c b/ir/tv/tv.c
index b8858d9cf..a26ed24c7 100644
--- a/ir/tv/tv.c
+++ b/ir/tv/tv.c
@@ -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.
  *
@@ -94,16 +94,29 @@ enum float_to_int_mode {
 
 #define fail_verify(a) _fail_verify((a), __FILE__, __LINE__)
 
-/****************************************************************************
- *   private variables
- ****************************************************************************/
-static struct set *tarvals = NULL;   /* container for tarval structs */
-static struct set *values = NULL;    /* container for values */
+/** A set containing all existing tarvals. */
+static struct set *tarvals = NULL;
+/** A set containing all existing values. */
+static struct set *values = NULL;
+
+/** The integer overflow mode. */
 static tarval_int_overflow_mode_t int_overflow_mode = TV_OVERFLOW_WRAP;
 
 /** if this is set non-zero, the constant folding for floating point is OFF */
 static int no_float = 0;
 
+/** IEEE-754r half precision */
+static const ieee_descriptor_t half_desc     = {  5,  10, 0, NORMAL };
+/** IEEE-754 single precision */
+static const ieee_descriptor_t single_desc   = {  8,  23, 0, NORMAL };
+/** IEEE-754 double precision */
+static const ieee_descriptor_t double_desc   = { 11,  52, 0, NORMAL };
+/** Intel x87 extended precision */
+static const ieee_descriptor_t extended_desc = { 15,  63, 1, NORMAL };
+
+/** IEEE-754r quad precision */
+static const ieee_descriptor_t quad_desc     = { 15, 112, 0, NORMAL };
+
 /****************************************************************************
  *   private functions
  ****************************************************************************/
@@ -114,10 +127,9 @@ static void _fail_verify(tarval *tv, const char* file, int line)
 {
 	/* print a memory image of the tarval and throw an assertion */
 	if (tv)
-		printf("%s:%d: Invalid tarval:\n  mode: %s\n value: [%p]\n", file, line, get_mode_name(tv->mode), tv->value);
+		panic("%s:%d: Invalid tarval: mode: %F\n value: [%p]", file, line, tv->mode, tv->value);
 	else
-		printf("%s:%d: Invalid tarval (null)", file, line);
-	assert(0);
+		panic("%s:%d: Invalid tarval (null)", file, line);
 }
 #ifdef __GNUC__
 INLINE static void tarval_verify(tarval *tv) __attribute__ ((unused));
@@ -134,8 +146,6 @@ INLINE static void tarval_verify(tarval *tv)
 
 	if (!FIND_TARVAL(tv)) fail_verify(tv);
 	if (tv->length > 0 && !FIND_VALUE(tv->value, tv->length)) fail_verify(tv);
-
-	return;
 }
 #endif /* NDEBUG */
 
@@ -158,6 +168,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;
@@ -256,12 +289,31 @@ static tarval *get_tarval_overflow(const void *value, int length, ir_mode *mode)
 /*
  *   public variables declared in tv.h
  */
-static tarval reserved_tv[4];
+static tarval reserved_tv[6];
 
-tarval *tarval_bad       = &reserved_tv[0];
-tarval *tarval_undefined = &reserved_tv[1];
-tarval *tarval_b_false   = &reserved_tv[2];
-tarval *tarval_b_true    = &reserved_tv[3];
+tarval *tarval_b_false     = &reserved_tv[0];
+tarval *tarval_b_true      = &reserved_tv[1];
+tarval *tarval_bad         = &reserved_tv[2];
+tarval *tarval_undefined   = &reserved_tv[3];
+tarval *tarval_reachable   = &reserved_tv[4];
+tarval *tarval_unreachable = &reserved_tv[5];
+
+/**
+ * get the float descriptor for given mode.
+ */
+static const ieee_descriptor_t *get_descriptor(const ir_mode *mode) {
+	switch (get_mode_size_bits(mode)) {
+	case 16:  return &half_desc;
+	case 32:  return &single_desc;
+	case 64:  return &double_desc;
+	case 80:
+	case 96:  return &extended_desc;
+	case 128: return &quad_desc;
+	default:
+		panic("Unsupported mode in get_descriptor()");
+		return NULL;
+	}
+}
 
 /*
  *   public functions declared in tv.h
@@ -272,6 +324,8 @@ tarval *tarval_b_true    = &reserved_tv[3];
  */
 tarval *new_tarval_from_str(const char *str, size_t len, ir_mode *mode)
 {
+	const ieee_descriptor_t *desc;
+
 	assert(str);
 	assert(len);
 	assert(mode);
@@ -280,29 +334,21 @@ tarval *new_tarval_from_str(const char *str, size_t len, ir_mode *mode)
 	case irms_control_flow:
 	case irms_memory:
 	case irms_auxiliary:
-		assert(0);
-		break;
+		panic("Unsupported tarval creation with mode %F", mode);
 
 	case irms_internal_boolean:
 		/* match [tT][rR][uU][eE]|[fF][aA][lL][sS][eE] */
-		if (strcasecmp(str, "true")) return tarval_b_true;
-		else if (strcasecmp(str, "false")) return tarval_b_true;
+		if (strcasecmp(str, "true"))
+			return tarval_b_true;
+		else if (strcasecmp(str, "false"))
+			return tarval_b_true;
 		else
 			/* XXX This is C semantics */
 			return atoi(str) ? tarval_b_true : tarval_b_false;
 
 	case irms_float_number:
-		switch (get_mode_size_bits(mode)) {
-		case 32:
-			fc_val_from_str(str, len, 8, 23, NULL);
-			break;
-		case 64:
-			fc_val_from_str(str, len, 11, 52, NULL);
-			break;
-		case 80:
-			fc_val_from_str(str, len, 15, 64, NULL);
-			break;
-		}
+		desc = get_descriptor(mode);
+		fc_val_from_str(str, len, desc, NULL);
 		return get_tarval(fc_get_buffer(), fc_get_buffer_length(), mode);
 
 	case irms_reference:
@@ -311,9 +357,7 @@ tarval *new_tarval_from_str(const char *str, size_t len, ir_mode *mode)
 		sc_val_from_str(str, len, NULL, mode);
 		return get_tarval(sc_get_buffer(), sc_get_buffer_length(), mode);
 	}
-
-	assert(0);  /* can't be reached, can it? */
-	return NULL;
+	panic("Unsupported tarval creation with mode %F", mode);
 }
 
 /*
@@ -366,19 +410,11 @@ long get_tarval_long(tarval* tv) {
 }
 
 tarval *new_tarval_from_double(long double d, ir_mode *mode) {
-	assert(mode && (get_mode_sort(mode) == irms_float_number));
+	const ieee_descriptor_t *desc;
 
-	switch (get_mode_size_bits(mode)) {
-	case 32:
-		fc_val_from_ieee754(d, 8, 23, NULL);
-		break;
-	case 64:
-		fc_val_from_ieee754(d, 11, 52, NULL);
-		break;
-	case 80:
-		fc_val_from_ieee754(d, 15, 64, NULL);
-		break;
-	}
+	assert(mode && (get_mode_sort(mode) == irms_float_number));
+	desc = get_descriptor(mode);
+	fc_val_from_ieee754(d, desc, NULL);
 	return get_tarval(fc_get_buffer(), fc_get_buffer_length(), mode);
 }
 
@@ -431,36 +467,35 @@ tarval *(get_tarval_b_true)(void) {
 	return _get_tarval_b_true();
 }
 
+tarval *(get_tarval_reachable)(void) {
+	return _get_tarval_reachable();
+}
+
+tarval *(get_tarval_unreachable)(void) {
+	return _get_tarval_unreachable();
+}
+
 tarval *get_tarval_max(ir_mode *mode) {
-	assert(mode);
+	const ieee_descriptor_t *desc;
 
+	assert(mode);
 	if (get_mode_n_vector_elems(mode) > 1) {
 		/* vector arithmetic not implemented yet */
 		return tarval_bad;
 	}
 
-	switch(get_mode_sort(mode)) {
+	switch (get_mode_sort(mode)) {
 	case irms_control_flow:
 	case irms_memory:
 	case irms_auxiliary:
-		assert(0);
-		break;
+		panic("mode %F does not support maximum value", mode);
 
 	case irms_internal_boolean:
 		return tarval_b_true;
 
 	case irms_float_number:
-		switch(get_mode_size_bits(mode)) {
-		case 32:
-			fc_get_max(8, 23, NULL);
-			break;
-		case 64:
-			fc_get_max(11, 52, NULL);
-			break;
-		case 80:
-			fc_get_max(15, 64, NULL);
-			break;
-		}
+		desc = get_descriptor(mode);
+		fc_get_max(desc, NULL);
 		return get_tarval(fc_get_buffer(), fc_get_buffer_length(), mode);
 
 	case irms_reference:
@@ -472,35 +507,26 @@ tarval *get_tarval_max(ir_mode *mode) {
 }
 
 tarval *get_tarval_min(ir_mode *mode) {
-	assert(mode);
+	const ieee_descriptor_t *desc;
 
+	assert(mode);
 	if (get_mode_n_vector_elems(mode) > 1) {
 		/* vector arithmetic not implemented yet */
 		return tarval_bad;
 	}
 
-	switch(get_mode_sort(mode)) {
+	switch (get_mode_sort(mode)) {
 	case irms_control_flow:
 	case irms_memory:
 	case irms_auxiliary:
-		assert(0);
-		break;
+		panic("mode %F does not support minimum value", mode);
 
 	case irms_internal_boolean:
 		return tarval_b_false;
 
 	case irms_float_number:
-		switch(get_mode_size_bits(mode)) {
-		case 32:
-			fc_get_min(8, 23, NULL);
-			break;
-		case 64:
-			fc_get_min(11, 52, NULL);
-			break;
-		case 80:
-			fc_get_min(15, 64, NULL);
-			break;
-		}
+		desc = get_descriptor(mode);
+		fc_get_min(desc, NULL);
 		return get_tarval(fc_get_buffer(), fc_get_buffer_length(), mode);
 
 	case irms_reference:
@@ -522,12 +548,11 @@ tarval *get_tarval_null(ir_mode *mode) {
 		return tarval_bad;
 	}
 
-	switch(get_mode_sort(mode)) {
+	switch (get_mode_sort(mode)) {
 	case irms_control_flow:
 	case irms_memory:
 	case irms_auxiliary:
-		assert(0);
-		break;
+		panic("mode %F does not support null value", mode);
 
 	case irms_float_number:
 		return new_tarval_from_double(0.0, mode);
@@ -545,18 +570,14 @@ tarval *get_tarval_null(ir_mode *mode) {
 tarval *get_tarval_one(ir_mode *mode) {
 	assert(mode);
 
-	if (get_mode_n_vector_elems(mode) > 1) {
-		/* vector arithmetic not implemented yet */
-		assert(0);
-		return tarval_bad;
-	}
+	if (get_mode_n_vector_elems(mode) > 1)
+		panic("vector arithmetic not implemented yet");
 
-	switch(get_mode_sort(mode)) {
+	switch (get_mode_sort(mode)) {
 	case irms_control_flow:
 	case irms_memory:
 	case irms_auxiliary:
-		assert(0);
-		break;
+		panic("mode %F does not support one value", mode);
 
 	case irms_internal_boolean:
 		return tarval_b_true;
@@ -574,18 +595,14 @@ tarval *get_tarval_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;
-	}
+	if (get_mode_n_vector_elems(mode) > 1)
+		panic("vector arithmetic not implemented yet");
 
-	switch(get_mode_sort(mode)) {
+	switch (get_mode_sort(mode)) {
 	case irms_control_flow:
 	case irms_memory:
 	case irms_auxiliary:
-		assert(0);
-		return tarval_bad;
+		panic("mode %F does not support all-one value", mode);
 
 	case irms_int_number:
 	case irms_internal_boolean:
@@ -599,21 +616,26 @@ tarval *get_tarval_all_one(ir_mode *mode) {
 	return tarval_bad;
 }
 
+int tarval_is_constant(tarval *tv) {
+	int num_res = sizeof(reserved_tv) / sizeof(reserved_tv[0]);
+
+	/* reserved tarvals are NOT constants. Note that although
+	   tarval_b_true and tarval_b_false are reserved, they are constants of course. */
+	return (tv < &reserved_tv[2] || tv > &reserved_tv[num_res - 1]);
+}
+
 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;
-	}
+	if (get_mode_n_vector_elems(mode) > 1)
+		panic("vector arithmetic not implemented yet");
 
-	switch(get_mode_sort(mode)) {
+	switch (get_mode_sort(mode)) {
 	case irms_control_flow:
 	case irms_memory:
 	case irms_auxiliary:
 	case irms_internal_boolean:
-		assert(0);
-		break;
+		panic("mode %F does not support minus one value", mode);
 
 	case irms_reference:
 		return tarval_bad;
@@ -628,88 +650,49 @@ tarval *get_tarval_minus_one(ir_mode *mode) {
 }
 
 tarval *get_tarval_nan(ir_mode *mode) {
-	assert(mode);
+	const ieee_descriptor_t *desc;
 
-	if (get_mode_n_vector_elems(mode) > 1) {
-		/* vector arithmetic not implemented yet */
-		return tarval_bad;
-	}
+	assert(mode);
+	if (get_mode_n_vector_elems(mode) > 1)
+		panic("vector arithmetic not implemented yet");
 
 	if (get_mode_sort(mode) == irms_float_number) {
-		switch(get_mode_size_bits(mode)) {
-		case 32:
-			fc_get_qnan(8, 23, NULL);
-			break;
-		case 64:
-			fc_get_qnan(11, 52, NULL);
-			break;
-		case 80:
-			fc_get_qnan(15, 64, NULL);
-			break;
-		}
+		desc = get_descriptor(mode);
+		fc_get_qnan(desc, NULL);
 		return get_tarval(fc_get_buffer(), fc_get_buffer_length(), mode);
-	} else {
-		assert(0 && "tarval is not floating point");
-		return tarval_bad;
-	}
+	} else
+		panic("mode %F does not support NaN value", mode);
 }
 
 tarval *get_tarval_plus_inf(ir_mode *mode) {
 	assert(mode);
-
-	if (get_mode_n_vector_elems(mode) > 1) {
-		/* vector arithmetic not implemented yet */
-		return tarval_bad;
-	}
+	if (get_mode_n_vector_elems(mode) > 1)
+		panic("vector arithmetic not implemented yet");
 
 	if (get_mode_sort(mode) == irms_float_number) {
-		switch(get_mode_size_bits(mode)) {
-		case 32:
-			fc_get_plusinf(8, 23, NULL);
-			break;
-		case 64:
-			fc_get_plusinf(11, 52, NULL);
-			break;
-		case 80:
-			fc_get_plusinf(15, 64, NULL);
-			break;
-		}
+		const ieee_descriptor_t *desc = get_descriptor(mode);
+		fc_get_plusinf(desc, NULL);
 		return get_tarval(fc_get_buffer(), fc_get_buffer_length(), mode);
-	} else {
-		assert(0 && "tarval is not floating point");
-		return tarval_bad;
-	}
+	} else
+		panic("mode %F does not support +inf value", mode);
 }
 
 tarval *get_tarval_minus_inf(ir_mode *mode) {
 	assert(mode);
 
-	if (get_mode_n_vector_elems(mode) > 1) {
-		/* vector arithmetic not implemented yet */
-		return tarval_bad;
-	}
+	if (get_mode_n_vector_elems(mode) > 1)
+		panic("vector arithmetic not implemented yet");
 
 	if (get_mode_sort(mode) == irms_float_number) {
-		switch(get_mode_size_bits(mode)) {
-		case 32:
-			fc_get_minusinf(8, 23, NULL);
-			break;
-		case 64:
-			fc_get_minusinf(11, 52, NULL);
-			break;
-		case 80:
-			fc_get_minusinf(15, 64, NULL);
-			break;
-		}
+		const ieee_descriptor_t *desc = get_descriptor(mode);
+		fc_get_minusinf(desc, NULL);
 		return get_tarval(fc_get_buffer(), fc_get_buffer_length(), mode);
-	} else {
-		assert(0 && "tarval is not floating point");
-		return tarval_bad;
-	}
+	} else
+		panic("mode %F does not support -inf value", mode);
 }
 
 /*
- * Arithmethic operations on tarvals ========================================
+ * Arithmetic operations on tarvals ========================================
  */
 
 /*
@@ -718,11 +701,8 @@ tarval *get_tarval_minus_inf(ir_mode *mode) {
 int tarval_is_negative(tarval *a) {
 	assert(a);
 
-	if (get_mode_n_vector_elems(a->mode) > 1) {
-		/* vector arithmetic not implemented yet */
-		assert(0 && "tarval_is_negative is not allowed for vector modes");
-		return 0;
-	}
+	if (get_mode_n_vector_elems(a->mode) > 1)
+		panic("vector arithmetic not implemented yet");
 
 	switch (get_mode_sort(a->mode)) {
 	case irms_int_number:
@@ -734,8 +714,7 @@ int tarval_is_negative(tarval *a) {
 		return fc_is_negative(a->value);
 
 	default:
-		assert(0 && "not implemented");
-		return 0;
+		panic("mode %F does not support negation value", a->mode);
 	}
 }
 
@@ -780,7 +759,7 @@ pn_Cmp tarval_cmp(tarval *a, tarval *b) {
 	assert(b);
 
 	if (a == tarval_bad || b == tarval_bad) {
-		assert(0 && "Comparison with tarval_bad");
+		panic("Comparison with tarval_bad");
 		return pn_Cmp_False;
 	}
 
@@ -805,8 +784,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!
@@ -836,13 +816,15 @@ pn_Cmp tarval_cmp(tarval *a, tarval *b) {
  * convert to other mode
  */
 tarval *tarval_convert_to(tarval *src, ir_mode *dst_mode) {
-	char *buffer;
-	fp_value *res;
+	char                    *buffer;
+	fp_value                *res;
+	const ieee_descriptor_t *desc;
 
 	assert(src);
 	assert(dst_mode);
 
-	if (src->mode == dst_mode) return src;
+	if (src->mode == dst_mode)
+		return src;
 
 	if (get_mode_n_vector_elems(src->mode) > 1) {
 		/* vector arithmetic not implemented yet */
@@ -859,19 +841,8 @@ tarval *tarval_convert_to(tarval *src, ir_mode *dst_mode) {
 	case irms_float_number:
 		switch (get_mode_sort(dst_mode)) {
 		case irms_float_number:
-			switch (get_mode_size_bits(dst_mode)) {
-			case 32:
-				fc_cast(src->value, 8, 23, NULL);
-				break;
-			case 64:
-				fc_cast(src->value, 11, 52, NULL);
-				break;
-			case 80:
-				fc_cast(src->value, 15, 64, NULL);
-				break;
-			default:
-				break;
-			}
+			desc = get_descriptor(dst_mode);
+			fc_cast(src->value, desc, NULL);
 			return get_tarval(fc_get_buffer(), fc_get_buffer_length(), dst_mode);
 
 		case irms_int_number:
@@ -883,7 +854,7 @@ tarval *tarval_convert_to(tarval *src, ir_mode *dst_mode) {
 				res = fc_rnd(src->value, NULL);
 				break;
 			default:
-				assert(0);
+				panic("Unsupported float to int conversion mode in tarval_convert_to()");
 				break;
 			}
 			buffer = alloca(sc_get_buffer_length());
@@ -923,17 +894,8 @@ tarval *tarval_convert_to(tarval *src, ir_mode *dst_mode) {
 			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(dst_mode)) {
-			case 32:
-				fc_val_from_str(buffer, 0, 8, 23, NULL);
-				break;
-			case 64:
-				fc_val_from_str(buffer, 0, 11, 52, NULL);
-				break;
-			case 80:
-				fc_val_from_str(buffer, 0, 15, 64, NULL);
-				break;
-			}
+			desc = get_descriptor(dst_mode);
+			fc_val_from_str(buffer, 0, desc, NULL);
 			return get_tarval(fc_get_buffer(), fc_get_buffer_length(), dst_mode);
 
 		default:
@@ -1013,8 +975,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);
@@ -1032,14 +995,22 @@ tarval *tarval_add(tarval *a, tarval *b) {
 
 	assert(a);
 	assert(b);
-	assert(a->mode == b->mode);
 
 	if (get_mode_n_vector_elems(a->mode) > 1 || get_mode_n_vector_elems(b->mode) > 1) {
 		/* vector arithmetic not implemented yet */
 		return tarval_bad;
 	}
 
+	if (mode_is_reference(a->mode) && a->mode != b->mode) {
+		b = tarval_convert_to(b, a->mode);
+	} else if (mode_is_reference(b->mode) && b->mode != a->mode) {
+		a = tarval_convert_to(a, b->mode);
+	}
+
+	assert(a->mode == b->mode);
+
 	switch (get_mode_sort(a->mode)) {
+	case irms_reference:
 	case irms_int_number:
 		/* modes of a,b are equal, so result has mode of a as this might be the character */
 		buffer = alloca(sc_get_buffer_length());
@@ -1047,7 +1018,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);
@@ -1061,18 +1032,27 @@ tarval *tarval_add(tarval *a, tarval *b) {
 /*
  * subtraction
  */
-tarval *tarval_sub(tarval *a, tarval *b) {
-	char *buffer;
+tarval *tarval_sub(tarval *a, tarval *b, ir_mode *dst_mode) {
+	char    *buffer;
 
 	assert(a);
 	assert(b);
-	assert(a->mode == b->mode);
 
 	if (get_mode_n_vector_elems(a->mode) > 1 || get_mode_n_vector_elems(b->mode) > 1) {
 		/* vector arithmetic not implemented yet */
 		return tarval_bad;
 	}
+
+	if (dst_mode != NULL) {
+		if (a->mode != dst_mode)
+			a = tarval_convert_to(a, dst_mode);
+		if (b->mode != dst_mode)
+			b = tarval_convert_to(b, dst_mode);
+	}
+	assert(a->mode == b->mode);
+
 	switch (get_mode_sort(a->mode)) {
+	case irms_reference:
 	case irms_int_number:
 		/* modes of a,b are equal, so result has mode of a as this might be the character */
 		buffer = alloca(sc_get_buffer_length());
@@ -1080,7 +1060,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);
@@ -1114,7 +1094,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);
@@ -1133,7 +1113,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) {
@@ -1238,8 +1218,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);
@@ -1263,7 +1244,7 @@ tarval *tarval_and(tarval *a, tarval *b) {
 
 	/* works even for vector modes */
 
-	switch(get_mode_sort(a->mode)) {
+	switch (get_mode_sort(a->mode)) {
 	case irms_internal_boolean:
 		return (a == tarval_b_false) ? a : b;
 
@@ -1280,7 +1261,7 @@ tarval *tarval_and(tarval *a, tarval *b) {
 /*
  * bitwise or
  */
-tarval *tarval_or (tarval *a, tarval *b) {
+tarval *tarval_or(tarval *a, tarval *b) {
 	assert(a);
 	assert(b);
 	assert(a->mode == b->mode);
@@ -1407,9 +1388,9 @@ tarval *tarval_shrs(tarval *a, tarval *b) {
 }
 
 /*
- * bitwise rotation
+ * bitwise rotation to left
  */
-tarval *tarval_rot(tarval *a, tarval *b) {
+tarval *tarval_rotl(tarval *a, tarval *b) {
 	char *temp_val = NULL;
 
 	assert(a);
@@ -1429,7 +1410,7 @@ tarval *tarval_rot(tarval *a, tarval *b) {
 	} else
 		temp_val = (char*)b->value;
 
-	sc_rot(a->value, temp_val, get_mode_size_bits(a->mode), mode_is_signed(a->mode), NULL);
+	sc_rotl(a->value, temp_val, get_mode_size_bits(a->mode), mode_is_signed(a->mode), NULL);
 	return get_tarval(sc_get_buffer(), sc_get_buffer_length(), a->mode);
 }
 
@@ -1437,6 +1418,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();
 }
 
@@ -1512,7 +1494,15 @@ int tarval_snprintf(char *buf, size_t len, tarval *tv) {
 	case irms_control_flow:
 	case irms_memory:
 	case irms_auxiliary:
-		return snprintf(buf, len, "<TV_OVERFLOW_BAD>");
+		if (tv == tarval_bad)
+			return snprintf(buf, len, "<TV_BAD>");
+		if (tv == tarval_undefined)
+			return snprintf(buf, len, "<TV_UNDEF>");
+		if (tv == tarval_unreachable)
+			return snprintf(buf, len, "<TV_UNREACHABLE>");
+		if (tv == tarval_reachable)
+			return snprintf(buf, len, "<TV_REACHABLE>");
+		return snprintf(buf, len, "<TV_??""?>");
 	}
 
 	return 0;
@@ -1630,6 +1620,15 @@ 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) {
+	const ieee_descriptor_t *desc = get_descriptor(mode);
+	return fc_can_lossless_conv_to(tv->value, desc);
+}
+
 /* Set the immediate precision for IEEE-754 results. */
 unsigned tarval_ieee754_set_immediate_precision(unsigned bits) {
 	return fc_set_immediate_precision(bits);
@@ -1640,6 +1639,17 @@ unsigned tarval_ieee754_get_exact(void) {
 	return fc_is_exact();
 }
 
+/* Return the size of the mantissa in bits (including possible
+   implicit bits) for the given mode. */
+unsigned tarval_ieee754_get_mantissa_size(const ir_mode *mode) {
+	const ieee_descriptor_t *desc;
+
+	assert(get_mode_arithmetic(mode) == irma_ieee754);
+	desc = get_descriptor(mode);
+
+	return desc->mantissa_size + desc->explicit_one;
+}
+
 /* check if its the a floating point NaN */
 int tarval_is_NaN(tarval *tv) {
 	if (! mode_is_float(tv->mode))
@@ -1700,38 +1710,48 @@ static const tarval_mode_info hex_output = {
 /*
  * Initialization of the tarval module: called before init_mode()
  */
-void init_tarval_1(long null_value) {
+void init_tarval_1(long null_value, int support_quad_precision) {
+	/* if these assertion fail, tarval_is_constant() will follow ... */
+	assert(tarval_b_false == &reserved_tv[0] && "b_false MUST be the first reserved tarval!");
+	assert(tarval_b_true  == &reserved_tv[1] && "b_true MUST be the second reserved tarval!");
+
 	_null_value = 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) */
-	init_strcalc(68);
-	init_fltcalc(0);
+	/* calls init_strcalc() with needed size */
+	init_fltcalc(support_quad_precision ? 112 : 64);
 }
 
 /*
  * Initialization of the tarval module: called after init_mode()
  */
 void init_tarval_2(void) {
-	tarval_bad->kind        = k_tarval;
-	tarval_bad->mode        = mode_BAD;
-	tarval_bad->value       = INT_TO_PTR(resid_tarval_bad);
+	tarval_bad->kind          = k_tarval;
+	tarval_bad->mode          = mode_BAD;
+	tarval_bad->value         = INT_TO_PTR(resid_tarval_bad);
+
+	tarval_undefined->kind    = k_tarval;
+	tarval_undefined->mode    = mode_ANY;
+	tarval_undefined->value   = INT_TO_PTR(resid_tarval_undefined);
+
+	tarval_b_true->kind       = k_tarval;
+	tarval_b_true->mode       = mode_b;
+	tarval_b_true->value      = INT_TO_PTR(resid_tarval_b_true);
 
-	tarval_undefined->kind  = k_tarval;
-	tarval_undefined->mode  = mode_ANY;
-	tarval_undefined->value = INT_TO_PTR(resid_tarval_undefined);
+	tarval_b_false->kind      = k_tarval;
+	tarval_b_false->mode      = mode_b;
+	tarval_b_false->value     = INT_TO_PTR(resid_tarval_b_false);
 
-	tarval_b_true->kind     = k_tarval;
-	tarval_b_true->mode     = mode_b;
-	tarval_b_true->value    = INT_TO_PTR(resid_tarval_b_true);
+	tarval_unreachable->kind  = k_tarval;
+	tarval_unreachable->mode  = mode_X;
+	tarval_unreachable->value = INT_TO_PTR(resid_tarval_unreachable);
 
-	tarval_b_false->kind    = k_tarval;
-	tarval_b_false->mode    = mode_b;
-	tarval_b_false->value   = INT_TO_PTR(resid_tarval_b_false);
+	tarval_reachable->kind    = k_tarval;
+	tarval_reachable->mode    = mode_X;
+	tarval_reachable->value   = INT_TO_PTR(resid_tarval_reachable);
 
 	/*
 	 * assign output modes that are compatible with the
@@ -1756,6 +1776,10 @@ void finish_tarval(void) {
 	del_set(values);  values = NULL;
 }
 
+int (is_tarval)(const void *thing) {
+	return _is_tarval(thing);
+}
+
 /****************************************************************************
  *   end of tv.c
  ****************************************************************************/