From: Danielle Bolan <n02702451@hawkmail.newpaltz.edu>
Date: Tue, 1 Jul 2014 14:03:38 +0000 (+0200)
Subject: Project_out working
X-Git-Tag: 1.0~183
X-Git-Url: https://svn.cri.ensmp.fr/git/linpy.git/commitdiff_plain/a8bb1d373098fbaaac285cc7aeea6b9b380f2fd7?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