You can try to compute over ZZ with p added as a generator, because
J in Z^n/p ~= J + p in Z^n
sage: p = next_prime(2^100)
sage: R.x,y,z = ZZ[]
sage: I = Ideal(R.random_element() for _ in range(R.ngens()))
sage: I += [p]
sage: I.groebner_basis()
[x^2 + 380295180068468820449010961696*x*z
Hi!
On 2012-12-13, Santanu Sarkar sarkar.santanu@gmail.com wrote:
When I want to calculate
Groebner basis, I have following error.
verbose 0 (3292: multi_polynomial_ideal.py, groebner_basis) Warning:
falling back to very slow toy implementation.
This is not an error but a warning.