equal_polyhedral for non-pointed cones

Questions and problems about using polymake go here.
simon
Posts: 13
Joined: 18 Mar 2012, 19:47

equal_polyhedral for non-pointed cones

Postby simon » 18 Mar 2012, 20:23

Hi,

I just tried

$c = new Cone<Rational>(RAYS=>[[1,0],[0,1],[-1,-1]]);
print equal_polyhedra($c, $c);

but this produces no output. Is this a bug?

paffenholz
Developer
Posts: 212
Joined: 24 Dec 2010, 13:47

Re: equal_polyhedral for non-pointed cones

Postby paffenholz » 20 Mar 2012, 10:20

You input is not valid, the section RAYS contains only rays for the pointed part of the cone, while your rays span a two dimensional linear space. This would be recorded in the property LINEALITY_SPACE (by giving a basis of the space).

There is a distinction between these reduced properties RAYS/LINEALITY_SPACE and the input properties INPUT_RAYS and INPUT_LINEALITY. The latter need not be reduced, and the section INPUT_RAYS may contain rays that lie in the lineality space. So you could define your cone using INPUT_RAYS instead of RAYS. It is generally advisable to use these input properties and let polymake do the reduction, as polymake never checks consistency of the input.

Nevertheless, even with the correct definition it might not work, depending on your current configuration. You indeed discovered a bug, polymake doesn't always properly handle the case of cones that have no pointed part.

simon
Posts: 13
Joined: 18 Mar 2012, 19:47

Re: equal_polyhedral for non-pointed cones

Postby simon » 20 Mar 2012, 15:38

thanks! I'm still learning.


Return to “Helpdesk”