From fcc8d1c4ad2ad54dbe1deb6bdf050618bc3c7fc5 Mon Sep 17 00:00:00 2001
From: Sven Verdoolaege <skimo@kotnet.org>
Date: Sat, 15 Mar 2008 13:22:11 +0100
Subject: [PATCH] evalue_split_domains_into_orthants: don't reduce evalue after
 splitting

reduce_evalue will remove any domain where the value is zero, but
some applications, most notably those dealing with bounds, depend
on these domains remaining in the input evalues.
---
 evalue.c | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/evalue.c b/evalue.c
index e18740e..10d45a5 100644
--- a/evalue.c
+++ b/evalue.c
@@ -3577,6 +3577,7 @@ evalue *evalue_sum(evalue *e, int nvar, unsigned MaxRays)
     }
 
     evalue_split_domains_into_orthants(e, MaxRays);
+    reduce_evalue(e);
     evalue_frac2floor2(e, 0);
     evalue_set_si(res, 0, 1);
 
@@ -3901,7 +3902,6 @@ void evalue_split_domains_into_orthants(evalue *e, unsigned MaxRays)
 	free_evalue_refs(&split);
 	Matrix_Free(C);
     }
-    reduce_evalue(e);
 }
 
 static evalue *find_fractional_with_max_periods(evalue *e, Polyhedron *D,
-- 
2.11.4.GIT