diff --git a/islpy/__init__.py b/islpy/__init__.py
index d89a6e1eacb24f06ea5bac4896fb3eafc91015b1..ee5d5a9de9839e4098a18d408690f46a724714ee 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