X-Git-Url: http://review.tizen.org/git/?a=blobdiff_plain;f=isl_lp_piplib.c;h=c95af4c161e5b70d8442ec47fc1e7eb59bc36236;hb=de51a9bc4da5dd3f1f9f57c2362da6f9752c44e0;hp=b4e9bd9647c94476885bb635301c86004dce4914;hpb=67497edd0ff3657eab38d3b71b3fff7131c0f412;p=platform%2Fupstream%2Fisl.git

diff --git a/isl_lp_piplib.c b/isl_lp_piplib.c
index b4e9bd9..c95af4c 100644
--- a/isl_lp_piplib.c
+++ b/isl_lp_piplib.c
@@ -1,11 +1,63 @@
-#include "isl_map.h"
-#include "isl_lp.h"
+/*
+ * Copyright 2008-2009 Katholieke Universiteit Leuven
+ *
+ * Use of this software is governed by the MIT license
+ *
+ * Written by Sven Verdoolaege, K.U.Leuven, Departement
+ * Computerwetenschappen, Celestijnenlaan 200A, B-3001 Leuven, Belgium
+ */
+
+#include <isl/map.h>
+#include <isl/vec.h>
+#include <isl/lp.h>
 #include "isl_piplib.h"
 #include "isl_map_piplib.h"
 
+static void copy_solution(struct isl_vec *vec, int maximize, isl_int *opt,
+	isl_int *opt_denom, PipQuast *sol)
+{
+	int i;
+	PipList *list;
+	isl_int tmp;
+
+	if (opt) {
+		if (opt_denom) {
+			isl_seq_cpy_from_pip(opt,
+				 &sol->list->vector->the_vector[0], 1);
+			isl_seq_cpy_from_pip(opt_denom,
+				 &sol->list->vector->the_deno[0], 1);
+		} else if (maximize)
+			mpz_fdiv_q(*opt, sol->list->vector->the_vector[0],
+					 sol->list->vector->the_deno[0]);
+		else
+			mpz_cdiv_q(*opt, sol->list->vector->the_vector[0],
+					 sol->list->vector->the_deno[0]);
+	}
+
+	if (!vec)
+		return;
+
+	isl_int_init(tmp);
+	isl_int_set_si(vec->el[0], 1);
+	for (i = 0, list = sol->list->next; list; ++i, list = list->next) {
+		isl_seq_cpy_from_pip(&vec->el[1 + i],
+			&list->vector->the_deno[0], 1);
+		isl_int_lcm(vec->el[0], vec->el[0], vec->el[1 + i]);
+	}
+	for (i = 0, list = sol->list->next; list; ++i, list = list->next) {
+		isl_seq_cpy_from_pip(&tmp, &list->vector->the_deno[0], 1);
+		isl_int_divexact(tmp, vec->el[0], tmp);
+		isl_seq_cpy_from_pip(&vec->el[1 + i],
+			&list->vector->the_vector[0], 1);
+		isl_int_mul(vec->el[1 + i], vec->el[1 + i], tmp);
+	}
+	isl_int_clear(tmp);
+}
+
 enum isl_lp_result isl_pip_solve_lp(struct isl_basic_map *bmap, int maximize,
 				      isl_int *f, isl_int denom, isl_int *opt,
-				      isl_int *opt_denom)
+				      isl_int *opt_denom,
+				      struct isl_vec **vec)
 {
 	enum isl_lp_result res = isl_lp_ok;
 	PipMatrix	*domain = NULL;
@@ -32,23 +84,18 @@ enum isl_lp_result isl_pip_solve_lp(struct isl_basic_map *bmap, int maximize,
 	if (!sol)
 		goto error;
 
-	if (!sol->list)
+	if (vec) {
+		isl_ctx *ctx = isl_basic_map_get_ctx(bmap);
+		*vec = isl_vec_alloc(ctx, 1 + total);
+	}
+	if (vec && !*vec)
+		res = isl_lp_error;
+	else if (!sol->list)
 		res = isl_lp_empty;
 	else if (entier_zero_p(sol->list->vector->the_deno[0]))
 		res = isl_lp_unbounded;
-	else {
-		if (opt_denom) {
-			isl_seq_cpy_from_pip(opt,
-				 &sol->list->vector->the_vector[0], 1);
-			isl_seq_cpy_from_pip(opt_denom,
-				 &sol->list->vector->the_deno[0], 1);
-		} else if (maximize)
-			mpz_fdiv_q(*opt, sol->list->vector->the_vector[0],
-					 sol->list->vector->the_deno[0]);
-		else
-			mpz_cdiv_q(*opt, sol->list->vector->the_vector[0],
-					 sol->list->vector->the_deno[0]);
-	}
+	else
+		copy_solution(*vec, maximize, opt, opt_denom, sol);
 	pip_matrix_free(domain);
 	pip_quast_free(sol);
 	return res;