From 692f773acbb544c808b3447553b9d8391a58f50b Mon Sep 17 00:00:00 2001
From: skimo <skimo>
Date: Thu, 30 Sep 2004 12:58:51 +0000
Subject: [PATCH] ensure all rays are positive

---
 barvinok.cc | 17 +++++++++++++++++
 1 file changed, 17 insertions(+)

diff --git a/barvinok.cc b/barvinok.cc
index ee03289..279efec 100644
--- a/barvinok.cc
+++ b/barvinok.cc
@@ -1567,6 +1567,22 @@ static bool Polyhedron_is_infinite(Polyhedron *P, unsigned nparam)
     return r <  P->NbRays;
 }
 
+/* Check whether all rays point in the positive directions
+ * for the parameters
+ */
+static bool Polyhedron_has_positive_rays(Polyhedron *P, unsigned nparam)
+{
+    int r;
+    for (r = 0; r < P->NbRays; ++r)
+	if (value_zero_p(P->Ray[r][P->Dimension+1])) {
+	    int i;
+	    for (i = P->Dimension - nparam; i < P->Dimension; ++i)
+		if (value_neg_p(P->Ray[r][i+1]))
+		    return false;
+	}
+    return true;
+}
+
 evalue* barvinok_enumerate_ev(Polyhedron *P, Polyhedron* C, unsigned MaxRays)
 {
     //P = unfringe(P, MaxRays);
@@ -3343,6 +3359,7 @@ gen_fun * barvinok_series(Polyhedron *P, Polyhedron* C, unsigned MaxRays)
 
     assert(!Polyhedron_is_infinite(P, nparam));
     assert(P->NbBid == 0);
+    assert(Polyhedron_has_positive_rays(P, nparam));
     assert(P->NbEq == 0);
 
     dim = P->Dimension;
-- 
2.11.4.GIT