recognize a rational function is a polynomial

asked 2022-02-25 17:28:49 +0100

rue82
113 ●1 ●4 ●13

updated 2022-02-25 18:42:08 +0100

slelievre
17704 ●22 ●160 ●350 http://carva.org/samue...

Is it possible to make sage recognize that a certain rational function is a Laurent polynomial, and treat it as such?

A simple example:

sage: R = LaurentPolynomialRing(ZZ, ['q1'])
sage: R.inject_variables(verbose=False)
sage: f = (1/(1-q1) + 1/(1-q1^-1))
sage: parent(f)
Fraction Field of Univariate Polynomial Ring in q1 over Integer Ring

Here, f lives correctly in the fraction field. However, it is in particular a polynomial, to which I'd like to apply the polynomial methods, e.g. get its monomial list.

Is this possible to achieve?

answered 2022-02-25 18:44:52 +0100

slelievre
17704 ●22 ●160 ●350 http://carva.org/samue...

updated 2022-02-27 14:52:37 +0100

Good thing you gave the Laurent polynomial ring a name.

That makes it easy to call it on an element of the fraction field.

Here is an example using the element f in the question.

sage: ff = R(f)
sage: parent(f)
Univariate Laurent Polynomial Ring in q1 over Integer Ring

In more complicated cases, it helps to use the numerator and denominator.

sage: g = (1-q1) * (1-q2) * (q2/(1-q1)+q2^-1/(1-q1^-1))
sage: g
(q2^3 - q1*q2 - q2^2 + q1)/(-q2)

sage: gg = R(g)
Traceback (most recent call last)
...
TypeError: fraction must have unit denominator

sage: gg = R(g.numerator()) / R(g.denominator())
sage: gg
-q2^2 + q1 + q2 - q1*q2^-1
sage: parent(gg)
Multivariate Laurent Polynomial Ring in q1, q2 over Integer Ring

edit flag offensive delete link

Comments

Thanks, that's also what I thought, but why is it not working e.g. here?

R.<q1,q2> = LaurentPolynomialRing(ZZ,2)
def test():
    x = (1-q1) * (1-q2) * (q2/(1-q1)+q2^-1/(1-q1^-1))
    return R(x)

rue82 ( 2022-02-26 15:51:53 +0100 )edit

the error it gives is: fraction must have unit denominator.

rue82 ( 2022-02-26 15:53:12 +0100 )edit

Thank you!

rue82 ( 2022-02-26 22:14:14 +0100 )edit

add a comment

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

recognize a rational function is a polynomial

1 Answer

Comments

Your Answer

Question Tools

Stats

Related questions

recognize a rational function is a polynomial edit

1 Answer

Comments

Your Answer

Question Tools

Stats

Related questions

recognize a rational function is a polynomial