examples/diamond.py

   1 #!/usr/bin/env python3
   2
   3 from pypol import *
   4
   5 x, y, z = symbols('x y z')
   6 diam = Ge(y, x - 1) & Le(y, x + 1) & Ge(y, -x - 1) & Le(y, -x + 1)
   7 print('diamond:', diam)
   8 print()
   9 rhom1 = Le(0, x) & Le(x, 3) & Le(0, y) & Le(y, 3) & Le(0, z) & Le(z, 3) \
  10 & Le(z - 2, x) & Ge(z + 2, x) & Ge(z - 1, -x) & Le(z - 5, -x) \
  11 & Le(z - 2, y) & Ge(z + 2, y) & Ge(z - 1, -y) & Le(z - 5, -y) \
  12 & Le(y - 2, x) & Ge(y + 2, x) & Ge(y - 1, -x) & Le(y - 5, -x)
  13 rhom1.plot()
  14 rhom2 = rhom1 & Le(x + y + z, 7) & Ge(-2, -x - y - z ) \
  15 & Le(x + y - z, 4) & Ge(x + y - z, -1) \
  16 & Le(x - y + z, 4) & Ge(x - y + z, -1) \
  17 & Le(-x + y + z, 4) & Ge(-x + y + z, -1)
  18 rhom2.plot()