The required lattice basis can be computed using the LLL-algorithm ( lll_lib, see see [Coh93]). For the computation of the saturation, there are various possibilities described in the
C.6.2 Algorithms Various algorithms for computing toric ideals. C.6.3 The Buchberger algorithm for toric ideals Specializing it for toric ideals.