Ask Your Question
0

Series Reversion [closed]

asked 2014-07-20 05:43:52 -0500

jjack gravatar image

updated 2014-07-22 22:52:04 -0500

I have a rational polynomial f(z) = P(z)/Q(z) which I want to revert (I'm trying to find the inverse of f).

I'll pick a simple example, the function f(z) = (1/z-1)^4

sage: z = var('z')
sage: z = FractionField(PolynomialRing(QQ, 'z')).objgen()
sage: f(z) = (1/z-1)^4
sage: g = f(z).taylor(x,2,4); #I expand f in a Taylor series just to make sure it's a series
sage: f.expand().reversion()
AttributeError: 'sage.symbolic.expression.Expression' object has no attribute 'reversion'
sage: g1 = g.power_series(QQ)
TypeError: denominator must be a unit
sage: g.reversion()
AttributeError: 'sage.symbolic.expression.Expression' object has no attribute 'reversion'

A power series P(z) is not a rational function P(z)/Q(z). That's why the TypeError: denominator must be a unit appears.

Is it possible that Sage cannot invert rational polynomials? In the Sage help text the computation is passed to pari first, before using Lagrangian inversion.

I thought Sage uses FLINT. Why is it not used for computing the inverse of rational polynomials?


EDIT (after answer by slelievre):

My problem now is how to avoid using the Symbolic Ring. I'm trying to do the following:

I want to invert an approximation formula for the Gamma-function. I am using Spouge's approximation. It has terms in the numerator and the denominator. Ithink I can compute them separetely per se, but what if I wanted to do a Pade approximation of the function

sage: f = (z+13)^(z+1/2)*exp(z)

sage: g = f.pade(5,5)

(it's the term ahead of the sum term in Spouge's approximation)

How do I get separate polynomials for the numerator and the denominator that are defined in the right algebraic structure?

edit retag flag offensive reopen merge delete

Closed for the following reason the question is answered, right answer was accepted by jjack
close date 2014-07-24 13:36:14.735427

Comments

I added a line with what y(x) looks like.

jjack gravatar imagejjack ( 2014-07-20 09:43:13 -0500 )edit

This is not sufficient for me to get the error. Please give all instructions necessary to reproduce the error from a fresh Sage session.

rws gravatar imagerws ( 2014-07-21 00:50:26 -0500 )edit

see above for the instructions

jjack gravatar imagejjack ( 2014-07-21 14:50:23 -0500 )edit

Note: (a) using Sage in a terminal will allow you to post your code by an easy copy-paste (where the input will have the `sage: ` prompt and the output will be mixed in with the input). (b) I edited the code in your original question to make it display better. (c) I added your follow-up question as an edit to your main question rather than an answer to your original question.

slelievre gravatar imageslelievre ( 2014-07-22 08:27:49 -0500 )edit

1 answer

Sort by ยป oldest newest most voted
1

answered 2014-07-21 22:48:40 -0500

updated 2014-07-22 22:53:55 -0500

Composition inverse of a power series

Avoid the symbolic ring whenever you can, and work in the appropriate structures for the task at hand.

A worked out example

Manipulate polynomials and rational fractions in the corresponding structures. Define

sage: R.<z> = PolynomialRing(QQ)
sage: F = FractionField(R)
sage: S = PowerSeriesRing(QQ,'z')

and check the result:

sage: R
Univariate Polynomial Ring in z over Rational Field
sage: F
Fraction Field of Univariate Polynomial Ring in z over Rational Field
sage: S
Power Series Ring in z over Rational Field

Define your polynomials and rational fractions.

sage: a = 2 - 3*z + z^2
sage: b = z
sage: f = a / b

Check what you get.

sage: a
z^2 - 3*z + 2
sage: b
z
sage: f
(z^2 - 3*z + 2)/z

Check out where a, b, f live.

sage: a.parent()
Univariate Polynomial Ring in z over Rational Field
sage: b.parent()
Univariate Polynomial Ring in z over Rational Field
sage: f.parent()
Fraction Field of Univariate Polynomial Ring in z over Rational Field

Want the Taylor series at 1 of f, and its inverse under composition.

First shift the function to move to the point 1.

sage: g = f(1+z)
sage: g
(z^2 - z)/(z + 1)

Move to the ring of power series.

sage: gg = S(g)

Check the first few terms.

sage: gg + O(z^8)
-z + 2*z^2 - 2*z^3 + 2*z^4 - 2*z^5 + 2*z^6 - 2*z^7 + O(z^8)

The constant term is zero and the coefficient of z is nonzero: we can compute the composition inverse.

sage: hh = gg.reversion()

Examine the first few terms.

sage: hh + O(z^8)
-z + 2*z^2 - 6*z^3 + 22*z^4 - 90*z^5 + 394*z^6 - 1806*z^7 + O(z^8)

Check that hh is the composition inverse of gg:

sage: hh(gg)
z + O(z^21)
sage: gg(hh)
z + O(z^21)

EDIT (to answer your follow-up question).

If you have a rational fraction g, then gnum = g.numerator() and gden = g.denominator() will get you the polynomials for the numerator and the denominator of g. You can check where they live by gnum.parent() and gden.parent(). If you need them to be in R, F or S, just do R(gnum), F(gnum), S(gnum), etc.

edit flag offensive delete link more

Comments

Thanks a lot for your help. I think I got the hang of it.

jjack gravatar imagejjack ( 2014-07-23 13:38:52 -0500 )edit

@jjack: Happy to help. If this answers your question, click the tick-mark (below score at top-left of answer) to indicate that. This marks the question as "answered" in the list of questions. When you reach 15 karma, you can also upvote answers and questions.

slelievre gravatar imageslelievre ( 2014-07-23 15:04:29 -0500 )edit

Question Tools

1 follower

Stats

Asked: 2014-07-20 05:43:52 -0500

Seen: 392 times

Last updated: Jul 23 '14