From: Danielle Bolan Date: Tue, 1 Jul 2014 14:03:38 +0000 (+0200) Subject: Project_out working X-Git-Tag: 1.0~183 X-Git-Url: https://scm.cri.mines-paristech.fr/git/linpy.git/commitdiff_plain/a8bb1d373098fbaaac285cc7aeea6b9b380f2fd7?ds=sidebyside;hp=c774b95a5e3c664d7adaeac83aace2b045d38a9f Project_out working --- diff --git a/pypol/domains.py b/pypol/domains.py index 4a1f9aa..6c10de3 100644 --- a/pypol/domains.py +++ b/pypol/domains.py @@ -154,15 +154,19 @@ class Domain: def drop_dims(self, dims): # use to remove certain variables use isl_set_drop_constraints_involving_dims instead? - from .polyhedra import Polyhedron + from .polyhedra import Polyhedron dims = list(dims) symbols = list(self.symbols) + print(symbols) islset = self._toislset(self.polyhedra, self.symbols) for dim in dims: if dim in symbols: - num = symbols.index(dim) - islbset = libisl.isl_set_drop_constraints_involving_dims(islset, libisl.isl_dim_set, num, 1) - return Polyhedron._fromislset(islbset, self.symbols) + first = symbols.index(dim) + islbset = libisl.isl_set_project_out(islset, libisl.isl_dim_set, first, 1) + symbols.__delitem__(first) + else: + islbset = libisl.isl_set_project_out(islset, libisl.isl_dim_set, 0, 0) + return Polyhedron._fromislset(islbset, symbols) def sample(self): from .polyhedra import Polyhedron