From bd13bed9199c993106eee07179bd84b3032cc7d5 Mon Sep 17 00:00:00 2001
From: Andreas Kloeckner <inform@tiker.net>
Date: Wed, 5 Jul 2017 04:25:25 -0500
Subject: [PATCH] PwQPolynomial.eval_from_dict: Check whether eval point is in
 domain

---
 islpy/__init__.py | 4 ++++
 1 file changed, 4 insertions(+)

diff --git a/islpy/__init__.py b/islpy/__init__.py
index d89a6e1..ee5d5a9 100644
--- a/islpy/__init__.py
+++ b/islpy/__init__.py
@@ -687,6 +687,10 @@ def _add_functionality():
             pt = pt.set_coordinate_val(
                 dim_type.param, i, value_dict[par_name])
 
+        if not (Set.from_point(pt) <= self.domain()):
+            raise ValueError("evaluation point outside of domain of "
+                    "definition of piecewise quasipolynomial")
+
         return self.eval(pt).to_python()
 
     PwQPolynomial.eval_with_dict = pwqpolynomial_eval_with_dict
-- 
GitLab