I am interested in factoring (sums of) Schubert polynomials using Sage (which uses Singular). Here's my input attempt:
> X = SchubertPolynomialRing(QQ)
> f = X([5,4,2,6,1,3]).expand(); f
> f.factor()
However, while the second line gives a correct output, x0^4*x1^3*x2^2*x3 + x0^4*x1^3*x2*x3^2, when I try to factor this output, I get the following error message: "NotImplementedError: Factorization of multivariate polynomials over non-fields is not implemented."
I don't know how to get Sage to recognize a SchubertPolynomialRing as being over a field. Any suggestions?
