Is it enough to use the primary_decomposition method of ideals?  E.g.:

R.<x,y,z> = PolynomialRing(QQ,3)
i1 = R.ideal([x*z,y*z])
[list(_.gens()) for _ in i1.primary_decomposition()]

gives

[[z], [y, x]]

-Marshall Hampton


On Nov 24, 10:26 am, andrew ewart <[email protected]> wrote:
> I want to be able to right a little code that would take any
> projective variety and break it down into its irreducible variety
> components
> any suggestions.

-- 
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to 
[email protected]
For more options, visit this group at 
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org

Reply via email to