X-Git-Url: https://scm.cri.ensmp.fr/git/linpy.git/blobdiff_plain/49b7b1b00f4cbbf5110ec77e72df004764f2b4bd..8fe086b5ae03df5fadf1ec7de1aa41c0ace4a3ae:/pypol/domains.py

diff --git a/pypol/domains.py b/pypol/domains.py
index be87ffb..bbb0946 100644
--- a/pypol/domains.py
+++ b/pypol/domains.py
@@ -154,14 +154,22 @@ 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 = 1     
         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:
-                first = symbols.index(dim)
-                islbset = libisl.isl_set_project_out(islset, libisl.isl_dim_set, first, 1)  
+                first = symbols.index(dim) 
+                try:
+                    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.__delitem__(first)
+                except:
+                    islbset = libisl.isl_set_project_out(islset, libisl.isl_dim_set, first, n)  
                 symbols.__delitem__(first)     
             else:
                 islbset = libisl.isl_set_project_out(islset, libisl.isl_dim_set, 0, 0)