X-Git-Url: https://scm.cri.ensmp.fr/git/linpy.git/blobdiff_plain/90d7f21ee1486524f499f66ef546a24331b2009c..51e97eade63b2f4c7b500feb503436cc4a886e59:/examples/nsad2010.py

diff --git a/examples/nsad2010.py b/examples/nsad2010.py
index 8194618..f310f9d 100755
--- a/examples/nsad2010.py
+++ b/examples/nsad2010.py
@@ -1,5 +1,22 @@
 #!/usr/bin/env python3
 
+"""
+    This file is part of Linpy.
+
+    Linpy is free software: you can redistribute it and/or modify
+    it under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    Linpy is distributed in the hope that it will be useful,
+    but WITHOUT ANY WARRANTY; without even the implied warranty of
+    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+    GNU General Public License for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with Linpy.  If not, see <http://www.gnu.org/licenses/>.
+"""
+
 from pypol import *
 
 
@@ -22,7 +39,7 @@ class Transformer:
         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
@@ -31,10 +48,10 @@ class Transformer:
             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)
 
 
@@ -44,3 +61,5 @@ if __name__ == '__main__':
         [i, j], [iprime, jprime])
     print('T  =', transformer.polyhedron)
     print('T* =', transformer.star().polyhedron)
+    
+# Copyright 2014 MINES ParisTech