
D.4.9 ffmodstd_lib
 Library:
 ffmodstd.lib
 Purpose:
 Groebner bases of ideals in polynomial rings
over rational function fields
 Authors:
 D.K. Boku boku@mathematik.unikl.de
W. Decker decker@mathematik.unikl.dei
C. Fieker fieker@mathematik.unikl.de
A. Steenpass steenpass@mathematik.unikl.de
 Overview:
 A library for computing a Groebner basis of an ideal in a polynomial
ring over an algebraic function field Q(T):=Q(t_1,...,t_m) using modular
methods and sparse multivariate rational interpolation, where the
t_i are transcendental over Q. The idea is as follows:
Given an ideal I in Q(T)[X], we map I to J via the map sending T to
Tz:=(t_1z+s_1,..., t_mz+s_m) for a suitable point s in Q^m\{(0,...,0)} and for some
extra variable z so that J is an ideal in Q(Tz)[X]. For a suitable point b in
Z^m\{(0,...,0)}, we map J to K via the map sending (T,z) to (b,z), where
b:=(b_1,...,b_m) (usually the b_i's are distinct primes), so that K is an ideal in
Q(z)[X]. For such a rational point b, we compute a Groebner basis G_b of K using
modular algorithms [1], where prime numbers are replaced by maximal ideals of the form
<zz_i>, and univariate rational interpolation [2,7]. Note that since Q[z]/<zz_i> = Q
we also use (if required) modular algorithms [1] over Q. The procedure is repeated for
many rational points b until their number is sufficiently large to recover the correct
coefficients in Q(T). Once we have these points, we obtain a set of polynomials G by
applying the sparse multivariate rational interpolation algorithm from [4] coefficient
wise to the list of Groebner bases G_b in Q(z)[X], where this algorithm makes use of
the following algorithms: univariate polynomial interpolation [2], univariate rational
function reconstruction [7], and multivariate polynomial interpolation [3]. The last
algorithm uses the wellknown Berlekamp/Massey algorithm [5] and its early termination
version [6]. The set G is then a Groebner basis of I with high probability.
 References:
 [1] E. A. Arnold: Modular algorithms for computing Groebner bases.
J. Symb. Comput. 35, 403419 (2003).
[2] R. L. Burden and J. D. Faires: Numerical analysis. 9th ed. (1993).
[3] M. BenOr and P. Tiwari: A deterministic algorithm for sparse multivariate
polynomial interpolation. Proc. of the 20th Annual ACM Symposium on
Theory of Computing, 301309 (1988).
[4] A. Cuyt and W.s. Lee: Sparse interpolation of multivariate rational functions.
Theor. Comput. Sci. 412, 14451456 (2011).
[5] E. Kaltofen and W.s. Lee: Early termination in sparse interpolation algorithms.
J. Symb. Comput. 36, 365400 (2003).
[6] E. Kaltofen, W.s. Lee and A. A. Lobo: Early termination in BenOr/Tiwari
sparse interpolation and a hybrid of Zippel's algorithm. Proc. ISSAC
(ISSAC '00), 192201 (2000).
[7] K. Sara and M. Monagan: Fast Rational Function Reconstruction. Proc. ISSAC
(ISSAC '06), 184190 (2006).
Procedures:
