Update drop_dims, now works with more than 2 dims passed
[linpy.git] / pypol / domains.py
index 4a1f9aa..3bcfa24 100644 (file)
@@ -154,15 +154,27 @@ 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 
+        n = 0 
         dims = list(dims)
         symbols = list(self.symbols)
         islset = self._toislset(self.polyhedra, self.symbols)
         for dim in dims:
+            dim_index = dims.index(dim)
             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) 
+                try:
+                    for dim in dims:
+                        if symbols[first+1] is dims[dim_index+1]: #check if next value in symbols is same as next value in dims
+                            n += 1
+                        islbset = libisl.isl_set_project_out(islset, libisl.isl_dim_set, first, n) 
+                        symbols.remove(dim)
+                except:
+                    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