Ask Your Question

find quotient of two multivariate polynomials (which are divisible)

asked 2020-07-22 15:35:21 -0600

Daniel Friedan gravatar image

updated 2020-07-23 08:22:18 -0600

slelievre gravatar image

I have two multivariate polynomials, num and denom such that denom divides num. But I have not been able to get SageMath to simplify the quotient num/denom.

The worksheet and the two polynomials are available at

Sample code is below.

sage: dR.<d1, d2, d3, d4, d12, d13, d14, d23, d24, d34, d123, d124, d134, d234, d1234> = PolynomialRing(ZZ, 15, order='lex')

sage: denom = load('denom')
sage: num = load('num')
sage: denom.divides(num)
sage: num.reduce(Ideal([denom]))
sage: F = num/denom
sage: F.denominator()/denom
sage: F.reduce()
sage: F.denominator() / denom
sage: num.number_of_terms()
sage: denom.number_of_terms()
edit retag flag offensive close merge delete



As pointed out in the answers, the data was created with one term ordering, while the workbook uses another term ordering.

However, why does the term ordering affect the truth of 'num is contained in the ideal generated by denom' i.e., the value of num.reduce(Ideal[denom]) ? One of the answers refers to "This 'local' term ordering ``. I don't understand how the term ordering affects the structure of the ideals in a polynomial ring over the integers.

Daniel Friedan gravatar imageDaniel Friedan ( 2020-07-25 16:14:09 -0600 )edit

This is an unfortunate consequence of the fact that PolynomialRing uses Singular's ring in the background, without restricting the term ordering. Singular's ring is an ordinary polynomial ring when given a "global" term ordering, but it is rather a localization when given a "local" term ordering. These local orderings should be disallowed on a PolynomialRing, and the Singular functionality for localizations should be exposed in another way.

rburing gravatar imagerburing ( 2020-07-27 05:13:27 -0600 )edit

2 answers

Sort by ยป oldest newest most voted

answered 2020-07-23 06:42:38 -0600

rburing gravatar image

updated 2020-07-23 07:00:48 -0600

Your data is fishy:

sage: denom.parent().term_order()
Negative weighted degree reverse lexicographic term order with weights (1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 4)

This 'local' term ordering causes reduce and (hence) divides to behave differently. In particular, it does not give the same result as with a well-ordering such as 'lex' in your ring dR:

sage: dR(denom).divides(dR(num))

So denom does not divide num in the ordinary sense.

Edit: In the 'local' sense, we have that denom is invertible, so the ideal generated by denom is the unit ideal, so of course num.reduce(Ideal([denom])) is zero.

edit flag offensive delete link more

answered 2020-07-22 18:36:27 -0600

nbruin gravatar image

There really should be a version of reduce somewhere that gives you which combination of the generators allowed for the reduction, and if someone can find where then that really should be the preferred solution, but at the expense of some extra variables, you can fake it:

Given two polynomials f and g in variables [z1,...,zn], embed the ring into a larger one with variables [z1,...,zn,F,G], with an elimination order for z1,...zn (i.e., those variables should be larger than F and G). Then reduce f-F with respect to the principal ideal generated by g-G. The result should be a polynomial that is free of z1,...,zn, and expresses F as a polynomial in G.

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: 2020-07-22 15:35:21 -0600

Seen: 202 times

Last updated: Jul 23