X-Git-Url: https://scm.cri.ensmp.fr/git/linpy.git/blobdiff_plain/bf0ba2bac20a5c3f288001a78275e83a0106a02a..1b9c772ffc21a71d146a093537d371d63be16c75:/examples/squares.py diff --git a/examples/squares.py b/examples/squares.py index 0f4e1a5..e140e30 100755 --- a/examples/squares.py +++ b/examples/squares.py @@ -2,15 +2,15 @@ from pypol import * -x, y = symbols('x y') +a, x, y, z = symbols('a x y z') sq1 = Le(0, x) & Le(x, 2) & Le(0, y) & Le(y, 2) sq2 = Le(2, x) & Le(x, 4) & Le(2, y) & Le(y, 4) - sq3 = Le(0, x) & Le(x, 3) & Le(0, y) & Le(y, 3) sq4 = Le(1, x) & Le(x, 2) & Le(1, y) & Le(y, 2) sq5 = Le(1, x) & Le(x, 2) & Le(1, y) sq6 = Le(1, x) & Le(x, 2) & Le(1, y) & Eq(y, 3) +sq7 = Le(0, x) & Le(x, 2) & Le(0, y) & Eq(z, 2) & Le(a, 3) u = Polyhedron([]) x = sq1 - sq2 @@ -62,4 +62,8 @@ print('is sq5 bounded?', sq5.isbounded()) #unbounded should return False print() print('sq6:', sq6) print('sq6 simplified:', sq6.sample()) +print() +#print(u.drop_dims(' ')) +print('sq7 with out constraints involving y and a', sq7.drop_dims('y a')) #drops dims that are passed +