Try this:
sage: I3_red_I2 = R.ideal([p.reduce(I2gb) for p in I3gb])
regards
john perry
On Sep 28, 12:24 am, Vinay Wagh <wagh...@gmail.com> wrote:
> Suppose I have two ideals I & J in k[X_1,\cdots,x_n], where k is a
> field. How do I reduce an ideal I wrt ideal J.
>
> e.g. Singular provides me a command
>
> singular > reduce(I,std(J));
>
> Without moving back and forth to Singular, is it possible to implement
> this in sage?
>
> I tried the following code:
>
> sage: R.<X,Y,Z> = PolynomialRing(QQ,3,order=TermOrder('wdeglex',[4,6,11]));
> sage: I = R.ideal(X^2-Y^3+Z^4, Y^5-Z^6+X^7, Z^13-X^12+Y^11);
> sage: I2 = I*I;
> sage: I3 = I2*I;
> sage: I2gb = I2.groebner_basis();
> sage: I3gb = I3.groebner_basis();
> sage: I2gb
> sage: I3_red_I2 = reduce(I3, I2gb);
>
> The last command (redece) is giving me an error. I am not getting what
> wrong I am doing...
>
> Thanks and regards
>
> -- VInay Wagh
--
To post to this group, send email to sage-support@googlegroups.com
To unsubscribe from this group, send email to
sage-support+unsubscr...@googlegroups.com
For more options, visit this group at
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org
