fix serious error in isl_mat_parameter_compression

The old version would sometimes remove valid solutions.
Since isl_mat_parameter_compression is used during simplification
on practically any set containing existentially quantified variables
involved in equalities, any such set could get corrupted.

1 files changed, 30 insertions, 0 deletions

diff --git a/isl_vec.c b/isl_vec.c
index 215e2398..59c26a9e 100644
--- a/isl_vec.c
+++ b/isl_vec.c
@@ -30,6 +30,31 @@ struct isl_vec *isl_vec_copy(struct isl_ctx *ctx, struct isl_vec *vec)
 	return vec;
 }
 
+struct isl_vec *isl_vec_dup(struct isl_ctx *ctx, struct isl_vec *vec)
+{
+	struct isl_vec *vec2;
+
+	if (!vec)
+		return NULL;
+	vec2 = isl_vec_alloc(ctx, vec->size);
+	isl_seq_cpy(vec2->block.data, vec->block.data, vec->size);
+	return vec2;
+}
+
+struct isl_vec *isl_vec_cow(struct isl_ctx *ctx, struct isl_vec *vec)
+{
+	struct isl_vec *vec2;
+	if (!vec)
+		return NULL;
+
+	if (vec->ref == 1)
+		return vec;
+
+	vec2 = isl_vec_dup(ctx, vec);
+	isl_vec_free(ctx, vec);
+	return vec2;
+}
+
 void isl_vec_free(struct isl_ctx *ctx, struct isl_vec *vec)
 {
 	if (!vec)
@@ -54,3 +79,8 @@ void isl_vec_dump(struct isl_ctx *ctx, struct isl_vec *vec,
 	}
 	fprintf(out, "]\n");
 }
+
+void isl_vec_lcm(struct isl_ctx *ctx, struct isl_vec *vec, isl_int *lcm)
+{
+	isl_seq_lcm(vec->block.data, vec->size, lcm);
+}