polyhedron = self.polyhedron
         for x, xprime, dx in zip(self.range_symbols, self.domain_symbols, delta_symbols):
             polyhedron &= Eq(dx, xprime - x)
-        polyhedron = polyhedron.project_out(self.symbols)
+        polyhedron = polyhedron.project(self.symbols)
         equalities, inequalities = [], []
         for equality in polyhedron.equalities:
             equality += (k-1) * equality.constant
             inequality += (k-1) * inequality.constant
             inequalities.append(inequality)
         polyhedron = Polyhedron(equalities, inequalities) & Ge(k, 0)
-        polyhedron = polyhedron.project_out([k])
+        polyhedron = polyhedron.project([k])
         for x, xprime, dx in zip(self.range_symbols, self.domain_symbols, delta_symbols):
             polyhedron &= Eq(dx, xprime - x)
-        polyhedron = polyhedron.project_out(delta_symbols)
+        polyhedron = polyhedron.project(delta_symbols)
         return Transformer(polyhedron, self.range_symbols, self.domain_symbols)
 
 
 
 print('sq6:', sq6)
 print('sq6 simplified:', sq6.sample())
 print()
-print(universe.project_out([x]))
-print('sq7 with out constraints involving y and a', sq7.project_out([a, z, x, y])) #drops dims that are passed
+print(universe.project([x]))
+print('sq7 with out constraints involving y and a', sq7.project([a, z, x, y])) #drops dims that are passed
 print()
 print('sq1 has {} parameters'.format(sq1.num_parameters()))
 print()
 
         self.assertFalse(self.empty == self.universe)
 
     def test_isdisjoint(self):
-        self.assertFalse(self.square1.isdisjoint(self.square2))        
+        self.assertFalse(self.square1.isdisjoint(self.square2))
         self.assertFalse(self.universe.isdisjoint(self.square1))
         self.assertTrue(self.square1.isdisjoint(self.square5))
         self.assertTrue(self.empty.isdisjoint(self.square1))
         self.assertFalse(self.square3 <= self.square4)
         self.assertTrue(self.empty <= self.square1)
         self.assertTrue(self.square1 <= self.universe)
-        
+
     def test_lt(self):
         self.assertTrue(self.square4 < self.square3)
         self.assertFalse(self.square3 < self.square4)
         self.assertEqual(self.universe.polyhedral_hull(), self.universe)
         self.assertEqual(self.empty.polyhedral_hull(), self.empty)
 
-    def test_project_out(self):
-        self.assertEqual(self.square1.project_out(symbols('x')), self.dropped)
-        self.assertEqual(self.square1.project_out(symbols('x y')), self.universe)
-        self.assertEqual(self.universe.project_out([]), self.universe)
-        self.assertEqual(self.empty.project_out([]), Empty)
+    def test_project(self):
+        self.assertEqual(self.square1.project(symbols('x')), self.dropped)
+        self.assertEqual(self.square1.project(symbols('x y')), self.universe)
+        self.assertEqual(self.universe.project([]), self.universe)
+        self.assertEqual(self.empty.project([]), Empty)
 
     def test_simplify(self):
         self.assertEqual(self.universe.simplify(), self.universe)
         self.assertEqual(self.square6.sample(), self.sample)
         self.assertEqual(self.empty.sample(), Empty)
         self.assertEqual(self.universe.sample(), self.universe)
-        
+
     def test_intersection(self):
         self.assertEqual(self.square1.intersection(self.square2), self.intersection)
 
         self.assertEqual(self.square1.lexmax(), self.lexmax)
         self.assertEqual(self.universe.lexmax(), self.universe)
         self.assertEqual(self.empty.lexmax(), Empty)
-        
+
     def test_num_parameters(self):
         self.assertEqual(self.square1.num_parameters(), 2)
         self.assertEqual(self.empty.num_parameters(), 0)
         self.assertEqual(self.universe.num_parameters(), 0)
-        
+
     def test_involves_dims(self):
         self.assertTrue(self.square1.involves_dims(symbols('x y')))
         self.assertFalse(self.empty.involves_dims(symbols('x')))