Ask Your Question

Implementing different groebner_basis() algorithms

asked 2010-08-23 15:28:46 +0200

Chris Aholt gravatar image

updated 2011-04-28 16:33:49 +0200

Kelvin Li gravatar image

I'm not sure if this is the correct forum to ask this type of question, but I'll try anyway.

I am curious if there is a way to do the following in Sage. I would like a way to make a Groebner basis calculation to not worry about calculations involving polynomials of too high degree. Macaulay2's gb command has the ability to specify this with an option called HardDegreeLimit, for instance. The Sage documentation on groebner_basis seems to imply that I can use Macaulay2's gb command with this option, but I don't know how to make it work. Singular's std command also has a degBound option, but I again don't know how to use this in my Sage implementation.

edit retag flag offensive close merge delete

1 Answer

Sort by ยป oldest newest most voted

answered 2010-08-23 17:02:27 +0200

Mike Hansen gravatar image

With the current Groebner basis commands in Sage, there is not a good way to do this. None of the underlying methods like _groebner_basis_macaulay2 support a degree limit.

As a workaround, you can do something like the following:

sage: P.<a,b,c> = PolynomialRing(QQ,3, order='lex')
sage: I = sage.rings.ideal.Katsura(P,3)
sage: I.groebner_basis()  #compute the normal Groebner basis
[a - 60*c^3 + 158/7*c^2 + 8/7*c - 1, b + 30*c^3 - 79/7*c^2 + 3/7*c, c^4 - 10/21*c^3 + 1/84*c^2 + 1/84*c]
sage: singular.eval('degBound = 2;')  #set the degree bound in Singular
'degBound = 2;'
sage: gb = Sequence(map(P, singular(I).std())); gb
[10*b*c - b + 12*c^2 - 4*c, 4*b^2 + 2*b*c - b, a + 2*b + 2*c - 1]
sage: singular.eval('degBound = 0;')  #reset it back to 0 (unlimited)
'degBound = 0;'

I've made this ticket #9789.

edit flag offensive delete link more


I agree with your answer, but I think your example is misleading. The Singular manual states that "degBound should not be used for a global ordering with inhomogeneous input". I'm pretty sure that the output for a non-homogeneous ideal is potentially not the degree truncation of an actual GB.

Volker Braun gravatar imageVolker Braun ( 2010-08-23 21:57:44 +0200 )edit

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools


Asked: 2010-08-23 15:28:46 +0200

Seen: 545 times

Last updated: Aug 23 '10