Probably this answer comes too late, but anyway:
What you want to do can be divided into two tasks:
- Find a Hilbert basis of (the dual of) your cone.
- Compute the toric ideal for the Hilbert basis.
I.e. let $M=\mathbb Z^n$ be the dual lattice of $N\times\mathbb Z$. You start with your cone $\sigma$, dualize it and obtain $\sigma^\vee$. Now you compute a minimal generating set $H:=\{h_1,\ldots,h_m\}\subseteq M$ of the semigroup $\sigma^\vee\cap M$. Then you get a surjection
$$
k[x_1,\ldots, x_m]\to k[\sigma^\vee\cap M],\ x_i\mapsto x^{h_i}
$$
The kernel of this morphism is the ideal $I$ you want. If you denote $H$ as a matrix, it can be described as
$$
I=( x^u - x^v\ | u,v\in\mathbb Z^m_{\ge 0},\ H\cdot(u-v)=0).
$$
Thus you can obtain generators of $I$ by just taking all vectors of $\ker H$, splitting them into positive and negative part, and then take the corresponding binomials.
This you probably knew already, but now for the software:
- I know two (three) frameworks for computing Hilbert bases: 4ti2 and Normaliz. You can enter your cone via rays or give the facet inequalities, i.e. the rays of $\sigma$. Thus, dualization can be done on the fly. Macaulay2 is also able to compute Hilbert bases via the Polyhedra package, but the first two are probably faster.
- To compute the toric ideal you can use the groebner method of 4ti2 (see Usage->groebner on the 4ti2 page). Singular provides some algorithms for computing toric ideals on its own, via the toric.lib package. If performance is an issue, 4ti2 is probably faster.
Both Macaulay2 and Singular provide packages interfacing 4ti2, which makes dealing with the algebraic aspect easier.
