X-Git-Url: https://scm.cri.ensmp.fr/git/linpy.git/blobdiff_plain/114cd46edf08bcf13c5f489c449db98f6098c468..2de96003a48ae1ed91c853f38c72374a832824de:/pypol/domains.py diff --git a/pypol/domains.py b/pypol/domains.py index be9dd4b..122d428 100644 --- a/pypol/domains.py +++ b/pypol/domains.py @@ -340,6 +340,8 @@ class Domain(GeometricObject): Return a list of vertices for this Polygon. """ from .polyhedra import Polyhedron + if not self.isbounded(): + raise ValueError('domain must be bounded') islbset = self._toislbasicset(self.equalities, self.inequalities, self.symbols) vertices = libisl.isl_basic_set_compute_vertices(islbset); vertices = islhelper.isl_vertices_vertices(vertices) @@ -443,6 +445,9 @@ class Domain(GeometricObject): return sorted(points, key=angles.get) def faces(self): + """ + Returns the vertices of the faces of a polyhedra. + """ faces = [] for polyhedron in self.polyhedra: vertices = polyhedron.vertices() @@ -503,6 +508,7 @@ class Domain(GeometricObject): axes.set_zlim(zmin, zmax) return axes + def plot(self, plot=None, **kwargs): """ Display plot of this set. @@ -523,6 +529,10 @@ class Domain(GeometricObject): return False def subs(self, symbol, expression=None): + """ + Subsitute the given value into an expression and return the resulting + expression. + """ polyhedra = [polyhedron.subs(symbol, expression) for polyhedron in self.polyhedra] return Domain(*polyhedra)