Thanks to the MUCH easier install now of things like pynormaliz and latte
(thanks to all who worked on those!), I can now do the following and
related computations nicely.
sage: n=1
sage: P =
Polyhedron(ieqs=[[-(n)/2,1,0,0],[-(n)/2,0,1,0],[(3*n)/2,-1,-1,-1],[0,1
....:
,0,0],[0,0,1,0],[0,0,0,1],[n,-1,0,0],[n,0,-1,0],[n,0,0,-1]],backend='norma
....: liz')
sage: [p.factor() for p in P.ehrhart_quasipolynomial()]
[(1/48) * (t + 2) * (t + 4) * (t + 6), (1/48) * (t - 1) * (t + 1) * (t + 3)]
However, what I really need is an Ehrhart quasi-polynomial for some of the
above inequalities to be *strict* inequalities, and I'm not sure how to do
that without tedious finding of some (not all) faces and subtracting them
off (which could be a nightmare and/or wrong in any case). Unfortunately
changing the non-strict inequalities "by hand" to other numbers gives the
wrong answers (really unsurprising, since it's a different polytope).
Any thoughts? Thanks!
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to sage-support+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/sage-support/c2e8811a-f6da-41f0-a294-31e5167c1336n%40googlegroups.com.
