Ask Your Question

Find polynomial in terms of ideal

asked 2011-09-26 19:24:55 -0500

bmanden gravatar image


I have been reading the tutorials and experimenting with the polynomial ring code but I have not been able to find out how to express a given polynomial in terms of the generators of an ideal. For instance, if I define a few polynomial generators of an ideal, say 'f=xyz-y^2' and 'g=x^3 + z', and have a polynomial p, I would like to express p as 'p = uf + vg', where I need to find the polynomials u,v. I tried a bunch of solve([],u,v) commands and such but I haven't found the magic solution. Thanks.

edit retag flag offensive close merge delete

2 answers

Sort by » oldest newest most voted

answered 2011-09-27 16:01:41 -0500

Mike Hansen gravatar image

You'll want to use the lift method on the polynomial p. For example,

sage: R.<x,y,z> = QQ[]
sage: f = x*y*z-y^2
sage: g = x^3 + z
sage: I = R.ideal([f,g])
sage: p = 3*f-x*g
sage: p.lift(I)
[3, -x]

That gives you the coefficients in front of each of the generators of the ideal. You'll want to make sure that p is actually in the ideal first; otherwise, you might get something like:

sage: x in I
sage: x.lift(I)
[0, 0]
edit flag offensive delete link more

answered 2011-09-27 09:39:31 -0500

Dirk Danckaert gravatar image

I have a hunch that you need to specify the needed polynomials u, v better. Surely there's an infinity of solutions for them? Or am I missing something? I would be very surprised if a solve() algorithm could handle this infinity of solutions. In the case of polynomials with a single variable these coefficients are usually found with the help of Bezout's theorem. I guess the multivariate case is handled analogously.

edit flag offensive delete link more

Your Answer

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

Add Answer

Question Tools

1 follower


Asked: 2011-09-26 19:24:55 -0500

Seen: 271 times

Last updated: Sep 27 '11