Implement method Point.aspolyhedron()

diff --git a/pypol/coordinates.py b/pypol/coordinates.py
index 78e8d4a151a3fea215593d29a8409bc49dd1772f..792364862a2fbd77dea9ecb5c30f9ca78963c411 100644 (file)
--- a/pypol/coordinates.py
+++ b/pypol/coordinates.py
@@ -112,6 +112,13 @@ class Point(Coordinates):
         return isinstance(other, Point) and \
             self._coordinates == other._coordinates
 
         return isinstance(other, Point) and \
             self._coordinates == other._coordinates
 

+    def aspolyhedron(self):
+        from .polyhedra import Polyhedron
+        equalities = []
+        for symbol, coordinate in self.coordinates():
+            equalities.append(symbol - coordinate)
+        return Polyhedron(equalities)
+

 
 class Vector(Coordinates):
     """
 
 class Vector(Coordinates):
     """

diff --git a/pypol/tests/test_coordinates.py b/pypol/tests/test_coordinates.py
index 2e14032a750785661f0e4c20321b28bbde10400b..044d7736342faaef63f5b129e178837265dab884 100644 (file)
--- a/pypol/tests/test_coordinates.py
+++ b/pypol/tests/test_coordinates.py
@@ -1,8 +1,9 @@
 import math
 import unittest
 
 import math
 import unittest
 

-from ..linexprs import Symbol

 from ..coordinates import *
 from ..coordinates import *

+from ..linexprs import Symbol
+from ..polyhedra import Eq

 
 
 class TestPoint(unittest.TestCase):
 
 
 class TestPoint(unittest.TestCase):


@@ -29,6 +30,9 @@ class TestPoint(unittest.TestCase):
         self.assertEqual(self.pt1 - self.pt2, Vector({self.x: -5, self.y: -35, self.z: -59}))
         self.assertEqual(self.pt1 - self.vec1, Point({self.x: -10, self.y: -25, self.z: -39}))
 
         self.assertEqual(self.pt1 - self.pt2, Vector({self.x: -5, self.y: -35, self.z: -59}))
         self.assertEqual(self.pt1 - self.vec1, Point({self.x: -10, self.y: -25, self.z: -39}))
 

+    def test_aspolyhedron(self):
+        self.assertEqual(self.pt1.aspolyhedron(), Eq(self.x, 10) & Eq(self.y, 5) & Eq(self.z, 1))
+

 
 class TestVector(unittest.TestCase):
 
 
 class TestVector(unittest.TestCase):