hilbertBasis -- computes the Hilbert basis of a Cone

Synopsis

• Usage:

  HB = hilbertBasis C

• Inputs:
  • C, a convex rational cone
• Outputs:
  • L, a list, containing the elements of the Hilbert basis

Description

The Hilbert basis of the cone C is computed by the Project-and-Lift-algorithm by Raymond Hemmecke (see below). It computes a Hilbert basis of the cone modulo the lineality space, so it returns the list of one column matrices that give the Hilbert basis of the Cone if one adds the basis of the lineality space and its negative. For the Project-and-Lift-algorithm see:

Raymond Hemmecke's On the computation of Hilbert bases of cones, in A. M. Cohen, X.-S. Gao, and N. Takayama, editors, Mathematical Software, ICMS 2002, pages 307317. World Scientific, 2002.

i1 : C = posHull matrix {{1,2},{2,1}}

o1 = {ambient dimension => 2           }
      dimension of lineality space => 0
      dimension of the cone => 2
      number of facets => 2
      number of rays => 2

o1 : Cone

i2 : hilbertBasis C

o2 = {| 1 |, | 2 |, | 1 |}
      | 1 |  | 1 |  | 2 |

o2 : List