next | previous | forward | backward | up | top | index | toc | Macaulay2 website
StatePolytope :: polymakeStatePolytope(ZZ,Ideal)

polymakeStatePolytope(ZZ,Ideal) -- computes the mth state polytope of an ideal

Synopsis

Description

See Sturmfels's book Groebner bases and convex polytopes, page 14 for the definition of Statem(I). (The difference between this and State(I) is that for all sufficiently large m, Statem(I) does not distinguish between initial ideals which have the same saturation with regard to the irrelevant ideal, whereas in State(I), these are separated.)
i1 : R = QQ[a..d];
i2 : I = ideal(a*c-b^2,a*d-b*c,b*d-c^2);

o2 : Ideal of R
i3 : polymakeStatePolytope(3,I)

o3 = {{9, 6, 6, 9}, {6, 12, 3, 9}, {4, 14, 5, 7}, {3, 14, 8, 5}, {9, 3, 12,
     ------------------------------------------------------------------------
     6}, {7, 5, 14, 4}, {5, 8, 14, 3}, {3, 12, 12, 3}}

o3 : List