#!/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 *
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)
[i, j], [iprime, jprime])
print('T =', transformer.polyhedron)
print('T* =', transformer.star().polyhedron)
+
+# Copyright 2014 MINES ParisTech