Ideal(PolynomialRing(RationalField(),['x','y']),-1).primary_decomposition() returns [] as expected but Ideal(PolynomialRing(RationalField(),['x','y']),1).primary_decomposition() returns [Ideal (1) of Multivariate Polynomial Ring in x, y over Rational Field].
Clearly the ideals are identical, so they should have the same primary decomposition (independent of the unit it is generated by). Either, both should return [] or [Ideal(1)]. One would also expect the output of primary decomposition to not include (1) similarly to the prime factorization of a natural number not including 1, e.g., list(factor(1)) returns []. -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-support+unsubscr...@googlegroups.com. To post to this group, send email to sage-support@googlegroups.com. Visit this group at http://groups.google.com/group/sage-support. For more options, visit https://groups.google.com/groups/opt_out.