toricBlowUp -- calculates the stellar subdivision of a polytope at a given face.

• Usage:

  toricBlowUp(P,Q),

  toricBlowUp(P,Q,k),

  toricBlowUp(M,N),

  toricBlowUp(M,N,k)

• Inputs:
  • P, a convex polyhedron,
  • Q, a convex polyhedron,
  • M, a matrix,
  • N, a matrix,
  • k, an integer,
• Outputs:
  • a convex polyhedron,

i1 : P=cayley(matrix{{0,2,0}},matrix{{0,0,2}})

o1 = P

o1 : Polyhedron

i2 : vertices oo

o2 = | 0 2 0 2 |
     | 0 0 1 1 |

              2       4
o2 : Matrix QQ  <-- QQ

i3 : Q=convexHull(matrix{(vertices P)_0})

o3 = Q

o3 : Polyhedron

i4 : toricBlowUp(P,Q,1)

o4 = Polyhedron{...1...}

o4 : Polyhedron

i5 : vertices oo

o5 = | 1 2 0 2 |
     | 0 0 1 1 |

              2       4
o5 : Matrix QQ  <-- QQ